Ecology

1.
Aschoff, J. Circadian activity pattern with two peaks. Ecology 47, 657–662 (1966).
Article  Google Scholar 
2.
Daan, S. & Aschoff, J. Circadian rhythms of locomotor activity in captive birds and mammals: Their variations with season and latitude. Oecologia 18, 269–316 (1975).
ADS  Article  PubMed  PubMed Central  Google Scholar 

3.
Katandukila, J. V., Bennett, N. C., Chimimba, C. T., Faulkes, C. G. & Oosthuizen, M. K. Locomotor activity patterns of captive East African root rats, Tachyoryctes splendens (Rodentia: Spalacidae), from Tanzania, East Africa. J. Mamm. 94, 1393–1400 (2013).
Article  Google Scholar 

4.
Bennie, J. J., Duffy, J. P., Inger, R. & Gaston, K. J. Biogeography of time partitioning in mammals. Proc. Natl. Acad. Sci. U.S.A. 111, 13727–13732 (2014).
ADS  CAS  Article  PubMed  PubMed Central  Google Scholar 

5.
Cloudsley-Thompson, J. L. Rhythmic Activity in Animal Physiology and Behaviour (Academic Press, Cambridge, 1961).
Google Scholar 

6.
Aschoff, J., Gercke, H., Pohl, P., Rieger, P. V. & Wever, S. P. U. R. Interdependent parameters of circadian activity rhythms in birds and man. In Biochronometry (ed. Menaker, M.) 3–29 (National Academy of Science, Washington, DC, 1971).
Google Scholar 

7.
Risenhoover, K. L. Winter activity patterns of moose in interior Alaska. J. Wildl. Manage. 50, 727–734 (1986).
Article  Google Scholar 

8.
Castillo-Ruiz, A., Paul, M. J. & Schwartz, W. J. In search of a temporal niche: Social interactions. Prog. Brain Res. 199, 267–280 (2012).
Article  PubMed  PubMed Central  Google Scholar 

9.
Kronfeld-Schor, N., Visser, M. E., Salis, L. & van Gils, J. A. Chronobiology of interspecific interactions in a changing world. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. https://doi.org/10.1098/rstb.2016.0248 (2017).
Article  Google Scholar 

10.
Farsi, H. et al. Validation of locomotion scoring as a new and inexpensive technique to record circadian locomotor activity in large mammals. Heliyon 4, e00980–e00980 (2018).
CAS  Article  PubMed  PubMed Central  Google Scholar 

11.
El Allali, K. et al. Smartphone and a freely available application as a new tool to record locomotor activity rhythm in large mammals and humans. Chronobiol. Int. 36, 1047–1057 (2019).
Article  PubMed  Google Scholar 

12.
Schmidt-Nielsen, K., Schmidt-Nielsen, B., Jarnum, S. A. & Houpt, T. R. Body temperature of the camel and its relation to water economy. Am. J. Physiol. 188, 103–112 (1957).
CAS  Article  PubMed  Google Scholar 

13.
Bouâouda, H. et al. Daily regulation of body temperature rhythm in the camel (Camelus dromedarius) exposed to experimental desert conditions. Physiol. Rep. 2, e12151 (2014).
Article  PubMed  PubMed Central  Google Scholar 

14.
El Allali, K. et al. Entrainment of the circadian clock by daily ambient temperature cycles in the camel (Camelus dromedarius). Am. J. Physiol. Regul. Integr. Comp. Physiol. 304, R1044–R1052 (2013).
Article  CAS  PubMed  Google Scholar 

15.
Farsi, H. et al. Melatonin rhythm and other outputs of the master circadian clock in the desert goat (Capra hircus) are entrained by daily cycles of ambient temperature. J. Pineal Res. 68, e12634 (2020).
CAS  Article  PubMed  Google Scholar 

16.
Ebling, F. J., Lincoln, G. A., Wollnik, F. & Anderson, N. Effects of constant darkness and constant light on circadian organization and reproductive responses in the ram. J. Biol. Rhythms 3, 365–384 (1988).
CAS  Article  PubMed  Google Scholar 

17.
Johnson, R. F., Randall, S. & Randall, W. Freerunning and entrained circadian rhythms in activity, eating and drinking in the cat. J. Interdiscipl. Cycle Res. 14, 315–327 (1983).
Article  Google Scholar 

18.
Jilge, B., Hörnicke, H. & Stähle, H. Circadian rhythms of rabbits during restrictive feeding. Am. J. Physiol. 253, R46–R54 (1987).
CAS  PubMed  Google Scholar 

19.
Decoursey, G. & Decoursey, P. J. Adaptive aspects of activity rhythms in bats. Biol. Bull. 126, 14–27 (1964).
Article  Google Scholar 

20.
Erkert, H. G., Nagel, B. & Stephani, I. Light and social effects on the free-running circadian activity rhythm in common marmosets (Callithrix jacchus; Primates): Social masking, pseudo-splitting, and relative coordination. Behav. Ecol. Sociobiol. 18, 443–452 (1986).
Article  Google Scholar 

21.
O’Reilly, H., Armstrong, S. M. & Coleman, G. J. Restricted feeding and circadian activity rhythms of a predatory marsupial, Dasyuroides byrnei. Physiol. Behav. 38, 471–476 (1986).
CAS  Article  PubMed  Google Scholar 

22.
Boulos, Z., Frim, D. M., Dewey, L. K. & Moore-Ede, M. C. Effects of restricted feeding schedules on circadian organization in squirrel monkeys. Physiol. Behav. 45, 507–515 (1989).
CAS  Article  PubMed  Google Scholar 

23.
Mahoney, M., Bult, A. & Smale, L. Phase response curve and light-induced fos expression in the suprachiasmatic nucleus and adjacent hypothalamus of Arvicanthis niloticus. J. Biol. Rhythms 16, 149–162 (2001).
CAS  Article  PubMed  Google Scholar 

24.
Alagaili, A. N., Bennett, N. C., Amor, N. M. & Hart, D. W. The locomotory activity patterns of the arid-dwelling desert hedgehog, Paraechinus aethiopicus, from Saudi Arabia. J. Arid Environ. 177, 104141 (2020).
ADS  Article  Google Scholar 

25.
Verwey, M., Robinson, B. & Amir, S. Recording and analysis of circadian rhythms in running-wheel activity in rodents. J. Vis. Exp. https://doi.org/10.3791/50186 (2013).
Article  PubMed  PubMed Central  Google Scholar 

26.
Refinetti, R. Early research on circadian rhythms. In Circadian Physiology 2nd edn (ed. Refinetti, R.) 1–667 (CRC Taylor and Frabcis Group, Boca Raton, 2006).
Google Scholar 

27.
Goldman, B. D., Goldman, S. L., Riccio, A. P. & Terkel, J. Circadian patterns of locomotor activity and body temperature in blind mole-rats, Spalax ehrenbergi. J. Biol. Rhythms 12, 348–361 (1997).
CAS  Article  PubMed  Google Scholar 

28.
Kopp, C. et al. Effects of a daylight cycle reversal on locomotor activity in several inbred strains of mice. Physiol. Behav. 63, 577–585 (1998).
CAS  Article  PubMed  Google Scholar 

29.
Giannetto, C., Casella, S., Caola, G. & Piccione, G. Photic and non-photic entrainment on daily rhythm of locomotor activity in goats. Anim. Sci. J. 81, 122–128 (2010).
Article  PubMed  Google Scholar 

30.
Piccione, G., Giannetto, C., Casella, S. & Caola, G. Daily locomotor activity in five domestic animals. Anim. Biol. 60, 15–24 (2010).
Article  Google Scholar 

31.
Challet, E. Minireview: Entrainment of the suprachiasmatic clockwork in diurnal and nocturnal mammals. Endocrinology 148, 5648–5655 (2007).
CAS  Article  PubMed  PubMed Central  Google Scholar 

32.
Dibner, C., Schibler, U. & Albrecht, U. The mammalian circadian timing system: Organization and coordination of central and peripheral clocks. Annu. Rev. Physiol. 72, 517–549 (2010).
CAS  Article  PubMed  PubMed Central  Google Scholar 

33.
Tanaka, M., Ichitani, Y., Okamura, H., Tanaka, Y. & Ibata, Y. The direct retinal projection to VIP neuronal elements in the rat SCN. Brain Res. Bull. 31, 637–640 (1993).
CAS  Article  PubMed  PubMed Central  Google Scholar 

34.
Jacomy, H., Burlet, A. & Bosler, O. Vasoactive intestinal peptide neurons as synaptic targets for vasopressin neurons in the suprachiasmatic nucleus. Double-label immunocytochemical demonstration in the rat. Neuroscience 88, 859–870 (1999).
CAS  Article  PubMed  PubMed Central  Google Scholar 

35.
Aton, S. J., Colwell, C. S., Harmar, A. J., Waschek, J. & Herzog, E. D. Vasoactive intestinal polypeptide mediates circadian rhythmicity and synchrony in mammalian clock neurons. Nat. Neurosci. 8, 476–483 (2005).
CAS  Article  PubMed  PubMed Central  Google Scholar 

36.
Reppert, S. M. & Weaver, D. R. Comparing clockworks: Mouse versus fly. J. Biol. Rhythms 15, 357–364 (2000).
CAS  Article  PubMed  PubMed Central  Google Scholar 

37.
Reppert, S. M. & Weaver, D. R. Molecular analysis of mammalian circadian rhythms. Annu. Rev. Physiol. 63, 647–676 (2001).
CAS  Article  PubMed  PubMed Central  Google Scholar 

38.
Shearman, L. P. et al. Interacting molecular loops in the mammalian circadian clock. Science 288, 1013–1019 (2000).
ADS  CAS  Article  PubMed  PubMed Central  Google Scholar 

39.
Okamura, H., Yamaguchi, S. & Yagita, K. Molecular machinery of the circadian clock in mammals. Cell Tissue Res. 309, 47–56 (2002).
CAS  Article  Google Scholar 

40.
Takahashi, J. S., Hong, H. K., Ko, C. H. & McDearmon, E. L. The genetics of mammalian circadian order and disorder: Implications for physiology and disease. Nat. Rev. Genet. 9, 764–775 (2008).
CAS  Article  PubMed  PubMed Central  Google Scholar 

41.
Mohawk, J. A., Green, C. B. & Takahashi, J. S. Central and peripheral circadian clocks in mammals. Annu. Rev. Neurosci. 35, 445–462 (2012).
CAS  Article  PubMed  PubMed Central  Google Scholar 

42.
Rensing, L. & Ruoff, P. Temperature effect on entrainment, phase shifting, and amplitude of circadian clocks and its molecular bases. Chronobiol. Int. 19, 807–864 (2002).
CAS  Article  PubMed  Google Scholar 

43.
Aschoff, J. & Tokura, H. Circadian activity rhythms in squirrel monkeys: Entrainment by temperature cycles 1. J. Biol. Rhythms 1, 91–99 (1986).
CAS  Article  PubMed  Google Scholar 

44.
Pálková, M., Sigmund, L. & Erkert, H. G. Effect of ambient temperature on the circadian activity rhythm in common marmosets, Callithrix j jacchus (primates). Chronobiol. Int. 16, 149–161 (1999).
Article  PubMed  Google Scholar 

45.
Rajaratnam, S. M. W. & Redman, J. R. Entrainment of activity rhythms to temperature cycles in diurnal palm squirrels. Physiol. Behav. 63, 271–277 (1998).
CAS  Article  PubMed  Google Scholar 

46.
Refinetti, R. Entrainment of circadian rhythm by ambient temperature cycles in mice. J. Biol. Rhythms 25, 247–256 (2010).
Article  Google Scholar 

47.
van Jaarsveld, B., Bennett, N. C., Hart, D. W. & Oosthuizen, M. K. Locomotor activity and body temperature rhythms in the Mahali mole-rat (C. h. mahali): The effect of light and ambient temperature variations. J. Therm. Biol. 79, 24–32 (2019).
Article  Google Scholar 

48.
Schmidt-Nielsen, K. The physiology of the camel. Sci. Am. 201, 140–151 (1959).
CAS  Article  PubMed  PubMed Central  Google Scholar 

49.
Wu, H. et al. Camelid genomes reveal evolution and adaptation to desert environments. Nat. Commun. 5, 5188 (2014).
ADS  CAS  Article  PubMed  PubMed Central  Google Scholar 

50.
Samara, E. M. Unraveling the relationship between the topographic distribution patterns of skin temperature and perspiration response in dromedary camels. J. Therm. Biol 84, 311–315 (2019).
Article  PubMed  Google Scholar 

51.
Tibary, A. & El Allali, K. Dromedary camel: A model of heat resistant livestock animal. Theriogenology 154, 203–211 (2020).
Article  PubMed  Google Scholar 

52.
Lindberg, R. G. & Hayden, P. Thermoperiodic entrainment of arousal from torpor in the little pocket mouse, Perognathus longimembris. Chronobiologia 1, 356–361 (1974).
CAS  PubMed  Google Scholar 

53.
Erkert, H. G. & Rothmund, E. Differences in temperature sensitivity of the circadian systems of homoeothermic and heterothermic neotropical bats. Comp. Biochem. Physiol. 68A, 383–390 (1980).
Google Scholar 

54.
Pohl, H. Temperature cycles as zeitgeber for the circadian clock of two burrowing rodents, the normothermic antelope ground squirrel and the heterothermic Syrian Hamster. Biol. Rhythm Res. 29, 311–325 (1998).
Article  Google Scholar 

55.
Cain, J. W., Krausman, P. R., Rosenstock, S. S. & Turner, J. C. Mechanisms of thermoregulation and water balance in Desert Ungulates. Wildl. Soc. Bull. 1973–2006(34), 570–581 (2006).
Article  Google Scholar 

