Ecology

Araújo MS, Perez SI, Magazoni MJC, Petry AC (2014) Body size and allometric shape variation in the molly Poecilia vivipara along a gradient of salinity and predation. BMC Evol Biol 14:251PubMed 

Google Scholar 
Arendt JD (2010) Morphological correlates of sprint swimming speed in five species of spadefoot toad tadpoles: comparison of morphometric methods. J Morphol 271:1044–1052PubMed 

Google Scholar 
Arendt JD, Reznick DN (2008) Convergence and parallelism reconsidered: what have we learned about the genetics of adaptation? Trends Ecol Evol 23:26–32PubMed 

Google Scholar 
Arnett HA, Kinnison MT (2017) Predator-induced phenotypic plasticity of shape and behavior: parallel and unique patterns across sexes and species. Curr Zool 63:369–378PubMed 

Google Scholar 
Arnett HA (2016) Sources of ecologically important trait variation in mosquitofish (Gambusia affinis and Gambusia holbrooki). Thesis, University of MaineArnold SJ (1983) Morphology, performance and fitness. Am Zool 23:347–361
Google Scholar 
Avise JC (1989) Gene trees and organismal histories: a phylogenetic approach to population biology. Evolution 43:1192–1208PubMed 

Google Scholar 
Baldwin BG (1997) Adaptive radiation of the Hawaiian silversword alliance: congruence and conflict of phylogenetic evidence from molecular and non-molecular investigations. In: Givnish TJ, Sytsma KJ (eds.) Molecular evolution and adaptive radiation. Cambridge University Press, Cambridge, UK, p 103–128Belk MC, Tuckfield RC (2010) Changing costs of reproduction: age‐based differences in reproductive allocation and escape performance in a livebearing fish. Oikos 119:163–169
Google Scholar 
Blount ZD, Lenski RE, Losos JB (2018) Contingency and determinism in evolution: replaying life’s tape. Science 362:eaam5979.PubMed 

Google Scholar 
Bryant EH, Meffert LM (1993) The effect of serial founder-flush cycles on quantitative genetic variation in the housefly. Heredity 70:122–129
Google Scholar 
Calsbeek R, Kuchta S (2011) Predator mediated selection and the impact of developmental stage on viability in wood frog tadpoles (Rana sylvatica). BMC Evol Biol 11:353PubMed 

Google Scholar 
Chapuis M-P, Estoup A (2007) Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol 24:621–631CAS 
PubMed 

Google Scholar 
Chenoweth SF, Blows MW (2008) QST meets the G matrix: the dimensionality of adaptive divergence in multiple correlated quantitative traits. Evolution 62:1437–1449PubMed 

Google Scholar 
Constantz GD (1989) Reproductive biology of poeciliid fishes. In: Meffe GK Jr, Snelson FF (eds.) Ecology and evolution of livebearing fishes (Poeciliidae). Prentice Hall, Englewood Cliffs, NJ, p 33–50Cunha RL, Tenorio MJ, Afonso C, Castilho R, Zardoya R (2008) Replaying the tape: recurring biogeographical patterns in Cape Verde Conus after 12 million years. Mol Ecol 17:885–901PubMed 

Google Scholar 
Dale MR, Fortin M-J (2014) Spatial analysis: a guide for ecologists. Cambridge University Press, Cambridge, UKDarwin CE (1859) The origin of species and the descent of man. The Modern Library, New York, NYDay T, Pritchard J, Schluter D (1994) A comparison of two sticklebacks. Evolution 48:1723–1734PubMed 

Google Scholar 
Dayton GH, Saenz D, Baum KA, Langerhans RB, DeWitt TJ (2005) Body shape, burst speed and escape behavior of larval anurans. Oikos 111:582–591
Google Scholar 
DeWitt TJ, Fuentes JI, Ioerger TR, Bishop MP (2021) Rectifying I: three point and continuous fit of the spatial autocorrelation metric, Moran’s I, to ideal form. Landsc Ecol 36:2897–2918
Google Scholar 
DeWitt TJ, Scheiner SM (2004) Phenotypic variation from single genotypes: a primer. In: DeWitt TJ, Scheiner SM (eds.) Phenotypic plasticity: functional and conceptual approaches. Oxford University Press, New York, NY, p 1–9DeWoody J, Avise J (2000) Microsatellite variation in marine, freshwater and anadromous fishes compared with other animals. J Fish Biol 56:461–473CAS 

Google Scholar 
Dobzhansky T (1955) A review of some fundamental concepts and problems of population genetics. Cold Spring Harb Symp Quant Biol 20:1–15CAS 
PubMed 

Google Scholar 
Endler JA (1986) Natural selection in the wild. Princeton University Press, Princeton, NJEroukhmanoff F, Hargeby A, Arnberg NN, Hellgren O, Bensch S, Svensson EI (2009) Parallelism and historical contingency during rapid ecotype divergence in an isopod. J Evol Biol 22:1098–1110CAS 
PubMed 

Google Scholar 
Fisher RA (1930) The genetical theory of natural selection. Oxford University Press, Oxford, UKFranssen NR (2011) Anthropogenic habitat alteration induces rapid morphological divergence in a native stream fish. Evol Appl 4:791–804PubMed 

Google Scholar 
Futuyma DJ, Moreno G (1988) The evolution of ecological specialization. Annu Rev Ecol Syst 19:207–233
Google Scholar 
Futuyma DJ (2021) How does phenotypic plasticity fit into evolutionary theory? In: Pfennig DW (ed) Phenotypic plasticity & evolution. CRC Press, Boca Raton, FL, p 349–366Ghalambor CK, Reznick DN, Walker JA (2004) Constraints on adaptive evolution: the functional trade-off between reproduction and fast-start swimming performance in the Trinidadian guppy (Poecilia reticulata). Am Nat 164:38–50PubMed 

Google Scholar 
Givnish TJ, Knox E, Patterson TB, Hapeman JR, Palmer JB, Sytsma KJ (1996) The Hawaiian lobelioids are monophyletic and underwent a rapid initial radiation roughly 15 million years ago. Am J Bot 83:159
Google Scholar 
Gomes JL, Montiero L (2008) Morphological divergence patterns among populations of Poecilia vivipara (Teleostei Poeciliidae): test of an ecomorphological paradigm Biol J Linn Soc 93:799–812
Google Scholar 
Gompel N, Prud’homme B (2009) The causes of repeated genetic evolution. Dev Biol 332:36–47CAS 
PubMed 

Google Scholar 
Grant PR, Grant BR (2014) 40 years of evolution: Darwin’s finches on Daphne Major Island. Princeton University Press, Princeton, NJGreenway R, Barts N, Henpita C, Brown AP, Rodriguez LA, Peña CMR et al. (2020) Convergent evolution of conserved mitochondrial pathways underlies repeated adaptation to extreme environments. Proc Natl Acad Sci USA 117:16424–16430CAS 
PubMed 

Google Scholar 
Hartl DL, Clark AG (1997) Principles of population genetics. Sinauer, Sunderland, MA
Google Scholar 
Haynes JL (1993) Annual reestablishment of mosquitofish populations in Nebraska. Copeia 1993:232–235
Google Scholar 
Holsinger KE, Weir BS (2009) Genetics in geographically structured populations: defining, estimating and interpreting FST. Nat Rev Genet 10:639–650CAS 
PubMed 

Google Scholar 
Ingley SJ, Johnson JB (2016) Divergent natural selection promotes immigrant inviability at early and late stages of evolutionary divergence. Evolution 70:600–616PubMed 

Google Scholar 
Ingley SJ, Billman EJ, Belk MC, Johnson JB (2014) Morphological divergence driven by predation environment within and between species of Brachyrhaphis fishes. PloS One 9:90274
Google Scholar 
James ME, Wilkinson MJ, Bernal DM, Liu H, North HL, Engelstädter J, et al. (2021) Phenotypic and genotypic parallel evolution in parapatric ecotypes of Senecio. Evolution (In press). Online version of record https://doi.org/10.1111/evo.14387Jiang S, Zhu K, Han L, Chen C, Wang M, Wang X (2021) Genetic variation and phylogeographic structure of Laodelphax striatellus in China based on microsatellite markers. J Appl Entomol 145:336–347CAS 

Google Scholar 
Johnson JB, Burt DB, DeWitt TJ (2008) Form, function, fitness-pathways to survival. Evolution 62:1243–1251PubMed 

Google Scholar 
Kobza RM, Trexler JC, Loftus WF, Perry SA (2004) Community structure of fishes inhabiting aquatic refuges in a threatened Karst wetland and its implications for ecosystem management. Biol Conserv 116:153–165
Google Scholar 
Langerhans RB (2009) Morphology, performance, fitness: functional insight into a post-Pleistocene radiation of mosquitofish. Biol Lett 5:488–491PubMed 

Google Scholar 
Langerhans RB, DeWitt TJ (2004) Shared and unique features of evolutionary diversification. Am Nat 164:335–349PubMed 

Google Scholar 
Langerhans RB, Makowicz AM (2009) Shared and unique features of morphological differentiation between predator regimes in Gambusia caymanensis. J Evol Biol 22:2231–2242CAS 
PubMed 

Google Scholar 
Langerhans RB, Layman CA, DeWitt TJ (2005) Male genital size reflects a tradeoff between attracting mates and avoiding predators in two live-bearing fish species. Proc Natl Acad Sci USA 102:7618–7623CAS 
PubMed 

Google Scholar 
Langerhans RB, Gifford ME, Joseph EO (2007) Ecological speciation in Gambusia fishes. Evolution 61:2056–2074CAS 
PubMed 

Google Scholar 
Langerhans RB, Layman CA, Shokrollahi AM, DeWitt TJ (2004) Predator-driven phenotypic diversification in Gambusia affinis. Evolution 58:2305–2318PubMed 

Google Scholar 
Leger EA, Rice KJ (2007) Assessing the speed and predictability of local adaptation in invasive California poppies (Eschscholzia californica). J Evol Biol 20:1090–1103CAS 
PubMed 

Google Scholar 
Levins R (1968) Evolution in changing environments. Princeton University Press, Princeton, NJLoera-Pérez J, Hernández-Stefanoni JL, Chiappa-Carrara X (2020) How do spatial and environmental factors affect the fish community structure in seasonally flooded karst systems? Lat Am J Aquat Res 48:268–279
Google Scholar 
Losos JB, Jackman TR, Larson A, de Queiroz K, Rodrı́guez-Schettino L (1998) Contingency and determinism in replicated adaptive radiations of island lizards. Science 279:2115–2118CAS 
PubMed 

Google Scholar 
Losos JB (2009) Lizards in an evolutionary tree: ecology and adaptive radiation of anoles. University of California Press, Berkeley, CAMaglio VJ, Rosen DE (1969) Changing preference for substrate color by reproductively active mosquitofish, Gambusia affinis (Baird and Girard)(Poeciliidae, Atheriniformes). American Museum novitates; no. 2397Martin RG(1975) Sexual and aggressive behavior, density and social structure in a natural population of mosquitofish, Gambusia affinis holbrooki. Copeia 1975:445–454
Google Scholar 
Maruyama T, Fuerst PA (1985) Population bottlenecks and nonequilibrium models in population genetics. II. Number of alleles in a small population that was formed by a recent bottleneck. Genetics 111:675–689CAS 
PubMed 

Google Scholar 
Matthews WJ, Marsh-Matthews E (2011) An invasive fish species within its native range: community effects and population dynamics of Gambusia affinis in the central United States. Freshw Biol 56:2609–2619
Google Scholar 
Mays JR, DeWitt TJ, Dharampal P, Andrus FT, Findlay RH (2019) Frequent habitat migration, phenotypic plasticity, and vestigial ecophenotypy revealed by isotope-based natal habitat inference in bluegill sunfish, Lepomis macrochirus. Evol Ecol Res 20. Available from https://evolutionary-ecology.com/abstracts/v20/3235.htmlMilano D, Ruzzante DE, Cussac VE, Macchi PJ, Ferriz RA, Barriga JP et al. (2006) Latitudinal and ecological correlates of morphological variation in Galaxias platei (Pisces, Galaxiidae) in Patagonia. Biol J Linn Soc 87:69–82
Google Scholar 
Moen DS, Morlon H, Wiens JJ (2016) Testing convergence versus history: convergence dominates phenotypic evolution for over 150 million years in frogs. Syst Biol 65:146–160PubMed 

Google Scholar 
Moody EK, Lozano-Vilano ML (2018) Predation drives morphological convergence in the Gambusia panuco species group among lotic andlentic habitats. J Evol Biol 31:491–501CAS 
PubMed 

Google Scholar 
Nei M (1978) Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89:583–590CAS 
PubMed 

Google Scholar 
Nepokroeff M, Sytsma KJ (1996) Systematics and patterns of speciation and colonization in Hawaiian Psychotria and relatives based on phylogenetic analysis of ITS sequence data. Am J Bot 83:181–182
Google Scholar 
Nievergelt CM, Libiger O, Schork NJ (2007) Generalized analysis of molecular variance. PLoS Genet 3:e51PubMed 

Google Scholar 
Oke KB, Rolshausen G, LeBlond C, Hendry AP (2017) How parallel is parallel evolution? A comparative analysis in fishes. Am Nat 190:1–16PubMed 

Google Scholar 
Ord TJ, Summers TC (2015) Repeated evolution and the impact of evolutionary history on adaptation. BMC Evol Biol 15:1–12
Google Scholar 
Pandey P, Ramegowda V, Senthil-Kumar M (2015) Shared and unique responses of plants to multiple individual stresses and stress combinations: physiological and molecular mechanisms. Front Plant Sci 6:723PubMed 

Google Scholar 
Pazmino SD, Kent MI, Ward AJ (2020) Locomotion and habituation to novel experimental environments in a social fish species. Behaviour 1:1–17
Google Scholar 
Peakall R, Smouse PE (2006) GenAlEx 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes 6:288–295
Google Scholar 
Piry S, Luikart G, Cornuet J-M (1999) Computer note. BOTTLENECK: a computer program for detecting recent reductions in the effective size using allele frequency data. J Hered 90:502–503
Google Scholar 
Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155:945–959CAS 
PubMed 

Google Scholar 
Purcell KM, Lance SL, Jones KLStockwell CA (2011) Ten novel microsatellite markers for the western mosquitofish Gambusia affinis Conserv Genet Resour 3:361–363
Google Scholar 
Putman AI, Carbone I (2014) Challenges in analysis and interpretation of microsatellite data for population genetic studies. Ecol Evol 4:4399–4428PubMed 

Google Scholar 
Pyke GH (2005) A review of the biology of Gambusia affinis and G. holbrooki. Rev Fish Biol Fish 15:339–365
Google Scholar 
Reed DH, Frankham R (2001) How closely correlated are molecular and quantitative measures of genetic variation? A meta-analysis. Evolution 55:1095–1103CAS 
PubMed 

Google Scholar 
Reznick DN, Shaw FH, Rodd FH, Shaw RG (1997) Evaluation of the rate of evolution in natural populations of guppies (Poecilia reticulata). Science 275:1934–1937CAS 
PubMed 

Google Scholar 
Richards TJ, Walter GM, McGuigan K, Ortiz-Barrientos D (2016) Divergent natural selection drives the evolution of reproductive isolation in an Australian wildflower Evolution 70:1993–2003PubMed 

Google Scholar 
Rivera G (2008) Ecomorphological variation in shell shape of the freshwater turtle Pseudemys concinna inhabiting different aquatic flow regimes. Integr Comp Biol 48:769–787PubMed 

Google Scholar 
Robinson BW, Wilson DS (1996) Genetic variation and phenotypic plasticity in a trophically polymorphic population of pumpkinseed sunfish (Lepomis gibbosus). Evol Ecol 10:631–652
Google Scholar 
Ruehl CB, DeWitt TJ (2005) Trophic plasticity and fine-grained resource variation in populations of western mosquitofish, Gambusia affinis. Evol Ecol Res 7:801–819
Google Scholar 
Ruehl CB, Shervette V, DeWitt TJ (2011) Replicated shape variation between simple and complex habitats in two estuarine fishes. Biol J Linn Soc 103:147–158
Google Scholar 
Santi F, Petry AC, Plath M, Riesch R (2020) Phenotypic differentiation in a heterogeneous environment: morphological and life-history responses to ecological gradients in a livebearing fish. J Zool 310:10–23
Google Scholar 
Schluter D, Clifford EA, Nemethy M, McKinnon JS (2004) Parallel evolution and inheritance of quantitative traits. Am Nat 163:809–822PubMed 

Google Scholar 
Schluter D (2000) The ecology of adaptive radiation. Oxford University Press, Oxford, UKSharpe DM, Langerhans RB, Low-Décarie E, Chapman LJ (2015) Little evidence for morphological change in a resilient endemic species following the introduction of a novel predator. J Evol Biol 28:2054–2067CAS 
PubMed 

Google Scholar 
Slatkin M (1995) A measure of population subdivision based on microsatellite allele frequencies. Genetics 139:457–462CAS 
PubMed 

Google Scholar 
Spencer CC, Chlan CA, Neigel JE, Scribner KT, Wooten MC, Leberg PL (1999) Polymorphic microsatellite markers in the western mosquitofish, Gambusia affinis. Mol Ecol 8:157–168CAS 
PubMed 

Google Scholar 
Spitze K (1993) Population structure in Daphnia obtusa: quantitative genetic and allozymic variation. Genetics 135:367–374CAS 
PubMed 

Google Scholar 
Thibault RE, Schultz RJ (1978) Reproductive adaptations among viviparous fishes (Cyprinodontiformes: Poeciliidae). Evolution 32:320–333PubMed 

Google Scholar 
Tobler M, DeWitt TJ, Schlupp I, García de León FJ, Herrmann R, Feulner PG et al. (2008) Toxic hydrogen sulfide and dark caves: phenotypic and genetic divergence across two abiotic environmental gradients in Poecilia mexicana. Evolution 62:2643–2659PubMed 

Google Scholar 
Tobler M, Palacios M, Chapman LJ, Mitrofanov I, Bierbach D, Plath M et al. (2011) Evolution in extreme environments: replicated phenotypic differentiation in livebearing fish inhabiting sulfidic springs. Evolution 65:2213–2228PubMed 

Google Scholar 
Van Oosterhout C, Weetman D, Hutchinson WF (2006) Estimation and adjustment of microsatellite null alleles in nonequilibrium populations. Mol Ecol Notes 6:255–256
Google Scholar 
Vázquez-Domínguez E, Hernández-Valdés A, Rojas-Santoyo A, Zambrano L (2009) Contrasting genetic structure in two codistributed freshwater fish species of highly seasonal systems. Rev Mex Biodivers 80:181–192
Google Scholar 
Via S, Lande R (1985) Genotype-environment interaction and the evolution of phenotypic plasticity. Evolution 39:505–522PubMed 

Google Scholar 
Waddington CH (1957) The strategy of the genes. Allen & Unwin, London
Google Scholar 
Walker JA (1997) Ecological morphology of lacustrine threespine stickleback Gasterosteus aculeatus L. (Gasterosteidae) body shape. Biol J Linn Soc 61:3–50
Google Scholar 
Walker JA, Bell MA (2000) Net evolutionary trajectories of body shape evolution within a microgeographic radiation of threespine sticklebacks (Gasterosteus aculeatus). J Zool 252:293–302
Google Scholar 
Wang X, Zorraquino V, Kim M, Tsoukalas A, Tagkopoulos I (2018) Predicting the evolution of Escherichia coli by a data-driven approach. Nat Commun 9:1–12
Google Scholar 
Ward RD, Woodwark M, Skibinski DOF (1994) A comparison of genetic diversity levels in marine, freshwater, and anadromous fishes. J Fish Biol 44:213–232
Google Scholar 
Waters JM, McCulloch GA (2021) Reinventing the wheel? Reassessing the roles of gene flow, sorting and convergence in repeated evolution. Mol Ecol 30:4162–4172PubMed 

Google Scholar 
Zambrano L, Vázquez-Domínguez E, García-Bedoya D, Loftus WF, Trexler JC (2006) Fish community structure in freshwater karstic water bodies of the Sian Ka’an Reserve in the Yucatan peninsula, Mexico. Ichthyol Explor Freshw 17:193–206Zane L, Nelson WS, Jones AG, Avise JC (1999) Microsatellite assessment of multiple paternity in natural populations of a live-bearing fish, Gambusia holbrooki. J Evol Biol 12:61–69
Google Scholar  More

