Novel Culturing Techniques Select for Heterotrophs and Hydrocarbon Degraders in a Subantarctic Soil

The soil substrate membrane system (SSMS) is a novel micro-culturing technique targeted at terrestrial soil systems. We applied the SSMS to pristine and diesel fuel spiked polar soils, along with traditional solid media culturing and culture independent 454 tag pyrosequencing to elucidate the effects of diesel fuel on the soil community. The SSMS enriched for up to 76% of the total soil diversity within high diesel fuel concentration soils, in contrast to only 26% of the total diversity for the control soils. The majority of organisms originally recovered with the SSMS were lost in the transfer to solid media, with all 300 isolates belonging to Proteobacteria, Firmicutes, Actinobacteria or Bacteroidetes, the four phyla most frequently associated with soil culturing efforts. The soils spiked with high diesel fuel concentrations exhibited reduced species richness, diversity and a selection towards heterotrophs and hydrocarbon degraders in comparison to the control soils. Based on these observations and the unusually high level of overlap in microbial taxa observed between methods, we suggest the SSMS holds potential to exploit hydrocarbon degraders and other targets within simplified bacterial systems, yet is inadequate for soil ecology and ecotoxicology studies where identifying rare oligotrophic species is paramount.

Microbial populations have been proposed as valuable indicators of polar soil health [1][2][3] . They constitute the bulk of the biomass in polar soils and as such, are primarily responsible for ecosystem services such as decomposition, mineralisation, inorganic nutrient turnover and pollutant removal 4 . Rapid sequencing through the use of next generation sequencing (NGS) platforms have evaluated the impact of petroleum hydrocarbons on bacterial community diversity in mangroves 5 , oil fields 6 , the Gulf of Mexico following the Deepwater Horizon oil spill 7 and contaminated soils in the sub-Antarctic 8 . It has been consistently demonstrated that the presence of petroleum hydrocarbons shift the phylogenetic make up of communities [9][10][11] . However, understanding the functional relevance of these phylogenetic shifts remains problematic. Cultured representatives provide greater insight into individual and community microbial functions and have the potential to validate sequence-based observations through functional experimentation.
The uncultured majority are thought to have oligotrophic growth strategies 12 and grow predominantly as micro-colonies [13][14][15] . By mimicking the oligotrophic conditions found in a specific environment, novel micro-culturing techniques have begun to bridge the gap between the cultured, and uncultured majority [16][17][18][19][20][21] . Novel culturing strategies are aimed at limiting nutrient availability, increasing incubation times and providing direct contact with the source environment 13,14,20 . Approaches designed to mimic aquatic environments have resulted in the successful growth of bacteria from many divisions, including the candidate divisions SAR11 17 and OP10 13,16 . The soil substrate membrane system (SSMS) is a novel micro-culturing technique targeted at terrestrial systems 22 . The technique mimics in-situ growth conditions of the soil environment by using environmental soil as the growth substrate, replacing the excessive nutrients characteristic of artificial media. The SSMS has been demonstrated to substantially increase the culturability of soil bacteria with reports of up to 40% of culture-independent diversity captured with the membrane enrichments [21][22][23] .
Macquarie Island is a remote Australian sub-Antarctic territory located in the Southern Ocean (54°37′ 53″ S, 158°52′ 15″ E) (Fig. 1). After a tumultuous history of exploitation from sealers, the island was declared a wildlife sanctuary in 1933 and has been operating as permanent Australian Antarctic research station since 1948. The continual human occupation has led to the presence of contaminants in the air, sea and soil, with petroleum hydrocarbons identified as one of the most prolific 24 . Three major petroleum hydrocarbon contaminated sites have been identified on Macquarie Island surrounding the fuel storage and power generation facilities ( Fig. 1) 25,26 . A bioremediation program began in 2009 to remediate the historical and recent fuel spills, utilising nutrient addition and bioventing to stimulate the petroleum hydrocarbon degradation potential of the indigenous microbial population 25,27,28 .
In order to inform bioremediation, ecotoxicology, and microbial ecology investigations, our aim was to utilise the SSMS to increase the culturability of the bacterial communities present in sub-Antarctic soils. We compared the results from SSMS culturing against direct culturing with artificial media, and both culturing techniques against a culture independent method; pyrosequencing targeting the 16S rRNA gene. Furthermore, by spiking the soil across a range of diesel fuel concentrations, we evaluated how the presence of diesel fuel altered the sub-Antarctic soil bacterial community.