56.
Mengistu, U., Dahlborn, K. & Olsson, K. Mechanisms of water economy in lactating Ethiopian Somali goats during repeated cycles of intermittent watering. Anim. Int. J. Anim. Biosci. 1, 1009–1017 (2007).
CAS  Article  Google Scholar 

57.
Gauthier-Pilters, H. Aspects of dromedary ecology and ethology. In The Camelid (ed. Cockrill, W. R.) (Scandinavian Institute of African Studies, Uppsala, 1984).
Google Scholar 

58.
Miller, G. D., Cochran, M. H. & Smith, E. L. Nighttime activity of desert bighorn sheep. Desert Bighorn Council Trans. 28, 23–25 (1984).
Google Scholar 

59.
Hayes, C. L. & Krausman, P. R. Nocturnal activity of female desert mule deer. J. Wildl. Manage. 57, 897–904 (1993).
Article  Google Scholar 

60.
Davimes, J. G. et al. Temporal niche switching in Arabian oryx (Oryx leucoryx): Seasonal plasticity of 24h activity patterns in a large desert mammal. Physiol. Behav. 177, 148–154 (2017).
CAS  Article  PubMed  Google Scholar 

61.
Davimes, J. G. et al. Seasonal variations in sleep of free-ranging Arabian oryx (Oryx leucoryx) under natural hyperarid conditions. Sleep https://doi.org/10.1093/sleep/zsy038 (2018).
Article  PubMed  Google Scholar 

62.
El Allali, K. et al. Seasonal variations in the nycthemeral rhythm of plasma melatonin in the camel (Camelus dromedarius). J. Pineal Res. 39, 121–128 (2005).
Article  CAS  PubMed  Google Scholar 

63.
Mrosovsky, N. Circannual cycles in golden-mantled ground squirrels: Phase shift produced by low temperature. J. Comp. Physiol. 136, 349–353 (1980).
Article  Google Scholar 

64.
Mrosovsky, N. Circannual cycles in golden-mantled ground squirrels: Experiments with food deprivation and effects of temperature on periodicity. J. Comp. Physiol. 136, 355–360 (1980).
Article  Google Scholar 

65.
Mrosovsky, N. Thermal effects on the periodicity, phasing and peristance of circannual cycles. In Living in the Cold (eds Heller, H. C. et al.) 403–410 (Elsevier, New York, 1986).
Google Scholar 

66.
Mrosovsky, N. Circannual cycles in golden-mantled ground squirrels: fall and spring cold pulses. J. Comp. Physiol. 167, 683–689 (1990).
Article  Google Scholar 

67.
Canguilhem, B., Schieber, J. P. & Koch, A. Circannual weight rhythm of the European hamster (Cricetus cricetus). Respective influence of the photoperiod and external temperature during its course. Arch. Sci. Physiol. 27, 67–90 (1973).
CAS  Google Scholar 

68.
Jallageas, M. & Assenmacher, I. External factors controlling annual testosterone and thyroxine cycles in the edible dormouse Glis glis. Comp. Biochem. Physiol. A Physiol. 77, 161–167 (1984).
CAS  Article  Google Scholar 

69.
Touitou, Y., Smolensky, M. H. & Portaluppi, F. Ethics, standards, and procedures of animal and human chronobiology research. Chronobiol. Int. 23, 1083–1096 (2006).
Article  PubMed  Google Scholar  More

1.
Claussnitzer, M. et al. A brief history of human disease genetics. Nature 577, 179–189 (2020).
ADS  CAS  PubMed  PubMed Central  Google Scholar 
2.
Hiller, M. et al. A “forward genomics” approach links genotype to phenotype using independent phenotypic losses among related species. Cell Rep. 2, 817–823 (2012).
CAS  PubMed  PubMed Central  Google Scholar 

3.
Wasser, S. K. et al. Genetic assignment of large seizures of elephant ivory reveals Africa’s major poaching hotspots. Science 349, 84–87 (2015).
ADS  CAS  PubMed  PubMed Central  Google Scholar 

4.
Wright, B. et al. Development of a SNP-based assay for measuring genetic diversity in the Tasmanian devil insurance population. BMC Genomics 16, 791 (2015).
PubMed  PubMed Central  Google Scholar 

5.
Lappalainen, T., Scott, A. J., Brandt, M. & Hall, I. M. Genomic analysis in the age of human genome sequencing. Cell 177, 70–84 (2019).
CAS  PubMed  PubMed Central  Google Scholar 

6.
Kircher, M. et al. A general framework for estimating the relative pathogenicity of human genetic variants. Nat. Genet. 46, 310–315 (2014).
CAS  PubMed  PubMed Central  Google Scholar 

7.
Lindblad-Toh, K. et al. A high-resolution map of human evolutionary constraint using 29 mammals. Nature 478, 476–482 (2011).
CAS  PubMed  PubMed Central  Google Scholar 

8.
Finucane, H. K. et al. Partitioning heritability by functional annotation using genome-wide association summary statistics. Nat. Genet. 47, 1228–1235 (2015).
CAS  PubMed  PubMed Central  Google Scholar 

9.
IUCN. The IUCN Red List of Threatened Species. Version 2019-2 http://www.iucnredlist.org (2019).

10.
Ryder, O. A. & Onuma, M. Viable cell culture banking for biodiversity characterization and conservation. Annu. Rev. Anim. Biosci. 6, 83–98 (2018).
PubMed  Google Scholar 

11.
Weisenfeld, N. I. et al. Comprehensive variation discovery in single human genomes. Nat. Genet. 46, 1350–1355 (2014).
CAS  PubMed  PubMed Central  Google Scholar 

12.
Putnam, N. H. et al. Chromosome-scale shotgun assembly using an in vitro method for long-range linkage. Genome Res. 26, 342–350 (2016).
CAS  PubMed  PubMed Central  Google Scholar 

13.
Kim, J. et al. Reconstruction and evolutionary history of eutherian chromosomes. Proc. Natl Acad. Sci. USA 114, E5379–E5388 (2017).
CAS  PubMed  Google Scholar 

14.
Lek, M. et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature 536, 285–291 (2016).
CAS  PubMed  PubMed Central  Google Scholar 

15.
Balasubramanian, S. et al. Using ALoFT to determine the impact of putative loss-of-function variants in protein-coding genes. Nat. Commun. 8, 382 (2017).
ADS  PubMed  PubMed Central  Google Scholar 

16.
Meadows, J. R. S. & Lindblad-Toh, K. Dissecting evolution and disease using comparative vertebrate genomics. Nat. Rev. Genet. 18, 624–636 (2017).
CAS  PubMed  Google Scholar 

17.
Cooper, G. M. & Shendure, J. Needles in stacks of needles: finding disease-causal variants in a wealth of genomic data. Nat. Rev. Genet. 12, 628–640 (2011).
CAS  PubMed  Google Scholar 

18.
Baiz, M. D., Tucker, P. K., Mueller, J. L. & Cortés-Ortiz, L. X-linked signature of reproductive isolation in humans is mirrored in a howler monkey hybrid zone. J. Hered. 111, 419–428 (2020).
PubMed  PubMed Central  Google Scholar 

19.
Dobzhansky, T. & Dobzhansky, T. G. Genetics and the Origin of Species (Columbia Univ. Press, 1937).

20.
Herrera-Álvarez, S., Karlsson, E., Ryder, O. A., Lindblad-Toh, K. & Crawford, A. J. How to make a rodent giant: genomic basis and tradeoffs of gigantism in the capybara, the world’s largest rodent. Preprint at https://doi.org/10.1101/424606 (2018).

21.
Abegglen, L. M. et al. Potential mechanisms for cancer resistance in elephants and comparative cellular response to DNA damage in humans. J. Am. Med. Assoc. 314, 1850–1860 (2015).
CAS  Google Scholar 

22.
Casewell, N. R. et al. Solenodon genome reveals convergent evolution of venom in eulipotyphlan mammals. Proc. Natl Acad. Sci. USA 116, 25745–25755 (2019).
CAS  PubMed  Google Scholar 

23.
Beichman, A. C. et al. Aquatic adaptation and depleted diversity: a deep dive into the genomes of the sea otter and giant otter. Mol. Biol. Evol. 36, 2631–2655 (2019).
CAS  PubMed  Google Scholar 

24.
Damas, J. et al. Broad host range of SARS-CoV-2 predicted by comparative and structural analysis of ACE2 in vertebrates. Proc. Natl Acad. Sci. USA 117, 22311–22322 (2020).
CAS  PubMed  Google Scholar 

25.
Xue, Y. et al. Mountain gorilla genomes reveal the impact of long-term population decline and inbreeding. Science 348, 242–245 (2015).
ADS  CAS  PubMed  PubMed Central  Google Scholar 

26.
Ceballos, F. C., Joshi, P. K., Clark, D. W., Ramsay, M. & Wilson, J. F. Runs of homozygosity: windows into population history and trait architecture. Nat. Rev. Genet. 19, 220–234 (2018).
CAS  PubMed  Google Scholar 

27.
Spielman, D., Brook, B. W. & Frankham, R. Most species are not driven to extinction before genetic factors impact them. Proc. Natl Acad. Sci. USA 101, 15261–15264 (2004).
ADS  CAS  PubMed  Google Scholar 

28.
Vinson, J. P. et al. Assembly of polymorphic genomes: algorithms and application to Ciona savignyi. Genome Res. 15, 1127–1135 (2005).
PubMed  PubMed Central  Google Scholar 

29.
MacManes, M. D. & Lacey, E. A. The social brain: transcriptome assembly and characterization of the hippocampus from a social subterranean rodent, the colonial tuco-tuco (Ctenomys sociabilis). PLoS ONE 7, e45524 (2012).
ADS  CAS  PubMed  PubMed Central  Google Scholar 

30.
Jones, K. E. et al. PanTHERIA: a species-level database of life history, ecology, and geography of extant and recently extinct mammals. Ecology 90, 2648 (2009).
Google Scholar 

31.
Cardillo, M. Biological determinants of extinction risk: why are smaller species less vulnerable? Anim. Conserv. 6, 63–69 (2003).
Google Scholar 

32.
Natesh, M. et al. Empowering conservation practice with efficient and economical genotyping from poor quality samples. Methods Ecol. Evol. 10, 853–859 (2019).
PubMed  PubMed Central  Google Scholar 

33.
Lowry, D. B. et al. Breaking RAD: an evaluation of the utility of restriction site-associated DNA sequencing for genome scans of adaptation. Mol. Ecol. Resour. 17, 142–152 (2017).
CAS  PubMed  Google Scholar 

34.
Shapiro, B. Pathways to de-extinction: how close can we get to resurrection of an extinct species? Funct. Ecol. 31, 996–1002 (2017).
Google Scholar 

35.
Benazzo, A. et al. Survival and divergence in a small group: the extraordinary genomic history of the endangered Apennine brown bear stragglers. Proc. Natl Acad. Sci. USA 114, E9589–E9597 (2017).
CAS  PubMed  Google Scholar 

36.
Saremi, N. F. et al. Puma genomes from North and South America provide insights into the genomic consequences of inbreeding. Nat. Commun. 10, 4769 (2019).
ADS  PubMed  PubMed Central  Google Scholar 

37.
Armstrong, J. et al. Progressive Cactus is a multiple-genome aligner for the thousand-genome era. Nature https://doi.org/10.1038/s41586-020-2871-y (2020).

38.
Haeussler, M. et al. The UCSC genome browser database: 2019 update. Nucleic Acids Res. 47, D853–D858 (2019).
CAS  PubMed  Google Scholar 

39.
Rands, C. M., Meader, S., Ponting, C. P. & Lunter, G. 8.2% of the human genome is constrained: variation in rates of turnover across functional element classes in the human lineage. PLoS Genet. 10, e1004525 (2014).
PubMed  PubMed Central  Google Scholar 

40.
ENCODE Project Consortium. An integrated encyclopedia of DNA elements in the human genome. Nature 489, 57–74 (2012).
ADS  Google Scholar 

41.
GTEx Consortium. Genetic effects on gene expression across human tissues. Nature 550, 204–213 (2017).
PubMed Central  Google Scholar 

42.
Regev, A. et al. The human cell atlas. eLife 6, e27041 (2017).
PubMed  PubMed Central  Google Scholar 

43.
Lewin, H. A. et al. Earth BioGenome project: sequencing life for the future of life. Proc. Natl Acad. Sci. USA 115, 4325–4333 (2018).
CAS  PubMed  Google Scholar 

44.
Koepfli, K.-P., Paten, B., the Genome 10K Community of Scientists & O’Brien, S. J. The Genome 10K project: a way forward. Annu. Rev. Anim. Biosci. 3, 57–111 (2015).
CAS  PubMed  PubMed Central  Google Scholar 

45.
Teeling, E. C. et al. Bat biology, genomes, and the Bat1K project: to generate chromosome-level genomes for all living bat species. Annu. Rev. Anim. Biosci. 6, 23–46 (2018).
PubMed  Google Scholar 

46.
Feng, S. et al. Dense sampling of bird diversity increases power of comparative genomics. Nature https://doi.org/10.1038/s41586-020-2873-9 (2020).

47.
Kumar, S., Stecher, G., Suleski, M. & Hedges, S. B. TimeTree: a resource for timelines, timetrees, and divergence times. Mol. Biol. Evol. 34, 1812–1819 (2017).
CAS  PubMed  Google Scholar 

48.
Wilson, D. E. & Reeder, D. M. (eds) Mammal Species of the World. A Taxonomic and Geographic Reference 3rd edn (Johns Hopkins Univ. Press, 2005).