1.Jarau, S. & Hrncir, M. Food Exploitation by Social Insects: Ecological, Behavioral, and Theoretical Approaches. (CRC Press, 2009).2.Detrain, C. & Deneubourg, J.-L. Collective decision-making and foraging patterns in ants and honeybees. Adv. Insect Physiol. 35, 123–173 (2008) (Elsevier).
Google Scholar 
3.Hölldobler, B. & Wilson, E. O. The Ants (Harvard University Press, 1990).
Google Scholar 
4.Beckers, R., Deneubourg, J.-L. & Goss, S. Modulation of trail laying in the ant Lasius niger (Hymenoptera: Formicidae) and its role in the collective selection of a food source. J. Insect Behav. 6, 751–759 (1993).
Google Scholar 
5.Detrain, C. & Prieur, J. Sensitivity and feeding efficiency of the black garden ant Lasius niger to sugar resources. J. Insect Physiol. 64, 74–80 (2014).CAS 
PubMed 

Google Scholar 
6.Jackson, D. E. & Châline, N. Modulation of pheromone trail strength with food quality in Pharaoh’s ant, Monomorium pharaonic. Animal Behav. 74, 463–470 (2007).
Google Scholar 
7.Sumpter, D. J. T. & Beekman, M. From nonlinearity to optimality: Pheromone trail foraging by ants. Anim. Behav. 66, 273–280 (2003).
Google Scholar 
8.Cerdá, X., Angulo, E., Boulay, R. & Lenoir, A. Individual and collective foraging decisions: A field study of worker recruitment in the gypsy ant Aphaenogaster senilis. Behav. Ecol. Sociobiol. 63, 551–562 (2009).
Google Scholar 
9.Detrain, C. & Deneubourg, J.-L. Scavenging by Pheidole pallidula key for understanding decision-making systems in ants. Anim. Behav. 53, 537–547 (1997).
Google Scholar 
10.Mailleux, A.-C., Deneubourg, J. L. & Detrain, C. How do ants assess food volume?. Anim. Behav. 59, 1061–1069 (2000).CAS 
PubMed 

Google Scholar 
11.Breed, M. D., Fewell, J. H., Moore, A. J. & Williams, K. R. Graded recruitment in a ponerine ant. Behav. Ecol. Sociobiol. 20, 407–411 (1987).
Google Scholar 
12.Cammaerts, M.-C. & Cammaerts, R. Food recruitment strategies of the ants Myrmica sabuleti and Myrmica ruginodis. Behav. Proc. 5, 251–270 (1980).CAS 

Google Scholar 
13.Portha, S., Deneubourg, J.-L. & Detrain, C. Self-organized asymmetries in ant foraging: A functional response to food type and colony needs. Behav. Ecol. 13, 776–781 (2002).
Google Scholar 
14.Devigne, C. & Detrain, C. How does food distance influence foraging in the ant Lasius niger: The importance of home-range marking. Insect. Soc. 53, 46–55 (2006).
Google Scholar 
15.Fewell, J. H. Directional fidelity as a foraging constraint in the western harvester ant, Pogonomyrmex occidentalis. Oecologia 82, 45–51 (1990).ADS 
PubMed 

Google Scholar 
16.Howard, D. F. & Tschinkel, W. R. The effect of colony size and starvation on food flow in the fire ant, Solenopsis invicta (Hymenoptera: Formicidae). Behav. Ecol. Sociobiol. 7, 293–300 (1980).
Google Scholar 
17.Mailleux, A.-C., Devigne, C., Deneubourg, J.-L. & Detrain, C. Impact of starvation on Lasius niger’ exploration. Ethology 116, 248–256 (2010).
Google Scholar 
18.Portha, S., Deneubourg, J.-L. & Detrain, C. How food type and brood influence foraging decisions of Lasius niger scouts. Anim. Behav. 68, 115–122 (2004).
Google Scholar 
19.Deneubourg, J.-L., Goss, S., Pasteels, J. M. & Beckers, R. Collective decision making through food recruitment. Insectes Soc. 37, 258–267 (1990).
Google Scholar 
20.Czaczkes, T. J., Salmane, A. K., Klampfleuthner, F. A. M. & Heinze, J. Private information alone can trigger trapping of ant colonies in local feeding optima. J. Exp. Biol. 219, 744–751 (2016).PubMed 

Google Scholar 
21.Collett, T. S. & Collett, M. Memory use in insect visual navigation. Nat. Rev. Neurosci. 3, 542–552 (2002).CAS 
PubMed 
MATH 

Google Scholar 
22.Azevedo, D. L. O., Medeiros, J. C. & Araújo, A. Adjustments in the time, distance and direction of foraging in Dinoponera quadriceps Workers. J. Insect. Behav. 27, 177–191 (2014).
Google Scholar 
23.Beverly, B. D., McLendon, H., Nacu, S., Holmes, S. & Gordon, D. M. How site fidelity leads to individual differences in the foraging activity of harvester ants. Behav. Ecol. 20, 633–638 (2009).
Google Scholar 
24.Fourcassié, V. & Traniello, J. F. A. Food searching behaviour in the ant Formica schaufussi (Hymenoptera, Formicidae): Response of naive foragers to protein and carbohydrate food. Anim. Behav. 48, 69–79 (1994).
Google Scholar 
25.Aron, S., Beckers, R., Deneubourg, J. L. & Pasteels, J. M. Memory and chemical communication in the orientation of two mass-recruiting ant species. Ins. Soc. 40, 369–380 (1993).
Google Scholar 
26.Lehue, M., Detrain, C. & Collignon, B. Nest entrances, spatial fidelity, and foraging patterns in the red ant Myrmica rubra: A field and theoretical study. Insects 11, 317 (2020).PubMed Central 

Google Scholar 
27.Bolek, S., Wittlinger, M. & Wolf, H. What counts for ants? How return behaviour and food search of Cataglyphis ants are modified by variations in food quantity and experience. J. Exp. Biol. 215, 3218–3222 (2012).PubMed 

Google Scholar 
28.Detrain, C., Natan, C. & Deneubourg, J. L. The influence of the physical environment on the self-organised foraging patterns of ants. Naturwissenschaften 88, 171–174 (2001).ADS 
CAS 
PubMed 

Google Scholar 
29.Traniello, J. F. A., Fujita, M. S. & Bowen, R. V. Ant foraging behavior: Ambient temperature influences prey selection. Behav. Ecol. Sociobiol. 15, 65–68 (1984).
Google Scholar 
30.Nonacs, P. & Dill, L. M. Mortality risk versus food quality trade-offs in ants: Patch use over time. Ecol. Entomol. 16, 73–80 (1991).
Google Scholar 
31.Tanner, C. J. Individual experience-based foraging can generate community territorial structure for competing ant species. Behav. Ecol. Sociobiol. 63, 591–603 (2009).
Google Scholar 
32.Brown, M. J. F. & Gordon, D. M. How resources and encounters affect the distribution of foraging activity in a seed-harvesting ant. Behav. Ecol. Sociobiol. 47, 195–203 (2000).
Google Scholar 
33.Fourcassié, V., Schmitt, T. & Detrain, C. Impact of interference competition on exploration and food exploitation in the ant Lasius niger. Psyche J. Entomol. 2012, 1–8 (2012).
Google Scholar 
34.Mehdiabadi, N. & Gilbert, L. Colony-level impacts of parasitoid flies on fire ants. Proceedings of the Royal Society of London. Series B: Biological Sciences. 269, 1695–1699 (2002).35.Feener, D. H. Competition between ant species: Outcome controlled by parasitic flies. Science 214, 815–817 (1981).ADS 
PubMed 

Google Scholar 
36.Schmid-Hempel, P. Parasites in Social Insects (Princeton University Press, 1998).
Google Scholar 
37.Cremer, S., Armitage, S. A. O. & Schmid-Hempel, P. Social immunity. Curr. Biol. 17, 693–702 (2007).
Google Scholar 
38.Zhukovskaya, M., Yanagawa, A. & Forschler, B. Grooming behavior as a mechanism of insect disease defense. Insects 4, 609–630 (2013).PubMed 
PubMed Central 

Google Scholar 
39.Ortius-Lechner, D., Maile, R., Morgan, E. D. & Boomsma, J. J. Metapleural gland secretion of the leaf-cutter ant Acromyrmex octospinosus: New compounds and their functional significance. J. Chem. Ecol. 26, 1667–1683 (2000).CAS 

Google Scholar 
40.Ballari, S., Farji-Brener, A. G. & Tadey, M. Waste management in the leaf-cutting ant Acromyrmex lobicornis: Division of labour, aggressive behaviour, and location of external refuse dumps. J. Insect Behav. 20, 87–98 (2007).
Google Scholar 
41.Leclerc, J.-B. & Detrain, C. Impact of colony size on survival and sanitary strategies in fungus-infected ant colonies. Behav. Ecol. Sociobiol. 72, 1–10 (2018).42.Pereira, H., Jossart, M. & Detrain, C. Waste management by ants: The enhancing role of larvae. Anim. Behav. 168, 187–198 (2020).
Google Scholar 
43.Diez, L., Urbain, L., Lejeune, P. & Detrain, C. Emergency measures: Adaptive response to pathogen intrusion in the ant nest. Behav. Proc. 116, 80–86 (2015).
Google Scholar 
44.López-Riquelme, G. O. & Fanjul-Moles, M. L. The funeral ways of social insects. Social strategies for corpse disposal. Trends Entomol. 9, 71–129 (2013).45.Heinze, J. & Walter, B. Moribund ants leave their nests to die in social isolation. Curr. Biol. 20, 249–252 (2010).CAS 
PubMed 

Google Scholar 
46.Leclerc, J.-B. & Detrain, C. Loss of attraction for social cues leads to fungal-infected Myrmica rubra ants withdrawing from the nest. Anim. Behav. 129, 133–141 (2017).
Google Scholar 
47.Fouks, B. & Lattorff, H. M. G. Recognition and avoidance of contaminated flowers by foraging bumblebees (Bombus terrestris). PLoS ONE 6, e26328 (2011).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
48.Pereira, H. & Detrain, C. Pathogen avoidance and prey discrimination in ants. R. Soc. Open Sci. 7, 191705 (2020).ADS 
PubMed 
PubMed Central 

Google Scholar 
49.Pereira, H. & Detrain, C. Prophylactic avoidance of hazardous prey by the ant host Myrmica rubra. Insects 11, 444 (2020).PubMed Central 

Google Scholar 
50.Tranter, C., LeFevre, L., Evison, S. E. F. & Hughes, W. O. H. Threat detection: Contextual recognition and response to parasites by ants. ISBE 26, 396–405 (2015).
Google Scholar 
51.Marikovsky, P. I. On some features of behavior of the ants Formica rufa L. infected with fungous disease. Insectes Soc. 9, 173–179 (1962).
Google Scholar 
52.Lehue, M., Detrain, C. & Collignon, B. Nest entrances, spatial fidelity, and foraging patterns in the red ant Myrmica rubra: A field and theoretical study. Insects 11, 317 (2020).PubMed Central 

Google Scholar 
53.Diez, L., Deneubourg, J.-L., Hoebeke, L. & Detrain, C. Orientation in corpse-carrying ants: Memory or chemical cues?. Anim. Behav. 81, 1171–1176 (2011).
Google Scholar 
54.Brütsch, T., Felden, A., Reber, A. & Chapuisat, M. Ant queens (Hymenoptera: Formicidae) are attracted to fungal pathogens during the initial stage of colony founding. Myrmecol. News 20, 71–76 (2014).
Google Scholar 
55.Pontieri, L., Vojvodic, S., Graham, R., Pedersen, J. S. & Linksvayer, T. A. Ant colonies prefer infected over uninfected nest sites. PLoS ONE 9, e111961 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
56.Leclerc, J.-B., Silva, J. P. & Detrain, C. Impact of soil contamination on the growth and shape of ant nests. R. Soc. Open Sci. 5, 180267 (2018).ADS 
PubMed 
PubMed Central 

Google Scholar 
57.Loreto, R. G. & Hughes, D. P. Disease dynamics in ants. Adv. Genet. 94, 287–306 (2016) (Elsevier).CAS 
PubMed 

Google Scholar 
58.Angelone, S. & Bidochka, M. J. Diversity and abundance of entomopathogenic fungi at ant colonies. J. Invertebr. Pathol. 156, 73–76 (2018).PubMed 

Google Scholar 
59.Evans, H. C., Groden, E. & Bischoff, J. F. New fungal pathogens of the red ant, Myrmica rubra, from the UK and implications for ant invasions in the USA. J. Funbiol. 114, 451–466 (2010).CAS 

Google Scholar 
60.Bos, N., Kankaanpää-Kukkonen, V., Freitak, D., Stucki, D. & Sundström, L. Comparison of twelve ant species and their susceptibility to fungal infection. Insects 10, 271 (2019).PubMed Central 

Google Scholar 
61.Theis, F. J., Ugelvig, L. V., Marr, C. & Cremer, S. Opposing effects of allogrooming on disease transmission in ant societies. Philos. Trans. R. Soc. B Biol. Sci. 370, 20140108–20140108 (2015).
Google Scholar 
62.Okuno, M., Tsuji, K., Sato, H. & Fujisaki, K. Plasticity of grooming behavior against entomopathogenic fungus Metarhizium anisopliae in the ant Lasius japonicus. J. Ethol. 30, 23–27 (2012).
Google Scholar 
63.Pereira, R. M. & Stimac, J. L. Transmission of Beauveria bassiana within nests of Solenopsis invicta (Hymenoptera: Formicidae) in the laboratory. Environ. Entomol. 21, 1427–1432 (1992).
Google Scholar 
64.Hughes, W. O. H., Thomsen, L., Eilenberg, J. & Boomsma, J. J. Diversity of entomopathogenic fungi near leaf-cutting ant nests in a Neotropical forest, with particular reference to Metarhizium anisopliae var. anisopliae. J. Invertebr. Pathol. 85, 46–53 (2004).CAS 
PubMed 

Google Scholar 
65.Novak, S. & Cremer, S. Fungal disease dynamics in insect societies: Optimal killing rates and the ambivalent effect of high social interaction rates. J. Theor. Biol. 372, 54–64 (2015).ADS 
MathSciNet 
PubMed 
MATH 

Google Scholar 
66.Qiu, H., Lu, L., Shi, Q. & He, Y. Fungus exposed Solenopsis invicta ants benefit from grooming. J. Insect. Behav. 27, 678–691 (2014).
Google Scholar 
67.Reber, A., Purcell, J., Buechel, S. D., Buri, P. & Chapuisat, M. The expression and impact of antifungal grooming in ants. J. Evol. Biol. 24, 954–964 (2011).CAS 
PubMed 

Google Scholar 
68.Sadd, B. M. & Schmid-Hempel, P. Insect immunity shows specificity in protection upon secondary pathogen exposure. Curr. Biol. 16, 1206–1210 (2006).CAS 
PubMed 

Google Scholar 
69.Traniello, J. F. A., Rosengaus, R. B. & Savoie, K. The development of immunity in a social insect: Evidence for the group facilitation of disease resistance. Proc. Natl. Acad. Sci. 99, 6838–6842 (2002).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
70.Konrad, M. et al. Social transfer of pathogenic fungus promotes active immunisation in ant colonies. PLoS Biol. 10, e1001300 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
71.Ugelvig, L. V. & Cremer, S. Social prophylaxis: Group interaction promotes collective immunity in ant colonies. Curr. Biol. 17, 1967–1971 (2007).CAS 
PubMed 

Google Scholar 
72.Bordoni, A. et al. No evidence of queen immunisation despite transgenerational immunisation in Crematogaster scutellaris ants. J. Insect Physiol. 120, 103998 (2020).CAS 
PubMed 

Google Scholar 
73.Reber, A. & Chapuisat, M. No evidence for immune priming in ants exposed to a fungal pathogen. PLoS ONE 7, e35372 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
74.Groden, E., Drummond, F. A., Garnas, J. & Franceour, A. Distribution of an invasive ant, Myrmica rubra (Hymenoptera: Formicidae), Maine. J. Econ. Entomol. 98, 1774–1784 (2005).PubMed 

Google Scholar 
75.Radchenko, A. & Elmes, G. W. Myrmica Ants (Hymenoptera: Formicidae) of the Old World (Natura optima dux Foundation, 2010).
Google Scholar 
76.Hänel, H. The life cycle of the insect pathogenic fungus Metarhizium anisopliae in the termite Nasutitermes exitiosus. Mycopathologia 80, 137–145 (1982).
Google Scholar 
77.Leclerc, J.-B. & Detrain, C. Ants detect but do not discriminate diseased workers within their nest. Sci. Nat. 103, 1–12 (2016).78.Lacey, L. A. Manual of Techniques in Invertebrate Pathology (Academic Press imprint of Elsevier Science, 2012).
Google Scholar 
79.R Core Team. R: A Language and Environment for Statistical Computing. (2020).80.Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer-Verlag, 2016).MATH 

Google Scholar 
81.Venables, W. N. & Ripley, B. D. Modern Applied Statistics with S (Springer, 2002).MATH 

Google Scholar 
82.Bates, D., Maechler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using {lme4}. J. Stat. Softw. 67, 1–48 (2015).
Google Scholar 
83.Kuznetsova, A., Brockhoff, P. B. & Christensen, R. H. B. lmerTest package: Tests in linear mixed effects models. J. Stat. Soft. 82, 1–26. https://doi.org/10.18637/JSS.V082.I13 (2017).84.Lenth, R. Emmeans: Estimated marginal means, aka least-squares means. R-package version 1.4.8. https://CRAN.R-project.org/package=emmeans (2020).85.Hothorn, T., Bretz, F. & Westfall, P. Simultaneous inference in general parametric models. Biom. J. 50, 346–363 (2008).MathSciNet 
PubMed 
MATH 

Google Scholar 
86.Therneau, T. A Package for Survival Analysis in R. R-package version 3.2-3. https://CRAN.R-project.org/package=survival (2020) More

1.Melbourne, B. A. & Hastings, A. Highly variable spread rates in replicated biological invasions: Fundamental limits to predictability. Science 325, 1536–1539 (2009).ADS 
CAS 
PubMed 

Google Scholar 
2.Lewis, M. A., Petrovskii, S. V. & Potts, J. R. The Mathematics Behind Biological Invasions (Springer, 2016).MATH 

Google Scholar 
3.Phillips, B. L. Evolutionary processes make invasion speed difficult to predict. Biol. Invasions 17, 1949–1960 (2015).
Google Scholar 
4.Peischl, S., Kirkpatrick, M. & Excoffier, L. Expansion load and the evolutionary dynamics of a species range. Am. Nat. 185, E81–E93 (2015).PubMed 

Google Scholar 
5.Burton, O. J., Travis, J. M. J. & Phillips, B. L. Trade-offs and the evolution of life-histories during range expansion. Ecol. Lett. 13, 1210–1220 (2010).PubMed 

Google Scholar 
6.Phillips, B. L. & Perkins, T. A. Spatial sorting as the spatial analogue of natural selection. Theor. Ecol. 12, 155–163 (2019).
Google Scholar 
7.Deforet, M., Carmona-Fontaine, C., Korolev, K. S. & Xavier, J. B. Evolution at the edge of expanding populations. Am. Nat. 194, 291–305 (2019).PubMed 
PubMed Central 

Google Scholar 
8.Travis, J. M. J. & Dytham, C. Dispersal evolution during invasions. Evol. Ecol. Res. 4, 1119–1129 (2002).
Google Scholar 
9.Bouin, E. et al. Invasion fronts with variable motility: Phenotype selection, spatial sorting and wave acceleration. C. R. Math. 350, 761–766 (2012).MathSciNet 
MATH 

Google Scholar 
10.Shine, R., Brown, G. P. & Phillips, B. L. An evolutionary process that assembles phenotypes through space rather than through time. Proc. Natl. Acad. Sci. USA 108, 5708–5711 (2011).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
11.Williams, J. L., Kendall, B. E. & Levine, J. M. Rapid evolution accelerates plant population spread in fragmented experimental landscapes. Science 353, 482–485 (2016).ADS 
CAS 
PubMed 

Google Scholar 
12.Weiss-Lehman, C., Hufbauer, R. A. & Melbourne, B. A. Rapid trait evolution drives increased speed and variance in experimental range expansions. Nat. Commun. 8, 14303 (2017).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
13.Ochocki, B. M. & Miller, T. E. Rapid evolution of dispersal ability makes biological invasions faster and more variable. Nat. Commun. 8, 14315 (2017).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
14.Urban, M. C., Phillips, B. L., Skelly, D. K. & Shine, R. The cane toad’s (Chaunus [Bufo] marinus) increasing ability to invade Australia is revealed by a dynamically updated range model. Proc. R. Soc. B 274, 1413–1419 (2007).PubMed 
PubMed Central 

