Ecology

Wang, L.C.H. Ecological, physiological, and biochemical aspects of torpor in mammals and birds in Advances in comparative and environmental physiology (ed. Wang, L.C.H.) 361–401 (Springer, 1989).Humphries, M.M., Speakman, J.R., & Thomas, D.W. Temperature, hibernation energetics, and the cave and continental distributions of little brown myotis in Functional and Evolutionary Ecology of Bats (ed. Zubaid, A., McCracken, G.F., & Kunz, T.H.) 23–37 (Oxford Press, 2005).Hudson, J.W. Torpidity in mammals in Comparative physiology of thermoregulation. (ed. Whittow, G.C., Hudson, J.W., & Deavers, D.R.) 97–165 (Academic Press, 1973).Geiser, F. & Ruf, T. Hibernation versus daily torpor in mammals and birds: Physiological variables and classification of torpor patterns. Phys. Zool. 68, 935–966 (1995).Article 

Google Scholar 
Davis, W. H. & Hitchcock, H. B. Biology and migration of the bat, Myotis lucifugus, New England. J. Mamm. 46, 296–313. https://doi.org/10.2307/1377850 (1965).Article 

Google Scholar 
Speakman, J. R. & Rowland, A. Preparing for inactivity: How insectivorous bats deposit a fat store for hibernation. Proc. Nutr. Soc. 58, 123–131. https://doi.org/10.1079/PNS19990017 (1999).CAS 
Article 
PubMed 

Google Scholar 
Boyles, J. G., Johnson, J. S., Blomberg, A. & Lilley, T. M. Optimal hibernation theory. Mammal Rev. 50, 91–100. https://doi.org/10.1111/mam.12181 (2020).Article 

Google Scholar 
Britzke, E. R., Sewell, P., Hohmann, M. G., Smith, R. & Scott, R. Use of temperature-sensitive transmitters to monitor the temperature profiles of hibernating bats affected with white-nose syndrome. Northeast. Nat. 17, 239–246. https://doi.org/10.1656/045.017.0207 (2010).Article 

Google Scholar 
Halsall, A. L., Boyles, J. G. & Whitaker, J. O. Jr. Body temperature patterns of big browns during winter in a building hibernaculum. J. Mamm. 93, 497–503. https://doi.org/10.1644/11-MAMM-A-262.1 (2012).Article 

Google Scholar 
Johnson, J. S., Lacki, M. J., Thomas, S. C. & Grider, J. F. Frequent arousals from winter torpor in Rafinesque’s big-eared bat (Corynorhinus rafinesquii). PLoS ONE 7, e49754. https://doi.org/10.1371/journal.pone.0049754 (2012).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Jonasson, K. A. & Willis, C. K. R. Hibernation energetics of free-ranging little brown bats. J. Exp. Bio. 215, 2141–2149. https://doi.org/10.1242/jeb.066514 (2012).Article 

Google Scholar 
Day, K. M. & Tomasi, T. E. Winter energetics of female Indiana bats Myotis sodalis. Physiol. Biochem. Zool. 87, 56–64. https://doi.org/10.1086/671563 (2014).Article 
PubMed 

Google Scholar 
Meierhofer, M. B. et al. Winter habitats of bats in Texas. PLoS ONE 14, e0220839. https://doi.org/10.1371/journal.pone.0220839 (2019).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Boyles, J. G. Benefits of knowing the costs of disturbance to hibernating bats. Wildl. Soc. Bull. 41, 388–392. https://doi.org/10.1002/wsb.755 (2017).Article 

Google Scholar 
Frick, W. F. et al. Pathogen dynamics during invasion and establishment of white-nose syndrome explain mechanisms of host persistence. Ecology 98, 624–631. https://doi.org/10.1002/ecy.1706 (2017).Article 
PubMed 

Google Scholar 
Cheng, T. L. et al. The scope and severity of white-nose syndrome on hibernating bats in North America. Conserv. Biol. 35, 1586–1597. https://doi.org/10.1111/cobi.13739 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Cryan, P. M., Meteyer, C. U., Boyles, J. G. & Blehert, D. S. Wing pathology of white-nose syndrome in bats suggests life-threatening disruption of physiology. BMC Biol. 8, 1–8 (2010).Article 

Google Scholar 
Cryan, P. M. et al. Electrolyte depletion in white-nose syndrome bats. J. Wild. Dis. 49, 398–402 (2013).CAS 
Article 

Google Scholar 
Frick, W. F. et al. Disease alters macroecological patterns of North American bats. Glob. Ecol. Biogeogr. 24, 741–749. https://doi.org/10.1111/geb.12290 (2015).Article 

Google Scholar 
Bernard, R.F., Willcox, E.V., Parise, K.L., Foster, J.T., & McCracken, G.F. White-nose syndrome fungus, Pseudogymnoascus destructans, on bats captured emerging from caves during winter in the southeastern United States. BMC Zool. https://doi.org/10.1186/s40850-017-0021-2 (2017).Davy, C. M. et al. The other white-nose syndrome transcriptome: Tolerant and susceptible hosts respond differently to the pathogen Pseudogymnoascus destructans. Ecol. Evol. 7, 7161–7170. https://doi.org/10.1002/ece3.3234 (2015).Article 

Google Scholar 
Lilley, T. M. et al. Resistance is futile: RNA-sequencing reveals differing responses to bat fungal pathogen in Nearctic Myotis lucifugus and Palearctic Myotis myotis. Oecologia 191, 295–309. https://doi.org/10.1007/s00442-019-04499-6 (2019).ADS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Bernard, R. F. & McCracken, G. F. Winter behavior of bats and the progression of white-nose syndrome in the southeastern United States. Ecol. Evol. 7, 1487–1496. https://doi.org/10.1002/ece3.2772 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Moosman, P.R., Warner, D.P., Hendren, R.H., & Hosler, M.J. Potential for monitoring eastern small-footed bats on talus slopes. Northeast. Nat. https://www.jstor.org/stable/26453719 (2015).Bernard, R. F. et al. Identifying research needs to inform white-nose syndrome management decisions. Cons. Sci. Prac. 2(e220), 2020. https://doi.org/10.1111/csp2.220 (2020).Article 

Google Scholar 
Reynolds, D. S., Shoemaker, K., Oettingen, S. V. & Najjar, S. High rates of winter activity and arousals in two New England bat species: implications for a reduced white-nose syndrome impact?. Northeast. Nat. 24, B188–B208. https://doi.org/10.1656/045.024.s720 (2017).Article 

Google Scholar 
Bernard, R. F., Willcox, E. V., Jackson, R. T., Brown, V. A. & McCracken, G. F. Feasting, not fasting: Winter diets of cave hibernating bats in the United States. Front. Zool. 18, 1–13. https://doi.org/10.1186/s12983-021-00434-9 (2021).Article 

Google Scholar 
Prendergast, B. J., Freeman, D. A., Zucker, I. & Nelson, R. J. Periodic arousal from hibernation is necessary for initiation of immune responses in ground squirrels. Am. J. Phys. 282, 1054–1062. https://doi.org/10.1152/ajpregu.00562.2001 (2002).Article 

Google Scholar 
Dobony, C. A. et al. Little brown myotis persist despite exposure to white-nose syndrome. J. Fish Wild. 2, 190–195. https://doi.org/10.3996/022011-JFWM-014 (2011).Article 

Google Scholar 
Rowley, J. J. & Alford, R. A. Hot bodies protect amphibians against chytrid infection in nature. Sci. Rep. 3, 1–4. https://doi.org/10.1038/srep01515 (2013).CAS 
Article 

Google Scholar 
Verant, M. L. et al. White-nose syndrome initiates a cascade of physiologic disturbances in the hibernating bat host. BMC Physiol. 14, 10 (2014).Article 

Google Scholar 
Verant, M. L. et al. Temperature-dependent growth of Geomyces destructans, the fungus that causes bat white-nose syndrome. PLoS ONE 7, e46280. https://doi.org/10.1371/journal.pone.0046280 (2012).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Brownlee-Bouboulis, S. A. & Reeder, D. M. White-nose syndrome-affected little brown Myotis (Myotis lucifugus) increase grooming and other active behaviors during arousals from hibernation. J. Wildl. Dis. 49, 850–859. https://doi.org/10.7589/2012-10-242 (2013).Article 
PubMed 

Google Scholar 
Campbell, J. Tennessee winter bat population and white-nose syndrome monitoring report for 2018–2019. TWRA Wildlife Technical Report 16-4. http://www.tnbwg.org/2019%20Annual%20Monitoring%20Report.pdf (2019).Langwig, K. E. et al. Sociality, density-dependence and microclimates determine the persistence of populations suffering from a novel fungal disease, white-nose syndrome. Eco. Let. 15, 1050–1057. https://doi.org/10.1111/j.1461-0248.2012.01829.x (2012).Article 

Google Scholar 
Langwig, K.E., et al. Drivers of variation in species impacts for a multi-host fungal disease of bats. Philos. Trans. R. Soc. B Biol. Sci. 371, 20150456. https://doi.org/10.1098/rstb.2015.0456 (2016).Aldridge, H.D.J.N., & Brigham, R.M. Load carrying and maneuverability in an insectivorous bat: A test of the 5% ‘rule’ of radio-telemetry. J. Mamm. 69, 379–382. https://doi.org/10.2307/1381393 (1988)Sikes, R. S. et al. Guidelines of the American Society of Mammalogists for the use of wild mammals in research and education. J. Mamm. 97, 663–688. https://doi.org/10.1093/jmammal/gyw078 (2016).Article 

Google Scholar 
Barclay, R. M. R. et al. Can external radiotransmitters be used to assess body temperature and torpor in bats?. J. Mammal. 77, 1102–1106. https://doi.org/10.2307/1382791 (1996).Article 

Google Scholar 
Turbill, C., & Geiser, F. Hibernation by tree-roosting bats. J. Comp. Physiol. B. 178, 597–605. https://doi.org/10.1007/s00360-007-0249-1 (2008)Park, K.J., Jones, G., & Ransome, R.D. Torpor, arousal and activity of hibernating greater horseshoe bats (Rhinolophus ferrumequinum). Funct. Ecol. 14, 580–588 (2000).Reeder, D. M. et al. Frequent arousal from hibernation linked to severity of infection and mortality in bats with white-nose syndrome. PLoS ONE 7, e38920. https://doi.org/10.1371/journal.pone.0038920 (2012).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Sirajuddin, P. Vulnerability of tri-colored bats (Perimyotis subflavus) to white-nose syndrome in the southeastern United States. M.S thesis at https://www.proquest.com/dissertations-theses/vulnerability-tri-colored-bats-i-perimyotis/docview/2185672601/se-2?accountid=8361 (2018).R Development Core Team 3.6.1. A Language and Environment for Statistical Computing. http://www.r-project.org (2019).Bates, D., Maechler, M., Bolker, B., & Walker, S. lme4: Linear mixed-effects models using ‘Eigen’ and S4. R Packag. version 1.1-13. ftp://cran.r-project.org/pub/R/web/packages/lme4/lme4.pdf. (2017)Czenze, Z. J. & Willis, C. K. R. Warming up and shipping out: Arousal and emergence timing in hibernating little brown bats (Myotis lucifugus). J. Comp. Physiol. B-Biochem. Syst. Environ. Physiol. 185, 575–586. https://doi.org/10.1007/s00360-015-0900-1 (2015).Article 