Methods
Site and sampling methods. At Macquarie Island, a soil sample (500 g) was collected from the top 30 cm of an uncontaminated plot (− 54.4974°, 158.9395°) approximately 170 m away from the largest contaminated site at the Main Power House (MPH) (Fig. 1) 25 . The soil was sub-sampled and spiked with Special Antarctic Blend (SAB) diesel fuel over the range of 0-20 000 mg kg −1 ( Table 1). The soil samples were homogenised and stored for 3 months at ambient temperature (10 °C). To determine total petroleum hydrocarbon (TPH) concentrations after spiking, sub-samples (10 g) from each spiked soil were extracted with hexane and assessed by gas chromatography flame ionisation detection (GC-FID), according to the method described in Schafer et al. 29 . As organisms usually respond logarithmically to toxicants, the spiked soils were grouped into the following biologically relevant diesel fuel concentration ranges; control (no fuel additions), low (DL-400 mg kg −1 ), medium (400-5000 mg kg −1 ), and high (> 5000 mg kg −1 ) ( Table 1). An adjacent core was collected to determine the nutrient and chemistry  parameters of the soil according to the methods described in Mooney et al. 30 . After the 3 month incubation period, soil samples were sieved aseptically through a 2 mm metal sieve to remove large particulate matter in preparation for cultivation and micro-cultivation.
Traditional cultivation on artificial media. To 22 . Plates were sealed with parafilm to prevent moisture loss and then incubated aerobically and anaerobically at 10 °C. For anaerobic conditions the plates were placed in an anaerobic chamber (Becton Dickinson, Australia) with anaerobic gas packs (BD). The resulting macro-colonies were sub-cultured (at least 3X) until pure isolates were obtained, paying particular attention to unique morphologies.
Enrichments with the SSMS. For micro-cultivation, 50 μ l of the 1:200 inoculums from above, was added to 10 ml of distilled H 2 O, then filtered onto a 0.2 μ m, white, Isopore PC membrane (Millipore, Australia) with a vacuum pump filtration manifold as described in Ferrari et al. 31 . The top of the sterile TCI membrane served as a barrier between the non-sterile soil slurry and the underside of the sterile PC membrane. For each individually spiked soil sample (Table 1), six replicate SSMS were set up in 6 well plates: three for aerobic and three for anaerobic incubation conditions and incubated at 10 °C. The optimal incubation time of 16 days was determined by evaluating growth membranes for micro-colonies, without overcrowding, as determined through epi-fluorescent microscopy (EFM) 22 . For aerobic conditions, the final culture vessels (6 well plates) were hydrated and sealed to prevent moisture loss during incubation 22 . For anaerobic conditions, the vessels were hydrated and placed in an anaerobic chamber with anaerobic gas packs. After 16 days of incubation, the membranes were removed and microbial growth was confirmed with EFM. After growth confirmation, the remaining membrane quarters were combined for each sample and used for DNA extraction and secondary culturing on artificial media. Secondary cultivation from the SSMS enrichments onto artificial media. Following micro-cultivation, the tips of the remaining quarters of growth membrane were secured in the lid of 1.5 ml tubes and vortexed with 1 ml of filtered Tris-EDTA buffer to remove the cells from the membrane. The buffer was centrifuged for 60 s at 13 000 g to create a pellet and the supernatant was removed. The pellet was resuspended in 300 μ l TE buffer and 100 μ l was used as an inoculum to spread in triplicate on low nutrient media plates (0.1 × RAVAN). This process was repeated for all the samples. The resulting macro-colonies were sub-cultured until pure cultures were obtained.
DNA extraction and RFLP of pure isolates. DNA was extracted from pure isolates by boiling cells at 100 °C for 10 min, followed by centrifugation at 13,000 g for 15 min. To amplify the 16S rRNA gene the resulting lysate was transferred to a clean tube for PCR and 10 ng was added to 50 μ l PCR reactions using primers F27 and R1492 32 . The PCR protocol consisted of initial denaturation at 94 °C for 5 min, followed by 35 cycles of 94 °C for 30 sec, 60 °C for 30 sec and 72 °C for 30 and terminated with a final step of 72 °C for 5 min. A portion (15 μ l) of the 16S rRNA PCR amplicons was digested with the restriction enzymes HinFl and Rsal (Promega, USA) following the manufacturer's instructions. Digested products were visualised on a 2% agarose gel after staining with SYBR safe. Unique RFLPs were selected for 16S rRNA gene sequencing. For RFLPs observed > 5 times, PCR products were sourced from more than one isolate. Sanger sequencing of pure cultured isolates. The PCR products of selected isolates were purified using a PCR purification kit (Qiagen, USA). Approximately 30 ng of the purified PCR product was utilised in a sequencing reaction with 1 μ l BigDye Terminator (Applied Biosystems) and 1.5 μ l 5X buffer, per 20 μ l reaction. The sequencing reaction protocol consisted of 25 cycles of 96 °C for 10 sec, 50 °C for 5 sec and 60 °C for 4 min. Reaction products were cleaned with ethanol and EDTA (Ethylene diamine tetra-acetic acid) and analysed with an ABI 3730 sequence scanner at the Ramaciotti Centre for Gene Function Analysis (UNSW Australia). All sequence data from the ABI sequencer was checked for chimeric artefacts using the program bellapheron http:// foo.maths.uq.edu.au/~huber/doc/doc/bellerophon.pdf. The data was then analysed using the NCBI nucleotide Basic Local Alignment Search Tool (BLAST) for nearest matches and closest cultured isolates in the Genbank Database.
DNA extraction from SSMS enrichments. The mixed bacterial community present on the enriched SSMS growth (PC) membranes were used for DNA extraction. The remaining membrane quarters were added to a 2 ml lysing matrix tube from the FastDNA SPIN kit for soil (MP Biomedicals, USA) and extracted according to the manufacturer's instructions. To confirm the presence of DNA in the extracts a portion of each sample (6 μ l) was electrophoresed on a 2% agarose gel. The DNA extracted from individual growth membranes was combined to generate one aerobic and one anaerobic sample for pyrosequencing at each concentration, for a total of 14 samples (Table 1).
DNA extraction from soil. After preliminary optimisation on the extraction of gDNA from polar soils, the FastDNA SPIN kit for soil (MP Biomedicals, Seven Hills, NSW, Australia) was selected over other commercial soil DNA extraction kits based on higher DNA yields, and the amount of soil used for the extractions was reduced from 0.5 g to approximately 0.25-0.30 g (data not shown). Each soil sample was extracted in triplicate and gDNA was quantified using a Picogreen assay (Life Technologies, Australia), with absorbance measured using a fluorescence plate reader (SpectraMax M3 Multi-Mode Microplate Reader, Molecular Devices, USA). In order to account for extraction, amplification and sequencing variability within samples, all three independently extracted DNA replicates were sequenced for the un-spiked control soil (Table 1).
Pyrosequencing targeting the 16S rRNA gene. The extracted DNA were purified and amplified for 454 FLX titanium pyrosequencing with the 16S rRNA gene universal primers; 28F and 519R according to the methods outlined in Dowd et al. 33 . Raw sequence data was provided in the form of standard flowgram format (sff) files. The sequences and flowgrams were extracted from the sff files, de-multiplexed and error-checked via the Pyronoise algorithm 34 in MOTHUR 35 . Further quality screening of the sequence data was completed according to the pipeline in van Dorst et al. Bacterial seed sequences were aligned to the curated SILVA secondary structure alignment 36 . Aligned 16S sequences were then clustered into OTUs based on 96% sequence similarity as the best definition for species-level OTUs in this region of the 16S gene 37 . Taxonomic assignment of the identified bacterial OTUs was performed using the Greengenes 2011 database trimmed to the same region as our amplicons (V1-V3) 38 . An OTU abundance-by-sample matrix was generated from the bacterial dataset with MOTHUR 35 . As an additional stage of quality control, global 16S rRNA gene sequence singletons were removed from the dataset and sequences were subsampled to the lowest number of reads pre sample (1950). The final processed sequences were deposited into NCBI database under the project accession number PRJNA338201.
Calculation of diversity estimates from pyrosequencing data. To evaluate the overall bacterial diversity recovered from the SSMS following pyrosequencing samples, aerobic and anaerobic sequences were combined for each fuel concentration analysed. Multivariate data analysis was conducted with the software packages PRIMER v6 and Permanova+ . To account for the log-normal distribution of the data, the abundance-by-sample matrix was square root transformed. After transformation the skewness and kurtosis was reduced closer to 0. For sample comparison the abundance-by-sample matrix was then standardised to express the OTU abundances as relative abundances. The community diversity estimates of richness, Pielou's evenness coefficient (J), Shannon (H') and Simpson (λ ) were calculated in PRIMER. To interrogate the shared diversity between the SSMS enrichments and the soil samples, the consistent OTU's were identified and multiplied by their relative abundance. The resulting shared diversity was calculated for each diesel fuel concentration range and expressed as a percentage of the total soil diversity. A two-tailed student's T-test was used to test the null hypothesis that diesel fuel concentration made no difference in the shared diversity between methods.
Response of individual genera to diesel fuel. The log abundance of individual genera was plotted against the log TPH concentration. Linear equations were fitted to the observed dose response and solved to determine the slope of the line. The fit of the linear model was assessed by calculating a sum of squares test with significance set at P < 0.05. The results were summarised in a volcano plot where the x-axis indicated the slope of the line (negative values infer inhibition and positive values infer stimulation), and the y-axis indicated the − log (P value). Calculating the -log of the P values allowed the genera with the greatest linear relationship to TPH to be positioned towards the top of the plot, and individual genera with the most substantial shifts in relative abundances positioned furthest from the middle of each plot.

