Philippa Kaur

1.
Zenner, E. Does old-growth condition imply high live-tree structural complexity?. For. Ecol. Manag. 195, 243–258 (2004).
Article  Google Scholar 
2.
Forest Ecosystem Management Assessment Team (FEMAT). Draft Supplemental Environmental Impact Statement on Management of Habitat for Late Successional and Oldgrowth Forest Related Species within the Range of the Northern Spotted Owl (US Government Printing Office, Washington, DC, 1993).
Google Scholar 

3.
Wan, P. et al. Impacts of different forest management methods on the stand spatial structure of a natural Quercus aliena var. acuteserrata forest in Xiaolongshan, China. Ecol. Inform. 50, 86–94 (2019).
Article  Google Scholar 

4.
Carrer, M., Castagneri, D., Popa, I., Pividori, M. & Lingua, E. Tree spatial patterns and stand attributes in temperate forests: The importance of plot size, sampling design, and null model. For. Ecol. Manag. 407, 125–134 (2018).
Article  Google Scholar 

5.
Bauhus, J., Puettmann, K. & Messier, C. Silviculture for old-growth attributes. For. Ecol. Manag. 258, 525–537 (2009).
Article  Google Scholar 

6.
Messier, C., Puettmann, K. J. & Coates, D. K. Managing Forests as Complex Adaptive Systems: Building Resilience to the Challenge of Global Change (Routledge, Abingdon, 2013).
Google Scholar 

7.
McElhinny, C., Gibbons, P., Brack, C. & Bauhus, J. Forest and woodland stand structural complexity: Its definition and measurement. For. Ecol. Manage. 218, 1–24 (2005).
Article  Google Scholar 

8.
Di Filippo, A., Biondi, F., Piovesan, G. & Ziaco, E. Tree ring-based metrics for assessing old-growth forest naturalness. J. Appl. Ecol. 54, 737–749 (2017).
Article  Google Scholar 

9.
Parrotta, J. A., Turnbull, J. W. & Jones, N. Catalyzing native forest regeneration on degraded tropical lands. For. Ecol. Manag. 99, 1–7 (1997).
Article  Google Scholar 

10.
Neumann, M. & Starlinger, F. The significance of different indices for stand structure and diversity in forests. For. Ecol. Manag. 145, 91–106 (2001).
Article  Google Scholar 

11.
McCleary, K. & Mowat, G. Using forest structural diversity to inventory habitat diversity of forest-dwelling wildlife in the West Kootenay region of British Columbia 2 1–13 (2002).

12.
Ishii, H. T., Tanabe, S.-I. & Hiura, T. Exploring the relationships among canopy structure, stand productivity, and biodiversity of temperate forest ecosystems. For. Sci. 50, 342–355 (2004).
Google Scholar 

13.
Tews, J. et al. Animal species diversity driven by habitat heterogeneity/diversity: the importance of keystone structures. J. Biogeogr. 31, 79–92 (2004).
Article  Google Scholar 

14.
Long, J. N. & Shaw, J. D. The influence of compositional and structural diversity on forest productivity. Forestry 83, 121–128 (2010).
Article  Google Scholar 

15.
Dănescu, A., Albrecht, A. & Bauhus, J. Structural diversity promotes productivity of mixed, uneven-aged forests in southwestern Germany. Oecologia 182, 319–333 (2016).
ADS  PubMed  Article  Google Scholar 

16.
Ehbrecht, M., Schall, P., Ammer, C. & Seidel, D. Quantifying stand structural complexity and its relationship with forest management, tree species diversity and microclimate. Agric. For. Meteorol. 242, 1–9 (2017).
ADS  Article  Google Scholar 

17.
Zenner, E. K. Do residual trees increase structural complexity in pacific northwest?. Ecol. Appl. 10, 800–810 (2000).
Article  Google Scholar 

18.
Hardiman, B. S., Bohrer, G., Gough, C. M., Vogel, C. S. & Curtisi, P. S. The role of canopy structural complexity in wood net primary production of a maturing northern deciduous forest. Ecology 92, 1818–1827 (2011).
PubMed  Article  Google Scholar 

19.
Puettmann, K. J., Coates, K. D. & Messier, C. C. A Critique of Silviculture: Managing for Complexity (Island Press, Washington, D.C., 2012).
Google Scholar 

20.
Robertson, G. P. & Tiedje, J. Spatial variability in a successional plant community: patterns of nitrogen availability. Ecology 69, 0–1524 (1988).

21.
Palmer, M. W. Spatial scale and patterns of species-environment relationships in hardwood forest of the North Carolina piedmont. Coenoses, 79–87 (1990).

22.
Lechowicz, M. & Bell, G. The ecology and genetics of fitness in forest plants. II. Microspatial heterogeneity of the edaphic environment. J. Ecol. 79, 687 (1991).

23.
Song, B. et al. Modeling canopy structure and heterogeneity across scales: from crowns to canopy. For. Ecol. Manage. 96, 217–229 (1997).
Article  Google Scholar 

24.
Zenner, E. & Peck, J. Characterizing structural conditions in mature managed red pine: spatial dependency of metrics and adequacy of plot size. For. Ecol. Manag. 257, 311–320 (2009).
Article  Google Scholar 

25.
Pommerening, A. & Uria-Diez, J. Do large forest trees tend towards high species mingling? Ecol. Inform. 42 (2017).

26.
Wang, H., Peng, H., Hui, G., Hu, Y. & Zhao, Z. Large trees are surrounded by more heterospecific neighboring trees in Korean pine broad-leaved natural forests. Sci. Rep. 8, 9149 (2018).
ADS  PubMed  PubMed Central  Article  CAS  Google Scholar 

27.
Hubbell, S. P., Ahumada, J. A., Condit, R. & Foster, R. B. Local neighborhood effects on long-term survival of individual trees in a neotropical forest. Ecol. Res. 16, 859–875 (2001).
Article  Google Scholar 

28.
Stoll, P. & Newbery, D. M. Evidence of species-specific neighborhood effects in the dipterocarpaceae of a bornean rain forest. Ecology 86, 3048–3062 (2005).
Article  Google Scholar 

29.
Pillay, T. & Ward, D. Spatial pattern analysis and competition between Acacia karroo trees in humid savannas. Plant Ecol. 213 (2012).

30.
Fueldner, K., Sattler, S., Zucchini, W. & Gadow, K. V. Modelling person-specific tree selection probabilities in a thinning. Allgemeine Forst Und Jagdzeitung (1996).

31.
Zenner, E. & Hibbs, D. A new method for modeling the heterogeneity of forest structure. For. Ecol. Manag. 129 (2000).

32.
Pommerening, A. Approaches to quantifying forest structures. Forestry 75(3), 305–324 (2002).
Article  Google Scholar 

33.
Beckschäfer, P. et al. Enhanced structural complexity index: an improved index for describing forest structural complexity. Open J. For. 3, 23–29 (2013).
Google Scholar 

34.
Kint, V., van Meirvenne, M., Nachtergale, L., Geudens, G. & Lust, N. Spatial methods for quantifying forest stand structure development: a comparison between nearest-neighbor indices and variogram analysis. For. Sci. 49, 36–49 (2003).
Google Scholar 

35.
Clark, P. J. & Evans, F. C. Distance to nearest neighbor as a measure of spatial relationships in populations. Ecology 35, 445–453 (1954).
Article  Google Scholar 

36.
Ripley, B. D. Spatial Statistics (Wiley, New York, 1981).

37.
Ripley, B. D. Modelling spatial patterns. J. R. Stat. Soc. 39(2), 172–212 (1977).
MathSciNet  Google Scholar 

38.
Pommerening, A. & Grabarnik, P. Individual-Based Methods in Forest Ecology and Management (Springer, Berlin, 2019).

39.
Gadow, K., Albert, M. & Hui, G. Das Winkelmaß – ein Strukturparameter zur beschreibung der Individualverteilung in Waldbeständen. Centralblatt für das gesamte Forstwesen 115(1), 1–10 (1998).
Google Scholar 

40.
Aguirre, O., Hui, G., Gadow, K. v. & Jiménez, J. An analysis of spatial forest structure using neighbourhood-based variables. For. Ecol. Manag. 183, 137–145 (2003).

41.
Hui, G. & Gadow, K. Das Winkelmass – Theoretische Überlegungen zum optimalen Standardwinkel. Allgemeine Forst u. Jagdzeitung 173(9), 173–177 (2002).
Google Scholar 

42.
Pommerening, A. Evaluating structural indices by reversing forest structural analysis. For. Ecol. Manage. 224, 266–277 (2006).
Article  Google Scholar 

43.
Li, Y., Hui, G., Zhao, Z., Hu, Y. & Adler, P. The bivariate distribution characteristics of spatial structure in natural Korean pine broad-leaved forest. Journal of Vegetation Science 23 (2012).

44.
Graz, F. P. Spatial diversity of dry savanna woodlands. Assessing the spatial diversity of a dry savanna woodland stand in northern Namibia using neighbourhood-based measures. Biodivers. Conserv. 00, 1–16 (2004).

45.
Pastorella, F. & Paletto, A. Stand structure indices as tools to support forest management: an application in Trentino forests (Italy). J. For. Sci. 59, 159–168 (2013).
Article  Google Scholar 

46.
Zhao, Z. et al. Testing the significance of different tree spatial distribution patterns based on the Uniform Angle Index. Can. J. For. Res. 44(11), 1417–1425 (2014).
Article  Google Scholar 

47.
Zhang, G. et al. Composition of basal area in natural forests based on the uniform angle index. Ecol. Inform. 45, 1–8 (2018).
Article  Google Scholar 

48.
Stiell, W. How uniformity of tree distribution affects stand growth. For. Chron. 54, 156–158 (1978).
Article  Google Scholar 

49.
Jay, A., Nichols, J. & Vanclay, J. Social and ecological issues for private native forestry in north-eastern New South Wales Australia. Small Scale For. 6, 115–126 (2007).
Article  Google Scholar 

50.
Zhang, G. et al. Designing near-natural planting patterns for plantation forests in China. For. Ecosyst. 6, 137 (2019).
Article  Google Scholar 

51.
Moeur, M. Characterising spatial patterns of trees using stem-mapped data. For. Sci. 39, 756–775 (1993).
ADS  Google Scholar 