Google Scholar 
Conover, W.J. Practical Nonparametric Statistics. (Wiley, 1999).Crawley, M.J. The R Book: Second Edition. (Wiley, 2013)Lenth, R., Singmann, H., Love, J., Buerkner, P., & Herve, M. Package ‘emmeans’. R Packag. version 1.7.0. https://cran.r-project.org/web/packages/emmeans/emmeans.pdf (2021)Best, T.L. & Jennings, J.B. Myotis leibii. Mammalian Species 547, 1–6. https://doi.org/10.2307/3504255 (1997)Frank, C.L., Herzog, C.,. Laske, J.P., & Cardino V. The role of skin temperature in the resistance of Myotis leibii to white-nose syndrome presented at the 49th Symposium of the North American Society for Bat Research, Kalamazoo, Michigan https://www.nasbr.org/resources/docs/meetings/year/2019/NASBR-2019-Abstracts-Draft-191001.pdf (2019)Johnson, J.S., Scafini, M.R., Sewall, B.J., & Turner, G.G. Hibernating bat species in Pennsylvania use colder winter habitats following the arrival of white-nose syndrome in Conservation and Ecology of Pennsylvania’s Bats (ed. Butchkoski, C.M., Reeder, D.M., Turner, G.G., & Whidden, H.P) 181–199 (The Pennsylvania Academy of Science, 2016)Haase, C.G., et al. Body mass and hibernation microclimate may predict bat susceptibility to white‐nose syndrome. Eco. Evo. 11, 506–515 https://doi.org/10.1002/ece3.7070 (2021)Veilleux, J.P. A Noteworthy Hibernation Record of Myotis leibii (Eastern Small-footed Bat) in Massachusetts. Northeast. Nat. 14, 501–502 https://doi.org/10.1656/1092-6194(2007)14[501:ANHROM]2.0.CO;2 (2007)Boyles, J. G., Dunbar, M. B. & Whitaker, J. O. Activity following arousal in winter in North American vespertilionid bats. Mamm. Rev. 36, 267–280. https://doi.org/10.1111/j.1365-2907.2006.00095.x (2006).Article 

Google Scholar 
Webb, P. I., Speakman, J. R. & Racey, P. A. How hot is a hibernaculum? A review of the temperatures at which bats hibernate. Can. J. Zool. 74, 761–765. https://doi.org/10.1139/z96-087 (1996).Article 

Google Scholar 
Jackson, R.T., Willcox, E.V., Zobel, J.M., & Bernard, R.F. Hibernation behavior of four bat species with differing susceptibility to white-nose syndrome. In review (2021).Fujita, M. S. & Kunz, T. H. Pipstrellus subflavus. Mamm. Spec. 228, 1–6. https://doi.org/10.2307/3504021 (1984).Article 

Google Scholar 
Bohn, S. J. et al. Evidence of ‘sickness behaviour’ in bats with white-nose syndrome. Behaviour 152, 981–1003. https://doi.org/10.1163/1568539X-00003384 (2016).Article 

Google Scholar 
Tuttle, M. D. Status, causes of decline, and management of endangered gray bats. J. Wild. Manag. 43, 1–17. https://doi.org/10.2307/3800631 (1979).Article 

Google Scholar 
Harvey, M.J., Altenbach, J.S., & Best, T.L. Bats of the Eastern United States. (JHU Press, 2011).Klüg-Baerwald, B. J., Lausen, C. L., Willis, C. K. & Brigham, R. M. Home is where you hang your bat: Winter roost selection by prairie-living big brown bats. J. Mamm. 98, 752–760. https://doi.org/10.1093/jmammal/gyx039 (2017).Article 

Google Scholar 
Sandel, J. K. et al. Use and selection of winter hibernacula by the eastern pipistrelle (Pipistrellus subflavus) in Texas. J. Mamm. 82, 173–178. https://doi.org/10.1644/1545-1542(2001)082%3c0173:UASOWH%3e2.0.CO;2 (2001).Article 

Google Scholar 
Hayes, J. P., Ober, H. K., Sherwin, R. E. Survey and monitoring of bats in Ecological and behavioral methods for the study of bats (ed. Kunz, T. H., Parsons, S.). 120–129 (JHU Press, 2009). More

