Temporal dynamics in microbial soil communities at anthrax carcass sites
© The Author(s). 2017
Received: 9 June 2017
Accepted: 13 September 2017
Published: 26 September 2017
Anthrax is a globally distributed disease affecting primarily herbivorous mammals. It is caused by the soil-dwelling and spore-forming bacterium Bacillus anthracis. The dormant B. anthracis spores become vegetative after ingestion by grazing mammals. After killing the host, B. anthracis cells return to the soil where they sporulate, completing the lifecycle of the bacterium. Here we present the first study describing temporal microbial soil community changes in Etosha National Park, Namibia, after decomposition of two plains zebra (Equus quagga) anthrax carcasses. To circumvent state-associated-challenges (i.e. vegetative cells/spores) we monitored B. anthracis throughout the period using cultivation, qPCR and shotgun metagenomic sequencing.
The combined results suggest that abundance estimation of spore-forming bacteria in their natural habitat by DNA-based approaches alone is insufficient due to poor recovery of DNA from spores. However, our combined approached allowed us to follow B. anthracis population dynamics (vegetative cells and spores) in the soil, along with closely related organisms from the B. cereus group, despite their high sequence similarity. Vegetative B. anthracis abundance peaked early in the time-series and then dropped when cells either sporulated or died. The time-series revealed that after carcass deposition, the typical semi-arid soil community (e.g. Frankiales and Rhizobiales species) becomes temporarily dominated by the orders Bacillales and Pseudomonadales, known to contain plant growth-promoting species.
Our work indicates that complementing DNA based approaches with cultivation may give a more complete picture of the ecology of spore forming pathogens. Furthermore, the results suggests that the increased vegetation biomass production found at carcass sites is due to both added nutrients and the proliferation of microbial taxa that can be beneficial for plant growth. Thus, future B. anthracis transmission events at carcass sites may be indirectly facilitated by the recruitment of plant-beneficial bacteria.
KeywordsBacillus anthracis Metabolism Metagenomics Semi-arid Shotgun sequencing Taphonomy Time-series analysis Sporulation Microbial diversity
The microbial composition of arid soils across the globe is distinct from other soil environments [1, 2]. The forces that shape arid soil microbial community composition include low water availability, temperature and UV radiation [3–8]. Especially, water restriction has a large influence since it affects several environmental factors, such as salinity, pH, and the availability of (in-) organic matter, which further modulate soil microbial diversity and activity [4, 5, 9–11]. The combination of the above factors may explain why soils in arid environments, such as deserts and arid savannahs, are taxonomically distinct from other soil types .
Most research on arid soil microbial diversity is directed to understand the community dynamics under changing environmental conditions, e.g. precipitation changes. However, soil communities can also be affected by the deposition of animal carcasses and their decomposition. This may have a profound impact on the soil microbiome as the carcass influences both biotic (adding new microbes) and abiotic (adding nutrients and moisture) factors. Moreover, after the introduction of a carcass, a succession within the microbial community will occur, changing abundances in accordance with nutrients being released during carcass decomposition . For instance, bacteria belonging to Proteobacteria and Acidobacteria are the most common in soils during the initial stages of carcass decomposition, while Firmicutes are more prominent during active decomposition [13, 14]. The described succession is however, found under experimentally controlled conditions, where carcasses were secured at one location for the entire experiment. In contrast, in natural ecosystems carcasses are often consumed and/or dragged away from the site of death by scavengers . Thus microbes and nutrients may only transiently enter the soil at the site of death , which might induce different microbial soil dynamics at natural carcass sites.
The present study investigates the effects of animal carcasses on soil microbial communities after an animal has died of an anthrax infection. The disease anthrax is caused by Bacillus anthracis, a gram-positive, rod shaped, sporulating bacterium . This species belongs together with Bacillus cereus, Bacillus thuringiensis and several other Bacillus spp., to the B. cereus group , which is commonly found in soil as vegetative cells or as spores . The B. cereus group bacteria are indistinguishable from each other using 16S rRNA gene sequences and show high genetic identity (> 99.6%) for certain housekeeping genes [18, 20, 21]. However, B. anthracis can genetically be distinguished from the other ‘species’ based on single nucleotide polymorphisms (SNPs) (e.g. in the plcR gene [22, 23]) and in most cases by the presence of two virulence plasmids specific to B. anthracis pXO1 and pXO2 [24–26]. The lifecycle of B. anthracis is different from other B. cereus group bacteria in that it predominantly targets, infects and kills mammalian herbivores, instead of insects as found for B. thuringiensis. Grazing by herbivores seems to be the main route for transmission in natural settings such as found in the semi-arid savannah in Etosha National Park (ENP), Namibia . After causing the death of its host, B. anthracis returns to the soil where it sporulates. Hence, B. anthracis spores will often be found in high densities in the top layer of soil where haemorrhagic fluids have leaked from an anthrax carcass [27–29].
Long term measurements at carcass sites in the ENP have identified several processes occurring in and above the soil. It was found that B. anthracis cell counts in rhizosphere soils increased in the second year after carcass deposition, but not in surface soils . Those results contrast experimental work where no multiplication of cells was observed in the second year . Furthermore, it was observed that grass biomass and quality increased at the localised area of anthrax carcass sites . This increase in localised plant growth is mostly likely due to nutrient release and/or rhizosphere plant-microbe interactions. The increase in above ground plant biomass results in attraction of grazers, which increases the potential exposure to the pathogen through ingestion of grasses and potentially roots / soil at carcass sites [16, 27]. It is thought that the bacterium sporulates shortly after entry into the soil and that it will stay dormant until favourable conditions arrive . This suggests that there is only a short period, after the entry of B. anthracis into the soil, where the pathogen could interact with other microbes or plants to enhance its transmission [16, 30–33]. It is unclear if such interactions influence further transmission of the pathogen. In particular there is a lack of understanding of the dynamics of the (arid) soil microbial community after the influx of animal fluids (e.g. blood, gut contents) with high densities of B. anthracis vegetative cells. In order to address how interactions between microbes / plants and B. anthracis benefits pathogen transmission, it is important to understand the dynamics of the arid soil microbial community after carcass deposition and influx of B. anthracis vegetative cells into the soil.
There are several technical challenges in following both the B. anthracis and microbial community dynamics in soil. Currently, DNA-based methods such as 16S rRNA amplicon sequencing and shotgun metagenomics are established methods for the study of complete microbial communities . As mentioned above B. anthracis is indistinguishable from other B. cereus group species on 16S rRNA gene sequences. Therefore, 16S rRNA amplicons are not suitable for studies aiming to distinguish members from this group. B. anthracis and many other microbes can be detected with shotgun metagenomics . Again, B. anthracis is genetically highly similar to B. cereus group bacteria, which makes bioinformatic identification of these species challenging . Therefore methods able to detect and distinguish various B. cereus group species are needed to complement and confirm the shotgun metagenomic results. Such methods include cultivation using selective media (e.g. PLET) and specific qPCR assays [23, 36]. The combination of these methods can then be used to monitor environmental B. anthracis populations and to generate hypotheses about possible species interactions, which subsequently can be tested with specifically designed experiments.
Here we present an analysis of soil microbial community composition following leaching of fluids into the soil from two spatially and temporally proximate zebra anthrax carcasses in ENP. The ENP is a semi-arid savannah environment, which for much of the year has sparse water availability, meagre vegetation and limited nutrient resources [37, 38]. Hence, carcass nutrients will likely have great influence on the localised soil microbial community . Large herds of ungulates ensure that every year there is an abundance of carcasses in ENP, often through predation or diseases such as anthrax. It is estimated that up to 400 plains zebras (Equus quagga) per year die from anthrax infections in the ENP . Due to the unpredictability of the presence of disease cases, and our interest in controlling for temporal and spatial variability in soil microbiota among study sites, our study was limited to two carcass sites situated proximately in space and time. Our aim was to describe the microbial community succession taking place in the soil after the influx of haemorrhagic fluids containing B. anthracis vegetative cells. Shotgun metagenomic sequencing of samples obtained through the first month of decomposition was employed to investigate temporal dynamics of the community structure and function. We investigate the temporal change of the taxonomic composition and the metabolic potential by identifying major metabolic pathways. Finally, we use a combination of techniques (qPCR, cultivation and metagenomic sequencing) to circumvent state (vegetative cell vs. spore) associated challenges to accurately track B. anthracis abundances.
Carcass information and rainfall recording
On 03.03.2014 we identified two plains zebra (Equus quagga) carcasses less than 1 km apart in ENP. Hereafter referred to as Carcass 1 (Ca1) and Carcass 2 (Ca2). Ca1 was intact when sampled on day 0, while Ca2 was minimally scavenged (some intestine was dragged out the anus by vultures). The time of collection (≈ 14:00), the state of both carcasses and only the presence of a few avian scavengers and not mammalian, indicates that both animals were likely dead for less than 12 h and certainly fewer than 24 h. Our study contrasts with other carcass decomposition experiments, since scavenging was not restricted [12–14, 39, 41]. As such, the carcass nutrients and fluids will only leak onto the soil for a short time before scavengers consume soft tissue and move the remains off the site.
After three days both carcasses were completely consumed by scavengers and the bones were found approximately 5 m away from each sampling site. The sampling sites were visited at days: 0, 3, 7, 14, 21 and 30, to collect material for cultivation and DNA extraction. At day zero an uncontaminated control sample (Ctrl0) was taken at both sites to function as a reference sample. The soil in the study area has a alkaline pH around 8.7 ± 0.4 (Additional file 1: Table S1), low moisture and dominance of bacteria (Additional file 2: Figure S1), which is characteristic of arid soils .
There was some rainfall at Okaukuejo in the days prior to day 0 followed by little rain until day 18 when heavy rainfalls occurred (day 18–21, Additional file 3: Figure S2).
Detection of B. anthracis
DNA extraction efficiency of bacterial gram-positive cells and/or spores can be poor . To control for extraction efficiency we spiked soil samples with 3.4 X 106 B. anthracis vaccine strain Sterne 34F2 spores (Onderstepoort Biological Products).
In order to distinguish B. anthracis metagenome reads from those of B. cereus and B. thuringiensis, mapping of reads was performed using the aln algorithm (Burrows-Wheeler Aligner (BWA-aln))  with very strict mapping parameters (Additional file 4: Methods). In addition, we added closely related Bacillus spp. reference strains as bait for sequences not unique to B. anthracis . (Note that without the usage of closely related reference strains, mapping of metagenomic reads against B. anthracis becomes highly unreliable). Such mapping against the genomes of two B. anthracis strains, K1 and K2 isolated from the carcasses studied here , resulted in a similar pattern as observed in B. anthracis specific qPCR experiments (Fig. 1a, d) with abundances peaking at day 0 for Ca1 and at day 3 for Ca2. For both carcasses there are no significant differences between the mapping results when using the K1 (Fig. 1b, e) or K2 (Fig. 1c, f) genomes. There are few reads (2–94) mapping to B. cereus, B. thuringiensis or B. subtilis in the Ctrl0 samples (Fig. 1, Additional file 5: Table S2). The frequencies for these three species peak at day 3 at both carcass sites, but abundances are lower than for B. anthracis (Fig. 1b–f). After day 3 the numbers of reads mapping to all Bacillus spp. decrease gradually, except for a slight increase in Ca1 at day 30. A close inspection of the mapping process revealed that many reads with low mapping qualities, e.g. poor alignments due to mismatches, were removed in our final filtering step (Additional file 4: Methods). Removal of reads with mismatches did not change the abundance pattern of our time-series for B. anthracis, while there was a significant change for the other species used (Additional file 6: Figure S3). This difference is likely due to the presence of DNA sequences in the metagenomes related, but not identical to the B. cereus group genomes available and used as references.
Culturing from soil samples that had been heated to kill vegetative cells, showed no B. anthracis spores present at days 0 and 3 (Fig. 1g) in contrast to the qPCR results. Ca1 had >100,000 spores on day 7 and 14, followed by a sharp drop in spores at day 21, and an increase to >50,000 spores at day 30. Ca2 had >2000 spores on day 7 and >60,000 spores on day 14; at the rest of the time-points Ca2 only had <130 spores. Some of these aberrations from the general trend in counts may be products of sampling biases. In addition, the reduction of B. anthracis levels on day 21 (Fig. 1g) could also be due to a response to the rainfall prior to day 21 (Additional file 3: Figure S2).
Temporal changes of GC content and average genome size
Overview metagenome shotgun samples
Raw PE reads
Cleaned PE reads
Cleaned Singleton reads
Average GC-content (%)
Average Genome size (Mb)
The drop in GC content corresponded with a drop in average genome size (AGS) and an increase in genome equivalents (Table 1) . The AGS drop and the shift in average % GC suggests that the carcass soil community composition is changing due to actively growing microbes responding to the carcass nutrients. In addition, the drop in AGS may indicate that the metabolic capacity of the soils changes over the time-series, since AGS can be used as a measure for the metabolic complexity of an ecosystem .
Taxonomic classification of shotgun sequences
Many tools exist for the taxonomic classification of metagenomic shotgun sequences. Our aim was to use a tool that captures most of the diversity present in the samples. We therefore compared metaxa2 , which classifies rRNA sequences in the metagenomes using a modified Silva SSU /LSU database, with metaBIT (wrapper around Metaphlan2) [51, 52], Kraken  and MEGAN . Metaphlan2 and Kraken uses a database derived from microbial whole genome sequences, with the difference that Metaphlan2 only uses signature sequences to identify taxa, instead of whole genomes. MEGAN uses the Non-Redundant protein database from NCBI (NR) for classification.
By comparing the relative abundance of reads classified by each tool we found that the tools differed in the amount of reads classified per sample (Additional file 7: Table S3; Additional file 8: Figure S4) . Interestingly, Kraken showed different temporal dynamics with respect to the change in relative abundance of classified reads compared to the other tools. MetaBIT shows an especially large change between days 0 and 3 compared to the other tools (Additional file 8: Figure S4). We assume that in both cases these results are due to the databases used by metaBIT and Kraken. This is also reflected in the number of prokaryotic taxa identified, with MEGAN detecting a maximum of 180 prokaryotic orders, while the three other tools detected fewer taxa (metaBIT:18, Kraken: 121 and metaxa2: 159). This shows that detection of environmental prokaryotes is significantly influenced by database choice.
The results of the above comparison suggest that metaxa2 is the best available choice among the four tools tested here to describe the bacterial and eukaryotic composition of the soil communities in our study (Additional file 7: Table S3), since the database it relies on is both well curated and contains the widest taxonomic range available. Moreover, MEGAN is a good complement to metaxa2, since it uses most reads for taxonomic and functional classification and can provide diversity measures that support community comparisons (e.g. PCoA, clustering etc).
Temporal dynamics of the soil microbial communities
Several orders show a sharp abundance increase after day 0 (log5 average fold change >1) at both carcass sites. These include Bacillales and Pseudomonodales species. In contrast, bacteria belonging to orders such as Frankiales, Rhizobiales and Solirubrobacteriales do not show large changes in their read abundances after the blood/nutrient influx at either carcass site. This is despite the large increase in total reads classified (Additional file 2: Figure S1). Finally, there are also orders that show limited increases. The abundance increase can be early (T = 0 days) and then decline (e.g. Clostridiales, Bacteriodales), or increase only later (T ≥ 7 days) in the time-series (e.g. Burkholderiales, Corynebacteriales).
The orders Bacillales and Pseudomonadales have the highest abundances at days 3 and 7 (Additional file 10: Figure S5), with genera such as Acinetobacter, Lysinibacillus and Kurthia being dominant (Additional file 9: Table S4). These orders are present at all time-points for both carcasses, but are especially abundant on days 3–14. Many of the genomes in the NCBI database from these genera are small (Size <3 Mbp) and their presence can explain the drop in average genome size (AGS) that we observed at day 3 (Table 1).
Between days 14 and 21 the relative rRNA read abundance for Bacillales drops at Ca1 and Ca2 by 33 and 23%, respectively. For Pseudomonadales, we identify a drop of 17 and 13% for Ca1 and Ca2. In total 9 orders show a relative rRNA read abundance drop (>10%) by day 21 in Ca1 and 12 orders in Ca2. At day 30 the abundances of both Bacillales and Pseudomondales increased again in Ca1, but not in Ca2. Interestingly, for the “stable” orders there was an overall increase in relative rRNA read abundances at day 21 and a subsequent drop at day 30.
The abundance changes observed at days 21 and 30 are likely influenced by heavy rainfall (16–103 mm) recorded at all the weather stations in the ENP in the three days prior to sampling day 21 (Additional file 3: Figure S2). Thus, we assume that the abundance increase at day 30 for many orders could be a response to this heavy rainfall. Moreover, the precipitation could explain the species abundance profile changes since it can result in short-term changes in the microbial community .
In addition to bacterial orders that are normal for soils, we also observed representatives of several bacterial orders likely introduced from the zebra. For instance, Fusobacteriales are present in the first days for both Ca1 and Ca2, where they are most abundant on day 0 for both samples. They are completely absent on day 14 and 30 in Ca1 and from day 7 onwards in Ca2. The increase in abundance for this order on day 0 is most likely due to the introduction of Fusobacterium spp. such as Fusobacterium equinum, which is a known inhabitant of the oral cavity and lower respiratory tract of horses . In Ca1 sequences classified to the order Pasteurellales are highly abundant on day 0, but completely disappear after day 7. Their abundance at day 0 is mainly caused by Actinobacillus spp. and Pasteurella caballi (Additional file 9: Table S4). The later species is a commensal of the upper respiratory tract of horses .
Microbial community metabolism
Shotgun metagenomic sequencing, which allows detection of seemingly all species in a community based on the presence of their DNA, is generally considered a superior method to culturing [1, 62, 63]. This is because it addresses “the great plate count anomaly”, which claims that only about 1% of the microorganisms seen during microscopy can be cultured [64, 65]. However, DNA-based analyses of microbial communities also have technical pitfalls, particularly when studying spore-forming microorganisms such as B. anthracis, which can result in large discrepancies between actual and estimated abundances . For B. anthracis this is further confounded by the high levels of DNA sequence identity shared with other members of the B. cereus group.
The technical pitfalls of DNA-based analyses are evident from our data, where B. anthracis spores could readily be cultured from our samples, but were not as easily detected in metagenomes without using a very strict mapping approach. The onset of B. anthracis sporulation occurs within the first 72 h after carcass deposition and can continue up to eight days post-mortem [29, 66]. The qPCR data (Fig. 1) indicates that over time B. anthracis abundance peaks before almost completely disappearing from the samples. These results, together with the cultivation data, suggest that the reduced quantities of B. anthracis observed in the qPCR experiments are due to reduced DNA extraction efficiency (because of sporulation) . This also suggests that abundances of B. anthracis and similar organisms are likely underestimated in the metagenomes at several time-points due to low extraction efficiency. Nonetheless, the metagenomes still captured a large part of the microbial community present in the carcass site soils, and can therefore be used to study the temporal changes taking place.
Why does a carcass promote plant growth?
The data generated reflects the temporal dynamics of the soil microbial community in a natural ecosystem after inoculation of a pathogen (B. anthracis), other host-associated microbes and nutrient influx from zebra carcasses. Turner et al.  showed that carcass sites in the ENP have higher quality and more abundant vegetation, which could be related to higher levels of phosphate, nitrogen and lower pH than the surrounding soils.
The microbial communities at both carcass sites show a clear response to the influx of bodily fluids, with similar community composition (at order level) of the dominant taxa for the two carcasses with only minor variation in abundances (Fig. 4). Shifts in the community are indicated by several observations: changes in relative abundance of classified reads between the samples (Additional file 3: Figure S2), changes of the AGS of the community studied (Fig. 5) and by variation in OTU diversity and abundance for each of the time-points (Fig. 2). The AGS decreases in the first week of the sample period before increasing again, suggesting a community shift towards copiothrophic bacteria in the first week (Table 1). This is expected as copiothrophic bacteria are better adapted to local increases in nutrient availability than typical soil bacteria and therefore can increase drastically in a short amount of time as observed here .
In contrast to the copiothrophic bacteria, the typical soil microbial community consisted of orders like Frankiales, Rhizobiales and Solirubrobacterales (Fig. 4), which remain relatively stable throughout the sample period. This suggests that these bacteria do not need, or are unable to use, the nutrients from the carcass. Moreover, Rhizobiales are nitrogen-fixing bacteria living in symbiosis with plants that thrive under nitrogen poor conditions . The extra nitrogen provided by an animal carcass would make them less competitive in soils. Frankiales spp. can be plant symbionts, but have also been found to grow on rock beds, which indicate a lifestyle specialized for oligothrophic or extreme conditions [70, 71]. In that respect, it is reasonable that these oligothrophic taxa do not react dramatically to the nutrients of a decomposing animal.
Among the copiotrophic bacteria we see an increase in genera like Acinetobacter, Lysinibacillus and Kurthia in the first week, which is similar to observations in other decomposition studies . These genera belong to the orders Pseudomonadales and Bacillales and are known to contain species that can be either pathogenic or beneficial for both animals and plants [72–76].
Interestingly, many of the orders reacting quickly to the nutrient influx remain at a relatively high level throughout the sampling period with variations in abundance of genera within orders (Fig. 4). For example, Xanthomonadales is abundant at relatively stable levels from day 3 onwards, but the genus Wohlfahrtiimonas within this order is only abundant on day 7 at Ca1. Wohlfahrtiimonas is a known parasite of different Wohlfahrtia spp. (flesh flies) and other such insects. The abundance of Wohlfahrtiimonas seen on day 7 for Ca1 may be a result of an increased number of insect hosts appearing through the decomposition process. This example indicates that short-term species composition changes might be due to local conditions in the soil and are not necessarily due to direct competition or sample variation.
Nonetheless, after day 7 other bacterial orders, such as Burkholderiales, increase in abundance. This order is typically found in higher abundances in the rhizosphere of plants than in the surrounding soils because they actively feed on plant root exudates [77–79]. Furthermore, Burkholderiales are known to have members that are plant growth-promoting bacteria (PGPB), like the Bacillales and Pseudomonadales genera described above [80–82]. Turner et al.  showed that decomposition of carcasses improves soil fertility and enhances plant growth. Our results suggest that this process is promoted by the soil microbial community, which shows a sharp increase in the abundance of orders known to have PGPBs. The promoting functions of PGPBs are many, such as nitrogen fixation, siderophore and phytohormone production, phosphorus solubilisation, and the suppression of plant diseases . Bacillus anthracis is an obligate-killer pathogen, able to transmit only by killing its host. Interestingly, by taking advantage of the nutrient- and microbial-driven stimulation of plant growth occurring at carcass sites, transmission of B. anthracis is enhanced since this plant growth is an attractive food source for herbivores .
Microbial metabolism in arid environments
To get a functional overview of the microbial communities in our study we investigated how metabolic pathways correlate with changes in AGS over time (Fig. 3). Several of the positively correlating pathways can be linked to environmental stressors typical for arid soils, such as desiccation, high temperatures and (UV) radiation. For instance, starch and sucrose (S&S) metabolism is essential for the synthesis of the disaccharides sucrose and trehalose, which can play a role in prokaryotic desiccation resistance . The importance of DNA repair in the surface soil exposed to high radiation levels is reflected in the positive correlation of the non-homologous end joining (NHEJ) repair system, which is also needed for survival during quiescent states such as in the spore stage [84–86].
The presence of other positively correlating pathways can be explained by their role in nutrient acquisition. Arid soils have low organic matter content and most of it comes from decaying plant material . The positive correlation of the aminobenzoate degradation pathway probably reflects its involvement in the degradation of recalcitrant compounds such as lignin.
Many of the pathways negatively correlated with AGS are essential pathways involved in membrane transport, DNA metabolism (transcription, translation, replication, repair and recombination), amino-acid, lipid, carbohydrate and cofactor metabolism (Additional file 13: Figure S6). This reflects the dominance of opportunistic bacteria, especially at day 3–7, with small genomes that lack the arid soil ‘specialist’ pathways discussed above.
Overall, the analysis of the functional metabolism shows a clear change between the different time-points. The “normal” microbial soil community metabolism in ENP is represented by species that have genes that are involved in dealing with stress encountered in arid soils. The carcass fluids and nutrients activate many dormant or low abundant copiothrophic lineages present in the soil, which is reflected by the pathways needed for proliferation especially present in the first week of our experiment. These two states of the soil community, dormant vs. active, are typical for arid environments where water input stimulates in a pulse-like manner microbial activity and if sufficient even plant growth . Rainfall gives enough moisture to activate both plant growth and the microbial decomposition compartment of the soil, providing plants with necessary nutrients [56, 87, 88]. However, in contrast to rainfall, animal carcasses not only provide moisture to the soil, but also additional nutrients. This pulse of nutrients and moisture allows the microbial community at very local scales to break down carcass nutrients, laying the foundation for enhanced plant growth directly and in the future [16, 56].
Our work is the first study to use shotgun metagenomics to describe the temporal changes in the soil microbial community after deposition of anthrax carcasses. We observed that the metagenomic approach is not the best way to study B. anthracis in a natural setting as it sporulates quickly making detection via DNA extraction difficult. In addition, it is too genetically similar to its close relatives such that most currently available tools to analyse (meta-) genomic sequences are not able to differentiate between such similar organisms. Nevertheless, the metagenomic data allows us to assess how this pathogen influences the microbial soil community.
The metagenomic time-series analysis indicated that carcass introduction stimulated the increase in abundance of part of the microbial soil community. This could play an indirect role in the future transmission of B. anthracis by stimulating the growth of taxa known to have PGPBs. The taxonomic shift in the microbial community was accompanied by shifts in the metabolic potential. The typical semi-arid soil community was involved in stress evasion, while the carcass introduction increased pathways involved in proliferation. Nonetheless, the time-series showed that within a month the semi-arid bacterial soil community was similar to the community at T = 0 with respect to taxonomic and metabolic composition.
Our work provides a background to study the ecology of spore-forming pathogens in a natural setting using culture-independent methods. It highlights the difficulty of using DNA based approaches to study spore-forming organisms such as B. anthracis in arid soils where they are dormant and poorly detectable.
Carcass and sample information
Sample collection in ENP was done after receiving permission from the Ministry of Environment and Tourism of Namibia (permit number: 1857/2013). Two fresh (likely <24 h since death) plains zebra (Equus quagga) anthrax carcasses (Ca1 & Ca2) were found on 03.03.2014 between Sprokieswoud and Charl Marais Dam, 835 m apart at coordinates S19.031/E015.548 (Ca1) and S19.037/E015.553 (Ca2), respectively. Information about ENP, carcass sites and sample collection can be found in the Additional file 4. Soil samples were taken from the area of blood-spill at six time-points (T), starting at T = 0 days (03.03.2014), T = 3, 7, 14, 21, 30 days. A control sample (Ctrl0) was also taken from the grids at T = 0 days, avoiding areas covered in blood (Additional file 8: Figure S4).
DNA isolation and purification
A total of 14 soil samples were collected (7 per carcass). DNA was isolated in two replicates from each sample (P1 and P2) resulting in a total of 14 DNA samples per carcass site. DNA was isolated using the FastDNA® spin kit for soil (MP Biomedicals, Santa Ana, California, USA), following the manufacturer’s protocol with adjustments specified in Additional file 4. The DNA samples were filter sterilised using an Ultrafree® Durapore PVDF 0.1 μM spinfilter (Millipore, Darmstadt, Germany).
DNA was shipped dry, following an ethanol precipitation, from ENP to Norway (see Additional file 4: for details). Upon arrival DNA was re-suspended in DES buffer from FastDNA® spin kit for soil. Samples were concentrated and purified using Agencourt® AMPure® XP beads (Beckman Coulter, Beverly, Massachusetts, USA), and treated with a PowerClean® Pro DNA clean-Up Kit, (MO BIO Laboratories, Carlsbad, California, USA) to remove any remaining inhibitors.
Metagenome sequencing and quality control
Environmental DNA samples were sequenced at the Norwegian Sequencing Centre (NSC) using Illumina MiSeq® (Illumina Inc., San Diego, California, USA). A 250 bp paired-end sequencing library (400 bp insert length) was generated using a Regular TruSeq® adapter ligation kit (Illumina Inc., San Diego, California, USA). Fastq files were processed using cut-adapt (v1.8)  and prinseq-lite (v0.20.4)  (See Additional file 4: for settings). Metagenomic read GC-content was calculated using infoseq (EMBOSS v. 6.5.7 ). The infoseq output was used to calculate average GC-content and GC-content skewness (timeData package version: 3012.100) using R-Studio (v3.3.0).
The average genome size (AGS) of each metagenome was estimated using MicrobeCensus on clean paired-end files (Additional file 4: Methods). The likelihood of sampling a universal single copy is inversely correlated with the AGS of a metagenomes [48, 92]. Thus universal single copy abundance will show significant differences when AGS between metagenomes differ.
Bacillus anthracis quantification from soils
Two methods were used to quantify B. anthracis in the soil samples; serial dilution culturing and quantitative PCR (qPCR) (see Additional file 4: for details).
The genome coverages of the B. anthracis isolate K1 and K2 genomes (Accession numbers: LBBZ00000000 and LBCA00000000 ) isolated from the Ca1 and Ca2 sites respectively, were investigated by mapping metagenomic sequences using the Burrows-Wheeler Aligner (BWA v0.7.8) . The genomes of B. cereus E33L, B. thuringiensis HD-771 and Bacillus subtilis 168 were also used as references (see Additional file 4: for details).
Taxonomic and functional classification of metagenomic reads
Taxonomic composition of the metagenomes was determined using metaxa2 (v2.0.1) which extracts SSU rRNA sequences , MEGAN (v5.10.15) , Kraken (0.10.5-beta)  and metaBIT . For the details on the comparison of the taxonomic classifiers see Additional file 4. The metaxa2 taxonomic profiles were visualised with a heatmap using ggplot2 in R-studio (see Additional file 4: for settings). Metaxa2 extracted SSU reads were run through the MetaAmp 1.1 pipeline to obtain operational taxonomical unit (OTU) abundances.
Shotgun sequences from each dataset were compared to the NCBI nr database (accessed on 20.01.2016) using Diamond (v0.7.11)  and the output was used for taxonomic and functional classification in MEGAN (v5.10.15) . Temporal species composition dynamics of metagenomic samples was examined by principle coordinate analysis (PCoA) in R-Studio (v3.3.0). Changes in the metabolic potential of the community over time were analysed by correlation of KEGG KO-terms and pathways [94, 95] abundances with AGS size (For details see Additional file 4).
The datasets supporting the conclusions of this article are available in the SRA repository, SRP076706. All sequence data was submitted to Genbank as part of bioproject: PRJNA281298.
We thank Namibia’s Ministry of Environment and Tourism for permission to conduct this research (permit number: 1857/2013). In addition, we thank Dr. Mari Espelund for help optimizing DNA extractions, Claudine C. Cloete and Zoe Barandongo for help with sample collecting and DNA extractions, the Etosha Ecological Institute and the Norwegian Defence Research Institute for providing laboratory facilities and the Department of Biosciences student funding at University of Oslo for travel money. We also thank Eric de Muinck and Pål Trosvik for critical feedback on our manuscript.
This work was supported by the Research Council of Norway (Project No. 180444/V40) awarded to CLN and NSF OISE-1103054 awarded to WCT.
Availability of data and materials
The datasets supporting the conclusions of this article are available in the SRA repository, SRP076706. All sequence data was submitted to Genbank as part of bioproject: PRJNA281298.
Additional file 4 accompanies the paper on the BMC Microbiology journal website.
All authors contributed to the design of the experiment. KV, WRE, WCT conducted the environmental sampling and qPCR experiments, bacterial cultivation and DNA extractions. KV prepared samples for metagenomic sequencing and performed bioinformatic quality control. KV, CLN and THA executed bioinformatics analysis of shotgun data. KV prepared the manuscript. All authors contributed to critical revision of the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
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.
- Fierer N, Leff JW, Adams BJ, Nielsen UN, Bates ST, Lauber CL, Owens S, Gilbert JA, Wall DH, Caporaso JG. Cross-biome metagenomic analyses of soil microbial communities and their functional attributes. Proc Natl Acad Sci U S A. 2012;109(52):21390–5.View ArticlePubMedPubMed CentralGoogle Scholar
- O'Brien SL, Gibbons SM, Owens SM, Hampton-Marcell J, Johnston ER, Jastrow JD, Gilbert JA, Meyer F, Antonopoulos DA. Spatial scale drives patterns in soil bacterial diversity. Environ Microbiol. 2016;18(6):2039–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Rousk J, Bååth E, Brookes PC, Lauber CL, Lozupone C, Caporaso JG, Knight R, Fierer N. Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J. 2010;4(10):1340–51.View ArticlePubMedGoogle Scholar
- Bell CW, Tissue DT, Loik ME, Wallenstein MD. Acosta - Martinez V, Erickson RA, Zak JC: soil microbial and nutrient responses to 7 years of seasonally altered precipitation in a Chihuahuan Desert grassland. Glob Chang Biol. 2014;20(5):1657–73.View ArticlePubMedGoogle Scholar
- Johnson SL, Kuske CR, Carney TD, Housman DC, Gallegos-Graves LV, Belnap J. Increased temperature and altered summer precipitation have differential effects on biological soil crusts in a dryland ecosystem. Glob Chang Biol. 2012;18(8):2583–93.View ArticleGoogle Scholar
- Garcia-Pichel F, Loza V, Marusenko Y, Mateo P, Potrafka RM. Temperature drives the continental-scale distribution of key microbes in topsoil communities. Science. 2013;340(6140):1574–7.View ArticlePubMedGoogle Scholar
- Niu F, He J, Zhang G, Liu X, Liu W, Dong M, Wu F, Liu Y, Ma X, An L, et al. Effects of enhanced UV-B radiation on the diversity and activity of soil microorganism of alpine meadow ecosystem in Qinghai–Tibet plateau. Ecotoxicology. 2014;23(10):1833–41.View ArticlePubMedGoogle Scholar
- Bell CW, Acosta-Martinez V, McIntyre NE, Cox S, Tissue DT, Zak JC. Linking microbial community structure and function to seasonal differences in soil moisture and temperature in a Chihuahuan desert grassland. Microb Ecol. 2009;58(4):827–42.View ArticlePubMedGoogle Scholar
- Bi J, Zhang NL, Liang Y, Yang HJ, Ma KP. Interactive effects of water and nitrogen addition on soil microbial communities in a semiarid steppe. J Plant Ecol. 2012;5(3):320–9.View ArticleGoogle Scholar
- Rath KM, Maheshwari A, Rousk J. The impact of salinity on the microbial response to drying and rewetting in soil. Soil Biol Biochem. 2017;108:17–26.View ArticleGoogle Scholar
- Van Horn DJ, Okie JG, Buelow HN, Gooseff MN, Barrett JE, Takacs-Vesbach CD. Soil microbial responses to increased moisture and organic resources along a salinity gradient in a polar desert. Appl Environ Microbiol. 2014;80(10):3034–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Howard GT, Duos B, Watson-Horzelski EJ. Characterization of the soil microbial community associated with the decomposition of a swine carcass. Int Biodeterior Biodegradation. 2010;64(4):300–4.View ArticleGoogle Scholar
- Cobaugh KL, Schaeffer SM, DeBruyn JM. Functional and structural succession of soil microbial communities below decomposing human cadavers. PLoS One. 2015;10(6):e0130201.View ArticlePubMedPubMed CentralGoogle Scholar
- Weiss S, Carter DO, Metcalf JL, Knight R. Carcass mass has little influence on the structure of gravesoil microbial communities. Int J Legal Med. 2016;130(1):253–63.View ArticlePubMedGoogle Scholar
- Bump JK, Peterson RO, Vucetich JA. Wolves modulate soil nutrient heterogeneity and foliar nitrogen by configuring the distribution of ungulate carcasses. Ecology. 2009;90(11):3159–67.View ArticlePubMedGoogle Scholar
- Turner WC, Kausrud KL, Krishnappa YS, Cromsigt JP, Ganz HH, Mapaure I, Cloete CC, Havarua Z, Küsters M, Getz WM. Fatal attraction: vegetation responses to nutrient inputs attract herbivores to infectious anthrax carcass sites. Proc R Soc Lond B Biol Sci. 2014;281(1795):20141785.View ArticleGoogle Scholar
- Sharp RJ, Roberts AG. Anthrax: the challenges for decontamination. J Chem Technol Biotechnol. 2006;81(10):1612–25.View ArticleGoogle Scholar
- Liu Y, Lai Q, Göker M, Meier-Kolthoff JP, Wang M, Sun Y, Wang L, Shao Z. Genomic insights into the taxonomic status of the Bacillus Cereus group. Sci Rep. 2015;5:14082.View ArticlePubMedPubMed CentralGoogle Scholar
- Jensen GB, Hansen BM, Eilenberg J, Mahillon J. The hidden lifestyles of Bacillus cereus and relatives. Environ Microbiol. 2003;5(8):631–40.View ArticlePubMedGoogle Scholar
- Helgason E, Økstad OA, Caugant DA, Johansen HA, Fouet A, Mock M, Hegna I, Kolstø A-B. Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis—one species on the basis of genetic evidence. Appl Environ Microbiol. 2000;66(6):2627–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Zwick ME, Joseph SJ, Didelot X, Chen PE, Bishop-Lilly KA, Stewart AC, Willner K, Nolan N, Lentz S, Thomason MK. Genomic characterization of the Bacillus cereus sensu lato species: backdrop to the evolution of Bacillus anthracis. Genome Res. 2012;22(8):1512–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Mignot T, Mock M, Robichon D, Landier A, Lereclus D, Fouet A. The incompatibility between the PlcR-and AtxA-controlled regulons may have selected a nonsense mutation in Bacillus anthracis. Mol Microbiol. 2001;42(5):1189–98.View ArticlePubMedGoogle Scholar
- Easterday WR, Ert MN, Simonson TS, Wagner DM, Kenefic LJ, Allender CJ, Keim P. Use of single nucleotide polymorphisms in the plcR gene for specific identification of Bacillus anthracis. J Clin Microbiol. 2005;43:1995–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Green BD, Battisti L, Koehler TM, Thorne CB, Ivins BE. Demonstration of a capsule plasmid in Bacillus anthracis. Infect Immun. 1985;49(2):291–7.PubMedPubMed CentralGoogle Scholar
- Mikesell P, Ivins BE, Ristroph JD, Dreier TM. Evidence for plasmid-mediated toxin production in Bacillus anthracis. Infect Immun. 1983;39(1):371–6.PubMedPubMed CentralGoogle Scholar
- Toby IT, Widmer J, Dyer DW. Divergence of protein-coding capacity and regulation in the Bacillus cereus sensu lato group. BMC Bioinformatics. 2014;15(11):1.Google Scholar
- Turner WC, Kausrud KL, Beyer W, Easterday WR, Barandongo ZR, Blaschke E, Cloete CC, Lazak J, Van Ert MN, Ganz HH. Lethal exposure: an integrated approach to pathogen transmission via environmental reservoirs. Sci Rep. 2016;6:27311.View ArticlePubMedPubMed CentralGoogle Scholar
- Dragon DC, Bader DE, Mitchell J, Woollen N. Natural dissemination of Bacillus anthracis spores in northern Canada. Appl Environ Microbiol. 2005;71(3):1610–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Bellan SE, Turnbull PC, Beyer W, Getz WM. Effects of experimental exclusion of scavengers from carcasses of anthrax-infected herbivores on Bacillus anthracis sporulation, survival, and distribution. Appl Environ Microbiol. 2013;79(12):3756–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Ganz HH, Turner WC, Brodie EL, Kusters M, Shi Y, Sibanda H, Torok T, Getz WM. Interactions between Bacillus anthracis and plants may promote anthrax transmission. PLoS Negl Trop Dis. 2014;8(6):e2903.View ArticlePubMedPubMed CentralGoogle Scholar
- Schuch R, Fischetti VA. The secret life of the anthrax agent Bacillus anthracis: bacteriophage-mediated ecological adaptations. PLoS One. 2009;4(8):e6532.View ArticlePubMedPubMed CentralGoogle Scholar
- Saile E, Koehler TM. Bacillus anthracis multiplication, persistence, and genetic exchange in the rhizosphere of grass plants. Appl Environ Microbiol. 2006;72(5):3168–74.View ArticlePubMedPubMed CentralGoogle Scholar
- Dey R, Hoffman PS, Glomski IJ. Germination and amplification of anthrax spores by soil-dwelling amoebas. Appl Environ Microbiol. 2012;78(22):8075–81.View ArticlePubMedPubMed CentralGoogle Scholar
- Logares R, Haverkamp TH, Kumar S, Lanzén A, Nederbragt AJ, Quince C, Kauserud H. Environmental microbiology through the lens of high-throughput DNA sequencing: synopsis of current platforms and bioinformatics approaches. J Microbiol Methods. 2012;91(1):106–13.View ArticlePubMedGoogle Scholar
- Be NA, Thissen JB, Gardner SN, McLoughlin KS, Fofanov VY, Koshinsky H, Ellingson SR, Brettin TS, Jackson PJ, Jaing CJ. Detection of bacillus anthracis DNA in complex soil and air samples using next-generation sequencing. PLoS One. 2013;8(9):e73455.View ArticlePubMedPubMed CentralGoogle Scholar
- Dragon DC, Rennie RP. Evaluation of spore extraction and purification methods for selective recovery of viable Bacillus anthracis spores. Lett Appl Microbiol. 2001;33(2):100–5.View ArticlePubMedGoogle Scholar
- Le Roux C, Grunow J, Morris J, Bredenkamp G, Scheepers J. A classification of the vegetation of the Etosha national park. S Afr J Bot. 1988;54(1):1–10.View ArticleGoogle Scholar
- Engert S. Spatial variability and temporal periodicity of rainfall in the Etosha National Park and surrounding areas in northern Namibia. Modoqua. 1997;20(1):115–20.Google Scholar
- Lauber CL, Metcalf JL, Keepers K, Ackermann G, Carter DO, Knight R. Vertebrate decomposition is accelerated by soil microbes. Appl Environ Microbiol. 2014;80(16):4920–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Bellan SE, Gimenez O, Choquet R, Getz WM. A hierarchical distance sampling approach to estimating mortality rates from opportunistic carcass surveillance data. Methods Ecol Evol. 2013;4(4):361–9.View ArticleGoogle Scholar
- Hyde ER, Haarmann DP, Petrosino JF, Lynne AM, Bucheli SR. Initial insights into bacterial succession during human decomposition. Int J Legal Med. 2015;129(3):661–71.View ArticlePubMedGoogle Scholar
- Lauber CL, Hamady M, Knight R, Fierer N. Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl Environ Microbiol. 2009;75(15):5111–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Albertsen M, Karst SM, Ziegler AS, Kirkegaard RH, Nielsen PH. Back to basics–the influence of DNA extraction and primer choice on phylogenetic analysis of activated sludge communities. PLoS One. 2015;10(7):e0132783.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Durbin R. Fast and accurate long-read alignment with burrows–wheeler transform. Bioinformatics. 2010;26(5):589–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Rasmussen S, Allentoft ME, Nielsen K, Orlando L, Sikora M, Sjögren K-G, Pedersen AG, Schubert M, Van Dam A, Kapel CMO. Early divergent strains of Yersinia pestis in Eurasia 5,000 years ago. Cell. 2015;163(3):571–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Valseth K, Nesbø CL, Easterday WR, Turner WC, Olsen JS, Stenseth NC, Haverkamp THA. Draft genome sequences of two bacillus anthracis strains from Etosha National Park, Namibia. Genome Announc. 2016;4(4):e00861–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Reichenberger ER, Rosen G, Hershberg U, Hershberg R. Prokaryotic nucleotide composition is shaped by both phylogeny and the environment. Genome Biol Evol. 2015;7(5):1380–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Nayfach S, Pollard KS. Average genome size estimation improves comparative metagenomics and sheds light on the functional ecology of the human microbiome. Genome Biol. 2015;16(1):51.View ArticlePubMedPubMed CentralGoogle Scholar
- Raes J, Korbel JO, Lercher MJ, von Mering C, Bork P. Prediction of effective genome size in metagenomic samples. Genome Biol. 2007;8(1):R10.View ArticlePubMedPubMed CentralGoogle Scholar
- Bengtsson-Palme J, Hartmann M, Eriksson KM, Pal C, Thorell K, Larsson DGJ, Nilsson RH. Metaxa2: improved identification and taxonomic classification of small and large subunit rRNA in metagenomic data. Mol Ecol Resour. 2015;15(6):1403–14.View ArticlePubMedGoogle Scholar
- Louvel G, Der Sarkissian C, Hanghøj K, Orlando L. metaBIT, an integrative and automated metagenomic pipeline for analysing microbial profiles from high-throughput sequencing shotgun data. Mol Ecol Resour. 2016;16(6):1415–27.View ArticlePubMedGoogle Scholar
- Truong DT, Franzosa EA, Tickle TL, Scholz M, Weingart G, Pasolli E, Tett A, Huttenhower C, Segata N. MetaPhlAn2 for enhanced metagenomic taxonomic profiling. Nat Methods. 2015;12(10):902–3.View ArticlePubMedGoogle Scholar
- Wood DE, Salzberg SL. Kraken: ultrafast metagenomic sequence classification using exact alignments. Genome Biol. 2014;15(3):R46.View ArticlePubMedPubMed CentralGoogle Scholar
- Huson DH, Auch AF, Qi J, Schuster SC. MEGAN analysis of metagenomic data. Genome Res. 2007;17(3):377–86.View ArticlePubMedPubMed CentralGoogle Scholar
- Lindgreen S, Adair KL, Gardner PP. An evaluation of the accuracy and speed of metagenome analysis tools. Sci Rep. 2016;6:19233.View ArticlePubMedPubMed CentralGoogle Scholar
- Collins SL, Sinsabaugh RL, Crenshaw C, Green L, Porras-Alfaro A, Stursova M, Zeglin LH. Pulse dynamics and microbial processes in aridland ecosystems. J Ecol. 2008;96(3):413–20.View ArticleGoogle Scholar
- Dorsch M, Lovet D, Bailey GD. Fusobacterium equinum sp. nov., from the oral cavity of horses. Int J Syst Evol Microbiol. 2001;51(6):1959–63.View ArticlePubMedGoogle Scholar
- Schlater LK, Brenner D, Steigerwalt A, Moss CW, Lambert M, Packer R. Pasteurella caballi, a new species from equine clinical specimens. J Clin Microbiol. 1989;27(10):2169–74.PubMedPubMed CentralGoogle Scholar
- Konstantinidis KT, Tiedje JM. Trends between gene content and genome size in prokaryotic species with larger genomes. Proc Natl Acad Sci U S A. 2004;101(9):3160–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Prosser JI. Dispersing misconceptions and identifying opportunities for the use of 'omics' in soil microbial ecology. Nat Rev Microbiol. 2015;13(7):439–46.View ArticlePubMedGoogle Scholar
- Teeling H, Glöckner FO. Current opportunities and challenges in microbial metagenome analysis--a bioinformatic perspective. Brief Bioinform. 2012;13(6):728–42.View ArticlePubMedPubMed CentralGoogle Scholar
- Fierer N, Bradford MA, Jackson RB. Towards an ecological classification of soil bacteria. Ecology. 2007;88(6):1354–64.View ArticlePubMedGoogle Scholar
- Tan J, Zuniga C, Zengler K. Unraveling interactions in microbial communities - from co-cultures to microbiomes. J Microbiol. 2015;53(5):295–305.View ArticlePubMedGoogle Scholar
- Staley JT, Konopka A. Measurement of in situ activities of nonphotosynthetic microorganisms in aquatic and terrestrial habitats. Annu Rev Microbiol. 1985;39(1):321–46.View ArticlePubMedGoogle Scholar
- Torsvik V, Øvreås L. Microbial diversity and function in soil: from genes to ecosystems. Curr Opin Microbiol. 2002;5(3):240–5.View ArticlePubMedGoogle Scholar
- Lindeque PM, Turnbull PCB. Ecology and epidemiology of anthrax in Etosha National Park, Namibia. Onderstepoort J Vet Res. 1994;61(1):71–83.PubMedGoogle Scholar
- Sedlackova V, Dziedzinska R, Babak V, Kralik P. The detection and quantification of Bacillus thuringiensis spores from soil and swabs using quantitative PCR as a model system for routine diagnostics of Bacillus anthracis. J Appl Microbiol. 2017; https://doi.org/10.1111/jam.13445.
- Ketola T, Mikonranta L, Laakso J, Mappes J. Different food sources elicit fast changes to bacterial virulence. Biol Lett. 2016;12(1):20150660.View ArticlePubMedPubMed CentralGoogle Scholar
- Erlacher A, Cernava T, Cardinale M, Soh J, Sensen CW, Grube M, Berg G. Rhizobiales as functional and endosymbiontic members in the lichen symbiosis of Lobaria pulmonaria L. Front Microbiol. 2015;6:53.View ArticlePubMedPubMed CentralGoogle Scholar
- Goodfellow M, Kämpfer P, Busse H-J, Trujillo ME, Suzuki K-I, Ludwig W, Whitman WB. Bergey's manual of systematic bacteriology : Volume 5: The Actinobacteria. In: vol. 5, 2 edn. New York: Springer Science & Business Media; 2012: 509–546.Google Scholar
- Foesel B, Geppert A, Rohde M, Overmann J. Parviterribacter kavangonensis and Parviterribacter multiflagellatus a novel genus and two novel species within the order Solirubrobacterales and emended description of the classes Thermoleophilia and Rubrobacteria and its orders and families. Int J Syst Evol Microbiol. 2016;66:652–65.View ArticleGoogle Scholar
- Weber BS, Harding CM, Feldman MF. Pathogenic Acinetobacter: from the cell surface to infinity and beyond. J Bacteriol. 2015;198(6):880–7.View ArticlePubMedGoogle Scholar
- Berry C. The bacterium, Lysinibacillus sphaericus, as an insect pathogen. J Invertebr Pathol. 2012;109(1):1–10.View ArticlePubMedGoogle Scholar
- Vos P, Garrity G, Jones D, Krieg NR, Ludwig W, Rainey FA, Schleifer K-H, Whitman W. Bergey's manual of systematic bacteriology: Volume 3: The Firmicutes. In: vol. 3, 2 edn. New York: Springer Science & Business Media; 2011: 364–370.Google Scholar
- Andrade LF, de Souza G, Nietsche S, Xavier AA, Costa MR, Cardoso AMS, Pereira MCT, Pereira D. Analysis of the abilities of endophytic bacteria associated with banana tree roots to promote plant growth. J Microbiol. 2014;52(1):27–34.View ArticlePubMedGoogle Scholar
- Batool F, Rehman Y, Hasnain S. Phylloplane associated plant bacteria of commercially superior wheat varieties exhibit superior plant growth promoting abilities. Front Life Sci. 2016;9(4):313–22.View ArticleGoogle Scholar
- Uroz S, Buée M, Murat C, Frey-Klett P, Martin F. Pyrosequencing reveals a contrasted bacterial diversity between oak rhizosphere and surrounding soil. Environ Microbiol Rep. 2010;2(2):281–8.View ArticlePubMedGoogle Scholar
- Vandenkoornhuyse P, Mahé S, Ineson P, Staddon P, Ostle N, Cliquet JB, Francez AJ, Fitter AH, Young JP. Active root-inhabiting microbes identified by rapid incorporation of plant-derived carbon into RNA. Proc Natl Acad Sci U S A. 2007;104(43):16970–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Bulgarelli D, Schlaeppi K, Spaepen S, Ver Loren van Themaat E, Schulze-Lefert P. Structure and functions of the bacterial microbiota of plants. Annu Rev Plant Biol. 2013;64:807–38.View ArticlePubMedGoogle Scholar
- Baldani JI, Rouws L, Cruz LM, Olivares FL, Schmid M, Hartmann A. The family Oxalobacteraceae. In: The prokaryotes: Alphaproteobacteria and Betaproteobacteria. Vol. 9783642301971; 2014. p. 919–74.View ArticleGoogle Scholar
- Pieterse CM, Zamioudis C, Berendsen RL, Weller DM, Van Wees SC, Bakker PA. Induced systemic resistance by beneficial microbes. Annu Rev Phytopathol. 2014;52:347–75.View ArticlePubMedGoogle Scholar
- Choudhary DK, Johri BN. Interactions of Bacillus spp. and plants--with special reference to induced systemic resistance (ISR). Microbiol Res. 2009;164(5):493–513.View ArticlePubMedGoogle Scholar
- Potts M. Desiccation tolerance of prokaryotes. Microbiol Rev. 1994;58(4):755–805.PubMedPubMed CentralGoogle Scholar
- Shuman S, Glickman MS. Bacterial DNA repair by non-homologous end joining. Nat Rev Microbiol. 2007;5(11):852–61.View ArticlePubMedGoogle Scholar
- Moeller R, Stackebrandt E, Reitz G, Berger T, Rettberg P, Doherty AJ, Horneck G, Nicholson WL. Role of DNA repair by nonhomologous-end joining in Bacillus subtilis spore resistance to extreme dryness, mono-and polychromatic UV, and ionizing radiation. J Bacteriol. 2007;189(8):3306–11.View ArticlePubMedPubMed CentralGoogle Scholar
- Garcia-Gonzalez A, Vicens L, Alicea M, Massey S. The distribution of recombination repair genes is linked to information content in bacteria. Gene. 2013;528(2):295–303.View ArticlePubMedGoogle Scholar
- Du Plessis W. Effective rainfall defined using measurements of grass growth in the Etosha National Park, Namibia. J Arid Environ. 2001;48(3):397–417.View ArticleGoogle Scholar
- Clark JS, Campbell JH, Grizzle H, Acosta-Martìnez V, Zak JC. Soil microbial community response to drought and precipitation variability in the Chihuahuan Desert. Microb Ecol. 2009;57(2):248–60.View ArticlePubMedGoogle Scholar
- Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet journal. 2011;17(1):10–2.View ArticleGoogle Scholar
- Schmieder R, Edwards R. Quality control and preprocessing of metagenomic datasets. Bioinformatics. 2011;27(6):863–4.View ArticlePubMedPubMed CentralGoogle Scholar
- Rice P, Longden I, Bleasby A. EMBOSS: the European molecular biology open software suite. Trends Genet. 2000;16(6):276–7.View ArticlePubMedGoogle Scholar
- Beszteri B, Temperton B, Frickenhaus S, Giovannoni SJ. Average genome size: a potential source of bias in comparative metagenomics. ISME J. 2010;4(8):1075–7.View ArticlePubMedGoogle Scholar
- Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12(1):59–60.View ArticlePubMedGoogle Scholar
- Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016;44(D1):D457–62.View ArticlePubMedGoogle Scholar
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.View ArticlePubMedPubMed CentralGoogle Scholar