52.
Stohlgren, T. Spatial patterns of giant sequoia (Sequoiadendrongiganteum) in two sequoia groves in Sequoia National Park California. Can. J. For. Res. 23, 120–132 (2011).
Article  Google Scholar 

53.
Pommerening, A. & Grabarnik, P. Individual-based Methods in Forest Ecology and Management (2019).

54.
Clark, P. & Evans, F. Distance to nearest neighbor as a measure of spatial relations. Ecology 35, 445–453 (1954).
Article  Google Scholar 

55.
Assunçáo, R. Testing spatial randomness by means of angle. Biometrics 50, 531–537 (1994).
MATH  Article  Google Scholar 

56.
Corral-Rivas JJ. PhD thesis. University of Göttingen (2006).

57.
Hui, G., Zhang, G., Zhao, Z. & Yang, A. Methods of forest structure research: a review. Curr. For. Rep. 5(3), 142–154. https://doi.org/10.1007/S40725-019-00090-7 (2019).
Article  Google Scholar 

58.
Gadow, K., Hui, G. & Albert, M. Das Winkelmaß – Ein Strukturparameter zur Beschreibung der Individualverteilung in Waldbeständen. Centralblatt für das Gesamte Forstwesen 115, 1–10 (1998).
Google Scholar 

59.
Wang, H. et al. The influence of sampling unit size and spatial arrangement patterns on neighborhood-based spatial structure analyses of forest stands. For. Syst. 25, e056 (2016).
Google Scholar 

60.
Kraft, G. Beiträge zur Lehre von den Durchforstungen, Schlagstellungen und Lichtungshieben, Vol. 154 (Klindworth’s Verlag, Hanover, 1884).

61.
Röhrig, E. & Gussone, H. A. Waldbau auf Ökologischer Grundlage: Zweiter Band (Hamburg, Paul Parey, 1982).
Google Scholar 

62.
Hawley, R. C. & Smith, M. D. The practice of silviculture. Ecology 17(1), 172 (1936).
Article  Google Scholar 

63.
Larsen, J. B. & Nielsen, A. B. Nature-based forest management—Where are we going?. For. Ecol. Manag. 238, 107–117 (2007).
Article  Google Scholar 

64.
Ajani, J. The Forest Wars (Melbourne University, Melbourne, 2007).
Google Scholar 

65.
Nichols, J. D., Bristow, M. & Vanclay, J. K. Mixed-species plantations: prospects and challenges. For. Ecol. Manag. 233, 383–390 (2006).
Article  Google Scholar 

66.
Carnus, J.-M. et al. Planted forests and biodiversity. J. For. 104, 65–77 (2006).
Google Scholar 

67.
Gadow, K. V. & Hui, G. Y. Characterizing forest spatial structure and diversity Institute of Forest Management, Georg-August-University Göttingen, Büsgenweg 5, D-37077 Göttingen, Germany Published in: Sustainable Forestry in Temperate Regions; Proc. of an international workshop organized at the University of Lund, Sweden: 20–30. More

1.
Clutton-Brock, T. H. The Evolution of Parental Care (Princeton University Press, Princeton, 1991).
Google Scholar 
2.
Trivers, R. L. Parental investment and sexual selection. In Sexual Selection and the Descent of Man (ed. Campbell, B.) 136–179 (John Murray, Aldine, 1972).
Google Scholar 

3.
Alonzo, S. H. & Klug, H. Maternity, paternity and parental care. In The Evolution of Parental Care (eds Royle, N. J. et al.) 189–203 (Oxford University Press, Oxford, 2012).
Google Scholar 

4.
Møller, A. P. & Cuervo, J. J. The evolution of paternity and paternal care in birds. Behav. Ecol. 11, 472–485 (2000).
Article  Google Scholar 

5.
Neff, B. D. Paternity and condition affect cannibalistic behavior in nest-tending bluegill sunfish. Behav. Ecol. Sociobiol. 54, 377–384 (2003).
Article  Google Scholar 

6.
Benowitz, K. M., Head, M. L., Williams, C. A., Moore, A. J. & , Royle, N.J. ,. Male age mediates reproductive investment and response to paternity assurance. Proc. R. Soc. B 280, 20131124 (2013).
PubMed  Article  Google Scholar 

7.
Wisenden, B. D. Alloparental care in fishes. Rev. Fish Biol. Fish. 9, 45–70 (1999).
Article  Google Scholar 

8.
Griffin, A. S., Alonzo, S. H. & Cornwallis, C. K. Why do cuckolded males provide paternal care? PLoS ONE 11, e1001520 (2013).
CAS  Article  Google Scholar 

9.
Stevens, M. Bird brood parasitism. Curr. Biol. 23, R909–R913 (2013).
CAS  PubMed  Article  Google Scholar 

10.
Cohen, M. S., Hawkins, M. B., Stock, D. W. & Cruz, A. Early life-history features associated with brood parasitism in the cuckoo catfish, Synodontis multipunctatus (Siluriformes: Mochokidae). Philos. Trans. R. Soc. B 374, 20180205 (2019).
Article  Google Scholar 

11.
Taborsky, M. Sneakers, satellites, and helpers: Parasitic and cooperative behavior in fish reproduction. Adv. Stud. Behav. 23, 1–100 (1994).
Article  Google Scholar 

12.
Zahavi, A. Mate selection: A selection for handicap. J. Theor. Biol. 53, 205–214 (1975).
CAS  PubMed  Article  Google Scholar 

13.
Price, T., Schluter, D. & Heckman, N. E. Sexual selection when the female directly benefits. Biol. J. Linn. Soc. 48, 187–211 (1993).
Article  Google Scholar 

14.
Arnold, S. J. & Duvall, D. Animal mating systems: A synthesis based on selection theory. Am. Nat. 143, 317–348 (1994).
Article  Google Scholar 

15.
Klug, H., Alonzo, S. H. & Bonsall, M. B. Theoretical foundations of parental care. In The Evolution of Parental Care (eds Royle, N. J. et al.) 21–39 (Oxford University Press, Oxford, 2012).
Google Scholar 

16.
Nazareth, T. M. & Machado, G. Mating system and exclusive postzygotic paternal care in a Neotropical harvestman (Arachnida: Opiliones). Anim. Behav. 79, 547–554 (2010).
Article  Google Scholar 

17.
Matsumoto, Y., Tawa, A. & Takegaki, T. Female mate choice in a paternal brooding blenny: the process and benefits of mating with males tending young eggs. Ethology 117, 227–235 (2011).
Article  Google Scholar 

18.
Rohwer, S. Selection for adoption versus infanticide by replacement “mates” in birds. In Current Ornithology (ed. Johnston, R. F.) 353–395 (Plenum Press, New York, 1986).
Google Scholar 

19.
Valencia-Aguilar, A., Zamudio, K. R., Haddad, C. F. B., Bogdanowicz, S. M. & Prado, C. P. A. Show me you care: Female mate choice based on egg attendance rather than male or territorial traits. Behav. Ecol. 31, 1054–1064 (2020).
Article  Google Scholar 

20.
Schulte, L. M., Ringler, E., Rojas, B. & Stynoski, J. L. Developments in amphibian parental care research: History, present advances, and future perspectives. Herpetol. Monogr. 34, 71–97 (2020).
Article  Google Scholar 

21.
Vági, B., Végvári, Z., Liker, A., Freckleton, R. P. & Székely, T. Parental care and the evolution of terrestriality in frogs. Proc. R. Soc. B 286, 20182737 (2019).
PubMed  Article  Google Scholar 

22.
Guayasamin, J. M., Cisneros-Heredia, D. F., McDiarmid, R. W., Peña, P. & Hutter, C. R. Glassfrogs of ecuador: Diversity, evolution, and conservation. Diversity 12, 222 (2020).
CAS  Article  Google Scholar 

23.
Stynoski, J. L. Discrimination of offspring by indirect recognition in an egg-feeding dendrobatid frog, Oophaga pumilio. Anim. Behav. 78, 1351–1356 (2009).
Article  Google Scholar 

24.
Ringler, E., Beck, K. B., Weinlein, S., Huber, L. & Ringler, M. Adopt, ignore, or kill? Male poison frogs adjust parental decisions according to their territorial status. Sci. Rep. 7, 43544 (2017).
ADS  PubMed  PubMed Central  Article  Google Scholar 

25.
Waldman, B. Mechanisms of kin recognition. J. Theor. Biol. 128, 159–185 (1987).
Article  Google Scholar 

26.
Penn, D. & Frommen, J. Kin recognition: An overview of conceptual issues, mechanisms and evolutionary theory. In Animal Behaviour: Evolution and Mechanisms (ed. Kappeler, P.) 55–86 (Springer, Heidelberg, 2010).
Google Scholar 

27.
Delia, J. R., Bravo-Valencia, L. & Warkentin, K. The evolution of extended parental care in glassfrogs: Do egg-clutch phenotypes mediate coevolution between the sexes? Ecol. Monogr. 90, e01411 (2020).
Article  Google Scholar 

28.
Pašukonis, A. et al. Induced parental care in a poison frog: A tadpole cross-fostering experiment. J. Exp. Biol. 220, 3949–3954 (2017).
PubMed  PubMed Central  Article  Google Scholar 

29.
Townsend, D. & Moger, W. H. Plasma androgen levels during male parental care in a tropical frog (Eleutherodactylus). Horm. Behav. 21, 93–99 (1987).
CAS  PubMed  Article  Google Scholar 

30.
Knapp, R., Wingfield, J. C. & Bass, A. H. Steroid hormones and paternal care in the plainfin midshipman fish (Porichthys notatus). Horm. Behav. 35, 81–89 (1999).
CAS  PubMed  Article  Google Scholar 

31.
Pikus, A. E., Guindre-Parker, S. & Rubenstein, D. R. Testosterone, social status and parental care in a cooperatively breeding bird. Horm. Behav. 97, 85–93 (2018).
CAS  PubMed  Article  Google Scholar 

32.
Fischer, E. K. & O’Connell, L. A. Hormonal and neural correlates of care in active versus observing poison frog parents. BioRxiv 27, 765503 (2019).
Google Scholar 

33.
Goymann, W. & Dávila, P. F. Acute peaks of testosterone suppress paternal care: evidence from individual hormonal reaction norms. Proc. R. Soc. B 284, 20170632 (2017).
PubMed  Article  CAS  Google Scholar 