Study areaThe study area covers the total land area of Thailand. Where it is useful, we divided Thailand into six regions (Fig. 2a), the names and boundaries of which are widely used, although they have no official administrative status. We focused on the elements of Thailand’s protected area system that were concerned principally with the in-situ conservation of biodiversity: existing and proposed National Parks, Wildlife Sanctuaries, Non-hunting Areas, and Forest Parks, covering 111, 201 km2 or 21.7% of the country’s land area37 (Fig. 1).Environmental dataA set of environmental variables that were expected to be directly or indirectly related to species distributions in Thailand was used to model suitable habitat in the present and future (Supplementary Material Table S1). These variables were chosen to encompass ecologically relevant variables and enable consistent comparison between species, regardless of species-specific preferences. GIS layers for the whole of the study area were compiled using a variety of data sources at 1-km2 resolution. For variables originally at higher than 1-km resolutions, we used the plus function in ArcMap to combine them with a mask of the study area to use the mask dimensions for all cells.The physical variables, altitude, slope, aspect, and soil pH are widely used in species distribution modeling. Slope and aspect have biologically significant impacts on both temperature and rainfall at these latitudes8 and are particularly important at the poleward margins of species ranges where species may be confined to one aspect. Slope also affects soil maturity and depth. Soil pH is a consistently measured soil variable that broadly correlates with fertility in tropical soils8. Additional soil variables, particularly soil phosphorus, have been shown to be important filters of plant species distributions in the tropics38, but they are not available for Thailand with a useful accuracy and spatial resolution. Altitudes were downloaded from the CGIAR-Consortium for Spatial Information, CGIAR-CSI version 4.1. Slope and aspect were generated by using surface tools in ArcGIS. Soil pH was extracted from ISRIC-World Soil Information version 2.0.Unlike the temperate zone, where tolerances of winter cold and requirements for summer warmth dominate plant and animal distributions, our understanding of how tropical climates filter species distributions is still weak38,39. In Thailand, as in most of the tropics, there are two major climatic gradients which correlate with changes in species composition: a rainfall gradient in the lowlands, along which total rainfall declines and the length of the dry season increases, and a gradient of steadily declining temperature with elevation7. There is no simple relationship between elevation, and thus temperature, and rainfall. An additional complication is that temperature seasonality may be significant in northern Thailand (north of c.18° N), where cooler winters reduce dry-season water stress and extreme low temperatures at high altitudes may exceed physiological tolerances. We therefore chose 8 bioclimatic variables (Supplementary Material Table S1) related to precipitation and temperature, and their seasonality, all of which have previously been used in species distribution modelling in this region9,40. These are available at a resolution of 30 arc sec (approximately 1 km at the equator) from WorldClim ver. 1.4 based on averages of 1970–1990. These variables are available from the same source (and downscaled using the same methods) for the future climate projections.Vegetation structure is an additional major influence on plant and animal distributions in the tropics, both in intact natural vegetation38,39 and when the original vegetation has been degraded or cleared8. Vegetation structure was represented through the inclusion of two continuous variables, percentage forest cover and tree density, as most of the modelled species are known to be sensitive to both the presence of forest and the degree of intactness of the tree cover9. Mean tree density per km2 was extracted from Crowther et al.41 version 2 and percentage coverage of forest per km2 was extracted from the European Space Agency (ESA) GlobCover Version 2.3.Note that the mechanistic basis of the correlations between all these variables and the current distributions of tropical plants and animals are rarely known. Temperature has a direct physiological impact on all organisms, and water supply may be seasonally limiting for plants and some amphibians, but indirect links through biotic interactions are expected to be more important in the tropics, including pest pressure on plants38 and food supply for animals39. Competition is probably also important in shaping local species assemblies. For future projections, we assumed that temperature and precipitation were changing, and that other variables (topography, soil, and vegetation) were stable, so our analysis represents the impacts of climate alone. For 2070, we used the same variables projected by three CMIP5 Earth System Models, CNRM-CM5, GFDL-CM3 and HadGEM2-ES, which have been previously used in Southeast Asia9,42 and in Thailand7. We used two Representative Concentration Pathways, RCP2.6 and RCP8.5, representing low and high greenhouse-gas concentration scenarios, respectively, and thus the potential range of radiative forcing by the end of the century43. RCP2.6 is consistent with meeting the Paris Agreement’s 2 °C global warming target.Species occurrence dataMany locality records for vertebrates were supplied by the Department of National Parks, Wildlife and Plant Conservation (DNP). Trained DNP staff walked along trails throughout the protected areas in Thailand during 2017–2018. They recorded 271,695 locations for 70 mammal species, 18 locations for 3 amphibian species, 318 locations for 18 reptile species, and 43,057 locations for 65 bird species44. We supplemented this with data downloaded from the Global Biodiversity Information Facility (GBIF, https://www.gbif.org/) for 1960–2019 for amphibians (2063 localities for 86 species)45, reptiles (1722 localities from 196 species)46, mammals (2508 localities from 191 species)47, and birds (1,559,222 localities from 884 species)48. More than 95% of the bird records from GBIF were identified as coming from eBird49, which is popular among birders in Thailand. For plants, we used occurrence data from the DNP’s forest resource inventory project from 221 plots, including 24,605 localities for 363 species, the DNP’s Forest Herbarium, including 227 localities for 141 species, and locations for 12 rare and endangered forest species collected from all over Thailand. We also downloaded data from the Botanical Information and Ecology Network (BIEN, https://bien.nceas.ucsb.edu/bien/), including 7209 localities for 1422 species.We removed suspect records (coordinate issues, name problems, etc.), duplicates from the same locality (i.e., more than one individual of the same species recorded in a cell), and species with  0.5 as adequate, but since only five SDMs out of the 1457 generated in this study had values lower than this (0.3–0.5), we retained all the models.Assessment of climate change impactsThe estimated current distribution for each species from Maxent was used as the baseline for comparison with projected distributions of suitable habitat for these species by 2070, under the two emission scenarios and three ESMs, and with and without unlimited dispersal into newly available habitat. We then assessed the impacts of climate change, both on the spatial distribution of individual species and on the pattern of species richness. To generate a species richness map, the binary habitat suitability maps for all species were stacked to produce a consolidated map, which showed the number of species for each 1 km grid cell, and then classified them into five classes (lowest, low, moderate, high, and highest), using the mean ± standard deviation as a break class40.Current and future maps were then compared for each species to calculate the change in species richness, and contingency tables were generated containing the numbers of cells (each of 1 km2) in each richness class. Suitable habitat areas were calculated for the current climate and projected for the future climate. For each species we estimated gained habitat as the areas that will become suitable for a species in future under that scenario, lost habitat as the areas currently predicted as suitable now but projected to become unsuitable under future climatic change, and stable habitat as the areas predicted as suitable now which will remain suitable into the future.We then assessed the vulnerability of each species by estimating the projected change in its range over the next 50 years and using a criteria-based approach, which combined the mean of the suitable habitat area (interpreted as equivalent to extent of occurrence) in the three models and a simplified version of the IUCN Red List criteria51. For 2070, we modified criterion A3(c) as follows; Extinct (Ex) species are projected to lose 100% of suitable habitat by 2070, Critically Endangered (CR) species are projected to lose over 80%, Endangered (EN) species are projected to lose 50–80%; Vulnerable (VU) species are projected to lose 30–50%, Near Threatened (NT) species are projected to lose  More

Collection of wolf hair samplesHair samples were collected by researchers from opportunistically found-dead wolves upon standard necropsy (all the Alpine and part of the Iberian samples) or in the field (all the Dinaric-Balkan and most of the Iberian samples), or from legally harvested wolves (only in the Scandinavian population). At the time of sample collection, wolves were legally harvested in Sweden, Slovenia, and Spain, and under total protection in Portugal and Italy. Hair samples were collected from four body regions, when possible: lumbar (n = 133), dorsal cervical (n = 66), tail (n = 33) and ventral thorax (n = 27) (Tables S1 and S2). The hair was cut as close as possible to the skin with scissors to avoid collecting hair follicles, but in some samples, hairs were pulled from the carcass. Samples were stored at room temperature in paper envelopes. Age, sex, date, and cause of death/capture, geographical location, body mass, and total length were obtained for most of the wolves.Age was estimated by the dental eruption and wear or cementum age analysis and classified as ‘juveniles’ ( 2 years)40, or ‘unknown’. Sex was assessed by inspection of genitalia. Causes of death were classified as ‘acute’, likely lasting minutes to hours (vehicle accident and legal or illegal shooting); ‘subacute’, likely lasting hours to days (drowning, poisoning, trapping and intraspecific aggression); ‘chronic’, likely lasting several weeks (infectious diseases—canine distemper, canine parvovirosis, leptospirosis; sarcoptic mange; or neoplastic diseases) or ‘unknown’. Total length was obtained by measuring with metric tape (1 mm precision) the distance from snout to the distal end of the last tail vertebrae. The body mass was measured with 100 g precision with scales.The detailed protocol for the handling of wolves live trapped in the scope of ecological and conservation studies (n = 7, all from the Iberian population) has been previously described5. Traps were monitored twice every day, in the early morning and late afternoon, hence the duration of restraint after capture was unknown for 8 wolves, potentially up to 12 h. Trap-alarms were deployed in the capture of 2 wolves, with 41 and 70 min intervals between activation of the alarm and administration of the drugs. Live trapping was conducted under permits issued by the nature conservation authorities of Portugal (Instituto de Conservação da Natureza e das Florestas: 338/2007/CAPT, 258/2008/CAPT, 286/2008/CAPT, 260/2009/CAPT, 332/2010/MANU, 333/2010/CAPT, 336/2010/MANU, 26/2012/MANU, and 72/2014/CAPT) and Spain (Dirección Xeral de Conservación da Naturaleza, Xunta de Galicia: E-0020/13-PNPE, 095/2013; Consejería de Medio Ambiente, Principado de Asturias: 31/08/2017-BOPA 05/09/17) and according to European Union directives on the protection of animals used for scientific purposes (Directive 2010/63/EU) and international wildlife standards41,42. The study was undertaken in compliance with the ARRIVE guidelines43.Cortisol extractionThe protocol for the extraction of cortisol from the hair was adapted from previously described procedures15,27. Forty mg of guard hairs were separated from the undercoat and placed in 15 ml falcon tubes. Hair follicles were cut whenever found in the sample. For each sample, the length of three intact hairs was recorded. The samples were washed twice with 40 µl of distilled water/mg hair and three times with the same amount of isopropanol. In each washing step, the samples and washing solution were vortexed, the supernatant discarded, and the hair dried using clean paper towels. After the final wash, samples were dried overnight at room temperature and 30 mg of hair cut into a 2 ml polypropylene screw cap plastic tube with five 4 mm steel beads added to each tube.The hair was ground to a fine powder in a FastPrep sample homogenizer (MP Biomedicals, USA) for four times 1 min at 6.0 m/s. 50 µl methanol/mg hair were added to each sample and sonicated for 30 min at 50 Hz at 50 °C. The samples were incubated for 18 h at 50 °C in an orbital shaker at 160 rpm, centrifuged for 15 min at 14,000g at 20 °C, and 1000 µl of supernatant was collected to a screw cap glass chromatography vial and dried at room temperature in a gentle stream of nitrogen gas. Due to restrictions on laboratory use during the SARS-Cov-2 pandemic, some batches of samples were instead evaporated overnight on a suction hood. This unexpected change in the methanol evaporation protocol was recorded and accounted for in the statistical analysis.Cortisol quantificationA commercial competitive ELISA kit (Cortisol free in Saliva ELISA, Demeditec, Germany) was used to quantify the concentration of cortisol, following the manufacturer’s instructions. The kit plate wells are provided coated with polyclonal rabbit antibody against cortisol, and cortisol-horseradish peroxidase was used as conjugate. According to the manufacturer, the cross-reactivity of the test to selected steroids is low (Table S3), the intra-assay variation is 3.8–5.8% and the inter-assay variation is 6.2–6.4%. Samples, standards, and controls were tested in duplicate.The 4-parameter standard curve was calculated from the log-transformed cortisol concentration of the standard solutions and their measured OD45044. Standard curves were estimated using the software GraphPad Prism 6.04 (GraphPad Software, La Jolla, California USA), and yielded an average R2adjusted = 0.991 (range 0.968–0.999). The cortisol concentration of the reconstituted samples was estimated from the standard curve and converted to cortisol concentration as picograms (pg) of cortisol/mg of guard hair.Intra and inter-assay coefficients of variation were estimated for six ELISA assays of 37–40 samples each. The low and high controls included in the kit were used to estimate the inter-assay coefficient of variation and the duplicate runs of each sample were used to estimate the intra-assay coefficient of variation. Linearity was assessed by two-fold dilutions (1:1, 1:2, 1:4 and 1:8) of 4 extracted samples, comparing the expected and observed concentrations. Recovery was assessed by spiking 6 ground hair samples with known concentrations of cortisol (50, 25, 12.5, 6.25 pg/mg, and no spiking), comparing the expected and observed concentrations.The intra-assay coefficient of variation of the ELISA assays ranged from 6.50 to 9.97% (average 7.66%). The inter-assay coefficient of variation was 11.54% for the low concentration controls and 9.08% for the high concentration controls (average 10.31%). Assay linearity was 91% for the 1:2 dilution, 103% for 1:4, and 117% for 1:8 (average 103%). The recovery of cortisol averaged 94%, being 73% for the 50 pg/mg spiked samples, 74% for 25 pg/mg, 95% for 12.5 pg/mg, and 113% for 6.25 pg/mg.Determinants of hair cortisol concentrationThe potential determinants of HCC investigated included wolf intrinsic variables: sex, age, body condition, body structural size, month of death/capture, and wolf population. The scaled mass index was selected as a measure of body condition45 and estimated from the log-transformed body weight (g) and total length (mm). Log-transformed total length was used as an indicator of body structural size46. Samples were assigned to the Iberian, Alpine, Dinaric-Balkan, or Scandinavian wolf populations16 from the geographical location of the death or live-trapping sites (Fig. 1).The relationship between HCC and additional variables related to the sampling procedure or to the work conducted in the laboratory (length of hair used for cortisol extraction, sample storage time, body region, cause of death/capture, and methanol evaporation protocol), herein referred to as methodological variables, was also investigated as potential confounding variables. Sample storage time was the period in months between death/capture and cortisol extraction. In those samples for which only the year of death was available, 30 June was assigned as the date of death, solely to estimate storage time. All continuous variables were standardized to their z-scores.Statistical analysisFirst, the effect of body region was investigated by a linear mixed model with HCC as the dependent variable, and the independent variables body region, as a categorical fixed effect, and individual wolf, as a random effect. The lumbar region was set as the reference class as it was the most represented in our sample (Table S1). Data from 27 wolves for which samples were available from all 4 body regions were used in this analysis. Four outliers in the dataset violated the assumption of normality in the residuals of the model comparing HCC across body regions (Fig. S1A) and were excluded from this model’s dataset (Fig. S1B).Second, the effect of intrinsic and methodological variables on HCC from the lumbar body region was investigated by another linear mixed model with sex, age, body condition, body structural size (standardized log-transformed total length), cause of death/capture, wolf population, hair length, sample storage time, and methanol evaporation protocol as fixed effect independent variables. The month of death/capture was included as a random effect. Reference classes for the categorical variables were set as female, adult, acute death, Iberian population, and methanol evaporation by nitrogen gas stream. Two outliers in the dataset violated the assumption of normality in the residuals of the model (Full model, Table S4) and were excluded from this analysis (Fig. S1C,D).The goal of this analysis was to assess the relationship between HCC and wolf intrinsic variables, controlling for the potential confounding effect of the methodological variables. Starting from the full model (Table S4), models including all possible combinations of variables were ranked by their AICc using the package “MuMIn”47 in R 3.6.148. The most supported model was selected for inference and models with ΔAICc  More

Lozupone CA, Knight R. Global patterns in bacterial diversity. Proc Natl Acad Sci USA. 2007;104:11436–40.CAS 
PubMed 
PubMed Central 

Google Scholar 
Logares R, Bråte J, Bertilsson S, Clasen JL, Shalchian-Tabrizi K, Rengefors K. Infrequent marine-freshwater transitions in the microbial world. Trends Microbiol. 2009;17:414–22.CAS 
PubMed 

Google Scholar 
Cavalier-Smith T. Megaphylogeny, cell body plans, adaptive zones: causes and timing of eukaryote basal radiations. J Eukaryot Microbiol. 2009;56:26–33.PubMed 

Google Scholar 
Nakov T, Beaulieu JM, Alverson AJ. Diatoms diversify and turn over faster in freshwater than marine environments. Evolution. 2019;73:2497–511.PubMed 

Google Scholar 
Dittami SM, Heesch S, Olsen JL, Collén J. Transitions between marine and freshwater environments provide new clues about the origins of multicellular plants and algae. J Phycol. 2017;53:731–45.PubMed 

Google Scholar 
Dickson B, Yashayaev I, Meincke J, Turrell B, Dye S, Holfort J. Rapid freshening of the deep North Atlantic Ocean over the past four decades. Nature. 2002;416:832–7.CAS 
PubMed 

Google Scholar 
Aretxabaleta AL, Smith KW, Kalra TS. Regime changes in global sea surface salinity trend. J Mar Sci Eng. 2017;5:57.
Google Scholar 
López-Maury L, Marguerat S, Bähler J. Tuning gene expression to changing environments: from rapid responses to evolutionary adaptation. Nat Rev Genet. 2008;9:583–93.PubMed 

Google Scholar 
Björck S. A review of the history of the Baltic Sea, 13.0-8.0 ka BP. Quat Int. 1995;27:19–40.
Google Scholar 
Krauss W. Baltic sea circulation. In: Steele JH, editor. Encyclopedia of ocean sciences. Oxford: Academic Press; 2001. p. 236–44.Telesh I, Schubert H, Skarlato S. Life in the salinity gradient: discovering mechanisms behind a new biodiversity pattern. Estuar Coast Shelf Sci. 2013;135:317–27.
Google Scholar 
Johannesson K, Le Moan A, Perini S, André C. A Darwinian laboratory of multiple contact zones. Trends Ecol Evol. 2020;35:1021–36.PubMed 

Google Scholar 
Olofsson M, Hagan JG, Karlson B, Gamfeldt L. Large seasonal and spatial variation in nano- and microphytoplankton diversity along a Baltic Sea-North Sea salinity gradient. Sci Rep. 2020;10:17666.CAS 
PubMed 
PubMed Central 

Google Scholar 
Sjöqvist C, Godhe A, Jonsson PR, Sundqvist L, Kremp A. Local adaptation and oceanographic connectivity patterns explain genetic differentiation of a marine diatom across the North Sea-Baltic Sea salinity gradient. Mol Ecol. 2015;24:2871–85.PubMed 
PubMed Central 

Google Scholar 
Jochem F. Distribution and importance of autotrophic ultraplankton in a boreal inshore area (Kiel Bight, Western Baltic). Mar Ecol Prog Ser. 1989;53:153–68.
Google Scholar 
Wasmund N, Nausch G, Gerth M, Busch S, Burmeister C, Hansen R, et al. Extension of the growing season of phytoplankton in the western Baltic Sea in response to climate change. Mar Ecol Prog Ser. 2019;622:1–16.CAS 

Google Scholar 
van Wirdum F, Andrén E, Wienholz D, Kotthoff U, Moros M, Fanget A-S, et al. Middle to Late Holocene variations in salinity and primary productivity in the Central Baltic Sea: a multiproxy study from the Landsort Deep. Front Mar Sci. 2019;6:51.
Google Scholar 
Alverson AJ. Timing marine–freshwater transitions in the diatom order Thalassiosirales. Paleobiology. 2014;40:91–101.
Google Scholar 
Nakov T, Beaulieu JM, Alverson AJ. Insights into global planktonic diatom diversity: the importance of comparisons between phylogenetically equivalent units that account for time. ISME J. 2018;12:2807–10.PubMed 
PubMed Central 

Google Scholar 
Kremp A, Godhe A, Egardt J, Dupont S, Suikkanen S, Casabianca S, et al. Intraspecific variability in the response of bloom-forming marine microalgae to changed climate conditions. Ecol Evol. 2012;2:1195–207.PubMed 
PubMed Central 

Google Scholar 
Olofsson M, Kourtchenko O, Zetsche E-M, Marchant HK, Whitehouse MJ, Godhe A, et al. High single-cell diversity in carbon and nitrogen assimilations by a chain-forming diatom across a century. Environ Microbiol. 2019;21:142–51.CAS 
PubMed 

Google Scholar 
Olofsson M, Almén A-K, Jaatinen K, Scheinin M. Temporal escape – adaptation to eutrophication by Skeletonema marinoi. FEMS Microbiol Lett. 2022;fnac011. https://pubmed.ncbi.nlm.nih.gov/35137038/.Godhe A, Härnström K. Linking the planktonic and benthic habitat: genetic structure of the marine diatom Skeletonema marinoi. Mol Ecol. 2010;19:4478–90.PubMed 

Google Scholar 
Dobin A, Gingeras TR. Mapping RNA-seq reads with STAR. Curr Protoc Bioinform. 2015;51:11.14.1–11.14.19.
Google Scholar 
Anders S, Pyl PT, Huber W. HTSeq-a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31:166–9.CAS 
PubMed 
PubMed Central 

Google Scholar 
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 
Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40.CAS 
PubMed 
PubMed Central 

Google Scholar 
Aramaki T, Blanc-Mathieu R, Endo H, Ohkubo K, Kanehisa M, Goto S, et al. KofamKOALA: KEGG Ortholog assignment based on profile HMM and adaptive score threshold. Bioinformatics. 2020;36:2251–2.CAS 
PubMed 

Google Scholar 
Emms DM, Kelly S. OrthoFinder: phylogenetic orthology inference for comparative genomics. Genome Biol. 2019;20:238.PubMed 
PubMed Central 

Google Scholar 
Almagro Armenteros JJ, Salvatore M, Emanuelsson O, Winther O, von Heijne G, Elofsson A, et al. Detecting sequence signals in targeting peptides using deep learning. Life Sci Alliance. 2019;2:e201900429.PubMed 
PubMed Central 

Google Scholar 
Gruber A, Rocap G, Kroth PG, Armbrust EV, Mock T. Plastid proteome prediction for diatoms and other algae with secondary plastids of the red lineage. Plant J. 2015;81:519–28.CAS 
PubMed 
PubMed Central 

Google Scholar 
Bendtsen JD, Nielsen H, von Heijne G, Brunak S. Improved prediction of signal peptides: SignalP 3.0. J Mol Biol. 2004;340:783–95.PubMed 

Google Scholar 
Gschloessl B, Guermeur Y, Cock JM. HECTAR: a method to predict subcellular targeting in heterokonts. BMC Bioinforma. 2008;9:393.
Google Scholar 
Claros MG. MitoProt, a Macintosh application for studying mitochondrial proteins. Comput Appl Biosci. 1995;11:441–7.CAS 
PubMed 

Google Scholar 
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.CAS 

Google Scholar 
Van den Berge K, Soneson C, Robinson MD, Clement L. stageR: a general stage-wise method for controlling the gene-level false discovery rate in differential expression and differential transcript usage. Genome Biol. 2017;18:151.PubMed 
PubMed Central 

Google Scholar 
Heller R, Manduchi E, Grant GR, Ewens WJ. A flexible two-stage procedure for identifying gene sets that are differentially expressed. Bioinformatics. 2009;25:1019–25.CAS 
PubMed 

Google Scholar 
Alexa A, and Rahnenfuhrer J. topGO: Enrichment Analysis for GeneOntology. R package version 2.44.0. 2021. https://bioconductor.org/packages/release/bioc/html/topGO.html.Wu D, Smyth GK. Camera: a competitive gene set test accounting for inter-gene correlation. Nucleic Acids Res. 2012;40:e133.CAS 
PubMed 
PubMed Central 

Google Scholar 
Supek F, Bošnjak M, Škunca N, Šmuc T. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS ONE. 2011;6:e21800.CAS 
PubMed 
PubMed Central 

Google Scholar 
Bussard A, Corre E, Hubas C, Duvernois-Berthet E, Le Corguillé G, Jourdren L, et al. Physiological adjustments and transcriptome reprogramming are involved in the acclimation to salinity gradients in diatoms. Environ Microbiol. 2017;19:909–25.CAS 
PubMed 

Google Scholar 
Matthijs M, Fabris M, Obata T, Foubert I, Franco-Zorrilla JM, Solano R, et al. The transcription factor bZIP14 regulates the TCA cycle in the diatom Phaeodactylum tricornutum. EMBO J. 2017;36:1559–76.CAS 
PubMed 
PubMed Central 

Google Scholar 
Kong L, Price NM. Transcriptomes of an oceanic diatom reveal the initial and final stages of acclimation to copper deficiency. Environ Microbiol. 2021;24:951–66.Amato A, Sabatino V, Nylund GM, Bergkvist J, Basu S, Andersson MX, et al. Grazer-induced transcriptomic and metabolomic response of the chain-forming diatom Skeletonema marinoi. ISME J. 2018;12:1594–604.CAS 
PubMed 
PubMed Central 

Google Scholar 
Maumus F, Allen AE, Mhiri C, Hu H, Jabbari K, Vardi A, et al. Potential impact of stress activated retrotransposons on genome evolution in a marine diatom. BMC Genomics. 2009;10:624.PubMed 
PubMed Central 

Google Scholar 
Pargana A, Musacchia F, Sanges R, Russo MT, Ferrante MI, Bowler C, et al. Intraspecific diversity in the cold stress response of transposable elements in the diatom Leptocylindrus aporus. Genes. 2019;11:9.PubMed Central 

Google Scholar 
Smith SR, Dupont CL, McCarthy JK, Broddrick JT, Oborník M, Horák A, et al. Evolution and regulation of nitrogen flux through compartmentalized metabolic networks in a marine diatom. Nat Commun. 2019;10:4552.PubMed 
PubMed Central 

Google Scholar 
Kageyama H, Tanaka Y, Shibata A, Waditee-Sirisattha R, Takabe T. Dimethylsulfoniopropionate biosynthesis in a diatom Thalassiosira pseudonana: Identification of a gene encoding MTHB-methyltransferase. Arch Biochem Biophys. 2018;645:100–6.CAS 
PubMed 

Google Scholar 
Nakov T, Judy KJ, Downey KM, Ruck EC, Alverson AJ. Transcriptional response of osmolyte synthetic pathways and membrane transporters in a euryhaline diatom during long-term acclimation to a salinity gradient. J Phycol. 2020;56:1712–28.CAS 
PubMed 

Google Scholar 
Kageyama H, Tanaka Y, Takabe T. Biosynthetic pathways of glycinebetaine in Thalassiosira pseudonana; functional characterization of enzyme catalyzing three-step methylation of glycine. Plant Physiol Biochem. 2018;127:248–55.CAS 
PubMed 

Google Scholar 
Krell A, Funck D, Plettner I, John U, Dieckmann G. Regulation of proline metabolism under salt stress in the psychrophilic diatom Fragilariopsis cylindrus (Bacillariophyceae). J Phycol. 2007;43:753–62.CAS 

Google Scholar 
Latta LC, Weider LJ, Colbourne JK, Pfrender ME. The evolution of salinity tolerance in Daphnia: a functional genomics approach. Ecol Lett. 2012;15:794–802.PubMed 

Google Scholar 
Ferrante MI, Entrambasaguas L, Johansson M, Töpel M, Kremp A, Montresor M, et al. Exploring molecular signs of sex in the marine diatom Skeletonema marinoi. Genes. 2019;10:494.Kroth PG. The biodiversity of carbon assimilation. J Plant Physiol. 2015;172:76–81.CAS 
PubMed 

Google Scholar 
Obata T, Fernie AR, Nunes-Nesi A. The central carbon and energy metabolism of marine diatoms. Metabolites. 2013;3:325–46.CAS 
PubMed 
PubMed Central 

Google Scholar 
Smith SR, Abbriano RM, Hildebrand M. Comparative analysis of diatom genomes reveals substantial differences in the organization of carbon partitioning pathways. Algal Res. 2012;1:2–16.CAS 

Google Scholar 
Kroth PG, Chiovitti A, Gruber A, Martin-Jezequel V, Mock T, Parker MS, et al. A model for carbohydrate metabolism in the diatom Phaeodactylum tricornutum deduced from comparative whole genome analysis. PLoS ONE. 2008;3:e1426.PubMed 
PubMed Central 

Google Scholar 
Furumoto T, Yamaguchi T, Ohshima-Ichie Y, Nakamura M, Tsuchida-Iwata Y, Shimamura M, et al. A plastidial sodium-dependent pyruvate transporter. Nature. 2011;476:472–5.CAS 
PubMed 

Google Scholar 
Chen G-Q, Jiang Y, Chen F. Salt-induced alterations in lipid composition of diatom Nitzschia laevis (Bacillariophyceae) under heterotrophic culture condition. J Phycol. 2008;44:1309–14.CAS 
PubMed 

Google Scholar 
Sayanova O, Mimouni V, Ulmann L, Morant-Manceau A, Pasquet V, Schoefs B, et al. Modulation of lipid biosynthesis by stress in diatoms. Philos Trans R Soc Lond B Biol Sci. 2017;372:20160407.PubMed 
PubMed Central 

Google Scholar 
Vårum KM, Myklestad S. Effects of light, salinity and nutrient limitation on the production of β-1,3-d-glucan and exo-d-glucanase activity in Skeletonema costatum (Grev.) Cleve. J Exp Mar Bio Ecol. 1984;83:13–25.
Google Scholar 
Radchenko IG, Il’yash LV. Growth and photosynthetic activity of diatom Thalassiosira weissflogii at decreasing salinity. Biol Bull. 2006;33:242–7.CAS 

Google Scholar 
Adams C, Bugbee B. Enhancing lipid production of the marine diatom Chaetoceros gracilis: synergistic interactions of sodium chloride and silicon. J Appl Phycol. 2014;26:1351–7.CAS 

Google Scholar 
Shetty P, Gitau MM, Maróti G. Salinity stress responses and adaptation mechanisms in eukaryotic green microalgae. Cells. 2019;8:1657.Jacob A, Kirst GO, Wiencke C, Lehmann H. Physiological responses of the Antarctic green alga Prasiola crispa ssp. antarctica to salinity stress. J Plant Physiol. 1991;139:57–62.CAS 

Google Scholar 
Bazzani E, Lauritano C, Mangoni O, Bolinesi F, Saggiomo M. Chlamydomonas responses to salinity stress and possible biotechnological exploitation. J Mar Sci Eng. 2021;9:1242.
Google Scholar 
Cheng R-L, Feng J, Zhang B-X, Huang Y, Cheng J, Zhang C-X. Transcriptome and gene expression analysis of an oleaginous diatom under different salinity conditions. Bioenergy Res. 2014;7:192–205.CAS 

Google Scholar 
Stock W, Blommaert L, Daveloose I, Vyverman W, Sabbe K. Assessing the suitability of imaging-PAM fluorometry for monitoring growth of benthic diatoms. J Exp Mar Bio Ecol. 2019;513:35–41.
Google Scholar 
Reichmann D, Voth W, Jakob U. Maintaining a healthy proteome during oxidative stress. Mol Cell. 2018;69:203–13.CAS 
PubMed 
PubMed Central 

Google Scholar 
Latowski D, Kuczyńska P, Strzałka K. Xanthophyll cycle-a mechanism protecting plants against oxidative stress. Redox Rep. 2011;16:78–90.CAS 
PubMed 

Google Scholar 
Chen D, Shao Q, Yin L, Younis A, Zheng B. Polyamine function in plants: metabolism, regulation on development, and roles in abiotic stress responses. Front Plant Sci. 2018;9:1945.PubMed 

Google Scholar 
Liu Q, Nishibori N, Imai I, Hollibaugh JT. Response of polyamine pools in marine phytoplankton to nutrient limitation and variation in temperature and salinity. Mar Ecol Prog Ser. 2016;544:93–105.CAS 

Google Scholar 
Scoccianti V, Penna A, Penna N, Magnani M. Effect of heat stress on polyamine content and protein pattern in Skeletonema costatum. Mar Biol. 1995;121:549–54.CAS 

Google Scholar 
Alscher RG, Erturk N, Heath LS. Role of superoxide dismutases (SODs) in controlling oxidative stress in plants. J Exp Bot. 2002;53:1331–41.CAS 
PubMed 

Google Scholar 
Kumar M, Kumari P, Gupta V, Reddy CRK, Jha B. Biochemical responses of red alga Gracilaria corticata (Gracilariales, Rhodophyta) to salinity induced oxidative stress. J Exp Mar Bio Ecol. 2010;391:27–34.CAS 

Google Scholar 
von Alvensleben N, Magnusson M, Heimann K. Salinity tolerance of four freshwater microalgal species and the effects of salinity and nutrient limitation on biochemical profiles. J Appl Phycol. 2016;28:861–76.
Google Scholar 
Rijstenbil JW, Wijnholds JA, Sinke JJ. Implications of salinity fluctuation for growth and nitrogen metabolism of the marine diatom Ditylum brightwellii in comparison with Skeletonema costatum. Mar Biol. 1989;101:131–41.CAS 

Google Scholar 
Mansour MMF. Nitrogen containing compounds and adaptation of plants to salinity stress. Biol Plant. 2000;43:491–500.CAS 

Google Scholar 
Garcia N, Lopez Elias JA, Miranda A, Martinez Porchas M, Huerta N, Garcia A. Effect of salinity on growth and chemical composition of the diatom Thalassiosira weissflogii at three culture phases. Lat Am J Aquat Res. 2012;40:435–40.
Google Scholar 
Van den Berge K, Hembach KM, Soneson C, Tiberi S, Clement L, Love MI, et al. RNA sequencing data: Hitchhiker’s guide to expression analysis. Annu Rev Biomed Data Sci. 2019;2:139–73.
Google Scholar 
Kremp A. Effects of cyst resuspension on germination and seeding of two bloom-forming dinoflagellates in the Baltic Sea. Mar Ecol Prog Ser. 2001;216:57–66.
Google Scholar 
Juneau P, Barnett A, Méléder V, Dupuy C, Lavaud J. Combined effect of high light and high salinity on the regulation of photosynthesis in three diatom species belonging to the main growth forms of intertidal flat inhabiting microphytobenthos. J Exp Mar Bio Ecol. 2015;463:95–104.CAS 

Google Scholar 
Vargas C, Argandoña M, Reina-Bueno M, Rodríguez-Moya J, Fernández-Aunión C, Nieto JJ. Unravelling the adaptation responses to osmotic and temperature stress in Chromohalobacter salexigens, a bacterium with broad salinity tolerance. Saline Syst. 2008;4:14.PubMed 
PubMed Central 

Google Scholar 
Khmelenina VN, Sakharovskii VG, Reshetnikov AS, Trotsenko YA. Synthesis of osmoprotectants by halophilic and alkaliphilic methanotrophs. Microbiology. 2000;69:381–6.CAS 

Google Scholar 
Fenizia S, Thume K, Wirgenings M, Pohnert G. Ectoine from bacterial and algal origin is a compatible solute in microalgae. Mar Drugs. 2020;18:42.CAS 
PubMed Central 

Google Scholar 
Amin SA, Hmelo LR, van Tol HM, Durham BP, Carlson LT, Heal KR, et al. Interaction and signalling between a cosmopolitan phytoplankton and associated bacteria. Nature. 2015;522:98–101.CAS 
PubMed 

Google Scholar 
Krell A, Beszteri B, Dieckmann G, Glöckner G, Valentin K, Mock T. A new class of ice-binding proteins discovered in a salt-stress-induced cDNA library of the psychrophilic diatom Fragilariopsis cylindrus (Bacillariophyceae). Eur J Phycol. 2008;43:423–33.CAS 

Google Scholar 
Helliwell KE, Kleiner FH, Hardstaff H, Chrachri A, Gaikwad T, Salmon D, et al. Spatiotemporal patterns of intracellular Ca2+ signalling govern hypo-osmotic stress resilience in marine diatoms. N Phytol. 2021;230:155–70.CAS 

Google Scholar 
Kaczmarska I, Poulíčková A, Sato S, Edlund MB, Idei M, Watanabe T, et al. Proposals for a terminology for diatom sexual reproduction, auxospores and resting stages. Diatom Res. 2013;28:263–94.
Google Scholar 
Godhe A, Kremp A, Montresor M. Genetic and microscopic evidence for sexual reproduction in the centric diatom Skeletonema marinoi. Protist. 2014;165:401–16.PubMed 

Google Scholar 
Annunziata R, Mele BH, Marotta P, Volpe M, Entrambasaguas L, Mager S, et al. Trade-off between sex and growth in diatoms: Molecular mechanisms and demographic implications. Sci Adv. 2022;8:eabj9466.CAS 
PubMed 
PubMed Central 

Google Scholar 
Ajani PA, Petrou K, Larsson ME, Nielsen DA, Burke J, Murray SA. Phenotypic trait variability as an indication of adaptive capacity in a cosmopolitan marine diatom. Environ Microbiol. 2021;23:207–23.CAS 
PubMed 

Google Scholar 
Sjöqvist CO, Kremp A. Genetic diversity affects ecological performance and stress response of marine diatom populations. ISME J. 2016;10:2755–66.PubMed 
PubMed Central 

Google Scholar 
Godhe A, Rynearson T. The role of intraspecific variation in the ecological and evolutionary success of diatoms in changing environments. Philos Trans R Soc Lond B Biol Sci. 2017;372:20160399.PubMed 
PubMed Central 

Google Scholar 
Bulankova P, Sekulić M, Jallet D, Nef C, van Oosterhout C, Delmont TO, et al. Mitotic recombination between homologous chromosomes drives genomic diversity in diatoms. Curr Biol. 2021;31:3221–32. e9CAS 
PubMed 

Google Scholar 
Pinseel E, Janssens SB, Verleyen E, Vanormelingen P, Kohler TJ, Biersma EM, et al. Global radiation in a rare biosphere soil diatom. Nat Commun. 2020;11:2382.CAS 
PubMed 
PubMed Central 

Google Scholar 
Savchuk OP. Large-scale nutrient dynamics in the Baltic sea, 1970–2016. Front Mar Sci. 2018;5:95.
Google Scholar 
Gomez-Mestre I, Jovani R. A heuristic model on the role of plasticity in adaptive evolution: plasticity increases adaptation, population viability and genetic variation. Proc Biol Sci. 2013;280:20131869.PubMed 
PubMed Central 

Google Scholar 
Lambert BS, Groussman RD, Schatz MJ, Coesel SN, Durham BP, Alverson AJ, et al. The dynamic trophic architecture of open-ocean protist communities revealed through machine-guided metatranscriptomics. Proc Natl Acad Sci USA. 2022;119:e2100916119.Harrison PF, Pattison AD, Powell DR, Beilharz TH. Topconfects: a package for confident effect sizes in differential expression analysis provides a more biologically useful ranked gene list. Genome Biol. 2019;20:67.PubMed 
PubMed Central 

Google Scholar  More

Cury, P. et al. Small pelagics in upwelling systems: patterns of interaction and structural changes in “wasp-waist” ecosystems. ICES J. Mar. Sci. 57, 603–618 (2000).
Google Scholar 
Fréon, P., Cury, P., Shannon, L. & Roy, C. Sustainable exploitation of small pelagic fish stocks challenged by environmental and ecosystem changes: a review. Bull. Mar. Sci. 76, 385–462 (2005).
Google Scholar 
The State of Mediterranean and Black Sea Fisheries 2020. (FAO, 2020). doi:https://doi.org/10.4060/cb2429en.Checkley, D. M., Asch, R. G. & Rykaczewski, R. R. Climate, Anchovy, and Sardine. Annu. Rev. Mar. Sci. 9, 469–493 (2017).ADS 

Google Scholar 
Bakun, A. Patterns in the ocean: Ocean processes and marine population dynamics. Oceanogr. Lit. Rev. 5, 530 (1997).
Google Scholar 
Houde, E. D. Recruitment variability. Fish. Reprod. Biol. 1, 91–171 (2009).
Google Scholar 
Cole, J. & McGlade, J. Clupeoid population variability, the environment and satellite imagery in coastal upwelling systems. Rev. Fish Biol. Fish. 8, 445–471 (1998).
Google Scholar 
Pörtner, H. O. & Peck, M. A. Climate change effects on fishes and fisheries: towards a cause-and-effect understanding. J. Fish Biol. 77, 1745–1779 (2010).PubMed 

Google Scholar 
European Commission. Joint Research Centre. Institute for the Protection and the Security of the Citizen. Scientific, technical and economic committee for fisheries (STECF) : Assessment of mediterranean stocks (Part I). (Publications Office, 2010).García, A. & Palomera, I. Anchovy early life history and its relation to its surrounding environment in the Western Mediterranean basin. 12 (1996).Tsikliras, A. C., Antonopoulou, E. & Stergiou, K. I. Spawning period of Mediterranean marine fishes. 40 (2010).Fernández-Corredor, E., Albo-Puigserver, M., Pennino, M. G., Bellido, J. M. & Coll, M. Influence of environmental factors on different life stages of European anchovy (Engraulis encrasicolus) and European sardine (Sardina pilchardus) from the Mediterranean Sea: A literature review. Reg. Stud. Mar. Sci. 41, 1006 (2021).
Google Scholar 
Agostini, V. N. & Bakun, A. Ocean triads’ in the Mediterranean Sea: physical mechanisms potentially structuring reproductive habitat suitability (with example application to European anchovy, Engraulis encrasicolus). Fish. Oceanogr. 11, 129–142 (2002).
Google Scholar 
Lasker, R. The relation between oceanographic conditions and larval anchovy food in the California Current: identification of factors contributing to recruitment failure. Rapp. P-V Reun. Cons. Int. Explor. Mer. 173, 212–230 (1978).
Google Scholar 
Cushing, D. H. Plankton Production and Year-class Strength in Fish Populations: an Update of the Match/Mismatch Hypothesis. in Advances in Marine Biology (eds. Blaxter, J. H. S. & Southward, A. J.) vol. 26 249–293 (Academic Press, 1990).Patti, B. et al. Anchovy (Engraulis encrasicolus) early life stages in the Central Mediterranean Sea: connectivity issues emerging among adjacent sub-areas across the Strait of Sicily. Hydrobiologia 821, 25–40 (2018).
Google Scholar 
Maynou, F., Sabatés, A. & Raya, V. Changes in the spawning habitat of two small pelagic fish in the Northwestern Mediterranean. Fish. Oceanogr. 29, 201–213 (2020).
Google Scholar 
Basilone, G. et al. Spawning site selection by European anchovy ( Engraulis encrasicolus ) in relation to oceanographic conditions in the Strait of Sicily. Fish. Oceanogr. 22, 309–323 (2013).
Google Scholar 
Malavolti, S. et al. Distribution of Engraulis encrasicolus eggs and larvae in relation to coastal oceanographic conditions (the South-western Adriatic Sea case study). Mediterr. Mar. Sci. 19, 180 (2018).
Google Scholar 
Zorica, B. et al. Spatiotemporal distribution of anchovy early life stages in the eastern part of the Adriatic Sea in relation to some oceanographic features. J. Mar. Biol. Assoc. UK 99, 1205–1211 (2019).
Google Scholar 
Punt, A. E. et al. Fisheries management under climate and environmental uncertainty: control rules and performance simulation. ICES J. Mar. Sci. 71, 2208–2220 (2014).
Google Scholar 
Szuwalski, C. S. & Hilborn, R. Environment drives forage fish productivity. Proc. Natl. Acad. Sci. 112, E3314–E3315 (2015).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ganias, K. Biology and Ecology of Sardines and Anchovies. 391.Palomera, I. et al. Small pelagic fish in the NW Mediterranean Sea: An ecological review. Prog. Oceanogr. 74, 377–396 (2007).ADS 

Google Scholar 
Patti, B., Torri, M. & Cuttitta, A. General surface circulation controls the interannual fluctuations of anchovy stock biomass in the Central Mediterranean Sea. Sci. Rep. 10, 1554 (2020).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Pisano, A. et al. New evidence of mediterranean climate change and variability from sea surface temperature observations. Rem. Sens. 12, 132 (2020).ADS 

Google Scholar 
Regner, S. Effects of environmental changes on early stages and reproduction of anchovy in the Adriatic Sea. Sci. Mar. 60, 167–177 (1996).
Google Scholar 
Lavigne, H. et al. Enhancing the comprehension of mixed layer depth control on the Mediterranean phytoplankton phenology: Mediterranean Phytoplankton Phenology. J. Geophys. Res. Oceans 118, 3416–3430 (2013).ADS 

Google Scholar 
Morote, E., Olivar, M. P., Villate, F. & Uriarte, I. A comparison of anchovy (Engraulis encrasicolus) and sardine (Sardina pilchardus) larvae feeding in the Northwest Mediterranean: Influence of prey availability and ontogeny. ICES J. Mar. Sci. 67, 897–908 (2010).
Google Scholar 
Basilone, G. et al. Spawning ecology of the European anchovy (Engraulis encrasicolus) in the Strait of Sicily: Linking variations of zooplankton prey, fish density, growth, and reproduction in an upwelling system. Prog. Oceanogr. 184, 102330 (2020).
Google Scholar 
McBride, R. S. et al. Energy acquisition and allocation to egg production in relation to fish reproductive strategies. Fish Fish. 16, 23–57 (2015).
Google Scholar 
d’Ortenzio, F. & Ribera d’Alcalà, M. On the trophic regimes of the Mediterranean Sea: a satellite analysis. Biogeosciences 6, 139–148 (2009).ADS 

Google Scholar 
Escribano, A., Aldanondo, N., Cotano, U., Boyra, G. & Urtizberea, A. Size- and density-dependent overwinter mortality of anchovy juveniles in the Bay of Biscay. Cont. Shelf Res. 183, 28–37 (2019).ADS 

Google Scholar 
Bellido, J. M. et al. Identifying essential fish habitat for small pelagic species in Spanish Mediterranean waters. Hydrobiologia 612, 171–184 (2008).
Google Scholar 
Van Beveren, E. et al. The fisheries history of small pelagics in the Northern Mediterranean. ICES J. Mar. Sci. 73, 1474–1484 (2016).
Google Scholar 
Edwards, M., Beaugrand, G., Helaouët, P., Alheit, J. & Coombs, S. Marine ecosystem response to the Atlantic Multidecadal Oscillation. PLoS ONE 8, e57212 (2013).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Erauskin-Extramiana, M. et al. Historical trends and future distribution of anchovy spawning in the Bay of Biscay. Deep Sea Res. Part II Top. Stud. Oceanogr. 159, 169–182 (2019).
Google Scholar 
Bertrand, A. et al. Schooling behaviour and environmental forcing in relation to anchoveta distribution: An analysis across multiple spatial scales. Prog. Oceanogr. 79, 264–277 (2008).ADS 

Google Scholar 
Giannoulaki, M. et al. Characterizing the potential habitat of European anchovy Engraulis encrasicolus in the Mediterranean Sea, at different life stages: Habitat of anchovy in the Mediterranean. Fish. Oceanogr. 22, 69–89 (2013).
Google Scholar 
Durant, J., Hjermann, D., Ottersen, G. & Stenseth, N. Climate and the match or mismatch between predator requirements and resource availability. Clim. Res. 33, 271–283 (2007).
Google Scholar 
Kharouba, H. M. & Wolkovich, E. M. Disconnects between ecological theory and data in phenological mismatch research. Nat. Clim. Change 10, 406–415 (2020).ADS 

Google Scholar 
Kodama, T. et al. Improvement in recruitment of Japanese sardine with delays of the spring phytoplankton bloom in the Sea of Japan. Fish. Oceanogr. 27, 289–301 (2018).
Google Scholar 
Platt, T., Fuentes-Yaco, C. & Frank, K. T. Spring algal bloom and larval fish survival. Nature 423, 398–399 (2003).ADS 
CAS 
PubMed 

Google Scholar 
Schweigert, J. F. et al. Factors linking Pacific herring (Clupea pallasi) productivity and the spring plankton bloom in the Strait of Georgia, British Columbia. Canada. Prog. Oceanogr. 115, 103–110 (2013).ADS 

Google Scholar 
Laurel, B. J. et al. Regional warming exacerbates match/mismatch vulnerability for cod larvae in Alaska. Progress Oceanogr. 193, 1055 (2021).
Google Scholar 
Brosset, P. et al. A fine-scale multi-step approach to understand fish recruitment variability. Sci. Rep. 10, 16064 (2020).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Ferreira, A. et al. Match-mismatch dynamics in the Norwegian-Barents Sea system. Mar. Ecol. Prog. Ser. 650, 81–94 (2020).ADS 

Google Scholar 
Ferreira, A. S. A., Hátún, H., Counillon, F., Payne, M. R. & Visser, A. W. Synoptic-scale analysis of mechanisms driving surface chlorophyll dynamics in the North Atlantic. Biogeosciences 12, 3641–3653 (2015).ADS 

Google Scholar 
Taboada, F. G. & Anadón, R. Determining the causes behind the collapse of a small pelagic fishery using Bayesian population modeling. Ecol. Appl. 26, 886–898 (2016).
Google Scholar 
Katara, I., Pierce, G. J., Illian, J. & Scott, B. E. Environmental drivers of the anchovy/sardine complex in the Eastern Mediterranean. Hydrobiologia 670, 49–65 (2011).
Google Scholar 
Robinson, A. R. et al. The Atlantic Ionian Stream. J. Mar. Syst. 20, 129–156 (1999).
Google Scholar 
Di Lorenzo, M., Sinerchia, M. & Colloca, F. The North sector of the Strait of Sicily: A priority area for conservation in the Mediterranean Sea. Hydrobiologia 821, 235–253 (2018).
Google Scholar 
García Lafuente, J. et al. Hydrographic phenomena influencing early life stages of the Sicilian Channel anchovy: Hydrographic relations with the Sicilian Channel anchovy. Fish. Oceanogr. 11, 31–44 (2002).
Google Scholar 
Patti, B. et al. Role of physical forcings and nutrient availability on the control of satellite-based chlorophyll a concentration in the coastal upwelling area of the Sicilian Channel. Sci. Mar. 74, 577–588 (2010).
Google Scholar 
Marini, M., Jones, B. H., Campanelli, A., Grilli, F. & Lee, C. M. Seasonal variability and Po River plume influence on biochemical properties along western Adriatic coast. J. Geophys. Res. 113, C0590 (2008).
Google Scholar 
Artegiani et al. The Adriatic Sea General Circulation. Part I Air–.pdf. (1997)Artegiani, A. et al. The Adriatic Sea general circulation. Part II: Baroclinic circulation structure. J. Phys. Oceanogr. 27, 18 (1997).
Google Scholar 
Millot, C. The Gulf of Lions’ hydrodynamics. Cont. Shelf Res. 10, 885–894 (1990).ADS 

Google Scholar 
Petrenko, A., Leredde, Y. & Marsaleix, P. Circulation in a stratified and wind-forced Gulf of Lions, NW Mediterranean Sea: In situ and modeling data. Cont. Shelf Res. 25, 7–27 (2005).ADS 

Google Scholar 
Sverdrup, H. U. On conditions for the vernal blooming of phytoplankton. J. Cons. Int. Explor. Mer. 18, 287–295 (1953).
Google Scholar 
Wilcox, R. R. Introduction to robust estimation and hypothesis testing. (Elsevier, 2016).Pernet, C. R., Wilcox, R. R. & Rousselet, G. A. Robust correlation analyses: False positive and power validation using a new open source Matlab toolbox. Front. Psychol. 1, 1 (2013).
Google Scholar 
Mair, P. & Wilcox, R. Robust statistical methods in R using the WRS2 package. Behav. Res. Methods 52, 464–488 (2020).PubMed 

Google Scholar 
Fox, J., & Weisberg, S. “Robust regression.” An R and S-Plus companion to applied regression 91 (2002).Crawley, M. J. The R book. John Wiley & Sons (2012).Fox, J. Bootstrapping Regression Models Appendix to An R and S-PLUS Companion to Applied Regression (2002). More

Seasonal conditions during crop growthThe groundnut crop produces optimum yield in the regions receiving rainfall between 200 to 1000 mm8. The total rainfall during groundnut growing season was 257.10 mm and 403.10 mm at Baljigapade in 2018 and 2019, respectively, wherein Pavagada (2018) total rainfall was 53.90 mm (Fig. 1). In 2018, both Pavagada and Baljigapade received very low and negligible rainfall during the reproduction and harvest stage of groundnut. Optimum temperature for groundnut production ranges between 20 to 30 °C and growth and pod formation limited below 16 °C and above 32 °C9. The monthly mean atmospheric temperature was ranged from 23 to 27 °C and 21 to 26 °C at Baljigapade in 2018 and 2019, respectively, wherein Pavagada ranged from 25 to 28 °C. At all three locations, the monthly mean atmospheric temperature was slightly high during the early vegetative growth of groundnut and it was progressively decreased as the crop reaches its maturity stage (Fig. 1). All three locations recorded higher and lower mean monthly sunshine hours (hours day−1) during peg initiation to pod filling stage (September and October) and early vegetative growth of groundnut (July and August), respectively.Growth parameters of groundnutAnalysis of variance revealed that treatment, location, and their interaction had a significant effect on plant height and number of branches at harvest (P  More

Download PDF

I’m from a fishing family. My grandfather was a fisherman when he was a young man, working out of Fano, the Italian town where I grew up and still live. I’m used to the smell of fish.I’m pictured during an overnight shift on the fishing boat RIMAS. I work from 5 p.m. until 9 a.m. with fishermen from nearby Cesenatico on the north Adriatic Sea. It’s a small boat: there’s only six or so of us on board. At night, the fish are most active and we can avoid other vessels.The nets scrape the sea bed for the catch but sometimes they also catch turtles who often die in the nets or on board. That’s where I come in. The net I’m holding is designed to allow turtles to escape: it has a hole at the top they can swim out of. We call it TED — short for ‘turtle excluder device’. The TED is made from a high-strength plastic, and is based on decades of work and research aimed at reducing the bycatch of turtles from trawling. Turtles and some larger fish can leave through the escape hatch, but the current holds most of the catch in the net.I ensure that the net is working, and that the fishermen we’re collaborating with can still catch enough for their livelihoods while protecting turtles. The work is part of research by the Cetacea Foundation, based in Riccione, Italy, in collaboration with the University of Pisa, where I’m a field researcher. It is financed by the LIFE programme, the European Union’s funding instrument for the environment and climate action.I love this work. It means I’m not stuck in an office all day and instead can enjoy the ocean and work closely with people who live by the sea. I get to be a researcher who works outside, rather than being hunched over a microscope.When my grandfather was fishing in the 1970s, there were more fish and more turtles around. At the foundation, we save 50–60 turtles a year, most of them harmed because of fishing. If we can protect turtles by rolling out this device to fishermen all across the Adriatic, I’d see this work as a success.

Nature 604, 210 (2022)
doi: https://doi.org/10.1038/d41586-022-00930-w

Related Articles

Protector of giant salamander

A photo celebration of scientists at work

Gut feeling: building a picture of Latin American microbiomes

Subjects

Careers

Climate change

Conservation biology

Latest on:

Careers

How and why to say ‘no’ to colleagues and collaborators
Career Column 28 MAR 22

Afghanistan’s girls’ schools can — and must — stay open. There is no alternative
Editorial 28 MAR 22

The marine biologist whose photography pastime became a profession
Career Column 25 MAR 22

Climate change

Tropical forests have big climate benefits beyond carbon storage
News 01 APR 22

Funding battles stymie ambitious plan to protect global biodiversity
News 31 MAR 22

Trends in Europe storm surge extremes match the rate of sea-level rise
Article 30 MAR 22

Jobs

Postdoctoral Fellow

LSU Health Sciences Center New Orleans (LSUHSC-NO)
New Orleans, LA, United States

Associate Research Scientist

Columbia University Medical Center (CUMC), CU
New York, NY, United States

Experimental Postdoctoral Research Associate in Optical/Quantum Metasurfaces

Los Alamos National Laboratory (LANL)
Los Alamos, NM, United States

Theory Postdoctoral Research Associate in Optical/Quantum Metasurfaces

Los Alamos National Laboratory (LANL)
Los Alamos, NM, United States More

Mabhaudhi, T. et al. Prospects of orphan crops in climate change. Planta 250, 695–708. https://doi.org/10.1007/s00425-019-03129-y (2019).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Singh, M., Bisht, I. S., Dutta, M., Springer. India, 221. https://doi.org/10.1007/978-81-322-2023-7(2014).Litrico, I. & Violle, C. Diversity in plant breeding: A new conceptual framework. Trends Plant Sci. 20, 604–613. https://doi.org/10.1016/j.tplants.2015.07.007 (2015).CAS 
Article 
PubMed 

Google Scholar 
Govindaraj, M., Vetriventhan, M. & Srinivasan, M. Importance of genetic diversity assessment in crop plants and its recent advances: An overview of its analytical perspectives. Genet. Res. Intern. 431–487, 2015. https://doi.org/10.1155/2015/431487 (2015).Article 

Google Scholar 
Akohoué, F., Sibiya, J. & Achigan-Dako, E. G. On-farm practices, mapping, and uses of genetic resources of Kersting’s groundnut [Macrotyloma geocarpum (Harms) Maréchal et Baudet] across ecological zones in Benin and Togo. Genet. Resour. Crop. Evol. 66, 195–214. https://doi.org/10.1007/s10722-018-0705-7 (2018).CAS 
Article 

Google Scholar 
Assogba, P. et al. Indigenous knowledge and agro-morphological evaluation of the minor crop Kersting’s groundnut [Macrotyloma geocarpum (Harms) Maréchal et Baudet] cultivars of Benin. Genet. Resour. Crop. Evol. 63, 513–529. https://doi.org/10.1007/s10722-015-0268-9 (2015).Article 

Google Scholar 
Adu-Gyamfi, R., Fearon, J., Bayorbor, T. B., Dzomeku, I. K. & Avornyo, V. K. The Status of Kersting’s groundnut [Macrotyloma geocarpum (Harms) Marechal and Baudet]. Outlook Agric. 40, 259–262. https://doi.org/10.5367/oa.2011.0050 (2011).Article 

Google Scholar 
Obasi, M. O. & Agbatse, A. Evaluation of nutritive value and some functional properties of Kersting’s groundnut seeds for optimum utilisation as a food and feed resource. E. Afr. Agric. For. J. 68, 173–181. https://doi.org/10.4314/eaafj.v68i4.1794 (2003).Article 

Google Scholar 
Ajayi, O. B. & Oyetayo, F. L. Potentials of Kerstingiella geocarpa as a health food. J. Med. Food 12, 184–187. https://doi.org/10.1089/jmf.2008.0100 (2009).CAS 
Article 
PubMed 

Google Scholar 
Mohammed, M., Jaiswal, S. K., Sowley, E. N. K., Ahiabor, B. D. K. & Dakora, F. D. Symbiotic N2 fixation and grain yield of endangered Kersting’s groundnut landraces in response to soil and plant associated bradyrhizobium inoculation to promote ecological resource-use efficiency. Front. Microbiol. 9, 1–14. https://doi.org/10.3389/fmicb.2018.02105 (2018).CAS 
Article 

Google Scholar 
Tamini, Z. Étude ethnobotanique de la Lentille de Terre [Macrotyloma geocarpum Maréchal & Baudet] au Burkina Faso. J. Agric. Trad. Bot. Appl. 37, 187–199. https://doi.org/10.3406/jatba.1995.3569 (1995).Article 

Google Scholar 
Coulibaly, M., Agossou, C. O. A., Akohoué, F., Sawadogo, M. & Achigan-Dako, E. G. Farmers’ preferences for genetic resources of Kersting’s groundnut [Macrotyloma geocarpum (Harms) Maréchal and Baudet] in the production systems of Burkina Faso and Ghana. Agronomy 10, 1–20. https://doi.org/10.3390/agronomy10030371 (2020).CAS 
Article 

Google Scholar 
AchiganDako, E. G. & Vodouhe, S. R. Macrotyloma geocarpum (Harms) Marechal & Baudet. In Plant Resources of Tropical Africa 1: Cereals and Pulses (ed. Brink, M. B. G.) 111–114 (Backhuys Publishers CTA, PROTA, 2006).
Google Scholar 
Mergeai, G. Influence des facteurs sociologiques sur la conservation des ressources phytogenetiques: Le cas de la lentille de terre [Macrotyloma geocarpum (Harms) Marechal et Baudet] au Togo. Bull Rech Agron 28, 487–500 (1993).
Google Scholar 
Long, S. P., Marshall-Colon, A. & Zhu, X. G. Meeting the global food demand of the future by engineering crop photosynthesis and yield potential. Cell 161, 56–66. https://doi.org/10.1016/j.cell.2015.03.019 (2015).CAS 
Article 
PubMed 

Google Scholar 
Akohoue, F., Achigan-Dako, E. G., Sneller, C., Van Deynze, A. & Sibiya, J. Genetic diversity, SNP-trait associations and genomic selection accuracy in a west African collection of Kersting’s groundnut [Macrotyloma geocarpum (Harms) Marechal & Baudet]. PLoS ONE 15, 1–24. https://doi.org/10.1371/journal.pone.0234769 (2020).CAS 
Article 

Google Scholar 
Schierenbeck, K. A. Population-level genetic variation and climate change in a biodiversity hotspot. Ann. Bot. 119, 215–228. https://doi.org/10.1093/aob/mcw214 (2017).Article 
PubMed 

Google Scholar 
Araújo, M. B., Whittaker, R. J., Ladle, R. J. & Erhard, M. Reducing uncertainty in projections of extinction risk from climate change. Glob. Ecol. Biogeogr. 14, 529–538. https://doi.org/10.1111/j.1466-822X.2005.00182.x (2005).Article 

Google Scholar 
Araujo, M. B. & Peterson, A. T. Uses and misuses of bioclimatic envelope modeling. Ecology 93, 1527–1539. https://doi.org/10.1890/11-1930.1 (2012).Article 
PubMed 

Google Scholar 
Martínez-Meyer, E. Climate change and biodiversity: Some considerations in forecasting shifts in species’ potential distributions. Biodiv. Inform. 2, 42–55. https://doi.org/10.17161/bi.v2i0.8 (2005).Article 

Google Scholar 
Elith, J. et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography 29, 129–151 (2006).Article 

Google Scholar 
Pironon, S. et al. Potential adaptive strategies for 29 sub-Saharan crops under future climate change. Nat. Clim. Chang. 9, 758–763. https://doi.org/10.1038/s41558-019-0585-7 (2019).ADS 
Article 

Google Scholar 
Ramirez-Cabral, N. Y. Z., Kumar, L. & Taylor, S. Crop niche modeling projects major shifts in common bean growing areas. Agric. For. Meteor. 218–219, 102–113. https://doi.org/10.1016/j.agrformet.2015.12.002 (2016).Article 

Google Scholar 
Syfert, M. M. et al. Crop wild relatives of the brinjal eggplant (Solanum melongena): Poorly represented in genebanks and many species at risk of extinction. Am. J. Bot. 103, 1–17. https://doi.org/10.3732/ajb.1500539 (2016).CAS 
Article 

Google Scholar 
Blonder, B. Hypervolume concepts in niche- and trait-based ecology. Ecography 41, 1441–1455. https://doi.org/10.1111/ecog.03187 (2018).Article 

Google Scholar 
Hampe, A. & Petit, R. J. Conserving biodiversity under climate change: The rear edge matters. Ecol. Lett. 8, 461–467. https://doi.org/10.1111/j.1461-0248.2005.00739.x (2005).Article 
PubMed 

Google Scholar 
Leimu, R. & Fischer, M. A meta-analysis of local adaptation in plants. PLoS ONE 3, 1–8. https://doi.org/10.1371/journal.pone.0004010 (2008).CAS 
Article 

Google Scholar 
Hereford, J. A quantitative survey of local adaptation and fitness trade-offs. Am. Naturalist. 173, 579–588. https://doi.org/10.1086/597611 (2009).Article 

Google Scholar 
Shaw, R. G. & Etterson, J. R. Rapid climate change and the rate of adaptation: Insight from experimental quantitative genetics. New Phytol. 195, 752–765. https://doi.org/10.1111/j.1469-8137.2012.04230.x (2012).Article 
PubMed 

Google Scholar 
Gotelli, N. J. & Stanton-Geddes, J. Climate change, genetic markers and species distribution modelling. J. Biogeogr. 42, 1577–1585. https://doi.org/10.1111/jbi.12562 (2015).Article 

Google Scholar 
Ikeda, D. H. et al. Genetically informed ecological niche models improve climate change predictions. Glob. Change Biol. 23, 164–176. https://doi.org/10.1111/gcb.13470 (2016).ADS 
Article 

Google Scholar 
Alvarado-Serrano, D. F. & Knowles, L. L. Ecological niche models in phylogeographic studies: Applications, advances and precautions. Mol. Ecol. Resour. 14, 233–248. https://doi.org/10.1111/1755-0998.12184 (2014).Article 
PubMed 

Google Scholar 
Thoen, M. P. et al. Genetic architecture of plant stress resistance: Multi-trait genome-wide association mapping. New Phytol. 213, 1346–1362. https://doi.org/10.1111/nph.14220 (2017).CAS 
Article 
PubMed 

Google Scholar 
Kafoutchoni, K. M., Agoyi, E. E., Agbahoungba, S., Assogbadjo, A. E. & Agbangla, C. Genetic diversity and population structure in a regional collection of Kersting’s groundnut [Macrotyloma geocarpum (Harms) Maréchal & Baudet]. Genet. Resour. Crop. Evol. https://doi.org/10.1007/s10722-021-01187-4 (2021).Article 

Google Scholar 
Brown, J. L. & Carnaval, A. C. A tale of two niches: Methods, concepts, and evolution. Front. Biogeogr. https://doi.org/10.21425/f5fbg44158 (2019).Article 

Google Scholar 
Marcer, A., Mendez-Vigo, B., Alonso-Blanco, C. & Pico, F. X. Tackling intraspecific genetic structure in distribution models better reflects species geographical range. Ecol. Evol. 6, 2084–2097. https://doi.org/10.1002/ece3.2010 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Oney, B., Reineking, B., O’Neill, G. & Kreyling, J. Intraspecific variation buffers projected climate change impacts on Pinus contorta. Ecol. Evol. 3, 437–449. https://doi.org/10.1002/ece3.426 (2013).Article 
PubMed 
PubMed Central 

Google Scholar 
Tamini, Z. Etude ethnobotanique et analyses morphophysiologiques du développement de la lentille de terre [Macrotyloma geocarpum (harms) Maréchal et Baudet] (Université de Ouagadougou, 1997).
Google Scholar 
Yohannes, H. A review on relationship between climate change and agriculture. J. Earth Sci. Clim. Change 7, 1–8. https://doi.org/10.4172/2157-7617.1000335 (2015).Article 

Google Scholar 
Sileshi, G. et al. Variation in maize yield gaps with plant nutrient inputs, soil type and climate across sub-Saharan Africa. Field Crops Res. 116, 1–13. https://doi.org/10.1016/j.fcr.2009.11.014 (2010).Article 

Google Scholar 
Padi, F. K. & Ehlers, J. D. Effectiveness of early generation selection in cowpea for grain yield and agronomic characteristics in Semiarid West Africa. Crop Sci. 48, 533–540. https://doi.org/10.2135/cropsci2007.05.0265 (2008).Article 

Google Scholar 
Kouelo, K. A. F. et al. Impact du travail du sol et de la fertilisation minérale sur la productivité de [Macrotyloma geocarpum (Harms) Maréchal et Baudet] au centre du Bénin. J. Appl. Biosci. 51, 3625–3632 (2012).
Google Scholar 
Wellenreuther, M., Larson, K. W. & Svensson, E. I. Climatic niche divergence or conservatism? Environmental niches and range limits in ecologically similar damselflie. Ecology 93, 1353–1366. https://doi.org/10.1890/11-1181.1 (2012).Article 
PubMed 

Google Scholar 
Akohoue, F., Achigan-Dako, E. G., Coulibaly, M. & Sibiya, J. Correlations, path coefficient analysis and phenotypic diversity of a West African germplasm of Kersting’s groundnut [Macrotyloma geocarpum (Harms) Maréchal & Baudet]. Genet. Resour. Crop Evol. 66, 1825–1842. https://doi.org/10.1007/s10722-019-00839-w (2019).Article 

Google Scholar 
Assogba, P. et al. Indigenous knowledge and agro-morphological evaluation of the minor crop Kersting’s groundnut [Macrotyloma geocarpum (Harms) Maréchal et Baudet] cultivars of Benin. Genet. Resour. Crop Evol. 63, 513–529. https://doi.org/10.1007/s10722-015-0268-9 (2015).Article 

Google Scholar 
Adu-Gyamfi, R., Dzomeku, I. K. & Lardi, J. Evaluation of growth and yield potential of genotypes of Kersting’s groundnut (Macrotyloma geocarpum Harms) in Northern Ghana. Int. Res. J. Agric. Sci. Soil Sci. 2, 509–515 (2012).
Google Scholar 
Burke, M. B., Lobell, D. B. & Guarino, L. Shifts in African crop climates by 2050, and the implications for crop improvement and genetic resources conservation. Glob. Environ. Change 19, 317–325. https://doi.org/10.1016/j.gloenvcha.2009.04.003 (2009).Article 

Google Scholar 
Ramirez-Cabral, N. Y. Z., Kumar, L. & Shabani, F. Global alterations in areas of suitability for maize production from climate change and using a mechanistic species distribution model (CLIMEX). Sci. Rep. 7, 1–13. https://doi.org/10.1038/s41598-017-05804-0 (2017).CAS 
Article 

Google Scholar 
Hancock, A. M. et al. Adaptation to climate across the Arabidopsis thaliana genome. Science 334, 83–86. https://doi.org/10.1126/science.1209244 (2011).ADS 
CAS 
Article 
PubMed 

Google Scholar 
Bellon, M. R. & van Etten, J. Climate change and on-farm conservation of crop landraces in centres of diversity. In Plant Genetic Resources and Climate Change Vol. 30 (eds Jackson, M. et al.) (CAB International, 2014).
Google Scholar 
Lane, A. & Jarvis, A. Changes in climate will modify the geography of crop suitability: Agricultural biodiversity can help with adaptation. ISI J 4, 12 (2007).
Google Scholar 
Vigouroux, Y., Barnaud, A., Scarcelli, N. & Thuillet, A. C. Biodiversity, evolution and adaptation of cultivated crops. C.R. Biol. 334, 450–457. https://doi.org/10.1016/j.crvi.2011.03.003 (2011).Article 
PubMed 

Google Scholar 
Coulibaly, M., Agossou, C. O. A., Akohoué, F., Sawadogo, M. & Achigan-Dako, E. G. Farmers’ preferences for genetic resources of Kersting’s groundnut [Macrotyloma geocarpum (Harms) Maréchal and Baudet] in the production systems of Burkina Faso and Ghana. Agronomy 10, 371. https://doi.org/10.3390/agronomy10030371 (2020).CAS 
Article 

Google Scholar 
Sohn, N., Fernandez, M. H., Papes, M. & Anciães, M. Ecological Niche modeling in practice: Flagship species and regional conservation planning. Oecol. Aust. 17, 429–440. https://doi.org/10.4257/oeco.2013.1703.11 (2013).Article 

Google Scholar 
Amujoyegbe, B., Obisesan, I., Ajayi, A. & Aderanti, F. Disappearance of Kersting’s groundnut [Macrotyloma geocarpum (harms) Maréchal et Baudet] in South-Western Nigeria: An indicator of genetic erosion. Plant Gen Res News 152, 45–50 (2007).
Google Scholar 
Banta, J. A. et al. Climate envelope modelling reveals intraspecific relationships among flowering phenology, niche breadth and potential range size in Arabidopsis thaliana. Ecol. Lett. 15, 769–777. https://doi.org/10.1111/j.1461-0248.2012.01796.x (2012).Article 
PubMed 

Google Scholar 
Kumar, J., Choudhary, A. K., Gupta, D. S. & Kumar, S. Towards exploitation of adaptive traits for climate-resilient smart pulses. Int. J. Mol. Sci. 20, 1–30. https://doi.org/10.3390/ijms20122971 (2019).ADS 
CAS 
Article 

Google Scholar 
Bohra, A., Mir, R. R., Jha, R., Maurya, A. K. & Varshney, R. K. Advances in genomics and molecular breeding for legume improvement. In Advancement in Crop Improvement Techniques (eds Bohra, A. et al.) 129–139 (Elsevier Inc, 2020).Chapter 

Google Scholar 
Gobu, R. et al. Accelerated crop breeding towards development of climate resilient varieties. In Climate Change and Indian Agriculture: Challenges and Adaptation Strategies (eds Srinivasarao, C. et al.) 49–69 (ICAR-National Academy of Agricultural Research Management, 2020).
Google Scholar 
Aliyu, S., Massawe, F. & Mayes, S. Genetic diversity and population structure of Bambara groundnut [Vigna subterranea (L.) Verdc.]: Synopsis of the past two decades of analysis and implications for crop improvement programmes. Genet. Resour. Crop. Evol. 63, 925–943. https://doi.org/10.1007/s10722-016-0406-z (2016).Article 

Google Scholar 
Al-Khayri, J. M., Jain, S. M., Johnson, D. V. Springer Nature Switzerland AG. Switzerland. https://doi.org/10.1007/978-3-030-23400-3(2019).Kilian, A. et al. Diversity arrays technology: a generic genome profiling technology on open platforms. Methods Mol. Biol. 888, 67–89. https://doi.org/10.1007/978-1-61779-870-2_5 (2012).Article 
PubMed 

Google Scholar 
Illumina, I. HiSeq®2500 Sequencing System: Unsurpassed power and efficiency for production scale sequencing. System Specification Sheet: Sequencing, 1–4. https://www.illumina.com/documents/products/datasheets/datasheet_hiseq2500.pdf (2015).Buckler, E. et al. User Manual for TASSEL Trait Analysis by association, Evolution and Linkage Version 5.0. The Buckler Lab at Cornell University. 1–70. https://www.maizegenetics.net/tassel (2014).Pritchard, J. K., Wen, X., Falush, D. Department of Human Genetics, University of Chicago. (2010).Francis, R. M. pophelper: An R package and web app to analyse and visualize population structure. Mol. Ecol. Resour. 17, 27–32. https://doi.org/10.1111/1755-0998.12509 (2016).CAS 
Article 
PubMed 

Google Scholar 
Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software structure: A simulation study. Mol. Ecol. 14, 2611–2620. https://doi.org/10.1111/j.1365-294X.2005.02553.x (2005).CAS 
Article 
PubMed 

Google Scholar 
Jombart, T. & Ahmed, I. adegenet version 1.3-1: New tools for the analysis of genome-wide SNP data. Bioinformatics https://doi.org/10.1093/bioinformatics/btr521 (2011).Article 
PubMed 
PubMed Central 

Google Scholar 
Boria, R. A., Olson, L. E., Goodman, S. M. & Anderson, R. P. Spatial filtering to reduce sampling bias can improve the performance of ecological niche models. Ecol. Model 275, 73–77. https://doi.org/10.1016/j.ecolmodel.2013.12.012 (2014).Article 

Google Scholar 
Peterson, A. T. et al. NicheBook (Princeton University Press, 2011).
Google Scholar 
Kass, J. M., Pinilla-Buitrago, G. E., Vilela, B., Aiello-Lammens, M. E., Muscarella, R., Merow, C., Anderson, R. P., Wallace: A Modular Platform for Reproducible Modeling of Species Niches and Distributions. R package version 1.0.6.3. (2020).Aiello-Lammens, M. E., Boria, R. A., Radosavljevic, A., Vilela, B. & Anderson, R. P. spThin: An R package for spatial thinning of species occurrence records for use in ecological niche models. Ecography 38, 541–545 (2015).Article 

Google Scholar 
Platts, P. J., Omeny, P. A. & Marchant, R. AFRICLIM: High-resolution climate projections for ecological applications in Africa. Afr. J. Ecol. 53, 103–108 (2014).Article 

Google Scholar 
Meinshausen, M. et al. The RCP greenhouse gas concentrations and their extensions from 1765 to 2300. Clim. Change 109, 213–241. https://doi.org/10.1007/s10584-011-0156-z (2011).ADS 
CAS 
Article 

Google Scholar 
IPCC, I. P. o. C. C. Cambridge University Press. https://doi.org/10.1017/CBO9781107415324(2013).Ramirez-Villegas, J. & Jarvis, A. Downscaling global circulation model outputs: The delta method decision and policy analysis, Working Paper No. 1. Policy Anal. Manag. 1, 1–18 (2010).
Google Scholar 
Hengl, T. et al. Mapping soil properties of Africa at 250 m resolution: Random forests significantly improve current predictions. PLoS ONE 10, 1–26. https://doi.org/10.1371/journal.pone.0125814 (2015).CAS 
Article 

Google Scholar 
Osorio-Olvera, L. et al. ntbox: An R package with graphical user interface for modeling and evaluating multidimensional ecological niches Methods. Ecol. Evol. 11, 1199–1206. https://doi.org/10.1111/2041-210X.13452 (2020).Article 

Google Scholar 
Phillips, S. J., Anderson, R. P. & Schapire, R. E. Maximum Entropy modeling of species geographic distributions. Ecol. Model 190, 231–259. https://doi.org/10.1016/j.ecolmodel.2005.03.026 (2005).Article 

Google Scholar 
Elith, J. et al. A statistical explanation of MaxEnt for ecologists. Diversity Distrib. 17, 43–57. https://doi.org/10.1111/j.1472-4642.2010.00725.x (2010).Article 

Google Scholar 
Richards, C. L., Carstens, B. C. & Lacey Knowles, L. Distribution modelling and statistical phylogeography: An integrative framework for generating and testing alternative biogeographical hypotheses. J. Biogeogr. 34, 1833–1845. https://doi.org/10.1111/j.1365-2699.2007.01814.x (2007).Article 

Google Scholar 
Phillips, S. J. & Dudík, M. Modeling of species distributions with Maxent: New extensions and a comprehensive evaluation. Ecography 31, 161–175. https://doi.org/10.1111/j.2007.0906-7590.05203.x (2008).Article 

Google Scholar 
Merow, C., Smith, M. J. & Silander, J. A. A practical guide to MaxEnt for modeling species’ distributions: What it does, and why inputs and settings matter. Ecography 36, 1058–1069. https://doi.org/10.1111/j.1600-0587.2013.07872.x (2013).Article 

Google Scholar 
Phillips, S. J. et al. Sample selection bias and presence-only distribution models: Implications for background and pseudo-absence data. Ecol. Appl. 19, 181–197 (2009).Article 

Google Scholar 
Barve, N., Barve, V. ENMGadgets: Pre and post processing in ENM Workflow. R package version 0.1.0.1. (2019).Kass, J. M. et al. ENMeval 20: Redesigned for customizable and reproducible modeling of species’ niches and distributions. Ecol. Evolut. https://doi.org/10.1111/2041-210X.13628 (2021).Article 

Google Scholar 
Warren, D. L. & Seifert, S. N. Ecological niche modeling in Maxent: The importance of model complexity and the performance of model selection criteria. Ecol. Appl. 21, 335–342. https://doi.org/10.1890/10-1171.1 (2011).Article 
PubMed 

Google Scholar 
Warren, D. L., Glor, R. E. & Turelli, M. Environmental Niche equivalency versus conservatism: Quantitative approaches to Niche evolution. Evolution 62, 2868–2883. https://doi.org/10.1111/j.1558-5646.2008.00482.x (2008).Article 
PubMed 

Google Scholar 
Broennimann, O. et al. Measuring ecological niche overlap from occurrence and spatial environmental data. Glob. Ecol. Biogeogr. 21, 481–497. https://doi.org/10.1111/j.1466-8238.2011.00698.x (2012).Article 

Google Scholar 
Benhamou, S. & Cornélis, D. Incorporating movement behavior and barriers to improve kernel home range space use estimates. J. Wildl. Manag. 74, 1353–1360. https://doi.org/10.2193/2009-441 (2010).Article 

Google Scholar  More