Results
Comparison of culture dependent verses independent methods to capture bacterial diversity.
As expected, molecular approaches sequencing total gDNA (directly from the soils or SSMS enrichments) recovered far greater numbers of species, genera and phyla than the pure culture dependent methods alone (Fig. 2, Table 2). Of the pure culture dependent approaches, the SSMS and artificial media combination recovered more genera and species than culturing directly from the soil alone. Yet both cultivation approaches recovered members from the same four phyla; Proteobacteria, Firmicutes, Actinobacteria and Bacteroidetes.
Diversity of pure cultured isolates. A total of 300 isolates were recovered from the soil using artificial media and from the combination of the SSMS and artificial media. From the 300 isolates, 69 RFLPs were identified as unique and following sequencing and a total of 46 species from 24 genera were identified ( Table 3). The SSMS and artificial media combination recovered a greater number of species in comparison to the traditional culturing directly from soil, with 32 species from 18 genera isolated compared to 18 species and 11 genera respectively. Of the 69 unique RFLPs identified, 10 were unable to be successfully sequenced. The majority of isolates recovered belonged to the genera Arthrobacter, Sphingomonas and Pseudomonas with 89, 60 and 43 isolates respectively. There were 5 novel isolates unable to be accurately identified below genus level with the following closest matches; Achromobacter xylosoxidans 96%, Arthrobacter psychrochitiniphilus 96%, Arthrobacter stackebrandtii 96% Hymenobacter ocellatus 96% Sphingomonas faeni 95% and Yersinia enterocolitica 96%, and 2 novel isolates were unable to be identified below family level with the following closest matches; Adhaeribacter terreus 91%, Arthrobacter sulfureus 94%.
Of the pure cultured genera, 4 were more dominant in the low TPH contamination samples including Acidovorax, Brevundimonas, Microbacterium and Sporasarcina, while the genera Rhodococcus were particularly prevalent in the mid-range TPH samples ( Table 3). The genera Arthrobacter, Burkholderia, and Sphingomonas were dominant across all TPH ranges. Pseudomonas and Rhodanobacter were also present in all TPH ranges but increased their prevalence in the samples with high TPH, along with Achromobacter. Of the genera recovered more than once, Microbacterium, Paenibacillus and Burkholderia were recovered only with the SSMS to artificial media combination, and Achromobacter was recovered only with the direct culturing from soil.
Pyrosequencing diversity from the SSMS enrichments and soil samples. For the SSMS enrichments a total of 34,308 16S rRNA gene sequences across the 14 samples (Table 1) (Table 1) were obtained. After filtering out the low quality sequences a total of 32,358 sequences (average 3,595 per sample) remained. For further community diversity analysis, all samples were subsampled to 1950 sequence reads.
As expected the species richness was higher in the soil community in comparison to the enriched SSMS community (Fig. 3A). On average, the species richness in the control SSMS enrichments was 17.9% of the richness found in the control soil samples, which increased to an average of 42.16% for the high diesel fuel concentration range. The species evenness (J'), Shannon (H') and Simpson (λ ) diversity estimates were also greater in the soil communities. Within the soil samples, all community estimates measured declined with increasing TPH. In contrast, the response within the SSMS enrichments was negligible, with small increases in the evenness and Simpson estimates and a decline in the species richness. Expressed as a percentage of the total soil diversity, the shared diversity ranged from an average of 25.6-67.1% (Fig. 3B). The percentage of shared OTUs was found to significantly increase from the control soils to the medium and high diesel fuel concentration ranges. Overall, the samples within the high TPH concentration range had the highest percentage of shared diversity at the OTU level (67.1%) and the greatest similarity across all measured diversity estimates (Fig. 3A).

