Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
Applied and Environmental Microbiology
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
Environmental Microbiology

Expanding the Diversity of Bacterioplankton Isolates and Modeling Isolation Efficacy with Large-Scale Dilution-to-Extinction Cultivation

Michael W. Henson, V. Celeste Lanclos, David M. Pitre, Jessica Lee Weckhorst, Anna M. Lucchesi, Chuankai Cheng, Ben Temperton, J. Cameron Thrash
Robert M. Kelly, Editor
Michael W. Henson
aDepartment of Biological Sciences, University of Southern California, Los Angeles, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Michael W. Henson
V. Celeste Lanclos
aDepartment of Biological Sciences, University of Southern California, Los Angeles, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
David M. Pitre
bDepartment of Biological Sciences, Louisiana State University, Baton Rouge, Louisiana, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jessica Lee Weckhorst
bDepartment of Biological Sciences, Louisiana State University, Baton Rouge, Louisiana, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Anna M. Lucchesi
bDepartment of Biological Sciences, Louisiana State University, Baton Rouge, Louisiana, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Chuankai Cheng
aDepartment of Biological Sciences, University of Southern California, Los Angeles, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Ben Temperton
cSchool of Biosciences, University of Exeter, Exeter, United Kingdom
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ben Temperton
J. Cameron Thrash
aDepartment of Biological Sciences, University of Southern California, Los Angeles, California, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for J. Cameron Thrash
Robert M. Kelly
North Carolina State University
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/AEM.00943-20
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Cultivated bacterioplankton representatives from diverse lineages and locations are essential for microbiology, but the large majority of taxa either remain uncultivated or lack isolates from diverse geographic locales. We paired large-scale dilution-to-extinction (DTE) cultivation with microbial community analysis and modeling to expand the phylogenetic and geographic diversity of cultivated bacterioplankton and to evaluate DTE cultivation success. Here, we report results from 17 DTE experiments totaling 7,820 individual incubations over 3 years, yielding 328 repeatably transferable isolates. Comparison of isolates to microbial community data for source waters indicated that we successfully isolated 5% of the observed bacterioplankton community throughout the study; 43% and 26% of our isolates matched operational taxonomic units and amplicon single-nucleotide variants, respectively, within the top 50 most abundant taxa. Isolates included those from previously uncultivated clades such as SAR11 LD12 and Actinobacteria acIV, as well as geographically novel members from other ecologically important groups like SAR11 subclade IIIa, SAR116, and others, providing isolates in eight putatively new genera and seven putatively new species. Using a newly developed DTE cultivation model, we evaluated taxon viability by comparing relative abundance with cultivation success. The model (i) revealed the minimum attempts required for successful isolation of taxa amenable to growth on our media and (ii) identified possible subpopulation viability variation in abundant taxa such as SAR11 that likely impacts cultivation success. By incorporating viability in experimental design, we can now statistically constrain the effort necessary for successful cultivation of specific taxa on a defined medium.

IMPORTANCE Even before the coining of the term “great plate count anomaly” in the 1980s, scientists had noted the discrepancy between the number of microorganisms observed under the microscope and the number of colonies that grew on traditional agar media. New cultivation approaches have reduced this disparity, resulting in the isolation of some of the “most wanted” bacterial lineages. Nevertheless, the vast majority of microorganisms remain uncultured, hampering progress toward answering fundamental biological questions about many important microorganisms. Furthermore, few studies have evaluated the underlying factors influencing cultivation success, limiting our ability to improve cultivation efficacy. Our work details the use of dilution-to-extinction (DTE) cultivation to expand the phylogenetic and geographic diversity of available axenic cultures. We also provide a new model of the DTE approach that uses cultivation results and natural abundance information to predict taxon-specific viability and iteratively constrain DTE experimental design to improve cultivation success.

INTRODUCTION

Axenic cultures of environmentally important microorganisms are critical for fundamental microbiological investigation, including generating physiological information about environmental tolerances, determining organism-specific metabolic and growth rates, testing hypotheses generated from in situ omics observations, and experimentally examining microbial interactions. Research using important microbial isolates has been critical to a number of discoveries, such as defining microorganisms involved in surface ocean methane saturation (1–3), the role of proteorhodopsin in maintaining cellular functions during states of carbon starvation (4, 5), the complete nitrification of ammonia within a single organism (6), and identifying novel metabolites and antibiotics (7–10). However, the vast majority of taxa remain uncultivated (11–13), restricting valuable experimentation on such topics as genes of unknown function, the role of analogous gene substitutions in overcoming auxotrophy, and the multifaceted interactions occurring in the environment inferred from sequence data (11, 14–16).

The quest to bring new microorganisms into culture and the recognition that traditional agar plate-based approaches have limited success (17–19) have compelled numerous methodological advances spanning a wide variety of techniques, such as diffusion chambers, microdroplet encapsulation, and slow acclimatization of cells to artificial media (20–25). Dilution-to-extinction (DTE) cultivation using sterile seawater as the medium has also proven highly successful for isolating bacterioplankton (26–32). Pioneered by Don Button and colleagues for the cultivation of oligotrophic bacteria, this method essentially preisolates organisms after serial dilution by separating individual or small groups of cells into their own incubation vessel (32, 33). This prevents slow-growing, obligately oligotrophic bacterioplankton from being outcompeted by faster-growing organisms, as would occur in enrichment-based isolation methods like those that would target aerobic heterotrophs. DTE is also a practical method for taxa that cannot grow on solid media. Natural seawater media provide these taxa with the same chemical surroundings from which they are collected, reducing the burden of anticipating all the relevant compounds required for growth (33).

Improvements to DTE cultivation in multiple labs have increased the number of inoculated wells and decreased the time needed to detect growth (26, 28, 34), thereby earning the moniker “high-throughput culturing” (26, 28). We (35) and others (30) have also adapted DTE culturing by incorporating artificial media in place of natural seawater media to successfully isolate abundant bacterioplankton. Thus far, DTE culturing has led to isolation of many important marine and freshwater groups, such as marine SAR11 Alphaproteobacteria (28, 29, 34–36), the freshwater SAR11 LD12 clade (29), SUP05/Arctic96BD-19 Gammaproteobacteria (37–39), OM43 Betaproteobacteria (26, 27, 31, 40, 41), HIMB11-type Roseobacter spp. (35, 42), numerous so-called oligotrophic marine Gammaproteobacteria (43), and acI Actinobacteria (44).

Despite the success of DTE cultivation, many taxa continue to elude domestication (11–13, 16). Explanations include a lack of required nutrients or growth factors in media (20, 45–49) and biological phenomena such as dormancy and/or phenotypic heterogeneity within populations (47, 48, 50–56). However, there have been few studies empirically examining the factors underlying isolation success in DTE cultivation experiments (34, 57, 58), restricting our ability to determine the relative importance of methodological versus biological influences on cultivation reliability for any given organism. Moreover, even for those taxa that we have successfully cultivated, in many cases we lack geographically diverse strains, restricting comparisons of the phenotypic and genomic diversity that may influence taxon-specific cultivability.

We undertook a 3-year cultivation effort in the coastal northern Gulf of Mexico (nGOM), from which we lack representatives of many common bacterioplankton groups, to provide new model organisms for investigating microbial function, ecology, biogeography, and evolution. Simultaneously, we paired our cultivation efforts with 16S rRNA gene amplicon analyses to compare cultivation results with the microbial communities in the source waters. We have previously reported on the success of our artificial media in obtaining abundant taxa over the course of the first seven experiments from this campaign (35). Here, we expand our report to include cultivation results from a total of 17 experiments, and we update the classic viability calculations of Button et al. (33) with a new model to estimate the viability of individual taxa using relative-abundance information. New isolates belonged to cultivated groups in eight putatively novel genera and seven putatively novel species in previously cultivated genera and expanded cultured geographic representation for many important clades like SAR11. Additionally, using model-based predictions, we identified possible taxon-specific viability variation that can influence cultivation success. By incorporating these new viability estimates into the model, our method facilitates statistically informed experimental design for targeting individual taxa, thereby reducing uncertainty for future culturing work (59).

RESULTS