Google Scholar 
15.Phillips, B. L., Brown, G. P. & Shine, R. Evolutionarily accelerated invasions: The rate of dispersal evolves upwards during the range advance of cane toads. J. Evol. Biol. 23, 2595–2601 (2010).CAS 
PubMed 

Google Scholar 
16.Chuang, A. & Peterson, C. R. Expanding population edges: Theories, traits, and trade-offs. Glob. Change Biol. 22, 494–512 (2016).ADS 

Google Scholar 
17.Phillips, B. L., Brown, G. P., Travis, J. M. & Shine, R. Reid’s paradox revisited: The evolution of dispersal kernels during range expansion. Am. Nat. 172, S34–S48 (2008).PubMed 

Google Scholar 
18.Alford, R. A., Brown, G. P., Schwarzkopf, L., Phillips, B. L. & Shine, R. Comparisons through time and space suggest rapid evolution of dispersal behaviour in an invasive species. Wildl. Res. 36, 23–28 (2009).
Google Scholar 
19.Lindström, T., Brown, G. P., Sisson, S. A., Phillips, B. L. & Shine, R. Rapid shifts in dispersal behavior on an expanding range edge. Proc. Natl. Acad. Sci. USA 110, 13452–13456 (2013).ADS 
PubMed 
PubMed Central 

Google Scholar 
20.Brown, G. P., Phillips, B. L. & Shine, R. The straight and narrow path: The evolution of straight-line dispersal at a cane toad invasion front. Proc. R. Soc. B 281, 20141385 (2014).PubMed 
PubMed Central 

Google Scholar 
21.DeVore, J., Ducatez, S. & Shine, R. Spatial ecology of cane toads (Rhinella marina) in their native range: A study from French Guiana. Sci. Rep. 11, 11817 (2021).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
22.Brattstrom, B. H. Homing in the giant toad, Bufo marinus. Herpetologica 18, 176–180 (1962).
Google Scholar 
23.Zug, G. R. & Zug, P. B. The marine toad Bufo marinus: A natural history resumé of native populations. Smithson. Contrib. Zool. 284, 1–58 (1979).
Google Scholar 
24.Bayliss, P. The ecology of post-metamorphic Bufo marinus in central Amazonian savanna, Brazil. Unpublished Ph.D. thesis (The University of Queensland, 1995).25.Turvey, N. Cane Toads: A Tale of Sugar, Politics and Flawed Science (Sydney University Press, 2013).
Google Scholar 
26.Carpenter, C. C. & Gillingham, J. C. Water hole fidelity in the marine toad, Bufo marinus. J. Herpetol. 21, 158–161 (1987).
Google Scholar 
27.Ward-Fear, G., Greenlees, M. J. & Shine, R. Toads on lava: Spatial ecology and habitat use of invasive cane toads (Rhinella marina) in Hawai’i. PLoS ONE 11, e0151700 (2016).PubMed 
PubMed Central 

Google Scholar 
28.Hastings, A. Can spatial variation alone lead to selection for dispersal? Theor. Popul. Biol. 24, 244–251 (1983).MATH 

Google Scholar 
29.Möbius, W., et al. The collective effect of finite-sized inhomogeneities on the spatial spread of populations in two dimensions. Preprint at http://arxiv.org/abs/1910.05332 (2019).30.Urban, M. C., Phillips, B. L., Skelly, D. K. & Shine, R. A toad more traveled: The heterogeneous invasion dynamics of cane toads in Australia. Am. Nat. 171, E134–E148 (2008).PubMed 

Google Scholar 
31.Macgregor, L. F., Greenlees, M., de Bruyn, M. & Shine, R. An invasion in slow motion: The spread of invasive cane toads (Rhinella marina) into cooler climates in southern Australia. Biol. Invasions 23(11), 3565–3581 (2021).
Google Scholar 
32.Perkins, A. T., Phillips, B. L., Baskett, M. L. & Hastings, A. Evolution of dispersal and life history interact to drive accelerating spread of an invasive species. Ecol. Lett. 16, 1079–1087 (2013).PubMed 

Google Scholar 
33.Seabrook, W. Range expansion of the introduced cane toad Bufo marinus in New South Wales. Aust. Zool. 27, 58–62 (1991).
Google Scholar 
34.Kearney, M. R. et al. Modelling species distributions without using species distributions: The cane toad in Australia under current and future climates. Ecography 31, 423–434 (2008).
Google Scholar 
35.McCann, S. M., Kosmala, G. K., Greenlees, M. J. & Shine, R. Physiological plasticity in a successful invader: Rapid acclimation to cold occurs only in cool-climate populations of cane toads (Rhinella marina). Conserv. Physiol. 6, cox072 (2018).PubMed 
PubMed Central 

Google Scholar 
36.Schwarzkopf, L. & Alford, R. A. Nomadic movement in tropical toads. Oikos 96, 492–506 (2002).
Google Scholar 
37.Seebacher, F. & Alford, R. A. Movement and microhabitat use of a terrestrial amphibian (Bufo marinus) on a tropical island: Seasonal variation and environmental correlates. J. Herpetol. 33, 208–214 (1999).
Google Scholar 
38.Phillips, B. L., Brown, G. P., Greenlees, M., Webb, J. K. & Shine, R. Rapid expansion of the cane toad (Bufo marinus) invasion front in tropical Australia. Austral Ecol. 32, 169–176 (2007).
Google Scholar 
39.Tingley, R. & Shine, R. Desiccation risk drives the spatial ecology of an invasive anuran (Rhinella marina) in the Australian semi-desert. PLoS ONE 6, e25979 (2011).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
40.Brown, G. P., Phillips, B. L., Webb, J. K. & Shine, R. Toad on the road: Use of roads as dispersal corridors by cane toads (Bufo marinus) at an invasion front in tropical Australia. Biol. Conserv. 133, 88–94 (2006).
Google Scholar 
41.Pettit, L. J., Greenlees, M. J. & Shine, R. Is the enhanced dispersal rate seen at invasion fronts a behaviourally plastic response to encountering novel ecological conditions? Biol. Lett. 12, 20160539 (2016).PubMed 
PubMed Central 

Google Scholar 
42.Jessop, T. S. et al. Exploring mechanisms and origins of reduced dispersal in island Komodo Dragons. Proc. R. Soc. B 285, 20181829 (2018).PubMed 
PubMed Central 

Google Scholar 
43.Mayr, E. Animal Species and Evolution (Harvard University Press, 1963).
Google Scholar 
44.Duckworth, R. A. The role of behavior in evolution: A search for mechanism. Evol. Ecol. 23, 513–531 (2009).
Google Scholar 
45.Muñoz, M. M. & Losos, J. B. Thermoregulatory behavior simultaneously promotes and forestalls evolution in a tropical lizard. Am. Nat. 191, E15–E26 (2017).PubMed 

Google Scholar 
46.Carroll, S. P. et al. And the beak shall inherit–evolution in response to invasion. Ecol. Lett. 8, 944–951 (2005).PubMed 

Google Scholar 
47.Stuart, Y. E. et al. Rapid evolution of a native species following invasion by a congener. Science 346, 463–466 (2014).ADS 
CAS 
PubMed 

Google Scholar 
48.Acevedo, A. A., Lampo, M. & Cipriani, R. The cane or marine toad, Rhinella marina (Anura, Bufonidae): Two genetically and morphologically distinct species. Zootaxa 4103, 574–586 (2016).PubMed 

Google Scholar 
49.Reilly, S. M. et al. Conquering the world in leaps and bounds: Hopping locomotion in toads is actually bounding. Funct. Ecol. 29, 1308–1316 (2015).
Google Scholar 
50.Griffis-Kyle, K. L., Kyle, S. & Jungels, J. Use of breeding sites by arid-land toads in rangelands: Landscape-level factors. Southwest. Nat. 56, 251–255 (2011).
Google Scholar 
51.Sinsch, U. Movement ecology of amphibians: From individual migratory behaviour to spatially structured populations in heterogeneous landscapes. Can. J. Zool. 92, 491–502 (2014).
Google Scholar 
52.Cayuela, H. et al. Determinants and consequences of dispersal in vertebrates with complex life cycles: A review of pond-breeding amphibians. Q. Rev. Biol. 95, 1–36 (2020).
Google Scholar 
53.Child, T., Phillips, B. L., Brown, G. P. & Shine, R. The spatial ecology of cane toads (Bufo marinus) in tropical Australia: Why do metamorph toads stay near the water? Austral Ecol. 33, 630–640 (2008).
Google Scholar 
54.Pettit, L., Ducatez, S., DeVore, J. L., Ward-Fear, G. & Shine, R. Diurnal activity in cane toads (Rhinella marina) is geographically widespread. Sci. Rep. 10, 5723 (2020).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
55.Shine, R., Ward-Fear, G. & Brown, G. P. A famous failure: Why were cane toads an ineffective biocontrol in Australia? Conserv. Sci. Pract. 2, e296 (2020).
Google Scholar 
56.Shine, R., Everitt, C., Woods, D. & Pearson, D. J. An evaluation of methods used to cull invasive cane toads in tropical Australia. J. Pest Sci. 91, 1081–1091 (2018).
Google Scholar 
57.Silvester, R., Greenlees, M., Shine, R. & Oldroyd, B. Behavioural tactics used by invasive cane toads (Rhinella marina) to exploit apiaries in Australia. Austral Ecol. 44, 237–244 (2019).
Google Scholar 
58.Finnerty, P. B., Shine, R. & Brown, G. P. The costs of parasite infection: Effects of removing lungworms on performance, growth and survival of free-ranging cane toads. Funct. Ecol. 32, 402–415 (2018).
Google Scholar 
59.Pettit, L., Greenlees, M. & Shine, R. The impact of transportation and translocation on dispersal behaviour in the invasive cane toad. Oecologia 184, 411–422 (2017).ADS 
PubMed 

Google Scholar 
60.Kraus, F. Alien Reptiles and Amphibians: A Scientific Compendium and Analysis (Springer, 2008).
Google Scholar 
61.McCann, S., Greenlees, M. J. & Shine, R. On the fringe of the invasion: The ecological impact of cane toads in marginally suitable habitats. Biol. Invasions 19, 2729–2737 (2017).
Google Scholar 
62.S. Kaiser et al., unpubl. Data.63.Finnerty, P., Shine, R. & Brown, G. P. Survival of the faeces: Does a nematode lungworm adaptively manipulate the behaviour of its cane toad host? Ecol. Evol. 8, 4606–4618 (2018).PubMed 
PubMed Central 

Google Scholar 
64.Brown, G. P., Kelehear, C., Pizzatto, L. & Shine, R. The impact of lungworm parasites on rates of dispersal of their anuran host, the invasive cane toad. Biol. Invasions 18, 103–114 (2016).
Google Scholar 
65.G. Ward-Fear et al., unpubl. Data. More

1.Emerson JB, Thomas BC, Alvarez W, Banfield JF. Metagenomic analysis of a high carbon dioxide subsurface microbial community populated by chemolithoautotrophs and bacteria and archaea from candidate phyla. Environ Microbiol. 2016;18:1686–703.CAS 
PubMed 

Google Scholar 
2.Hug LA, Thomas BC, Sharon I, Brown CT, Sharma R, Hettich RL, et al. Critical biogeochemical functions in the subsurface are associated with bacteria from new phyla and little studied lineages. Environ Microbiol. 2016;18:159–73.CAS 
PubMed 

Google Scholar 
3.Anantharaman K, Brown CT, Hug LA, Sharon I, Castelle CJ, Probst AJ, et al. Thousands of microbial genomes shed light on interconnected biogeochemical processes in an aquifer system. Nat Commun. 2016;7:1–11.
Google Scholar 
4.Lu X, Seuradge BJ, Neufeld JD. Biogeography of soil Thaumarchaeota in relation to soil depth and land usage. FEMS Microbiol Ecol. 2017;93:fiw246.PubMed 

Google Scholar 
5.Cardarelli EL, Bargar JR, Francis CA. Diverse Thaumarchaeota dominate subsurface ammonia-oxidizing communities in semi-arid floodplains in the western United States. Micro Ecol. 2020;80:778–92.CAS 

Google Scholar 
6.Tolar BB, Boye K, Bobb C, Maher K, Bargar JR, Francis CA. Stability of floodplain subsurface microbial communities through seasonal hydrological and geochemical cycles. Front Earth Sci. 2020;8:338.
Google Scholar 
7.Francis CA, Roberts KJ, Beman JM, Santoro AE, Oakley BB. Ubiquity and diversity of ammonia-oxidizing archaea in water columns and sediments of the ocean. PNAS. 2005;102:14683–8.CAS 
PubMed 
PubMed Central 

Google Scholar 
8.Treusch AH, Leininger S, Kletzin A, Schuster SC, Klenk H-P, Schleper C. Novel genes for nitrite reductase and Amo-related proteins indicate a role of uncultivated mesophilic crenarchaeota in nitrogen cycling. Environ Microbiol. 2005;7:1985–95.CAS 
PubMed 

Google Scholar 
9.Leininger S, Urich T, Schloter M, Schwark L, Qi J, Nicol GW, et al. Archaea predominate among ammonia-oxidizing prokaryotes in soils. Nature. 2006;442:806–9.CAS 
PubMed 

Google Scholar 
10.Wuchter C, Abbas B, Coolen MJL, Herfort L, van Bleijswijk J, Timmers P, et al. Archaeal nitrification in the ocean. PNAS. 2006;103:12317–22.CAS 
PubMed 
PubMed Central 

Google Scholar 
11.Prosser JI, Nicol GW. Archaeal and bacterial ammonia-oxidisers in soil:the quest for niche specialisation and differentiation. Trends Microbiol. 2012;20:523–31.CAS 
PubMed 

Google Scholar 
12.Mußmann M, Brito I, Pitcher A, Damste JSS, Hatzenpichler R, Richter A, et al. Thaumarchaeotes abundant in refinery nitrifying sludges express amoA but are not obligate autotrophic ammonia oxidizers. PNAS. 2011;108:16771–6.PubMed 
PubMed Central 

Google Scholar 
13.Weber EB, Lehtovirta-Morley LE, Prosser JI, Gubry-Rangin C, Laanbroek R. Ammonia oxidation is not required for growth of Group 1.1c soil Thaumarchaeota. FEMS Microbiol Ecol. 2015;91:fiv001.PubMed 
PubMed Central 

Google Scholar 
14.Lin X, Handley KM, Gilbert JA, Kostka JE. Metabolic potential of fatty acid oxidation and anaerobic respiration by abundant members of Thaumarchaeota and Thermoplasmata in deep anoxic peat. ISME J. 2015;9:2740–4.CAS 
PubMed 
PubMed Central 

Google Scholar 
15.Kato S, Itoh T, Yuki M, Nagamori M, Ohnishi M, Uematsu K, et al. Isolation and characterization of a thermophilic sulfur- and iron-reducing thaumarchaeote from a terrestrial acidic hot spring. ISME J. 2019;13:2465–74.CAS 
PubMed 
PubMed Central 

Google Scholar 
16.Aylward FO, Santoro AE. Heterotrophic Thaumarchaea with small genomes are widespread in the dark ocean. mSystems. 2020;5:e00415–20.CAS 
PubMed 
PubMed Central 

Google Scholar 
17.Reji L, Francis CA. Metagenome-assembled genomes reveal unique metabolic adaptations of a basal marine Thaumarchaeota lineage. ISME J. 2020;14:2105–15.CAS 
PubMed 
PubMed Central 

Google Scholar 
18.Ren M, Feng X, Huang Y, Wang H, Hu Z, Clingenpeel S, et al. Phylogenomics suggests oxygen availability as a driving force in Thaumarchaeota evolution. ISME J. 2019;13:2150–61.PubMed 
PubMed Central 

Google Scholar 
19.Kerou M, Alves RJE, Schleper C. Nitrososphaerales. In: Bergeys manual of systematics of archaea and bacteria ed. Bergey’s Manual Trust (Hoboken, NJ: John Wiley & Sons). 2016. https://doi.org/10.1002/9781118960608.cbm00055.20.Qin W, Martens-Habbena W, Kobelt JN, Stahl DA. Candidatus nitrosopumilales. In: Bergeys manual of systematics of archaea and bacteria ed. Bergey’s Manual Trust (Hoboken, NJ: John Wiley & Sons). 2016. https://doi.org/10.1002/9781118960608.gbm01290.21.Prosser JI, Nicol GW. Candidatus Nitrosotaleales. In: Bergeys manual of systematics of archaea and bacteria ed. Bergey’s Manual Trust (Hoboken, NJ: John Wiley & Sons). 2016. https://doi.org/10.1002/9781118960608.obm00123.22.Gubry-Rangin C, Kratsch C, Williams TA, McHardy AC, Embley TM, Prosser JI, et al. Coupling of diversification and pH adaptation during the evolution of terrestrial Thaumarchaeota. PNAS. 2015;112:9370–5.CAS 
PubMed 
PubMed Central 

Google Scholar 
23.Nicol GW, Leininger S, Schleper C, Prosser JI. The influence of soil pH on the diversity, abundance and transcriptional activity of ammonia oxidizing archaea and bacteria. Environ Microbiol. 2008;10:2966–78.CAS 
PubMed 

Google Scholar 
24.Szukics U, Abell GCJ, Hödl V, Mitter B, Sessitsch A, Hackl E, et al. Nitrifiers and denitrifiers respond rapidly to changed moisture and increasing temperature in a pristine forest soil. FEMS Microbiol Ecol. 2010;72:395–406.CAS 
PubMed 

Google Scholar 
25.Höfferle Š, Nicol GW, Pal L, Hacin J, Prosser JI, Mandić-Mulec I. Ammonium supply rate influences archaeal and bacterial ammonia oxidizers in a wetland soil vertical profile. FEMS Microbiol Ecol. 2010;74:302–15.PubMed 

Google Scholar 
26.Tourna M, Freitag TE, Nicol GW, Prosser JI. Growth, activity and temperature responses of ammonia-oxidizing archaea and bacteria in soil microcosms. Environ Microbiol. 2008;10:1357–64.CAS 
PubMed 

Google Scholar 
27.He J-Z, Shen J-P, Zhang L-M, Zhu Y-G, Zheng Y-M, Xu M-G, et al. Quantitative analyses of the abundance and composition of ammonia-oxidizing bacteria and ammonia-oxidizing archaea of a Chinese upland red soil under long-term fertilization practices. Environ Microbiol. 2007;9:2364–74.CAS 
PubMed 

Google Scholar 
28.Marusenko Y, Bates ST, Anderson I, Johnson SL, Soule T, Garcia-Pichel F. Ammonia-oxidizing archaea and bacteria are structured by geography in biological soil crusts across North American arid lands. Ecol Process. 2013;2:9.
Google Scholar 
29.Opitz S, Küsel K, Spott O, Totsche KU, Herrmann M. Oxygen availability and distance to surface environments determine community composition and abundance of ammonia-oxidizing prokaroytes in two superimposed pristine limestone aquifers in the Hainich region, Germany. FEMS Microbiol Ecol. 2014;90:39–53.CAS 
PubMed 

Google Scholar 
30.Purkamo L, Kietäväinen R, Miettinen H, Sohlberg E, Kukkonen I, Itävaara M, et al. Diversity and functionality of archaeal, bacterial and fungal communities in deep Archaean bedrock groundwater. FEMS Microbiol Ecol. 2018;94.31.Bushnell B BBTools software package. 2014. http://bbtools.jgi.doe.gov.32.Li H. BFC:correcting Illumina sequencing errors. Bioinformatics. 2015;31:2885–7.CAS 
PubMed 
PubMed Central 

Google Scholar 
33.Li D, Liu C-M, Luo R, Sadakane K, Lam T-W. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31:1674–6.CAS 
PubMed 

Google Scholar 
34.Li D, Luo R, Liu C-M, Leung C-M, Ting H-F, Sadakane K, et al. MEGAHIT v1.0: A fast and scalable metagenome assembler driven by advanced methodologies and community practices. Methods. 2016;102:3–11.CAS 
PubMed 

Google Scholar 
35.Kang D, Li F, Kirton ES, Thomas A, Egan RS, An H, et al. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ. 2019;7:e7359.PubMed 
PubMed Central 

Google Scholar 
36.Wu Y-W, Tang Y-H, Tringe SG, Simmons BA, Singer SW. MaxBin: an automated binning method to recover individual genomes from metagenomes using an expectation-maximization algorithm. Microbiome. 2014;2:26.CAS 
PubMed 
PubMed Central 