34.
Butin, J. D. Parental behavior and hormones in non-mammalian vertebrates. In Encyclopedia of Animal Behavior (eds Breed, M. & Moore, J.) 664–671 (Elsevier, Amsterdam, 2010).
Google Scholar 

35.
Townsend, D. S., Palmer, B. & Guillette, L. G. The lack of influence of exogenous testosterone on male parental behavior in a neotropical frog (Eleutherodactylus): A field experiment. Horm. Behav. 25, 313–322 (1991).
CAS  PubMed  Article  Google Scholar 

36.
Magee, S. E., Neff, B. D. & Knapp, R. Plasma levels of androgens and cortisol in relation to breeding behavior in parental male bluegill sunfish, Lepomis macrochirus. Horm. Behav. 49, 598–609 (2006).
CAS  PubMed  Article  Google Scholar 

37.
Ouyang, J. Q., Sharp, P. J., Dawson, A., Quetting, M. & Hau, M. Hormone levels predict individual differences in reproductive success in a passerine bird. Proc. R. Soc. B 278, 2537–2545 (2011).
CAS  PubMed  Article  Google Scholar 

38.
Mota, M. T. S., Franci, C. R. & Sousa, M. B. C. Hormonal changes related to paternal and alloparental care in common marmosets (Callithrix jacchus). Horm. Behav. 49, 293–302 (2006).
CAS  Article  Google Scholar 

39.
Romero, L. M. Seasonal changes in plasma glucocorticoid concentrations in free-living vertebrates. Gen. Comp. Endocrinol. 128, 1–24 (2002).
CAS  Article  Google Scholar 

40.
Consolmagno, R. C., Requena, G. S., Machado, G. & Brasileiro, C. A. Costs and benefits of temporary egg desertion in a rocky shore frog with male-only care. Behav. Ecol. Sociobiol. 70, 785–795 (2016).
Article  Google Scholar 

41.
Kelly, N. B. & Alonzo, S. H. Will male advertisement be a reliable indicator of paternal care, if offspring survival depends on male care? Proc. R. Soc. B 276, 3175–3183 (2009).
PubMed  Article  Google Scholar 

42.
Stiver, K. A. & Alonzo, S. H. Alloparental care increases mating success. Behav. Ecol. 22, 206–211 (2011).
Article  Google Scholar 

43.
Roldán, M. & Soler, M. Parental-care parasitism: How do unrelated offspring attain acceptance by foster parents? Behav. Ecol. 22, 679–691 (2011).
Article  Google Scholar 

44.
Maynard-Smith, J. Parental investment: A prospective analysis. Anim. Behav. 25, 1–9 (1977).
Article  Google Scholar 

45.
Valencia-Aguilar, A., Rodrigues, D. & Prado, C. P. A. Male care status influences the risk-taking decisions in a glassfrog. Behav. Ecol. Sociobiol. 74, 1–11 (2020).
Article  Google Scholar 

46.
Delia, J., Bravo-Valencia, L. & Warkentin, K. M. Patterns of parental care in Neotropical glassfrogs: Fieldwork alters hypotheses of sex-role evolution. J. Evol. Biol. 30, 898–914 (2017).
CAS  PubMed  Article  Google Scholar 

47.
Noronha, J. C. & Rodrigues, D. J. Reproductive behaviour of the glass frog Hyalinobatrachium cappellei (Anura: Centrolenidae) in the Southern Amazon. J. Nat. Hist. 52, 207–224 (2018).
Article  Google Scholar 

48.
Drake, D. L. & Ranvestel, A. W. Hyalinobatrachium colymbihpyllum (glass frog). Egg mass defense. Herpetol. Rev. 36, 434 (2005).
Google Scholar 

49.
Vockenhuber, E. A., Hödl, W. & Amézquita, A. Glassy fathers do matter: Egg attendance enhances embryonic survivorship in the glass frog Hyalinobatrachium valerioi. J. Herpetol. 43, 340–344 (2009).
Article  Google Scholar 

50.
Salgado, A. L. & Guayasamin, J. M. Parental care and reproductive behavior of the minute dappled glassfrog (Centrolenidae: Centrolene peristictum). S. Am. J. Herpetol. 13, 211–219 (2018).
Article  Google Scholar 

51.
Foster, W. A. & Treherne, J. E. Evidence for the dilution effect in the selfish herd from fish predation on a marine insect. Nature 293, 466–467 (1981).
ADS  Article  Google Scholar 

52.
Lehtonen, J. & Jaatinen, K. Safety in numbers: The dilution effect and other drivers of group life in the face of danger. Behav. Ecol. Sociobiol. 70, 449–458 (2016).
Article  Google Scholar 

53.
Gloag, R., Fiorini, V. D., Reboreda, J. C. & Kacelnik, A. Brood parasite eggs enhance egg survivorship in a multiply parasitized host. Proc. Biol. Sci. 279, 1831–1839 (2012).
PubMed  Google Scholar 

54.
Schulte, L. M. et al. The smell of success: Choice of larval rearing sites by means of chemical cues in a Peruvian poison frog. Anim. Behav. 81, 1147–1154 (2011).
Article  Google Scholar 

55.
Kam, Y. C. & Yang, H. W. Female–offspring communication in a Taiwanese tree frog, Chirixalus eiffingeri (Anura: Rhacophoridae). Anim. Behav. 64, 881–886 (2002).
Article  Google Scholar 

56.
Riedman, M. The evolution of alloparental care and adoption in mammals and birds. Q. Rev. Biol. 57, 405–435 (1982).
Article  Google Scholar 

57.
Briga, M., Pen, I. & Wright, J. Care for kin: Within-group relatedness and allomaternal care are positively correlated and conserved throughout the mammalian phylogeny. Biol. Lett. 8, 533–536 (2012).
PubMed  PubMed Central  Article  Google Scholar 

58.
Phillips, E., DeAngelis, R., Gogola, J. V. & Rhodes, J. S. Spontaneous alloparental care of unrelated offspring by non-breeding Amphiprion ocellaris in absence of the biological parents. Sci. Rep. 10, 4610 (2020).
ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

59.
Lee, H. J., Heim, V. & Meyer, A. Genetic evidence for prevalence of alloparental care in a socially monogamous biparental cichlid fish, Perissodus microlepis, from Lake Tanganyika supports the “selfish shepherd effect” hypothesis. Ecol. Evol. 6, 2843–2853 (2016).
PubMed  PubMed Central  Article  Google Scholar 

60.
Gosner, K. L. A simplified table for staging anuran embryos an larvae with notes on identification. Herpetologica 16, 183–190 (1960).
Google Scholar  More