Bacterial phylogenetic diversity from the SSMS enrichments and soil samples.
At the phyla level the SSMS enrichments were dominated by Proteobacteria, contributing 88-91% of the total relative abundance, followed by Firmicutes (2-10% of the total relative abundance), and Actinobacteria, Bacteroidetes, and Acidobacteria (< 5% relative abundance) (Fig. 4A). The total soil diversity was also dominated by species within Figure 2. The number of unique phyla (A) and genera (B) recovered with culture dependent and culture independent techniques. Numbers recovered are grouped by diesel fuel concentration ranges; low (DL-400 mg kg −1 ), medium (401-5000 mg kg −1 ) and high (> 5000 mg kg −1 ). The total extraction of gDNA from SSMS enrichments and directly from soil retrieved considerably more species than the culturing techniques. The pyrosequencing results revealed a peak in unique phyla and genera at mid-level diesel fuel concentrations.  Proteobacteria. However, the detected prevalence was lower than that estimated for the SSMS enrichments, ranging between 49-51% of the total relative abundance for the control, low and medium diesel fuel concentrations (Fig. 4B). Other phyla that were abundant in the control, low and medium diesel fuel range soils included Firmicutes (7-11%), Actinobacteria (9-18%), Bacteroidetes (5-12%), Gemmatimonadetes (2-3%), Acidobacteria (2-5%), and Chloroflexi (4-6%). The phyla Planctomycetes, Verrucomicrobia, Cyanobacteria, Candidate Division TM7 and the remaining phyla denoted 'Others' were also present, contributing less than 5% each of the total relative abundance. In the high diesel fuel ranges, the composition for the soil samples changed dramatically with Proteobacteria contributing to over 80% of the total relative abundance. Apart from Firmicutes at 8% of the relative abundance in the highest fuel concentration, all remaining phyla contributed < 5% each to the total relative abundance observed. Within Proteobacteria, the control, low and medium diesel fuel concentration ranges for the SSMS enrichments were dominated by Rhodoferax (26-35%) and Simplicispira (12-25%), yet within the high diesel fuel range both genera were barely detected (Fig. 4C). Pseudomonas was prevalent across all fuel concentrations but dominated the high-range diesel fuel samples contributing over 60% of the total relative abundance observed. For the high diesel fuel range Dyella (3-10%), Parvibaculum (1-4%) and Rhodanobacter (1%) were also detected (Fig. 4C). Within the soil samples, the community composition of the control, low and medium diesel fuel concentration ranges was found to be vastly different to the equivalent SSMS enrichments at the genera level diversity (Fig. 4D). The genera Simplicispira and Rhodoferax were barely detected and instead Delfuviicoccus was the most abundant genera contributing between 8-20% of the total relative abundance. A greater number of genera were detected in relatively low abundances in soil; Thioflavicoccus (3-8%), Rudaea (2-4%), Pelagibius (3-10%), unclassified species within Xanthomonadaceae (2-5%), unclassified species within β-Proteobacteria (3%), Rhodanobacter (1-3%) and Dyella (1-4%). Additionally, a further 305 genera represented by 'Others' contributed to a combined level of 40-55% of the total relative abundance. In the highest diesel fuel concentration range, there was a substantial decline in the number of total detected genera within the soil samples (Fig. 4D). Although present in different proportions, the major genera recovered within the high TPH concentration range were consistent between the SSMS enrichments and soil samples, in particular Pseudomonas, Rhodanobacter, Parvibaculum and Dyella.
Response of individual genera to TPH. Analysis of the response of bacterial genera present in the SSMS enrichments to TPH revealed 12 genera were significantly inhibited and 6 genera were significantly stimulated (p < 0.05) (Fig. 5, Table 4). In contrast, the genera identified in soil had 71 genera significantly inhibited and only 5 genera significantly stimulated. All of the genera significantly stimulated with increasing diesel fuel concentrations were members of Proteobacteria (Table 4), and the majority of genera were consistent between the two methods (Pseudomonas, Parvibaculum, Herbaspirillum and Dyella). However, the genera significantly inhibited with TPH originated from a total of 6 different phyla, and despite the larger numbers, there were no organisms identified as significantly inhibited within both the soil and SSMS enrichments (Table 4).