Google Scholar 
37.Wu Y-W, Simmons BA, Singer SW. MaxBin 2.0:an automated binning algorithm to recover genomes from multiple metagenomic datasets. Bioinformatics. 2016;32:605–7.CAS 
PubMed 

Google Scholar 
38.Uritskiy GV, DiRuggiero J, Taylor J. MetaWRAP—a flexible pipeline for genome-resolved metagenomic data analysis. Microbiome. 2018;6:158.PubMed 
PubMed Central 

Google Scholar 
39.Nurk S, Bankevich A, Antipov D, Gurevich A, Korobeynikov A, Lapidus A, et al. Assembling genomes and mini-metagenomes from highly chimeric reads. In: Deng M, Jiang R, Sun F, Zhang X, editors. Research in Computational Molecular Biology (RECOMB), Lecture Notes in Computer Science, Springer; Berlin, Heidelberg. 2013;7821:158–70.40.Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.CAS 
PubMed 
PubMed Central 

Google Scholar 
41.Chaumeil P-A, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk:a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics. 2020;36:1925–7.CAS 

Google Scholar 
42.Parks DH, Chuvochina M, Chaumeil P-A, Rinke C, Mussig AJ, Hugenholtz P. A complete domain-to-species taxonomy for Bacteria and Archaea. Nat Biotechnol. 2020;38:1079–86.CAS 
PubMed 

Google Scholar 
43.Parks DH, Chuvochina M, Waite DW, Rinke C, Skarshewski A, Chaumeil P-A, et al. A standardized bacterial taxonomy based on genome phylogeny substantially revises the tree of life. Nat Biotechnol. 2018;36:996–1004.CAS 
PubMed 

Google Scholar 
44.Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.CAS 
PubMed 
PubMed Central 

Google Scholar 
45.Hyatt D, Chen G-L, LoCascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal:prokaryotic gene recognition and translation initiation site identification. BMC Bioinform. 2010;11:119.
Google Scholar 
46.Seemann T. Prokka:rapid prokaryotic genome annotation. Bioinformatics. 2014;30:2068–9.CAS 
PubMed 

Google Scholar 
47.Kanehisa M, Sato Y, Morishima K. BlastKOALA and GhostKOALA: KEGG tools for functional characterization of genome and metagenome sequences. J Mol Biol. 2016;428:726–31.CAS 
PubMed 

Google Scholar 
48.Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007;35:W182–5.PubMed 
PubMed Central 

Google Scholar 
49.Huerta-Cepas J, Forslund K, Coelho LP, Szklarczyk D, Jensen LJ, Mering von C, et al. Fast genome-wide functional annotation through orthology assignment by eggNOG-mapper. Mol Biol Evol. 2017;34:2115–22.CAS 
PubMed 
PubMed Central 

Google Scholar 
50.Huerta-Cepas J, Szklarczyk D, Heller D, Hernández-Plaza A, Forslund SK, Cook H, et al. eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 2019;47:D309–14.CAS 

Google Scholar 
51.Overbeek R, Olson R, Pusch GD, Olsen GJ, Davis JJ, Disz T, et al. The SEED and the Rapid Annotation of microbial genomes using Subsystems Technology (RAST). Nucleic Acids Res. 2013;42:D206–14.PubMed 
PubMed Central 

Google Scholar 
52.Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.CAS 
PubMed 
PubMed Central 

Google Scholar 
53.Elbourne LDH, Tetu SG, Hassan KA, Paulsen IT. TransportDB 2.0: a database for exploring membrane transporters in sequenced genomes from all domains of life. Nucleic Acids Res. 2016;45:D320–4.PubMed 
PubMed Central 

Google Scholar 
54.Nielsen H, Engelbrecht J, Brunak S, von Heijne G. Identification of prokaryotic and eukaryotic signal peptides and prediction of their cleavage sites. Protein Eng. 1997;10:1–6.CAS 
PubMed 

Google Scholar 
55.Armenteros JJA, Tsirigos KD, Sønderby CK, Petersen TN, Winther O, Brunak S, et al. SignalP 5.0 improves signal peptide predictions using deep neural networks. Nat Biotechnol. 2019;37:420–3.
Google Scholar 
56.Sonnhammer EL, Heijne G, von, Krogh A. A hidden Markov model for predicting transmembrane helices in protein sequences. Proc Int Conf Intell Syst Mol Biol. 1998;6:175–82.CAS 
PubMed 

Google Scholar 
57.Krogh A, Larsson B, Heijne G, von, Sonnhammer EL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305:567–80.CAS 
PubMed 
PubMed Central 

Google Scholar 
58.Eren AM, Esen ÖC, Quince C, Vineis JH, Morrison HG, Sogin ML, et al. Anvi’o: an advanced analysis and visualization platform for ‘omics data. PeerJ. 2015;3:e1319.PubMed 
PubMed Central 

Google Scholar 
59.Edgar RC. MUSCLE:multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.CAS 
PubMed 
PubMed Central 

Google Scholar 
60.Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses. Nucleic Acids Res. 2009;25:1972–3.
Google Scholar 
61.Nguyen L-T, Schmidt HA, Haeseler von A, Minh BQ. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Nucleic Acids Res. 2015;32:268–74.CAS 

Google Scholar 
62.Hoang DT, Chernomor O, Haeseler von A, Minh BQ. Le Sy Vinh. UFBoot2: improving the ultrafast bootstrap approximation. Nucleic Acids Res. 2017;35:518–22.
Google Scholar 
63.Kalyaanamoorthy S, Minh BQ, Wong TKF, Haeseler von A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14:587–9.CAS 
PubMed 
PubMed Central 

Google Scholar 
64.Price MN, Dehal PS, Arkin AP. FastTree 2 – approximately maximum-likelihood trees for large alignments. PLoS ONE. 2010;5:e9490.PubMed 
PubMed Central 

Google Scholar 
65.Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, et al. Fast, scalable generation of high‐quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011;7:539–9.PubMed 
PubMed Central 

Google Scholar 
66.Chen I-MA, Chu K, Palaniappan K, Ratner A, Huang J, Huntemann M, et al. The IMG/M data management and analysis system v.6.0: new tools and advanced capabilities. Nucleic Acids Res. 2021;49:D751–63.CAS 
PubMed 

Google Scholar 
67.Alves RJE, Minh BQ, Urich T, Haeseler A, Schleper C. Unifying the global phylogeny and environmental distribution of ammonia-oxidising archaea based on amoA genes. Nat Commun. 2018;9:1517.PubMed 
PubMed Central 

Google Scholar 
68.Tolar BB, Mosier AC, Lund MB, Francis CA. Nitrosarchaeum. In: Bergeys manual of systematics of archaea and bacteria ed. Bergey’s Manual Trust (Hoboken, NJ: John Wiley & Sons). 2019:1–9. https://doi.org/10.1002/9781118960608.gbm01289.69.Park S-J, Kim J-G, Jung M-Y, Kim S-J, Cha I-T, Ghai R, et al. Draft genome sequence of an ammonia-oxidizing archaeon, “Candidatus Nitrosopumilus sediminis” AR2, from Svalbard in the Arctic Circle. J Bacteriol. 2012;194:6948–9.CAS 
PubMed 
PubMed Central 

Google Scholar 
70.Kim BK, Jung M-Y, Yu DS, Park S-J, Oh TK, Rhee S-K, et al. Genome sequence of an ammonia-oxidizing soil archaeon, “Candidatus Nitrosoarchaeum koreensis” MY1. J Bacteriol. 2011;193:5539–40.CAS 
PubMed 
PubMed Central 

Google Scholar 
71.Ochsenreiter T, Selezi D, Quaiser A, Bonch-Osmolovskaya L, Schleper C. Diversity and abundance of Crenarchaeota in terrestrial habitats studied by 16S RNA surveys and real time PCR. Environ Microbiol. 2003;5:787–97.CAS 
PubMed 

Google Scholar 
72.Lehtovirta-Morley LE, Stoecker K, Vilcinskas A, Prosser JI, Prosse, Nicol GW. Cultivation of an obligate acidophilic ammonia oxidizer from a nitrifying acid soil. PNAS. 2011;108:15892–7.CAS 
PubMed 
PubMed Central 

Google Scholar 
73.Lehtovirta-Morley LE, Ross J, Hink L, Weber EB, Gubry-Rangin C, Thion C, et al. Isolation of “Candidatus Nitrosocosmicus franklandus,” a novel ureolytic soil archaeal ammonia oxidiser with tolerance to high ammonia concentration. FEMS Microbiol Ecol. 2016;92:fiw057.PubMed 
PubMed Central 

Google Scholar 
74.Könneke M, Bernhard AE, la Torre de JR, Walker CB, Waterbury JB, Stahl DA. Isolation of an autotrophic ammonia-oxidizing marine archaeon. Nature. 2005;437:543–6.PubMed 

Google Scholar 
75.Qin W, Amin SA, Martens-Habbena W, Walker CB, Urakawa H, Devol AH, et al. Marine ammonia-oxidizing archaeal isolates display obligate mixotrophy and wide ecotypic variation. PNAS. 2014;111:12504–9.CAS 
PubMed 
PubMed Central 

Google Scholar 
76.Santoro AE, Dupont CL, Richter RA, Craig MT, Carini P, McIlvin MR, et al. Genomic and proteomic characterization of “Candidatus Nitrosopelagicus brevis”: an ammonia-oxidizing archaeon from the open ocean. PNAS. 2015;112:1173–8.CAS 
PubMed 
PubMed Central 

Google Scholar 
77.Bayer B, Vojvoda J, Offre P, Alves RJE, Elisabeth NH, Garcia JA, et al. Physiological and genomic characterization of two novel marine thaumarchaeal strains indicates niche differentiation. ISME J. 2015;10:1051–63.PubMed 
PubMed Central 

Google Scholar 
78.Larentis M, Psenner R, Alfreider A. Prokaryotic community structure in deep bedrock aquifers of the Austrian Central Alps. Antonie van Leeuwenhoek. 2015;107:687–701.PubMed 

Google Scholar 
79.Lazar CS, Stoll W, Lehmann R, Herrmann M, Schwab VF, Akob DM, et al. Archaeal diversity and CO2 fixers in carbonate-/siliciclastic-rock groundwater ecosystems. Archaea. 2017;2136287.80.Sheridan PO, Raguideau S, Quince C, Holden J, Zhang L, Williams TA, et al. Gene duplication drives genome expansion in a major lineage of Thaumarchaeota. Nat Commun. 2020;11:1–12.
Google Scholar 
81.Könneke M, Schubert DM, Brown PC, Hügler M, Standfest S, Schwander T, et al. Ammonia-oxidizing archaea use the most energy-efficient aerobic pathway for CO2 fixation. PNAS. 2014;111:8239–44.PubMed 
PubMed Central 

Google Scholar 
82.Hallam SJ, Konstantinidis KT, Putnam N, Schleper C, Watanabe Y-I, Sugahara J, et al. Genomic analysis of the uncultivated marine crenarchaeote Cenarchaeum symbiosum. PNAS. 2006;103:18296–301.CAS 
PubMed 
PubMed Central 

Google Scholar 
83.Spang A, Poehlein A, Offre P, Zumbr a gel S, Haider S, Rychlik N, et al. The genome of the ammonia-oxidizing Candidatus Nitrososphaera gargensis: insights into metabolic versatility and environmental adaptations. Environ Microbiol. 2012;14:3122–45.CAS 
PubMed 

Google Scholar 
84.Kamanda Ngugi D, Blom J, Alam I, Rashid M, Ba-Alawi W, Zhang G, et al. Comparative genomics reveals adaptations of a halotolerant thaumarchaeon in the interfaces of brine pools in the Red Sea. ISME J. 2015;9:396–411.CAS 
PubMed 

Google Scholar 
85.Abby SS, Melcher M, Kerou M, Krupovic M, Stieglmeier M, Rossel C, et al. Candidatus Nitrosocaldus cavascurensis, an ammonia oxidizing, extremely thermophilic archaeon with a highly mobile genome. Front Microbiol. 2018;9:28.PubMed 
PubMed Central 

Google Scholar 
86.Tourna M, Stieglmeier M, Spang A, Konneke M, Schintlmeister A, Urich T, et al. Nitrososphaera viennensis, an ammonia oxidizing archaeon from soil. PNAS. 2011;108:8420–5.CAS 
PubMed 
PubMed Central 

Google Scholar 
87.Johnson WV, Anderson PM. Bicarbonate is a recycling substrate for cyanase. J Biol Chem. 1987;262:9021–5.CAS 
PubMed 

Google Scholar 
88.Palatinszky M, Herbold C, Jehmlich N, Pogoda M, Han P, Bergen von M, et al. Cyanate as an energy source for nitrifiers. Nature. 2015;524:105–8.CAS 
PubMed 
PubMed Central 

Google Scholar 
89.Kitzinger K, Padilla CC, Marchant HK, Hach PF, Herbold CW, Kidane AT, et al. Cyanate and urea are substrates for nitrification by Thaumarchaeota in the marine environment. Nat Microbiol. 2019;4:234–43.CAS 
PubMed 

Google Scholar 
90.Pace HC, Brenner C. The nitrilase superfamily: classification, structure and function. Genome Biol. 2001;2:REVIEWS0001. https://doi.org/10.1186/gb-2001-2-1-reviews0001.91.Ramteke PW, Maurice NG, Joseph B, Wadher BJ. Nitrile-converting enzymes: an eco-friendly tool for industrial biocatalysis. Biotechnol Appl Biochem. 2013;60:459–81.CAS 
PubMed 

Google Scholar 
92.Walker CB, la Torre de JR, Klotz MG, Urakawa H, Pinel N, Arp DJ, et al. Nitrosopumilus maritimus genome reveals unique mechanisms for nitrification and autotrophy in globally distributed marine crenarchaea. PNAS. 2010;107:8818–23.CAS 
PubMed 
PubMed Central 

Google Scholar 
93.Mosier AC, Lund MB, Francis CA. Ecophysiology of an ammonia-oxidizing archaeon adapted to low-salinity habitats. Micro Ecol. 2012;64:955–63.CAS 

Google Scholar 
94.Lebedeva EV, Hatzenpichler R, Pelletier E, Schuster N, Hauzmayer S, Bulaev A, et al. Enrichment and genome sequence of the group i.1a ammonia-oxidizing archaeon “Ca. Nitrosotenuis uzonensis” representing a clade globally distributed in thermal habitats. PLoS ONE. 2013;8:e80835.PubMed 
PubMed Central 

Google Scholar 
95.Daebeler A, Herbold C, Vierheilig J, Sedlacek CJ, Pjevac P, Albertsen M, et al. Cultivation and genomic analysis of “Candidatus Nitrosocaldus islandicus,” an obligately thermophilic, ammonia-oxidizing thaumarchaeon from a hot spring biofilm in Graendalur valley, Iceland. Front Microbiol. 2018;9:193.PubMed 
PubMed Central 

Google Scholar 
96.Beam JP, Jay ZJ, Kozubal MA, Inskeep WP. Niche specialization of novel Thaumarchaeota to oxic and hypoxic acidic geothermal springs of Yellowstone National Park. ISME J. 2014;8:938–51.CAS 
PubMed 

Google Scholar 
97.Kim J-G, Park S-J, Damste JSS, Schouten S, Rijpstra WIC, Jung M-Y, et al. Hydrogen peroxide detoxification is a key mechanism for growth of ammonia-oxidizing archaea. PNAS. 2016;113:7888–93.CAS 
PubMed 
PubMed Central 

Google Scholar 
98.Imlay JA. Cellular defenses against superoxide and hydrogen peroxide. Annu Rev Biochem. 2008;77:755–76.CAS 
PubMed 
PubMed Central 

Google Scholar 
99.Zhalnina KV, Dias R, Leonard MT, de Quadros PD, Camargo FAO, Drew JC, et al. Genome sequence of Candidatus Nitrososphaera evergladensis from group I.1b enriched from everglades soil reveals novel genomic features of the ammonia-oxidizing archaea. PLoS ONE. 2014;9:e101648.PubMed 
PubMed Central 

Google Scholar 
100.Sauder LA, Albertsen M, Engel K, Schwarz J, Nielsen PH, Wagner M, et al. Cultivation and characterization of Candidatus Nitrosocosmicus exaquare, an ammonia-oxidizing archaeon from a municipal wastewater treatment system. ISME J. 2017;11:1142–57.CAS 
PubMed 
PubMed Central 

Google Scholar 
101.Tolar BB, Powers LC, Miller WL, Wallsgrove NJ, Popp BN, Hollibaugh JT. Ammonia oxidation in the ocean can be inhibited by nanomolar concentrations of hydrogen peroxide. Front Mar Sci. 2016;3:237.
Google Scholar 
102.Bayer B, Pelikan C, Bittner MJ, Reinthaler T, Könneke M, Herndl GJ, et al. Proteomic response of three marine ammonia-oxidizing archaea to hydrogen peroxide and their metabolic interactions with a heterotrophic alphaproteobacterium. mSystems. 2019;4:e00181–19.PubMed 
PubMed Central 

Google Scholar 
103.Woodcroft BJ, Singleton CM, Boyd JA, Evans PN, Emerson JB, Zhayed AAF, et al. Genome-centric view of carbon processing in thawing permafrost. Nature. 2018;560:49–54.CAS 
PubMed 

Google Scholar 
104.Yang Y, Herbold CW, Jung M-Y, Qin W, Cai M, Du H, et al. Survival strategies of ammonia-oxidizing archaea (AOA) in a full-scale WWTP treating mixed landfill leachate containing copper ions and operating at low-intensity of aeration. Water Res. 2021;191:116798.CAS 
PubMed 

Google Scholar 
105.Greening C, Biswas A, Carere CR, Jackson CJ, Taylor MC, Stott MB, et al. Genomic and metagenomic surveys of hydrogenase distribution indicate H2 is a widely utilised energy source for microbial growth and survival. ISME J. 2016;10:761–77.CAS 
PubMed 

Google Scholar 
106.Ma K, Schicho RN, Kelly RM, Adams MW. Hydrogenase of the hyperthermophile Pyrococcus furiosus is an elemental sulfur reductase or sulfhydrogenase:evidence for a sulfur-reducing hydrogenase ancestor. PNAS. 1993;90:5341–4.CAS 
PubMed 
PubMed Central 

Google Scholar 
107.Finney AJ, Sargent F. Formate hydrogenlyase:A group 4 [NiFe]-hydrogenase in tandem with a formate dehydrogenase. Adv Micro Physiol. 2019;74:465–86.
Google Scholar 
108.Baker BJ, Saw JH, Lind AE, Lazar CS, Hinrichs KU, Teske AP, et al. Genomic inference of the metabolism of cosmopolitan subsurface archaea, Hadesarchaea. Nat Microbiol. 2016;1:1–9.
Google Scholar 
109.He Y, Li M, Perumal V, Feng X, Fang J, Xie J, et al. Genomic and enzymatic evidence for acetogenesis among multiple lineages of the archaeal phylum Bathyarchaeota widespread in marine sediments. Nat Microbiol. 2016;1:1–9.
Google Scholar 
110.Lazar CS, Baker BJ, Seitz KW, Teske AP. Genomic reconstruction of multiple lineages of uncultured benthic archaea suggests distinct biogeochemical roles and ecological niches. ISME J. 2017;11:1118–29.CAS 
PubMed 
PubMed Central 

Google Scholar 
111.Farag IF, Biddle JF, Zhao R, Martino AJ, House CH, León-Zayas RI. Metabolic potentials of archaeal lineages resolved from metagenomes of deep Costa Rica sediments. ISME J. 2020;14:1345–58.CAS 
PubMed 
PubMed Central 

Google Scholar 
112.Orsi WD, Vuillemin A, Rodriguez P, Coskun ÖK, Gomez-Saez GV, Lavik G, et al. Metabolic activity analyses demonstrate that Lokiarchaeon exhibits homoacetogenesis in sulfidic marine sediments. Nat Microbiol. 2020;5:248–55.CAS 
PubMed 

Google Scholar 
113.Adam PS, Borrel G, Gribaldo S. Evolutionary history of carbon monoxide dehydrogenase/acetyl-CoA synthase, one of the oldest enzymatic complexes. PNAS. 2018;115:E1166–73.CAS 
PubMed 
PubMed Central 

Google Scholar 
114.Köpke M, Held C, Hujer S, Liesegang H, Wiezer A, Wollherr A, et al. Clostridium ljungdahlii represents a microbial production platform based on syngas. PNAS. 2010;107:13087–92.PubMed 
PubMed Central 