2017 honey sampling
Beekeepers were invited to provide honey for analysis via a nationwide campaign publicised on the gardening programme, BBC Gardener’s World (broadcast July 2017). Participating beekeepers were asked to supply ~30 ml of honey from any date in 2017, reporting the date of sample collection and the location of the apiary, using a grid reference or postcode. In total 441 honey samples were processed from beekeepers.
Honey DNA extraction
Any wax was removed using sterile forceps and DNA was extracted from 10 g of honey using a modified version of the DNeasy Plant Mini extraction kit (Qiagen). Firstly, the 10 g of honey was made up to 30 ml with molecular grade water and incubated in a water bath at 65 °C for 30 min. Samples were then centrifuged (Sorvall RC-5B) for 30 min at 15,000 rpm, the supernatant was discarded, and the pellet resuspended in 400 μL of a buffer made from a mix of 400 μL AP1 from the DNeasy Plant Mini Kit (Qiagen), 80 μL proteinase K (1 mg/ml) (Sigma) and 1 μL RNase A (Qiagen). This was incubated again for 60 min at 65 °C in a water bath and then disrupted using a TissueLyser II (Qiagen) for 4 min at 30 Hz with 3 mm tungsten carbide beads. The remaining steps were carried out according to the manufacturer’s protocol, excluding the use of the QIAshredder and the second wash stage. The extracted DNA was purified using the OneStep PCR Inhibitor Removal Kit (Zymo Research) and diluted 1 in 10.
PCR and library preparation
Illumina MiSeq paired-end indexed amplicon libraries were created via a two-step PCR protocol. Two libraries were prepared for the DNA barcode regions, rbcL and ITS2. Initial amplification used the template specific primers rbcLaf and rbcLr50637, and ITS2F and ITS3R, with universal tails designed to attach custom indices in the second-round PCR. To improve clustering on the Illumina MiSeq, a 6N sequence was also added between the forward template specific primer and the universal tail.
Forward universal tail, 6N sequence and rbcLaf: [ACACTCTTTCCCTACACGACGCTCTTCCGATCT]NNNNNN[ATGTCACCACAAACAGAGACTAAAGC]
Reverse universal tail and rbcLr506: [GTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT][AGGGGACGACCATACTTGTTCA]
Forward universal tail, 6N sequence and ITS2F: [ACACTCTTTCCCTACACGACGCTCTTCCGATCT]NNNNNN[ATGCGATACTTGGTGTGAAT]
Reverse universal tail and ITS3R: [GTGACTGGAGTTCAGACGTGTGCTCTTCCGATCT][GACGCTTCTCCAGACTACAAT]
This first PCR used a final volume of 20 μl: 2 μl template DNA, 10 μl of 2× Phusion Hot Start II High-Fidelity Mastermix (New England Biolabs UK), 0.4 μl (2.5 µM) forward and reverse primers, and 7.2 μl of PCR grade water. Thermal cycling conditions for rbcL were: 98 °C for 30 s, 95 °C for 2 min; 95 °C for 30 s, 50 °C for 30 s, 72 °C for 40 s (40 cycles); 72 °C for 5 min, 30 °C for 10 s. Thermal cycling conditions for the first ITS2 PCR were: 98 °C for 30 s 94 °C for 5 min; 94 °C for 30 s, 56 °C for 30 s, 72 °C for 40 s (40 cycles); 72 °C for 10 min, 30 °C for 1 min. The initial PCR was carried out three times and pooled.
The pooled products from the first PCR were purified following Illumina’s 16S Metagenomic Sequencing Library Preparation protocol using Agencourt AMPure XP beads (Beckman Coulter). The purified PCR product from round one was followed by a second round of amplification to anneal custom unique and identical i5 and i7 indices to each sample (Ultramer, Integrated DNA Technologies).
This index PCR stage used a final volume of 25 μl reaction (12.5 μl of 2× Phusion Hot Start II High-Fidelity Mastermix, 1 μl of i7 Index Primer and i5 Index Primer, 6.5 μl of PCR grade water, and 5 μl of purified first-round PCR product). Thermal cycling conditions were: 98 °C for 30 s; 95 °C for 30 s, 55 °C for 30 s, 72 °C for 30 s (8 cycles); 72 °C for 5 min, 4 °C for 10 min. Following the index PCR, a 1% gel was run to verify its success. The index PCR product was then purified following the PCR clean-up two sections of the Illumina protocol. The purified products of the index PCR were quantified using a Qubit 3.0 fluorescence spectrophotometer (Thermo Fisher Scientific) and pooled at equal concentrations to produce the final library. Positive and negative controls were amplified and sequenced alongside honey samples. The positive control was made from a mixture of five tropical tree species that were not present in the survey site. The species Baccaurea stipulata, Colona serratifolia., Dillenia excelsa, Kleinhovia hospita, and Pterospermum macrocarpum were used, taking 5 μl from each separate DNA extraction and mixing, before following the protocol as with the honey samples. All five species were detected within the sequencing results.
Bioinformatic analysis
Sequence data were processed using a modified data analysis pipeline14,38. Raw reads were trimmed to remove low-quality regions (Trimmomatic v. 0.33), paired, and then merged (FLASH v. 1.2.11), with merged reads shorter than 450 bp discarded. Identical reads were dereplicated within samples and then clustered at 100% identity across all samples (vsearch v. 2.3.2), with singletons (sequence reads that occurred only once across all samples) discarded.
The Barcode Wales and Barcode UK projects provide 98% coverage for the native flowering plants and conifers of the UK37. This reference library was supplemented with a curated library of the non-native and horticultural species, downloaded from GenBank. This UK species list was generated using the list of native species of the UK from Stace (2010)39, 505 naturalised alien species (BSBI), and horticultural species from the IRIS BG database at the National Botanic Garden of Wales.
The sequence data from the honey samples were compared against the reference database using blastn, using the script vsearch-pipe.py. The top BLAST hits were then summarised using the script vsearch_blast_summary.py. Sequences with bit scores below the 1st percentile were excluded. If the top bit scores of a sequence matched to a single species, then the sequence was identified to that species. If the top bit scores matched to different species within the same genus, then the result was attributed to the genus level. If the top bit score belonged to multiple genera within the same family then a family level designation was made. Sequences that returned families from different clades were excluded. These automated identifications were then checked manually for botanical veracity. To check identified plant species against their availability across the UK, species records from the BSBI (Botanical Society of Britain and Ireland) were used for native species, while commercial availability for horticultural species was verified with the RHS Plant Finder40. Within each sample, the number of sequences returned from rbcL and ITS2 for each plant taxon was summed to combine the results of each marker.
The proportion of sequences was used in the analysis, which has been shown to be an appropriate method to control for differences in read number41. Alternatively, the sequencing data can be rarefied, but this has been criticised as a statistical technique, due to requiring the removal of valid data41. To investigate the impact of rarefying on the conclusions drawn from the data, all analyses were rerun with rarefied data (Supplementary Results).
1952 Honey sampling
In 1952, 855 honey samples were characterised from 66 counties across the UK and Ireland using melissopalynology15,16. The methods reported for the research conducted in 1952 are described here fully for comparison. Samples were obtained via a general appeal and were all collected during the honey season of 1952. For each honey sample, ~200 pollen grains were identified using the morphology of the pollen under the microscope, following a standardised protocol42. To extract the pollen, 10 g of honey was dissolved in 20 ml of distilled water, from which 10 ml was taken and centrifuged at ~2000 rpm for one minute. The supernatant was discarded, and the sediment retained, and then the process was repeated for the remaining liquid. From the sediment, a drop was transferred to a glass slide and spread out over an area of 1 cm2, before being stained with fuchsin and dried. Euparal vert was used as a final mounting medium. Pollen was identified by comparison with a reference library of pollen preparations and available pollen morphological data43,44. Each plant taxon found in the sampled honey was reported according to the proportion of pollen grains found and classed into predominant ( >45% of pollen grains), secondary (15–45% of pollen grains) and important minor (1–15% of pollen grains). The location data for the honey samples were restricted to the county level, and summary data tables were presented for each UK county that returned honey.
Comparing the 1952 and 2017 honey samples
The plants detected using DNA metabarcoding and melissopalynology have been compared in previous studies with concordance found between the two methods45,46,47,48. Both methods detect the same major taxa, but rarer species in a sample are less likely to be found consistently, both when comparing methods and also during replicates of the same method45,46,47. DNA metabarcoding is often able to detect more taxa when compared to melissopalynology, by identifying rarer species in the sample and by achieving higher taxonomic resolution in certain cases. While melissopalynology uses counts of pollen grains to provide a starting point for quantitative analysis, DNA metabarcoding as a process is semi-quantitative, with biases associated with the process of DNA extraction, PCR and sequencing33,45. To allow for these considerations we placed the proportion of DNA sequence reads and pollen counts into four broad abundance classes matching the classifications used in melissopalynology (predominant, secondary, important minor and minor) and focus our analyses and conclusions on changes in the frequency of occurrence of the major taxa, classed as predominant and secondary. Both methods capture information on both nectar and pollen plants within the honey, however, certain species can be over or under represented in pollen analysis compared to their relative nectar contribution49. Both pollen and nectar plants are required to meet the foraging requirements of pollinators.
Statistics and reproducibility
Statistical analysis of DNA metabarcoding data
To understand how the plant taxa composition within the honey sample was structured in space and time, the effect of time (measured as the calendar month number in 2017), latitude and longitude of sampling location were included in a single, two-tailed generalized linear model using the ‘manyglm’ function in the package ‘mvabund’50. Honey samples with missing metadata were excluded, giving a sample size of 428. An abundance table of taxa (number of sequence reads) found in each sample was set as the multivariate response variable and a common set of predictor variables (month, latitude and longitude) were fit using a negative binomial distribution. The number of sequence reads per sample was included as an “offset” in the model in order to control for differences in the number of sequence reads between samples. Monte Carlo resampling was used to test for significant community-level responses to our predictors. The strong mean-variance relationship in the data (Supplementary Fig. 6) and the distribution of the count data (Supplementary Figs. 7, 8) support the use of a negative binomial distribution in the model. The appropriateness of the models was checked by visual inspection of the residuals against predicted values from the models (Supplementary Figs. 9–11).
We completed a spatial eigenfunction analysis using distance-based Moran’s eigenvectors. Moran’s Eigenvector Maps were computed using the ‘mem’ function from the adespatial package. Moran’s I was computed for each taxa using the ‘moran.randtest’, with Bonferroni correction for multiple testing. The direction of autocorrelation (positive and negative) was tested using the ‘moranNP.randtest’ function, using the adespatial package in R.
Statistical analysis of the 1952 and 2017 honey samples
Abundance classes were assigned based on the percentage of reads returned for the two DNA regions rbcL and ITS2, matching the classifications used in melissopalynology. Plant taxa represented by over 45% of reads were designated predominant for that sample; between 15 and 45% were secondary; between 1 and 15% were important minor taxa, and More

1.
Roddy, K. M. & Arita-Tsutsumi, L. A history of honey bees in the Hawaiian islands. J. Hawaiian Pac. Agric. 8, 59–70 (1997).
Google Scholar 
2.
Danka, R. G., Hellmich, R. L., Rinderer, T. E. & Collins, A. M. Diet-selection ecology of tropically and temperately adapted honey-bees. Anim. Behav. 35, 1858–1863 (1987).
Article  Google Scholar 

3.
Roberts, J. M. K., Anderson, D. L. & Durr, P. A. Absence of deformed wing virus and Varroa destructor in Australia provides unique perspectives on honeybee viral landscapes and colony losses. Sci. Rep. https://doi.org/10.1038/s41598-017-07290-w (2017).
PubMed  PubMed Central  Article  Google Scholar 

4.
de Guzman, L. I., Rinderer, T. E. & Stelzer, J. A. DNA evidence of the origin of Varroa jacobsoni Oudemans in the Americas. Biochem. Genet. 35, 327–335. https://doi.org/10.1023/a:1021821821728 (1997).
PubMed  Article  Google Scholar 

5.
Ramsey, S. D. et al. Varroa destructor feeds primarily on honey bee fat body tissue and not hemolymph. Proc. Natl. Acad. Sci. USA. 116, 1792–1801. https://doi.org/10.1073/pnas.1818371116 (2019).
CAS  PubMed  Article  Google Scholar 

6.
Rosenkranz, P., Aumeier, P. & Ziegelmann, B. Biology and control of Varroa destructor. J. Invertebr. Pathol. 103, S96–S119 (2010).
PubMed  Article  Google Scholar 

7.
Sammataro, D., Gerson, U. & Needham, G. Parasitic mites of honey bees: Life history, implications, and impact. Annu. Rev. Entomol. 45, 519–548 (2000).
CAS  PubMed  Article  Google Scholar 

8.
Amdam, G. V., Hartfelder, K., Norberg, K., Hagen, A. & Omholt, S. W. Altered physiology in worker honey bees (Hymenoptera: Apidae) infested with the mite Varroa destructor (Acari: Varroidae): A factor in colony loss during overwintering?. J. Econ. Entomol. 97, 741–747 (2004).
PubMed  Article  Google Scholar 

9.
Dejong, D., Dejong, P. H. & Goncalves, L. S. Weight-loss and other damage to developing worker honeybees from infestation with Varroa jacobsoni. J. Apic. Res. 21, 165–167. https://doi.org/10.1080/00218839.1982.11100535 (1982).
Article  Google Scholar 

10.
Ramadan, M. M., Reimer, N. J., Oishi, D. E., Young, C. L. & Heu, R. A. Varroa Mite Varroa destructor Anderson and Trueman (Acari: Varroidae) (Springer, New York, 2019).
Google Scholar 

11.
Martin, S. J. et al. Global honey bee viral landscape altered by a parasitic mite. Science 336, 1304–1306. https://doi.org/10.1126/science.1220941 (2012).
ADS  CAS  PubMed  Article  Google Scholar 