Discussion
The SSMS enrichments combined with secondary cultivation onto solid media recovered greater numbers of genera and species, and a greater number of novel isolates than artificial media alone. However, neither of the cultivation-dependent methods recovered isolates belonging to phyla outside Proteobacteria, Firmicutes, Actinobacteria or Bacteroidetes (Table 2). Organisms belonging to these four phyla are the most commonly cultured species from the environment, including extreme polar environments 39 . However, culturing soil bacteria from phyla outside of these four well-characterised phyla is possible, and has been achieved through the use of dilute media and extended incubation times 40 . Novel cultivation methods that mimic the natural environment continue to reduce the gap between culture-dependent and culture-independent approaches 21,23,41 , expanding our knowledge Here, the SSMS enrichments generated a considerable proportion (18% to 42%) of the total species richness, which is consistent with reported recoveries 21 . Yet, only a limited percentage (~7%) of the OTU's were consistent  with the control soils and neither of the two most dominant Genera (Simplicispira and Rhodoferax) were recovered. Additionally, only a small proportion of this diversity was successfully transferred to solid media, comparable to results obtained when using artificial media alone. In the future, further dilutions of our low nutrient SAB media, combined with dilution-to-extinction approaches may help improve the success rate of transferring micro-colonies from the SSMS enrichments to isolated macro-colonies. Furthermore, we suggest that combining the SSMS with single-cell applications such as flow cytometry 42,43 and micro-manipulation 43 , is a more appropriate method of domesticating rare organisms from extreme environmental samples such as polar soils.
The addition of diesel fuel is known to reduce phylogenetic diversity and selects for hydrocarbon degraders and heterotrophs 11 . This trend was expected and observed across all diesel fuel concentration ranges and methods utilised. Within the high diesel fuel soils, both the cultivation dependent and independent methods recovered Pseudomonas and Rhodanobacter, while the SSMS and total soil communities also recovered Dyella and Parvibaculum [8][9][10][11] . The selective pressure from the addition of diesel fuel also reduced the total species richness, evenness and diversity within the soils, as calculated by Shannon, Simpson, and Pileu diversity indices (Fig. 3). Overall, the SSMS enrichments exhibited substantially lower species richness, evenness, and phylogenetic diversity than the corresponding soil samples. This was particularly true for the control and low diesel fuel concentration ranges. However, the reduction in community complexity, combined with the selection towards heterotrophs following diesel fuel additions, resulted in the SSMS enrichments more closely reflecting both the phylogenetic  diversity and community structure of the high diesel fuel concentration samples (Fig. 3). The three most dominant genera (Pseudomonas, Rhodanobacter and Dyella) were consistent between both methods and made up and average of 60% and 70% of the total relative abundance in the SSMS enrichments and soil samples respectively. Overall, 67% of the total diversity was consistent between methods at the highest fuel concentrations, in contrast to just 26% for control soils. Based on these observations, we believe the SSMS may be a valuable tool to investigate the bacterial community, and to target hydrocarbon degraders or other species of interest within simplified heterotrophic communities 21,23 . The majority of genera identified from the soil were present in very low abundances and were significantly inhibited (p < 0.05) with increasing diesel fuel concentrations (Fig. 5). This is consistent with previous observations 8 . There is evidence that vital soil ecosystem services are carried out by organisms present in low abundances and that these organisms may be particularly vulnerable to disturbances 8,44 . For example, when examining the ecophysiology of anaerobic phototrophic bacteria, Musat et al., revealed that species representing approximately 0.3% of the total cell number contributed more than 40% of the total uptake of ammonium and 70% of the total uptake of carbon in the system. While in van Dorst et al., species essential to nitrification present in sub-Antarctic soils at < 1% relative abundances, were shown to be the most sensitive indicator of toxicity to diesel fuel. Thus, in order to capture the breadth of microbial ecosystem functions and the ecological implications of contamination and other disturbances, it is essential to be able to identify the depth and breadth of organisms present in the ecosystem in question.
Consistent with traditional culturing techniques 12 , the SSMS appears to maintain a selective bias towards heterotrophic species. While the bias towards heterotrophs remains, we suggest the SSMS is unsuitable to characterise total soil bacterial communities. As discussed earlier, the capacity to recover species in low abundance and species with oligotrophic growth strategies is crucial for ecotoxicology 4,8,45,46 and soil ecology investigations [47][48][49][50] . To that end, extracting total gDNA directly from the soil and targeting taxonomic markers on one of the current sequencing platforms remains the ideal method to profile microbial communities. It is important to note that except in the case of functionally specific clades, phylogenetic diversity offers limited information on the functional capabilities and the ecological relevance of community shifts that occur in response to contamination events. At this stage, representative functional data remains limited (particularly in underrepresented environments) for the classification of microbial communities into biologically meaningful categories 51 . Alternatively, methods that target the functional potential or ecological significance of microbial dynamics such as qPCR 8 , network analysis 45 , and shotgun metagenomics are useful strategies to combine with pyrosequencing in order to further illuminate the ecological relevance of microbial interactions in soil.