General cultivation campaign results.We conducted a total of 17 DTE cultivation experiments to isolate bacterioplankton (sub-2.7-μm fraction), with paired microbial community characterization of source waters (0.22- to 2.7-μm fraction), from six coastal Louisiana sites over a 3-year period (see Table S1 in the supplemental material, available at https://doi.org/10.6084/m9.figshare.12142113). We inoculated 7,820 distinct cultivation wells (all experiments) with an estimated 1 to 3 cells · well−1 using overlapping suites of artificial seawater media, JW (years 1 and 2 [35]) and MWH (year 3), designed to match the natural environment (Table 1). The MWH suite of media was modified from the JW media to additionally include choline, glycerol, glycine betaine, cyanate, dimethyl sulfoxide (DMSO), dimethylsulfoniopropionate (DMSP), thiosulfate, and orthophosphate (Table S1). These compounds have been identified as important metabolites and osmolytes for marine and freshwater microorganisms and were absent in the first iteration (JW) of our media (60–66). A total of 1,463 wells were positive (>104 cells · ml−1), and 738 of these were transferred to 125-ml polycarbonate flasks. For four experiments (FWC, FWC2, JLB2, and JLB3) we transferred only a subset of positives (48/301, 60/403, 60/103, and 60/146) because the number of isolates exceeded our ability to maintain and identify them at that time (Table 1). The subset of positive wells for these four experiments was selected using flow cytometry signatures usually indicative of smaller oligotrophic cells like SAR11 strain HTCC1062 (49) using our settings. Of the 738 wells from which we transferred cells across all experiments, 328 of these yielded repeatably transferable isolates that we deemed pure cultures based on 16S rRNA gene PCR and Sanger sequencing.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 1

Cultivation statistics, including whole-community viability estimates

Phylogenetic and geographic novelty of our isolates.The 328 isolates belonged to three phyla: Proteobacteria (n = 319), Actinobacteria (n = 8), and Bacteroidetes (n = 1) (see Fig. S1 to S5 in the supplemental material). We placed these isolates into 55 groups based on their positions within 16S rRNA gene phylogenetic trees (Fig. S1 to S5) and as a result of having ≥94% 16S rRNA gene sequence identity to other isolates. We applied a nomenclature to each group based on previous 16S rRNA gene database designations and/or other cultured representatives (Fig. 1; Table S1 [available at https://doi.org/10.6084/m9.figshare.12142113]). Isolates represented eight putatively novel genera with ≤94.5% 16S rRNA gene identity to a previously cultured representative: the Actinobacteria acIV subclades A and B and one other unnamed Actinobacteria group, an undescribed Acetobacteraceae clade (Alphaproteobacteria), the freshwater SAR11 LD12 (“Candidatus Fonsibacter ubiquis” [29]), the MWH-UniPo and an unnamed Burkholderiaceae clade (Betaproteobacteria), and the OM241 Gammaproteobacteria (Fig. 1; Table S1). Seven additional putatively novel species in other genera were also isolated (between 94.6 and 96.9% 16S rRNA gene sequence identity) in unnamed Comamonadaceae and Burkholderiales clades (Betaproteobacteria), the SAR92 clade and Pseudohongiella genus (Gammaproteobacteria), and unnamed Rhodobacteraceae and Bradyrhizobiaceae clades, as well as Maricaulis spp. (Alphaproteobacteria) (Fig. 1). Louisiana State University Culture Collection (LSUCC) isolates belonging to the groups BAL58 Betaproteobacteria (Fig. S4), OM252 Gammaproteobacteria, HIMB59 Alphaproteobacteria, and what we designated the LSUCC0101-type Gammaproteobacteria (Fig. S5) had close 16S rRNA gene matches to other isolates at the species level; however, none of those previously cultivated organisms have been formally described (Fig. 1). The OM252, BAL58, and MWH-UniPo clades were the most frequently cultivated, with 124 of our 328 isolates belonging to these three groups (Table S1). In total, 73 and 10 of the 328 isolates belonged in putatively novel genera and novel species in previously cultivated genera, respectively. We estimated that at least 310 of these isolates were geographically novel, being the first of their type cultivated from the nGOM (Fig. 2). This included isolates from cosmopolitan groups like SAR11 subclade IIIa, OM43 Betaproteobacteria, SAR116, and HIMB11-type “Roseobacter” spp. Cultivars from Vibrio sp. and Alteromonas sp. were the only two groups with close relatives (species level) isolated from the GOM.

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Percent identity of LSUCC isolate 16S rRNA genes compared with those from other isolates in NCBI (“Other,” gray dots) or from the DTE culture collections IMCC (gold dots), HTCC (blue dots), and HIMB (green dots). Each dot represents a pairwise 16S rRNA gene comparison (via BLASTn). x-axis categories are groups designated according to ≥94% sequence identity and phylogenetic placement (see Fig. S1 to S5 in the supplemental material). Above the graph is the 16S rRNA gene sequence percent identity to the closest non-LSUCC isolate within a column. Groups colored in red are those where LSUCC isolates represent putatively novel genera, whereas orange indicates putatively novel species.

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

A global map of the isolation locations of isolates from selected important aquatic bacterioplankton clades. All depicted taxa were isolated from surface water (<1 to 20 m), or the depth of the sample was not reported (see Table S1 in the supplemental material [available at https://doi.org/10.6084/m9.figshare.12142113] for details). Circles represent LSUCC isolates, while triangles are non-LSUCC isolates. Inset, a zoomed view of the coastal Louisiana region where LSUCC bacterioplankton originated. The map was made with the R package ggplot2 using the command map_data (see https://github.com/thrash-lab/Cultivar-novelty-map).

Natural abundance of isolates.We matched LSUCC isolate 16S rRNA gene sequences with both operational taxonomic units (OTUs) and amplicon single-nucleotide variants (ASVs) from bacterioplankton communities to assess the relative abundances of our isolates in the coastal nGOM waters that served as inocula. While OTUs provide a broad group-level designation (97% sequence identity), this approach can artificially combine multiple ecologically distinct taxa (67). Due to higher stringency for defining a matching 16S rRNA gene, ASVs can increase the confidence that our isolates represent environmentally relevant organisms (68, 69). However, while many abundant oligotrophic bacterioplankton clades, such as SAR11 (29, 70), OM43 (40, 41), SAR116 (71), and Sphingomonas spp. (72), have a single copy of the rRNA gene operon, other taxa can have multiple rRNA gene copies (70, 73), complicating ASV analyses. Since we could not a priori rule out multiple rRNA gene operons for novel groups with no genome-sequenced representatives, we used both OTU and ASV approaches.

In total, we obtained at least one isolate from 40 of the 777 OTUs and 71 of the 1,323 ASVs observed throughout the 3-year data set. Forty-three percent and 26% of LSUCC isolates matched the top 50 most abundant OTUs (median relative abundances for all sites, from 8.1 to 0.11% [see Fig. S6A in the supplemental material]) and ASVs (mean relative abundances for all sites, from 3.8 to 0.11% [Fig. S6B]), respectively, across all sites and samples. Microbial communities from all collected samples clustered into two groups corresponding to those inhabiting salinities below 7 and above 12, and salinity was the primary environmental driver distinguishing community beta diversity (OTU, R2 = 0.88, P = 0.001; ASV, R2 = 0.89, P = 0.001). As part of the cultivation strategy after the first five experiments, we used a suite of five media differing by salinity and matched the experiment with the medium that most closely resembled the salinity at the sample site. Consequently, our isolates matched abundant environmental groups from both high- and low-salinity regimes. At salinities above 12, LSUCC isolates matched 13 and 14 of the 50 most abundant OTUs and ASVs, respectively (Fig. 3A and 4A; Table S1 [available at https://doi.org/10.6084/m9.figshare.12142113]). These taxa included the abundant SAR11 subclade IIIa.1, HIMB59, HIMB11-type “Roseobacter,” and SAR116 Alphaproteobacteria; the OM43 Betaproteobacteria; and the OM182 and LSUCC0101-type Gammaproteobacteria. At salinities below seven, 10 and 9 of the 50 most abundant OTUs and ASVs, respectively, were represented by LSUCC isolates, including one of the most abundant taxa in both cluster sets, SAR11 LD12 (Fig. 3B and 4B). Some taxa, such as SAR11 IIIa.1 and OM43, were among the top 15 most abundant taxa in both salinity regimes (Fig. 3 and 4; Table S1), suggesting a euryhaline lifestyle. In fact, our cultured SAR11 IIIa.1 ASV7471 was the most abundant ASV in the aggregate data set (Fig. S6).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Rank abundances of the 50 most abundant OTUs from all sites based on median relative abundance at salinities greater than 12 (A) and less than 7 (B). The boxes indicate the interquartile range (IQR) of the data, with vertical lines indicating the upper and lower extremes according to 1.5 × IQR. Horizontal lines within each box indicate the median. The data points comprising the distribution are plotted on top of the boxplots. The shade of the dot represents the salinity at the sample site (red to blue, lower to higher), while the color of the box indicates broad taxonomic identity. LSUCC labels indicate OTUs with at least one cultivated representative.

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Rank abundances of the 50 most abundant ASVs from all sites based on median relative abundance at salinities greater than 12 (A) and less than 7 (B). The boxes indicate the interquartile range (IQR) of the data, with vertical lines indicating the upper and lower extremes according to 1.5 × IQR. Horizontal lines within each box indicate the median. The data points comprising the distribution are plotted on top of the boxplots. The shade of the dot represents the salinity at the sample site (red to blue, lower to higher), while the color of the box indicates broad taxonomic identity. LSUCC labels indicate ASVs with at least one cultivated representative.

Overall, this effort isolated taxa representing 18 and 12 of the top 50 most abundant OTUs and ASVs, respectively (Table 2; Fig. S6). When categorized by different median relative abundance categories of >1%, 0.1% to 1%, and <0.1%, isolate OTUs were distributed across those categories at 15%, 20%, and 27%, and isolate ASVs were distributed accordingly at 4%, 26%, and 37% (Table 2). Isolates with median relative abundances of <0.1%, such as Pseudohongiella spp., Rhodobacter spp., and Bordetella spp., would canonically fall within the rare biosphere (74) (Table S1 [available at https://doi.org/10.6084/m9.figshare.12142113]). A number of isolates did not match any identified OTUs or ASVs (38% and 33% of LSUCC isolates compared to available OTUs and ASVs, respectively), either because their matching OTUs/ASVs were below our thresholds for inclusion (at least two reads from at least two sites) or because they were below the detection limit from our sequencing effort (Table 2). Thus, 43% and 30% of our isolates belonged to OTUs and ASVs, respectively, with median relative abundances of >0.1%.

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 2

Median relative abundances of cultured OTUs and ASVs across all samples

Modeling DTE cultivation.An enigma that became immediately apparent through a review of our data was the absence of an obvious relationship between the abundance of a given taxon in the inoculum and the frequency of obtaining an isolate of the same type from a DTE cultivation experiment (see Fig. S7 and S8 in the supplemental material). For example, although we could culture SAR11 LD12 over a range of medium conditions (29) and the matching ASV had relative abundances of >5% in six of our 17 experiments (Fig. 5), we isolated only one representative (LSUCC0530). In an ideal DTE cultivation experiment where cells are randomly subsampled from a Poisson-distributed population, if the medium is sufficient for a given microorganism’s growth, then the number of isolates should correlate with that microorganism’s abundance in the inoculum. However, a qualitative examination of several abundant taxa that grew in our media, some of which we cultured on multiple occasions, revealed no clear pattern between abundance and isolation success (Fig. 5). Considering that medium composition was sufficient for cultivation of these organisms on at least some occasions, we hypothesized that cultivation frequency may reflect differences in the capacity for growth within populations of a given taxon. Thus, we decided to model cultivation frequency in relationship to estimated abundances in a way that could generate estimates of cellular viability, defined herein as meaning “presently able to grow in defined medium,” as opposed to a broader definition equating viability with being alive more generally, since we evaluated only growth capacity in this study. We hoped that modeling might also help us inform experimental design and make DTE cultivation efforts more predictable (59).

FIG 5
  • Open in new tab
  • Download powerpoint
FIG 5

Relative abundance of ASVs within key taxonomic groups compared with salinity. ASV types are colored independently, and triangle points indicate experiments for which at least one isolate was obtained. Nonlinear regression lines are provided as a visual aid for abundance trends.

Previously, Don Button and colleagues developed a statistical model for viability (V) of cells in the entire population for a DTE experiment (33): V=−ln(1−p)λ(1)

Where p is the proportion of wells or tubes, n, with growth z (p = z/n) and λ is the estimated number of cells inoculated per well (the authors used X originally). The equation uses a Poisson distribution to account for the variability in cell distribution within the inoculum and therefore the variability in the number of wells or tubes receiving the expected number of cells. We and others have used this equation in the past (26, 28, 35) to evaluate the efficacy of our cultivation experiments in the context of commonly cited numbers for cultivability using agar plate-based methods (13, 17, 75).

While equation 1 was effective for its intended purpose, it has a number of drawbacks that limit its utility for taxon-specific application: (i) if p = 1, i.e., all wells are positive, then the equation is invalid; (ii) At high values of p and low values of λ, estimates of V can exceed 100% (Table 1); (iii) accuracy of viability, calculated by the asymptotic standard error (ASE) or the coefficient of variation (CV), was shown to be nonuniform across a range of λ, with greatest accuracy when true viability was ∼10% (33), and thus, low viability, low values of λ, and small values of n were found to produce unreliable results; (iv) if p = 0, i.e., no positive wells are observed, estimates of viability that could produce 0 positive wells cannot be calculated; and (v) Button’s original model assumes that a well will produce a pure culture only if the inoculated well contains one cell, but in contrast, in low-diversity samples, samples dominated by a single taxon, or experiments evaluating viability from axenic cultures across different media, a limitation that only wells with single cells are axenic will underestimate the expected number of pure wells.

To overcome these limitations, we developed a Monte Carlo simulation model that facilitates the incorporation of relative-abundance data from complementary community profiling data (e.g., 16S rRNA gene amplicons) to calculate the likelihood of positive wells, pure wells, and viability at a taxon-specific level, based on the observed number of wells for which we obtained an isolate of a particular taxon (Fig. 6). By employing a Monte Carlo approach, our model is robust across all values of p and n with uniform prediction accuracy, and we can estimate the accuracy of our prediction within 95% confidence intervals (CI). Furthermore, the width of 95% CI boundaries of viability, as well as the expected number of positive and pure wells, is entirely controllable and dependent only on available computational capacity for bootstrapping (i.e., these can be improved with more bootstrapping, but at greater computational cost). When zero positive wells are observed experimentally, our approach enables estimation of a maximum viability that could explain such an observation by identifying the range of variability values for which zero resides within the bootstrapped 95% CI. Finally, the ability to calculate the viability of the entire community, as in equation 1, is retained simply by estimating viability using a relative abundance of one.

FIG 6
  • Open in new tab
  • Download powerpoint
FIG 6

Graphical depiction of the viability model.

We compared our model to that of Button et al. for evaluating viability from whole-community experimental results, similarly to previous reports (26, 28, 35) (Table 1). Our viability estimates (Vest) generally agreed with those using equation 1, but we have now provided 95% CI to depict the maximum and minimum viabilities that would match the returned positive-well distribution, as well as maximum and minimum values for the number of wells that ought to have contained a single cell. Maximum Vest ranged from 1.1% to >92.3% depending on the experiment, with a median Vest across all our experiments of 8.6% (Table 1). In one case, the extremely high value (FWC2) was better handled by our model than by equation 1, because it did not lead to a viability estimation greater than 100%. FWC and FWC2 represent Vest outliers compared with the entire data set (maximums of 59.7% and >92.3%, respectively) (Table 1). We believe these high numbers most likely resulted from underestimating the number of cells inoculated into each well (because of the use of microscopy, the presence of clumped cells, or possible pipette error [described in reference 35]), thus increasing the estimated viability.

Isolate-specific viability estimates.Our new model also facilitates taxon-specific viability estimates. Cultivation efficacy was evaluated for 71 cultured taxa matching ASVs within our detection limits (219 isolates) across 17 sites (1,207 pairwise combinations) by comparing the number of observed pure wells to those predicted by the Monte Carlo simulation using 9,999 bootstraps, 460 wells per experiment, and an assumption that all cells were viable (i.e., V = 100%). In total, for 1,158 out of 1,207 pairwise combinations (95.9%), the observed number of pure wells fell within the 95% CI of data simulated at matching relative abundance and inoculum size, suggesting that these two parameters alone could explain the observed cultivation success for most taxa (Table S1 [available at https://doi.org/10.6084/m9.figshare.12142113]). A total of 1,059 out of these 1,158 combinations (91%) recorded zero observed wells, but with a maximum relative abundance of 2.8% within these combinations, a score of zero fell within predicted 95% CI of simulations with 460 wells. Sensitivity analysis showed that with 460 wells per experiment, an observation of zero pure wells falls below the 95% CI’s lower bound (and is thus significantly depleted to enable viability to be estimated) for taxa with relative abundances of 2.3%, 2.9%, and 4.5% for inoculum sizes of one, two, and three cells per well, respectively (see Fig. S9 in the supplemental material). In fact, modeling DTE experiments from 92 wells to 9,200 wells per experiment showed that for taxa comprising ∼1% of a microbial community 1,104 wells (or 12 plates at 92 wells per plate), 1,380 wells (15 plates), and 2,576 wells (28 plates) were required to be statistically likely to recover at least one positive, pure well using inocula of one, two, or three cells per well, respectively, with V = 100% (Fig. S9).

A small but taxonomically relevant minority (49 out of 1,207) of pairwise combinations had a number of observed pure wells that fell outside the simulated 95% CI with V = 100% (Fig. 7). Of these, 28 had either one, two, or three more observed pure wells than the upper 95% CI (Table S1 [available at https://doi.org/10.6084/m9.figshare.12142113]), suggesting cultivability higher than expected based purely a model capturing the interaction between a Poisson-distributed inoculum and a binomially distributed relative abundance, with V = 100%. However, the deviance from the expected number of positive wells for those above the 95% CI was limited to three or fewer wells, meaning that we obtained only 1 to 3 more isolates than expected (Table S1). Conversely, those organisms that we isolated less frequently than expected showed greater deviance. 21 out of the 49 outliers had lower-than-expected cultivability (Fig. 7). These taxa had relative abundances ranging from 2.7% to 14.5% but recorded only 0, 1, or 2 isolates. In the most extreme case, ASV7629 (SAR11 LD12) at site ARD2c comprised 14.5% of the community but recorded no observed pure wells, compared to expected number of 13 to 30 isolates (95% CI) predicted by the Monte Carlo simulation. All the examples of taxa that were isolated less frequently than expected given the assumption of V = 100% belonged to either SAR11 LD12, SAR11 IIIa.1, or one particular OM43 ASV (7241) (Fig. 6).

FIG 7
  • Open in new tab
  • Download powerpoint
FIG 7

Actual versus expected numbers of isolates. Each point represents the actual number of isolates for every ASV/experiment pair compared to the expected number of isolates based on our model assuming 100% viability. Colors represent the relationship to the model predictions: green, isolates within the 95% CI for expected frequency; orange, actual isolates > maximum 95% CI for expected isolates; blue, actual isolates < minimum 95% CI for expected isolates. Circle size is proportional to the deviation of the number of actual isolates from the maximum (for orange) or minimum (for blue) 95% CI for expected isolates. The dotted line is the 1:1 ratio. Notable taxa on the extremities of the actual and expected values are labeled. All data points are provided in Table S1 in the supplemental material (available at https://doi.org/10.6084/m9.figshare.12142113).

We used our model to calculate estimated viability (Vest) for these organisms based on their cultivation frequency at sites where the assumption of V = 100% appeared to be violated (Table 3). Using the extreme example of SAR11 LD12 ASV7629 at site ARD2c, simulations across a range of V indicated that a result of zero positive wells fell within 95% of simulated values when the associated taxon Vest was ≤15%. When considering all anomalous cultivation results, LD12 had estimated maximum viabilities that ranged up to 55% (Table 3). OM43 (ASV7241) estimated maximum viabilities ranged from 52 to 80%, depending on the site, and similarly, SAR11 IIIa.1 ranged between 22 and 82% maximum viability (Table 3).

View this table:
  • View inline
  • View popup
  • Download powerpoint
TABLE 3

Estimated viabilities for taxa cultivated less frequently than expected

DISCUSSION

This work paired 17 DTE cultivation experiments with cultivation-independent assessments of microbial community structure in source waters to evaluate cultivation efficacy. We generated 328 new bacterial isolates representing 40 of the 777 OTUs and 71 of the 1,323 ASVs observed across all samples from which we inoculated DTE experiments. Stated another way, we successfully cultivated 5% of the total 3-year bacterioplankton community observed via either OTU or ASV analyses. A large fraction of our isolates (43% of cultured OTUs and 30% of cultured ASVs) represented taxa present at median relative abundances of >0.1%, with 15% and 4% of cultured OTUs and ASVs, respectively, at median abundances of >1%. One hundred forty of our isolates matched the top 50 most abundant OTUs, and 84 isolates matched the top 50 most abundant ASVs.

This campaign led to the first isolations of the abundant SAR11 LD12 and Actinobacteria acIV, the second isolate of the HIMB59 Alphaproteobacteria, and new genera within the Acetobacteraceae, Burkholderiaceae, OM241 and LSUCC0101-type Gammaproteobacteria, and MWH-UniPo Betaproteobacteria, thereby demonstrating again that continued DTE experimentation leads to isolation of previously uncultured organisms with value for aquatic microbiology. We have also added a considerable collection of isolates to previously cultured groups such as OM252 Gammaproteobacteria, BAL58 Betaproteobacteria, and HIMB11-type “Roseobacter” spp., and the majority of our isolates represent the first versions of these types of taxa from the Gulf of Mexico, adding comparative biogeographic value to these cultures.

Our viability model improved upon the statistical equation developed by Button and colleagues (33) to extend viability estimates to individual taxa within a mixed community and provide 95% CI to constrain those estimates. We cultured several groups of organisms abundant enough to evaluate viability with 460 wells (Fig. 7; see Fig. S9 in the supplemental material). The fact that these organisms were successfully cultured at least once meant that we could reasonably assume that the medium was sufficient for growth.

Some taxa were cultivated more frequently than expected (Fig. 7). We explore two possible explanations for this phenomenon: errors in quantification and variation in microbial cell organization. Any systematic error that led to underestimating the abundance of an organism would have correspondingly resulted in our underestimating the number of wells in which we would expect to find a pure culture of that organism. Such underestimations could come from primer biases associated with amplicon sequencing (76–78), but we do not know if those protocols specifically underestimate the OM252, MWH-UniPo, and HIMB11-type taxa cultured more frequently than expected (Fig. 7). However, due the low number of expected isolates in these groups and the small deviances in actual isolates from those expected numbers (within 1 to 3 isolates compared to expected values), the biases inherent in the relative abundance estimations for these taxa were probably small. Furthermore, one of the microorganisms isolated more frequently than expected matched the OM43 ASV1389 (see Fig. S6 in the supplemental material), whereas another OM43 ASV (7241) was cultivated less frequently than expected (see below), meaning that if primer bias were the cause of this discrepancy, it would have to be operating differently on very closely related organisms.

One possible biological explanation for why some isolates might have been cultured more frequently than expected is clumped cells. If cells of any given taxon in nature grew in small clusters, then the number of cells we added to a well would have been greater than expected based on a Poisson distribution. Furthermore, the model assumes that each cell is independent and that the composition of a subset of cells is only a function of the relative abundance of the taxon in the community. Within a cluster of cells, this assumption is violated, as the probability of cells being from the same taxon is higher. Thus, the model will underestimate the probability of a well being pure and therefore underestimate the number of pure wells likely to be observed within an experiment, leading to a greater number of isolates than expected. Future microscopy work could examine whether microorganisms such as OM252 and MWH-UniPo form small clusters in situ and/or in pure culture and whether this phenomenon may be different for different ASVs of OM43, or if clumping may be a transient phenotype.

We also identified three taxa—SAR11 LD12, SAR11 subclade IIIa.1, and the aforementioned OM43 ASV7241—that were isolated much less frequently than expected based on their abundances (Fig. 7; Table 3). This could mean that our assumption of V = 100% was incorrect or that, in contrast to the taxa that were cultured more frequently than expected (see above), our methods had biases that overestimated the abundance of these organisms, thereby overinflating the expected number of isolates. We used the modified 515/806RB primers, which have been shown to be much more accurate in quantifying SAR11 compared to fluorescent in situ hybridization (FISH) than the original 515/806 primers (within 6% ± 4% [standard deviation]), and this protocol almost always underestimates SAR11 abundance (76). This suggests that our expected number of isolates may have actually been underestimated and our cultivation success poorer than we measured, and therefore we may be overestimating viability for the SAR11 taxa in this study. Other sources of systematic error that might impinge on successful transfers, and thereby reduce our recovery, include sensitivity to pipette tip and/or flask material. However, the fact that these taxa were sometimes successfully isolated means that if these mechanisms were impacting successful transfers, then their activity was less than 100% efficient, which implies variations in subpopulation vulnerability that would be very similar conceptually to variations in subpopulation viability.

Another possible source of error that could have resulted in lower-than-expected numbers of isolates was the subset of experiments for which we did not transfer all positive wells due to limitations in available personnel time (Tables 1 and 3). However, our selection criteria for the subset of wells to transfer was based on flow cytometric signatures that would have encompassed small cells like SAR11 (see Results), and in any case, there were many examples of lower-than-expected recovery from other experiments where we transferred all positive wells (Table 3). Thus, we believe that these four experiments were unlikely to contribute major errors biasing our estimates of viability for SAR11 LD12, SAR11 IIIa.1, and other small cells like OM43.

If we instead explore biological reasons for the lower-than-expected numbers of positive wells in DTE experiments, a plausible explanation supported by the literature is simply that a large fraction of the population is in some state of inactivity or at least not actively dividing (79). Studies using uptake of a variety of radiolabeled carbon and sulfur sources have demonstrated that substantial fractions of SAR11 cells may be inactive, depending on the population (80–83). SAR11 cells in the northwest Atlantic and Mediterranean showed variable uptake of labeled leucine (30 to 50% [80, 81] and ∼25 to 55% [83, 84]) and amino acids (34 to 61% [80, 81] and 34 to 66% [80, 81]). Taken in reverse, this means that up to 75% of the SAR11 population may be dormant at any given time. In another study focused on brackish communities, fewer than 10% of SAR11 LD12 cells took up labeled leucine and/or thymidine (82). While this was likely not the ideal habitat for LD12 based on salinities above six (29, 82), this study supports the others described above that show substantial proportions of inactive SAR11 cells, the fraction of which may depend on environmental conditions and other unknown factors. Biorthogonal noncanonical amino acid tagging (BONCAT) showed a similar trend for SAR11 (85). These results also match general data indicating prevalent inactivity among aquatic bacterioplankton (79, 86–88). Although labeled uptake methods do not directly measure rates of cell division, the incorporation of these compounds requires active DNA replication or translation, which represent an even more fundamental level of activity than cell division (89).

Why might selection favor high percentages of subpopulation dormancy? One possibility is as an effective defense mechanism against abundant viruses. Viruses infecting SAR11 have been shown to be extremely abundant in both marine (90, 91) and freshwater (92) systems. Indeed, the paradox of high viral abundances and high host abundances in SAR11 has led to a refining of negative density-dependent selection through Lokta-Volterra predator-prey dynamics (93) to include heterogeneous susceptibility at the strain level (94, 95) and positive density-dependent selection through intraspecific proliferation of defense mechanisms (96). Activity of lytic viruses infecting SAR11 in situ demonstrated that phages infecting SAR11 have lower ratios of viral transcripts to host cells than in other abundant taxa and that observed abrupt changes in these ratios suggest coexistence of several SAR11 strains with different life strategies and phage susceptibilities (97). Phenotypic stochasticity of phage receptor expression has been shown to maintain a small proportion of phage-insensitive hosts within a population, enabling coexistence of predator and prey without extinction (98). Phages adsorb to a vast array of receptor proteins on their hosts, with many well-characterized receptors (e.g., OmpC, TonB, BtuB, and LamB) associated with nutrient uptake or osmoregulation (99). Selection therefore favors phenotypes that limit receptor expression, with an associated fitness cost, particularly in nutrient-limited environments.

However, an alternative mechanism is possible if a population of cells comprised a small number of susceptible cells and a large number of either resistant or dormant cells where presentation of receptor proteins is retained. The majority of host-virus encounters would occur with resistant or dormant cells and would constrain viral propagation through inefficient or failed infection, effectively acting as a sink for infectious particles. Prevalent lysogeny in SAR11 populations would provide a mechanism for establishing resistant cells via superinfection immunity (100, 101), where integration of a temperate phage prevents infection by other closely related viruses. There is growing evidence that many viruses infecting SAR11 are temperate (102, 103) and that reversion to virulence can be triggered through nutrient limitation (103), in contrast to other systems where lysogeny is favored in nutrient-poor conditions (104). Viral infection may also trigger host dormancy, lowering cellular metabolism to minimize energy requirements under nutrient-limited conditions (105). Such cells would be selected against during cultivation experiments, potentially explaining the rarity of SAR11 isolate genomes found to contain prophages. Dormancy and/or lysogeny would also enable long-term costability between abundant phages and their hosts (106) and resolve the apparent paradox of high host and virus abundances (101).

Detailed measurements of dormancy in SAR11 and what types of cellular functions become inactivated are part of our ongoing work. In the meantime, it is prudent to examine the implications of a substantial proportion of nondividing cells for our understanding of basic growth dynamics. Studies attempting to measure SAR11 growth rates in nature have yielded a wide range of results, ranging from 0.03 to 1.8 day−1 (70, 80, 83, 107–109). These span wider growth rates than observed for axenic cultures of SAR11 (0.4 to 1.2 day−1), but isolate-specific growth ranges within that spread are much more constrained (29, 36, 49, 110, 111). Conversion factors for determining production from [3H]leucine incorporation (112) are accurate for at least some Ia subclade members of SAR11 (113), so variations in growth rate estimates from microradiography experiments likely have other explanations. It is possible that different strains of SAR11 simply have variations in growth rate not captured by existing isolates. Another, not mutually exclusive, possibility is that the differences in in situ growth rate estimates also reflect variations in the proportion of actively dividing cells within the population. A simple model of cell division with binary fission where only a subset of cells divide and nondividing cells persist, rather than die, can still yield logarithmic growth curves (see Fig. S10 in the supplemental material) like those observed for SAR11 in pure culture (29, 49, 114). However, this subpopulation variability means that the division rate for the subset of cells that are actively dividing is much higher than calculated when assuming 100% dividing cells in the population. Based on our estimated viability for SAR11 LD12 of 15 to 55%, to obtain our previously calculated maximum division rate (0.5 day−1) for the whole culture (29), the per-cell division rate for only a subpopulation would span 2.48 to 0.79 day−1 (Fig. S10 and supplemental text). Verifying the proportion of SAR11 cells actively dividing in a culture may be challenging. Time-lapse microscopy (115) offers an elegant solution if SAR11 can be maintained for the requisite time periods for accurate measurements in a microfluidic device.

In addition to identifying taxa whose isolation success suggested deviations from biological assumptions of single planktonic cells with 100% viability, the model also revealed the limitations of DTE cultivation in assessing viability depending on relative abundance (Fig. S9). We cannot ascertain whether any given taxon may violate an assumption of V = 100% unless we have enough wells to demonstrate that it grew in fewer wells than expected. For example, taxa at 1% of the microbial community require more than 1,000 wells before the lack of a cultured organism represents a significant negative event, rather than a taxon simply lacking sufficient abundance to ensure inclusion in a well within 95% CI. In our 460-well experiments, we could not resolve whether taxa may have had viabilities below 100% if they were less than 3% of the community for any given experiment (Fig. S9). Modeling DTE experiments showed that for experiments targeting rare taxa, lower inoculum sizes are favored where a selective medium for enrichment is either unknown or undesirable. The exponential increase in the number of required wells with respect to the inoculum size is a function of a pure well requiring all cells within it to belong to the same taxon, assuming all cells are equally and optimally viable.

By providing taxon-specific predictions of viability from cultivation data, our model now facilitates an iterative process to improve experimental design and make cultivation more reliable. First, we use the cultivation success rates to determine for which taxa the assumption of 100% viability was violated. Second, we use the model to estimate viability for those organisms. Third, we use the viability and relative-abundance data to determine, within 95% CI, the appropriate number of inoculation attempts required to isolate a new version of that organism. Using SAR11 LD12 as an example, given a relative abundance of 10% and a viability of 15%, 800 DTE wells should yield four pure, positive wells (95% CI, 1 to 8). This means that, for microorganisms that we know successfully grow in our media, we can now statistically constrain the appropriate number of wells required to culture a given taxon again. For organisms that were not abundant enough to estimate viability using the model, we can use a conservative viability assumption (e.g., 50% [86]) with which to base our cultivation strategy, thereby still reducing uncertainty about the experimental effort necessary to reisolate one of these microorganisms.

Conclusions.This work has provided hundreds of new cultures for microbiological research, many among the most abundant members of the nGOM coastal bacterioplankton community. It also provides another demonstration of the effectiveness of sustained cultivation efforts for bringing previously uncultivated strains into culture. Our modeled cultivation results have generated compelling evidence for low viability within subpopulations of SAR11 LD12 and IIIa.1, as well as OM43 Betaproteobacteria. The prevalence of and controls on dormancy in these clades deserve further study. We anticipate that future work with larger DTE experiments will yield similar viability data about other groups of taxa with lower abundance, highlighting a valuable diagnostic application of DTE cultivation/modeling beyond the primary role in isolating new microorganisms. The integration of cultivation results, natural-abundance data from inoculum communities, and DTE modeling represents an important step forward in quantifying the risk associated with DTE efforts to isolate valuable taxa from new sources or repeating isolation from the same locations. We hope variations of this approach will be incorporated into wider community efforts to invest in culturing the uncultured.

MATERIALS AND METHODS

Sampling.Surface water samples were collected at six different sites once a year for 3 years, except for Terrebonne Bay, which was collected twice. The sites sampled were Lake Borgne (LKB) (Shell Beach, LA), Bay Pomme d’Or (JLB) (Buras, LA), Terrebonne Bay (TBON) (Cocodrie, LA), Atchafalaya River Delta (ARD) (Franklin, LA), Freshwater City (FWC) (Kaplan, LA), and Calcasieu Jetties (CJ) (Cameron, LA) (latitude and longitude coordinates are provided in Table S1 in the supplemental material [available at https://doi.org/10.6084/m9.figshare.12142113]). Water collection for biogeochemical and biological analysis followed the protocol described previously (35). Briefly, we collected surface water in a sterile, acid-washed polycarbonate bottle. Duplicate 120-ml water samples were filtered serially through 2.7-μm Whatman GF/D (GE Healthcare, Little Chalfont, UK) and 0.22-μm Sterivex (Millipore, MA, USA) filters and placed on ice until being transferred to −20°C in the laboratory (maximum of 3 h on ice). The University of Washington Marine Chemistry Laboratory analyzed duplicate subsamples of 50 ml 0.22-μm-filtered water collected in sterile 50-ml Falcon tubes (VWR, PA, USA) for concentrations of SiO4, PO43−, NH4+, NO3−, and NO2−. Samples for cell counts were filtered through a 2.7-μm GF/D filter, fixed with 10% formaldehyde, and stored on ice until enumeration (maximum of 3 h). Temperature, salinity, pH, and dissolved oxygen were measured using a handheld YSI 556 multiprobe system (YSI Inc., OH, USA). All metadata are available in Table S1.

DTE culturing and propagation.Isolation, propagation, and identification of isolates were completed as previously reported (29, 35, 116). A subsample of 2.7-μm-filtered surface water was stained with 1× SYBR green (Lonza, Basel, Switzerland) using a repeat pipettor and disposal tip (Gilson, WI, USA) and enumerated using a Guava EasyCyte 5HT HPL flow cytometer (Millipore, MA, USA) as described previously (116). After serial dilution to a predicted 1 to 3 cells · μl−1, 2 μl water was inoculated into five 2-ml 96-well polytetrafluoroethylene (PTFE) plates (Radleys, Essex, UK) containing 1.7 ml artificial seawater medium (Table S1 [available at https://doi.org/10.6084/m9.figshare.12142113]) using a 20-μl multichannel pipette (Gilson, WI, USA) to achieve an estimated 1 to 3 cells · well−1 (Table 1). The salinity of the medium was chosen to match in situ salinity after experiment JLB (January 2015) (Tables 1 and S1). After year two, a second generation of media, designated MWH, was designed to incorporate additional important osmolytes, reduced sulfur compounds, and other constituents (Tables 1 and S1) potentially necessary for in vitro growth of uncultivated clades (49, 117–123). The four corner wells of each plate were left uninoculated as negative controls for every experiment. Plates were covered using sterile, PTFE-coated silicon 96-well plate mats (Thermo Scientific, MA, USA). Cultures were incubated at in situ temperatures (Table S1) in the dark for 3 to 6 weeks and evaluated for positive growth (>104 cells · ml−1) by flow cytometry. Two hundred microliters from positive wells was transferred using a 200-μl single-channel pipette (Gilson, WI, USA) to duplicate 125-ml polycarbonate flasks (Corning, NY, USA) containing 50 ml of medium (29, 35, 116). At FWC, FWC2, JLB2c, and JLB3, not all positive wells were transferred because of the large number of positive wells. At each site, 48/301, 60/403, 60/103, and 60/146 of the positive wells were transferred, respectively, selected using flow cytometry signatures with <102 green fluorescence and <102 side scatter, which maximized our chances of isolating small microorganisms that encompass many of the most abundant and most wanted taxa, like SAR11, using our settings (116).

Culture identification.Cultures reaching ≥1 × 105 cells · ml−1 had 35 ml of the 50-ml volume filtered for identification via 16S rRNA gene PCR onto 25-mm, 0.22-μm polycarbonate filters (Millipore, MA, USA). DNA extractions were performed using the MoBio PowerWater DNA kit (Qiagen, MA, USA) following the manufacturer’s instructions and eluted in sterile water. The 16S rRNA gene was amplified as previously reported by Henson et al. (35) and sequenced at the Michigan State University Research Technology Support Facility Genomics Core. Evaluation of Sanger sequence quality was performed with 4Peaks (v. 1.7.1) (https://nucleobytes.com/4peaks/index.html), and forward and reverse complement sequences (converted via http://www.bioinformatics.org/sms/rev_comp.html) were assembled where overlap was sufficient using the CAP3 web server (http://doua.prabi.fr/software/cap3).

Community iTag sequencing, OTUs, and single-nucleotide variants.Sequentially filtered (2.7 μm, 0.22 μm) duplicate samples were extracted and analyzed using our previously reported protocols and settings (35, 124). We sequenced the 2.7- to 0.22-μm fraction for this study because this fraction corresponded with the <2.7-μm communities that were used for the DTE experiments. To avoid batch sequencing effects, DNA from the first seven collections reported previously (35) was resequenced with the additional samples from this study (FWC2 and after) (Table 1). We targeted the 16S rRNA gene V4 region with the 515F-806RB primer set (which corrects for poor amplification of taxa like SAR11) (76, 77) using Illumina MiSeq 2× 250-bp paired-end sequencing at Argonne National Laboratories, resulting in 2,343,106 raw reads for the 2.7- to 0.22-μm fraction. Using mothur v1.33.3 (125), we clustered 16S rRNA gene amplicons into distinctive operational taxonomic units (OTUs) with a 0.03 dissimilarity threshold (OTU0.03) and classified them according to the Silva v119 database (126, 127). After these steps, 55,256 distinct OTUs0.03 remained. We also used minimum entropy decomposition (MED) to partition reads into fine-scale amplicon single-nucleotide variants (ASVs) (68). Reads were first analyzed using mothur as described above up to the screen.seqs() command. The cleaned-reads fasta file was converted to MED-compatible headers with the “mothur2oligo” tool renamer.pl from the functions in MicrobeMiseq (https://github.com/DenefLab/MicrobeMiseq) (128) using the fasta output from screen.seqs() and the mothur group file. These curated reads were analyzed using MED (v. 2.1) with the flags –M 60, and -d 1. MED resulted in 2,813 refined ASVs. ASVs were classified in mothur using classify.seqs(), the Silva v119 database, and a cutoff bootstrap value of 80% (129). After classification, we removed ASVs identified as “chloroplast,” “mitochondria,” or “unknown” from the data set.

Community analyses.OTU (OTU0.03) and ASV abundances were analyzed within the R statistical environment v.3.2.1 (130) following previously published protocols (29, 35, 124). Using the package PhyloSeq (131), OTUs and ASVs were rarefied using the command rarefy_even_depth(), and OTUs/ASVs without at least two reads in four of the 34 samples (2 sites; ∼11%) were removed. This cutoff was used to remove potentially spurious OTUs/ASVs resulting from sequencing errors. Our modified PhyloSeq script is available on our GitHub repository, https://github.com/thrash-lab/Modified-Phyloseq. After filtering, the data sets contained 777 unique OTUs and 1,323 unique ASVs (Table S1 [available at https://doi.org/10.6084/m9.figshare.12142113]). For site-specific community comparisons, beta diversity between sites was examined using Bray-Curtis distances via ordination with nonmetric multidimensional scaling (NMDS) (Table S1). The nutrient data were normalized using the R function “scale,” which subtracts the mean and divides by the standard deviation for each column. The influence of the transformed environmental parameters on beta diversity was calculated in R with the envfit function. Relative abundances of an OTU or ASV from each sample were calculated as the number of reads over the sum of all the reads in that sample. The relative abundance was then averaged between biological duplicates for a given OTU or ASV. To determine the best-matching OTU or ASV for a given LSUCC isolate, the OTU representative fasta file, provided by mothur using get.oturep(), and the ASV fasta file were used to create a BLAST database (makeblastdb) against which the LSUCC isolate 16S rRNA genes could be searched via blastn (BLAST v 2.2.26) (“OTU_ASVrep_db”; available as supplemental information at https://doi.org/10.6084/m9.figshare.12142098). We designated an LSUCC isolate 16S rRNA gene match with an OTU or ASV sequence based on ≥97% or ≥99% sequence identity, respectively, as well as a ≥247-bp alignment.

16S rRNA gene phylogeny.Taxa in the Alpha-, Beta-, and Gammaproteobacteria phylogenies from reference 35 served as the backbone for the trees in the current work. For places in these trees with poor representation near isolate sequences, additional taxa were selected by searching the 16S rRNA genes of LSUCC isolates against the NCBI nucleotide database online with BLASTn (132) and selecting a variable number of best hits. The Bacteroidetes and Actinobacteria trees were composed entirely of nonredundant top 100 to 300 MegaBLAST hits to a local version of the NCBI nucleotide database, accessed August 2018. Sequences were aligned with MUSCLE v3.6 (133) using default settings and culled with TrimAl v1.4.rev22 (134) using the -automated1 flag, and the final alignment was inferred with IQ-TREE v1.6.11 (135) with default settings and -bb 1000 for ultrafast bootstrapping (136). Tips were edited with the nw_rename script within Newick Utilities v1.6 (137), and trees were visualized with Archaeopteryx (138). Fasta files for these trees and the naming keys are available as supplemental information at https://doi.org/10.6084/m9.figshare.12142098.

Assessment of isolate novelty.We quantified taxonomic novelty using BLASTn of our isolate 16S rRNA genes to those of other known isolates collected in three databases: (i) The NCBI nucleotide database (accessed August 2018) (NCBIdb), (ii) a custom database comprised of sequences from DTE experiments in other labs (DTEdb), and (iii) a database containing all of our isolate 16S rRNA genes (LSUCCdb). The DTEdb and LSUCCdb fasta files are available as supplemental information at https://doi.org/10.6084/m9.figshare.12142098. We compared our isolate sequences to these databases as follows.

(i) All representative sequences were searched against the nucleotide database using BLASTn (BLAST+v. 2.7.1) with the flags -perc_identity 84, -evalue 1E-6, -task blastn, -outfmt “6 qseqid sseqid pident length slen qlen mismatch evalue bitscore sscinames sblastnames stitle,” and -negative_gilist to remove uncultured and environmental sequences. The negative GI list was obtained by searching “environmental samples[organism] OR metagenomes[orgn]” in the NCBI nucleotide database (accessed 12 September 12 2018), and hits were downloaded in GI list format. This negative GI list is available as supplemental information at https://doi.org/10.6084/m9.figshare.12142098. The resultant hits from the NCBIdb search were further manually curated to remove sequences classified as single-cell genomes, clones, duplicates, and previously deposited LSUCC isolates.

(ii) We observed that many known HTCC, IMCC, and HIMB isolates that had previously been described as matching our clades (see Fig. S1 to S5 in the supplemental material) were missing from the resultant lists of nucleotide hits, so we extracted isolate accession numbers from numerous DTE experiments (26–28, 31, 34, 37, 44, 139, 140) from the nucleotide database via blastdbcmd and generated a separate DTEdb using makeblastdb. Duplicate accession numbers found in the NCBIdb were removed. The same BLASTn settings as for the first step were used to search our isolate sequences against DTEdb. Any match that fell below the lowest percent identity hit to the NCBIdb was removed from the DTEdb search, since the match would not have been present in the first NCBIdb search.

(iii) Finally, using the same BLASTn settings, we compared all pairwise identities of our 328 LSUCC isolate 16S rRNA gene sequences via the LSUCCdb.

The output from these searches is available in Table S1 (https://doi.org/10.6084/m9.figshare.12142113) under the “taxonomic novelty” tab.

We placed our LSUCC isolates into 55 taxonomic groups based on sharing ≥94% identity and/or their occurrence in monophyletic groups within our 16S rRNA gene trees (Fig. S1 to S5) (see above). For visualization purposes, in groups with multiple isolates we used our chronologically first cultivated isolate as the representative sequence for blastn searches, and these are the top point (100% identity to itself) in each group column of Fig. 1. Sequences from the other DTE culture collections were labeled with the corresponding collection name, while all other hits were labeled “other.”

Geographic novelty was assessed by manually screening the accession numbers from hits to LSUCC isolates with ≥99% 16S rRNA gene sequence identity for the latitude and longitude from a connected publication or location name (e.g., source, country, or site) in the NCBI description. LSUCC isolates in the Janibacter sp., Micrococcus sp., Altererythrobacter sp., Pseudomonas sp., and Phycicoccus sp. groups (16 total isolates) were not assessed because of missing isolation source information and no traceable publication. Isolation locations were plotted for a subset of important taxa (Table S1 [available at https://doi.org/10.6084/m9.figshare.12142113], “Map_cultivars” tab) using the “LSUCC_cultivar_map.R” available at our GitHub repository, https://github.com/thrash-lab/Cultivar-novelty-map.

Modeling DTE cultivation via Monte Carlo simulations.We developed a model using Monte Carlo simulation to estimate the median number of positive and pure wells (and associated 95% confidence intervals [CI]) expected from a DTE experiment for a given taxon at different inoculum sizes (λ), relative abundances (r), and viability (V) (Fig. 5). For each bootstrap, the number of cells added to each well was simulated using a Poisson distribution at a mean inoculum size of λ cells per well across n wells. The number of cells added to each well that belonged to a specific taxon was then estimated using a binomial distribution where the number of trials was set as the number of cells in a well and the probability of a cell belonging to a specific taxon, r, was the relative abundance of its representative ASV in the community analysis. Wells that contained at least one cell of a specific taxon were designated “positive.” Wells in which all the cells belonged to a specific taxon were designated “pure.” Finally, the influence of taxon-specific viability on recovery of pure wells was simulated using a second binomial distribution, where the number of cells within a pure well was used as the number of trials and the probability of growth was a viability score ranging from 0 to 1. For each simulation, 9,999 bootstraps were performed. Code for the model and all simulations is available at our GitHub repository, https://github.com/thrash-lab/viability_test.

Actual versus expected number of isolates.For each taxon in each DTE experiment, the Monte Carlo simulation was used to evaluate whether the number of recovered pure wells for each taxon was within 95% CI of simulated estimates, assuming optimum growth conditions (i.e., V = 100%). For each of 9,999 bootstraps, 460 wells were simulated with the inoculum size used for the experiment and the relative abundance of the ASV. For taxa where the number of expected wells fell outside the 95% CI of the model, a deviance score was calculated as the difference between the actual number of wells observed and median of the simulated data set. The results of this output are presented in Table S1 (available at https://doi.org/10.6084/m9.figshare.12142113) under the “Expected vs actual” tab, and the R script for visualizing this output (Fig. 7) is available at our GitHub repository, https://github.com/thrash-lab/EvsA-visualization.

Estimating viability in underrepresented taxa.For taxa where the observed number of positive wells was lower than the 95% CI lower limit within a given experiment, and because our analysis was restricted to only those organisms for which our medium was sufficient for growth at least once, the deviance was assumed to be a function of a viability term, V (ranging from 0 to 1), associated with suboptimal growth conditions, dormancy, persister cells, etc. To estimate a value of viability for a given taxon within a particular experiment, the Monte Carlo simulation was run using an experiment-appropriate inoculum size, relative abundance, and number of wells (460 for each experiment). Taxon-specific viability was tested across a range of decreasing values from 99% to 1% until such time as the observed number of pure wells for a given taxon fell between the 95% CI bounds of the simulated data. At this point, the viability value is the maximum viability of the taxon that enables the observed number of pure wells for a given taxon to be explained by the model. The results of this output are presented in Table 3 and Table S1 (available at https://doi.org/10.6084/m9.figshare.12142113) under the “Expected vs actual” tab.

Likelihood of recovering taxa at different relative abundances.To estimate the number of wells required in a DTE experiment to have a significant chance of recovering a taxon with a relative abundance of r, assuming optimum growth conditions (V = 100%), the Monte Carlo model was used to simulate experiments from 92 wells to 9,200 wells per experiment across a range of relative abundances from 0% to 100% in 0.1% increments and a range of inoculum sizes (cells per well of 1, 1.27, 1.96, 2, 3, 4, and 5). Each experiment was bootstrapped 999 times, and the number of bootstraps in which the lower bound of the 95% CI was ≥1 was recorded.

Data availability.All iTag sequences are available at the Sequence Read Archive with accession numbers SRR6235382 to SRR6235415 (29). PCR-generated 16S rRNA gene sequences from this study are accessible in NCBI GenBank under the accession numbers MK603525 to MK603769. Previously generated 16S rRNA gene sequences are accessible in NCBI GenBank under the accession numbers KU382357 to KU382438 (35). Table S1 is available at https://doi.org/10.6084/m9.figshare.12142113.

ACKNOWLEDGMENTS

We thank Nancy Rabalais for her comments and edits on a very early draft of the manuscript. Portions of this research were conducted with high-performance computing resources provided by Louisiana State University (http://www.hpc.lsu.edu) and the University of Southern California (https://hpcc.usc.edu/).

This work was supported by the Department of Biological Sciences at Louisiana State University, a Louisiana Board of Regents grant to J.C.T. [LEQSF(2014-17)-RD-A-06], a National Academies of Science, Engineering, and Medicine Gulf Research Program Early Career Fellowship to J.C.T., and the Dornsife College of Letters, Arts and Sciences at the University of Southern California. B.T. was partially funded by NERC award NE/R010935/1 and by the Simons Foundation BIOS-SCOPE program.

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

M.W.H. led sample collections, cultivation experiments, nucleic acid extraction, amplicon sequencing, and analyses; V.C.L., D.M.P., J.L.W., and A.M.L. supported sample collections, cultivation experiments, and nucleic acid extractions; M.W.H., and J.C.T. conducted cultivation comparisons and phylogenetic analyses; B.T. developed the viability model and led the statistical analyses; C.C. derived the cell-specific growth rate equations incorporating viability; J.C.T. designed the study and assisted with sample collections and model refinement; M.W.H., B.T., and J.C.T. led manuscript writing; and all authors contributed to and reviewed the text.

FOOTNOTES

    • Received 21 April 2020.
    • Accepted 13 June 2020.
    • Accepted manuscript posted online 19 June 2020.
  • Supplemental material is available online only.

  • Copyright © 2020 Henson et al.

This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

REFERENCES

  1. 1.↵
    1. Metcalf WW,
    2. Griffin BM,
    3. Cicchillo RM,
    4. Gao J,
    5. Janga SC,
    6. Cooke HA,
    7. Circello BT,
    8. Evans BS,
    9. Martens-Habbena W,
    10. Stahl DA,
    11. van der Donk WA
    . 2012. Synthesis of methylphosphonic acid by marine microbes: a source for methane in the aerobic ocean. Science 337:1104–1107. doi:10.1126/science.1219875.
    OpenUrlAbstract/FREE Full Text
  2. 2.↵
    1. Carini P,
    2. White AE,
    3. Campbell EO,
    4. Giovannoni SJ
    . 2014. Methane production by phosphate-starved SAR11 chemoheterotrophic marine bacteria. Nat Commun 5:4346. doi:10.1038/ncomms5346.
    OpenUrlCrossRefPubMed
  3. 3.↵
    1. Karl DM,
    2. Beversdorf L,
    3. Björkman KM,
    4. Church MJ,
    5. Martinez A,
    6. Delong EF
    . 2008. Aerobic production of methane in the sea. Nat Geosci 1:473–478. doi:10.1038/ngeo234.
    OpenUrlCrossRef
  4. 4.↵
    1. Steindler L,
    2. Schwalbach MS,
    3. Smith DP,
    4. Chan F,
    5. Giovannoni SJ
    . 2011. Energy starved Candidatus Pelagibacter ubique substitutes light-mediated ATP production for endogenous carbon respiration. PLoS One 6:e19725. doi:10.1371/journal.pone.0019725.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Gómez-Consarnau L,
    2. González JM,
    3. Coll-Lladó M,
    4. Gourdon P,
    5. Pascher T,
    6. Neutze R,
    7. Pedrós-Alió C,
    8. Pinhassi J
    . 2007. Light stimulates growth of proteorhodopsin-containing marine Flavobacteria. Nature 445:210–213. doi:10.1038/nature05381.
    OpenUrlCrossRefPubMedWeb of Science
  6. 6.↵
    1. Daims H,
    2. Lebedeva EV,
    3. Pjevac P,
    4. Han P,
    5. Herbold C,
    6. Albertsen M,
    7. Jehmlich N,
    8. Palatinszky M,
    9. Vierheilig J,
    10. Bulaev A,
    11. Kirkegaard RH,
    12. von Bergen M,
    13. Rattei T,
    14. Bendinger B,
    15. Nielsen PH,
    16. Wagner M
    . 2015. Complete nitrification by Nitrospira bacteria. Nature 528:504–509. doi:10.1038/nature16461.
    OpenUrlCrossRef
  7. 7.↵
    1. Lam KS
    . 2006. Discovery of novel metabolites from marine actinomycetes. Curr Opin Microbiol 9:245–251. doi:10.1016/j.mib.2006.03.004.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Jensen PR,
    2. Williams PG,
    3. Oh D-C,
    4. Zeigler L,
    5. Fenical W
    . 2007. Species-specific secondary metabolite production in marine actinomycetes of the genus Salinispora. Appl Environ Microbiol 73:1146–1152. doi:10.1128/AEM.01891-06.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    1. Nett M,
    2. Erol Ö,
    3. Kehraus S,
    4. Köck M,
    5. Krick A,
    6. Eguereva E,
    7. Neu E,
    8. König GM
    . 2006. Siphonazole, an unusual metabolite from Herpetosiphon sp. Angew Chem Int Ed Engl 45:3863–3867. doi:10.1002/anie.200504525.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Ling LL,
    2. Schneider T,
    3. Peoples AJ,
    4. Spoering AL,
    5. Engels I,
    6. Conlon BP,
    7. Mueller A,
    8. Schäberle TF,
    9. Hughes DE,
    10. Epstein S,
    11. Jones M,
    12. Lazarides L,
    13. Steadman VA,
    14. Cohen DR,
    15. Felix CR,
    16. Fetterman KA,
    17. Millett WP,
    18. Nitti AG,
    19. Zullo AM,
    20. Chen C,
    21. Lewis K
    . 2015. A new antibiotic kills pathogens without detectable resistance. Nature 517:455–459. doi:10.1038/nature14098.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    1. Lloyd KG,
    2. Steen AD,
    3. Ladau J,
    4. Yin J,
    5. Crosby L
    . 2018. Phylogenetically novel uncultured microbial cells dominate earth microbiomes. mSystems 3:e00055-18. doi:10.1128/mSystems.00055-18.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Hug LA
    . 2018. Sizing up the uncultured microbial majority. mSystems 3:e00185-18. doi:10.1128/mSystems.00185-18.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    1. Steen AD,
    2. Crits-Christoph A,
    3. Carini P,
    4. DeAngelis KM,
    5. Fierer N,
    6. Lloyd KG,
    7. Thrash JC
    . 2019. High proportions of bacteria and archaea across most biomes remain uncultured. ISME J 13:3126–3130. doi:10.1038/s41396-019-0484-y.
    OpenUrlCrossRef
  14. 14.↵
    1. Overmann J,
    2. Abt B,
    3. Sikorski J
    . 2017. Present and future of culturing bacteria. Annu Rev Microbiol 71:711–730. doi:10.1146/annurev-micro-090816-093449.
    OpenUrlCrossRef
  15. 15.↵
    1. Rappé MS
    . 2013. Stabilizing the foundation of the house that ’omics builds: the evolving value of cultured isolates to marine microbiology. Curr Opin Microbiol 16:618–624. doi:10.1016/j.mib.2013.09.009.
    OpenUrlCrossRefPubMed
  16. 16.↵
    1. Carini P
    . 2019. A “cultural” renaissance: genomics breathes new life into an old craft. mSystems 4:e00092-19. doi:10.1128/mSystems.00092-19.
    OpenUrlAbstract/FREE Full Text
  17. 17.↵
    1. Staley JT,
    2. Konopka A
    . 1985. Measurement of in situ activities of nonphotosynthetic microorganisms in aquatic and terrestrial habitats. Annu Rev Microbiol 39:321–346. doi:10.1146/annurev.mi.39.100185.001541.
    OpenUrlCrossRefPubMedWeb of Science
  18. 18.↵
    1. Amann RI,
    2. Ludwig W,
    3. Schleifer KH
    . 1995. Phylogenetic identification and in situ detection of individual microbial cells without cultivation. Microbiol Rev 59:143–169. doi:10.1128/MMBR.59.1.143-169.1995.
    OpenUrlAbstract/FREE Full Text
  19. 19.↵
    1. Hahn MW,
    2. Koll U,
    3. Schmidt J
    . 2019. Isolation and cultivation of bacteria, p 313–351. In Hurst CJ (ed), The structure and function of aquatic microbial communities. Springer International Publishing, Cham, Switzerland.
  20. 20.↵
    1. Kaeberlein T,
    2. Lewis K,
    3. Epstein SS
    . 2002. Isolating “uncultivable” microorganisms in pure culture in a simulated natural environment. Science 296:1127–1129. doi:10.1126/science.1070633.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    1. Zengler K,
    2. Toledo G,
    3. Rappe M,
    4. Elkins J,
    5. Mathur EJ,
    6. Short JM,
    7. Keller M
    . 2002. Cultivating the uncultured. Proc Natl Acad Sci U S A 99:15681–15686. doi:10.1073/pnas.252630999.
    OpenUrlAbstract/FREE Full Text
  22. 22.↵
    1. Steinert G,
    2. Whitfield S,
    3. Taylor MW,
    4. Thoms C,
    5. Schupp PJ
    . 2014. Application of diffusion growth chambers for the cultivation of marine sponge-associated bacteria. Mar Biotechnol (NY) 16:594–603. doi:10.1007/s10126-014-9575-y.
    OpenUrlCrossRef
  23. 23.↵
    1. Nichols D,
    2. Cahoon N,
    3. Trakhtenberg EM,
    4. Pham L,
    5. Mehta A,
    6. Belanger A,
    7. Kanigan T,
    8. Lewis K,
    9. Epstein SS
    . 2010. Use of Ichip for high-throughput in situ cultivation of “uncultivable” microbial Species. Appl Environ Microbiol 76:2445–2450. doi:10.1128/AEM.01754-09.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Tandogan N,
    2. Abadian PN,
    3. Epstein S,
    4. Aoi Y,
    5. Goluch ED
    . 2014. Isolation of microorganisms using sub-micrometer constrictions. PLoS One 9:e101429. doi:10.1371/journal.pone.0101429.
    OpenUrlCrossRef
  25. 25.↵
    1. Hahn MW,
    2. Stadler P,
    3. Wu QL,
    4. Pöckl M
    . 2004. The filtration–acclimatization method for isolation of an important fraction of the not readily cultivable bacteria. J Microbiol Methods 57:379–390. doi:10.1016/j.mimet.2004.02.004.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    1. Connon SA,
    2. Giovannoni SJ
    . 2002. High-throughput methods for culturing microorganisms in very-low-nutrient media yield diverse new marine isolates. Appl Environ Microbiol 68:3878–3885. doi:10.1128/aem.68.8.3878-3885.2002.
    OpenUrlAbstract/FREE Full Text
  27. 27.↵
    1. Yang SJ,
    2. Kang I,
    3. Cho JC
    . 2016. Expansion of cultured bacterial diversity by large-scale dilution-to-extinction culturing from a single seawater sample. Microb Ecol 71:29–43. doi:10.1007/s00248-015-0695-3.
    OpenUrlCrossRefPubMed
  28. 28.↵
    1. Stingl U,
    2. Tripp HJ,
    3. Giovannoni SJ
    . 2007. Improvements of high-throughput culturing yielded novel SAR11 strains and other abundant marine bacteria from the Oregon coast and the Bermuda Atlantic Time Series study site. ISME J 1:361–371. doi:10.1038/ismej.2007.49.
    OpenUrlCrossRefPubMedWeb of Science
  29. 29.↵
    1. Henson MW,
    2. Lanclos VC,
    3. Faircloth BC,
    4. Thrash JC
    . 2018. Cultivation and genomics of the first freshwater SAR11 (LD12) isolate. ISME J 12:1846–1860. doi:10.1038/s41396-018-0092-2.
    OpenUrlCrossRef
  30. 30.↵
    1. Hahnke RL,
    2. Bennke CM,
    3. Fuchs BM,
    4. Mann AJ,
    5. Rhiel E,
    6. Teeling H,
    7. Amann R,
    8. Harder J
    . 2015. Dilution cultivation of marine heterotrophic bacteria abundant after a spring phytoplankton bloom in the North Sea. Environ Microbiol 17:3515–3526. doi:10.1111/1462-2920.12479.
    OpenUrlCrossRef
  31. 31.↵
    1. Sosa OA,
    2. Gifford SM,
    3. Repeta DJ,
    4. DeLong EF
    . 2015. High molecular weight dissolved organic matter enrichment selects for methylotrophs in dilution to extinction cultures. ISME J 9:2725–2739. doi:10.1038/ismej.2015.68.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Schut F,
    2. de Vries EJ,
    3. Gottschal JC,
    4. Robertson BR,
    5. Harder W,
    6. Prins RA,
    7. Button DK
    . 1993. Isolation of typical marine bacteria by dilution culture: growth, maintenance, and characteristics of isolates under laboratory conditions. Appl Environ Microbiol 59:2150–2160. doi:10.1128/AEM.59.7.2150-2160.1993.
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    1. Button DK,
    2. Schut F,
    3. Quang P,
    4. Martin R,
    5. Robertson BR
    . 1993. Viability and isolation of marine bacteria by dilution culture: theory, procedures, and initial results. Appl Environ Microbiol 59:881–891. doi:10.1128/AEM.59.3.881-891.1993.
    OpenUrlAbstract/FREE Full Text
  34. 34.↵
    1. Song J,
    2. Oh HM,
    3. Cho JC
    . 2009. Improved culturability of SAR11 strains in dilution-to-extinction culturing from the East Sea, West Pacific Ocean. FEMS Microbiol Lett 295:141–147. doi:10.1111/j.1574-6968.2009.01623.x.
    OpenUrlCrossRefPubMedWeb of Science
  35. 35.↵
    1. Henson MW,
    2. Pitre DM,
    3. Weckhorst JL,
    4. Lanclos VC,
    5. Webber AT,
    6. Thrash JC
    . 2016. Artificial seawater media facilitate cultivating members of the microbial majority from the Gulf of Mexico. mSphere 1:e00028-16. doi:10.1128/mSphere.00124-16.
    OpenUrlAbstract/FREE Full Text
  36. 36.↵
    1. Rappé MS,
    2. Connon SA,
    3. Vergin KL,
    4. Giovannoni SJ
    . 2002. Cultivation of the ubiquitous SAR11 marine bacterioplankton clade. Nature 418:630–633. doi:10.1038/nature00917.
    OpenUrlCrossRefPubMedWeb of Science
  37. 37.↵
    1. Marshall KT,
    2. Morris RM
    . 2013. Isolation of an aerobic sulfur oxidizer from the SUP05/Arctic96BD-19 clade. ISME J 7:452–455. doi:10.1038/ismej.2012.78.
    OpenUrlCrossRefPubMedWeb of Science
  38. 38.↵
    1. Shah V,
    2. Chang BX,
    3. Morris RM
    . 2017. Cultivation of a chemoautotroph from the SUP05 clade of marine bacteria that produces nitrite and consumes ammonium. ISME J 11:263–271. doi:10.1038/ismej.2016.87.
    OpenUrlCrossRef
  39. 39.↵
    1. Spietz RL,
    2. Lundeen RA,
    3. Zhao X,
    4. Nicastro D,
    5. Ingalls AE,
    6. Morris RM
    . 2019. Heterotrophic carbon metabolism and energy acquisition in Candidatus Thioglobus singularis strain PS1, a member of the SUP05 clade of marine Gammaproteobacteria. Environ Microbiol 21:2391–2401. doi:10.1111/1462-2920.14623.
    OpenUrlCrossRef
  40. 40.↵
    1. Huggett MJ,
    2. Hayakawa DH,
    3. Rappé MS
    . 2012. Genome sequence of strain HIMB624, a cultured representative from the OM43 clade of marine Betaproteobacteria. Stand Genomic Sci 6:11–20. doi:10.4056/sigs.2305090.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Giovannoni SJ,
    2. Hayakawa DH,
    3. Tripp HJ,
    4. Stingl U,
    5. Givan SA,
    6. Cho J-C,
    7. Oh H-M,
    8. Kitner JB,
    9. Vergin KL,
    10. Rappé MS
    . 2008. The small genome of an abundant coastal ocean methylotroph. Environ Microbiol 10:1771–1782. doi:10.1111/j.1462-2920.2008.01598.x.
    OpenUrlCrossRefPubMedWeb of Science
  42. 42.↵
    1. Durham BP,
    2. Grote J,
    3. Whittaker KA,
    4. Bender SJ,
    5. Luo H,
    6. Grim SL,
    7. Brown JM,
    8. Casey JR,
    9. Dron A,
    10. Florez-Leiva L,
    11. Krupke A,
    12. Luria CM,
    13. Mine AH,
    14. Nigro OD,
    15. Pather S,
    16. Talarmin A,
    17. Wear EK,
    18. Weber TS,
    19. Wilson JM,
    20. Church MJ,
    21. DeLong EF,
    22. Karl DM,
    23. Steward GF,
    24. Eppley JM,
    25. Kyrpides NC,
    26. Schuster S,
    27. Rappé MS
    . 2014. Draft genome sequence of marine alphaproteobacterial strain HIMB11, the first cultivated representative of a unique lineage within the Roseobacter clade possessing an unusually small genome. Stand Genomic Sci 9:632–645. doi:10.4056/sigs.4998989.
    OpenUrlCrossRefPubMed
  43. 43.↵
    1. Cho JC,
    2. Giovannoni SJ
    . 2004. Cultivation and growth characteristics of a diverse group of oligotrophic marine gammaproteobacteria. Appl Environ Microbiol 70:432–440. doi:10.1128/AEM.70.1.432-440.2004.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Kim S,
    2. Kang I,
    3. Seo J-H,
    4. Cho J-C
    . 2019. Culturing the ubiquitous freshwater actinobacterial acI lineage by supplying a biochemical “helper”catalase. ISME J 13:2252–2263. doi:10.1038/s41396-019-0432-x.
    OpenUrlCrossRef
  45. 45.↵
    1. Tripp HJ
    . 2013. The unique metabolism of SAR11 aquatic bacteria. J Microbiol 51:147–153. doi:10.1007/s12275-013-2671-2.
    OpenUrlCrossRef
  46. 46.↵
    1. Nichols D,
    2. Lewis K,
    3. Orjala J,
    4. Mo S,
    5. Ortenberg R,
    6. O'Connor P,
    7. Zhao C,
    8. Vouros P,
    9. Kaeberlein T,
    10. Epstein SS
    . 2008. Short peptide induces an “uncultivable” microorganism to grow in vitro. Appl Environ Microbiol 74:4889–4897. doi:10.1128/AEM.00393-08.
    OpenUrlAbstract/FREE Full Text
  47. 47.↵
    1. Stewart EJ
    . 2012. Growing unculturable bacteria. J Bacteriol 194:4151–4160. doi:10.1128/JB.00345-12.
    OpenUrlAbstract/FREE Full Text
  48. 48.↵
    1. Kell DB,
    2. Young M
    . 2000. Bacterial dormancy and culturability: the role of autocrine growth factors. Curr Opin Microbiol 3:238–243. doi:10.1016/s1369-5274(00)00082-5.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    1. Carini P,
    2. Steindler L,
    3. Beszteri S,
    4. Giovannoni SJ
    . 2013. Nutrient requirements for growth of the extreme oligotroph “Candidatus Pelagibacter ubique” HTCC1062 on a defined medium. ISME J 7:592–602. doi:10.1038/ismej.2012.122.
    OpenUrlCrossRefPubMedWeb of Science
  50. 50.↵
    1. Epstein SS
    . 2013. The phenomenon of microbial uncultivability. Curr Opin Microbiol 16:636–642. doi:10.1016/j.mib.2013.08.003.
    OpenUrlCrossRefPubMed
  51. 51.↵
    1. Shah D,
    2. Zhang Z,
    3. Khodursky A,
    4. Kaldalu N,
    5. Kurg K,
    6. Lewis K
    . 2006. Persisters: a distinct physiological state of E. coli. BMC Microbiol 6:53. doi:10.1186/1471-2180-6-53.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Kell D,
    2. Potgieter M,
    3. Pretorius E
    . 2015. Individuality, phenotypic differentiation, dormancy and “persistence” in culturable bacterial systems: commonalities shared by environmental, laboratory, and clinical microbiology. F1000Res 4:179. doi:10.12688/f1000research.6709.2.
    OpenUrlCrossRef
  53. 53.↵
    1. Lennon JT,
    2. Jones SE
    . 2011. Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat Rev Microbiol 9:119–130. doi:10.1038/nrmicro2504.
    OpenUrlCrossRefPubMed
  54. 54.↵
    1. Grassi L,
    2. Di Luca M,
    3. Maisetta G,
    4. Rinaldi AC,
    5. Esin S,
    6. Trampuz A,
    7. Batoni G
    . 2017. Generation of persister cells of Pseudomonas aeruginosa and Staphylococcus aureus by chemical treatment and evaluation of their susceptibility to membrane-targeting agents. Front Microbiol 8:1917. doi:10.3389/fmicb.2017.01917.
    OpenUrlCrossRef
  55. 55.↵
    1. Bergkessel M,
    2. Basta DW,
    3. Newman DK
    . 2016. The physiology of growth arrest: uniting molecular and environmental microbiology. Nat Rev Microbiol 14:549–562. doi:10.1038/nrmicro.2016.107.
    OpenUrlCrossRef
  56. 56.↵
    1. Kussell E,
    2. Leibler S
    . 2005. Phenotypic diversity, population growth, and information in fluctuating environments. Science 309:2075–2078. doi:10.1126/science.1114383.
    OpenUrlAbstract/FREE Full Text
  57. 57.↵
    1. Kurm V,
    2. van der Putten WH,
    3. Hol WHG
    . 2019. Cultivation-success of rare soil bacteria is not influenced by incubation time and growth medium. PLoS One 14:e0210073. doi:10.1371/journal.pone.0210073.
    OpenUrlCrossRef
  58. 58.↵
    1. D'Onofrio A,
    2. Crawford JM,
    3. Stewart EJ,
    4. Witt K,
    5. Gavrish E,
    6. Epstein S,
    7. Clardy J,
    8. Lewis K
    . 2010. Siderophores from neighboring organisms promote the growth of uncultured bacteria. Chem Biol 17:254–264. doi:10.1016/j.chembiol.2010.02.010.
    OpenUrlCrossRefPubMed
  59. 59.↵
    1. Thrash JC
    . 2019. Culturing the uncultured: risk versus reward. mSystems 4:e00130-19. doi:10.1128/mSystems.00130-19.
    OpenUrlAbstract/FREE Full Text
  60. 60.↵
    1. Noell SE,
    2. Giovannoni SJ
    . 2019. SAR11 bacteria have a high affinity and multifunctional glycine betaine transporter. Environ Microbiol 21:2559–2575. doi:10.1111/1462-2920.14649.
    OpenUrlCrossRef
  61. 61.↵
    1. Thume K,
    2. Gebser B,
    3. Chen L,
    4. Meyer N,
    5. Kieber DJ,
    6. Pohnert G
    . 2018. The metabolite dimethylsulfoxonium propionate extends the marine organosulfur cycle. Nature 563:412–415. doi:10.1038/s41586-018-0675-0.
    OpenUrlCrossRef
  62. 62.↵
    1. Asher EC,
    2. Dacey JWH,
    3. Stukel M,
    4. Long MC,
    5. Tortell PD
    . 2017. Processes driving seasonal variability in DMS, DMSP, and DMSO concentrations and turnover in coastal Antarctic waters. Limnol Oceanogr 62:104–124. doi:10.1002/lno.10379.
    OpenUrlCrossRef
  63. 63.↵
    1. Mizuno CM,
    2. Rodriguez-Valera F,
    3. Ghai R
    . 2015. Genomes of planktonic Acidimicrobiales: widening horizons for marine Actinobacteria by metagenomics. mBio 6:e02083-14. doi:10.1128/mBio.02083-14.
    OpenUrlAbstract/FREE Full Text
  64. 64.↵
    1. Lee J,
    2. Kwon KK,
    3. Lim S-I,
    4. Song J,
    5. Choi AR,
    6. Yang S-H,
    7. Jung K-H,
    8. Lee J-H,
    9. Kang SG,
    10. Oh H-M,
    11. Cho JC
    . 2019. Isolation, cultivation, and genome analysis of proteorhodopsin-containing SAR116-clade strain Candidatus Puniceispirillum marinum IMCC1322. J Microbiol 57:676–687. doi:10.1007/s12275-019-9001-2.
    OpenUrlCrossRef
  65. 65.↵
    1. Kitzinger K,
    2. Padilla CC,
    3. Marchant HK,
    4. Hach PF,
    5. Herbold CW,
    6. Kidane AT,
    7. Könneke M,
    8. Littmann S,
    9. Mooshammer M,
    10. Niggemann J,
    11. Petrov S,
    12. Richter A,
    13. Stewart FJ,
    14. Wagner M,
    15. Kuypers MMM,
    16. Bristow LA
    . 2019. Cyanate and urea are substrates for nitrification by Thaumarchaeota in the marine environment. Nat Microbiol 4:234–243. doi:10.1038/s41564-018-0316-2.
    OpenUrlCrossRef
  66. 66.↵
    1. Kamennaya NA,
    2. Post AF
    . 2011. Characterization of cyanate metabolism in marine Synechococcus and Prochlorococcus spp. Appl Environ Microbiol 77:291–301. doi:10.1128/AEM.01272-10.
    OpenUrlAbstract/FREE Full Text
  67. 67.↵
    1. Eren AM,
    2. Maignien L,
    3. Sul WJ,
    4. Murphy LG,
    5. Grim SL,
    6. Morrison HG,
    7. Sogin ML
    . 2013. Oligotyping: differentiating between closely related microbial taxa using 16S rRNA gene data. Methods Ecol Evol 4:1111–1119. doi:10.1111/2041-210X.12114.
    OpenUrlCrossRef
  68. 68.↵
    1. Eren AM,
    2. Borisy GG,
    3. Huse SM,
    4. Mark Welch JL
    . 2014. Oligotyping analysis of the human oral microbiome. Proc Natl Acad Sci U S A 111:E2875–E2884. doi:10.1073/pnas.1409644111.
    OpenUrlAbstract/FREE Full Text
  69. 69.↵
    1. Needham DM,
    2. Fuhrman JA
    . 2016. Pronounced daily succession of phytoplankton, archaea and bacteria following a spring bloom. Nat Microbiol 1:16005. doi:10.1038/nmicrobiol.2016.5.
    OpenUrlCrossRef
  70. 70.↵
    1. Campbell BJ,
    2. Yu L,
    3. Heidelberg JF,
    4. Kirchman DL
    . 2011. Activity of abundant and rare bacteria in a coastal ocean. Proc Natl Acad Sci U S A 108:12776–12781. doi:10.1073/pnas.1101405108.
    OpenUrlAbstract/FREE Full Text
  71. 71.↵
    1. Oh HM,
    2. Kwon KK,
    3. Kang I,
    4. Kang SG,
    5. Lee JH,
    6. Kim SJ,
    7. Cho JC
    . 2010. Complete genome sequence of “Candidatus puniceispirillum marinum” IMCC1322, a representative of the SAR116 clade in the Alphaproteobacteria. J Bacteriol 192:3240–3241. doi:10.1128/JB.00347-10.
    OpenUrlAbstract/FREE Full Text
  72. 72.↵
    1. Fegatella F,
    2. Lim J,
    3. Kjelleberg S,
    4. Cavicchioli R
    . 1998. Implications of rRNA operon copy number and ribosome content in the marine oligotrophic ultramicrobacterium Sphingomonas sp. strain RB2256. Appl Environ Microbiol 64:4433–4438. doi:10.1128/AEM.64.11.4433-4438.1998.
    OpenUrlAbstract/FREE Full Text
  73. 73.↵
    1. Acinas SG,
    2. Marcelino LA,
    3. Klepac-Ceraj V,
    4. Polz MF
    . 2004. Divergence and redundancy of 16S rRNA sequences in genomes with multiple rrn operons. J Bacteriol 186:2629–2635. doi:10.1128/JB.186.9.2629-2635.2004.
    OpenUrlAbstract/FREE Full Text
  74. 74.↵
    1. Pedrós-Alió C
    . 2012. The rare bacterial biosphere. Annu Rev Mar Sci 4:449–466. doi:10.1146/annurev-marine-120710-100948.
    OpenUrlCrossRefPubMed
  75. 75.↵
    1. Martiny AC
    . 2019. High proportions of bacteria are culturable across major biomes. ISME J 13:2125–2128. doi:10.1038/s41396-019-0410-3.
    OpenUrlCrossRef
  76. 76.↵
    1. Apprill A,
    2. McNally S,
    3. Parsons R,
    4. Weber L
    . 2015. Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton. Aquat Microb Ecol 75:129–137. doi:10.3354/ame01753.
    OpenUrlCrossRef
  77. 77.↵
    1. Walters W,
    2. Hyde ER,
    3. Berg-Lyons D,
    4. Ackermann G,
    5. Humphrey G,
    6. Parada A,
    7. Gilbert JA,
    8. Jansson JK,
    9. Caporaso JG,
    10. Fuhrman JA,
    11. Apprill A,
    12. Knight R
    . 2016. Improved bacterial 16S rRNA gene (V4 and V4-5) and fungal internal transcribed spacer marker gene primers for microbial community surveys. mSystems 1:e00009-15. doi:10.1128/mSystems.00009-15.
    OpenUrlCrossRef
  78. 78.↵
    1. Parada AE,
    2. Needham DM,
    3. Fuhrman JA
    . 2016. Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ Microbiol 18:1403–1414. doi:10.1111/1462-2920.13023.
    OpenUrlCrossRefPubMed
  79. 79.↵
    1. del Giorgio PA,
    2. Gasol JM
    . 2008. Physiological structure and single-cell activity in marine bacterioplankton, p 243–298. In Kirchman DL (ed), Microbial ecology of the oceans, 2nd ed. Wiley, New York, NY.
  80. 80.↵
    1. Malmstrom RR,
    2. Kiene RP,
    3. Cottrell MT,
    4. Kirchman DL
    . 2004. Contribution of SAR11 bacteria to dissolved dimethylsulfoniopropionate and amino acid uptake in the North Atlantic Ocean. Appl Environ Microbiol 70:4129–4135. doi:10.1128/AEM.70.7.4129-4135.2004.
    OpenUrlAbstract/FREE Full Text
  81. 81.↵
    1. Malmstrom RR,
    2. Cottrell MT,
    3. Elifantz H,
    4. Kirchman DL
    . 2005. Biomass production and assimilation of dissolved organic matter by SAR11 bacteria in the Northwest Atlantic Ocean. Appl Environ Microbiol 71:2979–2986. doi:10.1128/AEM.71.6.2979-2986.2005.
    OpenUrlAbstract/FREE Full Text
  82. 82.↵
    1. Piwosz K,
    2. Salcher MM,
    3. Zeder M,
    4. Ameryk A,
    5. Pernthaler J
    . 2013. Seasonal dynamics and activity of typical freshwater bacteria in brackish waters of the Gulf of Gdańsk. Limnol Oceanogr 58:817–826. doi:10.4319/lo.2013.58.3.0817.
    OpenUrlCrossRef
  83. 83.↵
    1. Laghdass M,
    2. Catala P,
    3. Caparros J,
    4. Oriol L,
    5. Lebaron P,
    6. Obernosterer I
    . 2012. High contribution of SAR11 to microbial activity in the north west Mediterranean Sea. Microb Ecol 63:324–333. doi:10.1007/s00248-011-9915-7.
    OpenUrlCrossRefPubMed
  84. 84.↵
    1. Salter I,
    2. Galand PE,
    3. Fagervold SK,
    4. Lebaron P,
    5. Obernosterer I,
    6. Oliver MJ,
    7. Suzuki MT,
    8. Tricoire C
    . 2015. Seasonal dynamics of active SAR11 ecotypes in the oligotrophic Northwest Mediterranean Sea. ISME J 9:347–360. doi:10.1038/ismej.2014.129.
    OpenUrlCrossRefPubMed
  85. 85.↵
    1. Leizeaga A,
    2. Estrany M,
    3. Forn I,
    4. Sebastián M
    . 2017. Using click-chemistry for visualizing in situ changes of translational activity in planktonic marine bacteria. Front Microbiol 8:2360. doi:10.3389/fmicb.2017.02360.
    OpenUrlCrossRef
  86. 86.↵
    1. Kirchman DL
    . 2016. Growth rates of microbes in the oceans. Annu Rev Mar Sci 8:285–309. doi:10.1146/annurev-marine-122414-033938.
    OpenUrlCrossRef
  87. 87.↵
    1. Samo TJ,
    2. Smriga S,
    3. Malfatti F,
    4. Sherwood BP,
    5. Azam F
    . 2014. Broad distribution and high proportion of protein synthesis active marine bacteria revealed by click chemistry at the single cell level. Front Mar Sci 1:48. doi:10.3389/fmars.2014.00060.
    OpenUrlCrossRef
  88. 88.↵
    1. Smriga S,
    2. Samo TJ,
    3. Malfatti F,
    4. Villareal J,
    5. Azam F
    . 2014. Individual cell DNA synthesis within natural marine bacterial assemblages as detected by “click”chemistry. Aquat Microb Ecol 72:269–280. doi:10.3354/ame01698.
    OpenUrlCrossRef
  89. 89.↵
    1. Bradley JA,
    2. Amend JP,
    3. LaRowe DE
    . 2019. Survival of the fewest: microbial dormancy and maintenance in marine sediments through deep time. Geobiology 17:43–59. doi:10.1111/gbi.12313.
    OpenUrlCrossRef
  90. 90.↵
    1. Zhao Y,
    2. Temperton B,
    3. Thrash JC,
    4. Schwalbach MS,
    5. Vergin KL,
    6. Landry ZC,
    7. Ellisman M,
    8. Deerinck T,
    9. Sullivan MB,
    10. Giovannoni SJ
    . 2013. Abundant SAR11 viruses in the ocean. Nature 494:357–360. doi:10.1038/nature11921.
    OpenUrlCrossRefPubMedWeb of Science
  91. 91.↵
    1. Martinez-Hernandez F,
    2. Fornas Ò,
    3. Lluesma Gomez M,
    4. Garcia-Heredia I,
    5. Maestre-Carballa L,
    6. López-Pérez M,
    7. Haro-Moreno JM,
    8. Rodriguez-Valera F,
    9. Martinez-Garcia M
    . 2019. Single-cell genomics uncover Pelagibacter as the putative host of the extremely abundant uncultured 37-F6 viral population in the ocean. ISME J 13:232–236. doi:10.1038/s41396-018-0278-7.
    OpenUrlCrossRef
  92. 92.↵
    1. Chen L-X,
    2. Zhao Y,
    3. McMahon KD,
    4. Mori JF,
    5. Jessen GL,
    6. Nelson TC,
    7. Warren LA,
    8. Banfield JF
    . 2019. Wide distribution of phage that infect freshwater SAR11 bacteria. mSystems 4:e00410-19. doi:10.1128/mSystems.00410-19.
    OpenUrlAbstract/FREE Full Text
  93. 93.↵
    1. Thingstad TF,
    2. Lignell R
    . 1997. Theoretical models for the control of bacterial growth rate, abundance, diversity and carbon demand. Aquat Microb Ecol 13:19–27. doi:10.3354/ame013019.
    OpenUrlCrossRefWeb of Science
  94. 94.↵
    1. Våge S,
    2. Storesund JE,
    3. Thingstad TF
    . 2013. SAR11 viruses and defensive host strains. Nature 499:E3–E4. doi:10.1038/nature12387.
    OpenUrlCrossRefPubMedWeb of Science
  95. 95.↵
    1. Thingstad TF,
    2. Vage S,
    3. Storesund JE,
    4. Sandaa RA,
    5. Giske J
    . 2014. A theoretical analysis of how strain-specific viruses can control microbial species diversity. Proc Natl Acad Sci U S A 111:7813–7818. doi:10.1073/pnas.1400909111.
    OpenUrlAbstract/FREE Full Text
  96. 96.↵
    1. Giovannoni S,
    2. Temperton B,
    3. Zhao Y
    . 2013. Giovannoni et al. reply. Nature 499:E4–E5. doi:10.1038/nature12388.
    OpenUrlCrossRef
  97. 97.↵
    1. Alonso-Sáez L,
    2. Morán XAG,
    3. Clokie MR
    . 2018. Low activity of lytic pelagiphages in coastal marine waters. ISME J 12:2100–2102. doi:10.1038/s41396-018-0185-y.
    OpenUrlCrossRef
  98. 98.↵
    1. Chapman-McQuiston E,
    2. Wu XL
    . 2008. Stochastic receptor expression allows sensitive bacteria to evade phage attack. I. Experiments. Biophys J 94:4525–4536. doi:10.1529/biophysj.107.120212.
    OpenUrlCrossRefPubMedWeb of Science
  99. 99.↵
    1. Bertozzi Silva J,
    2. Storms Z,
    3. Sauvageau D
    . 2016. Host receptors for bacteriophage adsorption. FEMS Microbiol Lett 363:fnw002. doi:10.1093/femsle/fnw002.
    OpenUrlCrossRefPubMed
  100. 100.↵
    1. Williamson SJ,
    2. Houchin LA,
    3. McDaniel L,
    4. Paul JH
    . 2002. Seasonal variation in lysogeny as depicted by prophage induction in Tampa Bay, Florida. Appl Environ Microbiol 68:4307–4314. doi:10.1128/aem.68.9.4307-4314.2002.
    OpenUrlAbstract/FREE Full Text
  101. 101.↵
    1. Stough JMA,
    2. Tang X,
    3. Krausfeldt LE,
    4. Steffen MM,
    5. Gao G,
    6. Boyer GL,
    7. Wilhelm SW
    . 2017. Molecular prediction of lytic vs lysogenic states for Microcystis phage: metatranscriptomic evidence of lysogeny during large bloom events. PLoS One 12:e0184146. doi:10.1371/journal.pone.0184146.
    OpenUrlCrossRef
  102. 102.↵
    1. Zhao Y,
    2. Qin F,
    3. Zhang R,
    4. Giovannoni SJ,
    5. Zhang Z,
    6. Sun J,
    7. Du S,
    8. Rensing C
    . 2019. Pelagiphages in the Podoviridae family integrate into host genomes. Environ Microbiol 21:1989–2001. doi:10.1111/1462-2920.14487.
    OpenUrlCrossRef
  103. 103.↵
    1. Morris RM,
    2. Cain KR,
    3. Hvorecny KL,
    4. Kollman JM
    . 2020. Lysogenic host–virus interactions in SAR11 marine bacteria. Nat Microbiol doi:10.1038/s41564-020-0725-x.
    OpenUrlCrossRef
  104. 104.↵
    1. Warwick-Dugdale J,
    2. Buchholz HH,
    3. Allen MJ,
    4. Temperton B
    . 2019. Host-hijacking and planktonic piracy: how phages command the microbial high seas. Virol J 16:15. doi:10.1186/s12985-019-1120-1.
    OpenUrlCrossRef
  105. 105.↵
    1. Yu Z-C,
    2. Chen X-L,
    3. Shen Q-T,
    4. Zhao D-L,
    5. Tang B-L,
    6. Su H-N,
    7. Wu Z-Y,
    8. Qin Q-L,
    9. Xie B-B,
    10. Zhang X-Y,
    11. Yu Y,
    12. Zhou B-C,
    13. Chen B,
    14. Zhang Y-Z
    . 2015. Filamentous phages prevalent in Pseudoalteromonas spp. confer properties advantageous to host survival in Arctic sea ice. ISME J 9:871–881. doi:10.1038/ismej.2014.185.
    OpenUrlCrossRefPubMed
  106. 106.↵
    1. Shkoporov AN,
    2. Khokhlova EV,
    3. Fitzgerald CB,
    4. Stockdale SR,
    5. Draper LA,
    6. Ross RP,
    7. Hill C
    . 2018. ΦCrAss001 represents the most abundant bacteriophage family in the human gut and infects Bacteroides intestinalis. Nat Commun 9:4781. doi:10.1038/s41467-018-07225-7.
    OpenUrlCrossRef
  107. 107.↵
    1. Teira E,
    2. Martínez-García S,
    3. Lønborg C,
    4. Alvarez-Salgado XA
    . 2009. Growth rates of different phylogenetic bacterioplankton groups in a coastal upwelling system. Environ Microbiol Rep 1:545–554. doi:10.1111/j.1758-2229.2009.00079.x.
    OpenUrlCrossRefPubMedWeb of Science
  108. 108.↵
    1. Lankiewicz TS,
    2. Cottrell MT,
    3. Kirchman DL
    . 2016. Growth rates and rRNA content of four marine bacteria in pure cultures and in the Delaware estuary. ISME J 10:823–832. doi:10.1038/ismej.2015.156.
    OpenUrlCrossRefPubMed
  109. 109.↵
    1. Ferrera I,
    2. Gasol JM,
    3. Sebastián M,
    4. Hojerová E,
    5. Koblízek M
    . 2011. Comparison of growth rates of aerobic anoxygenic phototrophic bacteria and other bacterioplankton groups in coastal Mediterranean waters. Appl Environ Microbiol 77:7451–7458. doi:10.1128/AEM.00208-11.
    OpenUrlAbstract/FREE Full Text
  110. 110.↵
    1. Carini P,
    2. Campbell EO,
    3. Morré J,
    4. Sañudo-Wilhelmy SA,
    5. Thrash JC,
    6. Bennett SE,
    7. Temperton B,
    8. Begley T,
    9. Giovannoni SJ
    . 2014. Discovery of a SAR11 growth requirement for thiamin’s pyrimidine precursor and its distribution in the Sargasso Sea. ISME J 8:1727–1738. doi:10.1038/ismej.2014.61.
    OpenUrlCrossRefPubMedWeb of Science
  111. 111.↵
    1. Grant SR,
    2. Church MJ,
    3. Ferrón S,
    4. Laws EA,
    5. Rappé MS
    . 2019. Elemental composition, phosphorous uptake, and characteristics of growth of a SAR11 strain in batch and continuous culture mSystems 4:e00218-18. doi:10.1128/mSystems.00218-18.
    OpenUrlAbstract/FREE Full Text
  112. 112.↵
    1. Simon M,
    2. Azam F
    . 1989. Protein content and protein synthesis rates of planktonic marine bacteria. Mar Ecol Prog Ser 51:201–213. doi:10.3354/meps051201.
    OpenUrlCrossRefWeb of Science
  113. 113.↵
    1. White AE,
    2. Giovannoni SJ,
    3. Zhao Y,
    4. Vergin K,
    5. Carlson CA
    . 2019. Elemental content and stoichiometry of SAR11 chemoheterotrophic marine bacteria. Limnol Oceanogr 4:44–51. doi:10.1002/lol2.10103.
    OpenUrlCrossRef
  114. 114.↵
    1. Smith DP,
    2. Cameron Thrash J,
    3. Nicora CD,
    4. Lipton MS,
    5. Burnum-Johnson KE,
    6. Carini P,
    7. Smith RD,
    8. Giovannoni SJ
    . 2013. Proteomic and transcriptomic analyses of “Candidatus Pelagibacter ubique” describe the first PII-independent response to nitrogen limitation in a free-living alphaproteobacterium. mBio 4:e00133-12. doi:10.1128/mBio.00133-12.
    OpenUrlAbstract/FREE Full Text
  115. 115.↵
    1. van Vliet S,
    2. Dal Co A,
    3. Winkler AR,
    4. Spriewald S,
    5. Stecher B,
    6. Ackermann M
    . 2018. Spatially correlated gene expression in bacterial groups: the role of lineage history, spatial gradients, and cell-cell interactions. Cell Syst 6:496–507. doi:10.1016/j.cels.2018.03.009.
    OpenUrlCrossRef
  116. 116.↵
    1. Thrash JC,
    2. Weckhorst JL,
    3. Pitre DM
    . 2017. Cultivating fastidious microbes, p 57–78. In McGenity TJ, Timmis KN, Nogales B (ed), Hydrocarbon and lipid microbiology protocols. Springer, Berlin, Germany.
  117. 117.↵
    1. Vila-Costa M,
    2. Simó R,
    3. Harada H,
    4. Gasol JM,
    5. Slezak D,
    6. Kiene RP
    . 2006. Dimethylsulfoniopropionate uptake by marine phytoplankton. Science 314:652–654. doi:10.1126/science.1131043.
    OpenUrlAbstract/FREE Full Text
  118. 118.↵
    1. Mou X,
    2. Hodson RE,
    3. Moran MA
    . 2007. Bacterioplankton assemblages transforming dissolved organic compounds in coastal seawater. Environ Microbiol 9:2025–2037. doi:10.1111/j.1462-2920.2007.01318.x.
    OpenUrlCrossRefPubMedWeb of Science
  119. 119.↵
    1. Stein LY
    . 2015. Microbiology: cyanate fuels the nitrogen cycle. Nature 524:43–44. doi:10.1038/nature14639.
    OpenUrlCrossRef
  120. 120.↵
    1. Repeta DJ,
    2. Ferrón S,
    3. Sosa OA,
    4. Johnson CG,
    5. Repeta LD,
    6. Acker M,
    7. Delong EF,
    8. Karl DM
    . 2016. Marine methane paradox explained by bacterial degradation of dissolved organic matter. Nat Geosci 9:884–887. doi:10.1038/ngeo2837.
    OpenUrlCrossRef
  121. 121.↵
    1. Curson ARJ,
    2. Liu J,
    3. Bermejo Martínez A,
    4. Green RT,
    5. Chan Y,
    6. Carrión O,
    7. Williams BT,
    8. Zhang S-H,
    9. Yang G-P,
    10. Bulman Page PC,
    11. Zhang X-H,
    12. Todd JD
    . 2017. Dimethylsulfoniopropionate biosynthesis in marine bacteria and identification of the key gene in this process. Nat Microbiol 2:17009. doi:10.1038/nmicrobiol.2017.9.
    OpenUrlCrossRef
  122. 122.↵
    1. Widner B,
    2. Mulholland MR
    . 2017. Cyanate distribution and uptake in North Atlantic coastal waters. Limnol Oceanogr 62:2538–2549. doi:10.1002/lno.10588.
    OpenUrlCrossRef
  123. 123.↵
    1. Widner B,
    2. Mulholland MR,
    3. Mopper K
    . 2016. Distribution, sources, and sinks of cyanate in the coastal North Atlantic Ocean. Environ Sci Technol Lett 3:297–302. doi:10.1021/acs.estlett.6b00165.
    OpenUrlCrossRef
  124. 124.↵
    1. Henson MW,
    2. Hanssen J,
    3. Spooner G,
    4. Fleming P,
    5. Pukonen M,
    6. Stahr F,
    7. Thrash JC
    . 2018. Nutrient dynamics and stream order influence microbial community patterns along a 2914 kilometer transect of the Mississippi River. Limnol Oceanogr 63:1837–1855. doi:10.1002/lno.10811.
    OpenUrlCrossRef
  125. 125.↵
    1. Schloss PD,
    2. Westcott SL,
    3. Ryabin T,
    4. Hall JR,
    5. Hartmann M,
    6. Hollister EB,
    7. Lesniewski RA,
    8. Oakley BB,
    9. Parks DH,
    10. Robinson CJ,
    11. Sahl JW,
    12. Stres B,
    13. Thallinger GG,
    14. Van Horn DJ,
    15. Weber CF
    . 2009. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 75:7537–7541. doi:10.1128/AEM.01541-09.
    OpenUrlAbstract/FREE Full Text
  126. 126.↵
    1. Pruesse E,
    2. Quast C,
    3. Knittel K,
    4. Fuchs BM,
    5. Ludwig W,
    6. Peplies J,
    7. Glockner FO,
    8. Glöckner FO
    . 2007. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res 35:7188–7196. doi:10.1093/nar/gkm864.
    OpenUrlCrossRefPubMedWeb of Science
  127. 127.↵
    1. Quast C,
    2. Pruesse E,
    3. Yilmaz P,
    4. Gerken J,
    5. Schweer T,
    6. Yarza P,
    7. Peplies J,
    8. Glöckner FO
    . 2013. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res 41:590–596. doi:10.1093/nar/gks1219.
    OpenUrlCrossRef
  128. 128.↵
    1. Fujimoto M,
    2. Cavaletto J,
    3. Liebig JR,
    4. McCarthy A,
    5. Vanderploeg HA,
    6. Denef VJ
    . 2016. Spatiotemporal distribution of bacterioplankton functional groups along a freshwater estuary to pelagic gradient in Lake Michigan. J Great Lakes Res 42:1036–1048. doi:10.1016/j.jglr.2016.07.029.
    OpenUrlCrossRef
  129. 129.↵
    1. Cole JR,
    2. Wang Q,
    3. Fish JA,
    4. Chai B,
    5. McGarrell DM,
    6. Sun Y,
    7. Brown CT,
    8. Porras-Alfaro A,
    9. Kuske CR,
    10. Tiedje JM
    . 2014. Ribosomal Database Project: data and tools for high throughput rRNA analysis. Nucleic Acids Res 42:D633–D642. doi:10.1093/nar/gkt1244.
    OpenUrlCrossRefPubMedWeb of Science
  130. 130.↵
    R Core Team. 2017. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org.
  131. 131.↵
    1. McMurdie PJ,
    2. Holmes S
    . 2013. phyloseq: an R Package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi:10.1371/journal.pone.0061217.
    OpenUrlCrossRefPubMed
  132. 132.↵
    1. Camacho C,
    2. Coulouris G,
    3. Avagyan V,
    4. Ma N,
    5. Papadopoulos J,
    6. Bealer K,
    7. Madden TL
    . 2009. BLAST+: architecture and applications. BMC Bioinformatics 10:421. doi:10.1186/1471-2105-10-421.
    OpenUrlCrossRefPubMed
  133. 133.↵
    1. Edgar RC
    . 2004. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32:1792–1797. doi:10.1093/nar/gkh340.
    OpenUrlCrossRefPubMedWeb of Science
  134. 134.↵
    1. Capella-Gutiérrez S,
    2. Silla-Martínez JM,
    3. Gabaldón T
    . 2009. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Bioinformatics 25:1972–1973. doi:10.1093/bioinformatics/btp348.
    OpenUrlCrossRefPubMedWeb of Science
  135. 135.↵
    1. Nguyen L-T,
    2. Schmidt HA,
    3. von Haeseler A,
    4. Minh BQ
    . 2015. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Mol Biol Evol 32:268–274. doi:10.1093/molbev/msu300.
    OpenUrlCrossRefPubMed
  136. 136.↵
    1. Hoang DT,
    2. Chernomor O,
    3. von Haeseler A,
    4. Minh BQ,
    5. Vinh LS
    . 2018. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol 35:518–522. doi:10.1093/molbev/msx281.
    OpenUrlCrossRefPubMed
  137. 137.↵
    1. Junier T,
    2. Zdobnov EM
    . 2010. The Newick utilities: high-throughput phylogenetic tree processing in the UNIX shell. Bioinformatics 26:1669–1670. doi:10.1093/bioinformatics/btq243.
    OpenUrlCrossRefPubMedWeb of Science
  138. 138.↵
    1. Han MV,
    2. Zmasek CM
    . 2009. phyloXML: XML for evolutionary biology and comparative genomics. BMC Bioinformatics 10:356. doi:10.1186/1471-2105-10-356.
    OpenUrlCrossRefPubMed
  139. 139.↵
    1. Stingl U,
    2. Cho JC,
    3. Foo W,
    4. Vergin KL,
    5. Lanoil B,
    6. Giovannoni SJ
    . 2008. Dilution-to-extinction culturing of psychrotolerant planktonic bacteria from permanently ice-covered lakes in the McMurdo Dry Valleys, Antarctica. Microb Ecol 55:395–405. doi:10.1007/s00248-007-9284-4.
    OpenUrlCrossRefPubMedWeb of Science
  140. 140.↵
    1. Page KA,
    2. Connon SA,
    3. Giovannoni SJ
    . 2004. Representative freshwater bacterioplankton isolated from Crater Lake, Oregon. Appl Environ Microbiol 70:6542–6550. doi:10.1128/AEM.70.11.6542-6550.2004.
    OpenUrlAbstract/FREE Full Text
PreviousNext
Back to top
Download PDF
Citation Tools
Expanding the Diversity of Bacterioplankton Isolates and Modeling Isolation Efficacy with Large-Scale Dilution-to-Extinction Cultivation
Michael W. Henson, V. Celeste Lanclos, David M. Pitre, Jessica Lee Weckhorst, Anna M. Lucchesi, Chuankai Cheng, Ben Temperton, J. Cameron Thrash
Applied and Environmental Microbiology Aug 2020, 86 (17) e00943-20; DOI: 10.1128/AEM.00943-20

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this Applied and Environmental Microbiology article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Expanding the Diversity of Bacterioplankton Isolates and Modeling Isolation Efficacy with Large-Scale Dilution-to-Extinction Cultivation
(Your Name) has forwarded a page to you from Applied and Environmental Microbiology
(Your Name) thought you would be interested in this article in Applied and Environmental Microbiology.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Expanding the Diversity of Bacterioplankton Isolates and Modeling Isolation Efficacy with Large-Scale Dilution-to-Extinction Cultivation
Michael W. Henson, V. Celeste Lanclos, David M. Pitre, Jessica Lee Weckhorst, Anna M. Lucchesi, Chuankai Cheng, Ben Temperton, J. Cameron Thrash
Applied and Environmental Microbiology Aug 2020, 86 (17) e00943-20; DOI: 10.1128/AEM.00943-20
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

dilution to extinction
cultivation
bacterioplankton
LSUCC
microbial ecology
coastal microbiology

Related Articles

Cited By...

About

  • About AEM
  • Editor in Chief
  • Editorial Board
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Article Types
  • Ethics
  • Contact Us

Follow #AppEnvMicro

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

 

Print ISSN: 0099-2240; Online ISSN: 1098-5336