12.
Seeley, T. D. Honey bees of the Arnot Forest: A population of feral colonies persisting with Varroa destructor in the northeastern United States. Apidologie 38, 19–29 (2007).
Article  Google Scholar 

13.
Brettell, L. E. & Martin, S. J. Oldest Varroa tolerant honey bee population provides insight into the origins of the global decline of honey bees. Sci. Rep. https://doi.org/10.1038/srep45953 (2017).
PubMed  PubMed Central  Article  Google Scholar 

14.
Nielsen, D. I. Genetic structure of feral honey bee (Apis mellifera L.) populations in California Ph.D. thesis, University of California, Davis (2000).

15.
Doebler, S. A. The rise and fall of the honeybee: Mite infestations challenge the bee and the beekeeping industry. Bioscience 50, 738–742. https://doi.org/10.1641/0006-3568(2000)050[0738:Trafot]2.0.Co;2 (2000).
Article  Google Scholar 

16.
Fuchs, S. Preference for drone brood cells by Varroa jacobsoni oud in colonies of Apis-Mellifera-Carnica. Apidologie 21, 193–199 (1990).
Article  Google Scholar 

17.
Boot, W. J., Calis, J. N. M. & Beetsma, J. Differential periods of varroa mite invasion into worker and drone cells of honey-bees. Exp. Appl. Acarol. 16, 295–301 (1992).
Article  Google Scholar 

18.
Estoup, A., Solignac, M. & Cornuet, J. Precise assessment of the number of patrilines and of genetic relatedness in honeybee colonies. Proc. R. Soc. Lond. B 258, 1–7 (1994).
ADS  CAS  Article  Google Scholar 

19.
Tarpy, D. R., Nielsen, R. & Nielsen, D. I. A scientific note on the revised estimates of effective paternity frequency in Apis. Insectes Soc. 51, 203–204 (2004).
Article  Google Scholar 

20.
Akyol, E., Yeninar, H. & Kaftanoglu, O. Live weight of queen honey bees (Apis mellifera L.) predicts reproductive characteristics. J. Kansas Entomol. Soc. 81, 92–100. https://doi.org/10.2317/jkes-705.13.1 (2008).
Article  Google Scholar 

21.
Amiri, E., Strand, M. K., Rueppell, O. & Tarpy, D. R. Queen quality and the impact of honey bee diseases on queen health: Potential for interactions between two major threats to colony health. Insects. https://doi.org/10.3390/insects8020048 (2017).
PubMed  PubMed Central  Article  Google Scholar 

22.
Tarpy, D. R., Keller, J. J., Caren, J. R. & Delaney, D. A. Experimentally induced variation in the physical reproductive potential and mating success in honey bee queens. Insectes Soc. 58, 569–574. https://doi.org/10.1007/s00040-011-0180-z (2011).
Article  Google Scholar 

23.
Hatjina, F. et al. A review of methods used in some European countries for assessing the quality of honey bee queens through their physical characters and the performance of their colonies. J. Apic. Res. 53, 337–363. https://doi.org/10.3896/ibra.1.53.3.02 (2014).
Article  Google Scholar 

24.
De Souza, D. A. et al. Morphometric identification of queens, workers and intermediates in in vitro reared honey bees (Apis mellifera). PLoS ONE. https://doi.org/10.1371/journal.pone.0123663 (2015).
PubMed  PubMed Central  Article  Google Scholar 

25.
Woyke, J. Correlations between the age at which honeybee brood was grafted, characteristics of the resultant queens, and results of insemination. J. Apic. Res. 10, 45–55 (1971).
Article  Google Scholar 

26.
Dedej, S., Hartfelder, K., Aumeier, P., Rosenkranz, P. & Engels, W. Caste determination is a sequential process: Effect of larval age at grafting on ovariole number, hind leg size and cephalic volatiles in the honey bee (Apis mellifera carnica). J. Apic. Res. 37, 183–190 (1998).
Article  Google Scholar 

27.
Al-Lawati, H., Kamp, G. & Bienefeld, K. Characteristics of the spermathecal contents of old and young honeybee queens. J. Insect Physiol. 55, 116–121 (2009).
CAS  PubMed  Article  Google Scholar 

28.
Tarpy, D. R., Keller, J. J., Caren, J. R. & Delaney, D. A. Assessing the mating “health” of commercial honey bee queens. J. Econ. Entomol. 105, 20–25 (2012).
PubMed  Article  Google Scholar 

29.
Pettis, J. S., Rice, N., Joselow, K., vanEngelsdorp, D. & Chaimanee, V. Colony failure linked to low sperm viability in honey bee (Apis mellifera) queens and an exploration of potential causative factors. PLoS ONE. https://doi.org/10.1371/journal.pone.0147220 (2016).
PubMed  PubMed Central  Article  Google Scholar 

30.
Woyke, J. Natural and artificial insemination of queen honeybees. Bee World 43, 21–25 (1962).
Article  Google Scholar 

31.
Delaney, D. A., Keller, J. J., Caren, J. R. & Tarpy, D. R. The physical, insemination, and reproductive quality of honey bee queens (Apis mellifera). Apidologie 42, 1–13. https://doi.org/10.1051/apido/2010027 (2011).
Article  Google Scholar 

32.
McAfee, A. et al. Vulnerability of honey bee queens to heat-induced loss of fertility. Nat. Sustain. https://doi.org/10.1038/s41893-020-0493-x (2020).
Article  Google Scholar 

33.
Chaimanee, V., Evans, J. D., Chen, Y., Jackson, C. & Pettis, J. S. Sperm viability and gene expression in honey bee queens (Apis mellifera) following exposure to the neonicotinoid insecticide imidacloprid and the organophosphate acaricide coumaphos. J. Insect Physiol. 89, 1–8 (2016).
CAS  PubMed  Article  Google Scholar 

34.
Williams, G. R. et al. Neonicotinoid pesticides severely affect honey bee queens. Sci. Rep. https://doi.org/10.1038/srep14621 (2015).
PubMed  PubMed Central  Article  Google Scholar 

35.
Rangel, J. & Tarpy, D. R. The combined effects of miticides on the mating health of honey bee (Apis mellifera L.) queens. J. Apicult. Res. 54, 275–283. https://doi.org/10.1080/00218839.2016.1147218 (2015).
Article  Google Scholar 

36.
Buechler, R. et al. Standard methods for rearing and selection of Apis mellifera queens. J. Apicult. Res. https://doi.org/10.3896/ibra.1.52.1.07 (2013).
Article  Google Scholar 

37.
Lee, K. V., Goblirsch, M., McDermott, E., Tarpy, D. R. & Spivak, M. Is the brood pattern within a honey bee colony a reliable indicator of queen quality?. Insects. https://doi.org/10.3390/insects10010012 (2019).
PubMed  PubMed Central  Article  Google Scholar 

38.
de Miranda, J. R. et al. Standard methods for virus research in Apis mellifera. J. Apicult. Res. https://doi.org/10.3896/ibra.1.52.4.22 (2013).
Article  Google Scholar 

39.
Kevill, J. L. et al. The pathogen profiles of queen honey bees does not reflect those of their colonies workers. Insects 11, 382. https://doi.org/10.3390/insects11060382 (2020).
PubMed Central  Article  Google Scholar 

40.
Evans, J. D. et al. Standard methods for molecular research in Apis mellifera. J. Apicult. Res. https://doi.org/10.3896/ibra.1.52.4.11 (2013).
Article  Google Scholar 

41.
Jones, O. R. & Wang, J. COLONY: A program for parentage and sibship inference from multilocus genotype data. Mol. Ecol. Resour. 10, 551–555. https://doi.org/10.1111/j.1755-0998.2009.02787.x (2010).
PubMed  Article  Google Scholar 

42.
Nielsen, R., Tarpy, D. R. & Reeve, H. K. Estimating effective paternity number in social insects and the effective number of alleles in a population. Mol. Ecol. 12, 3157–3164 (2003).
PubMed  Article  Google Scholar 

43.
Neumann, P. & Carreck, N. L. Honey bee colony losses. J. Apic. Res. 49, 1–6 (2010).
Article  Google Scholar 

44.
van Engelsdorp, D. & Meixner, M. D. A historical review of managed honey bee populations in Europe and the United States and the factors that may affect them. J. Invertebr. Pathol. 103, S80–S95 (2010).
Article  Google Scholar 

45.
Ellis, J. D., Evans, J. D. & Pettis, J. Colony losses, managed colony population decline, and Colony Collapse Disorder in the United States. J. Apic. Res. 49, 134–136. https://doi.org/10.3896/ibra.1.49.1.30 (2010).
Article  Google Scholar 

46.
Lee, K. V. et al. A national survey of managed honey bee 2013–2014 annual colony losses in the USA: Results from the Bee Informed Partnership. Apidologie. https://doi.org/10.1007/s13592-015-0356-z (2015).
Article  Google Scholar 

47.
Steinhauer, N. et al. Drivers of colony losses. Curr. Opin. Insect Sci. 26, 142–148. https://doi.org/10.1016/j.cois.2018.02.004 (2018).
PubMed  Article  Google Scholar 

48.
Tarpy, D. R., Delaney, D. A. & Seeley, T. D. Mating frequencies of honey bee queens (Apis mellifera L.) in a population of feral colonies in the northeastern United States. PLoS ONE. https://doi.org/10.1371/journal.pone.0118734 (2015).
PubMed  PubMed Central  Article  Google Scholar 

49.
Lensky, Y. & Demter, M. Mating flights of the queen honeybee (Apis mellifera) in a subtropical climate. Comp. Biochem. Physiol. 81, 229–241 (1985).
Article  Google Scholar 

50.
USDA-NASS. (ed National Agricultural Statistics Service) (2018).

51.
DeGrandi-Hoffman, G. et al. Comparisons of pollen substitute diets for honey bees: Consumption rates by colonies and effects on brood and adult populations. J. Apic. Res. 47, 265–270. https://doi.org/10.3896/ibra.1.47.4.06 (2008).
Article  Google Scholar 

52.
Loftus, J. C., Smith, M. L. & Seeley, T. D. How honey bee colonies survive in the wild: Testing the importance of small nests and frequent swarming. PLoS ONE. https://doi.org/10.1371/journal.pone.0150362 (2016).
PubMed  PubMed Central  Article  Google Scholar 