Google Scholar 
115.Lazar CS, Baker BJ, Seitz KW, Hyde AS, Dick GJ, Hinrichs KU, et al. Genomic evidence for distinct carbon substrate preferences and ecological niches of Bathyarchaeota in estuarine sediments. Nucleic Acids Res. 2015;18:1200–11.
Google Scholar 
116.Debnar-Daumler C, Seubert A, Schmitt G, Heider J. Simultaneous involvement of a tungsten-containing aldehyde:ferredoxin oxidoreductase and a phenylacetaldehyde dehydrogenase in anaerobic phenylalanine metabolism. J Bacteriol. 2014;196:483–92.PubMed 
PubMed Central 

Google Scholar 
117.Kletzin A, Mukund S, Kelley-Crouse TL, Chan MK, Rees DC, Adams MW. Molecular characterization of the genes encoding the tungsten-containing aldehyde ferredoxin oxidoreductase from Pyrococcus furiosus and formaldehyde ferredoxin oxidoreductase from Thermococcus litoralis. J Bacteriol. 1995;177:4817–9.CAS 
PubMed 
PubMed Central 

Google Scholar 
118.Arndt F, Schmitt G, Winiarska A, Saft M, Seubert A, Kahnt J, et al. Characterization of an aldehyde oxidoreductase from the mesophilic bacterium Aromatoleum aromaticum ebn1, a member of a new subfamily of tungsten-containing enzymes. Front Microbiol. 2019;10. https://doi.org/10.3389/fmicb.2019.00071.119.Lloyd KG, Schreiber L, Petersen DG, Kjeldsen KU, Lever MA, Steen AD, et al. Predominant archaea in marine sediments degrade detrital proteins. Nature. 2013;496:215–8.CAS 
PubMed 

Google Scholar 
120.Dimapilis JRR. Tungsten is essential for long-term maintenance of members of candidate archaeal genus Aigarchaeota Group 4. [dissertation on the Internet]. San Bernardino, California State University; 2019. https://scholarworks.lib.csusb.edu/etd/927/.121.Anthony C. The quinoprotein dehydrogenases for methanol and glucose. Arch Biochem Biophys. 2004;428:2–9.CAS 
PubMed 

Google Scholar 
122.Jaffe AL, Castelle CJ, Dupont CL, Banfield JF. Lateral gene transfer shapes the distribution of rubisco among candidate phyla radiation bacteria and DPANN archaea. Nucleic Acids Res. 2019;36:435–46.CAS 

Google Scholar 
123.Herbold CW, Lehtovirta-Morley LE, Jung M-Y, Jehmlich N, Hausmann B, Han P, et al. Ammonia-oxidising archaea living at low pH: insights from comparative genomics. Environ Microbiol. 2017;19:4939–52.CAS 
PubMed 
PubMed Central 

Google Scholar 
124.Aono R, Sato T, Imanaka T, Atomi H. A pentose bisphosphate pathway for nucleoside degradation in Archaea. Nat Chem Biol. 2015;11:355–60.CAS 
PubMed 

Google Scholar 
125.Chadwick GL, Hemp J, Fischer WW, Orphan VJ. Convergent evolution of unusual complex I homologs with increased proton pumping capacity: energetic and ecological implications. ISME J. 2018;12:2668–80.CAS 
PubMed 
PubMed Central 

Google Scholar 
126.Cai C, Leu AO, Xie G-J, Guo J, Feng Y, Zhao J-X, et al. A methanotrophic archaeon couples anaerobic oxidation of methane to Fe(III) reduction. ISME J. 2018;12:1929–39.CAS 
PubMed 
PubMed Central 

Google Scholar 
127.Leu AO, McIlroy SJ, Ye J, Parks DH, Orphan VJ, Tyson GW. Lateral gene transfer drives metabolic flexibility in the anaerobic methane-oxidizing archaeal family Methanoperedenaceae. mBio. 2020;11:e01325–20.PubMed 
PubMed Central 

Google Scholar 
128.Zhou Z, L Y, Xu W, Pan J, Luo Z-H, Li M. Genome- and community-level interaction insights into carbon utilization and element cycling functions of Hydrothermarchaeota in hydrothermal sediment. mSystems. 2020;5:e00795–19.CAS 
PubMed 
PubMed Central 

Google Scholar 
129.Tully BJ, Graham ED, Heidelberg JF. The reconstruction of 2,631 draft metagenome-assembled genomes from the global oceans. Sci Data. 2018;5:170203.CAS 
PubMed 
PubMed Central 

Google Scholar  More

1.Wang, W., Li, Q. K. & Ma, K. P. Establishment and spatial distribution of Quercus liaotungensis Koidz. seedlings in Dongling Mountain. Acta Phytoecol. Sin. 24, 595 (2000).2.Han, H. R., He, S. Q. & Zhang, X. P. The effect of light intensity on the growth and development of Quercus liaotungensis seedlings. J. Beijing For. Univ. 22, 97–100 (2000).
Google Scholar 
3.Chen, Z. P., Wang, H. & Yuan, H. B. Studies on soil seed bank and seed fate of Quercus liaotungensis forest in the Ziwu Mountains. J. Gansu Agric. Univ. 40, 7–12 (2005).
Google Scholar 
4.Li, Y. Resource investigation and superior germplasm resources selection of woody energy plants Quercus mongolica Fisch and Quercus liaotungensis Koidz, Dissertation, Chinese Academy of Forestry, (2011).5.Yin, X., Zhou, G., Sui, X., He, Q. & Li, R. Dominant climatic factors of Quercus mongolica geographical distribution and their thresholds. Acta Ecol. Sin 33, 103–109 (2013).Article 

Google Scholar 
6.Takai, T. et al. A natural variant of NAL1, selected in high-yield rice breeding programs, pleiotropically increases photosynthesis rate. Sci. Rep. 3, 1–11 (2013).Article 

Google Scholar 
7.Yang, Y. J., Tong, Y. G., Yu, G. Y., Zhang, S. B. & Huang, W. Photosynthetic characteristics explain the high growth rate for Eucalyptus camaldulensis: Implications for breeding strategy. Ind. Crop. Prod. 124, 186–191 (2018).CAS 
Article 

Google Scholar 
8.Spyridaki, A., Psylinakis, E. & Ghanotakis, D. F. Photosystem II. In Biotechnological Applications of Photosynthetic Proteins: Biochips, Biosensors and Biodevices (ed. Giardi, M.T. & Piletska, E. V.) 11–13 (Springer, Boston, 2006).9.Dąbrowski, P. et al. Prompt chlorophyll a fluorescence as a rapid tool for diagnostic changes in PSII structure inhibited by salt stress in Perennial ryegrass. J. Photochem. Photobiol. B 157, 22–31 (2016).10.Van Rooijen, R. et al. Natural variation of YELLOW SEEDLING1 affects photosynthetic acclimation of Arabidopsis thaliana. Nat. Commun. 8, 1–9 (2017).Article 

Google Scholar 
11.Zushi, K., Kajiwara, S. & Matsuzoe, N. Chlorophyll a fluorescence OJIP transient as a tool to characterize and evaluate response to heat and chilling stress in tomato leaf and fruit. Sci. Hortic. 148, 39–46 (2012).CAS 
Article 

Google Scholar 
12.Fan, J. et al. Alleviation of cold damage to photosystem II and metabolisms by melatonin in Bermudagrass. Front. Plant Sci. 6, 925 (2015).Article 

Google Scholar 
13.Van Heerden, P., Swanepoel, J. & Krüger, G. Modulation of photosynthesis by drought in two desert scrub species exhibiting C3-mode CO2 assimilation. Environ. Exp. Bot. 61, 124–136 (2007).Article 

Google Scholar 
14.Živčák, M., Brestič, M., Olšovská, K. & Slamka, P. Performance index as a sensitive indicator of water stress in Triticum aestivum L. Plant Soil Environ. 54, 133–139 (2008).Article 

Google Scholar 
15.Kalaji, H. M., Bosa, K., Kościelniak, J. & Żuk-Gołaszewska, K. Effects of salt stress on photosystem II efficiency and CO2 assimilation of two Syrian barley landraces. Environ. Exp. Bot. 73, 64–72 (2011).CAS 
Article 

Google Scholar 
16.Singh, D. P. & Sarkar, R. K. Distinction and characterisation of salinity tolerant and sensitive rice cultivars as probed by the chlorophyll fluorescence characteristics and growth parameters. Funct. Plant Biol. 41, 727–736 (2014).CAS 
Article 

Google Scholar 
17.Song, X. L. et al. NaCl stress aggravates photoinhibition of photosystem II and photosystem I in Capsicum annuum leaves under high irradiance stress. Acta Phytoecol. Sin. 35, 681 (2011).18.Sun, Y. J., Du, Y. P. & Zhai, H. Effects of different light intensity on PSII activity and recovery of Vitis vinifera cv. cabernet sauvignon leaves under high temperature stress. Plant Physiol. J. 50, 1209–1215 (2014).
Google Scholar 
19.Chen, S., Strasser, R. J. & Qiang, S. In vivo assessment of effect of phytotoxin tenuazonic acid on PSII reaction centers. Plant Physiol. Biochem. 84, 10–21 (2014).Article 

Google Scholar 
20.Zorić, A. S. et al. Resource allocation in response to herbivory and gall formation in Linaria vulgaris. Plant Physiol. Biochem. 135, 224–232 (2019).Article 

Google Scholar 
21.Butler, W. & Kitajima, M. Fluorescence quenching in photosystem II of chloroplasts. Biochim. Biophys. Acta. 376, 116–125 (1975).CAS 
Article 

Google Scholar 
22.Baker, N. R. Chlorophyll fluorescence: A probe of photosynthesis in vivo. Annu. Rev. Plant Biol. 59, 89–113 (2008).CAS 
Article 

Google Scholar 
23.Strasser, R. J., Srivastava, A. & Tsimilli-Michael, M. Screening the vitality and photosynthetic activity of plants by fluorescence transient. In Crop Improvement for Food Security (ed. Behl, R. K., Punia, M. S. & Lather, B. P. S.) 72–115 (SSARM, Hisar, 1999).24.Appenroth, K. J., Stöckel, J., Srivastava, A. & Strasser, R. Multiple effects of chromate on the photosynthetic apparatus of Spirodela polyrhiza as probed by OJIP chlorophyll a fluorescence measurements. Environ. Pollut. 115, 49–64 (2001).CAS 
Article 

Google Scholar 
25.Stirbet, A., Lazár, D., Kromdijk, J. & Govindjee, G. Chlorophyll a fluorescence induction: Can just a one-second measurement be used to quantify abiotic stress responses?. Photosynthetica 56, 86–104. https://doi.org/10.1007/s11099-018-0770-3 (2018).CAS 
Article 

Google Scholar 
26.Tsimilli-Michael, M., Strasser, R. J. In vivo assessment of plants’ vitality: applications in detecting and evaluating the impact of mycorrhization on host plants. In Mycorrhiza: State of the Art. Genetics and Molecular Biology, Eco-Function, Biotechnology, Eco-Physiology, Structure and Systematics (ed. Varma, A.) 679–703 (Springer, Dordrecht, 2008).27.Albert, K. R., Mikkelsen, T. N., Michelsen, A., Ro-Poulsen, H. & van der Linden, L. Interactive effects of drought, elevated CO2 and warming on photosynthetic capacity and photosystem performance in temperate heath plants. J. Plant Physiol. 168, 1550–1561 (2011).CAS 
Article 

Google Scholar 
28.Chen, L. et al. Melatonin is involved in regulation of bermudagrass growth and development and response to low K+ stress. Front. Plant Sci. 8, 2038 (2017).Article 

Google Scholar 
29.Zhang, L. et al. The alleviation of heat damage to photosystem II and enzymatic antioxidants by exogenous spermidine in tall fescue. Front. Plant Sci. 8, 1747 (2017).Article 

Google Scholar 
30.Yao, X. et al. Effect of shade on leaf photosynthetic capacity, light-intercepting, electron transfer and energy distribution of soybeans. Plant Growth Regul. 83, 409–416 (2017).CAS 
Article 

Google Scholar 
31.Samborska, I. A. et al. Structural and functional disorder in the photosynthetic apparatus of radish plants under magnesium deficiency. Funct. Plant Biol. 45, 668–679 (2018).CAS 
Article 

Google Scholar 
32.dos Santos, V. A. H. F. & Ferreira, M. J. Are photosynthetic leaf traits related to the first-year growth of tropical tree seedlings? A light-induced plasticity test in a secondary forest enrichment planting. For. Ecol. Manage. 460, 7900 (2020).
Google Scholar 
33.Pavlović, I. et al. Early Brassica crops responses to salinity stress: A comparative analysis between Chinese cabbage, white cabbage, and kale. Front. Plant Sci. 10, 450 (2019).Article 

Google Scholar 
34.Xin, J., Ma, S., Li, Y., Zhao, C. & Tian, R. Pontederia cordata, an ornamental aquatic macrophyte with great potential in phytoremediation of heavy-metal-contaminated wetlands. Ecotox. Environ. Safe. 203, 111024 (2020).CAS 
Article 

Google Scholar 
35.Wang, M. X. Forest genetics and breeding (ed. Wang, M. X.) 130–137 (China Forestry Publishing House, Beijing, 2001).36.Kurjak, D. et al. Variation in the performance and thermostability of photosystem II in European beech (Fagus sylvatica L.) provenances is influenced more by acclimation than by adaptation. Eur. J. For. Res. 138, 79–92 (2019).CAS 
Article 

Google Scholar 
37.Navarro-Cerrillo, R. M. et al. Growth and physiological sapling responses of eleven Quercus ilex ecotypes under identical environmental conditions. For. Ecol. Manage. 415, 58–69 (2018).Article 

Google Scholar 
38.Guo, H., Wang, X. A., Zhu, Z. H., Wang, S. X. & Guo, J. C. Seed and microsite limitation for seedling recruitment of Quercus wutaishanica on Mt. Ziwuling, Loess Plateau, China. New For. 41, 127–137 (2011).39.Li, Z. S. et al. Tree-ring growth responses of Liaodong Oak (Quercus wutaishanica) to climate in the Beijing Dongling Mountain of China. Acta Phytoecol. Sin. 41, 11 (2021).
Google Scholar 
40.Holland, V., Koller, S. & Bruggemann, W. Insight into the photosynthetic apparatus in evergreen and deciduous European oaks during autumn senescence using OJIP fluorescence transient analysis. Plant Biol. 16, 801–808. https://doi.org/10.1111/plb.12105 (2014).CAS 
Article 
PubMed 

Google Scholar 
41.Ahammed, G. J., Xu, W., Liu, A. & Chen, S. COMT1 silencing aggravates heat stress-induced reduction in photosynthesis by decreasing chlorophyll content, photosystem II activity, and electron transport efficiency in tomato. Front. Plant Sci. 9, 998 (2018).Article 

Google Scholar 
42.Kalaji, H. M. et al. Chlorophyll a fluorescence as a tool to monitor physiological status of plants under abiotic stress conditions. Acta Physiol. Plant. 38, 102 (2016).Article 

Google Scholar 
43.Liu, J., Lu, Y., Hua, W. & Last, R. L. A new light on photosystem II maintenance in oxygenic photosynthesis. Front. Plant Sci. 10, 975 (2019).Article 

Google Scholar 
44.Shucun, S. & Lingzhi, C. Leaf growth and photosynthesis of Quercus liaotungensis in Dongling Mountain region. Acta Phytoecol. Sin. 20, 212–217 (2000).
Google Scholar 
45.Wu, A., Hammer, G. L., Doherty, A., von Caemmerer, S. & Farquhar, G. D. Quantifying impacts of enhancing photosynthesis on crop yield. Nat. Plants 5, 380–388 (2019).Article 

Google Scholar 
46.Pšidová, E. et al. Altitude of origin influences the responses of PSII photochemistry to heat waves in European beech (Fagus sylvatica L.). Environ. Exp. Bot. 152, 97–106 (2018).Article 

Google Scholar 
47.Liang, D. et al. Exogenous melatonin promotes biomass accumulation and photosynthesis of kiwifruit seedlings under drought stress. Sci. Hortic. 246, 34–43 (2019).CAS 
Article 

Google Scholar 
48.Panda, D., Ray, A. & Sarkar, R. K. Yield and photochemical activity of selected rice cultivars from Eastern India under medium depth stagnant flooding. Photosynthetica 57, 1084–1093 (2019).CAS 
Article 

Google Scholar 
49.Zhang, H. H. et al. Effects of flooding stress on the photosynthetic apparatus of leaves of two Physocarpus cultivars. J. For. Res. 29, 1049–1059. https://doi.org/10.1007/s11676-017-0496-2 (2018).ADS 
CAS 
Article 

Google Scholar 
50.Lu, W. J. Plant physiology (ed. Lu, W. J.) 88–89 (China Forestry Publishing House, Beijing, 2017).51.Xiao, C. W. & Zhou, G. S. Effect of simulated precipitation change on growth, gas exchange and chlorophyll fluorescence of Caragana intermedia in Manwusu sandland. Chin. J. Appl. Ecol. 5, 692–696 (2001).ADS 

Google Scholar  More

1.Gremillion, K. J., Barton, L. & Piperno, D. R. Particularism and the retreat from theory in the archaeology of agricultural origins. Proc. Natl. Acad. Sci. U.S.A. 111, 6171–6177 (2014).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
2.Piperno, D. R., Ranere, A. J., Dickau, R. & Aceituno, F. Niche construction and optimal foraging theory in Neotropical agricultural origins: A re-evaluation in consideration of the empirical evidence. J. Archaeol. Sci. 78, 214–220 (2017).
Google Scholar 
3.Piperno, D. R. The origins of plant cultivation and domestication in the Neotropics: A behavioral ecological perspective. In Behavioral Ecology and the Transition to Agriculture (eds Kennett, D. J. & Winterhalder, B.) 137–166 (University of California Press, 2006).
Google Scholar 
4.Zeder, M. A. Core questions in domestication research. Proc. Natl. Acad. Sci. U.S.A. 112, 3191–3198 (2015).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
5.Goldberg, A., Mychajliw, A. M. & Hadly, E. A. Post-invasion demography of prehistoric humans in South America. Nature 532, 232–235 (2016).ADS 
CAS 
PubMed 

Google Scholar 
6.Riris, P. & Arroyo-Kalin, M. Widespread population decline in South America correlates with mid-Holocene climate change. Sci. Rep. 9, 6850 (2019).ADS 
PubMed 
PubMed Central 

Google Scholar 
7.de Souza, J. G. & Riris, P. Delayed demographic transition following the adoption of cultivated plants in the eastern La Plata Basin and Atlantic coast, South America. J. Archaeol. Sci. 125, 105293 (2021).
Google Scholar 
8.Bocquet-Appel, J.-P. When the world’s population took off: The springboard of the Neolithic demographic transition. Science 333, 560–561 (2011).ADS 
CAS 
PubMed 

