Ecology

Nutrient, FOM, and physiochemical contents of the fecal manures
For the three main fertility nutrients (N, P and K) in the fecal manures (Table 1a) the E manure contained the highest N level at 80 µg/g, 16-fold more than the other samples. In addition, it had a high P level of 100 µg/g (along with PC1 and PC3) compared with an average of 34.66 µg/g for the others, and a moderate K level of 150 µg/g (others ranged from 125.00 to 350.00 µg/g. The FOM content ranged from 1.00 to 1.50% for all the animal manures (Table 1a). For the general physiochemical properties, chicken manures showed a relatively high salinity (Table 1b: 3.52 ± 0.28 to 4.01 ± 0.64 ppt), and these ionic salts caused a high conductivity (11.27 ± 0.69 to 12.79 ± 0.08 mS/cm) and high water saturation (73.29 ± 0.40 to 80.56 ± 3.17 mS/cm), giving less available water for plant roots. A moderate level of salinity, conductivity and water content, such as those in E, are generally preferred for farming soil. The water content of the manure was lowest in G2, goats fed with P. purpureum (23.75 ± 3.07%), but increased two-fold in G1 goats fed with D. eriantha (49.35 ± 6.08%), suggesting a different type of grass feed might confer different fecal liquid and dry contents, or the chemical composition influenced the water activity in the goat’s gut. The R manure also had a low water content (41.53 ± 5.57%).
Table 1 Characteristics of (a) NPK and FOM, and (b) physiochemical levels of the different animal fecal manures.
Full size table

The visual observation of each manure was in accord with the reported water content, with a dry and hard texture (insufficient moisture) for G1, G2, and R manures, and a wet texture for the chicken, cattle, pig, and buffalo manures. Too much water in the soil can adversely affect aeration, where the plant’s roots receive insufficient oxygen and rot. Moreover, most plants prefer a slightly acidic soil pH of 6.2–6.8, which was found in the R, MB and E composts (Table 1b). The manure pH can change through reactions such as organic matter decomposition that affects the availabilities of NPK39.
From the overall comparison of the NPK nutrients and physiochemical properties, the fecal manures from the Phuparn chicken species might be less suitable as a fertilizer because of the poor N level, and high level of salts and conductivity. Indeed, the manures from most species contained a low N level, and those from goat and rabbit also had a relatively low water content (insufficient water for a plant to absorb nutrients from the soil through the roots, to the trunk, leaves and fruits). A moderate water content in the manure was appropriate. But of course, a somewhat low water level could be adjusted by adding water. The biophysiochemical quality analyses of the manures suggested that the E manure would be appropriate as a fertilizer because of the enhanced NPK levels, low salinity, low conductivity, slightly acidic-neutral pH, and a moderate water content (Table 1b).
Quantification of total bacteria in the fecal manures
The total number of bacteria in copies/g DW fecal manure was derived from the 16S rRNA gene qPCR. Figure 1a showed the data were consistent between independent triplicate samples, and that the E and R manures had the lowest bacterial counts (at 2.1 × 109 and 2.5 × 109 copies/g DW, respectively). The other animal fecal manures contained more than 5 × 109 copies/g DW, and the greatest bacterial load was found in PC2 and MB at 2.6 × 1010 and 1.9 × 1010 copies/g DW, respectively, while PC1 and PC3 had a significantly lower level (p  98% of the sequencing coverage of taxonomic compositions at the genus and species levels, respectively (Table 2). The average Good’s coverage indices were 99.54% and 99.43% for the genus and species levels, respectively (Supplemental Table 1). This was consistent with the plateau rarefaction curves, showing the frequencies of OTUs become constant despite increasing sequencing reads, meaning a sufficient sequencing coverage was obtained (Supplemental Fig. 1).
Table 2 Good’s coverage indices (estimated sequencing coverage) and alpha diversity indices of bacterial taxonomic profiles by 16S rRNA gene sequences at the (a) genus and (b) species levels.
Full size table

Combining the count of total 16S rRNA gene copies with the OTU percent compositions gave the quantitative number of copies of each OTU in a community. These quantitative microbiota data were then used to compute the alpha and beta diversity measurements. The alpha diversity revealed that the microbiota of E had the relatively most diverse OTUs at both the genus and species levels (Chao richness, Fig. 2a,c), regardless of having the lowest total bacterial count (Fig. 1). The low total bacterial abundance but high diversity in E underlined that the various OTUs of bacteria in E might be present in small numbers when compared to the copy numbers of OTUs in the other animal manures. Note that the alpha diversity obtained by considering the distribution of general OTUs among the different animal manures were similar (Shannon diversity, Fig. 2b,d).
Figure 2

Alpha diversity measurements of OTU compositions at (a,b) genus and (c,d) species levels, by richness (Chao) and evenness (Shannon). Box plot with bar representing the mean from three sequencing replicates, and asterisk (*) indicates a statistically significant difference by one-way ANOVA at p  1% abundance). In (b), the OTUs where Mothur could not identify the genus name were denoted by small letters (p_ abbreviates phylum; o_, order; and f_, family) to the deepest taxonomic names that could be identified.

Full size image

Relative frequencies of plant symbiotic and pathogenic bacterial genera
Plant symbiosis and pathogenic bacteria were analysed across the different animal manures. Manures PC1-3, G1, G2 and CC demonstrated generally abundant symbiotic bacteria comprised of the genera Pseudomonas, Bacillus, Arthrobacter, Flavobacterium, Alcaligenes and Streptomyces (Fig. 4a). For examples, PC2 contained abundant Alcaligenes (1.16 × 1010 cells/g DW), Pseudomonas (3.35 × 109), Flavobacterium (2.82 × 108) and Arthrobacter (3.8 × 108). Streptomyces was only found in G1, G2 and E in moderate numbers (5.9–29.0 × 106 cells/g DW).
Figure 4

Comparative number of bacterial genera categorized as plant (a) symbiotic or (b) pathogenic genera across the different animal manures. Data represented mean ± S.D. List of bacteria categorized as plant symbionts and pathogens were downloaded from the Virulence Factor Database (VFDB).

Full size image

For the plant pathogenic bacteria, high levels were found in the manures that also contained high levels of plant symbiosis bacteria (PC1-3, G1, G2 and CC), except for E, plus the other manures, such as S, D and MB. The E had none of VFDB-listed plant pathogens. Thus, all manures except for E contained plant pathogenic bacteria, and these were from Escherichia, Shigella, Enterococcus, Clostridium, Acinetobacter, Treponema, Staphylococcus and Bacteroides (Fig. 4b). This finding correlated with the percent abundance of genera (Fig. 3b) where, examples, Treponema were relatively high in PCO and PC3 at 2.03 × 108 and 1.31 × 108 cells/g DW, respectively. Escherichia were relatively low in R and DC at 1.15 × 105 and 2.06 × 105 cells/g DW, respectively. In contrast, only E did not contain any VFDB-listed plant pathogenic bacteria suggesting that E offers a plant pathogen-free organic fertilizer.
Relationship among bacterial communities, and statistical correlation with the biophysiochemical properties
The NMDS demonstrated both the reproducibility of the data between the independent triplicate samples, except for R and PCO. The larger variation in bacterial communities were found generally between the manures from different animal species, while the minor variation in bacterial communities were found between the breeds, or the feeding diets. For examples, the variation in quantitative microbiota profiles between G1, G2 and E, compared with PC1-3 (p = 0.085) and MB (p = 0.59) (Supplemental Fig. 2). Indeed, the quantitative microbiota structures belonging to the six animal manures that had all the VFDB’s categorized plant pathogens (PC, S, R, D and MB, except PCO) were rather distant from E (p = 0.101, 0.052, 0.089, 0.104 and 0.099, respectively).
Seven parameters (NPK and four physiochemical properties) were analyzed for possible Pearson’s correlation with any of the quantitative microbiota structures. The N level was strongly correlated (p = 0.01) to the structures and in the same direction of the E (Fig. 5). The percent water, salinity and conductivity characteristics were significant and associated in the direction opposite to E, and also to the G1, G2 and PCO microbiota structures. On the other hand, the PC1, PC2 and MB microbiota structures were strongly associated with the water content, salinity and conductivity. These parameters are suggested to be important in controlling the diversity of microbiota structures.
Figure 5

Relationship among bacterial communities via NMDS constructed from thetayc distance coefficients among quantitative microbiota (stress value = 0.15, R2 = 0.86), and Pearson’s correlation with nutrients and physiochemical properties (AMOVA, p  More

1.
Würsig, B. & Clark, C. W. Behavior. in The bowhead whale (eds. Burns, J. J., Montague, J. & Cowles, C. J.) 157–199 (the Society for Marine Mammalogy, 1993).
2.
Stafford, K. M., Moore, S. E., Laidre, K. L. & Heide-Jørgensen, M. P. Bowhead whale springtime song off West Greenland. J. Acoust. Soc. Am. 124, 3315–3323 (2008).
PubMed  Article  ADS  Google Scholar 

3.
Delarue, J., Laurinolli, M. & Martin, B. Bowhead whale (Balaena mysticetus) songs in the Chukchi Sea between October 2007 and May 2008. J. Acoust. Soc. Am. 126, 3319–3328 (2009).
PubMed  Article  ADS  Google Scholar 

4.
Tervo, O. M., Christoffersen, M. F., Parks, S. E., Kristensen, R. M. & Madsen, P. T. Evidence for simultaneous sound production in the bowhead whale (Balaena mysticetus). J. Acoust. Soc. Am. 130, 2257–2262 (2011).
PubMed  Article  ADS  Google Scholar 

5.
Tervo, O. M., Parks, S. E., Christoffersen, M. F., Miller, L. A. & Kristensen, R. M. Annual changes in the winter song of bowhead whales (Balaena mysticetus) in Disko Bay, Western Greenland . Mar. Mammal Sci. 27, E241–E252 (2011).
Article  Google Scholar 

6.
Payne, R. S. & McVay, S. Songs of humpback whales. Science 173, 585–597 (1971).
CAS  PubMed  Article  ADS  Google Scholar 

7.
Collins, S. Vocal fighting and flirting: the functions of birdsong. in Nature’s Music 39–79 (Elsevier, Amsterdam, 2004).

8.
Johnson, H. D., Stafford, K. M., George, J. C., Ambrose, W. G. & Clark, C. W. Song sharing and diversity in the Bering-Chukchi-Beaufort population of bowhead whales (Balaena mysticetus), spring 2011. Mar. Mammal Sci. 31, 902–922 (2015).
Article  Google Scholar 

9.
Cholewiak, D. M., Cerchio, S., Jacobsen, J. K., Urbán-R, J. & Clark, C. W. Songbird dynamics under the sea: acoustic interactions between humpback whales suggest song mediates male interactions. R. Soc. Open Sci. 5, 171298 (2018).
PubMed  PubMed Central  Article  ADS  Google Scholar 

10.
Murray, A., Dunlop, R. A., Noad, M. J. & Goldizen, A. W. Stereotypic and complex phrase types provide structural evidence for a multi-message display in humpback whales (Megaptera novaeangliae). J. Acoust. Soc. Am. 143, 980–994 (2018).
PubMed  Article  ADS  Google Scholar 

11.
Garland, E. C. et al. Humpback whale song on the Southern Ocean feeding grounds: Implications for cultural transmission. PLoS ONE 8, e79422 (2013).
PubMed  PubMed Central  Article  ADS  CAS  Google Scholar 

12.
Rekdahl, M. L. et al. Culturally transmitted song exchange between humpback whales (Megaptera novaeangliae) in the southeast Atlantic and southwest Indian Ocean basins. 15 (2018).

13.
Darling, J. D., Acebes, J. M. V., Frey, O., Jorge Urbán, R. & Yamaguchi, M. Convergence and divergence of songs suggests ongoing, but annually variable, mixing of humpback whale populations throughout the North Pacific. Sci. Rep. 9, 7002 (2019).
PubMed  PubMed Central  Article  ADS  CAS  Google Scholar 

14.
Garland, E. C. et al. Quantifying humpback whale song sequences to understand the dynamics of song exchange at the ocean basin scale. J. Acoust. Soc. Am. 133, 560–569 (2013).
PubMed  Article  ADS  Google Scholar 

15.
Garland, E. C. et al. The devil is in the detail: Quantifying vocal variation in a complex, multi-levelled, and rapidly evolving display. J. Acoust. Soc. Am. 142, 460–472 (2017).
PubMed  Article  ADS  Google Scholar 

16.
Allen, J. A., Garland, E. C., Dunlop, R. A. & Noad, M. J. Cultural revolutions reduce complexity in the songs of humpback whales. Proc. R. Soc. B Biol. Sci. 285, 20182088 (2018).
Article  Google Scholar 

17.
Tervo, O. M., Parks, S. E. & Miller, L. A. Seasonal changes in the vocal behavior of bowhead whales (Balaena mysticetus) in Disko Bay, Western-Greenland. J. Acoust. Soc. Am. 126, 1570–1580 (2009).
PubMed  Article  ADS  Google Scholar 

18.
Payne, K. & Payne, R. Large scale changes over 19 years in songs of humpback whales in Bermuda. Ethology 68, 89–114 (1985).
Google Scholar 

19.
Cholewiak, D. M., Sousa-Lima, R. S. & Cerchio, S. Humpback whale song hierarchical structure: historical context and discussion of current classification issues. Mar. Mammal Sci. 29, E312–E332 (2013).
Article  Google Scholar 

20.
Kowarski, K., Moors-Murphy, H., Maxner, E. & Cerchio, S. Western North Atlantic humpback whale fall and spring acoustic repertoire: Insight into onset and cessation of singing behavior. J. Acoust. Soc. Am. 145, 2305–2316 (2019).
PubMed  Article  ADS  PubMed Central  Google Scholar 

21.
Stafford, K. M. et al. Spitsbergen’s endangered bowhead whales sing through the polar night. Endanger. Spec. Res. 18, 95–103 (2012).
Article  Google Scholar 

22.
Stafford, K. M., Lydersen, C., Wiig, Ø. & Kovacs, K. M. Extreme diversity in the songs of Spitsbergen’s bowhead whales. Biol. Lett. 14, 20180056 (2018).

23.
Boertmann, D., Kyhn, L. A., Witting, L. & Heide-Jørgensen, M. P. A hidden getaway for bowhead whales in the Greenland Sea. Polar Biol. 38, 1315–1319 (2015).
Article  Google Scholar 

24.
André, M. et al. Listening to the deep: live monitoring of ocean noise and cetacean acoustic signals. Mar. Pollut. Bull. 63, 18–26 (2011).
PubMed  Article  CAS  PubMed Central  Google Scholar 

25.
Bioacoustics Research Program. Raven Pro: Interactive sound analysis software (Version 1.5). (The Cornell Lab of Ornithology, 2014).

26.
Wilden, I., Herzel, H., Peters, G. & Tembrock, G. Subharmonics, biphonation, and deterministic chaos in mammal vocalization. Bioacoustics 9, 171–196 (1998).
Article  Google Scholar 

27.
Bonett, D. G. Confidence interval for a coefficient of quartile variation. Comput. Stat. Data Anal. 50, 2953–2957 (2006).
MathSciNet  MATH  Article  Google Scholar 

28.
Demšar, J. et al. Orange: data mining toolbox in python. J. Mach. Learn. Res. 14, 2349–2353 (2013).
MATH  Google Scholar 

29.
Kruskal, J. B. Multidimensional scaling by optimizing goodness of fit to a nonmetric hypothesis. Psychometrika 29, 1–27 (1964).
MathSciNet  MATH  Article  Google Scholar 

30.
Hair, J. F., Black, W. C., Babin, B. J., Anderson, R. E. & Tatham, R. L. Multidimensional scaling. in Multivariate data analysis: A global perspective 539–582 (Upper Saddle River, NJ: Pearson, 2010).

31.
Rousseeuw, P. J. Silhouettes: a graphical aid to the interpretation and validation of cluster analysis. J. Comput. Appl. Math. 20, 53–65 (1987).
MATH  Article  Google Scholar 

32.
Arbelaitz, O., Gurrutxaga, I., Muguerza, J., Pérez, J. M. & Perona, I. An extensive comparative study of cluster validity indices. Pattern Recognit. 46, 243–256 (2013).
Article  Google Scholar 

33.
Ahonen, H. et al. The underwater soundscape in western Fram Strait: breeding ground of Spitsbergen’s endangered bowhead whales. Mar. Pollut. Bull. 123, 97–112 (2017).
CAS  PubMed  Article  Google Scholar 

34.
Ljungblad, D. K., Thompson, P. O. & Moore, S. E. Underwater sounds recorded from migrating bowhead whales, Balaena mysticetus, in 1979. J. Acoust. Soc. Am. 71, 477–482 (1982).
Article  ADS  Google Scholar 

35.
Kroodsma, D. E. Reproductive development in a female songbird: differential stimulation by quality of male song. Science 192, 574–575 (1976).
CAS  PubMed  Article  ADS  Google Scholar 

36.
Byers, B. E. & Kroodsma, D. E. Female mate choice and songbird song repertoires. Anim. Behav. 77, 13–22 (2009).
Article  Google Scholar 

37.
Payne, K. The progressively changing songs of humpback whales: A window on the creative process in a wild animal. in The origins of music 135–150 (2000).

38.
Parsons, E. C. M., Wright, A. J. & Gore, M. A. The nature of humpback whale (Megaptera novaeangliae) song. J. Mar. Anim. Their Ecol. 1, 22–31 (2008).
Google Scholar 

39.
Van Parijs, S. M., Lydersen, C. & Kovacs, K. M. Vocalizations and movements suggest alternative mating tactics in male bearded seals. Anim. Behav. 65, 273–283 (2003).
Article  Google Scholar 

40.
Ballentine, B., Hyman, J. & Nowicki, S. Vocal performance influences female response to male bird song: an experimental test. Behav. Ecol. 15, 163–168 (2004).
Article  Google Scholar 

41.
Podos, J., Huber, S. K. & Taft, B. Bird song: the interface of evolution and mechanism. Annu. Rev. Ecol. Evol. Syst. 35, 55–87 (2004).
Article  Google Scholar 

42.
Podos, J. A performance constraint on the evolution of trilled vocalizations in a songbird family (Passeriformes: Emberizidae). Evolution 51, 537–551 (1997).
PubMed  Article  Google Scholar 

43.
Moseley, D. L., Lahti, D. C. & Podos, J. Responses to song playback vary with the vocal performance of both signal senders and receivers. Proc. R. Soc. B Biol. Sci. 280, 20131401 (2013).
Article  Google Scholar 

44.
Reby, D. et al. Evidence of biphonation and source–filter interactions in the bugles of male North American wapiti (Cervus canadensis). J. Exp. Biol. 219, 1224–1236 (2016).
CAS  PubMed  Article  Google Scholar 

45.
Aubin, T., Jouventin, P. & Hildebrand, C. Penguins use the two-voice system to recognize each other. Proc. R. Soc. B Biol. Sci. 267, 1081–1087 (2000).
CAS  Article  Google Scholar 

46.
Fitch, W. T., Neubauer, J. & Herzel, H. Calls out of chaos: the adaptive significance of nonlinear phenomena in mammalian vocal production. Anim. Behav. 63, 407–418 (2002).
Article  Google Scholar 

47.
Volodina, E. V., Volodin, I. A., Isaeva, I. V. & Unck, C. Biphonation may function to enhance individual recognition in the dhole, Cuon alpinus. Ethology 112, 815–825 (2006).
Article  Google Scholar 

48.
Volodin, I. A. & Nagaylik, M. M. Cues to orientation of a caller to a listener in biphonic and non-biphonic close range contact calls in the dhole (Cuon alpinus). 11 (2006).

49.
Miller, P. J. O., Samarra, F. I. P. & Perthuison, A. D. Caller sex and orientation influence spectral characteristics of “two-voice” stereotyped calls produced by free-ranging killer whales. J. Acoust. Soc. Am. 121, 3932 (2007).
PubMed  Article  ADS  Google Scholar 

50.
Schoenfuss, H. L. et al. The anatomy of the larynx of the bowhead whale, Balaena mysticetus, and its sound-producing functions. Anat. Rec. 297, 1316–1330 (2014).
Article  Google Scholar 

51.
Root-Gutteridge, H. et al. A lifetime of changing calls: North Atlantic right whales, Eubalaena glacialis, refine call production as they age. Anim. Behav. 137, 21–34 (2018).
Article  Google Scholar 

52.
Mercado, E., Herman, L. M. & Pack, A. A. Song copying by humpback whales: themes and variations. Anim. Cogn. 8, 93–102 (2005).
PubMed  Article  Google Scholar 

53.
Zoloth, S. & Green, S. Monkey vocalizations and human speech: Parallels in perception?. Brain. Behav. Evol. 16, 430–442 (1979).
CAS  PubMed  Article  Google Scholar 

54.
Soltis, J., Leong, K. & Savage, A. African elephant vocal communication II: rumble variation reflects the individual identity and emotional state of callers. Anim. Behav. 70, 589–599 (2005).
Article  Google Scholar 

55.
Brady, B., Hedwig, D., Trygonis, V. & Gerstein, E. Classification of Florida manatee (Trichechus manatus latirostris ) vocalizations. J. Acoust. Soc. Am. 147, 1597–1606 (2020).
PubMed  Article  ADS  Google Scholar 

56.
Murray, S. O., Mercado, E. & Roitblat, H. L. Characterizing the graded structure of false killer whale (Pseudorca crassidens ) vocalizations. J. Acoust. Soc. Am. 104, 1679–1688 (1998).
CAS  PubMed  Article  ADS  Google Scholar 

57.
Karlsen, J. D., Bisther, A., Lydersen, C., Haug, T. & Kovacs, K. M. Summer vocalisations of adult male white whales (Delphinapterus leucas) in Svalbard, Norway. 10 (2002).

58.
Nemiroff, L. & Whitehead, H. Structural characteristics of pulsed calls of long-finned pilot whales Globicephala melas. Bioacoustics 19, 67–92 (2009).
Article  Google Scholar 

59.
Dunlop, R. A., Noad, M. J., Cato, D. H. & Stokes, D. The social vocalization repertoire of east Australian migrating humpback whales (Megaptera novaeangliae). J. Acoust. Soc. Am. 122, 2893 (2007).
PubMed  Article  ADS  Google Scholar 

60.
Clark, C. W. The acoustic repertoire of the Southern right whale, a quantitative analysis. Anim. Behav. 30, 1060–1071 (1982).
Article  Google Scholar 

61.
Trygonis, V., Gerstein, E., Moir, J. & McCulloch, S. Vocalization characteristics of North Atlantic right whale surface active groups in the calving habitat, southeastern United States. J. Acoust. Soc. Am. 134, 4518–4531 (2013).
PubMed  Article  ADS  PubMed Central  Google Scholar 

62.
Janik, V. M. Pitfalls in the categorization of behaviour: a comparison of dolphin whistle classification methods. Anim. Behav. 57, 133–143 (1999).
CAS  PubMed  Article  PubMed Central  Google Scholar 

63.
Tougaard, J. & Eriksen, N. Analysing differences among animal songs quantitatively by means of the Levenshtein distance measure. Behaviour 143, 239–252 (2006).
Article  Google Scholar 

64.
Garland, E. C. et al. Improved versions of the Levenshtein distance method for comparing sequence information in animals’ vocalisations: tests using humpback whale song. Behaviour 149, 1413–1441 (2012).
Article  Google Scholar 

65.
Garland, E. C. et al. Population structure of humpback whales in the western and central South Pacific Ocean as determined by vocal exchange among populations. Conserv. Biol. 29, 1198–1207 (2015).
PubMed  Article  PubMed Central  Google Scholar 

66.
Koski, W. R., Davis, R. A., Miller, G. W. & Withrow, D. E. Reproduction. in The bowhead whale (eds. Burns, J. J., Montague, J. J. & Cowles, C. J.) 239–274 (The Society for Marine Mammalogy, 1993).

67.
Heide-Jørgensen, M. P., Laidre, K. L., Quakenbush, L. T. & Citta, J. J. The Northwest Passage opens for bowhead whales. Biol. Lett. 8, 270–273 (2012).
PubMed  Article  PubMed Central  Google Scholar 

68.
Laxon, S. W. et al. CryoSat-2 estimates of Arctic sea ice thickness and volume. Geophys. Res. Lett. 40, 732–737 (2013).
Article  ADS  Google Scholar 

69.
Stafford, K. M., Lydersen, C., Wiig, ø. & Kovacs, K. M. Data from: Extreme diversity in the songs of Spitsbergen’s bowhead whales. Dryad, Dataset. https://doi.org/https://doi.org/10.5061/dryad.1ck400f (2018). More

Percoll density gradient centrifugation
The removal of the associated bacteria from KMHK cells against four Percoll density gradients, 90%, 90–50%, 90–50–30% and 90–50–30–10% were shown in Fig. 2. There was no significant difference of the remaining total bacterial counts (with an initial total bacterial count of 6.36 ± 0.04 log10 CFU/mL) between algal samples centrifuged with 90% and 90%–50% Percoll gradients, but decreased significantly from 5.02 ± 0.2 log10 CFU/mL to 4.38 ± 0.05 log10 CFU/mL when dinoflagellate samples were centrifuged with 90–50–30% Percoll gradient with no further increase on adding another layer of 10% density to the gradient (i.e., 90–50–30–10%). These suggested that the highest bacterial removal capacity was achieved by centrifugation of the KMHK cells with the three-layer discontinuous (90–50–30%) gradient. This gradient was adopted in subsequent experiments in the present study, but it was different from Cho et al. who centrifuged the small Haptphyta, Isochrysis galbana (6–12 μm) with the five-layer discontinuous gradient (50%–40%–30%–20%–10%) and harvested the algal cell in between 40 and 30% Percoll5,27. As far as we know, this is the only previous study employing discontinuous gradient for algal culture, and it is obvious the optimized gradient composition varies among algal species, probably because of the diverse algal size and morphology. Vu et al. (2018) suggested that cells with a swimming ability may swim away from the concentrated zone after centrifugation, resulting in low cell recovery efficiency11. In this study, however, the swimming ability of KMHK cells was lost only temporarily for several minutes after centrifugation. This indicated that KMHK cell recovery would not be affected if the supernatant were removed immediately after centrifugation.
Figure 2

Total bacterial count in the algal sample after centrifugation with different Percoll density gradients. All data are presented as means ± standard deviations of three independent experiments (n = 3). Different letters on the top of the bar indicate that the means were significantly different among gradients at p ≤ 0.05 according to one-way analysis of variance followed by Tukey multiple comparison tests.

Full size image

Density gradient centrifugation provides an excellent alternative to filtration and micro-pipetting to physically separate bacterial cells from the algal cells, especially for fragile cell separation11. The density medium provides padding and thus protects the algal cells from the shearing force and enhance the separation efficiency17. Although filtration is one of the commonly adopted separating techniques because of its convenience, inexpensive and easy to use11, it is infeasible for use with dinoflagellate samples. However, membrane filters can be easily clogged by the algal cells and thus prolong the processing time of the filtration. Our previous experience found that more than 1 h was required to filter 10 mL of a dinoflagellate sample with 1 L of sterile medium for washing. More, algal cell recovery becomes extremely difficult when the cells are stuck firmly to the filter, and the algal cells may be easily damaged during cell harvest. This is a very important consideration for fragile cells of the unarmored dinoflagellates as they are vulnerable to the shear force generated during filtration. Micropipetting method applied in Coscinodiscus wailesii was laborious, and a skilled operator was required to pipette the algal cells from a culture droplet and sequential wash in droplets19.
In the gradient centrifugation, the commonly used density media include Ficoll, Ludox and Percoll. Of these, Ficoll is not suitable for marine samples because it is a polysaccharide-based medium that is non-isotonic at high medium concentrations and becomes viscous liquid when dissolved in seawater28,29. Ludox and Percoll are both silicon-based density media compatible with seawater. Percoll is preferred for cell separation because it is less cytotoxic to the algal cells than is Ludox30. Few studies have demonstrated the feasibility of using Percoll density gradient centrifugation for microalgal samples17,31. Two fragile dinoflagellate species, Heterosigma carterae and Cochlodinium polykrikoides, have been efficiently harvested and recovered through centrifugation using 90% Percoll17,31.
Bacterial removal by basic and extended protocols
In the present study, Percoll density gradient centrifugation was coupled with antibiotic treatment. It is because the use of antibiotics is one of the most common bacterial killing methods but antibiotics alone rarely achieve complete bacterial elimination from the algal culture11. Antibiotic susceptibility testing results in this study reveal that the bacteria associated with KMHK culture were sensitive to the antibiotic cocktail used, that is, a combination of 100 U of penicillin, 100 µg/mL streptomycin, 100 µg/mL gentamicin and 1 µg/mL tetracycline. Similarly, Ki and Han also reported that the combination of 100 mg/L streptomycin, 150 mg/L ampicillin, 150 mg/L penicillin G and 200 mg/L gentamicin effectively killed bacteria without having detrimental effects on the dinoflagellates Peridinium bipes and A. tamarense12. Guillard demonstrated that most algal species tolerated 100 mg/L penicillin, 25 mg/L streptomycin and 25 mg/L gentamicin reasonably well32. However, many red colonies, identified as those of Rhodopirellula baltica through 16 s rDNA sequencing analysis, were observed on the antibiotic susceptibility testing plate of KMHK cells containing only penicillin, streptomycin and gentamicin in our study. Tetracycline was therefore included in the present antibiotic treatment based on previous report that Rhodopirellula sp. was highly susceptible to 0.5 ppm of tetracycline33. It is common to modify the antibiotic cocktail for different algal cultures since the antibiotics depends on the microbiome. For instance, Su et al. treated Alexandrium cultures with a combination of gentamycin, streptomycin, cephalothin and rifampicin for 7 days to obtain axenic cultures13.
The effectiveness of the basic and extended protocols to remove bacteria in the algal samples is summarized in Fig. 3. The total bacterial count significantly reduced from initial 5.79 ± 0.22 log10 CFU/mL to 4.88 ± 0.05 log10 CFU/mL (p ≤ 0.05) after centrifuged with 90% Percoll (Step 1), even though the primary aim of this step was to condense the algal cells. This is probably due to the removal of numerous free-living bacteria in the supernatant, which was discarded. In the subsequent two gradient centrifugation using 90–50–30% density layers (Steps 2 and 3), approximately 12% of bacteria were further removed at p ≤ 0.05. The bacterial count did not show any additional reduction even when a step of a 90–50–30% density gradient centrifugation was added between Steps 3 and 4 (data not shown). After Step 4 with 48-h antibiotic treatment, the bacterial count was significantly reduced by 31% (Fig. 3a). These results indicated that numerous algae-associated bacteria could be effectively inhibited using the antibiotics but these steps could not completely eradicate the bacteria. No further decline in the bacterial count was found after Step 5 but decreased significant after Step 6, although both steps employed the same gradient centrifugation (90–50–30%). This could probably be the killing effect of the extended incubation of the bacteria with the remaining antibiotics, and/or the bacterial count significantly reduced by the gradient centrifugation in Step 5 was offset by the intracellular bacteria released from algal cell lysis during the antibiotic treatment. The effect of antibiotic exposure time on bacterial removal efficiency was shown in Table 1 below, while the existence of intracellular bacteria inside the dinoflagellate cells remains controversial6 and deserves more in-depth studies. The residue antibiotics, bacteria and algal cell debris were reported to suppress the algal cell growth17,31, these suppressions could be effectively removed in the present study as shown by algal regrowth (Fig. 3).
Figure 3

Bacterial removal in the KMHK sample using the basic and extended protocols. (a) Total bacterial count after different steps in the basic protocol. Initial: the initial bacterial count in KMHK samples at the beginning; Supplementary Fig. 1 illustrate the Step 1 to 6 of basic protocol. (b) Total bacterial count after the basic and extended protocols. Extended protocol refers to the descriptions in material and method. (c) Total bacterial count in the treated KMHK cultures at different days of cultivation. (d) Algal cell concentration during the regrowth of the treated KMHK cultures. All data are presented as mean ± standard deviations of three independent experiments (n = 3). Different letters on the top of the bar indicate that means were significantly different among samples at p ≤ 0.05 according to one-way analysis of variance followed by Tukey multiple comparison tests.

Full size image

After all six steps of the present protocol, the total bacterial count remaining in the KMHK culture was 1.13 ± 0.07 log10 CFU/mL (equivalent to approximately 13 CFU/mL), indicating that  > 99.9% of the bacteria were removed. Nevertheless, the bacterial count increased significantly to 5.75 ± 0.14 log10 CFU/mL (equivalent to approximately 5.86 × 105 CFU/mL) in day 7 KMHK regrow culture after the protocol (Fig. 3c). We used the extended protocol to remove the remaining bacteria by repeating the 48-h antibiotic incubation and Steps 4–6 after the basic protocol (Fig. 3b). However, the total bacterial count did not have any significant change between the basic and extended protocols (Fig. 3b), demonstrating that the additional steps, that is, repeating the 48-h antibiotic incubation and Steps 4–6, in the extended protocol failed to eliminate the few remaining bacteria in the KMHK culture. Similar to the basic protocol, the remaining bacteria regrew significantly in day 3 and day 7 KMHK cultures after the extended one (Fig. 3c). More, the additional steps of antibiotic treatment and centrifugation even damaged the algal cells, as reflected by the regrowth of KMHK cells was strongly inhibited by the extended protocol (Fig. 3d). After the extended protocol, the KMHK cell density at day 7 was only half of that after the basic protocol that regrew normally and reached the cell density of approximately 10,000 cells/mL after 7 days, comparable to that of a routine culture. The algal growth rates (µ) after the basic and extended protocols were 1.29 and 1.14, respectively. This implied that the additional steps of antibiotic treatment in the extended protocol damaged the algal cells. It has been reported that antibiotics treatment could interfere the peptidoglycan biosynthesis and eventually inhibited chloroplast division34. The extended protocol was also ineffective in removing bacteria, probably due to limited bacterial removal capacity and/or insufficient dose and exposure time of antibiotics.
Effect of initial algal cell density and antibiotic exposure times on bacterial removal
Bacterial cell concentration, antibiotic dose and antibiotic exposure time are critical factors affecting the efficiency of bacterial removal in algal samples. We hypothesized that the incomplete bacterial removal after the basic protocol was attributable to (1) the bacterial concentration in the algal culture exceeded the treatment capacity and (2) the dose and exposure time used in the antibiotic treatment were insufficient. To test insufficient dosing, a double antibiotic dose was used but there was no significant difference in the amount of bacterial removal between the normal and double doses of antibiotics used in the treatments (data not shown). Our observations also showed that the cell density of the 7-day KMHK culture treated with a double dose of antibiotics decreased from 10,000 to 489 cells/mL, clearly indicating a high antibiotic dose severely damaged the algal cells.
Both the amount and percentage of bacterial removal were independent of the initial algal cell density (Table 1). At least 94.49% of bacteria were removed in all treatments. With a 48-h antibiotic exposure time, the percentage of bacterial removal did not have any significant changes with increases of initial algal cell density, and were 94.49%, 99.84% and 99.93% at high, moderate and low densities, respectively. Similarly, no significant difference in the amount of bacterial removal (in terms of Log10 CFU/mL) was observed between low and moderate initial algal cell densities but were significantly higher than that at high initial algal cell density. With 96-h of antibiotic exposure, bacterial removal percentage was 100% at both high and low algal densities and 98.53% at moderate algal density. Similarly, bacterial removal of low and moderate initial algal cell densities were significantly higher than that at high initial algal cell density. These results indicated that the bacterial removal ability was independent of the initial algal cell density.
Table 1 Bacterial removal in the algal cultures with different initial algal cell densities and antibiotic exposure times in the basic protocol.
Full size table

When compare the antibiotic exposure time, bacterial removal with 48-h antibiotic treatment was significantly higher than that with other exposure times, regardless of the initial algal density (Table 1). A complete bacterial removal (100%) was detected under three conditions, that is, high and low initial densities with 96-h exposure time and moderate initial algal density with 72-h exposure time. The antibiotic exposure times ranging from several hours to 1 week have been reported in pervious researches, depending on the method used, and the bacterial cell count and composition in the algal sample13,32. The exposure time should be considered and controlled carefully because prolonged exposure to antibiotics may damage the algal cells and suppress their growth. However, the present study reveals antibiotic exposure time was not a critical factor for achieving 100% bacterial removal. On the other hand, the bacterial removal might be related to the bacterial cell count and composition present in the algal culture at the beginning. The microbiome can also change frequently during the routine cultivation of the algal cultures35.
Algal cells obtained from different protocols were then cultured for 7 days and then subjected to bacterial counting. A substantial number of bacteria was found in all treatments at the end of the algal regrow cultures, including the 72-h exposure achieving 100% bacterial removal, except two 96-h antibiotic exposure treatments (Table 1). Even in these two 96-h exposure treatments with 100% bacteria-free algal cultures, bacteria were detected after several generations of the algal culture (data not shown). These indicated that the axenicity of algal culture may not be guaranteed even when no bacteria are detected after treating the algal samples with our developed protocol.
The reappearance of bacterial growth may be due to the existence of a very small proportion of bacteria attached onto the algal cells after the protocol which was too little to be detected. Dinoflagellates are covered with complex and irregular cell surface structures36. Steric hindrance from parts of these algal surface areas may protect the firmly attached bacteria, making their detachment from the algal cells during gradient centrifugation difficult and decreasing antibiotic accessibility. Another possible explanation is that a few bacteria may have developed antibiotic resistance11, but this is unlikely in the present study as antibiotic resistance usually develops when the bacteria are continuously exposed to a nonlethal dose of an antibiotic37. Although why a few bacteria remained on the algal cells and regrew rapidly along with algal cell growth are poorly understood, it is necessary to have additional treatments such as serial dilution selection to ensure a true axenic algal culture is obtained.
Selection of axenic algal cells through serial dilution
After 7 days of cultivation, KMHK cells from all dilutions described in Fig. 1 survived. No bacterial growth was observed in 42 KMHK cultures (out of total 45 cultures), although bacterial colonies were observed in one of the day 7 KMHK cultures (one culture in the 100 dilution) and in another two of the day 21 KMHK cultures (one in 100 dilution and one in 10–1 dilution) in the first trial (Supplementary Table 1). The results reiterate that a high proportion of the KMHK cells in the population is in fact axenic and the ratio of axenic to non-axenic clones in the algal cell population is high after the basic protocol, it is therefore highly feasible to obtain the axenic clones from the population through such serial dilution approach. Similar approaches have been reported in previous studies5,12,38. For instance, Ki and Han dispersed the algal cells in a 96-well plate after filtration and antibiotic treatment12. Sena et al. also serially diluted the cyanobacterial sample, Arthrospira spp., after antibiotics treatment38. Although the axenic clone could be obtained by plating the algal cells on agar5, this approach was not feasible for marine dinoflagellates because the cells were unable to grow on a solid medium.
Verification of axenicity of algal cultures
It has been reported that some marine bacteria grow very slowly on agar, something like 50 days, and some of them are unculturable39,40. Therefore, confirming the axenic state of the algal cultures is paramount. The axenic state of the two selected KMHK cultures was tested and results of DAPI epifluorescence microscopy show that no bacteria were observed in the treated KMHK cells (Fig. 4b) while bacteria were found in the untreated cells (Fig. 4a). For the rDNA sequencing analysis, a 1500-bp PCR product was obtained after bacterial 16S rDNA amplification (Fig. 5a). The BLAST of the sequence reveal that it shared 99.59% similarity with 16 s rDNA sequence in the plastid gene of K. mikimotoi (accession no. AB027236). Similarly, the 600-bp amplicon observed in the amplification of fungal ITS shared 99.67% similarity with the ITS sequence of K. mikimotoi (accession no. KT733616; Fig. 5c). These results confirm the absence of both culturable and unculturable bacteria and fungi in the treated KMHK cultures. It has been reported that the algal cultures must continually be treated with antibiotics in order to maintain their axenic status15, but algal cells may die after several sub-cultures because of prolonged antibiotic exposure. When this happens, it is usually too late to recover the cultures. In the present study, regular monitoring of the axenicity of the cultures was performed through bacterial colony counting, DAPI epifluorescence microscopy and rDNA sequence analysis. The established axenic cultures were maintained generations after generations without adding any antibiotics, and no bacteria were found in any of the sub-cultures being tested even after 30 generations (data not shown). The axenic cultures of KMHK were established successfully and maintained sustainably, indicating this methodology was a promising approach applicable to other unarmored dinoflagellates. To the best of our knowledge, this is the first successful establishment of an axenic culture for the unarmored dinoflagellate K. mikimotoi.
Figure 4

DAPI epifluorescence microscopic images of KMHK and AT6 samples under ×1000 magnification: (a) untreated (control) and (b) treated KMHK samples; (c) untreated (control) and (d) treated AT6 samples.

Full size image

Figure 5

PCR amplification of bacterial 16 s rDNA for (a) KMHK and (b) AT6 samples and that of fungal ITS region for (c) KMHK and (d) AT6 samples obtained from basic protocol and serial dilution. +ve: positive control, -ve: negative control.

Full size image

Development of axenic culture of A. tamarense using our established methodology
The established method was applied to obtain the axenic cultures of another dinoflagellate species, A. tamarense (AT6), a well-known paralytic shellfish toxin–producing agent which has been extensively studied in the past decades41,42,43,44. The AT6 culture with an initial bacterial count of 7.9 ± 0.08 log10 CFU/mL was subjected to the basic protocol and the total bacterial counts recorded after Steps 3, 4 and 6 were shown in Fig. 6a. The result was generally similar to that of KMHK culture using the basic protocol (Fig. 3a), except no bacteria were detected in the AT6 culture after Step 6. Even though 100% bacterial removal was achieved, bacterial regrowth was observed on days 3 and 7 of the treated AT6 culture (Fig. 6b). The bacterial count regrew significantly to 6.71 ± 0.08 log10 CFU/mL after 7 days of cultivation (Fig. 6c). The regrowth of bacteria from the treated AT6 culture achieving 100% bacterial removal confirmed that few bacteria attached at some points onto the algal surface were shielded. Biegala et al. reported that associated bacteria were attached onto the cell surface within the sulci and cingula of A. tamarense6.
Figure 6

Bacterial removal in the Alexandrium tamarense (AT6) samples using the basic protocol. (a) Total bacterial count against different steps. Initial: the initial bacterial count present in AT6 at the beginning. (b) Total bacterial count in the AT6 culture obtained after the protocol at different days of cultivation. (c) Algal cell concentration during the regrowth of the AT6 culture obtained after the protocol. All data are presented as means ± standard deviations of three independent experiments (n = 3). N.D.: not detected. Different letters on the top of the bar indicate that means were significantly different among different samples at p ≤ 0.05 according to one-way analysis of variance followed by Tukey multiple comparison tests.

Full size image

For the serial dilution selection of AT6 cells after the basic protocol, culturable bacteria were observed in 11 of the 15 cultures in 100 dilution (5/5 in trial 1; 4/5 in trial 2 and 2/5 in trial 3) after 7 days of algal cultivation. All these 15 cultures showed bacterial regrowth after 21 days of algal cultivation, but the number of cultures with bacterial regrowth decreased to 1 and 4 of the 15 cultures in 10–1 dilutions after 7 and 21 days of algal cultivation, respectively. In 10−2 dilutions, no culturable bacteria were observed in all the cultures after both 7 and 21 days of algal cultivation (Supplementary Table 1). These results further demonstrate the feasibility of using the stepwise serial dilution method to select axenic algae, and 10–2 dilutions offer the highest probability in acquiring the axenic clones. The bacterial status of two of these potential axenic AT6 cultures was further assessed through DAPI epifluorescence microscopy and bacterial rDNA and fungal ITS sequencing analysis. No bacteria were observed in the DAPI epifluorescence image of the treated AT6 cultures compared to the untreated control cultures (Figs. 4c,d). Neither bacterial 16 s rDNA band (Fig. 5b) nor fungal ITS region (Fig. 5d) was amplified in the treated AT6 samples. These results confirmed the axenic status of the AT6 cultures.
Our established methodology
The present results demonstrate the potential of our methodology to be used in the establishment of axenic cultures for both armored and unarmored dinoflagellates. Figure 7 summarizes the workflow and procedures of our methodology. This promising approach combines three techniques, Percoll density gradient centrifugation, antibiotic treatment and serial dilution. Density gradient centrifugation considerably reduces the bacterial population by the physical separation between the associated bacteria, mainly the free-living and loosely attached bacteria, and the dinoflagellate cells on the basis of cell size. Percoll density layers not only provide a matrix for separating the two types of cells effectively but also cushion the dinoflagellate cells against the impact of the mechanical force. The Percoll density layers together with the bactericidal action of the antibiotic treatment typically eradicate  > 99% of the associated bacteria from the dinoflagellate culture. Our strategies not target at removing the remaining  More

1.
Mendil, D. et al. Determination of trace metal levels in seven fish species in lakes in Tokat, Turkey. Food Chem. 90, 175–179 (2005).
CAS  Article  Google Scholar 
2.
Kennish, M. J. Pollution in estuaries and coastal marine waters. J. Coast. Res. 12, 27–49 (1994).
Google Scholar 

3.
de Souza Machado, A. A., Spencer, K., Kloas, W., Toffolon, M. & Zarfl, C. Metal fate and effects in estuaries: A review and conceptual model for better understanding of toxicity. Sci. Total Environ. 541, 268–281 (2016).
ADS  PubMed  Article  CAS  PubMed Central  Google Scholar 

4.
Mamtani, R., Stern, P., Dawood, I. & Cheema, S. Metals and disease: A global primary health care perspective. J. Toxicol. 2011, 1–11 (2011).
Article  CAS  Google Scholar 

5.
van Dam, J. W., Negri, A. P., Uthicke, S. & Mueller, J. F. Chemical pollution on coral reefs: Exposure and ecological effects. In Ecological Impacts of Toxic Chemicals (eds Sanchez-Bayo, F. et al.) 187–211 (Bentham Science Publishers, Oak Park, 2011).
Google Scholar 

6.
Dong, W. et al. Developmental toxicity from exposure to various forms of mercury compounds in medaka fish (Oryzias latipes) embryos. PeerJ 4, e2282 (2016).
PubMed  PubMed Central  Article  CAS  Google Scholar 

7.
Morcillo, P., Esteban, M. A. & Cuesta, A. Mercury and its toxic effects on fish. AIMS Environ. Sci. 4, 386 (2017).
CAS  Article  Google Scholar 

8.
Bosch, A. C., O’Neill, B., Sigge, G. O., Kerwath, S. E. & Hoffman, L. C. Heavy metal accumulation and toxicity in smoothhound (Mustelus mustelus) shark from Langebaan Lagoon, South Africa. Food Chem. 190, 871–878 (2016).
CAS  PubMed  Article  PubMed Central  Google Scholar 

9.
Fisher, N. S. & Reinfelder, J. R. The trophic transfer of metals in marine systems. In Metal Speciation and Bioavailability in Aquatic Systems, Vol. 3 (eds. Tessier, A. & Turner, D. R.) 407–411, 363–406 (Wiley, Chichester, 1995).

10.
Ali, H. & Khan, E. Bioaccumulation of non-essential hazardous heavy metals and metalloids in freshwater fish. Risk to human health. Environ. Chem. Lett. 16, 903–917 (2018).
CAS  Article  Google Scholar 

11.
Domi, N., Bouquegneau, J. M. & Das, K. Feeding ecology of five commercial shark species of the Celtic Sea through stable isotope and trace metal analysis. Mar. Environ. Res. 60, 551–569 (2005).
CAS  PubMed  Article  PubMed Central  Google Scholar 

12.
Rumbold, D., Wasno, B., Hammerschlag, N. & Volety, A. Mercury accumulation in sharks from the coastal waters of Southwest Florida. Arch. Environ. Contam. Toxicol. 67, 402–412 (2014).
CAS  PubMed  Article  PubMed Central  Google Scholar 

13.
Lee, C. S. et al. Declining mercury concentrations in bluefin tuna reflect reduced emissions to the North Atlantic Ocean. Environ. Sci. Technol. 50, 12825–12830 (2016).
ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

14.
Merly, L. et al. Blood plasma levels of heavy metals and trace elements in white sharks (Carcharodon carcharias) and potential health consequences. Mar. Pollut. Bull. 142, 85–92 (2019).
CAS  PubMed  Article  PubMed Central  Google Scholar 

15.
Fisk, A. T., Tittlemier, S. A., Pranschke, J. L. & Norstrom, R. J. Using anthropogenic contaminants and stable isotopes to assess the feeding ecology of Greenland sharks. Ecology 83, 2162–2172 (2002).
Article  Google Scholar 

16.
Buck, D. G. et al. A global-scale assessment of fish mercury concentrations and the identification of biological hotspots. Sci. Total Environ. 687, 956–966 (2019).
ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 

17.
Tiktak, G. P. et al. Are concentrations of pollutants in sharks, rays and skates (Elasmobranchii) a cause for concern? A systematic review. Mar. Pollut. Bull. 160, 111701 (2020).
CAS  PubMed  Article  PubMed Central  Google Scholar 

18.
Boening, D. W. Ecological effects, transport, and fate of mercury: A general review. Chemosphere 40, 1335–1351 (2000).
ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 

19.
Sevcikova, M., Modra, H., Slaninova, A. & Svobodova, Z. Metals as a cause of oxidative stress in fish: A review. Vet. Med. 56, 537–546 (2011).
CAS  Article  Google Scholar 

20.
Bezerra, M. F., Lacerda, L. D. & Lai, C. T. Trace metals and persistent organic pollutants contamination in batoids (Chondrichthyes: Batoidea): A systematic review. Environ. Pollut. 248, 684–695 (2019).
CAS  PubMed  Article  PubMed Central  Google Scholar 

21.
Cortés, E. Standardized diet compositions and trophic levels of sharks. ICES J. Mar. Sci. 56, 707–717 (1999).
Article  Google Scholar 

22.
Hussey, N. E. et al. Rescaling the trophic structure of marine food webs. Ecol. Lett. 17, 239–250 (2014).
PubMed  Article  PubMed Central  Google Scholar 

23.
Shipley, O. N., Gallagher, A. J., Shiffman, D. S., Kaufman, L. & Hammerschlag, N. Diverse resource-use strategies in a large-bodied marine predator guild: Evidence from differential use of resource subsidies and intraspecific isotopic variation. Mar. Ecol. Prog. Ser. 623, 71–83 (2019).
ADS  Article  Google Scholar 

24.
Roff, G. et al. The ecological role of sharks on coral reefs. Trends Ecol. Evol. 31, 395–407 (2016).
PubMed  Article  PubMed Central  Google Scholar 

25.
Hammerschlag, N. et al. Ecosystem function and services of aquatic predators in the Anthropocene. Trends Ecol. Evol. 34, 369–383 (2019).
PubMed  Article  PubMed Central  Google Scholar 

26.
Dent, F. & Clarke, S. State of the global market for shark products. FAO Fisheries and Aquaculture Technical Paper, 590 (2015).

27.
Buchan, K. C. The Bahamas. Mar. Pollut. Bull. 41, 94–111 (2000).
CAS  Article  Google Scholar 

28.
Brooks, E. J., Sloman, K. A., Sims, D. W. & Danylchuk, A. J. Validating the use of baited remote underwater video surveys for assessing the diversity, distribution and abundance of sharks in the Bahamas. Endanger. Spec. Res. 13, 231–243 (2011).
Article  Google Scholar 

29.
Zhu, Y., Newman, S. P., Reid, W. D. & Polunin, N. V. Fish stable isotope community structure of a Bahamian coral reef. Mar. Biol. 166, 160 (2019).
Article  Google Scholar 

30.
Burgess, G. H. & Fordham, S. V. Regional overviews. In Sharks, Rays & Chimaeras: The Status of Chondrichthyan Fishes (eds. Fowler, S. L. et al.) 461 (IUCN/SSG Shark Specialist Group, Gland, 2005).

31.
Sherman, K. D. et al. Contemporary and emerging fisheries in The Bahamas—Conservation and management challenges, achievements and future directions. Fish. Manag. Ecol. 25, 319–331 (2018).
Article  Google Scholar 

32.
Ward-Paige, C. A. A global overview of shark sanctuary regulations and their impact on shark fisheries. Mar. Policy 82, 87–97 (2017).
Article  Google Scholar 

33.
Hammerschlag, N., Gallagher, A. J., Wester, J., Luo, J. & Ault, J. S. Don’t bite the hand that feeds: Assessing ecological impacts of provisioning ecotourism on an apex marine predator. Funct. Ecol. 26, 567–576 (2012).
Article  Google Scholar 

34.
Graham, F. et al. Use of marine protected areas and exclusive economic zones in the subtropical western North Atlantic Ocean by large highly mobile sharks. Divers. Distrib. 22, 534–546 (2016).
Article  Google Scholar 

35.
Chapman, D. D., Pikitch, E. K., Babcock, E. A. & Shivji, M. S. Deep-diving and diel changes in vertical habitat use by Caribbean reef sharks Carcharhinus perezi. Mar. Ecol. Prog. Ser. 344, 271–275 (2007).
ADS  Article  Google Scholar 

36.
Shipley, O. N. et al. Horizontal and vertical movements of Caribbean reef sharks (Carcharhinus perezi): Conservation implications of limited migration in a marine sanctuary. R. Soc. Open Sci. 4, 160611 (2017).
ADS  PubMed  PubMed Central  Article  Google Scholar 

37.
Shipley, O. N. et al. Fine-scale movement and activity patterns of Caribbean reef sharks (Carcharhinus perezi) in the Bahamas. Environ. Biol. Fish. 101, 1097–1104 (2018).
Article  Google Scholar 

38.
Rosa, R.S., Mancini, P., Caldas, J.P., Graham, R.T. Carcharhinus perezi. The IUCN Red List of Threatened Species 2006: e.T60217A12323052. https://doi.org/10.2305/IUCN.UK.2006.RLTS.T60217A12323052.en (2006). Accessed 23 Aug 19.

39.
Wei, T. & Simko, V. R package “corrplot”: Visualization of a Correlation Matrix (Version 0.84) (2017).

40.
Wood, S. N. Generalized Additive Models: An Introduction with R 2nd edn. (Chapman and Hall/CRC Press, Boca Raton, 2017).
Google Scholar 

41.
Adams, D. H. & McMichael, R. H. Jr. Mercury levels in four species of sharks from the Atlantic coast of Florida. Fish. Bull. 97, 372–379 (1999).
Google Scholar 

42.
Matulik, A. G. et al. Bioaccumulation and biomagnification of mercury and methylmercury in four sympatric coastal sharks in a protected subtropical lagoon. Mar. Polllut. Bull. 116, 357–364 (2017).
CAS  Article  Google Scholar 

43.
Gelsleichter, J., Sparkman, G., Howey, L. A., Brooks, E. J. & Shipley, O. N. Elevated accumulation of the toxic metal mercury in the Critically Endangered oceanic whitetip shark Carcharhinus longimanus from the northwestern Atlantic Ocean. Endanger. Spec. Res. 43, 267–279 (2020).
Article  Google Scholar 

44.
Hammerschlag, N., Skubel, R. A., Sulikowski, J., Irschick, D. J. & Gallagher, A. J. A comparison of reproductive and energetic states in a marine apex predator (the tiger shark, Galeocerdo cuvier). Physiol. Biochem. Zool. 91, 933–942 (2018).
PubMed  Article  PubMed Central  Google Scholar 

45.
Froese, R., & Pauly, D. FishBase, version (06/2017). World Wide Web Electronic Publication (2017). Accessed 23 Aug 19.

46.
Pikitch, E. K., Chapman, D. D., Babcock, E. A. & Shivji, M. S. Habitat use and demographic population structure of elasmobranchs at a Caribbean atoll (Glover’s Reef, Belize). Mar. Ecol. Prog. Ser. 302, 187–197 (2005).
ADS  Article  Google Scholar 

47.
Dulvy, N. K. et al. Extinction risk and conservation of the world’s sharks and rays. Elife 3, e00590 (2014).
PubMed  PubMed Central  Article  Google Scholar 

48.
Castro-González, M. I. & Méndez-Armenta, M. Heavy metals: Implications associated to fish consumption. Environ. Toxicol. Pharmacol. 26, 263–271 (2008).
PubMed  Article  CAS  PubMed Central  Google Scholar 

49.
Amirah, M. N., Afiza, A. S., Faizal, W. I. W., Nurliyana, M. H. & Laili, S. Human health risk assessment of metal contamination through consumption of fish. J. Environ. Pollut. Hum. Health 1, 1–5 (2013).
Google Scholar 

50.
Lyons, K. et al. Species-specific characteristics influence contaminant accumulation trajectories and signatures across ontogeny in three pelagic shark species. Environ. Sci. Technol. 53, 6997–7006 (2019).
ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 

51.
Rejomon, G., Nair, M. & Joseph, T. Trace metal dynamics in fishes from the southwest coast of India. Environ. Monit. Assess. 167, 243–255 (2010).
CAS  PubMed  Article  PubMed Central  Google Scholar 

52.
Shipley, O. N. et al. Trophodynamics and mercury bioaccumulation in reef and open-ocean fishes from The Bahamas with a focus on two teleost predators. Mar. Ecol. Prog Ser. 608, 221–232 (2019).
ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

53.
Karimi, R., Fisher, N. S. & Folt, C. L. Multielement stoichiometry in aquatic invertebrates: When growth dilution matters. Am. Nat. 176, 699–709 (2010).
PubMed  Article  PubMed Central  Google Scholar 

54.
Choy, C. A., Popp, B. N., Kaneko, J. J. & Drazen, J. C. The influence of depth on mercury levels in pelagic fishes and their prey. Proc. Natl. Acad. Sci. U. S. A. 106, 13865–13869 (2009).
ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

55.
Lacerda, L. D., Goyanna, F., Bezerra, M. F. & Silva, G. B. Mercury concentrations in tuna (Thunnus albacares and Thunnus obesus) from the Brazilian Equatorial Atlantic Ocean. B Eniviron. Contam. Toxicol. 98, 49–155 (2017).
Google Scholar 

56.
Lee, C. S. & Fisher, N. S. Methylmercury uptake by diverse marine phytoplankton. Limnol. Oceanogr. 61, 1626–1639 (2016).
ADS  PubMed  PubMed Central  Article  Google Scholar 

57.
Mohammed, A. & Mohammed, T. Mercury, arsenic, cadmium and lead in two commercial shark species (Sphyrna lewini and Caraharinus porosus) in Trinidad and Tobago. Mar. Pollut. Bull. 119, 214–218 (2017).
CAS  PubMed  Article  PubMed Central  Google Scholar 

58.
Burger, J. et al. Metal levels in fish from the Savannah River: potential hazards to fish and other receptors. Environ Res 89, 85–97 (2002).
ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 

59.
Kim, S. W. et al. Heavy metal accumulation in and food safety of shark meat from Jeju island, Republic of Korea. PLoS ONE 14, e0212410 (2019).
CAS  PubMed  PubMed Central  Article  Google Scholar 

60.
Roesijadi, G. Metallothioneins in metal regulation and toxicity in aquatic animals. Aquat. Toxicol. 22, 81–113 (1992).
CAS  Article  Google Scholar 

61.
Palmiter, R. D. The elusive function of metallothioneins. Proc. Natl. Acad. Sci. U. S. A. 95, 8428–8430 (1998).
ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

62.
Pagenkopf, G. K. Gill surface interaction model for trace-metal toxicity to fishes: Role of complexation, pH, and water hardness. Environ. Sci. Technol. 17, 342–347 (1983).
ADS  CAS  Article  Google Scholar 

63.
Playle, R. C. Modelling metal interactions at fish gills. Sci. Total Environ. 219, 147–163 (1998).
ADS  CAS  Article  Google Scholar 

64.
Barrera-García, A. et al. Trace elements and oxidative stress indicators in the liver and kidney of the blue shark (Prionace glauca). Comp. Biochem. Physiol. A 165, 483–490 (2013).
Article  CAS  Google Scholar 

65.
Dhanakumar, S., Solaraj, G. & Mohanraj, R. Heavy metal partitioning in sediments and bioaccumulation in commercial fish species of three major reservoirs of river Cauvery delta region, India. Ecotox. Environ. Safe 113, 145–151 (2015).
CAS  Article  Google Scholar 

66.
Buchanan, K. et al. Guidelines for the treatment of animals in behavioural research and teaching. Anim. Behav. 83, 301–309 (2012).
Article  Google Scholar 

67.
Bergés-Tiznado, M. E. et al. Mercury and selenium in muscle and target organs of Scalloped Hammerhead Sharks Sphyrna lewini of the SE Gulf of California: Dietary intake, molar ratios, loads, and human health risks. Arch. Environ. Contam. Toxicol. 69, 440–452 (2015).
PubMed  Article  CAS  PubMed Central  Google Scholar  More

1.
Batalha, L. S., Silva Filho, D. F. & Martius, C. Using termite nests as a source of organic matter in agrosilvicultural production systems in Amazonia. Scientia Agricola 52, 318–325 (1995).
Article  Google Scholar 
2.
Ayuke, F. O. et al. Soil fertility management: Impacts on soil macrofauna, soil aggregation and soil organic matter allocation. Appl. Soil. Ecol. 48, 53–62 (2011).
Article  Google Scholar 

3.
Jouquet, P., Chaudhary, E. & Kumar, A. R. V. Sustainable use of termite activity in agro-ecosystems with reference to earthworms A review. Agron. Sustain. Dev. 38, 3 (2018).
Article  Google Scholar 

4.
Deligne, J., Quennedey, A. & Blum, M. S. The enemies and defense mechanisms of termites. In Social Insects (ed. Hermann, H. R.) 1–76 (Academic Press, Cambridge, 1981).
Google Scholar 

5.
Grasse, P. P. Termitologia, Tome III (Masson, Paris, 1986).
Google Scholar 

6.
Prestwich, G. D. Defense mechanisms of termites. Annu. Rev. Entomol. 29(1), 201–232 (1984).
CAS  Article  Google Scholar 

7.
Chuah, C. H. Chemical Weapons and Defense Mechanism of Malaysian Termites. In Chemistry in Malaysia 4–11 (Institut Kimia Malaysia, 2010).

8.
Iida, M. & Akino, T. Defensive effect of soldier-specific secretion by Reticulitermes speratus (Isoptera: Rhinotermitidae) on the facultative termitophagous ponerine ant, Brachyponera chinensis (Hymenoptera: Ponerinae). Appl. Entomol. Zool. 51, 111–116 (2016).
Article  Google Scholar 

9.
Kori, N. S. M. & Arumugam, N. Termites of Agropark, Universiti Malaysia Kelantan, Jeli Campus: Diversity and pest composition. J. Trop. Resour. Sustain. Sci. 5, 104–108 (2017).
Google Scholar 

10.
Alia-Diyana, M. H., Appalasamy, S. & Arumugam, N. Termite species and structural pest identification in selected rural areas of Kelantan, Malaysia. IOP Conf. Ser. Earth and Environ. Sci. https://doi.org/10.1088/1755-1315/549/1/012053 (2020).
Article  Google Scholar 

11.
Sillam-Dussès, D. et al. Comparative Study of the Labial Gland Secretion in Termites (Isoptera). PLoS ONE 7(10), e46431 (2012).
ADS  PubMed  PubMed Central  Article  CAS  Google Scholar 

12.
Wyatt, T. D. Pheromones and other chemical communication in animals. In Encyclopedia of Neuroscience (ed. Squire, L. R.) 611–616 (Academic Press, Oxford, 2017).
Google Scholar 

13.
Ahmad, N. & Kamarudin, N. Pheromone Trapping in Controlling Key Insect Pests: Progress and Prospects (Malaysia Palm Oil Board, Kajang, 2016).
Google Scholar 

14.
Matthews, G. Pesticides: Health, Safety and the Environment (Wiley Blackwell, Chichester, 2015).
Google Scholar 

15.
Desneux, N., Decourtye, A. & Delpuech, J.-M. The sublethal effects of pesticides on beneficial arthropods. Ann. Rev. Entomol. 52(1), 81–106 (2007).
CAS  Article  Google Scholar 

16.
Rattan, R. Mechanism of action of insecticidal secondary metabolites of plant origin. Crop Prot. 29(9), 913–920 (2010).
CAS  Article  Google Scholar 

17.
Regnault-Roger, C., Vincent, C. & Arnason, J. T. Essential oils in insect control: Low-risk products in a high-stakes world. Annu. Rev. Entomol. 57, 405–424 (2012).
CAS  PubMed  Article  Google Scholar 

18.
Pavela, R. & Benelli, G. Essential oils as ecofriendly biopesticides? Challenges and constraints. Trends Plant Sci. 21(12), 1000–1007 (2016).
CAS  PubMed  Article  Google Scholar 

19.
Chouvenc, T., Su, N. & Kenneth, G. J. Fifty years of attempted biological control of termites—Analysis of a failure. Biol. Control 59(2), 69–82 (2011).
Article  Google Scholar 

20.
Meikle, W. G. et al. Evaluation of an entomopathogenic fungus, Paecilomyces fumosoroseus (Wize) Brown and Smith (Deuteromycota: Hyphomycetes) obtained from Formosan subterranean termites (Isop, Rhinotermitidae). J. Appl. Entomol. 129(6), 315–322 (2005).
Article  Google Scholar 

21.
Tho, Y. P. Termites of Peninsular Malaysia (Forest Research Institute Malaysia, Selangor, 1992).
Google Scholar 

22.
Krasulova, J. et al. Chemistry and anatomy of the frontal gland in soldiers of the sand termite Psammotermes hybostoma. J. Chem. Ecol. 38(5), 557–565 (2012).
CAS  PubMed  Article  Google Scholar 

23.
Bakaruddin, N. H., Dieng, H., Sulaiman, S. F. & Ab Majid, A. H. Evaluation of the toxicity and repellency of tropical plant extract against subterranean termites, Globitermes sulphureus and Coptotermes gestroi. Inf. Process. Agric. 5(3), 298–307 (2018).
Google Scholar 

24.
Lee, C. C. & Lee, C. Y. A laboratory maintenance regime for a fungus-growing termite Macrotermes gilvus (Blattodea: Termitidae). J. Econ. Entomol. 108(3), 1243–1250 (2015).
CAS  PubMed  Article  Google Scholar 

25.
Zibaee, I. & Pooya, B. K. Evaluation of repellent activity of two essential oils and their mixed formulation against cockroaches (Dictyoptera: Blattidae, Blattellidae) in Iran. J. Entomol. Zool. Stud. 4, 106 (2016).
Google Scholar 

26.
OECD. Guidance Document on Assays for Testing the Efficacy of Baits Against Cockroaches Health and Safety Publications (OECD Environment, Paris, 2013).
Google Scholar 

27.
Syed, R., Manzoor, F., Adalat, R., Abdul-Sattar, A. & Syed, A. Laboratory evaluation of toxicity of insecticide formulations from different classes against American cockroach (Dictyoptera: Blattidae). J. Arthropod-Borne Dis. 8(1), 21–34 (2014).
PubMed  Google Scholar 

28.
Ohta, M., Matsuura, F., Henderson, G. & Laine, R. A. Novel free ceramides as components of the soldier defense gland of the Formosan subterranean termite (Coptotermes formosanus). J. Lipid Res. 48(3), 656–664 (2007).
CAS  PubMed  Article  Google Scholar 

29.
McDonald, L. L., Guy, R. H. & Speirs, R. D. Preliminary evaluation of new candidate materials as toxicants, repellents, and attractants against stored-product insects-1 (Agriculture Research Service, 1970).

30.
Johnson, R. A., Thomas, R. J., Wood, T. G. & Swift, M. J. The inoculation of the fungus comb in newly founded colonies of some species of the Macrotermitinae (Isoptera) from Nigeria. J. Nat. Hist. 15(5), 751–756 (2007).
Article  Google Scholar 

31.
de Mello, A. P., Azevedo, N. R., da Silva, A. M. B. & Gusmão, M. A. B. Chemical composition and variability of the defensive secretion in Nasutitermes corniger (Motschulsky, 1885) in urban area in the Brazilian semiarid region. Entomotropica 31, 82–90 (2016).
Google Scholar 

32.
Bordereau, C., Robert, A., Van Tuyen, V. & Peppuy, A. Suicidal defensive behaviour by frontal gland dehiscence in Globitermes sulphureus Haviland soldiers (Isoptera). Insectes Soc. 44(3), 289–297 (1997).
Article  Google Scholar 

33.
Touchard, A., Dejean, A., Escoubas, P. & Orivel, J. Intraspecific variations in the venom peptidome of the ant Odontomachus haematodus (Formicidae: Ponerinae) from French Guiana. J. Hymenoptera Res. 47, 87–101 (2015).
Article  Google Scholar 

34.
Aguilera-Olivares, D., Burgos-Lefimil, C., Melendez, W., Flores-Prado, L. & Niemeyer, H. M. Chemical basis of nestmate recognition in a defense context in a one-piece nesting termite. Chemoecology 26(5), 163–172 (2016).
Article  Google Scholar 

35.
Kuswanto, E., Ahmad, I., Putra, R. E. & Harahap, I. S. Two novel volatile compounds as the key for intraspecific colony recognition in Macrotermes gilvus (Isoptera: Termitidae). J. Entomol. 12(2), 87–94 (2015).
CAS  Article  Google Scholar 

36.
Ismanto, A. & Baedowi, A. Efikasi ekstrak akar tuba dalam mengendalikan rayap tanah Macrotermes gilvus hagen pada pertanaman kayu putih. Jurnal Ecogreen. 5(1), 57–62 (2019).
Google Scholar 

37.
Jones, T. H. et al. The chemistry of exploding ants, Camponotus spp. (cylindricus complex). J. Chem. Ecol. 30(8), 1479–1492 (2004).
CAS  PubMed  Article  PubMed Central  Google Scholar 

38.
Laciny, A. et al. Colobopsis explodens sp. n., model species for studies on “exploding ants” (Hymenoptera, Formicidae), with biological notes and first illustrations of males of the Colobopsis cylindrica group. ZooKeys 751, 1–40 (2018).
Article  Google Scholar 

39.
Kuwahara, Y. Chemical Ecology of Astigmatid Mites (Cambridge University Press, Cambridge, 2004).
Google Scholar 

40.
Iqbal, N. & Saeed, S. Toxicity of six new chemical insecticides against the termite, Microtermes mycophagus D. (Isoptera: Termitidae: Macrotermitinae). Pak. J. Zool. 45(3), 709–713 (2013).
CAS  Google Scholar 

41.
Lihoreau, M. & Rivault, C. Kin recognition via cuticular hydrocarbons shapes cockroach social life. Behav. Ecol. 20, 46–53 (2009).
Article  Google Scholar 

42.
Emanuel, S. & Libersat, F. Nociceptive Pathway in the cockroach Periplaneta americana. Front. Physiol. https://doi.org/10.3389/fphys.2019.01100 (2019).
Article  PubMed  PubMed Central  Google Scholar 

43.
Deisig, N., Dupuy, F., Anton, S. & Renou, M. Responses to pheromones in a complex odor world: Sensory processing and behavior. Insects. 5(2), 399–422 (2014).
PubMed  PubMed Central  Article  Google Scholar 

44.
Nishino, H. et al. Spatial receptive fields for odor localization. Curr. Biol. 28(4), 600–608 (2018).
CAS  PubMed  Article  PubMed Central  Google Scholar 

45.
Sandoz, J. C., Pham, D. M., Renou, M. & Wadhams, L. Asymmetrical generalisation between pheromonal and floral odours in appetitive olfactory conditioning of the honey bee (Apis mellifera L.). J. Comp. Physiol. 187, 559–568 (2001).
CAS  Article  Google Scholar 

46.
Kreher, S. A., Kwon, J. Y. & Carlson, J. R. The molecular basis of odor coding in the Drosophila larva. Neuron 46(3), 445–456 (2005).
CAS  PubMed  Article  PubMed Central  Google Scholar 

47.
Grosjean, Y. et al. An olfactory receptor for food-derived odours promotes male courtship in Drosophila. Nature 478(7368), 236 (2011).
ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 

48.
Vosshall, L. B. & Hansson, B. S. A unified nomenclature system for the insect olfactory coreceptor. Chem. Senses 36(6), 497–498 (2011).
PubMed  Article  PubMed Central  Google Scholar 

49.
Peschke, K. & Eisner, T. Defensive secretion of the tenebrionid beetle, Blaps mucronata: Physical and chemical determinants of effectiveness. J. Comp. Physiol. 161(3), 377–388 (1987).
CAS  Article  Google Scholar 

50.
Dettner, K. Solvent-dependent variablity of effectiveness of quinone-defensive systems of Oxytelinae beetles (Coleoptera: Staphylinidae). Entomologia Generalis. 15, 275–292 (1991).
Article  Google Scholar 

51.
Roth, L. & Eisner, T. Chemical defenses of arthropods. Ann. Rev. Entomol. 7, 107–136 (2003).
Article  Google Scholar 

52.
Li, J. et al. Odoriferous defensive stink gland transcriptome to identify novel genes necessary for quinone synthesis in the red flour beetle Tribolium castaneum. PLoS Genet. 9(7), 1003596–1003596 (2013).
Article  CAS  Google Scholar 

53.
Delattre, O. et al. Complex alarm strategy in the most basal termite species. Behav. Ecol. Sociobiol. 69(12), 1945–1955 (2015).
Article  Google Scholar 

54.
Prestwich, G. D. & Chen, D. Soldier defense secretions of Trinervitermes bettonianus (Isoptera, Nasutitermitinae): Chemical variation in allopatric populations. J. Chem. Ecol. 7(1), 147–157 (1981).
CAS  PubMed  Article  Google Scholar 

55.
Piper, R. Extraordinary Animals: An Encyclopedia of Curious and Unusual Animals (Greenwood Publishing Group, Westport, 2007).
Google Scholar 

56.
Costa-Leonardo, A. A new interpretation of the defense glands of neotropical Ruptitermes (Isoptera, Termitidae, Apicotermitinae). Sociobiology 44, 391–402 (2004).
Google Scholar 

57.
Reinhard, J., Lacey, M. & Lenz, M. Application of the natural phagostimulant hydroquinone in bait systems for termite management (Isoptera). Sociobiology 39(2), 213–230 (2002).
Google Scholar 

58.
Hasyim, A., Istianto, M. & de Kogel, W. Male fruit fly, Bactrocera tau (Diptera; Tephritidae) attractants from Elsholtzia pubescens Bth. Asian J. Plant Sci. 6(1), 181–183 (2007).
Article  Google Scholar 

59.
Chen, Z. Y. et al. Insecticidal and repellent activity of essential oil from Amomum villosum Lour. and its main compounds against two stored-product insects. Int. J. Food Prop. 21(1), 2265–2275 (2018).
CAS  Article  Google Scholar 

60.
Reisenman, C. E., Lei, H. & Guerenstein, P. G. Neuroethology of olfactory-guided behavior and its potential application in the control of harmful insects. Front. Physiol. 7, 271–271 (2016).
PubMed  PubMed Central  Article  Google Scholar 

61.
Raina, A. K., Bland, J. M. & Osbrink, W. Hydroquinone is not a phagostimulant for the Formosan subterranean termite. J. Chem. Ecol. 31(3), 509–517 (2005).
CAS  PubMed  Article  Google Scholar 

62.
Bagnères, A.-G. & Hanus, R. Communication and social regulation in termites. In Social Recognition in Invertebrates: The Knowns and the Unknowns (eds Aquiloni, L. & Tricarico, E.) 193–248 (Springer, Cham, 2015).
Google Scholar 

63.
Alia Diyana, M. H., Appalasamy, S., Arumugam, N. & Boon, J. G. A study of a termite chemical defense fluid compound of Macrotermes carbonarius. IOP Conf. Ser. Earth Environ. Sci. https://doi.org/10.1088/1755-1315/269/1/012009 (2019).
Article  Google Scholar 

64.
Environmental Protection Agency. Furanone. Prevention P A T S (National Center for Environmental Publications and Information, 1993).

65.
Igwe, O. U. & Udofia, D. E. Secondary metabolites of the cuticular abdominal glands of variegated grasshopper (Zonocerus variegatus L.). Int. J. Spectrosc. 2015, 1–6 (2015).
Article  CAS  Google Scholar 

66.
Neoh, K. B., Yeap, B.-K., Tsunoda, K., Yoshimura, T. & Lee, C.-Y. Do termites avoid carcasses? Behavioral responses depend on the nature of the carcasses. PLoS ONE 7(4), 36375 (2012).
ADS  Article  CAS  Google Scholar 

67.
Blassioli-Moraes, M. C., Laumann, R. A., Michereff, M. F. F. & Borges, M. Semiochemicals for integrated pest management. In Sustainable Agrochemistry: A Compendium of Technologies (ed. Vaz, S., Jr.) 85–112 (Springer, Cham, 2019).
Google Scholar 

68.
de Melo, A. R. et al. Toxicity of different fatty acids and methyl esters on Culex quinquefasciatus larvae. Ecotoxicol. Environ. Saf. 154, 1–5 (2018).
ADS  PubMed  Article  CAS  Google Scholar 

69.
Xie, Y., Wang, K., Huang, Q. & Lei, C. Evaluation toxicity of monoterpenes to subterranean termite, Reticulitermes chinensis Snyder. Ind. Crops Prod. 53, 163–166 (2014).
CAS  Article  Google Scholar 

70.
Xie, Y. et al. Antitermitic and antifungal activities of eugenol and its congeners from the flower buds of Syzgium aromaticum (clove). Ind. Crops Prod. 77, 780–786 (2015).
CAS  Article  Google Scholar 

71.
Zhang, Z., Yang, T., Zhang, Y., Wang, L. & Xie, Y. Fumigant toxicity of monoterpenes against fruitfly, Drosophila melanogaster. Ind. Crops Prod. 81, 147–151 (2016).
CAS  Article  Google Scholar 

72.
Silva, L. N. D. et al. The influence of fatty acid methyl esters (FAMEs) in the biochemistry and the Na+/K+-ATPase activity of Culex quinquefasciatus Larvae. J. Membr. Biol. 249(4), 459–467 (2016).
CAS  PubMed  Article  PubMed Central  Google Scholar  More

1.
IPCC. Managing the risks of extreme events and disasters to advance climate change adaptation. A special report of working groups i and ii of the intergovernmental panel on climate change. Cambridge, UK, and New York, NY, USA: Cambridge University Press; 2012.
Google Scholar 
2.
Trenberth KE, Dai AG, van der Schrier G, Jones PD, Barichivich J, Briffa KR, et al. Global warming and changes in drought. Nat Clim Chang. 2014;4:17–22.
Article  Google Scholar 

3.
Kreyling J, Henry HAL. Vanishing winters in germany: soil frost dynamics and snow cover trends, and ecological implications. Clim Res. 2011;46:269–76.
Article  Google Scholar 

4.
Congreves KA, Wagner-Riddle C, Si BC, Clough TJ. Nitrous oxide emissions and biogeochemical responses to soil freezing-thawing and drying-wetting. Soil Biol Biochem. 2018;117:5–15.
CAS  Article  Google Scholar 

5.
Schimel J, Balser TC, Wallenstein M. Microbial stress-response physiology and its implications for ecosystem function. Ecology 2007;88:1386–94.
PubMed  Article  Google Scholar 

6.
Tecon R, Or D. Biophysical processes supporting the diversity of microbial life in soil. FEMS Microbiol Rev. 2017;41:599–623.
CAS  PubMed  PubMed Central  Article  Google Scholar 

7.
Koponen HT, Bååth E. Soil bacterial growth after a freezing/thawing event. Soil Biol Biochem. 2016;100:229–32.
CAS  Article  Google Scholar 

8.
Schostag M, Priemé A, Jacquiod S, Russel J, Ekelund F, Jacobsen CS. Bacterial and protozoan dynamics upon thawing and freezing of an active layer permafrost soil. ISME J. 2019;13:1345–59.
CAS  PubMed  PubMed Central  Article  Google Scholar 

9.
Barnard RL, Osborne CA, Firestone MK. Changing precipitation pattern alters soil microbial community response to wet-up under a mediterranean-type climate. ISME J. 2015;9:946–57.
CAS  PubMed  Article  Google Scholar 

10.
Drigo B, Nielsen UN, Jeffries TC, Curlevski NJA, Singh BK, Duursma RA, et al. Interactive effects of seasonal drought and elevated atmospheric carbon dioxide concentration on prokaryotic rhizosphere communities. Environ Microbiol. 2017;19:3175–85.
CAS  PubMed  Article  PubMed Central  Google Scholar 

11.
Priemé A, Christensen S. Natural perturbations, drying-wetting and freezing-thawing cycles, and the emission of nitrous oxide, carbon dioxide and methane from farmed organic soils. Soil Biol Biochem. 2001;33:2083–91.
Article  Google Scholar 

12.
Kim DG, Vargas R, Bond-Lamberty B, Turetsky MR. Effects of soil rewetting and thawing on soil gas fluxes: a review of current literature and suggestions for future research. Biogeosciences 2012;9:2459–83.
CAS  Article  Google Scholar 

13.
Meisner A, Leizeaga A, Rousk J, Bååth E. Partial drying accelerates bacterial growth recovery to rewetting. Soil Biol Biochem. 2017;112:269–76.
CAS  Article  Google Scholar 

14.
Blazewicz SJ, Schwartz E, Firestone MK. Growth and death of bacteria and fungi underlie rainfall-induced carbon dioxide pulses from seasonally dried soil. Ecology. 2014;95:1162–72.
PubMed  Article  PubMed Central  Google Scholar 

15.
Barnard RL, Osborne CA, Firestone MK. Responses of soil bacterial and fungal communities to extreme desiccation and rewetting. ISME J. 2013;7:2229–41.
CAS  PubMed  PubMed Central  Article  Google Scholar 

16.
Meisner A, Jacquiod S, Snoek BL, ten Hooven FC, van der Putten WH. Drought legacy effects on the composition of soil fungal and prokaryote communities. Front Microbiol. 2018;9:294.
PubMed  PubMed Central  Article  Google Scholar 

17.
de Vries FT, Griffiths RI, Bailey M, Craig H, Girlanda M, Gweon HS, et al. Soil bacterial networks are less stable under drought than fungal networks. Nat Commun. 2018;9:3033.
PubMed  PubMed Central  Article  CAS  Google Scholar 

18.
Bouskill NJ, Lim HC, Borglin S, Salve R, Wood TE, Silver WL, et al. Pre-exposure to drought increases the resistance of tropical forest soil bacterial communities to extended drought. ISME J. 2013;7:384–94.
CAS  PubMed  Article  Google Scholar 

19.
Ren J, Song C, Hou A, Song Y, Zhu X, Cagle GA. Shifts in soil bacterial and archaeal communities during freeze-thaw cycles in a seasonal frozen marsh, northeast china. Sci Total Environ. 2018;625:782–91.
CAS  PubMed  Article  Google Scholar 

20.
Mannisto MK, Tiirola M, Haggblom MM. Effect of freeze-thaw cycles on bacterial communities of arctic tundra soil. Micro Ecol. 2009;58:621–31.
Article  Google Scholar 

21.
Evans SE, Wallenstein MD. Climate change alters ecological strategies of soil bacteria. Ecol Lett. 2014;17:155–64.
PubMed  Article  PubMed Central  Google Scholar 

22.
Fuchslueger L, Bahn M, Hasibeder R, Kienzl S, Fritz K, Schmitt M, et al. Drought history affects grassland plant and microbial carbon turnover during and after a subsequent drought event. J Ecol. 2016;104:1453–65.
CAS  PubMed  PubMed Central  Article  Google Scholar 

23.
Hawkes CV, Waring BG, Rocca JD, Kivlin SN. Historical climate controls soil respiration responses to current soil moisture. Proc Natl Acad Sci. 2017;114:6322–7.
CAS  PubMed  Article  PubMed Central  Google Scholar 

24.
Preece C, Verbruggen E, Liu L, Weedon JT, Peñuelas J. Effects of past and current drought on the composition and diversity of soil microbial communities. Soil Biol Biochem. 2019;131:28–39.
CAS  Article  Google Scholar 

25.
de Nijs EA, Hicks LC, Leizeaga A, Tietema A, Rousk J. Soil microbial moisture dependences and responses to drying–rewetting: the legacy of 18 years drought. Glob Change Biol. 2019;25:1005–15.
Article  Google Scholar 

26.
Evans SE, Wallenstein MD. Soil microbial community response to drying and rewetting stress: does historical precipitation regime matter? Biogeochemistry 2012;109:101–16.
Article  Google Scholar 

27.
Butterly CR, Bunemann EK, McNeill AM, Baldock JA, Marschner P. Carbon pulses but not phosphorus pulses are related to decreases in microbial biomass during repeated drying and rewetting of soils. Soil Biol Biochem. 2009;41:1406–16.
CAS  Article  Google Scholar 

28.
Schimel JP, Clein JS. Microbial response to freeze-thaw cycles in tundra and taiga soils. Soil Biol Biochem. 1996;28:1061–6.
Article  Google Scholar 

29.
Jurburg SD, Nunes I, Brejnrod A, Jacquiod S, Priemé A, Sørensen SJ, et al. Legacy effects on the recovery of soil bacterial communities from extreme temperature perturbation. Front Microbiol. 2017;8:1832.
PubMed  PubMed Central  Article  Google Scholar 

30.
Calderón K, Philippot L, Bizouard F, Breuil M-C, Bru D, Spor A. Compounded disturbance chronology modulates the resilience of soil microbial communities and n-cycle related functions. Front Microbiol. 2018;9:2721.
PubMed  PubMed Central  Article  Google Scholar 

31.
Wilson SL, Grogan P, Walker VK. Prospecting for ice association: characterization of freeze–thaw selected enrichment cultures from latitudinally distant soils. Can J Microbiol. 2012;58:402–12.
CAS  PubMed  Article  Google Scholar 

32.
Evans SE, Wallenstein MD, Burke IC. Is bacterial moisture niche a good predictor of shifts in community composition under long-term drought? Ecology. 2014;95:110–22.
PubMed  Article  Google Scholar 

33.
Yergeau E, Kowalchuk GA. Responses of antarctic soil microbial communities and associated functions to temperature and freeze-thaw cycle frequency. Environ Microbiol. 2008;10:2223–35.
PubMed  Article  Google Scholar 

34.
Lennon JT, Aanderud ZT, Lehmkuhl BK, Schoolmaster DR. Mapping the niche space of soil microorganisms using taxonomy and traits. Ecology. 2012;93:1867–79.
PubMed  Article  Google Scholar 

35.
Mackey BM, Derrick CM. Conductance measurements of the lag phase of injured salmonella-typhimurium. J Appl Bacteriol. 1984;57:299–308.
CAS  PubMed  Article  Google Scholar 

36.
Nocker A, Fernandez PS, Montijn R, Schuren F. Effect of air drying on bacterial viability: a multiparameter viability assessment. J Microbiol Methods. 2012;90:86–95.
CAS  PubMed  Article  Google Scholar 

37.
Potts M. Desiccation tolerance of prokaryotes. Microbiological Rev. 1994;58:755–805.
CAS  Article  Google Scholar 

38.
Öquist MG, Sparrman T, Klemedtsson L, Drotz SH, Grip H, Schleucher J, et al. Water availability controls microbial temperature responses in frozen soil co2 production. Glob Change Biol. 2009;15:2715–22.
Article  Google Scholar 

39.
Manzoni S, Schimel JP, Porporato A. Responses of soil microbial communities to water-stress: results from a meta-analysis. Ecology. 2012;93:930–8.
PubMed  Article  Google Scholar 

40.
Iovieno P, Bååth E. Effect of drying and rewetting on bacterial growth rates in soil. FEMS Microbiol Ecol. 2008;65:400–7.
CAS  PubMed  Article  Google Scholar 

41.
Drotz SH, Sparrman T, Nilsson MB, Schleucher J, Öquist MG. Both catabolic and anabolic heterotrophic microbial activity proceed in frozen soils. Proc Natl Acad Sci USA. 2010;107:21046–51.
CAS  PubMed  Article  Google Scholar 

42.
Lennon JT, Jones SE. Microbial seed banks: the ecological and evolutionary implications of dormancy. Nat Rev Microbiol. 2011;9:119–30.
CAS  PubMed  Article  Google Scholar 

43.
Manzoni S, Schaeffer SM, Katul G, Porporato A, Schimel JP. A theoretical analysis of microbial eco-physiological and diffusion limitations to carbon cycling in drying soils. Soil Biol Biochem. 2014;73:69–83.
CAS  Article  Google Scholar 

44.
Larsen KS, Jonasson S, Michelsen A. Repeated freeze-thaw cycles and their effects on biological processes in two arctic ecosystem types. Appl Soil Ecol. 2002;21:187–95.
Article  Google Scholar 

45.
Williams MA, Xia K. Characterization of the water soluble soil organic pool following the rewetting of dry soil in a drought-prone tallgrass prairie. Soil Biol Biochem. 2009;41:21–8.
CAS  Article  Google Scholar 

46.
Placella SA, Brodie EL, Firestone MK. Rainfall-induced carbon dioxide pulses result from sequential resuscitation of phylogenetically clustered microbial groups. Proc Natl Acad Sci USA. 2012;109:10931–6.
CAS  PubMed  Article  Google Scholar 

47.
Fukami T. Historical contingency in community assembly: Integrating niches, species pools, and priority effects. Annu Rev Ecol Evol Syst. 2015;46:1–23.
Article  Google Scholar 

48.
Szekely AJ, Langenheder S. Dispersal timing and drought history influence the response of bacterioplankton to drying-rewetting stress. ISME J. 2017;11:1764–76.
PubMed  PubMed Central  Article  Google Scholar 

49.
Engelhardt IC, Welty A, Blazewicz SJ, Bru D, Rouard N, Breuil M-C, et al. Depth matters: effects of precipitation regime on soil microbial activity upon rewetting of a plant-soil system. ISME J. 2018;12:1061–71.
CAS  PubMed  PubMed Central  Article  Google Scholar 

50.
Guhr A, Borken W, Spohn M, Matzner E. Redistribution of soil water by a saprotrophic fungus enhances carbon mineralization. Proc Natl Acad Sci USA. 2015;112:14647–51.
CAS  PubMed  Article  PubMed Central  Google Scholar 

51.
Yergeau E, Hogues H, Whyte LG, Greer CW. The functional potential of high arctic permafrost revealed by metagenomic sequencing, qpcr and microarray analyses. ISME J. 2010;4:1206–14.
CAS  PubMed  Article  PubMed Central  Google Scholar 

52.
Gittel A, Barta J, Kohoutova I, Mikutta R, Owens S, Gilbert J, et al. Distinct microbial communities associated with buried soils in the siberian tundra. ISME J. 2014;8:841–53.
CAS  PubMed  Article  PubMed Central  Google Scholar 

53.
Hultman J, Waldrop MP, Mackelprang R, David MM, McFarland J, Blazewicz SJ, et al. Multi-omics of permafrost, active layer and thermokarst bog soil microbiomes. Nature 2015;521:208–12.
CAS  PubMed  Article  PubMed Central  Google Scholar 

54.
Feng X, Nielsen LL, Simpson MJ. Responses of soil organic matter and microorganisms to freeze–thaw cycles. Soil Biol Biochem. 2007;39:2027–37.
CAS  Article  Google Scholar 

55.
Carini P, Marsden PJ, Leff J, Morgan EE, Strickland MS, Fierer N. Relic DNA is abundant in soil and obscures estimates of soil microbial diversity. Nat Microbiol. 2017;2:6.
Article  CAS  Google Scholar 

56.
Lennon JT, Muscarella ME, Placella SA, Lehmkuhl BK. How, when, and where relic DNA affects microbial diversity. mBio 2018;9:e00637–18.
PubMed  PubMed Central  Google Scholar 

57.
Blazewicz SJ, Barnard RL, Daly RA, Firestone MK. Evaluating rrna as an indicator of microbial activity in environmental communities: limitations and uses. ISME J. 2013;7:2061–8.
CAS  PubMed  PubMed Central  Article  Google Scholar 

58.
Blagodatskaya E, Kuzyakov Y. Active microorganisms in soil: critical review of estimation criteria and approaches. Soil Biol Biochem. 2013;67:192–211.
CAS  Article  Google Scholar 

59.
van de Voorde TFJ, van der Putten WH, Martijn Bezemer T. Intra- and interspecific plant–soil interactions, soil legacies and priority effects during old-field succession. J Ecol. 2011;99:945–53.
Article  Google Scholar 

60.
Nunes I, Jurburg S, Jacquiod S, Brejnrod A, Falcão Salles J, Priemé A, et al. Soil bacteria show different tolerance ranges to an unprecedented disturbance. Biol Fertil Soils. 2018;54:189–202.
Article  Google Scholar 

61.
Christiansen CT, Haugwitz MS, Priemé A, Nielsen CS, Elberling B, Michelsen A, et al. Enhanced summer warming reduces fungal decomposer diversity and litter mass loss more strongly in dry than in wet tundra. Glob Change Biol. 2017;23:406–20.
Article  Google Scholar 

62.
Borg Dahl M, Brejnrod AD, Russel J, Sørensen SJ, Schnittler M. Different degrees of niche differentiation for bacteria, fungi, and myxomycetes within an elevational transect in the german alps. Micro Ecol. 2019;78:764–80.
CAS  Article  Google Scholar 

63.
RCoreTeam. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2018.
Google Scholar 

64.
McMurdie PJ, Holmes S. Phyloseq: an r package for reproducible interactive analysis and graphics of microbiome census data. PLOS ONE. 2013;8:e61217.
CAS  PubMed  PubMed Central  Article  Google Scholar 

65.
Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. Vegan: community ecology package. R package version 2.5-5. https://CRAN.R-project.org/package=vegan. 2019.

66.
Kuznetsova A, Brockhoff PB, Christensen RHB. lmerTest Package: Tests in Linear Mixed Effects Models. J Stat Softw. 2017;82:1–26.
Article  Google Scholar 

67.
Bates D, Maechler M, Bolker B, Walker S. Fitting Linear Mixed-Effects Models Using lme4. J Stat Softw. 2015;67:1–48.
Article  Google Scholar 

68.
Warnes GR, Bolker B, Bonebakker L, Gentleman R, Liaw WHA, Lumley T, et al. Package ‘gplots’: various r programming tools for plotting data. R package version 3.0.1. https://CRAN.R-project.org/package=gplots. 2019.

69.
Wickham H. Ggplot2: elegant graphics for data analysis. New York: Springer-Verlag; 2016.
Google Scholar 

70.
Sharma S, Szele Z, Schilling R, Munch JC, Schloter M. Influence of freeze-thaw stress on the structure and function of microbial communities and denitrifying populations in soil. Appl Environ Microbiol. 2006;72:2148–54.
CAS  PubMed  PubMed Central  Article  Google Scholar 

71.
Tilston EL, Sparrman T, Öquist MG. Unfrozen water content moderates temperature dependence of sub-zero microbial respiration. Soil Biol Biochem. 2010;42:1396–407.
CAS  Article  Google Scholar 

72.
Ho A, Di Lonardo DP, Bodelier PLE. Revisiting life strategy concepts in environmental microbial ecology. FEMS Microbiol Ecol. 2017;93:14.
Google Scholar 

73.
Yuste JC, Penuelas J, Estiarte M, Garcia-Mas J, Mattana S, Ogaya R, et al. Drought-resistant fungi control soil organic matter decomposition and its response to temperature. Glob Change Biol. 2011;17:1475–86.
Article  Google Scholar 

74.
Kaisermann A, Maron PA, Beaumelle L, Lata JC. Fungal communities are more sensitive indicators to non-extreme soil moisture variations than bacterial communities. Appl Soil Ecol. 2015;86:158–64.
Article  Google Scholar 

75.
Clemmensen KE, Bahr A, Ovaskainen O, Dahlberg A, Ekblad A, Wallander H, et al. Roots and associated fungi drive long-term carbon sequestration in boreal forest. Science. 2013;339:1615–8.
CAS  PubMed  Article  Google Scholar 

76.
Bates ST, Clemente JC, Flores GE, Walters WA, Parfrey LW, Knight R, et al. Global biogeography of highly diverse protistan communities in soil. ISME J. 2013;7:652–9.
CAS  PubMed  Article  Google Scholar 

77.
Geisen S, Mitchell EAD, Adl S, Bonkowski M, Dunthorn M, Ekelund F, et al. Soil protists: a fertile frontier in soil biology research. FEMS Microbiol Rev. 2018;42:293–323.
CAS  PubMed  Article  Google Scholar 

78.
Geisen S, Bandow C, Römbke J, Bonkowski M. Soil water availability strongly alters the community composition of soil protists. Pedobiologia 2014;57:205–13.
Article  Google Scholar 

79.
Clarholm M. Protozoan grazing of bacteria in soil—impact and importance. Micro Ecol. 1981;7:343–50.
CAS  Article  Google Scholar 

80.
Harder CB, Rønn R, Brejnrod A, Bass D, Al-Soud WA, Ekelund F. Local diversity of heathland cercozoa explored by in-depth sequencing. ISME J. 2016;10:2488.
CAS  PubMed  PubMed Central  Article  Google Scholar 

81.
Ekelund F, Frederiksen HB, Rønn R. Population dynamics of active and total ciliate populations in arable soil amended with wheat. Appl Environ Microbiol. 2002;68:1096–101.
CAS  PubMed  PubMed Central  Article  Google Scholar 

82.
Müller H, Achilles-Day UEM, Day JG. Tolerance of the resting cysts of colpoda inflata (ciliophora, colpodea) and meseres corlissi (ciliophora, spirotrichea) to desiccation and freezing. Eur J Protistol. 2010;46:133–42.
PubMed  Article  Google Scholar 

83.
Stieglmeier M, Alves RJE, Schleper C. The phylum thaumarchaeota. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. The prokaryotes: other major lineages of bacteria and the archaea. Berlin, Heidelberg: Springer Berlin Heidelberg; 2014. p. 347–62.

84.
Shi Y, Adams JM, Ni Y, Yang T, Jing X, Chen L, et al. The biogeography of soil archaeal communities on the eastern tibetan plateau. Sci Rep. 2016;6:38893.
CAS  PubMed  PubMed Central  Article  Google Scholar 

85.
Pedrós-Alió C. The rare bacterial biosphere. Ann Rev Mar Sci. 2012;4:449–66.
PubMed  Article  Google Scholar 

86.
Lynch MDJ, Neufeld JD. Ecology and exploration of the rare biosphere. Nat Rev Microbiol. 2015;13:217–29.
CAS  PubMed  Article  Google Scholar 

87.
Chowdhury N, Nakatani AS, Setia R, Marschner P. Microbial activity and community composition in saline and non-saline soils exposed to multiple drying and rewetting events. Plant Soil. 2011;348:103–13.
CAS  Article  Google Scholar 

88.
Fierer N, Schimel JP. Effects of drying-rewetting frequency on soil carbon and nitrogen transformations. Soil Biol Biochem. 2002;34:777–87.
CAS  Article  Google Scholar 

89.
Schöler A, Jacquiod S, Vestergaard G, Schulz S, Schloter M. Analysis of soil microbial communities based on amplicon sequencing of marker genes. Biol Fertil Soils. 2017;53:485–9.
Article  CAS  Google Scholar 

90.
Lau JA, Lennon JT. Rapid responses of soil microorganisms improve plant fitness in novel environments. Proc Natl Acad Sci USA. 2012;109:14058–62.
CAS  PubMed  Article  PubMed Central  Google Scholar 

91.
Meisner A, De Deyn GB, de Boer W, van der Putten WH. Soil biotic legacy effects of extreme weather events influence plant invasiveness. Proc Natl Acad Sci USA. 2013;110:9835–8.
CAS  PubMed  Article  PubMed Central  Google Scholar 

92.
Solyanikova IP, Suzina NE, Egozarian NS, Polivtseva VN, Prisyazhnaya NV, El-Registan GI, et al. The response of soil arthrobacter agilis lush13 to changing conditions: Transition between vegetative and dormant state. J Environ Sci Health B 2017;52:745–51.
CAS  PubMed  Article  PubMed Central  Google Scholar 

93.
Levinson HS, Hyatt MT. Correlation of respiratory activity with phases of spore germination and growth in bacillus-megaterium as influenced by manganese and l-alanine. J Bacteriol. 1956;72:176–83.
CAS  PubMed  PubMed Central  Article  Google Scholar 

94.
Williamson KE, Fuhrmann JJ, Wommack KE, Radosevich M. Viruses in soil ecosystems: an unknown quantity within an unexplored territory. Annu Rev Virol. 2017;4:201–19.
CAS  PubMed  Article  PubMed Central  Google Scholar 

95.
van Kruistum H, Bodelier PLE, Ho A, Meima-Franke M, Veraart AJ. Resistance and recovery of methane-oxidizing communities depends on stress regime and history; a microcosm study. Front Microbiol. 2018;9:1714.
PubMed  PubMed Central  Article  Google Scholar 

96.
Bérard A, Ben Sassi M, Renault P, Gros R. Severe drought-induced community tolerance to heat wave. An experimental study on soil microbial processes. J Soils Sed. 2012;12:513–8.
Article  Google Scholar  More