53.
Le Conte, Y. et al. Honey bee colonies that have survived Varroa destructor. Apidologie 38, 566–572 (2007).
Article  Google Scholar 

54.
Chen, Y. P., Evans, J. & Feldlaufer, M. Horizontal and vertical transmission of viruses in the honeybee, Apis mellifera. J. Invertebr. Pathol. 92, 152–159 (2006).
PubMed  Article  Google Scholar 

55.
de Miranda, J. R. & Fries, I. Venereal and vertical transmission of deformed wing virus in honeybees (Apis mellifera L.). J. Invertebr. Pathol. 98, 184–189 (2008).
PubMed  Article  Google Scholar 

56.
Yue, C., Schroder, M., Bienefeld, K. & Genersch, E. Detection of viral sequences in semen of honeybees (Apis mellifera): Evidence for vertical transmission of viruses through drones. J. Invertebr. Pathol. 92, 105–108 (2006).
CAS  PubMed  Article  Google Scholar 

57.
Amiri, E., Meixner, M. D. & Kryger, P. Deformed wing virus can be transmitted during natural mating in honey bees and infect the queens. Sci. Rep. https://doi.org/10.1038/srep33065 (2016).
PubMed  PubMed Central  Article  Google Scholar 

58.
Szalanski, A. L., Tripodi, A. D., Trammel, C. E. & Downey, D. Mitochondrial DNA genetic diversity of honey bees, Apis mellifera in Hawaii. Apidologie 47, 679–687. https://doi.org/10.1007/s13592-015-0416-4 (2016).
CAS  Article  Google Scholar 

59.
Danka, R. G., Harris, J. W., Villalobos, E. & Glenn, T. Varroa destructor resistance of honey bees in Hawaii, USA, with different genetic proportions of Varroa Sensitive Hygiene (VSH). J. Apic. Res. 51, 288–290. https://doi.org/10.3896/ibra.1.51.3.13 (2012).
Article  Google Scholar 

60.
Metz, B. N. & Tarpy, D. R. Reproductive senescence in drones of the honey bee (Apis mellifera). Insects. https://doi.org/10.3390/insects10010011 (2019).
PubMed  PubMed Central  Article  Google Scholar 

61.
Sturup, M., Baer-Imhoof, B., Nash, D. R., Boomsma, J. J. & Baer, B. When every sperm counts: Factors affecting male fertility in the honeybee Apis mellifera. Behav. Ecol. 24, 1192–1198. https://doi.org/10.1093/beheco/art049 (2013).
Article  Google Scholar  More

The keystone species concept is an important aspect of population ecology, community ecology, and conservation biology1,2, and its application is likely to be critical with ongoing climate change3. Keystone species can be identified because they have a larger effect on communities and ecosystems than would be predicted based on their abundance or dominance. Loss of keystone species within communities and ecosystems is likely to result in secondary extinction events, and in extreme cases these events can lead to community and ecosystem collapse4. The critical importance of keystone species is derived from the wide range of biotic interactions they engage in with other community members (predation, competition, herbivory, mutualism, facilitation, etc.) and their influence on abiotic environmental conditions2. Keystone species have been described in a range of ecosystems (e.g., marine, fresh water, terrestrial, etc.) and have included a variety of taxa (e.g., fungi, animals, and plants)1,3,5.
Plant communities consisting of isolated or scattered trees occur across the globe, and such trees have been described as keystone species, or “keystone structures”6. This certainly applies to trees and shrubs that are members of plant communities in arid and semi-arid habitat7. Many members of Acacia s.l. (Fabaceae: Mimosoideae8), which are broadly distributed around the world, are considered keystone species within the communities they reside. For example, they are considered keystone species in parts of Australia9, Pakistan10, the Kalahari Desert, Botswana11, Tunisia12,13,14, the Sinai Desert, Egypt15,16, and south-eastern Egypt16,17. As pointed out by Abdallah et al.12, isolated trees in arid habitats, including Vachellia species., have several characteristics that contribute to their keystone status: (1) shade from their canopies prevents extreme temperature fluctuations, increases soil moisture levels, and provides shelter for wildlife, (2) they improve soil conditions through biological nitrogen fixation and litter fall by increasing soil nitrogen content, organic carbon, and water-holding capacity, (3) they increase plant and animal biodiversity as a consequence of characteristics one and two, (4) they provide a source of food for wildlife, and (5) they provide a source of fuel, fodder, and medicines for local people and their domesticated animals. Because of their critical importance, a full characterization of keystone species and the roles they play within communities and ecosystems is urgently needed; especially as they are adversely impacted by various human activities.
The Gebel Elba mountain range is an extension of the Afromontane “biodiversity hotspot” and is at the northern limit of the Eritreo-Arabian province and the Sahel regional transition zone18. The relatively high abundance of moisture of this mountain range leads to higher plant biodiversity than reported elsewhere in Egypt, it consists of 458 species, which constitutes approximately 21% of the Egyptian flora19,20. According to the plant checklist provided by Boulos21, the flora of Egypt consists of 2100 taxa belonging to 755 genera and 129 families; including 45 genera and 228 taxa in the Fabaceae. Gebel Elba is one of the seven main phytogeographical regions in Egypt21. Additionally, the region’s tree and shrub species diversity is higher than in any other regions in Egypt19, with some Sahelian woody elements restricted to the Gebel Elba region and not reported elsewhere in Egypt. Of the 10 Vachellia (synonym: Acacia8) species reported in Egypt, seven are known to occur in the Gebel Elba region, with Vachellia asak (synonym: Acacia asak) and Vachellia oerfota subsp. oerfota (synonym: Acacia oerfota subsp. oerfota) restricted to this region.
An analysis of the plant communities of wadi Yahmib and three of its tributaries, on the north-western slopes of Gebel Elba, revealed the presence of seven plant communities, with these communities being arrayed across an elevational (environmental) gradient17. The Vachellia tortilis subsp. tortilis (synonym: Acacia tortilis subsp. tortilis) community was the main vegetation type on Gebel Elba. This community type occurred commonly in the water channels of wadis and gravel terraces from low to mid elevations (130–383 m), and the species was a member of all of the other six communities in the study area17. In addition, Vachellia tortilis subsp. raddiana (synonym: Acacia tortilis subsp. raddiana) was an overstory co-dominant species in another community on Gebel Elba. Finally, a third acacia species, Vachellia etbaica (synonym: Acacia etbaica), was also detected in this study.
Within arid and semi-arid ecosystems across north Africa and the Arabian Peninsula, plant ecologists have focused their attention on describing the vegetation of wadis that drain to the Red Sea, with these studies focusing on keystone Vachellia species12,13,14,15,16,17,22,23. The present study aimed to contribute to this body of knowledge by determining the distribution, abundance, and describing the growth characteristics of three Vachellia tree taxa in wadi Khoda and wadi Rahaba, in Gebel Elba National Park, south-eastern Egypt. These data will allow us to provide detailed descriptions of the characteristics of these three taxa. This study is essential at this moment because these tree taxa are keystone species within these ecosystems, and their presence and conservation are likely to be threatened by human activities and ongoing climate change. More