Google Scholar 
9.Bonomo, M., Politis, G. & Gianotti, C. Montículos, Jerarquía Social y Horticultura en Las Sociedades Indígenas Del Delta Del Río Paraná (Argentina). Latin Am. Antiq. 22, 297–333 (2011).
Google Scholar 
10.Milheira, R. G., Attorre, T. & Borges, C. Construtores de cerritos na Laguna Dos Patos, Pontal da Barra, sul do Brasil: Lugar persistente, território e ambiente construído no Holoceno recente. Latin Am. Antiq. 30, 35–54 (2019).
Google Scholar 
11.Gaspar, M. D. Considerations of the sambaquis of the Brazilian coast. Antiquity 72, 592–615 (1998).
Google Scholar 
12.De Blasis, P., Fish, S. K., Gaspar, M. D. & Fish, P. R. Some references for discussion of complexity among the Sambaqui moundbuilders from the southern shores of Brazil. Rev. Arqueol. Am. 15, 75–105 (1998).
Google Scholar 
13.Schaan, D. Long-term human induced impacts on Marajó Island Landscapes, Amazon Estuary. Diversity 2, 182–206 (2010).
Google Scholar 
14.Chanca, I. et al. Food and diet of the pre-Columbian mound builders of the Patos Lagoon region in southern Brazil with stable isotope analysis. J. Archaeol. Sci. 133, 105439 (2021).
Google Scholar 
15.Lima, T. A. Em busca dos frutos do mar: Os pescadores-coletores do litoral centro-sul do Brasil. Rev. Usp 44, 270–327 (2000).
Google Scholar 
16.Prous, A. Arqueologia brasileira (Editora Universidade de Brasília Brasília, 1991).
Google Scholar 
17.Rohr, J. A. Sítios Arqueológicos de Santa Catarina. Anais do Museu de Antropol. UFSC XVI, 77–167 (1984).
Google Scholar 
18.Schmitz, P. I. Considerações sobre a ocupação pré-histórica do litoral meridional do Brasil. Pesqui. Antropol. 63, 355–364 (2006).
Google Scholar 
19.Schmitz, P. I. Visão de conjunto dos sítios da Tapera, Armação do Sul, Laranjeiras I e II, Pântano do Sul e Cabeçudas. Pesqui. Antropol. 53, 183–190 (1996).
Google Scholar 
20.Fish, S. K., De Blasis, P., Gaspar, M. D. & Fish, P. R. Incremental events in the construction of sambaquis, southeastern Santa Catarina. Rev. Mus. Arqueol. Etnol. 1, 69–87 (2000).
Google Scholar 
21.Tenorio, M. C. Abandonment of Brazilian coastal sites: Why leave the Eden. In Explorations in American archaeology. Essays in Honor of Wesley R. Hurt (ed. Hurt, W. R.) 221–257 (University Press of America, 1998).
Google Scholar 
22.Fossile, T. et al. Pre-Columbian fisheries catch reconstruction for a subtropical estuary in South America. Fish Fish 47, 67 (2019).
Google Scholar 
23.Villagran, X. S., Klokler, D., Peixoto, S., DeBlasis, P. & Giannini, P. C. F. Building coastal landscapes: Zooarchaeology and geoarchaeology of Brazilian shell mounds. J. Island Coast. Archaeol. 6, 211–234 (2011).
Google Scholar 
24.Fish, P. R. et al. Monumental shell mounds as persistent places in southern coastal Brazil. In The Archaeology and Historical Ecology of Small Scale Economies 120–140 (2013).25.Villagran, X. S. A redefinition of waste: Deconstructing shell and fish mound formation among coastal groups of southern Brazil. J. Anthropol. Archaeol. 36, 211–227 (2014).
Google Scholar 
26.Schmitz, P. I. Acampamentos litorâneos em Içara, SC. Um exercício em padrão de assentamento. Clio 35, 99–118 (1996).
Google Scholar 
27.Kneip, A., Farias, D. & DeBlasis, P. Longa duração e territorialidade da ocupação sambaquieira na laguna de Santa Marta, Santa Catarina. Rev. Arqueol. 31, 25–51 (2018).
Google Scholar 
28.DeBlasis, P., Farias, D. S. & Kneip, A. Velhas tradições e gente nova no pedaço: Perspectivas longevas de arquitetura funerária na paisagem do litoral sul catarinense. Rev. Mus. Arqueol. Etnol. 24, 109 (2014).
Google Scholar 
29.Bandeira, D. R. Ceramistas Pre-coloniais da Baia da Babitonga, SC: Arqueologia e Etnicidade (Universidade Estadual de Campinas, 2004).
Google Scholar 
30.Colonese, A. C. et al. Long-term resilience of late Holocene coastal subsistence system in Southeastern South America. PLoS ONE 9, e93854 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
31.Pezo-Lanfranco, L. et al. Middle Holocene plant cultivation on the Atlantic Forest coast of Brazil?. R. Soc. Open Sci. 5, 180432 (2018).ADS 
PubMed 
PubMed Central 

Google Scholar 
32.Bastos, M. Q. R. et al. Isotopic evidences regarding migration at the archeological site of Praia da Tapera: New data to an old matter. J. Archaeol. Sci. Rep. 4, 588–595 (2015).
Google Scholar 
33.Bastos, M. Q. R., Lessa, A., Rodrigues-Carvalho, C., Tykot, R. H. & Santos, R. V. Carbon and nitrogen isotope analysis: Diet before and after the arrival of ceramic at Forte Marechal Luz Site. Re. Mus. Arqueol. Etnol. 24, 137–151 (2014).
Google Scholar 
34.Pezo-Lanfranco, L., DeBlasis, P. & Eggers, S. Weaning process and subadult diets in a monumental Brazilian shellmound. J. Archaeol. Sci. Rep. 22, 452–469 (2018).
Google Scholar 
35.De Masi, M. A. N. Aplicações de isótopos estáveis de O, C e N em estudos de sazonalidade, mobilidade e dieta de populações pré-históricas no sul do Brasil. Rev. Arqueol. 22, 55–76 (2009).
Google Scholar 
36.De Masi, M. Pescadores coletores da costa sul do Brasil. Pesquisas 57, 1–136 (2001).
Google Scholar 
37.Oppitz, G. et al. Pensando sobre mobilidade, dieta e mudança cultural: Análises isotópicas no sítio Armação do Sul, Florianópolis/SC. Cadernos do LEPAARQ (UFPEL) 15, 237–266 (2018).
Google Scholar 
38.Figuti, L. et al. Investigações arqueológicas e geofísicas dos sambaquis fluviais do vale do Ribeira de Iguape, Estado de São Paulo. In Museu de Arqueologia e Etnologia, USP. Relatório Final de Atividades de Projeto Temático, Processo FAPESP. 2 (2004).39.Figuti, L. Construindo o sambaqui: A ocupação e os processos de construção de sitio na bacia do Canal do Palmital, Santa Catarina. São Paulo: MAE/USP, 2009. Relatório. Processo FAPESP 08/01285-0 (2009).40.Crouch, M. S. P. Testing the Subsistence Model for the Adoption of Ceramic Technology Among Coastal Sambaqui Foragers of Southern Brazil (2013).41.Scheel-Ybert, R. & Boyadjian, C. Gardens on the coast: Considerations on food production by Brazilian shellmound builders. J. Anthropol. Archaeol. 60, 101211 (2020).
Google Scholar 
42.Pezo-Lanfranco, L. Evidence of variability in carbohydrate consumption in prehistoric fisher–hunter–gatherers of Southeastern Brazil: Spatiotemporal trends of oral health markers. Am. J. Phys. Anthropol. 167, 507–523 (2018).PubMed 

Google Scholar 
43.Boyadjian, C. H. C., Eggers, S., Reinhard, K. J. & Scheel-Ybert, R. Dieta no sambaqui Jabuticabeira-II (SC): Consumo de plantas revelado por microvestígios provenientes de cálculo dentário. Cadernos do LEPAARQ (UFPEL) 13, 131–161 (2016).
Google Scholar 
44.Merencio, F. T. & DeBlasis, P. Análises de mobilidade no litoral sul de Santa Catarina entre 2000–500 cal AP. Rev. Mus. Arqueol. Etnol. 36, 57–91 (2021).
Google Scholar 
45.Angulo, R. J., Lessa, G. C. & de Souza, M. C. A critical review of mid- to late-Holocene sea-level fluctuations on the eastern Brazilian coastline. Quat. Sci. Rev. 25, 486–506 (2006).ADS 

Google Scholar 
46.DeBlasis, P., Gaspar, M. & Kneip, A. Sambaquis from the Southern Brazilian coast: Landscape building and enduring heterarchical societies throughout the Holocene. Land 10, 757 (2021).
Google Scholar 
47.Wesolowski, V., Ferraz Mendonça de Souza, S. M., Reinhard, K. J. & Ceccantini, G. Evaluating microfossil content of dental calculus from Brazilian sambaquis. J. Archaeol. Sci. 37, 1326–1338 (2010).
Google Scholar 
48.da Rocha Bandeira, D. The use of wildlife by sambaquianos in prehistoric Babitonga Bay, North coast of Santa Catarina. Brazil. Rev. Chil. Antropol. https://doi.org/10.5354/0719-1472.2016.40613 (2015).Article 

Google Scholar 
49.Fossile, T., Ferreira, J., Bandeira, D. R., Dias-da-Silva, S. & Colonese, A. C. Integrating zooarchaeology in the conservation of coastal-marine ecosystems in Brazil. Quat. Int. https://doi.org/10.1016/j.quaint.2019.04.022 (2019).Article 

Google Scholar 
50.Ramsey, C. B. Methods for summarizing radiocarbon datasets. Radiocarbon 59, 1809–1833 (2017).CAS 

Google Scholar 
51.Crema, E. R., Bevan, A. & Shennan, S. Spatio-temporal approaches to archaeological radiocarbon dates. J. Archaeol. Sci. 87, 1–9 (2017).CAS 

Google Scholar 
52.Shennan, S. et al. Regional population collapse followed initial agriculture booms in mid-Holocene Europe. Nat. Commun. 4, 2486 (2013).ADS 
PubMed 

Google Scholar 
53.Toniolo, T. F. et al. Sea-level fall and coastal water cooling during the Late Holocene in Southeastern Brazil based on vermetid bioconstructions. Mar. Geol. 428, 106281 (2020).ADS 
CAS 

Google Scholar 
54.Cruz, F. W. et al. Orbitally driven east–west antiphasing of South American precipitation. Nat. Geosci. 2, 210–214 (2009).ADS 
CAS 

Google Scholar 
55.Carvalho do Amaral, P. G., Fonseca Giannini, P. C., Sylvestre, F. & Ruiz Pessenda, L. C. Paleoenvironmental reconstruction of a Late Quaternary lagoon system in southern Brazil (Jaguaruna region, Santa Catarina state) based on multi-proxy analysis. J. Quat. Sci. 27, 181–191 (2012).
Google Scholar 
56.Zular, A. et al. Late Holocene intensification of colds fronts in southern Brazil as indicated by dune development and provenance changes in the São Francisco do Sul coastal barrier. Mar. Geol. 335, 64–77 (2013).ADS 

Google Scholar 
57.Robinson, M. et al. Uncoupling human and climate drivers of late Holocene vegetation change in southern Brazil. Sci. Rep. 8, 7800 (2018).ADS 
PubMed 
PubMed Central 

Google Scholar 
58.Bryan, A. L. & Gruhn, R. The Sambaqui at Forte Marechal Luz, State of Santa Catarina, Brazil: Archaeological Research at Six Cave Or Rockshelter Sites in Interior Bahia, Brazil (Oregon State University, 1993).
Google Scholar 
59.Oppitz, G. Coisas que mudam: Os processos de mudança nos sítios conchíferos catarinenses e um olhar isotópico sobre o caso do sítio Armação do Sul, Florianópolis/SC (Universidade de São Paulo, 2015). https://doi.org/10.11606/D.71.2015.tde-11112015-105226.Book 

Google Scholar 
60.Garcia, A. M., Hoeinghaus, D. J., Vieira, J. P. & Winemiller, K. O. Isotopic variation of fishes in freshwater and estuarine zones of a large subtropical coastal lagoon. Estuar. Coast. Shelf Sci. 73, 399–408 (2007).ADS 

Google Scholar 
61.Stuthmann, L. E. & Castellanos-Galindo, G. A. Trophic position and isotopic niche of mangrove fish assemblages at both sides of the Isthmus of Panama. Bull. Mar. Sci. 96, 449–468 (2020).
Google Scholar 
62.Romanuk, T. N., Hayward, A. & Hutchings, J. A. Trophic level scales positively with body size in fishes: Trophic level and body size in fishes. Glob. Ecol. Biogeogr. 20, 231–240 (2011).
Google Scholar 
63.Guiry, E. Complexities of stable carbon and nitrogen isotope biogeochemistry in ancient freshwater ecosystems: Implications for the study of past subsistence and environmental change. Front. Ecol. Evol. 7, 313 (2019).
Google Scholar 
64.Gu, B. Variations and controls of nitrogen stable isotopes in particulate organic matter of lakes. Oecologia 160, 421–431 (2009).ADS 
CAS 
PubMed 

Google Scholar 
65.Kendall, C. Tracing nitrogen sources and cycling in catchments. In Isotope Tracers in Catchment Hydrology (eds Kendall, C. & McDonnell, J. J.) 519–576 (Elsevier, 1998).
Google Scholar 
66.Alves-Costa, C. P., da Fonseca, G. A. B. & Christófaro, C. Variation in the diet of the brown-nosed coati (Nasua nasua) in Southeastern Brazil. J. Mammal. 85, 478–482 (2004).
Google Scholar 
67.Beisiegel, B. M. Notes on the coati, Nasua nasua (Carnivora: Procyonidae) in an Atlantic forest area. Braz. J. Biol. 61, 689–692 (2001).CAS 
PubMed 

Google Scholar 
68.Norton, M. The chicken or the Iegue: Human–animal relationships and the Columbian exchange. Am. Hist. Rev. 120, 28–60 (2015).
Google Scholar 
69.Métraux, A. La Civilisation Matérielle Des Tribus Tupi-Guarani (Librairie Orientaliste Paul Geuthner, 1928).
Google Scholar 
70.de Azevedo, A. Q. et al. Hydrological influence on the evolution of a subtropical mangrove ecosystem during the late Holocene from Babitonga Bay, Brazil. Palaeogeogr. Palaeoclimatol. Palaeoecol. 574, 110 (2021).
Google Scholar 
71.França, M. C. et al. Late-Holocene subtropical mangrove dynamics in response to climate change during the last millennium. Holocene 29, 445–456 (2019).ADS 

Google Scholar 
72.Behling, H. & Negrelle, R. R. B. Tropical rain forest and climate dynamics of the Atlantic Lowland, Southern Brazil, during the late quaternary. Quat. Res. 56, 383–389 (2001).
Google Scholar 
73.Gaspar, M., DeBlasis, P., Fish, S. K. & Fish, P. R. Sambaquis (shell mound) societies of coastal Brazil. In Handbook of South American Archaeology (eds Silverman, H. & Isbell, W.) 319–335 (Springer, 2008).
Google Scholar 
74.Posth, C. et al. Reconstructing the deep population history of central and south America. Cell 175, 1185-1197.e22 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
75.Fidalgo, D., Hubbe, M. & Wesolowski, V. Population history of Brazilian south and southeast shellmound builders inferred through dental morphology. Am. J. Phys. Anthropol. https://doi.org/10.1002/ajpa.24342 (2021).Article 
PubMed 

Google Scholar 
76.Hubbe, M., Okumura, M., Bernardo, D. V. & Neves, W. A. Cranial morphological diversity of early, middle, and late Holocene Brazilian groups: Implications for human dispersion in Brazil. Am. J. Phys. Anthropol. 155, 546–558 (2014).PubMed 

Google Scholar 
77.Morgan, C. Is it intensification yet? Current archaeological perspectives on the evolution of hunter–gatherer economies. J. Archaeol. Res. 23, 163–213 (2015).
Google Scholar 
78.Zeder, M. A. The broad spectrum revolution at 40: Resource diversity, intensification, and an alternative to optimal foraging explanations. J. Anthropol. Archaeol. 31, 241–264 (2012).
Google Scholar 
79.Schmitz, P. E., Verardi, I., de Masi, M. A., Rogge, J. H. & Jacobus, A. L. Escavações Arqueológicas do Pe. João Alfredo Rohr. O Sítio da Praia das Laranjeiras II. Uma Aldeia de Tradição Ceramista Itararé. Pesqui. Antropol. 49, 181 (1993).
Google Scholar 
80.Tiburtius, G. & Bigarella, I. K. Nota sobre os anzóis de osso da jazida páleo-etnográfica de Itacoara, Santa Catarina. Rev. Mus. Paul. Nova Série 7, 381–387 (1953).
Google Scholar 
81.Lessa, A. & Medeiros, J. C. D. Preliminary thoughts about the occurence of violence among the Brazilian shellmound builders: Analysis of the skeletons from Cabeçuda (Santa Catarina) and Arapuan (Rio de Janeiro) sites. Rev. Mus. Arqueol. Etnol. 11, 77 (2001).
Google Scholar 
82.Lessa, A. Reflexões preliminares sobre paleoepidemiologia da violência em grupos ceramistas litorâneos: (I) Sítio Praia da Tapera—SC. Rev. Mus. Arqueol. Etnol. https://doi.org/10.11606/issn.2448-1750.revmae.2001.109411 (2006).Article 

Google Scholar 
83.García-Escárzaga, A. & Gutiérrez-Zugasti, I. The role of shellfish in human subsistence during the Mesolithic of Atlantic Europe: An approach from meat yield estimations. Quat. Int. 584, 9–19 (2021).
Google Scholar 
84.Erlandson, J. M. The role of shellfish in prehistoric economies: A protein perspective. Am. Antiq. 53, 102–109 (1988).
Google Scholar 
85.Bandeira, D. R. Ceramistas Pré-coloniais da Baía da Babitonga—Arqueologia e Etnicidade (Universidade Estadual de Campinas, 2004).
Google Scholar 
86.Gilson, S.-P. & Lessa, A. Arqueozoologia do sítio Rio do Meio (SC). rsab 34, 217–248 (2021).
Google Scholar 
87.Gilson, S.-P. & Lessa, A. Capture, processing and utilization of sharks in archaeological context: Its importance among fisher–hunter–gatherers from southern Brazil. J. Archaeol. Sci. Rep. 35, 102693 (2021).
Google Scholar 
88.Hayden, B. Nimrods, piscators, pluckers, and planters: The emergence of food production. J. Anthropol. Archaeol. 9, 31–69 (1990).
Google Scholar 
89.Figuti, L. & Klokler, D. Resultados preliminares dos vestígios zooarqueológicos do sambaqui Espinheiros II (Joinville, SC). Rev. Mus. Arqueol. Etnol. 6, 169–187 (1996).
Google Scholar 
90.Benz, D. M. Levantamento preliminar de algumas espécies de vertebrados pretéritos do sítio arqueológico Ilha dos Espinheiros II Joinville—SC (Universidade da Região de Joinville, 2000).
Google Scholar 
91.Gilson, S.-P. & Lessa, A. Ocupação tardia do litoral norte e central catarinense por grupos pescadores-caçadores-coletores. rsab 33, 55–77 (2020).
Google Scholar 
92.Silva, S. B., Schmitz, P. I., Rogge, J. H., de Masi, M. A. N. & Jacobus, A. L. Escavações arqueológicas do pe. João Alfredo Rohr, S. J. o sítio arqueológico da Praia da Tapera: Um assentamento Itararé e Tupiguarani (Pesquisas, 1990).
Google Scholar 
93.Cardoso, J. M. O sítio Costeiro Galheta IV: Uma Perspectiva Zooarqueológica (Museu de Arqueologia e Etnologia, 2019). https://doi.org/10.11606/d.71.2019.tde-27112018-142710.Book 

Google Scholar 
94.Allen, M. W., Bettinger, R. L., Codding, B. F., Jones, T. L. & Schwitalla, A. W. Resource scarcity drives lethal aggression among prehistoric hunter-gatherers in central California. Proc. Natl. Acad. Sci. U.S.A. 113, 12120–12125 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
95.Kelly, R. L. The Lifeways of Hunter–Gatherers: The Foraging Spectrum (Cambridge University Press, 2013).
Google Scholar 
96.Hansel, F. A. & Evershed, R. P. Formation of dihydroxy acids from Z-monounsaturated alkenoic acids and their use as biomarkers for the processing of marine commodities in archaeological pottery vessels. Tetrahedron Lett. 50, 5562–5564 (2009).CAS 

Google Scholar 
97.Hansel, F. A., Bull, I. D. & Evershed, R. P. Gas chromatographic mass spectrometric detection of dihydroxy fatty acids preserved in the ‘bound’ phase of organic residues of archaeological pottery vessels. Rapid Commun. Mass Spectrom. 25, 1893–1898 (2011).ADS 
CAS 
PubMed 

Google Scholar 
98.Hansel, F. A. & Schmitz, P. I. Classificação e interpretação dos resíduos orgânicos preservados em fragmentos de cerâmica arqueológica por cromatografia gasosa e cromatografia gasosa-espectrometria de massas. Pesqui. Antropol. 63, 81–112 (2006).
Google Scholar 
99.Bueno, L., & Gilson, S. Brazilian Radiocarbon Database. Brazilian Radiocarbon Database. https://brc14database.com.br/?page_id=32 (2021).100.Prous, A. As esculturas de pedra (zoólitos) e de osso dos sambaquis do Brasil meridional e do Uruguay. Rev. Mem. 5, 197–217 (2018).
Google Scholar 
101.Prous, A. Les sculptures préhistoriques du Sud-Brésilien. bspf 71, 210–217 (1974).
Google Scholar 
102.Ramsey, C. B. Bayesian analysis of radiocarbon dates. Radiocarbon 51, 337–360 (2009).CAS 

Google Scholar 
103.Carleton, W. C. & Groucutt, H. S. Sum things are not what they seem: Problems with point-wise interpretations and quantitative analyses of proxies based on aggregated radiocarbon dates. Holocene 31, 630–643 (2021).ADS 

Google Scholar 
104.Hogg, A. G. et al. SHCal20 Southern Hemisphere Calibration, 0–55,000 Years cal BP. Radiocarbon 62, 759–778 (2020).CAS 

