- Open Access
Metatranscriptomic and comparative genomic insights into resuscitation mechanisms during enrichment culturing
© The Author(s). 2018
- Received: 14 June 2018
- Accepted: 4 December 2018
- Published: 26 December 2018
The pure culture of prokaryotes remains essential to elucidating the role of these organisms. Scientists have reasoned that hard to cultivate microorganisms might grow in pure culture if provided with the chemical components of their natural environment. However, most microbial species in the biosphere that would otherwise be “culturable” may fail to grow because of their growth state in nature, such as dormancy. That means even if scientist would provide microorganisms with the natural environment, such dormant microorganisms probably still remain in a dormant state.
We constructed an enrichment culture system for high-efficiency isolation of uncultured strains from marine sediment. Degree of enrichment analysis, dormant and active taxa calculation, viable but non-culturable bacteria resuscitation analysis, combined with metatranscriptomic and comparative genomic analyses of the interactions between microbial communications during enrichment culture showed that the so-called enrichment method could culture the “uncultured” not only through enriching the abundance of “uncultured,” but also through the resuscitation mechanism. In addition, the enrichment culture was a complicated mixed culture system, which contains the competition, cooperation, or coordination among bacterial communities, compared with pure cultures.
Considering that cultivation techniques must evolve further—from axenic to mixed cultures—for us to fully understand the microbial world, we should redevelop an understanding of the classic enrichment culture method. Enrichment culture methods can be developed and used to construct a model for analyzing mixed cultures and exploring microbial dark matter. This study provides a new train of thought to mining marine microbial dark matter based on mixed cultures.
- Enrichment culture
- Mixed culture
- Uncultured bacteria
The ocean microbiome was one of the first microbiomes to be studied, but the diversity and distribution of its members are only now becoming familiar [1, 2]. However, key questions remain because most microbial species in the biosphere resist cultivation in the laboratory [3, 4]. These species are often referred to as microbial “dark matter” [5, 6]. Researchers have pointed out that if one could provide the chemical components of natural environment, uncultured microbes might grow under pure culture condition. To provide these components and allow the culturing of important microbial taxa, some cultivation methods such as culture cocktails , miniaturized cultures , diffusion chambers , and nature’s incubator  have been developed.
These have been great successes; however, many microbial species in the biosphere that would otherwise be “culturable” may fail to grow because of their growth state in nature, such as dormancy , which results in species that are referred to as viable but non-culturable (VBNC) . Sediments can be highly variable in terms of nutrient input and cells might be dormant because of the lack of a major nutrient and/or carbon source even though they naturally occur in these sediments. As a result, nature’s incubator  and diffusion chambers , in which diffusion provides microorganisms with their natural environment incubated in situ, cannot easily capture dormant microbes. Recent studies have shown that marine environments are a “seed bank,” where dormant bacteria are still present in the marine sediments but at much lower abundance . If one could predict the metabolism of an organism and provide the preferred carbon/nutrient sources, one should be theoretically able to switch a dormant to a non-dormant microbe.
Moreover, Lewis and Epstein reported that some isolates merely grew with mixed-culture , thereby demonstrating mixed-culture dependence for some certain isolates . Mixed cultures can provide a simple enough community to research the individual community members , to fully understand the communication of microorganisms, , and represent strategies to grow unculturable bacteria .
In fact, mixed cultures are not a recent idea. Clostridium pasteurianum was isolated through mixed culture by Sergei Winogradsky in 1895 , which were a type of mixed culture. Enrichment is a long-used and common practice and effectively increases the population of target organisms. By using the enrichment culturing method, several rare microorganisms belonging to infrequently isolated or recently described taxa were isolated in our previous study, including Bradymonadales , Marinilabiliales , and Draconibacteriaceae . Furthermore, the enrichment culturing method has been helpful for single-cell genomics analysis . However, the many questions that remain regarding enrichment culturing limit its application in the isolation of microbial dark matter. For example, how and why microorganism communities change during the enrichment culturing is unknown, as is, more importantly, whether “yet-to-be cultivable” strains can be isolated due to their enriched abundance. Clarifying these questions might be helpful for precisely culturing microorganisms in mixed cultures for drug discovery or screening novel natural products that may finally end the antibiotics discovery void.
Culturing the marine sediment microbiotas based on the enrichment method
We developed an enrichment culture method using a low-nutrient medium containing 10 mM sodium pyruvate in the culture (Fig. 1) because low-nutrient media facilitated the cultivation of previously uncultured marine bacterioplankton . Meanwhile, sodium pyruvate was reported to have a resuscitative effect on VBNC Salmonella  and shortened the lag phase of Nitrosomonas europaea during mixed culture stimulation .
The quantities of specific isolates differed at various stages in the enrichment culture, and 2- to 4-fold more different species could be obtained from the same sediment sample after enrichment culturing (Additional file 5: Table S5). The culturable bacterial species isolated from different marine sediment samples were different, and only approximately 20% of the species were isolated in two sediment samples (Additional file 5: Table S5). However, the different samples in an enrichment culture exhibited similarities: (a) Proteobacteria are most abundant in the samples at various times during enrichment culturing. (b) The number of culturable Bacteroidetes isolates increased after enrichment culturing. (c) Some novel species could be isolated only after enrichment culturing, including some members of Marinilabiliales, ε-Proteobacteria, and δ-Proteobacteria. These results suggested that the method can change the flora of cultivable microbial communities and can be very effective in the separation of marine sediment bacteria.
Composition and diversity of microbial communities in enrichment culture
To explore the compositional changes of microbial communities in different enrichment culture stages, the 16S rRNA gene PCR products of the V4-V5 region were sequenced. The barcoded high-throughput sequencing generated 914,059 quality sequences from 15 samples (Additional file 6: Table S6), with an average of 60,937 sequences per treatment. The number of OTUs detected in each sample ranged from 1285 to 2982 (Additional file 6: Table S6). Bacteria were determined to be more relatively abundant in samples from day 0 than later samples by analysis of the ACE (abundance-based coverage estimator), Chao1 estimator, and Shannon and Simpson indexes (Additional file 6: Table S6). The rarefaction curves, which were based on OTU numbers, almost reached the asymptote for different samples (Additional file 7: Figure S1a); the Shannon diversity index curves also reached plateau levels (Additional file 7: Figure S1b). These results suggested that the vast majority of prokaryotes were well represented in these reservoirs. However, new unique OTUs appeared after the enrichment treatment, and the number of novel OTUs was 16.8% greater after the enrichment cultivation than in the original sediment samples (Additional file 7: Figure S1c). Thus, sequencing analyses would miss some low-abundance sequences of rare strains, i.e., the real microbial dark matter. However, enrichment culture could change the abundance of some OTUs, allowing some low-abundant OTUs to reach the detection threshold. Principal components analysis (PCoA) based on weighted UniFrac metrics also showed that the OTUs found in each enrichment treatment were different in their abundance from each other and that the different sediment samples showed different OTUs under the same enrichment treatment (Additional file 7: Figure S1d).
Of the classifiable sequences, 53 phyla were identified, with Proteobacteria, Fusobacteria, Bacteroidetes, Firmicutes, and Chloroflexi representing the most dominant lineages in the samples receiving the 5-stage treatment (Additional file 7: Figure S1e). During enrichment cultivation, the abundance of the phyla Fusobacteria and Firmicutes increased significantly, whereas that of phyla Planctomycetes and Acidobacteria decreased (Additional file 7: Figure S1). At the family level, some anaerobes or facultative anaerobes, such as Marinilabiaceae, Fusobacteriaceae, and Clostridiaceae, were significantly enriched under enrichment culturing (Additional file 7: Figure S1f).
These results strongly indicate that the composition and diversity of microbial communities changed with environmental factors and that novel operational taxonomic units (OTUs) appeared after enrichment, suggesting that marine sediment may contain a “microbial seed bank”  or a set of rare microbes that can barely survive until the proper nutrients arrive [24, 25].
Mechanism of culturing the isolates
It is usually considered that enrichment culturing enhances the probability of isolating uncultured microorganisms by increasing their abundance. However, few systemic analyses of the mechanism responsible for culturing isolates during enrichment exist. Whether the “yet-to-be cultivable” strains are cultured due to their abundance enrichment thus remains unclear.
However, the abundances of some isolates were enriched at certain enrichment stages. The species most typical of this enrichment is Draconibacterium orientale (belonging to a novel genus and a novel family), which was isolated by this method . In our previous long-term studies, this species was recalcitrant to be isolated by directly culturing on MA medium, except after enrichment culture. In addition, strains of this species were repeatedly isolated in these three sediment samples after enrichment. Local BLAST analysis revealed that the abundance of this bacterium was very low in the original samples (usually undetected) but significantly increased with a prolonged enrichment culturing (to 251 reads) (Fig. 2, Additional file 3: Table S3 and Additional file 8: Table S7). Draconibacterium orientale belonged to the family Draconibacteriaceae (it was transferred to family Prolixibacteraceae in 2014) within the order Marinilabiliales  and the phylum Bacteroidetes, which are commonly assumed to be specialized in degrading high-molecular-weight compounds in marine environments . Marinilabiliales contains four families (Marinilabiliaceae, Prolixibacteraceae, Marinifilaceae, and Salinivirgaceae ) and has approximately 45 species (Additional file 8: Table S7). Approximately 70% of species were isolated after the so-called enrichment treatment (Additional file 8: Table S7).
In our previous study, we identified 12 novel species in the order Marinilabiliales , and 10 of those 12 species were isolated using enrichment methods (Additional file 8: Table S7). These results showed that enrichment culturing might be an efficient method for isolating some low-abundance uncultured bacteria by enrichment. Meanwhile, some species, such as Halomonas denitrificans and Marinobacter xestospongiae (Additional file 5: Table S5), were relatively abundant in the original sediment samples but still recalcitrant to be isolated by direct spread onto MA medium, though enrichment culture could aid their isolation. The results were also consistent with that in Salinivirga cyanobacteriivorans isolation .
However, during the enrichment culture, some bacteria could be enriched but could not be isolated on the MA medium. One possibility might be due to the different environmental conditions between enrichment and isolation, such as differences in oxygen concentration, liquid enrichment but solid culture, and nutrient composition changes after enrichment. An alternative possibility might have to do with the “Black Queen Hypothesis” proposed by Morris, Lenski, and Zinser . According to this hypothesis, some bacteria, enriched in the enrichment culture, might need to grow close to other bacteria on a petri plate.
We were interested in why some certain groups of bacteria could be enriched and cultured by this method and why such low-abundance strains could be isolated. When most bacteria in energy-limited environments are slow-growing bacteria (“k-strategists”), they may be prepared for a slow but steady existence under the nutrient limitations. Rapid growth is inhibited because the bacteria would not be ready for the end of a brief nutrient flush . Other bacteria, which are rarer in nutrient-limited environments, respond rapidly to nutrient flushes and have a rapid growth when nutrients are plentiful. However, these bacteria (“r-strategists,” such as most Marinilabiliales species) must re-enter a dormant state before such a nutrient flush is exhausted. Once such bacteria are in a dormant state, they wait for the next nutrient flush  or reaching a VBNC state .
As a result, we hypothesized that the enrichment culture in this study helped some species (“k-strategists”) adapt and grow on MA medium and could aid the resuscitation of dormant or VBNC bacteria (“r-strategists”) (Additional file 7: Figure S2). Those resuscitated bacteria were unable to divide during the enrichment culturing but might have been then able to grow on the MA medium. To assess the role of resuscitation, we calculated the proportions of dormant and active taxa in each sample using the (16S rRNA+ 1):(16S rRNA gene + 1) gene ratio and defined any taxon with a ratio > 1 as resuscitation or active , owing to the observed low abundance in some bacterial diversity . Active taxa remained 27% of the community in original sediments (Additional file 7: Figure S3), but during enrichment culture, the proportion of active taxa decreased over time (Additional file 7: Figure S3) to 13% of the community. However, order Marinilabiliales had considerably higher ratios in enrichment culture compared with original sediments (Fig. 3b). At the family and genus level, Marinilabiaceae, Tangfeifania, Draconibacterium, and Saccharicrinis also had considerably higher ratios in enrichment culture (Fig. 3b), suggesting that most bacteria of Marinilabiliales were dormant in original sediments and that they could be significantly active after enrichment culturing. In addition, these results were consistent with the idea that certain bacteria were more frequently isolated after enrichment culturing.
We also treated 10 strains that had been isolated but with low abundance in the corresponding enrichment culture stage for induction into a VBNC state (Additional file 7: Figure S4a). Five strains could not grow on MA medium after VBNC induction treatment, even after 10 days of culturing; however, these VBNC cells were resuscitated and formed a colony on MA medium after an enrichment culture treatment (Additional file 7: Figure S4b). The other five strains, which were not fully in the VBNC state under our experimental conditions, were also cultured after enrichment, and their numbers of CFUs were significantly increased by enrichment (Additional file 7: Figure S4b). Real-time PCRs showed that most of the resuscitated strains were not enriched during the enrichment culturing (Additional file 7: Figure S4c). These results suggested that the enrichment culture method included important mechanisms for “resuscitating” some marine sediment bacteria.
Metatranscriptomic and comparative genomic analyses of the resuscitation mechanism during enrichment culturing
To determine the in situ metabolic activities of the cultured microbes and to explore the resuscitation mechanism in the enrichment mixed culture system, we used metatranscriptomic and comparative genomic analyses of the bacterial metabolic activity in the enrichment system.
During enrichment culturing, most of these general stress response genes were downregulated (Figs. 4 and 5). In contrast, genes that participate in aspects of cellular metabolism such as the citrate cycle (TCA cycle), glycolysis/gluconeogenesis, the pentose phosphate pathway, nitrogen metabolism, oxidative phosphorylation, sulfur metabolism, and amino acid metabolism were upregulated, suggesting the microbes have high metabolic activity (Fig. 4). Carbohydrate metabolism was one of the most abundantly represented categories and has a key role in cell growth . Expression of genes involved in most of the important steps in glycolysis was detected, and the expression of several glycolysis genes increased significantly in all three sediment samples during enrichment culturing (Fig. 4). The gene encoding acetyl-CoA synthetase was significantly upregulated during the 5 days of enrichment culturing, probably because we added 0.2% CH3COONa, which becomes a substrate of acetyl-CoA synthetase, in the broth before enrichment. Another way to enhance acetyl-CoA production was using the oxidative decarboxylation of pyruvate, which is generated from glycolysis, and the expression of genes related to the pyruvate dehydrogenase system was also significantly upregulated in the first stage (5 days) of enrichment. Acetyl-CoA is a molecule that participates in many biochemical reactions in protein, carbohydrate, and lipid metabolisms . Its main function is to deliver the acetyl group to the TCA cycle to be oxidized for energy production. CH3COONa and acetyl-CoA might thus have potential roles in resuscitation; however, the detailed mechanisms of this process should be examined in our further studies.
In all three sediment samples, the expression of genes involved in the first half of the TCA cycle during different enrichment stages (range of 1.5-fold and 1000-fold) greatly increased, leading to synthesis of 2-oxoglutarate, a key metabolite in the regulation of carbon and nitrogen metabolism  (Fig. 4). The genes related to carbohydrate metabolism were also upregulated. The aerobic respiration pathway genes (the ubiquinol-cytochrome c reductase iron-sulfur subunit and the ubiquinol-cytochrome c reductase cytochrome b subunit) were highly expressed in the day 0 treatment samples, while anaerobic respiration pathways (the cytochrome bd-I oxidase subunit) were upregulated with prolonged enrichment culturing (Fig. 4), suggesting that the conversion of aerobic respiration to anaerobic respiration occurred in this enrichment culture. Genes related to the biosynthesis and degradation of L-serine (D-3-phosphoglycerate dehydrogenase, SerA, more than 100-fold after 5 days of enrichment culturing) and L-threonine (threonine synthase, ThrC, more than 100-fold) were upregulated in the enrichment culture (Fig. 4). Catabolism of l-serine and l-threonine was reported to provide E. coli cells with a source of energy under anaerobic conditions ; enhancement of production and degradation of l-serine and l-threonine is thus probably also a mode of energy generation used by microbes under the anaerobic conditions caused by this anaerobic enrichment culture method.
Likewise, glutamine synthetase gene (glnA) and glutamate synthase gene (gltBD), which encoding the enzymes permitted the incorporation of ammonium into glutamine and utilization of glutamate , were regulated (Fig. 4). Two-component system-related nitrogen assimilation (NtrC family) and peptide and amino acid transport systems (ABC transporters) were also highly upregulated (Fig. 4), suggesting that populating microorganisms had massive demand for nitrogen resources. However, genes encoding enzymes for the utilization of nitrate and nitrite (for example, nitrate reductase and nitrite reductase) exhibited relatively low transcriptional activities, and nitrate reductase was highly expressed during only the late state of enrichment culturing (Fig. 5). These results showed that the enrichment culture communities might acquire nitrogen resources via the direct transformation and utilization of organic nitrogen sources. In addition, genes involved in sulfate reduction, such as those encoding adenylyl-sulfate reductase , were also enriched, which suggests that the main sulfur metabolism in the microbiota might change during resuscitation.
Bacteria within the Marinilabiliales group also showed highly active cellular metabolism in the enrichment culture. Genes encoding pyruvate phosphate dikinase, glycerol-3-phosphate dehydrogenase, glyceraldehyde 3-phosphate dehydrogenase, glycerol kinase, pyruvate-flavodoxin oxidoreductase, and a cell division protein (ZapA) were upregulated (Fig. 5). Pyruvate phosphate dikinase is a key enzyme in gluconeogenesis and enhances the survival of bacteria . Glycerol-3-phosphate dehydrogenase serves as a major link between carbohydrate metabolism and lipid metabolism. Glyceraldehyde 3-phosphate dehydrogenase and glycerol kinase are enzymes that participate in glycolysis. Pyruvate-flavodoxin oxidoreductase participates in the oxidative decarboxylation of pyruvate to acetyl-CoA. Most of those genes affect pyruvate metabolism are related to carbohydrate and energy metabolism and aid resuscitation from the VBNC state , suggesting Marinilabiliales microbes were resuscitated from non-culturable (dormant) states, and the results were consistent with the 16S rRNA:16S rRNA gene analysis (Fig. 3b). We added 0.11% sodium pyruvate to the enrichment culture medium in this study, and the results was consistent with a study that showed that pyruvate might be a key molecule in resuscitation by converting bacteria from a non-culturable state to a growing and colony-forming state .
Genes involved in flagellar assembly and type IV and type VI secretion systems (T4SSs and T6SSs) were upregulated, especially during the late stages of enrichment culturing (21 and 30 days of treatment) (Fig. 4). With the dynamic changes in bacterial communities in this enrichment mix culture, the T6SSs are among the mechanisms that aid in the struggle against other bacterial species . T6SSs participate in a broad variety of functions, including virulence, antibacterial activity, and conquering new territory . Similarly, T4SSs represent a highly diverse superfamily of secretion systems found in many bacterial species and also can deliver a killing toxin to bacterial neighbors so the organism persists in an ecological niche , suggesting that an increase in cell density increases the concentration of T4SS and T6SS components in bacteria and the capacity of the cell to win the competition. These results also support the ability of novel species Bradymonas sediminis FA350 and Bradymonas sp. B210 (with T4SSs) to prey on other bacteria and to reach a relatively high growth rate after enrichment culturing (data not shown).
Cooperation or coordination strategy
Target bacterial culture based on the enrichment culture method
Because most of the bacteria whose abundance increased during enrichment culturing were strictly anaerobic bacteria (Additional file 7: Figure S1e), we also designed an anaerobic culture strategy using the broad-range medium MA with NO3− and vitamin B complex to specifically select for anaerobic or facultative anaerobic bacteria. The almost full-length 16S rRNA gene library sequencing was performed for all of the colonies on the plates under anaerobic culture conditions to analyze the microorganism groups in the anaerobic culture system. Many novel anaerobic or facultative anaerobic species were obtained by this method (Additional file 13: Table S12). Some novel species within the order Marinilabiliales could be isolated at high frequencies due to their abundant OTUs on the MA medium (Additional file 13: Table S12). However, the enrichment culture method is currently not perfect due to the unique culture conditions required by some taxa. To understand and manipulate the microbiome, researchers must dissect and engineer the interactions within these critical communities .
Our workflow enables large-scale culturing, archiving, genomic sequencing, and phenotyping of novel bacteria from marine sediment microbiotas, which were previously considered recalcitrant. Degree of enrichment analysis, dormant and active taxa calculation, VBNC bacteria resuscitation analysis, combined with metatranscriptomic and comparative genomic analyses of the culturability mechanism during enrichment culture showed that the so-called enrichment method could isolate “recalcitrant” strains not only through enriching their abundance but also through resuscitation. In addition, the enrichment culture is a complicated mixed culture system that utilizes a competition, cooperation, or coordination strategy among bacterial communities. Considering that cultivation techniques must evolve from axenic to mixed cultures, comprehending microbial world , we should redevelop an understanding of the classic enrichment culture method. In addition, enrichment culture methods can be developed and used to construct a model for analyzing mixed cultures and exploring microbial dark matter. This simple method still has great potential for cultivating recalcitrant microorganisms. Our streamlined, single-medium approach builds on the considerable efforts of others [4, 9, 11, 58] and provides a new approach to mining marine microbial dark matter based on resuscitation in mixed cultures. In future studies, we hope to develop a series of enrichment culture methods to form a multi-type “enrichment culturomics,” which will then be combined with the current “culturomics”  to reveal many more uncultured microbes.
Sediment collection and bacterial isolation
Three sediment samples were collected from the intertidal zone of Weihai, China (HGD 37°32′0″ N 122°03′59″ E, XSD 37°31′17″ N 122°01′8″ E, and S 37°28′45″ N 121°57′15″ E) on 10 May 2016. Sediment samples were collected by a surface sediment sampler. Following collection, samples were placed in sterile 50 ml plastic Whirl-Pak bags (NASCO, Modesto, CA, USA) and kept cool until they were processed (within 4 h) using the following procedures.
The medium for enrichment culture was a low-nutrient medium that consisted of the following ingredients in seawater: 0.1% NH4Cl, 0.2% CH3COONa, 0.02% MgSO4·7H2O, 0.02% yeast extract, 0.02% peptone, 0.1% EDTA-Na2, and 0.11% sodium pyruvate. The pH of the medium was adjusted to 7.5 and then autoclaved. A 10% (w/v) NaHCO3 solution was filtered, and a 2% (w/v) KH2PO4 solution was autoclaved. Each solution was added to the autoclaved media (10 ml per liter).
Enrichment culture incubation was performed at 25 °C for 0, 5, 12, 21, and 30 days in separate 500 ml sealed glass bottles (filled with medium and 20 g of sediment sample). A detailed protocol is described in the Extended Experimental Procedures.
Culture-independent microbial community composition
We extracted DNA and RNA from enriched samples (from 0 to 30 days of incubation) using commercially available kits and protocols (PowerSoil DNA Isolation Kit and further purification using a Power Clean ProDNA Clean-up Kit from MoBio for DNA, and an RNeasy isolation kit from Qiagen for RNA). Total RNA quality was checked with a Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). DNA contamination was analyzed by PCR with universal archaeal and bacterial 16S rRNA gene primers and ethidium bromide agarose gel electrophoresis. When DNA was present, it was removed by DNase I treatment (TaKaRa, China) for 30 min at 37 °C, followed by lithium chloride/EtOH precipitation overnight at − 20 °C . We synthesized cDNA from RNA using random hexamer primers and a SuperScript III first-strand synthesis kit (Invitrogen). PCR amplification, purification, and pooling and the subsequent pyrosequencing of a region of the 16S rRNA gene were performed following the procedure described by Fierer et al. . We used the primer set composed of 515F (5′-GTGCCAGCMGCCGCGG-3′) and 907R (5′-CCGTCAATTCMTTTRAGTTT-3′), which was designed to amplify the V4–V5 region and was demonstrated in silico to be universal for nearly all bacterial taxa. Sequencing was carried out on a MiSeq platform at the Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China). Detailed protocols are described in the Extended Experimental Procedures.
To assess the proportion of dormant and active taxa during enrichment culturing, we calculated the 16S rRNA+ 1:16S rRNA gene + 1 ratio for each taxa in each sample and defined a taxon as active with a ratio > 1 [30, 31]. Implicit in this calculation is the assumption that if a 16S rRNA sequence is present in the RNA portion of a sample, it must also be present at least one time on the DNA portion but was not detected due to incomplete sequencing. To assess the diversity of active taxa, we calculated the 16S rRNA:16S rRNA gene ratio for each taxon and defined a taxon as active when its ratio was > 1.
Degree of enrichment detection
The dataset was normalized by random resampling, based on smallest sample size (n = 50,443 sequences in this study). To detect the abundance of the isolated strains in each enrichment period, local BLAST analysis of the 16S rRNA gene of isolated strains and 16S rRNA gene libraries of different stages of enrichment from MiSeq results were conducted using NCBI local BLAST, and 16S rRNA gene libraries with 99% sequence similarity (E value<1e-50) were considered the same species.
To assess the degree of enrichment, we calculated the “enriched 16S rRNA gene (abundance of corresponding isolated sample)+1”:“16S rRNA gene (abundance of 0 day)+1” ratio for each taxon in each sample. When the ratio > 1, the species was enriched when it was isolated. In contrast, a ratio ≤ 1 indicates that the species was not enriched when it was isolated. The “+1” is a correction calculation to prevent the error caused by an initial abundance of 0 because an isolated species could have at least 1 read in the initial sample and the isolated sample.
rRNA removal, cDNA library construction, and transgenomics sequencing
The total RNA of samples was subjected to an rRNA removal procedure using a Ribo-zero Magnetic kit according to the manufacturer’s instructions (Epicentre, an Illumina® company). Next, cDNA libraries were constructed using a TruSeq™ RNA sample prep kit (Illumina). Barcoded libraries were paired-end sequenced on an Illumina HiSeq 2500 platform at the Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China) using a HiSeq 4000 PE Cluster Kit and a HiSeq 4000 SBS Kit according to the manufacturer’s instructions (www.illumina.com). Detailed protocols are described in the Extended Experimental Procedures.
Induction of VBNC state and resuscitation of growth
The selected cells were cultured in marine agar 2216 (MA; BD) broth at 30 °C for 24 h with mild shaking at 150 rpm and were harvested by centrifugation at 16,000×g for 3 min. The bacterial cells were washed twice with sterile seawater, harvested again by centrifugation at 16,000×g for 3 min and finally resuspended in sterile seawater to a density of 108 CFU/ml. The cells were then stored at 4 °C for 35 days. The bacterial cell state was monitored over time using the dilution plate count method on MA medium. The VBNC cells were the bacteria that could not form colonies on MA after such treatment. The experiments were conducted in triplicate.
The resuscitation of VBNC cells used the enrichment culture method. VBNC cells were incubated in enrichment culture medium at 25 °C for 5 days, spread onto MA and incubated at 28 °C for 1 day. Culturable cells were enumerated based on the formation of visible colonies on MA medium. VBNC cells directly spread onto MA, incubated at 28 °C for 6 days, were used as a control treatment.
Quantitative real-time PCR
The enumeration of each of the selected bacterial 16S rRNA gene copies of VBNC samples before and after enrichment culture of 5 days was performed via quantitative real-time PCR to evaluate the cell number during the enrichment culture. Gene copy numbers of VBNC sample enrichment cultures at 0 days were used as controls. We used the primer pair composed of 341F (5΄-CCTACGGGAGGCAGCAG-3΄) and 534R (5΄-ATTACCGCGGCTGCTGGCA-3΄) for detection and quantification . Quantification was based on the fluorescence intensity of SYBR green dye, and reactions for each sample were carried out in an ABI StepOnePlus thermal cycler. The reactions were performed in a total volume of 20 μl using Absolute QPCR SYBR Green Rox Abgene/Thermo (Bio-Medicine), 1 μM of each primer, and 10 ng of sample DNA, which was added to each reaction to reach a final concentration of 0.025 μg/μl. The bacterial 16S rRNA gene fragments were amplified using an initial denaturation step at 95 °C for 15 min and then 35 cycles of 15 s at 95 °C, 30 s at 60 °C, 30 s at 72 °C, and 30 s at 80 °C . Then, the Plasmid DNA Standard was constructed by introducing the 16S rDNA gene amplified from Bacillus subtilis into the pMD19 T-Vector (TaKaRa) according to the manufacturer’s instructions. The Plasmid DNA Standard was then gel-purified and quantified. DNA copy number was determined by the concentration and relative molecular weight of the Plasmid DNA. Then, dilutions were made with EASY dilution solution (TaKaRa) to produce concentrations ranging from 1 × 107 to 1 × 101 DNA copies/μl to generate calibration curves. In this study, results with Cq values greater than 35 were considered negative. Error bars in diagrams represent standard deviations from three independent experiments.
Generation of interaction networks
The network analyses were performed using the Molecular Ecological Network Approach (MENA)  pipeline, which was used to generate interaction networks. Read counts were filtered manually to exclude OTUs with > 2 “zero” occurrences across 15 samples and were uploaded to the pipeline. The majority was set to one (a parameter means that only OTUs > 1 could be kept for further analysis), missing data were kept blank, read counts were converted by logarithm, and the Pearson correlation coefficient was used as the similarity measure. To make the correlations more stringent and clearer, we set the threshold value ≥ |0.7| .
Comparative genomic analysis of vitamin-related cofactor biosynthesis pathways
To reveal the potential relationships between the Marinilabiliales species and the interworking groups, we used the taxonomy annotation of metatranscriptome sequences to search for target species, and selected the abundant taxonomy annotation as the reference species. It is hard to define all the exact interworking strains and the genomes during the enrichment culture. Therefore, we used some of the published genomes as the representatives. Along with 5 genomes obtained from this study, totally 20 genomes, which included 10 genomes of organisms within the Marinilabiliales group and 10 genomes of enriched organisms having positive connections with the Marinilabiliales group, were used to make the comparative genomics analysis. (Additional file 12: Table S11). Genome-based analysis of B vitamin-related cofactor biosynthesis pathways in 20 genomes of bacteria was performed using the subsystem-based comparative genomic approach  implemented in SEED/RAST [66, 67] combined with the genomic reconstruction of vitamin-specific transcriptional regulons and the identification of candidate vitamin transporters, as previously described . All of the genomes can be accessed via NCBI (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/genome/), and the detailed accession numbers are shown in Additional file 12: Table S11.
Extensive manual pathway curation was assisted by KEGG orthology assignments from the BlastKOALA annotation tool . Cofactor requirements were asserted by cataloging respective enzymes from annotated genomes and connecting them to metabolic pathways in KEGG  and SEED .
The 16S rRNA gene and 16S rRNA data sets have been deposited in the Sequence Read Archive under accession numbers SRP132130 and SRP133458 for all the samples. The accession numbers of full-length 16S rRNA gene data sets for the five anaerobic cultured plates were SRP132148. The metatranscriptomic sequences for all the samples have been deposited in the Sequence Read Archive under accession numbers SRP133988. All bacterial isolates have been deposited at the Shandong Infrastructure of Marine Microbial Resources hosted by the Laboratory of Marine Microbiology at Shandong University (http://www.sdum.wh.sdu.edu.cn/search.html?itemId=14). The accession numbers of the isolated strains are listed in Additional file 3: Table S3. Any isolates with accession numbers are available upon request.
This work was supported by the National Natural Science Foundation of China (31770002, 41876166, and 31370057), Science and Technology Basic Resources Investigation Program of China (2017FY100302).
DSM and ZJD developed culturing procedures. QYL, DCL, XMW, and MJS carried out culturing and bacterial isolation. DSM, QYL, and XMW performed bioinformatic analyses. DSM and ZJD designed the study. DSM, GJC, and ZJD analyzed data and wrote the paper. All authors read and approved the manuscript.
No conflict of interest exists in the submission of this manuscript, and the manuscript has been approved by all authors for publication. The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Moran MA. The global ocean microbiome. Science. 2015;350(6266):aac8455.View ArticleGoogle Scholar
- Sunagawa S, Coelho LP, Chaffron S, Kultima JR, Labadie K, Salazar G, Djahanschiri B, Zeller G, Mende DR, Alberti A, et al. Ocean plankton. Structure and function of the global ocean microbiome. Science. 2015;348(6237):1261359.View ArticleGoogle Scholar
- Rappe MS, Giovannoni SJ. The uncultured microbial majority. Annu Rev Microbiol. 2003;57:369–94.View ArticleGoogle Scholar
- Rappe MS, Connon SA, Vergin KL, Giovannoni SJ. Cultivation of the ubiquitous SAR11 marine bacterioplankton clade. Nature. 2002;418(6898):630–3.View ArticleGoogle Scholar
- Lok C. Mining the microbial dark matter. Nature. 2015;522(7556):270–3.View ArticleGoogle Scholar
- Rinke C, Schwientek P, Sczyrba A, Ivanova NN, Anderson IJ, Cheng JF, Darling A, Malfatti S, Swan BK, Gies EA, et al. Insights into the phylogeny and coding potential of microbial dark matter. Nature. 2013;499(7459):431–7.View ArticleGoogle Scholar
- Omsland A, Cockrell DC, Howe D, Fischer ER, Virtaneva K, Sturdevant DE, Porcella SF, Heinzen RA. Host cell-free growth of the Q fever bacterium Coxiella burnetii. Proc Natl Acad Sci U S A. 2009;106(11):4430–4.View ArticleGoogle Scholar
- Ma L, Kim J, Hatzenpichler R, Karymov MA, Hubert N, Hanan IM, Chang EB, Ismagilov RF. Gene-targeted microfluidic cultivation validated by isolation of a gut bacterium listed in human microbiome project’s most wanted taxa. Proc Natl Acad Sci U S A. 2014;111(27):9768–73.View ArticleGoogle Scholar
- Kaeberlein T, Lewis K, Epstein SS. Isolating “uncultivable” microorganisms in pure culture in a simulated natural environment. Science. 2002;296(5570):1127–9.View ArticleGoogle Scholar
- Nichols D, Cahoon N, Trakhtenberg EM, Pham L, Mehta A, Belanger A, Kanigan T, Lewis K, Epstein SS. Use of ichip for high-throughput in situ cultivation of “uncultivable” microbial species. Appl Environ Microb. 2010;76(8):2445–50.View ArticleGoogle Scholar
- Connon SA, Giovannoni SJ. High-throughput methods for culturing microorganisms in very-low-nutrient media yield diverse new marine isolates. Appl Environ Microbiol. 2002;68(8):3878–85.View ArticleGoogle Scholar
- Xu HS, Roberts N, Singleton FL, Attwell RW, Grimes DJ, Colwell RR. Survival and viability of nonculturable Escherichia coli and Vibrio cholerae in the estuarine and marine environment. Microb Ecol. 1982;8(4):313–23.View ArticleGoogle Scholar
- Gonnella G, Bohnke S, Indenbirken D, Garbe-Schonberg D, Seifert R, Mertens C, Kurtz S, Perner M. Endemic hydrothermal vent species identified in the open ocean seed bank. Nat Microbiol. 2016;1(8):16086.View ArticleGoogle Scholar
- Stewart EJ. Growing unculturable bacteria. J Bacteriol. 2012;194(16):4151–60.View ArticleGoogle Scholar
- Garcia SL. Mixed cultures as model communities: hunting for ubiquitous microorganisms, their partners, and interactions. Aquat Microb Ecol. 2016;77(2):79–85.View ArticleGoogle Scholar
- Nai C, Meyer V. From axenic to mixed cultures: technological advances accelerating a paradigm shift in microbiology. Trends Microbiol. 2018;26(6):538–54.View ArticleGoogle Scholar
- Little AEF, Robinson CJ, Peterson SB, Raffa KF, Handelsman J. Rules of engagement: interspecies interactions that regulate microbial communities. Annu Rev Microbiol. 2008;62(1):375–401.View ArticleGoogle Scholar
- Wang ZJ, Liu QQ, Zhao LH, Du ZJ, Chen GJ. Bradymonas sediminis gen. Nov., sp. nov., isolated from coastal sediment, and description of Bradymonadaceae fam. Nov. and Bradymonadales Ord. Nov. Int J Syst Evol Microbiol. 2015;65(Pt 5):1542–9.View ArticleGoogle Scholar
- Wu WJ, Zhao JX, Chen GJ, Du ZJ. Description of Ancylomarina subtilis gen. Nov., sp. nov., isolated from coastal sediment, proposal of Marinilabiliales Ord. Nov. and transfer of Marinilabiliaceae, Prolixibacteraceae and Marinifilaceae to the order Marinilabiliales. Int J Syst Evol Microbiol. 2016;66(10):4243–9.View ArticleGoogle Scholar
- Du ZJ, Wang Y, Dunlap C, Rooney AP, Chen GJ. Draconibacterium orientale gen. Nov., sp. nov., isolated from two distinct marine environments, and proposal of Draconibacteriaceae fam. Nov. Int J Syst Evol Microbiol. 2014;64(Pt 5):1690–6.View ArticleGoogle Scholar
- Utturkar SM, Cude WN, Robeson MS Jr, Yang ZK, Klingeman DM, Land ML, Allman SL, Lu TY, Brown SD, Schadt CW, et al. Enrichment of root endophytic bacteria from populus deltoides and single-cell-genomics analysis. Appl Environ Microbiol. 2016;82(18):5698–708.View ArticleGoogle Scholar
- Morishige Y, Fujimori K, Amano F. Differential resuscitative effect of pyruvate and its analogues on VBNC (viable but non-culturable) Salmonella. Microbes Environ. 2013;28(2):180–6.View ArticleGoogle Scholar
- Clark C, Schmidt EL. Effect of mixed culture on Nitrosomonas europaea simulated by uptake and utilization of pyruvate. J Bacteriol. 1966;91(1):367–73.PubMedPubMed CentralGoogle Scholar
- Metcalf JL, Xu ZZ, Weiss S, Lax S, Van Treuren W, Hyde ER, Song SJ, Amir A, Larsen P, Sangwan N, et al. Microbial community assembly and metabolic function during mammalian corpse decomposition. Science. 2016;351(6269):158–62.View ArticleGoogle Scholar
- Lennon JT, Jones SE. Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat Rev Microbiol. 2011;9(2):119–30.View ArticleGoogle Scholar
- Fernandez-Gomez B, Richter M, Schuler M, Pinhassi J, Acinas SG, Gonzalez JM, Pedros-Alio C. Ecology of marine Bacteroidetes: a comparative genomics approach. ISME J. 2013;7(5):1026–37.View ArticleGoogle Scholar
- Ben Hania W, Joseph M, Bunk B, Sproer C, Klenk HP, Fardeau ML, Spring S. Characterization of the first cultured representative of a Bacteroidetes clade specialized on the scavenging of cyanobacteria. Environ Microbiol. 2017;19(3):1134–48.View ArticleGoogle Scholar
- Morris JJ, Lenski RE, Zinser ER. The black queen hypothesis: evolution of dependencies through adaptive gene loss. Mbio. 2012;3(2):e00036-12Google Scholar
- Janssen PH. Dormant microbes: scouting ahead or plodding along? Nature. 2009;458(7240):831.View ArticleGoogle Scholar
- Jones SE, Lennon JT. Dormancy contributes to the maintenance of microbial diversity. Proc Natl Acad Sci U S A. 2010;107(13):5881–6.View ArticleGoogle Scholar
- Kearns PJ, Angell JH, Howard EM, Deegan LA, Stanley RH, Bowen JL. Nutrient enrichment induces dormancy and decreases diversity of active bacteria in salt marsh sediments. Nat Commun. 2016;7:12881.View ArticleGoogle Scholar
- Jones DS, Flood BE, Bailey JV. Metatranscriptomic insights into polyphosphate metabolism in marine sediments. ISME J. 2016;10(4):1015–9.View ArticleGoogle Scholar
- Chen LX, Hu M, Huang LN, Hua ZS, Kuang JL, Li SJ, Shu WS. Comparative metagenomic and metatranscriptomic analyses of microbial communities in acid mine drainage. ISME J. 2015;9(7):1579–92.View ArticleGoogle Scholar
- Abu-Ali GS, Mehta RS, Lloyd-Price J, Mallick H, Branck T, Ivey KL, Drew DA, DuLong C, Rimm E, Izard J, et al. Metatranscriptome of human faecal microbial communities in a cohort of adult men. Nat Microbiol. 2018;3(3):356–66.View ArticleGoogle Scholar
- Moran MA, Satinsky B, Gifford SM, Luo H, Rivers A, Chan LK, Meng J, Durham BP, Shen C, Varaljay VA, et al. Sizing up metatranscriptomics. ISME J. 2013;7(2):237–43.View ArticleGoogle Scholar
- Fujiwara K, Taguchi H. Mechanism of methionine synthase overexpression in chaperonin-depleted Escherichia coli. Microbiology. 2012;158(Pt 4):917–24.View ArticleGoogle Scholar
- Bepperling A, Alte F, Kriehuber T, Braun N, Weinkauf S, Groll M, Haslbeck M, Buchner J. Alternative bacterial two-component small heat shock protein systems. Proc Natl Acad Sci U S A. 2012;109(50):20407–12.View ArticleGoogle Scholar
- Guisbert E, Yura T, Rhodius VA, Gross CA. Convergence of molecular, modeling, and systems approaches for an understanding of the Escherichia coli heat shock response. Microbiol Mol Biol Rev. 2008;72(3):545–54.View ArticleGoogle Scholar
- Zhao F, Wang YT, An HR, Hao YL, Hu XS, Liao XJ. New insights into the formation of viable but nonculturable Escherichia coli O157:H7 induced by high-pressure CO2. Mbio. 2016;7(4):e00961-16.Google Scholar
- Salina EG, Waddell SJ, Hoffmann N, Rosenkrands I, Butcher PD, Kaprelyants AS. Potassium availability triggers Mycobacterium tuberculosis transition to, and resuscitation from, non-culturable (dormant) states. Open Biol. 2014;4(10). https://doi.org/10.1098/rsob.140106.
- Su X, Guo L, Ding L, Qu K, Shen C. Induction of viable but nonculturable state in Rhodococcus and transcriptome analysis using RNA-seq. PLoS One. 2016;11(1):e0147593.View ArticleGoogle Scholar
- Muthusamy S, Lundin D, Mamede Branca RM, Baltar F, Gonzalez JM, Lehtio J, Pinhassi J. Comparative proteomics reveals signature metabolisms of exponentially growing and stationary phase marine bacteria. Environ Microbiol. 2017;19(6):2301–19.View ArticleGoogle Scholar
- Huergo LF, Dixon R. The emergence of 2-oxoglutarate as a master regulator metabolite. Microbiol Mol Biol Rev. 2015;79(4):419–35.View ArticleGoogle Scholar
- Hesslinger C, Fairhurst SA, Sawers G. Novel keto acid formate-lyase and propionate kinase enzymes are components of an anaerobic pathway in Escherichia coli that degrades L-threonine to propionate. Mol Microbiol. 1998;27(2):477–92.View ArticleGoogle Scholar
- Leigh JA, Dodsworth JA. Nitrogen regulation in bacteria and archaea. Annu Rev Microbiol. 2007;61:349–77.View ArticleGoogle Scholar
- Muyzer G, Stams AJ. The ecology and biotechnology of sulphate-reducing bacteria. Nat Rev Microbiol. 2008;6(6):441–54.View ArticleGoogle Scholar
- Zuniga-Ripa A, Barbier T, Conde-Alvarez R, Martinez-Gomez E, Palacios-Chaves L, Gil-Ramirez Y, Grillo MJ, Letesson JJ, Iriarte M, Moriyon I. Brucella abortus depends on pyruvate phosphate dikinase and malic enzyme but not on Fbp and GlpX fructose-1,6-bisphosphatases for full virulence in laboratory models. J Bacteriol. 2014;196(16):3045–57.View ArticleGoogle Scholar
- Hood RD, Singh P, Hsu F, Guvener T, Carl MA, Trinidad RR, Silverman JM, Ohlson BB, Hicks KG, Plemel RL, et al. A type VI secretion system of Pseudomonas aeruginosa targets a toxin to bacteria. Cell Host Microbe. 2010;7(1):25–37.View ArticleGoogle Scholar
- Zhao W, Caro F, Robins W, Mekalanos JJ. Antagonism toward the intestinal microbiota and its effect on Vibrio cholerae virulence. Science. 2018;359(6372):210–3.View ArticleGoogle Scholar
- Grohmann E, Christie PJ, Waksman G, Backert S. Type IV secretion in Gram-negative and Gram-positive bacteria. Mol Microbiol. 2018;107(4):455–71.View ArticleGoogle Scholar
- Zengler K, Zaramela LS. The social network of microorganisms - how auxotrophies shape complex communities. Nat Rev Microbiol. 2018;16(6):383–90.View ArticleGoogle Scholar
- Degnan PH, Taga ME, Goodman AL. Vitamin B-12 as a modulator of gut microbial ecology. Cell Metab. 2014;20(5):769–78.View ArticleGoogle Scholar
- Wexler AG, Goodman AL. An insider’s perspective: bacteroides as a window into the microbiome. Nature Microb. 2017;2(5):17026.View ArticleGoogle Scholar
- Garcia SL, Buck M, Hamilton JJ, Wurzbacher C, Grossart HP, Mcmahon KD, Eiler A. Model communities hint at promiscuous metabolic linkages between ubiquitous free-living freshwater bacteria. mSphere. 2018;3:e00202–18.View ArticleGoogle Scholar
- Romine MF, Rodionov DA, Maezato Y, Osterman AL, Nelson WC. Underlying mechanisms for syntrophic metabolism of essential enzyme cofactors in microbial communities. The ISME J. 2017;11(6):1434–46.View ArticleGoogle Scholar
- Fredrickson JKECOLOGY. Ecological communities by design. Science. 2015;348(6242):1425–7.View ArticleGoogle Scholar
- Coyte KZ, Schluter J, Foster KR. The ecology of the microbiome: networks, competition, and stability. Science. 2015;350(6261):663–6.View ArticleGoogle Scholar
- Tripp HJ, Kitner JB, Schwalbach MS, Dacey JW, Wilhelm LJ, Giovannoni SJ. SAR11 marine bacteria require exogenous reduced sulphur for growth. Nature. 2008;452(7188):741–4.View ArticleGoogle Scholar
- Lagier JC, Dubourg G, Million M, Cadoret F, Bilen M, Fenollar F, Levasseur A, Rolain JM, Fournier PE, Raoult D. Culturing the human microbiota and culturomics. Nat Rev Microbiol. 2018;16(9):540–50.Google Scholar
- Beam JP, Jay ZJ, Schmid MC, Rusch DB, Romine MF, Jennings Rde M, Kozubal MA, Tringe SG, Wagner M, Inskeep WP. Ecophysiology of an uncultivated lineage of Aigarchaeota from an oxic, hot spring filamentous ‘streamer’ community. ISME J. 2016;10(1):210–24.View ArticleGoogle Scholar
- Fierer N, Hamady M, Lauber CL, Knight R. The influence of sex, handedness, and washing on the diversity of hand surface bacteria. Proc Natl Acad Sci U S A. 2008;105(46):17994–9.View ArticleGoogle Scholar
- Wessen E, Nyberg K, Jansson JK, Hallin S. Responses of bacterial and archaeal ammonia oxidizers to soil organic and fertilizer amendments under long-term management. Appl Soil Ecol. 2010;45(3):193–200.View ArticleGoogle Scholar
- Williams HN, Lymperopoulou DS, Athar R, Chauhan A, Dickerson TL, Chen H, Laws E, Berhane TK, Flowers AR, Bradley N, et al. Halobacteriovorax, an underestimated predator on bacteria: potential impact relative to viruses on bacterial mortality. ISME J. 2016;10(2):491–9.View ArticleGoogle Scholar
- Deng Y, Jiang YH, Yang Y, He Z, Luo F, Zhou J. Molecular ecological network analyses. BMC Bioinformatics. 2012;13:113.View ArticleGoogle Scholar
- Kurtz ZD, Muller CL, Miraldi ER, Littman DR, Blaser MJ, Bonneau RA. Sparse and compositionally robust inference of microbial ecological networks. PLoS Comput Biol. 2015;11(5):e1004226.View ArticleGoogle Scholar
- Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, Edwards RA, Gerdes S, Parrello B, Shukla M, et al. The SEED and the rapid annotation of microbial genomes using subsystems technology (RAST). Nucleic Acids Res. 2014;42(Database issue):D206–14.View ArticleGoogle Scholar
- Brettin T, Davis JJ, Disz T, Edwards RA, Gerdes S, Olsen GJ, Olson R, Overbeek R, Parrello B, Pusch GD, et al. RASTtk: a modular and extensible implementation of the RAST algorithm for building custom annotation pipelines and annotating batches of genomes. Sci Rep. 2015;5:8365.View ArticleGoogle Scholar
- Rodionova IA, Li X, Plymale AE, Motamedchaboki K, Konopka AE, Romine MF, Fredrickson JK, Osterman AL, Rodionov DA. Genomic distribution of B-vitamin auxotrophy and uptake transporters in environmental bacteria from the Chloroflexi phylum. Environ Microbiol Rep. 2015;7(2):204–10.View ArticleGoogle Scholar
- Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J Mol Biol. 2016;428(4):726–31.View ArticleGoogle Scholar