49.
Vlieghe, D. et al. A new generation of JASPAR, the open-access repository for transcription factor binding site profiles. Nucleic Acids Res. 34, D95–D97 (2006).
CAS  PubMed  Google Scholar 

50.
Simão, F. A., Waterhouse, R. M., Ioannidis, P., Kriventseva, E. V. & Zdobnov, E. M. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210–3212 (2015).
PubMed  PubMed Central  Google Scholar 

51.
Farré, M. et al. A near-chromosome-scale genome assembly of the gemsbok (Oryx gazella): an iconic antelope of the Kalahari desert. Gigascience 8, giy162 (2019).
PubMed Central  Google Scholar 

52.
McKenna, A. et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303 (2010).
CAS  PubMed  PubMed Central  Google Scholar 

53.
DePristo, M. A. et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 43, 491–498 (2011).
CAS  PubMed  PubMed Central  Google Scholar 

54.
Li, H. et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics 25, 2078–2079 (2009).
PubMed  PubMed Central  Google Scholar 

55.
Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics 25, 1754–1760 (2009).
CAS  PubMed  PubMed Central  Google Scholar 

56.
Benaglia, T., Chauveau, D., Hunter, D. & Young, D. mixtools: an R package for analyzing finite mixture models. J. Stat. Softw. 32, 1–29 (2009).
Google Scholar 

57.
R Core Team. R: A Language and Environment for Statistical Computing. https://www.R-project.org/ (2019).

58.
Paten, B. et al. Cactus: algorithms for genome multiple sequence alignment. Genome Res. 21, 1512–1528 (2011).
CAS  PubMed  PubMed Central  Google Scholar 

59.
Ondov, B. D. et al. Mash: fast genome and metagenome distance estimation using MinHash. Genome Biol. 17, 132 (2016).
PubMed  PubMed Central  Google Scholar 

60.
Smit, A. F. A., Hubley, R. & Green, P. RepeatMasker Open-4.0. http://www.repeatmasker.org/ (2013–2015).

61.
Hubisz, M. J., Pollard, K. S. & Siepel, A. PHAST and RPHAST: phylogenetic analysis with space/time models. Brief. Bioinform. 12, 41–51 (2011).
CAS  PubMed  Google Scholar 

62.
Benjamini, Y. & Hochberg, Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. B 57, 289–300 (1995).
MathSciNet  MATH  Google Scholar 

63.
Benson, D. A. et al. GenBank. Nucleic Acids Res. 41, D36–D42 (2013).
CAS  PubMed  Google Scholar 

64.
Nguyen, N. et al. Comparative assembly hubs: web-accessible browsers for comparative genomics. Bioinformatics 30, 3293–3301 (2014).
CAS  PubMed  PubMed Central  Google Scholar 

65.
Karczewski, K. J. et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581, 434–443 (2020).
ADS  CAS  PubMed  PubMed Central  Google Scholar 

66.
Pinheiro, E. C., Taddei, V. A., Migliorini, R. H. & Kettelhut, I. C. Effect of fasting on carbohydrate metabolism in frugivorous bats (Artibeus lituratus and Artibeus jamaicensis). Comp. Biochem. Physiol. B Biochem. Mol. Biol. 143, 279–284 (2006).
PubMed  Google Scholar 

67.
Gordon, L. M. et al. Amorphous intergranular phases control the properties of rodent tooth enamel. Science 347, 746–750 (2015).
ADS  CAS  PubMed  Google Scholar 

68.
Hindle, A. G. & Martin, S. L. Intrinsic circannual regulation of brown adipose tissue form and function in tune with hibernation. Am. J. Physiol. Endocrinol. Metab. 306, E284–E299 (2014).
CAS  PubMed  Google Scholar 

69.
Stanford, K. I. et al. Brown adipose tissue regulates glucose homeostasis and insulin sensitivity. J. Clin. Invest. 123, 215–223 (2013).
CAS  PubMed  Google Scholar 

70.
Chondronikola, M. et al. Brown adipose tissue improves whole-body glucose homeostasis and insulin sensitivity in humans. Diabetes 63, 4089–4099 (2014).
CAS  PubMed  PubMed Central  Google Scholar 

71.
Saito, M. et al. High incidence of metabolically active brown adipose tissue in healthy adult humans: effects of cold exposure and adiposity. Diabetes 58, 1526–1531 (2009).
CAS  PubMed  PubMed Central  Google Scholar  More

Listen to the latest from the world of science, with Benjamin Thompson and Nick Howe.
Your browser does not support the audio element.
Download MP3

In this episode:
00:46 Sensory pollution and bird reproduction
Light- and noise-pollution have been shown to affect the behaviour of birds. However, it’s been difficult to work out whether these behavioural changes have led to bird species thriving or declining. Now, researchers have assembled a massive dataset that can begin to give some answers. Research article: Senzaki et al.
10:17 Coronapod
Interim results from a phase III trial show compelling evidence that a coronavirus vaccine candidate can prevent COVID-19. However, amid the optimism there remain questions to be answered – we discuss these, and what the results might mean for other vaccines in development. News: What Pfizer’s landmark COVID vaccine results mean for the pandemic
23:29 Research Highlights
A tiny bat breaks a migration record, and researchers engineer a mouse’s sense of place. Research Highlight: The record-setting flight of a bat that weighs less than a toothbrush; Research Article: Robinson et al.
25:39 Organised crime in fisheries
When you think of fishing, organised crime probably isn’t the first thing that springs to mind. However, billions of dollars every year from the fishing industry are lost to criminal enterprises. We discuss some of the impacts and what can be done about it. Research Article: Witbooi et al.
32:13 Briefing Chat
We discuss some highlights from the Nature Briefing. This time, a time-capsule discovered on the Irish coast provides a damning indictment of Arctic warming, and some human remains challenge the idea of ‘man-the-hunter’. The Guardian: Arctic time capsule from 2018 washes up in Ireland as polar ice melts; Science: Woman the hunter: Ancient Andean remains challenge old ideas of who speared big game
Subscribe to Nature Briefing, an unmissable daily round-up of science news, opinion and analysis free in your inbox every weekday.
Never miss an episode: Subscribe to the Nature Podcast on Apple Podcasts, Google Podcasts, Spotify or your favourite podcast app. Head here for the Nature Podcast RSS feed. More

Species distribution modelling
A total of 1548 occurrence locations were obtained for B. fagacearum from the United States Forest Service (provided to us by Erin Bullas-Appleton of the Canadian Food Inspection Agency in July 2018) and the Global Biodiversity Information Facility17. These records were filtered at a 300 arcsecond (approximately 10-km) resolution to remove duplicates and reduce spatial clustering, leaving 1401 unique location records (Fig. 1a). Occurrence locations for two key insect vectors of oak wilt, C. truncatus (Fig. 2a) and C. sayi (Fig. 3a), were obtained from GBIF, publications10,18, and from specimens in the following collections: Atlantic Forestry Centre, Fredericton, NB, Canada; Canadian National Collection of Insects, Arachnids, and Nematodes, Agriculture and Agri-Food Canada Research Centre, Ottawa, ON, Canada; Ontario Forest Research Institute, Sault Ste. Marie, ON, Canada; Gareth S. Powell Collection, Nephi, UT, USA; Reginald Webster Collection, Charters Settlement, New Brunswick, Canada; and the Florida State Collection of Arthropods, Gainesville, FL, USA . After filtering at a 10-km resolution, there were 82 and 58 unique occurrence records for C. truncatus and C. sayi respectively.
Figure 1