Google Scholar 
105.Heaton, T. J. et al. Marine20—The marine radiocarbon age calibration curve (0–55,000 cal BP). Radiocarbon 62, 779–820 (2020).CAS 

Google Scholar 
106.Angulo, R. J., de Souza, M. C., Reimer, P. J. & Sasaoka, S. K. Reservoir effect of the southern and southeastern Brazilian coast. Radiocarbon 47, 67–73 (2008).
Google Scholar 
107.Alves, E. et al. Radiocarbon reservoir corrections on the Brazilian coast from pre-bomb marine shells. Quat. Geochronol. 29, 30–35 (2015).
Google Scholar 
108.De Masi, M. A. N. Prehistoric Hunter–Gatherer Mobility on the Southern Brazilian Coast: Santa Catarina Island (Unpublished PhD dissertation. Stanford University, 1999).109.DeNiro, M. J. Postmortem preservation and alteration of in vivo bone collagen isotope ratios in relation to palaeodietary reconstruction. Nature 317, 806 (1985).ADS 
CAS 

Google Scholar 
110.Ambrose, S. H. Preparation and characterization of bone and tooth collagen for isotopic analysis. J. Archaeol. Sci. 17, 431–451 (1990).
Google Scholar 
111.van Klinken, G. J. Bone collagen quality indicators for palaeodietary and radiocarbon measurements. J. Archaeol. Sci. 26, 687–695 (1999).
Google Scholar 
112.Szpak, P., Buckley, M., Darwent, C. M. & Richards, M. P. Long-term ecological changes in marine mammals driven by recent warming in northwestern Alaska. Glob. Change Biol. 24, 490–503 (2018).ADS 

Google Scholar 
113.Garcia, A. M., Vieira, J. P. & Winemiller, K. O. Effects of 1997–1998 El Niño on the dynamics of the shallow-water fish assemblage of the Patos Lagoon Estuary (Brazil). Estuar. Coast. Shelf Sci. 57, 489–500 (2003).ADS 

Google Scholar 
114.Wiley, A. E. et al. Millennial-scale isotope records from a wide-ranging predator show evidence of recent human impact to oceanic food webs. Proc. Natl. Acad. Sci. U.S.A. 110, 8972–8977 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
115.Buckley, M., Collins, M., Thomas-Oates, J. & Wilson, J. C. Species identification by analysis of bone collagen using matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry. Rapid Commun. Mass Spectrom. 23, 3843–3854 (2009).ADS 
CAS 
PubMed 

Google Scholar 
116.Strohalm, M., Hassman, M., Kosata, B. & Kodícek, M. mMass data miner: An open source alternative for mass spectrometric data analysis. Rapid Commun. Mass Spectrom. 22, 905–908 (2008).ADS 
PubMed 

Google Scholar 
117.Buckley, M. et al. Species identification of archaeological marine mammals using collagen fingerprinting. J. Archaeol. Sci. 41, 631–641 (2014).CAS 

Google Scholar 
118.Kirby, D. P., Buckley, M., Promise, E., Trauger, S. A. & Holdcraft, T. R. Identification of collagen-based materials in cultural heritage. Analyst 138, 4849–4858 (2013).ADS 
CAS 
PubMed 

Google Scholar 
119.Brown, T. A., Nelson, E. E., Vogel, S. J. & Southon, J. R. Improved collagen extraction by modified longin method. Radiocarbon 30, 171–177 (1988).CAS 

Google Scholar 
120.Craig, O. E. et al. Stable isotope analysis of Late Upper Palaeolithic human and faunal remains from Grotta del Romito (Cosenza), Italy. J. Archaeol. Sci. 37, 2504–2512 (2010).
Google Scholar 
121.Kragten, J. Calculating standard deviations and confidence intervals with a universally applicable spreadsheet technique. Analyst 119, 2161–2166 (1994).ADS 
CAS 

Google Scholar 
122.Good Practice Guide for Isotope Ratio Mass Spectrometry (FIRMS, 2018).123.Guiry, E. J. & Szpak, P. Improved quality control criteria for stable carbon and nitrogen isotope measurements of ancient bone collagen. J. Archaeol. Sci. 132, 105416 (2021).CAS 

Google Scholar 
124.Phillips, D. L. et al. Best practices for use of stable isotope mixing models in food-web studies. Can. J. Zool. 92, 823–835 (2014).
Google Scholar 
125.Fernandes, R., Grootes, P., Nadeau, M.-J. & Nehlich, O. Quantitative diet reconstruction of a Neolithic population using a Bayesian mixing model (FRUITS): The case study of Ostorf (Germany). Am. J. Phys. Anthropol. https://doi.org/10.1002/ajpa.22788 (2015).Article 
PubMed 

Google Scholar 
126.Jim, S., Jones, V., Ambrose, S. H. & Evershed, R. P. Quantifying dietary macronutrient sources of carbon for bone collagen biosynthesis using natural abundance stable carbon isotope analysis. Br. J. Nutr. 95, 1055–1062 (2006).CAS 
PubMed 

Google Scholar 
127.Colonese, A. C. et al. Stable isotope evidence for dietary diversification in the pre-Columbian Amazon. Sci. Rep. 10, 1–11 (2020).
Google Scholar 
128.Fernandes, R., Nadeau, M.-J. & Grootes, P. M. Macronutrient-based model for dietary carbon routing in bone collagen and bioapatite. Archaeol. Anthropol. Sci. 4, 291–301 (2012).
Google Scholar 
129.Galetti, M., Rodarte, R. R., Neves, C. L., Moreira, M. & Costa-Pereira, R. Trophic niche differentiation in rodents and marsupials revealed by stable isotopes. PLoS ONE 11, e0152494 (2016).PubMed 
PubMed Central 

Google Scholar 
130.Hellevang, H. & Aagaard, P. Constraints on natural global atmospheric CO2 fluxes from 1860 to 2010 using a simplified explicit forward model. Sci. Rep. 5, 17352 (2015).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
131.Fernandes, R. A Simple(R) model to predict the source of dietary carbon in individual consumers. Archaeometry 58, 500–512 (2016).CAS 

Google Scholar 
132.Fernandes, R., Millard, A. R., Brabec, M., Nadeau, M.-J. & Grootes, P. Food reconstruction using isotopic transferred signals (FRUITS): A Bayesian model for diet reconstruction. PLoS ONE 9, e87436 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
133.Wickham, H. ggplot2: Elegant graphics for data analysis (Springer, 2016).MATH 

Google Scholar 
134.Mcgill, R., Tukey, J. W. & Larsen, W. A. Variations of box plots. Am. Stat. 32, 12–16 (1978).
Google Scholar 
135.Kwak, S. G. & Kim, J. H. Central limit theorem: The cornerstone of modern statistics. Korean J. Anesthesiol. 70, 144–156 (2017).PubMed 
PubMed Central 

Google Scholar  More

Our novel approach to bioacoustic source separation involves an end-to-end pipeline consisting of multiple discrete steps, including (1) synthesizing a dataset, (2) developing and training a separator network to disentangle the input mixture, and (3) constructing and training a classifier model to employ as a downstream evaluation task. This workflow requires few hyperparameter modifications to account for unique vocal behavior across different biological taxa but is otherwise general and makes no species-level assumptions about the spectrotemporal structure of the source calls. We develop a complete framework for bioacoustic source separation in a permutation-invariant mode using overlapping waveforms drawn from the same class of signals. We apply BioCPPNet to macaques, dolphins, and Egyptian fruit bats, and we consider two or three concurrent “speakers”. Note that we henceforth refer to non-human animal signalers as “speakers” for consistency with the human speech separation literature2. We address both the closed speaker regime in which the training and evaluation data subsets contain calls produced by individuals drawn from the same distribution as well the open speaker regime in which the model is tested on calls generated by individuals not present in the training dataset.Bioacoustic dataWe investigate a set of species with dissimilar vocal behaviors in terms of spectral and temporal properties. We apply BioCPPNet to a macaque coo call dataset47 consisting of 7285 coos produced by 8 unique individuals; a bottlenose dolphin signature whistle dataset26 comprised of 400 signature whistles generated by 20 individuals, of which we randomly select 8 for the purposes of this study; and an Egyptian fruit bat vocalization dataset48 containing a heterogeneous distribution of individuals, call types, and call contexts. In the case of the bat dataset, we extract the data (31399 calls) corresponding to the 15 most heavily represented individual bats, reserving 12 individuals (27586 calls) to address the closed speaker regime and the remaining 3 individuals (3813 calls) to evaluate model performance in the open speaker scenario.DatasetsThe mixture dataset is generated from a species-specific corpus of bioacoustic recordings containing signals annotated according to the known identity of the signaller. Motivated by WSJ0-2mix2, a preeminent reference dataset used for human single-channel acoustic source separation, we adopt a similar approach of constructing bioacoustic datasets by temporally overlapping and summing ground truth individual-specific signals to enable supervised training of our model. For macaques and dolphins, the mixture waveforms contain discrete source calls that overlap in the time domain, by design. For bats, mixtures are constructed by adding signal streams, each of which may exhibit one or more temporally separated sequential vocal elements. In all cases, the mixtures operate under the assumption that, without loss of generality, the constituent sources vary in the degree of spectral overlap due to differential spectrotemporal properties of sources, in accordance with the DUET principle (i.e, the mixtures contain approximately disjoint sources that rarely coincide in dominant frequency)19. The resultant dataset consists of an input array of the composite mixture waveforms, a target array containing the separated ground truth waveforms corresponding to the respective mixtures, and a class label array denoting the identities of the vocalizing animals responsible for generating the signals. In the case of macaques, we here consider closed speaker set mixtures of two and three simultaneous speakers, but our method is functionally not limited in the number of sources (N) it can handle. For dolphins, we consider the closed speaker regime with two overlapping calls, and for bats, we consider the closed and open speaker scenarios with two sources.We first extract the labeled waveforms either by truncating or zero-padding the waveforms to ensure that all the samples are of fixed duration. We select the number of frames either by computing the mean plus three-sigma of the durations of the calls contained in the corpus from which we draw the signals, by selecting the maximum duration of all calls, or by choosing a fixed value. For macaques, dolphins, and bats, we use 23156 frames (0.95s), 290680 frames (3.03s), and 250000 frames (1.0s), respectively. We then randomly select vocalizations from N different speakers drawn from the distribution of individuals used in the study (8 macaques, 8 dolphins, 12 bats for the closed speaker regime, 3 bats for the open speaker regime) and mix them additively, ensuring to randomly shift the overlaps to simulate a more plausible scenario and to provide for asynchronicity of start times, an important acoustic cue that has been suggested as a mechanism with which the animal brain can solve the CPP1. Despite higher computational and memory costs, we opt to use native sampling rates, since certain animal vocalizations may reach frequencies near the native Nyquist frequency. With this in mind, however, our method does provide for resampling when the vocalizations of the particular species of interest are amenable to downsampling. Explicitly, for the three species we consider including macaques, dolphins, and bats, we use sampling rates of 24414 Hz, 96 kHz, and 250 kHz, respectively. For the closed speaker regime, the training and evaluation subsets contain calls produced by the same distribution of individuals to ensure a closed speaker set. We segment the original nonoverlapping vocalizations into 80/20 training/validation subsets. We generate the mixture training waveforms using 80% of the calls, and we construct the mixture validation subset using the remaining 20% of calls held out from the training data. In the case of overlapping bat calls (for which the corpus of bioacoustic recordings contains (mathscr {O}(10text { hours})) of data as opposed to (mathscr {O}(10^{-1}text { hours})) for macaques and dolphins), we also address the open speaker source separation problem by constructing a further testing data subset of mixtures of calls of additional vocalizers not contained in the training distribution. For macaques, we construct a training data subset comprised of 12k samples and a validation subset with 3k samples, all of which contain calls drawn from 8 animals. For dolphins, we randomly select 8 individuals and construct training/validations subsets with 8k and 2k samples, respectively. For bats, we select 15 individuals, randomly reserving 12 for the closed speaker problem and the remaining 3 for the open speaker situation. We train the bat separator model on 24k mixtures. We evaluate performance in both the closed and open speaker scenarios using data subsets consisting of 6k mixtures containing unseen vocalizations produced by the appropriate distribution of individuals according to the regime under consideration. We repeat the bat training using a larger mixture dataset (denoted by +) containing 72k samples. We here report validation metrics to ensure that we are evaluating model performance on unseen mixtures of unseen calls in the closed speaker regime and on unseen mixtures of unseen calls of unseen individuals in the open speaker regime.For the downstream classification task, we extract vocalizations annotated according to the individual identity, and we segment the calls into an 80/20 training/testing split to ensure that we are evaluating model performance on unseen calls. For both the training and evaluation data subsets, we employ an augmentation scheme in which we apply random temporal shifts to call onsets to better reflect more plausible real-world scenarios.Classification modelsIn an effort to provide a more physically interpretable evaluation metric to supplement the commonly-implemented SI-SDR used in human speech separation studies, we develop CNN-based classifier models to label the individual identity of the separated vocalizations as a downstream task. This requires training classification networks to predict the speaker class label of the original unmixed waveforms. For each species we consider, we design and train custom simple and lightweight CNN-based architectures largely motivated by previous work24, tailored to accommodate the unique vocal behavior of the given species.The first layer in the model is an optional high pass filter constructed using a nontrainable 1D convolution (Conv1D) layer with frozen weights determined by a windowed sinc function49,50 to eliminate low-frequency background noise. We omit this computationally intensive layer for macaques and Egyptian fruit bats, but we implement a high pass filter for the dolphin dataset, selecting an arbitrary cutoff frequency of 4.7 kHz and transition bandwidth 0.08 to remove background without impinging on the region of support for dolphin whistles. After the optional filter is an encoder layer to compute on-the-fly feature extraction. We experimented with a fully learnable free Conv1D filterbank, a spectrogram, and a log-magnitude spectrogram and observed optimal performance using a non-decibel (dB)-scaled STFT layer computed with a nfft window width, a hop window shift, and a Hann window where nfft and hop are species-dependent variables. For macaques, we select nfft=1024 and hop=64 corresponding to temporal scales on the order of 40ms and frequency resolutions on the order of 20 Hz. We choose nfft=1024 and hop=256 for dolphins and nfft=2048 and hop=512 for bats, corresponding to temporal resolutions of ~ 10 ms and ~ 8 ms and frequency resolutions of ~ 90 Hz and ~ 120 Hz, respectively.Following the built-in feature engineering, the architecture includes 4 convolutional blocks, which consist of two sequential 2D convolution (Conv2D) layers with leaky ReLU activation and a max pooling layer with pool size 4. Next is a dense fully connected layer with leaky ReLU activation followed by another linear layer with log softmax activation to output the V log probabilities (i.e. confidences) where V is the number of individual vocalizers used in the study (8, 8, 12 for macaques, dolphins, and bats, respectively). We also include dropout regularization with p=0.25 for the macaque classifier and p=0.5 for the dolphin and bat classifiers to address potential overfitting. With these architectures, the macaque, dolphin, and bat classifier models have 230k, 279k, and 247k trainable parameters, respectively.For all species, we minimize the negative log-likelihood objective loss function using the Adam optimizer51 with learning rate lr = 3e−4. For macaques, dolphins, and bats, respectively, we train for 100, 50, and 100 epochs with batch sizes 32, 8, and 8. We serialize the model after each epoch and select the top-performing models. We opt not to carry out hyperparameter optimization since the classification task is of secondary importance and is used solely as a downstream task.Figure 1(a) Schematic overview of the BioCPPNet pipeline. Source vocalization waveforms are overlapped in time and mixed additively. BioCPPNet operates on the mixture waveform, yielding predictions for the separated waveforms, which are compared to the source ground truths, up to a permutation. The estimated waveforms are classified by the identity classification model24 (ID) to compute the downstream classification accuracy metric. (b) Block diagram of the BioCPPNet architecture. The input mixture waveform is transformed to a learnable or handcrafted representation (Rep), which then passes through a 2-dimensional U-Net52 composed of a contracting encoder path and an expanding decoder path with skip connnections. The encoder path consists of sequential downsampling convolutional blocks, each of which is constructed using two convolutional layers (Conv2D) with leaky ReLU activation and batch normalization (BatchNorm) followed by a max pooling. The decoder path employs upsampling convolutional blocks, consisting of an upsampling and skip connection concatenation followed again by the Conv2D layers with leaky ReLU and BatchNorm. The U-Net predicts masks (Mask 0 and Mask 1), the number of which is determined by the number of sources (N), that are multiplicatively applied to the original mixture representation. The predicted time-frequency representations of the separated waveforms are inverted with learnable or handcrafted inverse transforms (iRep) to output raw waveforms. All schematic diagrams were created using Affinity Designer (version 1.8.1) https://affinity.serif.com/en-us/designer/.Full size imageSeparation modelsBioCPPNet (Fig. 1) is a lightweight and modular architecture with a modifiable representation encoder, a 2D U-Net core, and an inverse transform decoder, which acts directly on raw audio via on-the-fly learnable or handcrafted transforms. The structure of the network is designed to provide for extensive experimentation, optimization, and enhancement across a range of species with variable vocal behavior. We construct and train a separation model for each species and each number N of sources contained in the input mixture.Figure 2Schematic diagram demonstrating the application of BioCPPNet to dolphin signature whistles using handcrafted STFT-based encoders and decoders. The source waveforms produced by N speakers of unique identity (e.g. T. truncatus 0 and T. truncatus 1) are overlapped in time, summed, and transformed to time-frequency space using an STFT layer, resulting in the mixture time-frequency representation (Mixture TFR). The U-Net predicts masks (Mask 0 and Mask 1) that are applied to the mixture representation. The separated spectrogram estimations (TFR 0 and TFR 1) are inverted using an iSTFT layer to yield the model’s predictions for the separated raw waveforms, which are compared to the ground truth waveforms and classified according to predicted identity using the classification model.Full size imageModel architectureAs with the classifier model, the network’s encoder consists of a feature engineering block, the initial layer of which is an optional high pass filter. This is followed by the representation transform, which includes several options including the Conv1D free encoder, the STFT filterbank, and the log-magnitude (dB) STFT filterbank. We choose the same kernel size (nfft) and stride (hop) parameters defined in the classifier model. Sequentially following the feature extraction encoder is a 2D U-Net core. This architecture consists of B (4 for macaques, 3 for dolphins, and 4 for bats) downsampling convolutional blocks, a middle convolutional block, and B upsampling convolutional blocks. The downsampling blocks consist of two 2D convolutional layers with filter number that increases with model depth with leaky ReLU activation followed by a max pooling with pool size 2, 6, and 3 for macaques, dolphins, and bats. The middle block contains two 2D convolutional layers with leaky ReLU activation. The upsampling blocks include an upsampling using the bilinear algorithm and a scale factor corresponding to the pool size used during downsampling, followed by skip connections in which the corresponding levels of the contracting and expanding paths are concatenated before passing through two 2D convolutional layers with leaky ReLU activation. All convolutional layers in the downsampling, middle, and upsampling blocks include batch normalization after the activation function to stabilize and expedite training and to promote regularization. Though our default implementation is phase-unaware, we also offer the option for a parallel U-Net pathway working directly on phase information, which has been shown to improve performance in other applications53,54,55. The final layer in the U-Net core is a 2D convolutional layer with N channels, which are then split prior to entering the inverse transform decoder. For the inverse transform, we again provide numerous choices including a free filterbank decoder based on a 1D convolutional transpose (ConvTranspose1D) layer, an iSTFT layer, an iSTFT layer accepting dB-scaled inputs, and a multi-head convolutional neural network (MCNN) for fast spectrogram inversion56. In detail, the U-Net returns N masks that are then multiplied by the original encoded representation of the mixture waveform. The separated representations are then passed into the inverse transform layer in order to yield the raw waveforms corresponding to the model’s predictions for the separated vocalizations. We initialize all trainable weights using the Xavier uniform initialization. In the case of macaques, we experiment across all combinations of representation encoders and inverse transform decoders, and we find optimal performance using the handcrafted non-dB STFT/iSTFT layers operating in the time-frequency domain. Since the model with the fully learnable Conv1D-based encoder/decoder uniquely operates in the time domain, we report evaluation metrics for this model, as well. For dolphins and bats, we here report metrics using exclusively the non-dB STFT/iSTFT technique.BioCPPNet (Fig. 1) is designed as a lightweight fully convolutional model in order to efficiently process large amounts of bioacoustic data sampled at high sampling rates while simultaneously minimizing computational costs and limitations and the likelihood of overfitting. For the macaque separators, the networks consist of 1.2M parameters (for the STFT, iSTFT combination), 2.5M parameters (for the STFT, iSTFT combination with parallel phase pathway), or 2.8M parameters (for the Conv1D free filterbanks). For the dolphin separator (Fig. 2), the model has 304k parameters, while the bat separator model has 1.2M parameters. This is to be contrasted with the comparatively heavyweight default implementations of models commonly used in human speech separation problems, such as Conv-TasNet3, which has 5.1M parameters; DPTNet4 with 2.7M parameters; or Wavesplit5 with 29M parameters. Regardless of the lower complexity of BioCPPNet, the model achieves comparable performance or even outperforms reference human speech separator models while still being lightweight enough to train on a single NVIDIA P100 GPU.Model training objectiveThe model training objective aims to optimize the reconstruction of separated waveforms from the aggregated composite input signal. We adopt a permutation-invariant training (PIT)57 scheme in which the model’s predicted outputs are compared with the ground truth sources by searching over the space of permutations of source orderings. This fundamental property of our training objective reflects that the order of estimations and their corresponding labels from a mixture waveform is not expressly germane to the task of acoustic source separation, i.e. separation is a set prediction problem independent of speaker identity ordering5.Source separation involves training a separator model f to reconstruct the source single-channel waveforms given a mixture (x=sum _{i=1}^N s^i) of N sources, where each source signal (s^i) for (i in [1, N]) is a real-valued continuous vector with fixed length T, i.e., (s^i in mathbb {R}^{1 times T}). The model outputs the predicted waveforms ({hat{s}^i}_{i=1}^N) where (forall i, hat{s}^i = f^i(x)), and a loss function is evaluated by comparing the predictions to the ground truth sources ({s^i}_{i=1}^N) up to a permutation. Explicitly, we consider a permutation-invariant objective function5,$$begin{aligned} mathscr {L}(hat{s}, s) = min _{sigma in S_N} frac{1}{N} sum _{i=1}^N ell (hat{s}^{sigma (i)}, s^i) qquad text {where} forall i, hat{s}^i = f^i(x) end{aligned}$$Here, (ell (cdot , cdot )) represents the loss function computed on an (output, target) pair, (sigma) indicates a permutation, and (S_N) is the space of permutations. In certain scenarios, we include the L2 regularization term,$$begin{aligned} mathscr {L} mapsto mathscr {L} + lambda sum _{j=1}^P beta _j^2 end{aligned}$$where (beta _j) represent the model parameters, P denotes the model complexity, and (lambda) is a hyperparameter empirically selected to minimize overfitting (i.e. enhance convergence of training and evaluation losses and metrics).For the single-channel loss function (ell), we consider a linear combination of several loss terms that compute the error in estimated waveform reconstructions ({hat{s}^i}_{i=1}^N) relative to the ground truth waveforms ({s^i}_{i=1}^N).