Virus and cells
RVFV strain 35/74 was originally isolated from the liver of a sheep that died during a RVFV outbreak in the Free State province of South Africa in 197421. The strain was previously passaged four times in suckling mouse brain and three times in BHK cells. The virus used for IV inoculation of sheep was prepared by a further amplification in BHK-21 cells (ATCC CCL-10) cultured in CO2-independent medium (CIM, Invitrogen), supplemented with 5% FBS (Bodinco) and 1% Pen/Strep (Invitrogen).
To prepare a virus-spiked blood meal for membrane feeding of mosquitoes, the virus was amplified in Aedes albopictus C6/36 cells (ATCC CRL-1660). To this end, C6/36 cells were inoculated with a multiplicity of infection of 0.005 and cultured at 28 °C in absence of CO2 in L-15 medium (Sigma) supplemented with 10% fetal bovine serum (FBS), 2% Tryptose Phosphate Broth (TPB) and 1% MEM nonessential amino acids solution (MEMneaa). At 4 days post infection, culture medium was harvested, cleared by slow-speed centrifugation and titrated using Vero-E6 cells (ATCC CRL-1586), grown in DMEM supplemented with GlutaMAX, 3% FBS, 1% Pen/Strep and 1% Fungizone (DMEM +) at 37 °C and 5% CO2. Titers were determined using the Spearman-Kärber algorithm22,23.
Mosquito rearing and feeding on lambs
Rockefeller strain Ae. aegypti mosquitoes (Bayer AG, Monheim, Germany) were maintained at Wageningen University, Wageningen, the Netherlands, as described24. Briefly, mosquitoes were kept in Bugdorm-1 rearing cages at a temperature of 27 °C with a 12:12 light:dark cycle and a relative humidity of 70% with a 6% glucose solution provided ad libitum. Mosquitoes were subsequently transported to biosafety level three (BSL-3) facilities of Wageningen Bioveterinary Research (Lelystad, the Netherlands), where the mosquitoes were maintained with sugar water (6% sucrose in H2O), provided via soaked cotton pads covered with a lid to prevent evaporation in an insect incubator (KBWF 240, Binder) at 28 °C at a humidity of 70% and a 16:8 light:dark cycle.
Mosquito feeding on lambs was preceded by sedating the lambs with IV administration of medetomidine (Sedator). When fully sedated, cardboard boxes containing 40–50 female mosquitoes were placed on the shaved inner thigh of each hind leg (Fig. 1b,c). After 20 min of feeding, cardboard boxes were removed and atipamezol (Atipam) was administered via intramuscular (IM) route to wake up the animals. Fully engorged mosquitoes were collected using an automated insect aspirator and maintained with sugar water (6% sucrose in H2O), provided via soaked cotton pads covered with a lid to prevent evaporation, in an insect incubator (KBWF 240, Binder) at 28 °C at a humidity of 70% and a 16:8 light:dark cycle.
Feeding of mosquitoes using a Hemotek system
Blood meals to be used for Hemotek membrane feeding were prepared essentially as described before25. Briefly, erythrocytes were harvested from freshly collected bovine EDTA blood by slow-speed centrifugation (650 xg), followed by three wash steps with PBS. Washed erythrocytes were resuspended in L15 complete medium (L15 + 10% FBS, 2% TPB, 1% MEMneaa) to a concentration that is four times higher than found in blood. To prepare a blood meal, one part of the erythrocyte suspension was mixed with two parts of culture medium containing RVFV resulting in a final titer of 107.5 TCID50/ml as determined on Vero-E6 cells.
Mosquitoes were allowed to take a RVFV-spiked blood meal through a Parafilm M membrane using the Hemotek PS5 feeding system (Discovery Workshops, Lancashire, United Kingdom). Feeding was performed in plastic buckets (1 l) covered with mosquito netting. After blood feeding for approximately 1.5–2 h, fully engorged mosquitoes were collected using an automated insect aspirator and maintained with sugar water (6% sucrose in H2O), provided via soaked cotton pads covered with a lid to prevent evaporation in an insect incubator (KBWF 240, Binder) at 28 °C at a humidity of 70% and a 16:8 light:dark cycle.
Virus isolation
Virus isolation from plasma samples was performed using BHK-21 cells, seeded at a density of 20,000 cells/well in 96-wells plates. Serial dilutions of samples were incubated with the cells for 1.5 h before medium replacement. Cytopathic effect was evaluated after 5–7 days post infection. Virus titers (TCID50/ml) were determined using the Spearman-Kärber algorithm22,23.
To check for positive saliva, mosquitoes were sedated on a semi-permeable CO2-pad connected to 100% CO2 and wings and legs were removed. Saliva was collected by forced salivation using 20 µl filter tips containing 7 µl of a 1:1 mixture of FBS and 50% sucrose (capillary tube method). After 1–1.5 h, saliva samples were collected and used to inoculate Vero-E6 cell monolayers. Cytopathic effect (CPE) was scored 5–7 days later.
Serology
Weekly collected serum samples were used to detect RVFV-specific antibodies using the ID Screen Rift Valley Fever Competition Multi-species ELISA (ID-VET). This ELISA measures percentage competition between antibodies present in test sera and a monoclonal antibody. Neutralizing antibodies were detected using the RVFV-4 s-based virus neutralization test as described26.
RT-qPCR
Viral RNA was isolated with the NucliSENS easyMAG system according the manufacturer’s instructions (bioMerieux, France) from 0.5 ml plasma samples. Briefly, 5 µl RNA was used in a RVFV RT-qPCR using the LightCycler one-tube RNA Amplification Kit HybProbe (Roche, Almere, The Netherlands) in combination with a LightCycler 480 real-time PCR system (Roche) and the RVS forward primers (AAAGGAACAATGGACTCTGGTCA), the RVAs (CACTTCTTACTACCATGTCCTCCAAT) reverse primer and a FAM-labelled probe RVP (AAAGCTTTGATATCTCTCAGTGCCCCAA). Primers and probes were earlier described by Drosten et al.27. Virus isolations were performed on RT-qPCR positive samples with a threshold above 105 RNA copies/ml as this was previously shown to be a cut-off point below which no live virus can be isolated.
Pathology and (immuno)histopathology
Liver samples were placed on ice during the necropsies and subsequently stored at − 80 °C until virus isolations and RT-qPCR Tissue samples for histology and IHC were collected, placed in 10% neutral buffered formalin, embedded into paraffin and prepared for H&E staining or IHC staining for RVFV antigen using the RVFV Gn-specific 4-D4 mAb as described5.
Statistics
For statistical analysis, mosquito feeding and mosquito saliva positive rates per group were compared by fitting logistic regression mixed models where lamb or membrane were introduced as random effects. To compare viremia (based on virus isolation results) the area under the curve (AUC) representing the overall viremia during the infected period was calculated for each infected sheep. This AUC and peak of viremia was used for comparison between groups, which was done by fitting linear regression models.
Additionally we also assessed the variability observed between groups on the above mentioned variables (feeding and saliva positive rates, AUC and peak viremia). For these comparisons, data were first assessed for normality using the Shapiro–Wilk test. If data from all groups were normally distributed, the Bartlett’s test of homogeneity of variance was used. If the data did not have a normal distribution, the Fligner-Killeen test was applied.
Survival of infected lambs (time to death) was compared between experiment groups using Kaplan–Meier survival analysis and the mortality rates were compared fitting a logistic regression model.
For all comparisons, the threshold for significance was p  More

1.
Dansgaard, W. et al. A new Greenland deep ice core. Science 218, 1273–1277 (1982).
ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 
2.
Stuiver, M. & Grootes, P. M. GISP2 oxygen isotope ratios. Quatern. Res. 53, 277–284 (2000).
ADS  CAS  Article  Google Scholar 

3.
Genty, D. et al. Isotopic characterization of rapid climatic events during OIS3 and OIS4 in Villars Cave stalagmites (SW-France) and correlation with Atlantic and Mediterranean pollen records. Quatern. Sci. Rev. 29, 2799–2820 (2010).
ADS  Article  Google Scholar 

4.
Pérez-Mejías, C. et al. Orbital-to-millennial scale climate variability during Marine Isotope Stages 5 to 3 in northeast Iberia. Quatern. Sci. Rev. 224, 105946 (2019).
Article  Google Scholar 

5.
Sánchez Goñi, M. F., Bard, E., Landais, A., Rossignol, L. & D’Errico, F. Air–sea temperature decoupling in western Europe during the last interglacial–glacial transition. Nat. Geosci. 6, 837 (2013).
ADS  Article  CAS  Google Scholar 

6.
Fletcher, W. J. et al. Millennial-scale variability during the last glacial in vegetation records from Europe. Quatern. Sci. Rev. 29, 2839–2864 (2010).
ADS  Article  Google Scholar 

7.
Hofreiter, M. & Stewart, J. Ecological change, range fluctuations and population dynamics during the pleistocene. Curr. Biol. 19, R584–R594 (2009).
CAS  PubMed  Article  PubMed Central  Google Scholar 

8.
Hodgkins, J. et al. Climate-mediated shifts in Neandertal subsistence behaviors at Pech de l’Azé IV and Roc de Marsal (Dordogne Valley, France). J. Hum. Evol. 96, 1–18 (2016).
PubMed  Article  PubMed Central  Google Scholar 

9.
Rendu, W. et al. Subsistence strategy changes during the Middle to Upper Paleolithic transition reveals specific adaptations of human populations to their environment. Sci. Rep. 9, 15817 (2019).
ADS  PubMed  PubMed Central  Article  CAS  Google Scholar 

10.
Dibble, H. L. et al. How did hominins adapt to ice age Europe without fire?. Curr. Anthropol. 58, S278–S287 (2017).
Article  Google Scholar 

11.
Sorensen, A. C. On the relationship between climate and Neandertal fire use during the Last Glacial in south-west France. Quatern. Int. 436, 114–128 (2017).
Article  Google Scholar 

12.
Delagnes, A. & Rendu, W. Shifts in Neandertal mobility, technology and subsistence strategies in western France. J. Archaeol. Sci. 38, 1771–1783 (2011).
Article  Google Scholar 

13.
Discamps, E., Jaubert, J. & Bachellerie, F. Human choices and environmental constraints: Deciphering the variability of large game procurement from Mousterian to Aurignacian times (MIS 5–3) in southwestern France. Quatern. Sci. Rev. 30, 2755–2775 (2011).
ADS  Article  Google Scholar 

14.
Faivre, J.-P. et al. The contribution of lithic production systems to the interpretation of Mousterian industrial variability in south-western France: The example of Combe-Grenal (Dordogne, France). Quatern. Int. 350, 227–240 (2014).
Article  Google Scholar 

15.
Hublin, J. J. The modern human colonization of western Eurasia: When and where?. Quatern. Sci. Rev. 118, 194–210 (2015).
ADS  Article  Google Scholar 

16.
Higham, T. et al. The timing and spatiotemporal patterning of Neanderthal disappearance. Nature 512, 306–309 (2014).
ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 

17.
Discamps, E. & Royer, A. Reconstructing palaeoenvironmental conditions faced by Mousterian hunters during MIS 5 to 3 in southwestern France: A multi-scale approach using data from large and small mammal communities. Quatern. Int. 433, 64–87 (2016).
Article  Google Scholar 

18.
Bernard, A. et al. Pleistocene seasonal temperature variations recorded in the δ 18O of Bison priscus teeth. Earth Planet. Sci. Lett. 283, 133–143 (2009).
ADS  CAS  Article  Google Scholar 

19.
Fabre, M. et al. Late Pleistocene climatic change in the French Jura (Gigny) recorded in the δ18O of phosphate from ungulate tooth enamel. Quatern. Res. 75, 605–613 (2011).
ADS  CAS  Article  Google Scholar 

20.
Richards, M. P. et al. Temporal variations in Equus tooth isotope values (C, N, O) from the Middle Paleolithic site of Combe Grenal, France (ca. 150,000 to 50,000 BP). J. Archaeol. Sci. Rep. 14, 189–198 (2017).
Google Scholar 

21.
Scherler, L., Tütken, T. & Becker, D. Carbon and oxygen stable isotope compositions of late Pleistocene mammal teeth from dolines of Ajoie (Northwestern Switzerland). Quatern. Res. (United States) 82, 378–387 (2014).
ADS  CAS  Article  Google Scholar 

22.
Skrzypek, G., Winiewski, A. & Grierson, P. F. How cold was it for Neanderthals moving to Central Europe during warm phases of the last glaciation?. Quatern. Sci. Rev. 30, 481–487 (2011).
ADS  Article  Google Scholar 

23.
Capitan, L. & Peyrony, D. Découverte d’un sixième squelette moustérien à La Ferrassie (Dordogne). Rev. Anthropol. 31, 382–388 (1921).
Google Scholar 

24.
Peyrony, D. La Ferrassie. Moustérien, Périgordien, Aurignacien. Préhistoire III. Préhistoire (1934)

25.
Turq, A. et al. La Ferrassie: Rapport d’opération pour l’année 2012 (2012).

26.
Delporte, H. & Delibrias, G. Le grand abri de la Ferrassie: fouilles 1968–1973. (Ed. du Laboratoire de paléontologie humaine et de préhistoire, 1984).

27.
Guérin, G. et al. A multi-method luminescence dating of the Palaeolithic sequence of La Ferrassie based on new excavations adjacent to the La Ferrassie 1 and 2 skeletons. J. Archaeol. Sci. 58, 147–166 (2015).
Article  Google Scholar 