Occurrence data (a) used for generating climate suitability models for Bretziella fagacearum. Maps with colour gradients indicate Maxent-derived climate suitability for B. fagacearum for the: 1981–2010 period (b); 2011–2040 period (c); and 2041–2070 period (d). Stippling delineates the ANUCLIM-derived climate envelope for B. fagacearum in each time period. Hatching delineates the current distribution of Quercus in Canada. Climate projections are based on a composite of four climate models and the RCP 4.5 emissions scenario (see text for further details). Maps were generated using ARCGIS v.9.3 (ESRI, Redlands, CA, USA; https://www.esri.com/arcgis/about-arcgis).

Full size image

Figure 2

Occurrence data (a) used for generating climate suitability models for Colopterus truncatus. Maps with colour gradients indicate Maxent-derived climate suitability for C. truncatus for the: 1981–2010 period (b); 2011–2040 period (c); and 2041–2070 period (d). Stippling delineates the ANUCLIM-derived climate envelope for C. truncatus in each time period. Hatching delineates the current distribution of Quercus in Canada. Climate projections are based on a composite of four climate models and the RCP 4.5 emissions scenario (see text for further details). Maps were generated using ARCGIS v.9.3 (ESRI, Redlands, CA, USA; https://www.esri.com/arcgis/about-arcgis).

Full size image

Figure 3

Occurrence data (a) used for generating climate suitability models for Carpophilus sayi. Maps with colour gradients indicate Maxent-derived climate suitability for C. sayi for the: 1981–2010 period (b); 2011–2040 period (c); and 2041–2070 period (d). Stippling delineates the ANUCLIM-derived climate envelope for C. sayi in each time period. Hatching delineates the current distribution of Quercus in Canada. Climate projections are based on a composite of four climate models and the RCP 4.5 emissions scenario (see text for further details). Maps were generated using ARCGIS v.9.3 (ESRI, Redlands, CA, USA; https://www.esri.com/arcgis/about-arcgis).

Full size image

Climate estimates were obtained at each occurrence location by interrogating North American climate models (described in McKenney et al.19) of the 1981–2010 normal period for the following four variables: (1) fall (i.e., September–December) precipitation (FALLPCP); (2) average spring (i.e., March–June) temperature (SPRINGTMP); (3) annual climate moisture index (CMI; a climate-based moisture balance variable; see Hogg20 for details); and (4) Average Minimum Temperature of the Coldest Month (MINTCM). These climate variables were selected based on their reported influences on B. fagacearum (FALLPCP and CMI21,22) C. sayi and C. truncatus (SPRINGTMP10), and insect distributions in general (MINTCM23). Further, none of the selected variables were highly correlated (i.e., r  60 years) for each province in our study area. Focusing on current/near-term oak timber stocks (since we are not estimating B. fagacearum spread), we multiplied merchantable volume over 40 years old—roughly the age at which oak becomes harvestable in Ontario42—by average provincial stumpage values.
Estimating stumpage values was somewhat challenging due to inter-provincial variation in stumpage systems and reporting of stumpage fees. For the province of Québec, we obtained oak-specific stumpage fees for 191 harvest zones for the period April 1, 2019 to March 31, 2020 (Bureau de mise en marché43). Since stumpage fees in Québec vary by wood quality class (i.e., A, B, and C), we further obtained information on the proportion of wood harvested in each class over the same period (unpublished dataset, Bureau de mise en marché des Bois). We then calculated the average stumpage fee for the province by averaging across harvest zones and quality classes, while weighting by the proportion of wood in each quality class. For Ontario, we obtained stumpage fees for two quality classes of hardwoods (i.e., Class 1 and Class 2) and four oak-related product types (Veneer, Sawlogs, Composite, and Firewood) for January 1 to December 31, 2019 (Ontario Ministry of Natural Resources and Forestry44). Given that oak is typically considered a higher value hardwood, and in lieu of information on how oak is partitioned across product types in Ontario, we calculated oak stumpage fees as an average across product types for the Class 1 hardwood category. Note that the stumpage rates employed here include the Renewal and Forest Futures fees that are part of the Ontario stumpage system. Finally, stumpage values for the province of Nova Scotia were obtained for a single hardwood quality class for the period April 1, 2017 to March 31, 2018 (Province of Nova Scotia45). As in Ontario, stumpage fees were averaged across product types. The stumpage values for Nova Scotia were applied to the relatively small amount of oak in the neighbouring Maritime Provinces of New Brunswick and Prince Edward Island (PEI).
Alternatively, gross domestic product (GDP) can provide an estimate of the total economic activity associated with a given industry. Annual GDP estimates for broad categories (e.g., forestry and logging industry, and wood product manufacturing industry) are available for each province at Natural Resources Canada’s forestry statistics website (https://cfs.nrcan.gc.ca/statsprofile/overview/ca). In order to estimate GDP specifically for oak-related timber products, we first multiplied these provincial broad category GDP values by the proportion of the total provincial harvest that was composed of hardwoods (multipliers obtained from published provincial data sources as detailed in the Results section below). This value was then further refined by multiplying by the proportion of hardwoods in the province that was composed of oak species. These estimates were obtained from forest attribute grids13, by summing merchantable volume of (1) oak and (2) all broadleaf species within the industrial forestry limits of each province. Spatial summaries were carried out using the raster and rgdal packages in r. Though admittedly coarse, we felt this approach was the best available given the dearth of readily available economic data for individual tree species/genera; similar approaches have been used previously to estimate economic impacts of invasive species46,47.
These two approaches (i.e., stumpage-based and GDP-based) provide different perspectives on oak-related timber values at risk. The stumpage approach attaches a basic price to standing timber resources, but does not consider downstream economic activities associated with harvest, such as wages, equipment purchases, and capital expenditures. This approach implicitly assumes that substitution possibilities (e.g., other tree species) can fully replace oak-related contributions to the economy with minimal adjustment costs and, as such, is a conservative estimate of potential timber value losses. Alternatively, the GDP approach attempts to include all downstream economic contributions and assumes little or no opportunity for substitution, such that oak timber losses would be accompanied by a proportional reduction in economic activity. We present both estimates here to provide policy-makers with a range of possible impacts. The value of costs through time is generally arrived at using economic discounting; however, here we have no estimates of the timeline associated with oak wilt spread and hence have chosen to report gross, undiscounted values. See Aukema et al.48 for further discussion. More

1.
BirdLife International. Pygoscelis antarcticus. The IUCN Red List of Threatened Species 2018:e.T22697761A132601557. https://doi.org/10.2305/IUCN.UK.2018-2.RLTS.T22697761A132601557.en (2018).
2.
Sander, M., Balbao, T. C., Polito, M. J., Costa, E. S. & Carneiro, A. P. B. Recent decrease in Chinstrap penguin (Pygoscelis antarctica) populations at two of Admiralty Bay’s islets on King George Island, South Shetland Islands, Antarctica. Polar Biol. 30, 659–661 (2007).
Article  Google Scholar 

3.
Naveen, R., Lynch, H. J., Forrest, S., Mueller, T. & Polito, M. First direct, site-wide penguin survey at Deception Island, Antarctica, suggests significant declines in breeding Chinstrap penguins. Polar Biol. 35(12), 1879–1888 (2012).
Google Scholar 

4.
Lynch, H. J. et al. In stark contrast to widespread declines along the Scotia Arc, a survey of the South Sandwich Islands finds a robust seabird community. Polar Biol. 39, 1–11 (2016).
ADS  Article  Google Scholar 

5.
Dunn, M. et al. Population size and decadal trends of three penguin species nesting at Signy Island South Orkney Islands. PLoS ONE 11, e0164025 (2016).
Article  Google Scholar 

6.
Strycker, N. et al. Fifty-year change in penguin abundance on Elephant Island, South Shetland Islands, Antarctica: Results of the 2019–20 census (2020) (in review).

7.
Croxall, J. P. & Furse, J. R. Food of Chinstrap penguins (Pygoscelis antarctica) and Macaroni penguins (Eudyptes chrysolophus) at Elephant Island Group, South Shetland Islands. Ibis 122, 237–245 (1980).
Article  Google Scholar 

8.
Trivelpiece, W. Z. et al. Variability in krill biomass links harvesting and climate warming to penguin population changes in Antarctica. PNAS 108(18), 7625–7628 (2011).
ADS  CAS  Article  Google Scholar 

9.
Lynch, H. J., Naveen, R., Trathan, P. N. & Fagan, W. F. Spatially integrated assessment reveals widespread changes in penguin populations on the Antarctic Peninsula. Ecology 93(6), 1367–1377 (2012).
Article  Google Scholar 

10.
Casanovas, P., Naveen, R., Forrest, S., Poncet, J. & Lynch, H. J. A comprehensive coastal seabird survey maps out the front lines of ecological change on the Western Antarctic Peninsula. Polar Biol. 38, 927–940 (2015).
Article  Google Scholar 

11.
Croxall, J. P., Trathan, P. & Murphy, E. Environmental change and Antarctic seabird populations. Science 297, 1510–1514 (2002).
ADS  CAS  Article  Google Scholar 

12.
CCAMLR. Standard methods for monitoring parameters of predator species. CCAMLR Ecosystem Monitoring Programme, Hobart (2004).

13.
Hinke, J., Salwicka, K., Trivelpiece, S., Watters, G. & Trivelpiece, W. Divergent responses of Pygoscelis penguins reveal common environmental driver. Oecologia 153, 845–855 (2007).
ADS  Article  Google Scholar 

14.
Lynch, H. J., Fagan, W. F. & Naveen, R. Population trends and reproductive success at a frequently visited penguin colony on the western Antarctic Peninsula. Polar Biol. 33, 493–503 (2010).
Article  Google Scholar 

15.
Fraser, W. R., Trivelpiece, W. Z., Ainley, D. G. & Trivelpiece, S. G. Increases in Antarctic penguin populations: Reduced competition with whales or a loss of sea ice due to environmental warming?. Polar Biol. 11(8), 525–531 (1992).
Article  Google Scholar 

16.
Hinke, J. et al. Identifying risk: concurrent overlap of the Antarctic krill fishery with krill-dependent predators in the Scotia Sea. PLoS ONE 12, e0170132 (2017).
Article  Google Scholar 

17.
Hill, S. L. et al. Reference points for predators will progress ecosystem-based management of fisheries. Fish Fish. 21(2), 368–378 (2020).
Article  Google Scholar 

18.
Juáres, M. A. et al. Adélie penguin population changes at Stranger Point: 19 years of monitoring. Antarct. Sci. 27(5), 455–461 (2015).
ADS  Article  Google Scholar 

19.
Lynch, M., Youngflesh, C., Agha, N. M., Ottinger, M. A. & Lynch, H. Tourism and stress hormone measures in Gentoo Penguins on the Antarctic Peninsula. Polar Biol. 42, 1299–1306 (2019).
Article  Google Scholar 

20.
Croxall, J. P. & Kirkwood, E. D. The Distribution of Penguins on the Antarctic Peninsula and Islands of the Scotia Sea (British Antarctic Survey, Cambridge, 1979).
Google Scholar 

21.
Woehler, E. J. The Distribution and Abundance of Antarctic and Subantarctic Penguins (Scientific Committee on Antarctic Research, Cambridge, 1993).
Google Scholar 

22.
Conroy, J. W. H., White, M. G., Furse, J. R. & Bruce, G. Observations on the breeding biology of the Chinstrap penguin, Pygoscelis Antarctica, at Elephant Island, South Shetland Islands. Br. Antarct. Surv. Bull. 40, 23–32 (1975).
Google Scholar 

23.
Lowther, A. D. Antarctic marine mammals. Encyclopedia of Marine Mammals 3rd edn, 27–32 (Academic Press, London, 2018).
Google Scholar 

24.
Hinke, J. et al. Spatial and isotopic niche partitioning during winter in Chinstrap and Adélie penguins from the South Shetland Islands. Ecosphere 6, 125 (2015).
Article  Google Scholar 

25.
Dimitrijević, D. et al. Isotopic niches of sympatric Gentoo and Chinstrap penguins: evidence of competition for Antarctic krill?. Polar Biol. 41, 1655–1669 (2018).
Article  Google Scholar 

26.
Trathan, P. N., Croxall, J. P. & Murphy, E. J. Dynamics of Antarctic penguin populations in relation to inter-annual variability in sea ice distribution. Polar Biol. 16, 321–330 (1996).
Article  Google Scholar 

27.
He, H. et al. Aerial photography based census of Adélie penguin and its application in CH4 and N20 budget estimation in Victoria Land, Antarctic. Sci. Rep. 7, 12942 (2017).
ADS  Article  Google Scholar 

28.
Bromwich, D. H. et al. Central West Antarctica among the most rapidly warming regions on Earth. Nat. Geosci. 6, 139–145 (2013).
ADS  CAS  Article  Google Scholar 

29.
World Meteorological Organization. 2020.

30.
Flores, H. Impact of climate change on Antarctic krill. Mar. Ecol. Prog. Ser. 458, 1–19 (2012).
ADS  Article  Google Scholar 

31.
Reid, K., Sims, M., White, R. W. & Gillon, K. W. Spatial distribution of predator/prey interactions in the Scotia Sea: Implications for measuring predator/fisheries overlap. Deep Sea Res. Part II 51, 1383–1396 (2004).
ADS  Article  Google Scholar 

32.
Watters, G. M., Hinke, J. & Reiss, C. S. Long-term observations from Antarctica demonstrate that mismatched scales of fisheries management and predator–prey interaction lead to erroneous conclusions about precaution. Sci. Rep. 10, 2314 (2020).
ADS  CAS  Article  Google Scholar 

33.
Polito, M. J. et al. Stable isotope analyses of feather amino acids identify penguin migration strategies at ocean basin scales. Biol. Lett. 13, 20170241 (2017).
Article  Google Scholar 

34.
Hinke, J., Santos, M., Malgorzata Korczak-Abshire, G. M. & Watters, G. Individual variation in migratory movements of Chinstrap penguins leads to widespread occupancy of ice-free winter habitats over the continental shelf and deep ocean basins of the Southern Ocean. PLoS ONE 14(12), e0226207 (2019).
CAS  Article  Google Scholar 

35.
Hurlbert, A. H. & Liang, Z. Spatiotemporal variation in avian migration phenology: citizen science reveals effects of climate change. PLoS ONE 7, e31662 (2012).
ADS  CAS  Article  Google Scholar 

36.
Cohen, J. M., Lajeunesse, J. M. & Rohr, J. R. A global synthesis of animal phenological responses to climate change. Nat. Clim. Change 8, 224–228 (2018).
ADS  Article  Google Scholar 

37.
Samplonius, J. M. et al. Strengthening the evidence base for temperature-mediated phenological asynchrony and its impacts. EcoEvoRxiv Preprints. https://doi.org/10.32942/osf.io/jmy67. Published online ahead of print May 5, 2020.

38.
Robinson, R. A. et al. Travelling through a warming world: climate change and migratory species. Endanger. Species Res. 7, 87–99 (2009).
ADS  Article  Google Scholar 

39.
Atkinson, A. et al. Krill (Euphausia superba) distribution contracts southward during rapid regional warming. Nat. Clim. Change 9, 142–147 (2019).
ADS  Article  Google Scholar 

40.
Freer, J. et al. Limited genetic differentiation among Chinstrap penguin (Pygoscelis antarctica) colonies in the Scotia Arc and Western Antarctic Peninsula. Polar Biol. 38(9), 1493–1502 (2015).
Article  Google Scholar 

41.
Clucas, G. V. et al. Comparative population genomics reveals key barriers to dispersal in Southern Ocean penguins. Mol. Ecol. 27, 4680–4697 (2018).
CAS  Article  Google Scholar 

42.
Humphries, G. et al. Mapping Application for Penguin Populations and Projected Dynamics (MAPPPD): data and tools for dynamic management and decision support. Polar Rec. 53, 1–7 (2017).
Article  Google Scholar 

43.
Lynch, H. J. & LaRue, M. A. First global census of the Adélie penguin. The Auk Ornithol. Adv. 131, 457–466 (2014).
Google Scholar 

44.
Che-Castaldo, C. et al. Pan-Antarctic analysis aggregating spatial estimates of Adélie penguin abundance reveals robust dynamics despite stochastic noise. Nat. Commun. 8, 832 (2017).
ADS  Article  Google Scholar  More

Microorganisms drive biogeochemical cycling in the earth system1 (Fig. 1). Photoautotrophic microorganisms are responsible for about half of CO2 fixation and O2 production on earth, and heterotrophic microorganisms are responsible for much of the return reaction: the oxidation of organic matter back into CO2. The temporal and spatial separation of photoautotrophy and heterotrophy in the global environment drives the biological sequestration of carbon, the reduction of atmospheric CO2, and the maintenance of elevated atmospheric and oceanic O22,3,4,5. Chemoautotrophic microorganisms also fix CO2 and, together with anaerobic heterotrophic metabolisms, carry out diverse chemical transformations including the fluxes of nitrogen to and from biologically available states and the formation of the potent greenhouse gas nitrous oxide (N2O)6,7. Since these transformations respond to, and feedback on, changes in climate (Fig. 1), estimating microbial activity accurately at global scales is important for climate science.
Fig. 1: Key microbially driven redox transformations that mediate the atmospheric fluxes of climatically relevant gases.

Radiatively active gases are notated with red type. The processes in black type are represented in some way (though not necessarily with electron balancing) in both the marine and terrestrial biospheres in earth system models within the Coupled Model Intercomparison Project (land: NCAR Community Earth System Model103; ocean: GFDL COBALTv2104), which are used for projections of climate change in reports by the Intergovernmental Panel on Climate Change. Processes in green type are represented in only the terrestrial model. Current models do not yet include other relevant reactions, some of which are represented in gray type, such as anaerobic ammonia oxidation (anammox), the marine production and consumption of methane, the redox cycling of iron, manganese, and other metals, and the methane-relevant redox chemistry of phosphorus105. COBALTv2 does account for sulfate reduction in marine sediments, but sulfate is not represented. Image courtesy of NASA.

Full size image

However, understanding and projecting the impacts of microbial processes are limited in part due to oversimplified representation in earth system models. For example, in marine biogeochemical models, much attention is given to the complex impacts of phytoplankton—the photoautotrophic microorganisms responsible for primary production—and their small zooplankton predators8,9,10. The bacterial and archaeal activities responsible for other critical aspects of biogeochemical cycling in the land and ocean—remineralization, denitrification, nitrogen fixation, methanogenesis, etc.—are often crudely parameterized10. Such models have limited prognostic capability. For example, models typically prescribe the ecological niche of a given metabolism with imposed, empirically determined parameters that are site- or organism-specific. These parameterizations may or may not apply to other environments, including past and future ecosystems.
These simplistic approaches have been largely necessary due to the difficulties of characterizing the taxonomy and metabolic capabilities of natural microbial communities. However, the rapid expansion of genetic sequencing capabilities has enabled a clearer view of microbial biogeography and activity in the environment. In consequence, computational biogeochemistry is opening up the black box of remineralization and other microbially mediated processes in marine and terrestrial environments11,12,13,14,15,16,17,18.
As we expand models to include the full metabolic potential of microorganisms, how can we organize and reduce the complexity of the descriptions of metabolic diversity? Non-photosynthetic organisms oxidize chemical species for energy, and thus their respiration is biogeochemically significant19. Here, we explain how the key reduction-oxidation (redox) reactions that supply energy for metabolisms can provide an additional organizing principle for explicit descriptions of microbial populations in ecosystem models. This redox basis can be exploited to quantitatively resolve chemical transformations in terms of assimilatory and respiratory fluxes. While not yet incorporated into earth system models, this view has been advocated for such applications20, and has been embraced and employed in the field of environmental biotechnology, such as in the interpretation and modeling of wastewater bioreactors21. Just as models of ocean and atmospheric circulation are constrained by conservation of energy and potential vorticity, complementing mass balance with powerful redox and energetic constraints enables self-consistent descriptions of diverse microbial metabolisms.
This approach aims to advance ecological modeling beyond species-specific descriptions to those that matter for biogeochemical function, in line with trait-based modeling approaches9. In analogy to the use of redox chemistry, trait-based functional type models of phytoplankton have used cell size as an organizing principle for understanding phytoplankton biogeography, biodiversity, and impact on biogeochemistry9,22,23. These types of theoretical constraints allow for the inclusion of more functional types without introducing as many degrees of freedom as would be necessary if each were empirically described. The guiding perspective is that organizing complex biological behavior by its underlying chemical and physical constraints gives more universally applicable descriptions of large-scale biogeochemistry.
When incorporating a redox-balanced approach into ecosystem models, microbial function emerges from underlying chemistry as a consequence of interactions between populations modeled as metabolic functional types and their environment. Resulting theoretically grounded ecosystem models independently simulate microbial growth, respiration, and abundances in ways that we can compare with observations such as sequencing datasets. Thus, sequencing datasets are used as critical tests for the models, as external constraints rather than as input to the models, allowing for an iterative relationship between theory, observations, and models.
In contrast with empirically informed models, this approach involves constructing a model of microbial activity theoretically, and then comparing the results with the observations in order to gain an understanding of the system. The goal is to understand why biology functions as it does, in addition to anticipating global impacts. From a first-principles biogeochemical perspective with respect to physical and chemical forcing, genes are an intermediate step between forcing and function, with the detailed complexity of biological reality following the underlying chemical and physical constraints (analogous to the form follows function principle of architect Louis Sullivan). This does not equate to thinking that biology (or genetic information) does not matter or can be replaced entirely in models by physics and chemistry. Rather, we want to fundamentally understand biological activity as an integrated part of an ecosystem, and physics and chemistry become tools for doing so.
Here, we outline the basis for using redox chemistry as an organizing principle and its translation into quantitative descriptions of microbial activity that are simple enough for global earth system models. We then discuss the benefits of this approach in the context of their implications for improved understanding and projections of global change impacts. Finally, we discuss limitations and possible future developments.
Predicting microbial activity
From one perspective, microbial communities are characterized by interactions at the micro-scale: gene expression, enzymatic capabilities, metabolites, species-specific interdependencies, etc., as well as the physical and chemical environment surrounding small cells24,25,26,27,28. The information from sequencing in particular has allowed for a huge expansion of insight into the detailed in situ activity of uncultivated species. When investigating global-scale impacts, how do we decide which of these details may be bypassed for simplicity? Or, if this simplification is impossible, must we incrementally construct a microbial ecosystem model that incorporates all known micro-scale detail?
Another way forward arises from a macro-scale perspective, which examines how ecosystem function relates to the chemical potential utilized by organisms for energy20,29,30. For example, it is well known that microbial communities in sediments and anoxic zones organize according to the redox tower – the ranking of half-reactions by electrochemical potential14,31,32. Furthermore, respiration by living organisms increases the entropy of the environment by dissipating concentrated sources of chemical energy in accordance with the Second Law of Thermodynamics30,33,34.
This perspective suggests that chemical potential can be used to predict the activity of microbial communities and their biogeochemical impact. However, given the notorious complexity of microbial cells and systems, which is many steps away from governing chemical or physical equations, how can we be sure that this activity is indeed predictable? Frentz et al.35 demonstrated that external conditions cause the seemingly random fluctuations observed in microbial growth, rather than stochastic variation in gene expression. This provides direct evidence of deterministic behavior, and so the authors conclude that microbial systems can in principle be determined by macroscopic laws.
How is this determinism manifested? If microbial communities can respond relatively quickly to changes to their local environment, they may predictably optimize the exploitation of locally available resources. In the ocean, dispersal in microbes is thought to be a highly efficient process such that microbial communities can in effect draw from an extensive seed bank36,37, as captured in the phrase “everything is everywhere, the environment selects”38. Furthermore, recent evidence also shows that gene acquisitions and deletions happen quickly enough to allow for horizontal gene transfer to dominate bacterial adaptation39,40,41,42,43, implying that evolution can occur within few generations and thus on timescales similar to ecological interactions. Perhaps consequentially, similar geochemical environments have been demonstrated to have high-microbial functional redundancy despite different taxonomic compositions17,44. This may be interpreted with the hypothesis that physics and chemistry selects for metabolic traits, and that these traits can be housed in different organisms with taxonomic composition shaped by micro-scale or biotic interactions17,44,45.
The prediction of microbial activity from environmental chemical potential has a long history in microbiology20,21,46,47,48,49,50,51,52, and is conceptually similar to other redox-balanced approaches to understanding microbial activity in sediments, soils, subsurfaces, and aquatic systems13,31,53,54,55,56,57,58. Illustrating the power of these approaches, anticipating metabolism from chemical potential resulted in a prediction that anaerobic ammonia oxidation (anammox) should exist decades before it was observed59,60. Quantitatively understanding microbially mediated rates of conversion of substrates has practical implications for wastewater treatment, and thus the field of biotechnology has established methodologies for an approach in textbook form21. Flux balance analysis (FBA) models can be considered as much more highly detailed analogs of this approach that resolve the mass and electron balances among a multitude of chemical reactions within a cell61,62.
Redox-balanced metabolic functional types
We can resolve microbial activity in global ecosystem models using the underlying redox chemistry of diverse metabolisms as a constraint. One specific way forward is to model distinct metabolisms as populations of metabolic functional types. This systematically quantifies relative rates of substrate consumption, biomass synthesis, and excretions of transformed products associated with each metabolism. Coupled with estimates of substrate uptake, this replaces implicit parameterizations of processes such as organic matter consumption, oxygen depletion, and denitrification with electron-balanced respiratory fluxes of dynamic microbial populations. Box 1 provides a detailed description of this methodology for multi-dimensional models.
A particular set of redox reactions may distinguish a functional type, such as the oxidation of organic matter using oxygen (aerobic heterotrophy), or the oxidation of ammonia or nitrite using oxygen (chemoautotrophic nitrification) as exemplified in Table 1. For each metabolism, an electron-balanced description consists of multiple half-reactions: biomass synthesis, oxidation of an electron donor, and reduction of an electron acceptor21,48. The ratio of anabolism and catabolism can then be represented by the fraction f of electrons fueling cell synthesis vs. respiration for energy, following ref. 21. This provides a yield y (moles biomass synthesized per mole substrate utilized) of each required substrate that reflects two inputs: electron fraction f and the coefficients of the half-reactions (Fig. 2). The interlinked yields reflect the energy supplied by the redox reaction, the energy required for synthesis and other cellular demands, and the inefficiencies of energy conversion. Either f or y for any one of the substrates may be estimated theoretically with Gibbs free energies of reaction21 or with a combination of theoretical and empirical strategies63.
Table 1 Simplified equations describing two exemplary metabolic functional types.
Full size table

Fig. 2: Schematic of a single cell represented as a metabolic functional type carrying out the aerobic oxidation of ammonia.

The redox balance informs the elemental ratios of substrates utilized, biomass synthesized, and waste products excreted (Table 1).

Full size image

The result is a stoichiometric budget of the metabolism of the whole organism (Table 1). These descriptions quantify the elemental ratios of utilized substrates, biomass, and the excretion of waste products. For example, the descriptions account for the CO2 produced by heterotrophic metabolisms as well as the CO2 fixed by chemoautotrophic metabolisms (Table 1, Fig. 2, and Supplementary Fig. 2), linking microbial activity directly to global carbon cycling.
To estimate the growth rate of each functional type, the yields from the metabolic budgets are combined with the uptake rates of the required substrates (Box 1). Limiting uptake rates may rely on empirically derived uptake kinetic parameters, or they can be estimated theoretically from diffusive supply, cell size, membrane physiology, and other physical constraints64,65,66. If theoretical models of uptake are used, the physical constraints on substrate acquisition and the redox chemical constraints on energy acquisition can provide an entirely theoretical estimate of the growth of each metabolic functional type.
One strategy is to represent the populations carrying out each of these discrete metabolisms as one functional type population, which aggregates the diverse community of many species that are fueled by the same (or a similar) redox reaction (Fig. 3). Such aggregation has been deemed a useful strategy for representing the biogeochemical impacts of microbial communities for certain research questions67,68. However, for other questions this wipes out critical diversity among the aggregated populations. For example, diverse aerobic heterotrophic populations consume organic matter over a wide range of rates, and these rates dictate the amount of biologically sequestered carbon in the ocean. Redox chemistry and physical limitations alone may not inform the heterogeneity among similar metabolisms. One additional constraint is the limited capacity of the cell and thus its allocation of proteome towards different functions69. While the electrons supplied to the cell must be conserved following the redox balance, the electrons may be partitioned differently into machinery for substrate uptake vs. biomass synthesis, for instance, for different phenotypes. This partitioning can be quantitatively related to ecological fitness and biogeochemical impact via uptake kinetics, effective yields, and other traits9 (Supplementary Note 1 and Supplementary Fig. 1).
Fig. 3: Solutions from a global simulation resolving multiple metabolic functional types.

Net primary productivity (NPP), the biomasses of the metabolic functional types, and the sinking particulate organic carbon (POC) flux are resolved along a transect of a global microbial ecosystem model coupled with an estimate of the ocean circulation (Darwin-MITgcm18).

Full size image

Box 1 Incorporating metabolic functional types into ecosystem models

A metabolic functional type can be represented as a population with a growth rate that is limited or co-limited by multiple required substrates. If Liebig’s Law of the Minimum is employed, the limiting growth rate μ is described as

$$mu = min (V_{mathrm{i}}y_{mathrm{i}})$$
(1)

where Vi is the specific uptake rate of substrate i, and yield yi is the biomass yield with respect to that substrate. Yields for the different substrates and elements are interlinked in the metabolic budget derived from the underlying redox chemistry. Yields reflect Gibbs free energies of reaction among other factors. In the simplest model, non-limiting substrates are consumed in proportion to the limiting resource according to the metabolic budget, although in reality they may accumulate in the form of storage molecules.
Each metabolic functional type population can be incorporated into a multi-dimensional environmental model (e.g., an ocean simulation) with physical transport as

$$frac{{dB}}{{dt}} = mu B – Lleft( B right)B – underbrace {nabla cdot left( {bf{u}}{B} right)}_{{mathrm{advection}}} + underbrace {nabla cdot left( {{mathbf{kappa }}nabla B} right)}_{{mathrm{diffusion}}}$$
(2)

for biomass concentration B, loss rate L, velocity u, and diffusion coefficient κ. The loss rate function varies with biomass and represents a combination of processes, including predation, viral lysis, maintenance, and senescence. These processes remain largely unconstrained, although efforts have been made to relate losses to ecological dynamics107,108.
The yield partitions the amount of substrate taken up by the population into that used for growth, Viyi, versus that exiting the cell in modified form as a waste product, Vi(1 − yi) (Fig. 2 and Table 1). Equation 1 suggests a correlation between μ and y, but yields may be further modified by other factors. For example, accounting for maintenance energy decreases the ratio of growth to respiration, contributing to a decoupling between growth rate and yield particularly at low growth rates109. Furthermore, a trade-off between uptake rate and yield at the cellular level reflects the allocation of enzyme towards machinery for substrate uptake vs. biomass synthesis, among other factors. Considering a proteome constraint can incorporate this trade-off (Supplementary Note 1).

Benefits and implications for anticipating global change
Redox chemistry aids in reducing the number of degrees of freedom in descriptions of diverse microbial metabolisms. We next discuss the benefits of this electron-balanced approach, each contextualized by specific projected impacts of global change due to microbial activity and broad challenges in the fields of microbial ecology and biogeochemistry.
Flexible and broadly applicable metabolic thresholds: A key question for microbial biogeochemical studies, for which biogeochemical models are primed to answer, is how the biogeographies of diverse, active metabolisms vary with changes in the physical and chemical environment. What threshold determines the viability of a given metabolism?
Redox-balanced metabolic budgets obviate the need to impose critical concentrations or other thresholds that determine the presence of any given metabolism. Rather than being imposed following empirical relationships, metabolic biogeography emerges dynamically from ecological interactions and reflects environmental chemical potential. This flexibility aids in understanding metabolic thresholds more fundamentally, and it expands model applicability to diverse and unobserved environments. This is of particular importance for understanding global change, since past and future worlds may include very different ecosystems that do not reflect current empirical trends.
For example, the oceans are currently losing oxygen due to global warming5,70. If anoxic zones continue to expand, this will increase the habitat of anaerobic microorganisms, whose respiration results in emissions of N2 and N2O to the atmosphere7. Many biogeochemical models prescribe O2 concentrations that inhibit anaerobic activity in accordance with observations of specific organisms or communities in experimental conditions. This assumes that the same O2 concentrations limit metabolism similarly in all environments, and often trades mechanistic understanding of oxygen limitation for empirical correlations that may reflect a variety of natural and introduced biases, such as micro-scale heterogeneity, physical mixing in the ocean, and experimental bottle effects.
In contrast, a metabolic functional type model does not require imposed oxygen threshold concentrations (Supplementary Fig. 3). When oxygen supply is abundant, anaerobic types are competitively excluded because growth using alternative electron acceptors is lower than with oxygen. When oxygen supply is low, aerobic populations may persist and continue to deplete any available oxygen even as their growth is limited by oxygen, allowing for a steady state stable coexistence of aerobic and anaerobic metabolisms, which is consistent with a variety of observations71.
Descriptions of microbial growth that reflect underlying chemical potential can enable predictions of many other metabolic transitions, such as nitrogen fixation, nitrification, and the transition to sulfur oxidation and reduction13,18,72,73,74. As another example, this approach predicts the restriction of nitrification from the sunlit surface ocean as a consequence of competitive exclusion by phytoplankton in many environments (Fig. 4 and Supplementary Fig. 2), as well as active nitrification in some surface locations where phytoplankton are limited by another factor not affecting the chemoautotrophs, such as at high latitudes where phytoplankton are limited by light18. The emergent exclusion from most of the surface ocean anticipates that many clades of nitrifying microorganisms have adapted to long-term exclusion from the surface and consequentially lost (or did not develop) photoprotective cellular machinery.
Fig. 4: Model simulation and observations of the marine nitrification system.

Biogeochemistry is driven by microbial metabolic functional types in a vertical water column model18. Lines are model solutions, and marked points are observations from two stations in the Pacific Ocean75,106 (see Supplementary Fig. 2 for more detail) (a). Chlorophyll a concentrations and abundances of ammonia-oxidizing organisms (AOO) and nitrite-oxidizing organisms (NOO). Observed abundances are of the 16S rRNA abundances of archaeal Marine Group I and Nitrospina-like bacteria75,106. Model abundances are converted from biomass with 0.1 fmol N cell-1 for AOO, 0.2 fmol N cell-1 for NOO76, and one gene copy per cell. (b). Light (solar irradiance) and bulk nitrification rates.

Full size image

Replacing implicit descriptions of organic matter remineralization: The fate of organic matter dictates the amount of carbon sequestered in the marine and terrestrial biospheres. Microbial consumption mediates the carbon stored in soils, the carbon stored in the ocean as dissolved organic matter (DOM), and the sinking flux of organic carbon that constitutes the marine biological carbon pump4, without which atmospheric CO2 would be 100–200 ppm higher than current levels. We want to understand how these carbon reservoirs respond to changes in climate, such as increased temperatures and changes in precipitation patterns. However, in biogeochemical models, simple rate constants often dictate the remineralization of elements from organic back into inorganic constituents.
Replacing simplistic parameterizations with dynamic metabolic functional types means that electron-balanced descriptions of growth and respiration instead drive the fate of organic matter in earth system models (Fig. 3). In addition to a more sophisticated and responsive description of carbon sequestration, non-living organic matter is fully integrated into ecosystem frameworks, enabling theoretical studies of phytoplankton-bacteria interactions to complement observational and experimental approaches.
Much work remains in the development of these descriptions. As we discuss below, accurate estimates of organic matter turnover rates require more accurate descriptions of the complex processes governing microbial uptake rates of organic matter. However, the redox-informed yields are still useful for quantifying the relative amount of CO2 excreted and the absolute amount of biomass sustained on a given substrate, independent of uptake kinetics (Supplementary Note 2).
Relationships between abundances, rates, nutrient concentrations, and elemental ratios: An overarching puzzle challenging microbial ecology is to understand how chemical transformations in the environment are set by the ecological interactions at the organism level, among individual microscopic cells. It is clear that abundances of populations are not simply and directly correlated with biogeochemical impact (i.e., higher abundance does not necessarily imply an associated higher rate of chemical transformation). Untangling the relationship between abundances and biogeochemical function is also necessary for interpretation of genetic evidence that provide insight into this complex ecosystem structure.
Redox-balanced metabolic functional type modeling links rates of biomass synthesis associated with a particular metabolism to its rate of respiration as well as the standing stock of limiting nutrients. As functional type modeling is coupled with estimates of population loss rates due to grazing, viral lysis, or other mortality, simulations also resolve the standing stocks of functional biomass. This quantifies the relationship between biomass concentrations and volumetric rates of chemical transformations, emphasizing how relatively low biomass may be associated with relatively high bulk rates71.
For example, the approach has revealed a clear example of the signature of chemical potential in the ecology of marine nitrification18 (Fig. 4). In this model, the two steps of nitrification are represented by two functional type populations. This predicts about a three-fold difference in the abundances of the organisms responsible for each of the two steps of nitrification, despite the fact that the two populations carry out the same rate of subsurface N-cycling at steady state6,18. A three-fold or greater difference in abundance and associated ammonium (NH4+) and nitrite (NO2−) concentrations is consistent with observed differences18,75, and it reflects that the oxidation of one mole of NH4+ generates three times more electrons than the oxidation of one mole of NO2−, with differences in cell size further contributing to differences in abundances (Fig. 4 and Supplementary Fig. 2). Recent observations confirm the redox-based difference in NH4+ and NO2− biomass yield76,77, although measured rates from a nonsteady environment suggest that NO2−-oxidizing bacteria can partition electrons more efficiently than NH4+-oxidizing archaea76 (i.e., higher fraction f despite lower yield y; see Supplementary Note 3).
As redox-based descriptions resolve the stoichiometry of whole organism metabolism, they also link together elemental cycles. Explicit description of relative elemental flow through the ecosystem, and specifically their variation from average values, is critical for understanding climate-biogeochemical feedbacks78,79,80. For example, the nitrification model also estimates the CO2 fixation rates associated with nitrification rates (Supplementary Fig. 2), enabling global-scale, electron-balanced projections of the amount of carbon converted to organic form by chemoautotrophic nitrifying microorganisms.
Connections with sequencing datasets: How do we relate metabolic functional type models to sequencing datasets measuring genetic, transcriptomic, and proteomic diversity? Connecting biogeochemical models with sequencing data is critical because this data provides an enormous amount of information about ecosystem structure and function. Genes (or transcripts) themselves are not necessarily the most concise or useful currency given functional redundancies, unattributed function, and variation in gene dosage from horizontal gene transfer as well as growth rate40. Recent gene-centric models aim to resolve the abundances of key genes as proxies for a predetermined set of metabolic pathways13,14,17. However, the parameters used to describe metabolic pathways in these models are estimated similarly to the redox-balanced yields and efficiencies described here.
The innovation of gene-centric models is the sophisticated conversion of estimates of biogeochemical activity and biomass to genes. For example, the model of Coles et al.17 resolves biomass and nutrient concentrations prognostically, and then uses a three-part formula—representing constitutive, regulated, and steady state transcription—to diagnostically calculate transcription rates from modeled biomass and growth rates17. Thus, the two types of modeling are complimentary, with redox chemistry providing estimates of metabolic activity from fundamental principles, and the careful calibrations between activity and sequencing providing a comparative metric.
The examples here externalize the conversion between modeled activity and sequencing information as a transparent process. In Fig. 4, the predicted functional biomass of ammonia-oxidizing population is related to archaeal Marine Group I (MGI) and Nitrospina-like 16S rRNA genes with two conversion factors: the cell elemental quota (fmol N cell−1) and the number of cellular gene copies. Conversion error arises since cell mass and size vary with growth rate81,82. Maintaining transparency of the conversion from predicted microbial activity to genes and transcripts allows interdisciplinary audiences to understand and critique the models.
Limitations and possible extensions
Using chemical potential as a theoretically grounding organizing principle for the resolution of diverse metabolisms can greatly improve microbial descriptions in global biogeochemical models. However, the approach does have its limitations, which generally increase in significance with increased temporal or spatial resolution.
Modeling metabolic diversity with functional type populations requires choosing how metabolisms are distributed among the populations. This has consequences when interpreting time-varying states: model solutions become dependent on the partitioning of metabolism among the functional types as the timescales of physical change approach the timescales of microbial growth (see Supplementary Note 4, Supplementary Fig. 3, and Supplementary Fig. 4 for a detailed example). Other species-specific time-varying phenomena such as the lag response of organisms to substrate availability also become relevant83. On one hand, this is beneficial for resolution of microbial processes in fine-grained ocean circulation models where flow can vary on the order of days. However, incorporating another constraint, such as proteome allocation69, is necessary to inform these choices. For example, considering enzymatic allocation in combination with energetics allowed for the prediction of both the division of nitrification into a two-step process in mixed environments and the combined, complete pathway in one organism (comammox) in biofilms, which preceded observations of the latter84,85,86.
Uncertainty in distributions of metabolism lies not only in the length of a metabolic pathway, but also in the degree of metabolic versatility (metabolic mixotrophy). Such versatility characterizes key players in large-scale biogeochemistry, such as nitrite-oxidizing bacteria and photoheterotrophs87,88,89. Mixotrophic lifestyles can increase the fitness of populations in their environments, impacting overall ecosystem function90. In one sense, the approach here provides a prediction of where we might expect such mixotrophy by resolving stable coexistences of diverse metabolisms. In Fig. 3, for example, syntrophic coexistence occurs at depth among heterotrophs, ammonia oxidizers, and nitrite oxidizers, and future work could investigate what determines which combinations of these coexistences remain as passive interactions, which develop into mutualistic dependencies as active interactions91, and which evolve into mixotrophic phenotypes or endosymbionts. Additionally, by considering the potential to carry out a metabolism as a trait, we can use the current framework along with an additional constraint to investigate implications of metabolic mixotrophy. For example, Coles et al.17 impose a trade-off between the degree of metabolic diversity of a single functional type and growth rate, enabling the exploration the consequences of distribution of metabolism on the biogeochemical state.
Also, the metabolic functional type approach resolves only active functional biomass, while evidence suggests that less than 10% to more than 75% of the microbial community may be inactive92. Some seemingly inactive populations may slowly metabolize over long timescales, requiring longer model integration times and careful attention to their loss rates for resolution, while some populations are periodically active as revealed by high-resolution observations in time37.
The proposed modeling approach relies on estimates of the limiting uptake rates of required substrates. In lieu of suitable theoretical descriptions, the use of empirically derived uptake kinetic parameters still employs the benefits of the redox-informed yields (Supplementary Note 2). However, underlying physical constraints to substrate acquisition can in principle be exploited to develop more universally applicable descriptions for a variety of substrates and contexts. Uptake kinetics are complex, but for many limiting resources, encounter effectively controls the uptake, and the physics of encounter has been relatively well described. For example, uptake rates estimated from diffusive supply of substrate, cellular geometry, and membrane physiology64,65,66 have been empirically supported93. For organic matter, future work is needed to develop suitable descriptions of consumption rates, whether empirical or theoretical. For example, descriptions require attention to the hydrolysis of organic compounds by extracellular enzymes and the ecology of sinking marine particles—the diffusion of monomer away from the particle, within-particle transport, and dynamic ecological interactions on particle surfaces, among other processes11,74,94,95.
In Fig. 3, the electron-balanced description consists of an average stoichiometry and electron fraction for one sinking pool of organic matter in the ocean, which, as mentioned above, is not sufficient to accurately resolve the carbon storage that is shaped by a distribution of turnover rates. As one of the many factors impacting the rates, an energetics-based perspective can serve as a tool for further deciphering organic matter complexity. For example, organic matter may be partially organized by the nominal oxidation state of its carbon atoms, which relates to a measure of free energy and accessibility96,97. This could be used to improve the phenomenological description of organic matter in models as labile vs. non-labile, for example, with a more mechanistic underpinning.
Descriptions of phytoplankton are currently much more sophisticated than of bacteria and archaea in models, reflecting a longer history of comprehensive sets of observations. However, further work could develop simple descriptions of photoautotrophic metabolisms from underlying energetics by connecting the supply of photons to available energy for biosynthesis within the cell. Many biogeochemical models account for an inefficiency of phytoplankton metabolism with a parameter that dictates their excretions of dissolved organic matter98. Incorporating this excretion into an energetic framework would enhance studies of phytoplankton ecology, such as studies of photoautotrophic-heterotrophic interactions in the ocean surface or photoautotrophic-chemoautotrophic interactions at the base of the euphotic zone where some phytoplankton excrete nitrite due to incomplete reduction of nitrate99.
As a more radical extension, can we progress past population modeling and model microbial consortia as one aggregate community biomass34,100? This may improve resolution of time-varying metabolic versatility. However, if both steps of nitrification were a part of such a consortium, would the characteristic accumulation of nitrite be predicted (Supplementary Fig. 2)? We leave these questions for future research and conclude that the best choice for the degree of resolution of metabolism will depend on the specific research question and the available observations.
We have described a useful approach for understanding and anticipating microbial control of biogeochemical cycling that is suitable for global applications. The approach aims to represent microbial growth and respiration explicitly and consistently from knowledge of chemical gradients in the environment, towards a goal of building an independently constructed theoretical ecosystem model that can then be compared to observations. Describing microbial communities with underlying energetic constraints connects metabolisms dynamically with global geochemical distributions, such as those of carbon dioxide, oxygen, and biologically available nitrogen. This deepens our understanding of microbial ecosystems and enables the incorporation of the feedbacks of microbial activity to changes in global biogeochemistry and the climate system. More

Experimental limitations and implications
Many factors are known to contribute to mollusk shell strength including thickness, microstructure, previous shell damage, and ornamentation19,21,29,30,31. Thus, it is challenging to separate the influence of one feature from another when testing the strength of real shells. Fortunately, there are alternative methods by which to model shells22,32,33,34. 3D printing is an effective tool to create model shells that exhibit brittle behavior in compression, serving as an accurate first-order approximation for shell breakage (Fig. 2). 3D printed models are not created to replicate shell microstructure; rather they serve to normalize confounding factors encountered with real shells, like size, variations in thickness, and taphonomy/degradation of microstructure, to isolate variables of interest while maintaining the predominantly brittle behavior seen in real shells. As a result, the failure of different 3D printed shapes under bulk mechanical compression can be used for relative comparisons between morphologies without needing to reproduce the exact magnitude of load to failure of real shells22,33 (see “Methods”, Fig. 2). Natural shells ultimately break after cracking through microstructural layers. Analogously, in this study, 3D printed shells cracked predominantly through layers of printed materials, not along printing boundaries, demonstrating an appropriate proxy for shell failure (Figs. 1, 2).
Resistance to compression by a predator can be tested directly using a variety of loading tests, the effectiveness of which has been demonstrated repeatedly by investigators studying topics ranging from predation to climate change22,33,34. Compression experiments of bivalve shells, specifically, have been used extensively for this purpose35,36,37,38,39,40,41,42,43,44,45,46. In this study, the mechanical strength of model bivalve shell shapes was analyzed to understand potential defensive value against durophagous (crushing) predators using hypothetical flat teeth and jaws to crush prey (see Crofts and Summers (2014)22 for examples of flat crushing morphologies) (Figs. 1, 2). These experiments are most analogous to small shell-crushers with flat dentitions, like fishes (e.g., guitar fishes, stingrays, etc.), rather than large predators like the modern walrus or extinct marine reptiles (e.g., placodonts, mosasaurs) which are so much larger than their prey that the differences in shell strength resulting from shape are likely insignificant. Bivalves are also preyed upon by many other predators including asteroids, gastropods, birds, and mammals47. The ability to escape is likely dependent upon predator capabilities, where not all escape mechanisms are equally effective against different predators. The results of these experiments are not intended to be general proxies of predation resistance—they are only applicable as a proxy for predators that use flat crushing dentitions. Additionally, it is important to note that compression of a shell by vertebrate predators is a different mechanical process from compression by an arthropod, as claws localize forces differently than teeth and jaws48. Therefore, the experiments used in this study are not meant to model predation by invertebrate durophages. This experimental setup represents an idealized case of shell compression assuming consistent shell thickness and a common predator using a quasistatic loading regime to isolate the influence of shell shape. Therefore, this study makes no conclusions as to the effects of shell shape on strength under conditions of point loading by claws or impact. Furthermore, while bivalve shape is undoubtedly influenced by many factors including fabrication (shell growth and construction)49, phylogenetics, location of soft tissue, etc., these experiments were designed to study an idealized case in which defense against vertebrate shell crushing predators is the most important functional constraint on shell shape.
Advantageously, the use of mathematically generated theoretical bivalve shells, rather than real shells, enables testing a range of shapes—some of which can be or have been found in nature—and others that have yet to exist. For example, these methods enable physical testing of morphologies which are only found in the fossil record. Fossils cannot be used for mechanical experimentation due to changes in the integrity of the shell structure resulting from the fossilization process. Furthermore, testing shapes which do not exist due to biological constraints like fabrication, can provide valuable insights when combined with the study of shapes that have evolved naturally. For example, testing theoretical shapes can reveal morphologies that perform better than natural morphologies in specific functional settings. Thus, testing shapes which perform well, but do not exist, is one way to identify the potential influence of evolutionary constraints such as fabrication (shell growth patterns), phylogenetics, or other necessary biological functions. In this study, the extremely perpendicular-elongate shells represent theoretical shapes that are not seen in nature (likely due to space limitations for soft tissue). Developing theoretical physical models also allows for the elimination of many confounding variables that cannot be avoided when testing real shells. The effects of variations in shell thickness and ornamentation, which also play a role in shell strength, can thus be separated from the strength imparted by gross shell shape when theoretical models are generated. However, because overall shell strength is derived from a combination of these many factors (microstructure, thickness, etc.31.) the conclusions from these experiments represent a first-order approximation of the effects of shape on shell strength alone. Because the models for this study were generated mathematically, and are therefore not exact replicas of specific taxa, the discussion below aims to address potential tradeoffs based on general shell shapes that could be further studied for specific bivalve taxa in future experiments.
Shell shape and strength
Three parameters of shell shape were modified for this study. (1) Generating curve shape is both a modeling parameter and biological feature of mollusk shells27. In mathematical models, generating curve shape describes the shape that is rotated around an axis to create a surface that represents the overall shell shape. Biologically, the generating curve refers to the portion of the shell upon which new material is secreted by the mantle. In bivalves, this is the shape of the commissure. Thus, by changing the shape of the generating curve in a model bivalve, i.e., changing the shape of the ellipse that generates the surface of the shell, models can be perpendicular-elongated (a  > b) or parallel-elongated (a  More

1.
Turner, M. G. Landscape ecology: The effects of pattern on process. Annu. Rev. Ecol. Syst. 20, 171–197 (1989).
Article  Google Scholar 
2.
Noss, R. F. Wildlife corridors. in Ecology of Greenways (eds. Smith, D. & Hellmund, P.) 43–98 (University of Minesota Press, Minesota, 1993).

3.
Brooks, T. M. et al. Habitat loss and extinction in the hotspots of biodiversity. Conserv. Biol. 16, 909–923 (2002).
Article  Google Scholar 

4.
Hanski, I. The Shrinking World: Ecological Consequences of Habitat Loss, Vol. 14 (International Ecology Institute, Philadelphia, 2005).
Google Scholar 

5.
IUCN. The World Conservation Strategy. (IUCN, UNEP, 1980).

6.
IUCN. The IUCN Red List of Threatened Species. Version 2014.3. (2014).

7.
Dirnböck, T., Dullinger, S. & Grabherr, G. A regional impact assessment of climate and land-use change on alpine vegetation. J. Biogeogr. 30, 401–417 (2003).
Article  Google Scholar 

8.
Gonçalves, J., Honrado, J. P., Vicente, J. R. & Civantos, E. A model-based framework for assessing the vulnerability of low dispersal vertebrates to landscape fragmentation under environmental change. Ecol. Complex. 28, 174–186 (2016).
Article  Google Scholar 

9.
Saunders, D. A., Hobbs, R. J. & Margules, C. R. Biological consequences of ecosystem fragmentation: A review. Conserv. Biol. 5, 18–32 (1991).
Article  Google Scholar 

10.
Fahrig, L. & Merriam, G. Conservation of fragmented populations. Conserv. Biol. 8, 50–59 (1994).
Article  Google Scholar 

11.
Wiens, J. A. Habitat fragmentation: Island v landscape perspectives on bird conservation. Ibis 137, S97–S104 (1994).
Article  Google Scholar 

12.
Diamond, J. M. ‘Normal’ extinctions of isolated populations. In extinctions (ed. Nitecki, M. H.) 191–246 (University of Chicago Press, Chicago, 1984).
Google Scholar 

13.
Laurance, W. F. Comparative responses of five arboreal marsupials to tropical forest fragmentation. J. Mammal. 71, 641–653 (1990).
Article  Google Scholar 

14.
Bennett, A. F. Linkages in the Landscape: The Role of Corridors and Connectivity in Wildlife Conservation. (IUCN, 2003).

15.
Opdam, P. Metapopulation theory and habitat fragmentation: a review of holarctic breeding bird studies. Landsc. Ecol. 5, 93–106 (1991).
Article  Google Scholar 

16.
Thomas, C. D. & Jones, T. M. Partial recovery of a skipper butterfly (Hesperia comma) from population refuges: Lessons for conservation in a fragmented landscape. J. Anim. Ecol. 62, 472–481 (1993).
Article  Google Scholar 

17.
Haddad, N. M. et al. Corridor use by diverse taxa. Ecology 84, 609–615 (2003).
Article  Google Scholar 

18.
Grab, H. et al. Habitat enhancements rescue bee body size from the negative effects of landscape simplification. J. Appl. Ecol. 56, 2144–2154 (2019).
Article  Google Scholar 

19.
Smeraldo, S. et al. Modelling risks posed by wind turbines and power lines to soaring birds: The black stork (Ciconia nigra) in Italy as a case study. Biodivers. Conserv. 29, 1959–1976 (2020).
Article  Google Scholar 

20.
Noss, R. F. A regional landscape approach to maintain diversity. Bioscience 33, 700–706 (1983).
Article  Google Scholar 

21.
Noss, R. F. & Harris, L. D. Nodes, networks and MUMS: Preserving diversity at all scales. Environ. Manag. 10, 299–309 (1986).
ADS  Article  Google Scholar 

22.
Grumbine, R. E. What is ecosystem management?. Conserv. Biol. 8, 27–38 (1994).
Article  Google Scholar 

23.
Forman, R. T. T. Land Mosaics. The Ecology of Landscapes and Regions ( Cambridge University Press, Cambridge, 1995).
Google Scholar 

24.
Jongman, R. H. G. Nature conservation planning in Europe: Developing ecological networks. Landsc. Urban Plan. 32, 169–183 (1995).
Article  Google Scholar 

25.
Kubeš, J. Biocentres and corridors in a cultural landscape. A critical assessment of the ‘territorial system of ecological stability’. Landsc. Urban Plan. 35, 231–240 (1996).
Article  Google Scholar 

26.
Diamond, J. M. The island dilemma: lessons of modern biogeographic studies for the design of natural reserves. Biol. Cons. 7, 129–146 (1975).
Article  Google Scholar 

27.
Wilson, E. O. & Willis, E. O. Applied biogeography. In Ecology and Evolution of Communities (eds Cody, M. L. & Diamond, J. M.) 522–534 (Belknap Press, New York, 1975).
Google Scholar 

28.
Soulé, M. E. Land use planning and wildlife maintainance: Guidelines for conserving wildlife in an urban landscape. J. Am. Plan. Assoc. 3, 313–323 (1991).
Article  Google Scholar 

29.
Opdam, P., Van Apeldoorn, R., Schotman, A. & Kalkhoven, J. Population responses to landscape fragmentation. In Landscape Ecology of A Stressed Environment (eds Vos, C. C. & Opdam, P.) 147–171 (Chapman and Hall, London, 1993).
Google Scholar 

30.
Beier, P. & Noss, R. F. Do habitat corridors provide connectivity?. Conserv. Biol. 12, 1241–1252 (1998).
Article  Google Scholar 

31.
Brown, J. H. & Kodric-Brown, A. Turnover rates in insular biogeography: Effect of immigration on extinction. Ecology 58, 445–449 (1977).
Article  Google Scholar 

32.
Barrett, G. W. & Bohlen, P. J. Landscape Ecology Landscape Linkages and Biodiversity (Island Press, New York, 1991).
Google Scholar 

33.
Forman, R. T. T. & Godron, M. Landscape Ecology (Wiley, New York, 1986).
Google Scholar 

34.
Gilbert-Norton, L., Wilson, R., Stevens, J. R. & Beard, K. H. A meta-analytic review of corridor effectiveness. Conserv. Biol. 24, 660–668 (2010).
Article  Google Scholar 

35.
Mech, S. G. & Hallett, J. G. Evaluating the effectiveness of corridors: A genetic approach. Conserv. Biol. 15, 467–474 (2001).
Article  Google Scholar 

36.
Harris, L. D. & Scheck, J. From implications to applications: the dispersal corridor principle applied to the conservation of biological diversity. in Nature Conservation 2: The Role of Corridors (eds. Saunders, D. A. & Hobbs, R. J.) 189–220 (Surrey Beatty & Sons, 1991).

37.
Hobbs, R. J. & Hopkins, A. J. M. The role of conservation corridors in a changing climate. In The Role of Corridors (eds Saunders, D. A. & Hobbs, R. J.) 281–290 (Surrey Beaty & Sons, New York, 1991).
Google Scholar 

38.
McLaughlin, J. F., Hellmann, J. J., Boggs, C. L. & Ehrlich, P. R. Climate change hastens population extinctions. Proc. Natl. Acad. Sci. 99, 6070–6074 (2002).
ADS  CAS  Article  Google Scholar 

39.
Bennett, G. & Mulongoy, K. J. Review of Experience with Ecological Networks, Corridors and Buffer Zones (Secretariat of the Convention on Biological Diversity, 2006).

40.
MacArthur, R. H. & Wilson, E. O. An equilibrium theory of insular zoogeography. Evolution 17, 373–387 (1963).
Article  Google Scholar 

41.
MacArthur, R. H. & Wilson, E. O. The Theory of Island Biogeography (Princeton University Press, Princeton, 1967).
Google Scholar 

42.
Hanski, I. & Gilpin, M. Metapopulation dynamics: Brief history and conceptual domain. Biol. J. Lin. Soc. 42, 3–16 (1991).
Article  Google Scholar 

43.
Bosso, L., Mucedda, M., Fichera, G., Kiefer, A. & Russo, D. A gap analysis for threatened bat populations on Sardinia. Hystrix Ital. J. Mammal. 27, 11788. https://doi.org/10.4404/hystrix-27.2-11788 (2016).
Article  Google Scholar 

44.
Deus, E. et al. Current and future conflicts between eucalypt plantations and high biodiversity areas in the Iberian Peninsula. J. Nat. Conserv. 45, 107–117 (2018).
Article  Google Scholar 

45.
Johovic, I., Gama, M., Banha, F., Tricarico, E. & Anastácio, P. M. A potential threat to amphibians in the European Natura 2000 network: Forecasting the distribution of the American bullfrog Lithobates catesbeianus. Biol. Conserv. 245, 108–551. https://doi.org/10.1016/j.biocon.2020.108551 (2020).
Article  Google Scholar 

46.
van der Sluis, T. et al. How much Biodiversity is in Natura 2000? The “Umbrella Effect” of the European Natura 2000 protected area network. 147 (Alterra, Wageningen, 2016).

47.
Natura 2000 https://ec.europa.eu/environment/nature/natura2000/index_en.htm (2019).

48.
European Comission. Green Infrastructure (GI)—Enhancing Europe’s Natural Capital. 11 (Brussels, 2013).

49.
European Comission. Technical information on Green Infrastructure (GI). 24 (Brussels, 2013).

50.
Araújo, M. B., Thuiller, W. & Pearson, R. G. Climate warming and the decline of amphibians and reptiles in Europe. J. Biogeogr. 33, 1712–1728 (2006).
Article  Google Scholar 

51.
Bennett, A. F. & Saunders, D. A. Habitat fragmentation and landscape change. In Conservation Biology for All (eds Sodhi, N. S. & Ehrlich, P. R.) 88–106 (Oxford University Press, Oxford, 2010).
Google Scholar 

52.
Cushman, S. A. Effects of habitat loss and fragmentation on amphibians: A review and prospectus. Biol. Cons. 128, 231–240 (2006).
Article  Google Scholar 

53.
Keinath, D. A. et al. A global analysis of traits predicting species sensitivity to habitat fragmentation. Glob. Ecol. Biogeogr. 26, 115–127 (2017).
Article  Google Scholar 

54.
Blaustein, A. R. et al. Amphibian breeding and climate change. Conserv. Biol. 15, 1804–1809 (2001).
Article  Google Scholar 

55.
Gibbons, W. J. et al. The Global Decline of Reptiles Déjà Vu Amphibians. Bioscience 50, 653–666 (2000).
Article  Google Scholar 

56.
Rivera-Ortiz, F. A., Aguilar, R., Arizmendi, M. D. C., Quesada-Avendaño, M. & Oyama, K. Habitat fragmentation and genetic variability of tetrapod populations. Anim. Conserv. 18, 249–258 (2015).
Article  Google Scholar 

57.
Andrews, K. M., Gibbons, J. W. & Jochimsen, D. M. Ecological effects of roads on amphibians and reptiles: A literature review. In Urban Herpetology (eds Mitchell, J. C. et al.) 121–143 (Society for the Study of Amphibians & Reptiles, London, 2008).
Google Scholar 

58.
Hansen, N. A., Sato, C. F., Michael, D. R., Lindenmayer, D. B. & Driscoll, D. A. Predation risk for reptiles is highest at remnant edges in agricultural landscapes. J. Appl. Ecol. 56, 31–43 (2019).
Article  Google Scholar 

59.
McCallum, M. L. Tropical Herpetology: A drop in the bucket. Trends Ecol. Evol. 20, 289–290 (2005).
Article  Google Scholar 

60.
Bonnet, X., Shine, R. & Lourdais, O. Taxonomic chauvinism. Trends Ecol. Evol. 17, 1–3 (2002).
Article  Google Scholar 

61.
Tingley, R., Meiri, S. & Chapple, D. G. Addressing knowledge gaps in reptile conservation. Biol. Cons. 204, 1–5 (2016).
Article  Google Scholar 

62.
Beier, P., Majka, D. & Jenness, J. Conceptual Steps for Designing Wildlife Corridors. www.corridordesign.org (2007).

63.
Valencia-Aguilar, A., Cortés-Gómez, A. M. & Ruiz-Agudelo, C. A. Ecosystem services provided by amphibians and reptiles in Neotropical ecosystems. Int. J. Biodivers. Sci. Ecosyst. Serv. Manag. 9, 257–272 (2013).
Article  Google Scholar 

64.
Hager, H. A. Area-sensitivity of reptiles and amphibians: Are there indicator species for habitat fragmentation?. Écoscience 66, 139–147 (1998).
Article  Google Scholar 

65.
Almasieh, K., Mirghazanfari, S. M. & Mahmoodi, S. Biodiversity hotspots for modeled habitat patches and corridors of species richness and threatened species of reptiles in central Iran. Eur. J. Wildl. Res. 65, 92. https://doi.org/10.1007/s10344-019-1335-x (2019).
Article  Google Scholar 

66.
Albert, C., Luque, G. M. & Courchamp, F. The twenty most charismatic species. PLoS ONE 13, e0199149. https://doi.org/10.1371/journal.pone.0199149 (2018).
CAS  Article  Google Scholar 

67.
Brooke, Z. M., Bielby, J., Nambiar, K. & Carbone, C. correlates of research effort in carnivores: Body size, range size and diet matter. PLoS ONE 9, e93195. https://doi.org/10.1371/journal.pone.0093195 (2014).
ADS  CAS  Article  Google Scholar 

68.
Rozylowicz, L., Popescu, V. D., Pătroescu, M. & Chișamera, G. The potential of large carnivores as conservation surrogates in the Romanian Carpathians. Biodivers. Conserv. 20, 561–579 (2011).
Article  Google Scholar 

69.
Beier, P., Majka, D. R. & Spencer, W. D. Forks in the road choices in procedures for designing wildland linkages. Conserv. Biol. 22, 836–851 (2008).
Article  Google Scholar 

70.
Burbrink, F. T., Phillips, C. A. & Heske, E. J. A riparian zone in southern Illinois as a potential dispersal corridor for reptiles and amphibians. Biol. Cons. 86, 107–115 (1998).
Article  Google Scholar 

71.
Dixo, M. & Metzger, J. P. Are corridors, fragment size and forest structure important for the conservation of leaf-litter lizards in a fragmented landscape?. Oryx 43, 435–442 (2009).
Article  Google Scholar 

72.
ArcGIS Release 10.4 (Redlands, CA, 2013).

73.
Hamer, A. J. & McDonnell, M. J. The response of herpetofauna to urbanization: Inferring patterns of persistence from wildlife databases. Austral Ecol. 35, 568–580 (2010).
Article  Google Scholar 

74.
Vignoli, L., Mocaer, I., Luiselli, L. & Bologna, M. A. Can a large metropolis sustain complex herpetofauna communities? An analysis of the suitability of green space fragments in Rome. Anim. Conserv. 12, 456–466 (2009).
Article  Google Scholar 

75.
Strugariu, A., Gherghel, I., Huțuleac-Volosciuc, M. V. & Pușcașu, C. M. Preliminary aspects concerning the herpetofauna from urban and peri-urban environments from North-Eastern Romania: A case study in the city of Suceava. Herpetol. Roman. 1, 53–61 (2007).
Google Scholar 

76.
Gherghel, I., Strugariu, A., Sahlean, T. C. & Zamfirescu, O. Anthropogenic impact or anthropogenic accommodation? Distribution range expansion of the common wall lizard (Podarcis muralis) by means of artificial habitats in the north-eastern limits of its distribution range. Acta Herpetol. 4, 183–189 (2009).
Google Scholar 

77.
Gherghel, I. & Tedrow, R. Manmade structures are used by an invasive species to colonize new territory across a fragmented landscape. Acta Oecol. 101, 103479. https://doi.org/10.1016/j.actao.2019.103479 (2019).
Article  Google Scholar 

78.
Ward, M. et al. Just ten percent of the global terrestrial protected area network is structurally connected via intact land. Nat. Commun. 11, 4563. https://doi.org/10.1038/s41467-020-18457-x (2020).
CAS  Article  Google Scholar 

79.
Watson, J. E. M. et al. The exceptional value of intact forest ecosystems. Nat. Ecol. Evol. 2, 599–610 (2018).
Article  Google Scholar 

80.
Fitzgerald, L. A. et al. The future for reptiles: advances and challenges in the anthropocene. Encycl. Anthropocene 3, 163–174 (2018).
Article  Google Scholar 

81.
Hof, C., Araújo, M. B., Jetz, W. & Rahbek, C. Additive threats from pathogens, climate and land-use change for global amphibian diversity. Nature 480, 516–519 (2011).
ADS  CAS  Article  Google Scholar 

82.
Meta-Analysis, A. Rey Benayas, J. M., Newton, A. C., Diaz, A. & Bullock, J. M. Enhancement of biodiversity and ecosystem services by ecological restoration. Science 325, 1121–1124 (2009).
Article  CAS  Google Scholar 

83.
Van Der Windt, H. J. & Swart, J. A. A. Ecological corridors, connecting science and politics: the case of the Green River in the Netherlands. J. Appl. Ecol. 45, 124–132 (2008).
Article  Google Scholar 

84.
Hilty, J. et al. Guidelines for conserving connectivity through ecological networks and corridors (International Union for Conservation of Nature, 2020).

85.
Gregory, A. J. & Beier, P. Response variables for evaluation of the effectiveness of conservation corridors. Conserv. Biol. 28, 689–695 (2014).
Article  Google Scholar 

86.
Mráz, P. & Ronikier, M. Biogeography of the Carpathians: Evolutionary and spatial facets of biodiversity. Biol. J. Lin. Soc. 119, 528–559 (2016).
Article  Google Scholar 

87.
Deodatus, F. et al. Creation of ecological corridors in the Ukrainian Carpathians. In The Carpathians: Integrating Nature and Society Towards Sustainability Environmental Science and Engineering (eds Kozak, J. et al.) 701–717 (Springer, Berlin, 2013).
Google Scholar 

88.
Favilli, F., Hoffmann, C., Elmi, M., Ravazzoli, E. & Streifeneder, T. The BioREGIO Carpathians project: Aims, methodology and results from the “Continuity and Connectivity” analysis. Nat. Conserv. 11, 95–111 (2015).
Article  Google Scholar 

89.
Csagoly, P., Magnin, G. & Hulea, O. Lower Danube Green Corridor. in The Wetland Book: II: Distribution, Description and Conservation (eds. Finlayson, M. C., Milton, R. G., Prentice, C. R. & Davidson, N. C.) 1–6 (Springer, Netherlands, 2016).

90.
Belote, R. T. et al. Identifying corridors among large protected areas in the United States. PLoS ONE 11, e0154223. https://doi.org/10.1371/journal.pone.0154223 (2016).
CAS  Article  Google Scholar 

91.
Breckheimer, I. et al. Defining and evaluating the umbrella species concept for conserving and restoring landscape connectivity. Conserv. Biol. 28, 1584–1593 (2014).
Article  Google Scholar 

92.
Meurant, M., Gonzales, A., Doxa, A. & Albert, C. H. Selecting surrogate species for connectivity conservation. Biol. Cons. 227, 326–334 (2018).
Article  Google Scholar 

93.
Dondina, O., Orioli, V., Chiatante, G. & Bani, L. Practical insights to select focal species and design priority areas for conservation. Ecol. Indic. 108, 105767. https://doi.org/10.1016/j.ecolind.2019.105767 (2020).
Article  Google Scholar 

94.
Churko, G., Kienast, F. & Bolliger, J. A multispecies assessment to identify the functional connectivity of amphibians in a human-dominated landscape. Int. J. Geo-Inf. 9, 287. https://doi.org/10.3390/ijgi9050287 (2020).
Article  Google Scholar 

95.
Cushman, S. A. & Landguth, E. L. Multi-taxa population connectivity in the Northern Rocky Mountains. Ecol. Model. 231, 101–112 (2012).
Article  Google Scholar 

96.
Krosby, M. et al. Focal species and landscape “naturalness” corridor models offer complementary approaches for connectivity conservation planning. Landsc. Ecol. 30, 2121–2132 (2015).
Article  Google Scholar 

97.
Wiens, J. A., Hayward, G. D., Holthausen, R. S. & Wisdom, M. J. Using surrogate species and groups for conservation planning and management. Bioscience 58, 241–252 (2008).
Article  Google Scholar 

98.
Macdonald, E. A. et al. Identifying ambassador species for conservation marketing. Glob. Ecol. Conserv. 12, 204–214 (2017).
Article  Google Scholar 

99.
Fleury, A. M. & Brown, R. D. A framework for the design of wildlife conservation corridors with specific application to southwestern Ontario. Landsc. Urban Plan. 37, 163–186 (1997).
Article  Google Scholar 

100.
Cogălniceanu, D. et al. Diversity and distribution of amphibians in Romania. ZooKeys 296, 35–57 (2013).
Article  Google Scholar 

101.
Cogălniceanu, D. et al. Diversity and distribution of reptiles in Romania. ZooKeys 341, 49–76 (2013).
Article  Google Scholar 

102.
LaRue, M. A. & Nielsen, C. K. Modelling potential dispersal corridors for cougars in midwestern North America using least-cost path methods. Ecol. Model. 212, 372–381 (2008).
Article  Google Scholar 

103.
Adriaensen, F. et al. The application of ‘least-cost’ modelling as a functional landscape model. Landsc. Urban Plan. 64, 233–247 (2003).
Article  Google Scholar 

104.
Correa Ayram, C. A., Mendoza, M. E., Etter, A. & Salicrup, D. R. P. Habitat connectivity in biodiversity conservation: a review of recent studies and applications. Prog. Phys. Geogr. 1, 1–32 (2015).
Google Scholar 

105.
Ribeiro, J. W. et al. LandScape Corridors (LSCORRIDORS): A new softwarepackage for modelling ecological corridors based onlandscape patterns and species requirements. Methods Ecol. Evol. 8, 1425–1432 (2017).
Article  Google Scholar 

106.
Linkage Mapper Connectivity Analysis Software v. 2.0.0 (The Nature Conservancy, Seattle, 2011).

107.
Popescu, V. D., Rozylowicz, L., Cogălniceanu, D., Niculae, I. M. & Cucu, A. L. Moving into protected areas? Setting conservation priorities for Romanian Reptiles and Amphibians at risk from climate change. PLoS ONE 8, e79330. https://doi.org/10.1371/journal.pone.0079330 (2014).
ADS  CAS  Article  Google Scholar 

108.
Lambeck, R. J. Focal species: A multispecies umbrella for nature conservation. Conserv. Biol. 11, 849–856 (1997).
Article  Google Scholar  More