L1 Loss $$begin{aligned} |hat{s}^{sigma (i)} – s^{i}| end{aligned}$$ This represents the absolute error on raw time domain waveforms.

STFT L1 Loss $$begin{aligned} |text {STFT}(hat{s}^{sigma (i)}) – text {STFT}(s^{i})| end{aligned}$$ This term functions to minimize absolute error on time-frequency space representations. Empirically, the inclusion of this contribution enhances the reconstruction of signal harmonicity.

Spectral Convergence Loss $$begin{aligned} ||text {STFT}(hat{s}^{sigma (i)}) – text {STFT}(s^{i})||_F / ||text {STFT}(s^{i})||_F end{aligned}$$ where (||cdot ||_F) denotes the Frobenius norm over time and frequency. This term emphasizes high-magnitude spectral components56.

We also experimented with additional terms including L1 loss on log-magnitude spectrograms to address spectral valleys and negative SI-SDR (nSI-SDR), but the inclusion of these contributions did not yield empirical improvements in results.For macaques, we modify the training algorithm according to the representation transform and inverse transform built into the model. For the model with the fully learnable Conv1D encoder and decoder, we train using the AdamW58 optimizer with a learning rate 3e-4 and batch size 16 for 100 epochs. In order to stabilize training and avoid local minima when using handcrafted STFT and iSTFT filterbanks, we initially begin training the models for 3 epochs with batch size 16 using stochastic gradient descent (SGD) with Nesterov momentum 0.6 and learning rate 1e-3 before switching to the AdamW optimizer until reaching 100 epochs.For dolphins, we provide the model with the original mixture as input, but we use high pass-filtered source waveforms as the target, which means the separation model must additionally learn to denoise the input. We again initialize training with 3 epochs and batch size 8 using SGD with Nesterov momentum 0.6 and learning rate 1e-3 before switching to the AdamW optimizer with learning rate 3e-4 for the remaining 97 epochs. We use a similar training scheme for bats, initially training with SGD for 3 epochs before employing the optimizer switcher callback to switch to AdamW and to complete 100 epochs.Model evaluation metricsWe consider the reconstruction performance by computing evaluation metrics using an expression given by5,$$begin{aligned} mathscr {M}(hat{s}, s) = max _{sigma in S_N} frac{1}{N} sum _{i=1}^N m(hat{s}^{sigma (i)}, s^i) qquad text {where} forall i, hat{s}^i = f^i(x) end{aligned}$$where (m(cdot , cdot )) is the single-channel evaluation metric computed on permutations of (output, target) pairs.Specifically, we implement two evaluation metrics to assess reconstruction quality, including (1) SI-SDR and (2) downstream classification accuracy. We consider the signal-to-distortion ratio (SDR)2, defined as the negative log squared error normalized by reference signal energy5. However, as is commonly implemented in the human speech separation literature, we instead compute the scale-invariant SDR (SI-SDR), which disregards prediction scale by searching over gains5,40. Explicitly, SI-SDR((hat{s}, s) = -10log _{10}(|hat{s} – s|^2) + 10log _{10}(|alpha s|^2)) for optimal scaling factor (alpha = hat{s}^Ts / |s|^2).Additionally, to provide a physically interpretable metric, we evaluate the performance of the trained classifier models in labeling separated waveforms according to the predicted identity of the vocalizer. This metric assumes that the classification accuracy on a downstream task reflects the fidelity of the estimated signal relative to the ground truth source and thus serves as a proxy for reconstruction quality. More

Phylogenetic analysisFor both the mycobiont and photobiont molecular phylogenies from multi-locus sequence data (nrITS, mtSSU and RPB1 for the mycobiont (140 samples) and nrITS, psbJ-L and COX2 for the photobiont (139 samples) were inferred (Supplementary Figs. S1 and S3 online). Additionally, phylogenies based solely on the marker nrITS were calculated (Supplementary Figs. S2 and S4 online), to include samples where the additional markers were not available. Both analyses include only accessions from the study sites (Fig. 1; Table 1). The phylogenies based on the multi-locus data were congruent to the clades of the phylogenies based on the marker nrITS. Thus, in the following, the focus will be only on the latter.MycobiontThe final data matrix for the phylogeny based on the marker nrITS comprised 306 single sequences with a length of 550 bp. It included sequences of the families Lecanoraceae and Lecideaceae. The phylogenetic tree was midpoint rooted and shows a total of 19 strongly supported clades on species level, assigned to five genera. The backbone is not supported and therefore the topology will not be discussed. All genera are clearly assigned to their family level and are strongly supported. Only Lecanora physciella forms an extra clade as sister to the families Lecideaceae and Lecanoraeae, which is not the case at the multimarker phylogeny. L. physciella has still an uncertain status, because of morphological similarities to both sister families6. The clade of the genus Lecidea revealed seven species (L. andersonii, L. polypycnidophora, L. UCR1, L. sp. 5, L. lapicida, L. cancriformis and L. sp. 6), Lecanora five species (L. physciella, L. sp. 2, L. fuscobrunnea, L. cf. mons-nivis, L. sp. 3), Carbonea three species (C. sp. URm1, C. vorticosa, C. sp. 2), and Lecidella three species (L. greenii, L. siplei, L. sp. nov2). The samples allocated to the genus Rhizoplaca were monospecific (R. macleanii). The taxonomical assignment of the obtained sequences were based on the studies of Ruprecht et al.48 and Wagner et al.10.PhotobiontThe final data matrix for the phylogeny based on the marker nrITS comprised 281 single sequences with a length of 584 bp. The phylogenetic tree was midpoint rooted and shows six strongly supported clades, assigned to seven different OTU levels67, using the concept of Muggia et al.51 and Ruprecht et al.48. All of the OTUs belong to the genus Trebouxia (clades A, I, S), comprising Tr_A02, Tr_A04a, Tr_I01, Tr_I17, Tr_S02, Tr_S15 and Tr_S18. Photobiont sequences taken from Perez-Ortega et al.50, which were labelled only with numbers, were renamed to assign them to the appropriate OTUs48.Analysis of spatial distributionIn general, the most common mycobionts species were Lecidea cancriformis (94 of the 306 samples), Rhizoplaca macleanii (51 samples) and Lecidella greenii (37 samples), followed by Carbonea sp. 2 (13 samples), C. vorticosa (11 samples), Lecidea polypycnidophora (10 samples) and Lecidella siplei (10 samples; see Supplementary Fig. S5 online). Nine mycobiont species were found exclusively in area 5 (MDV, 78°S): Carbonea vorticosa, Lecanora cf. mons-nivis, L. sp. 2, Lecidea lapicida, L. polypycnidophora, L. sp. 5, L. sp. 6, L. UCR1 and Rhizoplaca macleanii. On the other hand, only Lecidea cancriformis was found in all the six areas; Lecanora fuscobrunnea was present in all the areas with the exception of area 2.The most common photobiont OTUs were Tr_A02 (165 of the 281 samples) and Tr_S02 (59 samples), both of them occurring in all the six different areas, followed by Tr_S18 (32 samples), Tr_S15 (10 samples, confined to area 5) and Tr_I01 (10 samples). However, of the 149 photobiont accessions of area 5, 134 (89.93%) were assigned to Tr_A02. This percentage is much higher than in the other areas (area 1: 44.44%, area 2: 69.23%, area 3: 21.74%, area 4a: 7.69%, area 4b: 6.67%), even if those samples with mycobionts occurring exclusively in area 5 (see above) were excluded (76.56% of the 64 remaining samples are assigned to Tr_A02).The alpha, beta and gamma diversity values are given in Table 2. For the mycobionts, the alpha diversity of the communities was the highest in area 5 (8.93, which results in nine species) and the lowest in area 4b (two species, 1.88). In contrast, for the photobionts, the lowest alpha diversity value was found in area 5 (two OTUs, 1.50) and the highest in area 4a (four OTUs, 4.06). Thus, referring to this, area 5 plays a remarkable role: compared to the other areas, it shows the highest diversity of mycobiont species on the one hand and the lowest diversity of photobiont OTUs on the other hand.Table 2 Number of lichen samples, number of identified mycobiont species and photobiont OTUs, as well as alpha, beta and gamma diversity values of mycobiont species/photobiont OTUs for the different areas.Full size tableThe beta diversity values (diversity of local assemblages) for mycobiont species and photobiont OTUs are quite similar (1.69 and 1.64, respectively). This is in contrast to gamma diversity values: the overall diversity for the different areas within the whole region is much higher for the mycobionts (ten species, 9.92) than for the photobionts (three OTUs, 3.35).For mycobionts, the overall sample coverage equals to 0.993. That means that the probability for an individual of the community to belong to a sampled species is 99.3%, or, from another point of view, the probability for an individual of the whole community to belong to a species that has not been sampled is 0.7%. The sample coverage is highest for area 4b (1.000) and lowest for area 2 (0.771). Sample coverage values of the other areas are in between (area 1: 0.895, area 3: 0.931, area 4a: 0.939, area 5: 0.981). The rarefaction/extrapolation curves for the mycobiont species (see Supplementary Fig. S6a) suggest that for any sample size up to the specified level of sample coverage of 0.95, alpha diversity within area 4b is significantly lower than alpha diversity within any other area, and alpha diversity within area 5 is significantly greater than that of area 4a and 4b (based on 95% confidence intervals).For photobionts, the overall sample coverage as well as the sample coverages of area 1, area 2, area 3, area 4b as well as area 5 is equal 1.000. Only the sample coverage of area 4a (0.951) differs. The rarefaction/extrapolation curves for the photobiont OTUs (see Supplementary Fig. S6b) suggest that for any sample size up to the specified level of sample coverage of 0.95, alpha diversity within area 1 is significantly lower than alpha diversity of area 3 and 4a and significantly greater than that of area 5. Alpha diversity of area 5 is significantly lower than that of area 1, area 3 and area 4a.Influence of environmental factors (elevation, precipitation and temperature)First, the proportion of the OTU Tr_A02 samples was significantly correlated to BIO10 means of the areas (R = 0.87, p = 0.022; see Supplementary Fig. S7 online): the higher the temperature mean values of the warmest quarter of an area, the higher the proportion of samples containing photobionts that are assigned to Tr_A02.The alpha diversity values of mycobiont species significantly positively correlated with BIO10 (R = 0.88, p = 0.021; see Supplementary Fig. S8 online): the higher the temperature mean values of the warmest quarter, the higher the mycobiont diversity within this particular area.Furthermore, the differences in mycobiont species community composition were significantly related to BIO10 (constrained principal coordinate analysis: F = 14.7137, p = 0.001, see Supplementary Fig. S9 online), BIO12 (F = 2.7535, p = 0.012), elevation (F = 2.5108, p = 0.025) and the geographic separation of the samples (Mantel statistic r = 0.1288, p = 0.0002).The differences in community composition of photobiont OTUs were related significantly to BIO10 (constrained principal coordinate analysis: F = 48.5952, p = 0.001, see Supplementary Fig. S10 online), BIO12 (F = 4.4848, p = 0.008), elevation (F = 6.8608, p = 0.002), and physical distance (Mantel statistic r = 0.4472, p = 0.0001).Haplotype analysisHaplotype networks were computed for the mycobiont species and photobiont OTUs with h ≥ 2 and at least one haplotype with n ≥ 3 (Carbonea sp. 2, Lecanora fuscobrunnea, Lecidea cancriformis, Lecidella greenii, L. siplei, L. sp. nov2 and Rhizoplaca macleanii, as well as Tr_A02, Tr_I01 and Tr_S02), in both cases based on nrITS sequence data (Figs. 2, 3). The samples of Carbonea vorticosa (11) were all assigned to a single haplotype, which was also true for Lecidea polypycnidophora (10 samples), Tr_S15 (10 samples) and Tr_S18 (32 samples). Figure 3b, c illustrate the subdivision of Tr_I0151 into Tr_I01j35,48 and Tr_I01k (in this study), and the subdivision of Tr_S02 into Tr_S0235, and Tr_S02b and Tr_S02c48.Figure 2Haplotype networks of mycobiont species with h ≥ 2 and at least one haplotype with n ≥ 3, showing the spatial distribution within the different areas, based on nrITS data. (a) Carbonea sp. 2, (b) Lecanora fuscobrunnea, (c) Lecidea cancriformis, (d) Lecidella greenii, (e) Lecidella siplei, (f) Lecidella sp. nov2, (g) Rhizoplaca macleanii. Roman numerals at the center of the pie charts refer to the haplotype IDs; the italic numbers next to the pie charts give the total number of samples per haplotype. The circle sizes reflect relative frequency within the species; the frequencies were clustered in ten (e.g. the circles of all haplotypes making up between 20 and 30% have the same size). Note: only complete sequences were included.Full size imageFigure 3Haplotype networks of photobiont OTUs with h ≥ 2 and at least one haplotype with n ≥ 3, showing the spatial distribution within the different areas, based on nrITS data. (a) Tr_A02, (b) Tr_I01, (c) Tr_S02. Roman numerals at the center of the pie charts refer to the haplotype IDs; the italic numbers next to the pie charts give the total number of samples per haplotype. The circle sizes reflect relative frequency within the species; the frequencies were clustered in ten (e.g. the circles of all haplotypes making up between 20 and 30% have the same size). Note: only complete sequences were included.Full size imageThe haplotype networks include pie charts showing the occurrence of the different haplotypes within the different areas. All haplotypes of Rhizoplaca macleanii are restricted to area 5, as well as Lecidella greenii mainly to area 5 and areas 1 and 4a, and Lecidella sp. 2 to areas 2 and 3. However, all other species do not suggest a spatial pattern with different haplotypes being specific for different areas. Moreover, the distribution turned out to be rather unspecific, with a great part of the haplotypes found in multiple areas. For the sake of completeness, additionally, haplotype networks based on multi-locus sequence data were computed for the most abundant mycobiont species and photobiont OTU with multi-locus data available (Lecidea cancriformis and Tr_S02). Not surprisingly, those networks show a greater number of different haplotypes, but they also do not allow conclusions concerning spatial patterns of area specific haplotypes (see Supplementary Fig. S11 online).Diversity and specificity indices of mycobiont species and photobiont OTUsThe diversity and specificity indices for the different mycobiont species and photobiont OTUs are given in Supplementary Table S8 online.For the sample locations of mycobiont species with n ≥ 10, BIO10 was strongly correlated to the specificity indices NRI (net relatedness index) and significantly correlated to PSR (phylogenetic species richness) and 1 – J′ (Pielou evenness index). BIO12 was significantly correlated to NRI, PSR and 1 – J′. Figure 4 illustrates these correlations: the higher the BIO10 and BIO12 mean values, the higher was the NRI (phylogenetic clustering of the photobiont symbiotic partners), the lower was the PSR (increased phylogenetically relatedness of photobiont symbiotic partners) and the higher was 1 – J′ (less numerically evenness of the photobiont symbiotic partners). Thus, for the mean values of the sample locations of a mycobiont species, a comparatively high temperature of the warmest quarter and high annual precipitation occurs with associated photobionts that are phylogenetically clustered and closer related to each other. The lowest values of NRI and the highest values of PSR were developed by Lecidea cancriformis and Lecanora fuscobrunnea, which also showed the lowest BIO10 and BIO12 mean values at their sample sites. On the contrary, the highest values of NRI and PSR were developed by Rhizoplaca macleanii, which also had the highest BIO10 and BIO12 means.Figure 4Correlation plots. Specificity indices NRI (net relatedness index), PSR (phylogenetic species richness and 1 – J′ (Pielou evenness index) against mean values of BIO10 (mean temperature of warmest quarter) and BIO12 (annual precipitation) for mycobiont species with n ≥ 10.Full size imageFor the sample locations of photobiont OTUs with n ≥ 10, elevation significantly negatively correlated with h (number of haplotypes) and Hd (haplotype diversity): the higher the mean elevation of sample sites, the lower the number of haplotypes and the lower the probability that two randomly chosen haplotypes are different (Fig. 5). The highest values of h and Hd were shown by Tr_A02, Tr_I01 and Tr_S02, which occurred at sample sites with comparatively low elevations. In contrast, Tr_S15 and Tr_S18 occurred at very high elevations and showed very low values of h and Hd.Figure 5Correlation plots. Diversity indices h (number of haplotypes) and Hd (haplotype diversity) against mean elevation of sample sites for photobiont OTUs with n ≥ 10.Full size imageAnalysis of mycobiont–photobiont associationsBipartite networks were calculated for all associations between mycobiont species (lower level) and the respective photobiont OTUs (higher level) for all areas (Fig. 6). The H2′ value (overall level of complementary specialization of all interacting species) was highest in area 2 (0.921), indicating a network with mostly specialized interactions: within this network, with the exception of Lecidea andersonii, the mycobiont species are associated exclusively with one single photobiont OTU. The second highest H2′ value was developed by area 4b (0.710); in contrast, area 4a showed the lowest H2′ value (0.260), with the most abundant mycobiont species Lecidea cancriformis showing associations with five different photobiont OTUs. The H2′ values of area 1, area 3 and area 5 indicate medium specification.Figure 6Bipartite networks showing the associations between mycobiont species and photobiont OTUs for the different areas. Rectangles represent species/OTUs, and the width is proportional to the number of samples. Associated species/OTUs are linked by lines whose width is proportional to the number of associations.Full size imageIn addition, the bipartite networks illustrate the different occurrence of mycobiont species and photobiont OTUs within the different areas: For example, in area 1 (and area 2), five (seven) different mycobiont species are associated with only three different photobiont OTUs. In contrast, in area 4b, only two different mycobiont species are associated with four different photobiont OTUs. In area 5, the number of associated photobiont OTUs is also four, but those four OTUs are associated with 16 different mycobiont species.The network matrix giving all the associations between the mycobiont species and photobiont OTUs is presented in Supplementary Table S9 online. More