28.
Talamo, S. et al. The new 14C chronology for the Palaeolithic site of La Ferrassie, France: The disappearance of Neanderthals and the arrival of Homo sapiens in France. J. Quatern. Sci. https://doi.org/10.1002/jqs.3236 (2020).
Article  Google Scholar 

29.
Britton, K. et al. Sampling plants and malacofauna in 87Sr/86Sr bioavailability studies: Implications for isoscape mapping and reconstructing of past mobility patterns. Front. Ecol. Evol. 8, 579473 (2020). 

30.
Willmes, M. et al. Mapping of bioavailable strontium isotope ratios in France for archaeological provenance studies. Appl. Geochem. 90, 75–86 (2018).
CAS  Article  Google Scholar 

31.
Deutscher Wetterdienst. Monthly mean air temperature of Gourdon, Dépt. Lot; Aquitaine/France (1996–2017) (2020).

32.
Hoppe, K. A. Correlation between the oxygen isotope ratio of North American bison teeth and local waters: Implication for paleoclimatic reconstructions. Earth Planet. Sci. Lett. 244, 408–417 (2006).
ADS  CAS  Article  Google Scholar 

33.
D’Angela, D. & Longinelli, A. Oxygen isotopes in living mammal’s bone phosphate: Further results. Chem. Geol. Isot. Geosci. Sect. 86, 75–82 (1990).
Article  Google Scholar 

34.
Rozanski, K., Araguás-Araguás, L. & Gonfiantini, R. Relation between long-term trends of oxygen-18 isotope composition of precipitation and source. Sci. New Ser. 258, 981–985 (1992).
CAS  Google Scholar 

35.
Levin, N. E., Cerling, T. E., Passey, B. H., Harris, J. M. & Ehleringer, J. R. A stable isotope aridity index for terrestrial environments. Proc. Natl. Acad. Sci. U.S.A. 103, 11201–11205 (2006).
ADS  CAS  PubMed  PubMed Central  Article  Google Scholar 

36.
Kohn, M. J., Schoeninger, M. J. & Valley, J. W. Herbivore tooth oxygen isotope compositions: Effects of diet and physiology. Geochim. Cosmochim. Acta 60, 3889–3896 (1996).
ADS  CAS  Article  Google Scholar 

37.
Bocherens, H., Koch, P. L., Mariotti, A., Geraads, D. & Jaeger, J.-J. Isotopic biogeochemistry (13C, 18O) of mammalian enamel from african pleistocene hominid sites. Palaios 11, 306–318 (1996).
ADS  Article  Google Scholar 

38.
Shackleton, N. Oxygen isotope analyses and Pleistocene temperatures re-assessed. Nature 215, 15–17 (1967).
ADS  CAS  Article  Google Scholar 

39.
Dansgaard, W. Stable isotopes in precipitation. Tellus 16, 436–468 (1964).
ADS  Article  Google Scholar 

40.
Gat, J. R. Isotope Hydrology: A Study of the Water Cycle Vol. 6 (Imperila College Press, London, 2010).
Google Scholar 

41.
Pryor, A. J. E., Stevens, R. E., O’Connell, T. C. & Lister, J. R. Quantification and propagation of errors when converting vertebrate biomineral oxygen isotope data to temperature for palaeoclimate reconstruction. Palaeogeogr. Palaeoclimatol. Palaeoecol. 412, 99–107 (2014).
Article  Google Scholar 

42.
Rozanski, K., Araguás-Araguás, L. & Gonfiantini, R. Isotopic patterns in modern global precipitation. Clim. Change Cont. Isot. Rec. 78, 1–36 (1993).
Google Scholar 

43.
Bocherens, H. et al. Direct isotopic evidence for subsistence variability in Middle Pleistocene Neanderthals (Payre, southeastern France). Quatern. Sci. Rev. 154, 226–236 (2016).
ADS  Article  Google Scholar 

44.
Ingraham, N. L., Caldwell, E. A. & Verhagen, B. T. Arid Catchments. in Isotope tracers in catchment hydrology, 435–465 (Elsevier, 1998). https://doi.org/10.1016/B978-0-444-81546-0.50020-3.

45.
Tütken, T., Furrer, H. & Walter Vennemann, T. Stable isotope compositions of mammoth teeth from Niederweningen, Switzerland: Implications for the Late Pleistocene climate, environment, and diet. Quatern. Int. 164–165, 139–150 (2007).
Article  Google Scholar 

46.
Ecker, M. et al. Middle pleistocene ecology and neanderthal subsistence: Insights from stable isotope analyses in Payre (Ardèche, southeastern France). J. Hum. Evol. 65, 363–373 (2013).
PubMed  Article  PubMed Central  Google Scholar 

47.
Stevens, R. E. et al. Nitrogen isotope analyses of reindeer (Rangifer tarandus), 45,000 BP to 9,000 BP: Palaeoenvironmental reconstructions. Palaeogeogr. Palaeoclimatol. Palaeoecol. 262, 32–45 (2008).
Article  Google Scholar 

48.
Stevens, R. E., Hermoso-Buxán, X. L., Marín-Arroyo, A. B., González-Morales, M. R. & Straus, L. G. Investigation of Late Pleistocene and Early Holocene palaeoenvironmental change at El Mirón cave (Cantabria, Spain): Insights from carbon and nitrogen isotope analyses of red deer. Palaeogeogr. Palaeoclimatol. Palaeoecol. 414, 46–60 (2014).
Article  Google Scholar 

49.
Drucker, D. G., Bridault, A., Hobson, K. A., Szuma, E. & Bocherens, H. Can carbon-13 in large herbivores reflect the canopy effect in temperate and boreal ecosystems? Evidence from modern and ancient ungulates. Palaeogeogr. Palaeoclimatol. Palaeoecol. 266, 69–82 (2008).
Article  Google Scholar 

50.
Diefendorf, A. F., Mueller, K. E., Wing, S. L., Koch, P. L. & Freeman, K. H. Global patterns in leaf 13C discrimination and implications for studies of past and future climate. Proc. Natl. Acad. Sci. 107, 5738–5743 (2010).
ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 

51.
Feranec, R. S., García, N., Díez, J. C. & Arsuaga, J. L. Understanding the ecology of mammalian carnivorans and herbivores from Valdegoba cave (Burgos, northern Spain) through stable isotope analysis. Palaeogeogr. Palaeoclimatol. Palaeoecol. 297, 263–272 (2010).
Article  Google Scholar 

52.
Bocherens, H., Drucker, D. G. & Madelaine, S. Evidence for a 15N positive excursion in terrestrial foodwebs at the Middle to Upper Palaeolithic transition in south-western France: Implications for early modern human palaeodiet and palaeoenvironment. J. Hum. Evol. 69, 31–43 (2014).
PubMed  Article  PubMed Central  Google Scholar 

53.
Sanchez Goni, M. F. et al. Contrasting impacts of Dansgaard-Oeschger events over a western European latitudinal transect modulated by orbital parameters. Quatern. Sci. Rev. 27, 1136–1151 (2008).
ADS  Article  Google Scholar 

54.
Ruddiman, W. F. & McIntyre, A. Oceanic mechanisms for amplification of the 23,000-year ice-volume cycle. Science 212, 617–627 (1981).
ADS  CAS  PubMed  Article  PubMed Central  Google Scholar 

55.
Kindler, P. et al. Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core. Clim. Past 10, 887–902 (2014).
Article  Google Scholar 

56.
Guiter, F. et al. The last climatic cycles in Western Europe: A comparison between long continuous lacustrine sequences from France and other terrestrial records. Quatern. Int. 111, 59–74 (2003).
Article  Google Scholar 

57.
de Beaulieu, J.-L. & Reille, M. The last climatic cycle at La Grande Pile (Vosges, France) a new pollen profile. Quatern. Sci. Rev. 11, 431–438 (1992).
ADS  Article  Google Scholar 

58.
Van Andel, T. H. & Tzedakis, P. C. Palaeolithic landscapes of Europe and environs, 150,000–25,000 years ago: An overview. Quatern. Sci. Rev. 15, 481–500 (1996).
ADS  Article  Google Scholar 

59.
Ponel, P. Rissian, Eemian and Würmian Coleoptera assemblages from La Grande Pile (Vosges, France). Palaeogeogr. Palaeoclimatol. Palaeoecol. 114, 1–41 (1995).
Article  Google Scholar 

60.
Royer, A. et al. Late Pleistocene (MIS 3–4) climate inferred from micromammal communities and δ 18O of rodents from Les Pradelles, France. Quatern. Res. (United States) 80, 113–124 (2013).
ADS  CAS  Article  Google Scholar 

61.
Barron, E., Andel, T. H. van & Pollard, D. Glacial environments II: Reconstructing the climate of Europe in the last glaciation. Neanderthals and modern humans in the European landscape during the last glaciation 57–78 (2003).

62.
Guérin, G. et al. Multi-method (TL and OSL), multi-material (quartz and flint) dating of the Mousterian site of Roc de Marsal (Dordogne, France): Correlating Neanderthal occupations with the climatic variability of MIS 5-3. J. Archaeol. Sci. 39, 3071–3084 (2012).
Article  CAS  Google Scholar 

63.
Copeland, S. R. et al. Strontium isotope ratios (87Sr/86Sr) of tooth enamel: A comparison of solution and laser ablation multicollector inductively coupled plasma mass spectrometry methods. Rapid Commun. Mass Spectrom. 22, 3187–3194 (2008).
ADS  CAS  PubMed  Article  PubMed Central  Google Scholar  More

Biodegradability of crude oil by two enrichment cultures
The enrichment culture H7S showed no obvious proliferation in the first five days because sample 7S was a sulphide rock, while H11S showed visible proliferation after the fourth day. After 14 days of cultivation, gravimetric analysis demonstrated that the enrichment cultures H7S and H11S exhibited similar oil-degrading abilities and degraded 54% and 56% of the crude oil, respectively (Fig. 1).
Figure 1

The oil degradation efficiency of the two enrichment cultures H7S and H11S.

Full size image

The biodegradation percentages for total n-alkanes (C10–C34) and polycyclic aromatic hydrocarbons (PAHs) were calculated by comparing the two enrichment cultures with the negative controls (Fig. 2). Based on evaluation with C17/pristane and C18/phytane, the degradation efficiencies of the two enrichment cultures were significantly better than those of the negative controls (P  More