More stories

  • in

    Hair cortisol concentration reflects the life cycle and management of grey wolves across four European populations

    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

  • in

    Winter torpor expression varies in four bat species with differential susceptibility to white-nose syndrome

    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

  • in

    Strain-specific transcriptional responses overshadow salinity effects in a marine diatom sampled along the Baltic Sea salinity cline

    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

  • in

    Vulnerability to climate change of species in protected areas in Thailand

    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

  • in

    Restructuring of plankton genomic biogeography in the surface ocean under climate change

    Field, C. B., Behrenfeld, M. J., Randerson, J. T. & Falkowski, P. Primary production of the biosphere: integrating terrestrial and oceanic components. Science https://doi.org/10.1126/science.281.5374.237 (1998).Guidi, L. et al. Plankton networks driving carbon export in the oligotrophic ocean. Nature https://doi.org/10.1038/nature16942 (2016).Henson, S. A., Sanders, R. & Madsen, E. Global patterns in efficiency of particulate organic carbon export and transfer to the deep ocean. Glob. Biogeochem. Cycles https://doi.org/10.1029/2011GB004099 (2012).Azam, F. et al. The ecological role of water-column microbes in the sea. Mar. Ecol. Prog. Ser. https://doi.org/10.3354/meps010257 (1983).Saab, M. A. Day-to-day variation in phytoplankton assemblages during spring blooming in a fixed station along the Lebanese coastline. J. Plankton Res. https://doi.org/10.1093/plankt/14.8.1099 (1992).Djurhuus, A. et al. Environmental DNA reveals seasonal shifts and potential interactions in a marine community. Nat. Commun. https://doi.org/10.1038/s41467-019-14105-1 (2020).Kavanaugh, M. T. et al. Seascapes as a new vernacular for pelagic ocean monitoring, management and conservation. ICES J. Mar. Sci. https://doi.org/10.1093/icesjms/fsw086 (2016).Longhurst, A. R. Ecological Geography of the Sea (Elsevier, 2007).Fay, A. R. & McKinley, G. A. Global open-ocean biomes: mean and temporal variability. Earth Syst. Sci. Data https://doi.org/10.5194/essd-6-273-2014 (2014).Reygondeau, G. et al. Dynamic biogeochemical provinces in the global ocean. Glob. Biogeochem. Cycles https://doi.org/10.1002/gbc.20089 (2013).Richter, D. J. et al. Genomic evidence for global ocean plankton biogeography shaped by large-scale current systems. Preprint at bioRxiv https://doi.org/10.1101/867739 (2020).Dutkiewicz, S. et al. Dimensions of marine phytoplankton diversity. Biogeosciences https://doi.org/10.5194/bg-17-609-2020 (2020).Hellweger, F. L., Van Sebille, E. & Fredrick, N. D. Biogeographic patterns in ocean microbes emerge in a neutral agent-based model. Science https://doi.org/10.1126/science.1254421 (2014).Laso-Jadart, R. et al. Investigating population-scale allelic differential expression in wild populations of Oithona similis (Cyclopoida, Claus, 1866). Ecol. Evol. https://doi.org/10.1002/ece3.6588 (2020).Delmont, T. O. et al. Single-amino acid variants reveal evolutionary processes that shape the biogeography of a global SAR11 subclade. eLife https://doi.org/10.7554/eLife.46497 (2019).Carradec, Q. et al. A global ocean atlas of eukaryotic genes. Nat. Commun. https://doi.org/10.1038/s41467-017-02342-1 (2018).Salazar, G. et al. Gene expression changes and community turnover differentially shape the global ocean metatranscriptome. Cell https://doi.org/10.1016/j.cell.2019.10.014 (2019).Alberti, A. et al. Viral to metazoan marine plankton nucleotide sequences from the Tara Oceans expedition. Sci. Data https://doi.org/10.1038/sdata.2017.93 (2017).Pesant, S. et al. Open science resources for the discovery and analysis of Tara Oceans data. Sci. Data https://doi.org/10.1038/sdata.2015.23 (2015).Karsenti, E. et al. A holistic approach to marine eco-systems biology. PLoS Biol. https://doi.org/10.1371/journal.pbio.1001177 (2011).Duarte, C. M. Seafaring in the 21st century: the Malaspina 2010 circumnavigation expedition. Limnol. Oceanogr. Bull. https://doi.org/10.1002/lob.10008 (2015).Barton, A. D., Irwin, A. J., Finkel, Z. V. & Stock, C. A. Anthropogenic climate change drives shift and shuffle in North Atlantic phytoplankton communities. Proc. Natl Acad. Sci. USA https://doi.org/10.1073/pnas.1519080113 (2016).Benedetti, F., Guilhaumon, F., Adloff, F. & Ayata, S. D. Investigating uncertainties in zooplankton composition shifts under climate change scenarios in the Mediterranean Sea. Ecography https://doi.org/10.1111/ecog.02434 (2018).Beaugrand, G. et al. Prediction of unprecedented biological shifts in the global ocean. Nat. Clim. Change 9, 237–243 (2019).Article 

    Google Scholar 
    Pinsky, M. L., Worm, B., Fogarty, M. J., Sarmiento, J. L. & Levin, S. A. Marine taxa track local climate velocities. Science https://doi.org/10.1126/science.1239352 (2013).Bopp, L. et al. Multiple stressors of ocean ecosystems in the 21st century: projections with CMIP5 models. Biogeosciences https://doi.org/10.5194/bg-10-6225-2013 (2013).Thomas, M. K., Kremer, C. T., Klausmeier, C. A. & Litchman, E. A global pattern of thermal adaptation in marine phytoplankton. Science https://doi.org/10.1126/science.1224836 (2012).Ibarbalz, F. M. et al. Global trends in marine plankton diversity across kingdoms of life. Cell https://doi.org/10.1016/j.cell.2019.10.008 (2019).Busseni, G. et al. Large scale patterns of marine diatom richness: drivers and trends in a changing ocean. Glob. Ecol. Biogeogr. https://doi.org/10.1111/geb.13161 (2020).Hutchinson, G. E. Concluding remarks. Cold Spring Harb. Symp. Quant. Biol. 22, 415–427 (1957).Article 

    Google Scholar 
    Delmont, T. O. et al. Functional repertoire convergence of distantly related eukaryotic plankton lineages revealed by genome-resolved metagenomics. Preprint at bioRxiv https://doi.org/10.1101/2020.10.15.341214 (2020).Delmont, T. O. et al. Heterotrophic bacterial diazotrophs are more abundant than their cyanobacterial counterparts in metagenomes covering most of the sunlit ocean. ISME J. https://doi.org/10.1038/s41396-021-01135-1 (2021).Boyer, et al. World Ocean Database 2013, NOAA Atlas NESDIS 72 (National Oceanic and Atmospheric Administration, 2013); https://doi.org/10.7289/V5NZ85MTSunagawa, S. et al. Tara Oceans: towards global ocean ecosystems biology. Nat. Rev. Microbiol. https://doi.org/10.1038/s41579-020-0364-5 (2020).Moon, K. R. et al. Visualizing structure and transitions in high-dimensional biological data. Nat. Biotechnol. https://doi.org/10.1038/s41587-019-0336-3 (2019).van Vuuren, D. P. et al. The representative concentration pathways: an overview. Climatic Change https://doi.org/10.1007/s10584-011-0148-z (2011).Polovina, J. J., Dunne, J. P., Woodworth, P. A. & Howell, E. A. Projected expansion of the subtropical biome and contraction of the temperate and equatorial upwelling biomes in the North Pacific under global warming. ICES J. Mar. Sci. https://doi.org/10.1093/icesjms/fsq198 (2011).Flombaum, P., Wang, W. L., Primeau, F. W. & Martiny, A. C. Global picophytoplankton niche partitioning predicts overall positive response to ocean warming. Nat. Geosci. https://doi.org/10.1038/s41561-019-0524-2 (2020).Richardson, A. J. In hot water: zooplankton and climate change. ICES J. Mar. Sci. 65, 279–295 (2008).Article 

    Google Scholar 
    Wrightson, L. & Tagliabue, A. Quantifying the impact of climate change on marine diazotrophy: insights from Earth system models. Front. Mar. Sci. 7, 635 (2020).Article 

    Google Scholar 
    Zehr, J. P. & Capone, D. G. Changing perspectives in marine nitrogen fixation. Science 368, eaay9514 (2020).CAS 
    Article 

    Google Scholar 
    Luo, Y.-W. et al. Database of diazotrophs in global ocean: abundance, biomass and nitrogen fixation rates. Earth Syst. Sci. Data 4, 47–73 (2012).Article 

    Google Scholar 
    Eppley, R. W. & Peterson, B. J. Particulate organic matter flux and planktonic new production in the deep ocean. Nature 282, 677–680 (1979).Article 

    Google Scholar 
    Laws, E. A., Falkowski, P. G., Smith, W. O., Ducklow, H. & McCarthy, J. J. Temperature effects on export production in the open ocean. Glob. Biogeochem. Cycles 14, 1231–1246 (2000).CAS 
    Article 

    Google Scholar 
    Agrawal, R. & Srikant, R. in Proceedings of the 20th International Conference on Very Large Data Bases (eds Bocca, J. B. et al.) 487–499 (Morgan Kaufmann, 1994).Laufkötter, C. et al. Projected decreases in future marine export production: the role of the carbon flux through the upper ocean ecosystem. Biogeosciences 13, 4023–4047 (2016).Article 

    Google Scholar 
    Iudicone, D. Some may like it hot. Nat. Geosci. https://doi.org/10.1038/s41561-020-0535-z (2020).Gorsky, G. et al. Expanding Tara Oceans protocols for underway, ecosystemic sampling of the ocean–atmosphere interface during Tara Pacific expedition (2016–2018). Front. Mar. Sci. https://doi.org/10.3389/fmars.2019.00750 (2019).Istace, B. et al. de novo assembly and population genomic survey of natural yeast isolates with the Oxford Nanopore MinION sequencer. Gigascience https://doi.org/10.1093/gigascience/giw018 (2017).Grand, M. M. et al. Developing autonomous observing systems for micronutrient trace metals. Front. Mar. Sci. https://doi.org/10.3389/fmars.2019.00035 (2019).Becker, R. A., Wilks, A. R., Brownrigg, R., Minka, T. P. & Deckmyn, A. maps: Draw geographical maps. R version 3.5.0 https://cran.r-project.org/web/packages/maps/index.html (2021).Jaccard, P. Distribution comparée de la flore alpine dans quelques régions des Alpes occidentales et orientales. Bull. Murith. 31, 81–92 (1902).
    Google Scholar 
    Watson, R. A. A database of global marine commercial, small-scale, illegal and unreported fisheries catch 1950–2014. Sci. Data https://doi.org/10.1038/sdata.2017.39 (2017).Maritime Boundaries Geodatabase: Maritime Boundaries and Exclusive Economic Zones (200NM), version 11 (Flanders Marine Institute, 2019); https://doi.org/10.14284/386Aumont, O., Ethé, C., Tagliabue, A., Bopp, L. & Gehlen, M. PISCES-v2: an ocean biogeochemical model for carbon and ecosystem studies. Geosci. Model Dev. https://doi.org/10.5194/gmd-8-2465-2015 (2015).Bibby, T. S. & Moore, C. M. Silicate:nitrate ratios of upwelled waters control the phytoplankton community sustained by mesoscale eddies in sub-tropical North Atlantic and Pacific. Biogeosciences https://doi.org/10.5194/bg-8-657-2011 (2011).Brun, P., Kiørboe, T., Licandro, P. & Payne, M. R. The predictive skill of species distribution models for plankton in a changing climate. Glob. Change Biol. https://doi.org/10.1111/gcb.13274 (2016).Redfield, A. C. in James Johnstone Memorial Volume (ed. Daniel, R. J.) 176–192 (Liverpool Univ. Press, 1934).Michelangeli, P. A., Vrac, M. & Loukos, H. Probabilistic downscaling approaches: application to wind cumulative distribution functions. Geophys. Res. Lett. https://doi.org/10.1029/2009GL038401 (2009).Ridgeway, G. gbm: Generalized boosted regression models. R version 1.6–3.1 https://cran.r-project.org/web/packages/gbm/gbm.pdf (2010).Breiman, L. & Cutler, A. randomForest: Breiman and Cutler’s random forests for classification and regression. R package 4.1.0 https://www.stat.berkeley.edu/~breiman/RandomForests/ (2012).Venables, W. N. & Ripley, B. D. Modern Applied Statistics with S 4th edn (Springer, 2002).Wood, S. N. Stable and efficient multiple smoothing parameter estimation for generalized additive models. J. Am. Stat. Assoc. https://doi.org/10.1198/016214504000000980 (2004).Fawcett, T. An introduction to ROC analysis. Pattern Recognit. Lett. https://doi.org/10.1016/j.patrec.2005.10.010 (2006).Biecek, P. DALEX: explainers for complex predictive models. J. Mach. Learn. Res. 19, 1–5 (2018).
    Google Scholar 
    Jones, M. C. & Cheung, W. W. L. Multi-model ensemble projections of climate change effects on global marine biodiversity. ICES J. Mar. Sci. https://doi.org/10.1093/icesjms/fsu172 (2015).Vallejos, C. A. Exploring a world of a thousand dimensions. Nat. Biotechnol. https://doi.org/10.1038/s41587-019-0330-9 (2019).Kaufman, L. and Rousseeuw, P.J. in Statistical Data Analysis Based on the L1 Norm and Related Methods (ed. Dodge, Y.) 405–416 (North-Holland, 1987).Rousseeuw, P. J. Silhouettes: a graphical aid to the interpretation and validation of cluster analysis. J. Comput. Appl. Math. https://doi.org/10.1016/0377-0427(87)90125-7 (1987).Orsi, A. H., Whitworth, T. & Nowlin, W. D. On the meridional extent and fronts of the Antarctic Circumpolar Current. Deep Sea Res. Part I https://doi.org/10.1016/0967-0637(95)00021-W (1995).Hubert, L. & Arabie, P. Comparing partitions. J. Classif. https://doi.org/10.1007/BF01908075 (1985).Somerfield, P. J. Identification of the Bray–Curtis similarity index: comment on Yoshioka (2008). Mar. Ecol. Prog. Ser. https://doi.org/10.3354/meps07841 (2008).Bloom, S. Similarity indices in community studies: potential pitfalls. Mar. Ecol. Prog. Ser. https://doi.org/10.3354/meps005125 (1981).Welch, B. L. The generalisation of student’s problems when several different population variances are involved. Biometrika 34, 28–35 (1947).CAS 

    Google Scholar 
    Holm, S. A simple sequentially rejective multiple test procedure. Scand. J. Stat. 6, 65–70 (1979).
    Google Scholar 
    Mann, H. B. & Whitney, D. R. On a test of whether one of two random variables is stochastically larger than the other. Ann. Math. Stat. 18, 50–60 (1947).Article 

    Google Scholar 
    Sthle, L. & Wold, S. Analysis of variance (ANOVA). Chemom. Intell. Lab. Syst. 6, 259–272 (1989).Article 

    Google Scholar 
    Bozdogan, H. Model selection and Akaike’s Information Criterion (AIC): the general theory and its analytical extensions. Psychometrika 52, 345–370 (1987).Article 

    Google Scholar 
    Frémont, P. et al. Biogeographies of genomic provinces from ‘Restructuring of plankton genomic biogeography in the surface ocean under climate change’. figshare. https://figshare.com/articles/dataset/Biogeographies_genomic_provinces/19071620 (2022). More

  • in

    Biological trade-offs underpin coral reef ecosystem functioning

    Welti, N. et al. Bridging food webs, ecosystem metabolism, and biogeochemistry using ecological stoichiometry theory. Front. Microbiol. 8, 1298 (2017).Article 

    Google Scholar 
    Ceballos, G. et al. Accelerated modern human-induced species losses: entering the sixth mass extinction. Sci. Adv. 1, e14002 (2015).Hughes, T. P. et al. Global warming and recurrent mass bleaching of corals. Nature 543, 373–377 (2017).CAS 
    Article 

    Google Scholar 
    Pauly, D. et al. Towards sustainability in world fisheries. Nature 418, 689–695 (2002).Bellwood, D. R., Streit, R. P., Brandl, S. J. & Tebbett, S. B. The meaning of the term ‘function’ in ecology: a coral reef perspective. Funct. Ecol. 33, 948–961 (2019).Williams, G. J. et al. Coral reef ecology in the Anthropocene. Funct. Ecol. 33, 1014–1022 (2019).Article 

    Google Scholar 
    Brandl, S. J. et al. Coral reef ecosystem functioning: eight core processes and the role of biodiversity. Front. Ecol. Environ. 17, 445–454 (2019).Article 

    Google Scholar 
    Cinner, J. E. et al. Meeting fisheries, ecosystem function, and biodiversity goals in a human-dominated world. Science 368, 307–311 (2020).CAS 
    Article 

    Google Scholar 
    Mouillot, D. et al. Functional over-redundancy and high functional vulnerability in global fish faunas on tropical reefs. Proc. Natl Acad. Sci. USA 111, 13757–13762 (2014).CAS 
    Article 

    Google Scholar 
    Mora, C. et al. Global human footprint on the linkage between biodiversity and ecosystem functioning in reef fishes. PLoS Biol. 9, e1000606 (2011).CAS 
    Article 

    Google Scholar 
    Barneche, D. R. et al. Scaling metabolism from individuals to reef-fish communities at broad spatial scales. Ecol. Lett. 17, 1067–1076 (2014).CAS 
    Article 

    Google Scholar 
    McIntyre, P. B. et al. Fish distributions and nutrient cycling in streams: can fish create biogeochemical hotspots? Ecology 89, 2335–2346 (2008).Article 

    Google Scholar 
    Allgeier, J. E., Layman, C. A., Mumby, P. J. & Rosemond, A. D. Consistent nutrient storage and supply mediated by diverse fish communities in coral reef ecosystems. Glob. Change Biol. 20, 2459–2472 (2014).Article 

    Google Scholar 
    Morais, R. A. & Bellwood, D. R. Pelagic subsidies underpin fish productivity on a degraded coral reef. Curr. Biol. 29, 1521–1527.e6 (2019).CAS 
    Article 

    Google Scholar 
    Morais, R. A., Connolly, S. R. & Bellwood, D. R. Human exploitation shapes productivity–biomass relationships on coral reefs. Glob. Change Biol. 26, 1295–1305 (2020).Article 

    Google Scholar 
    Barneche, D. R. et al. Body size, reef area and temperature predict global reef-fish species richness across spatial scales. Glob. Ecol. Biogeogr. 28, 315–327 (2019).Article 

    Google Scholar 
    Schiettekatte, N. M. D. et al. Nutrient limitation, bioenergetics and stoichiometry: a new model to predict elemental fluxes mediated by fishes. Funct. Ecol. 34, 1857–1869 (2020).Article 

    Google Scholar 
    Schramski, J. R., Dell, A. I., Grady, J. M., Sibly, R. M. & Brown, J. H. Metabolic theory predicts whole-ecosystem properties. Proc. Natl Acad. Sci. USA 112, 2617–2622 (2015).CAS 
    Article 

    Google Scholar 
    Morais, R. A. & Bellwood, D. R. Global drivers of reef fish growth. Fish Fish. 19, 874–889 (2018).Article 

    Google Scholar 
    Hood, J. M., Vanni, M. J. & Flecker, A. S. Nutrient recycling by two phosphorus-rich grazing catfish: the potential for phosphorus-limitation of fish growth. Oecologia 146, 247–257 (2005).Article 

    Google Scholar 
    Barneche, D. R. & Allen, A. P. The energetics of fish growth and how it constrains food-web trophic structure. Ecol. Lett. 21, 836–844 (2018).Article 

    Google Scholar 
    Brandl, S. J. et al. Demographic dynamics of the smallest marine vertebrates fuel coral reef ecosystem functioning. Science 364, 1189–1192 (2019).CAS 
    Article 

    Google Scholar 
    Lefcheck, J. S. et al. Tropical fish diversity enhances coral reef functioning across multiple scales. Sci. Adv. 5, eaav6420 (2019).Topor, Z. M., Rasher, D. B., Duffy, J. E. & Brandl, S. J. Marine protected areas enhance coral reef functioning by promoting fish biodiversity. Conserv. Lett. 12, e12638 (2019).Article 

    Google Scholar 
    Bellwood, D. R., Hughes, T. P. & Hoey, A. S. Sleeping functional group drives coral-reef recovery. Curr. Biol. 16, 2434–2439 (2006).CAS 
    Article 

    Google Scholar 
    Darling, E. S. & D’agata, S. Coral reefs: fishing for sustainability. Curr. Biol. 27, R65–R68 (2017).CAS 
    Article 

    Google Scholar 
    Graham, N. A. J. et al. Human disruption of coral reef trophic structure. Curr. Biol. 27, 231–236 (2017).CAS 
    Article 

    Google Scholar 
    Graham, N. A. J. et al. Dynamic fragility of oceanic coral reef ecosystems. Proc. Natl Acad. Sci. USA 103, 8425–8429 (2006).CAS 
    Article 

    Google Scholar 
    Stuart-Smith, R. D., Brown, C. J., Ceccarelli, D. M. & Edgar, G. J. Ecosystem restructuring along the great barrier reef following mass coral bleaching. Nature 560, 92–96 (2018).CAS 
    Article 

    Google Scholar 
    Burkepile, D. E. et al. Nutrient supply from fishes facilitates macroalgae and suppresses corals in a Caribbean coral reef ecosystem. Sci. Rep. 3, 1493 (2013).CAS 
    Article 

    Google Scholar 
    Graham, N. A. J. et al. Changing role of coral reef marine reserves in a warming climate. Nat. Commun. 11, 2000 (2020).Reynolds, R. W. et al. Daily high-resolution-blended analyses for sea surface temperature. J. Clim. 20, 5473–5496 (2007).Article 

    Google Scholar 
    Froese, R., Thorson, J. T. & Reyes, R. B. A Bayesian approach for estimating length–weight relationships in fishes. J. Appl. Ichthyol. 30, 78–85 (2014).Article 

    Google Scholar 
    Froese, R. & Pauly, D. FishBase (2018); https://www.fishbase.in/home.htmParravicini, V. et al. Delineating reef fish trophic guilds with global gut content data synthesis and phylogeny. PLoS Biol. 18, e3000702 (2020).CAS 
    Article 

    Google Scholar 
    Brown, J. H., Gillooly, J. F., Allen, A. P., Savage, V. M. & West, G. B. Toward a metabolic theory of ecology. Ecology 85, 1771–1789 (2004).Article 

    Google Scholar 
    Bürkner, P.-C. brms: an R package for Bayesian multilevel models using Stan. J. Stat. Softw. 80, 1–28 (2017).Article 

    Google Scholar 
    Carpenter, B. et al. Stan: a probabilistic programming language. J. Stat. Softw. 76, 1–31 (2017).Article 

    Google Scholar  More

  • in

    Latitudinal gradients in avian colourfulness

    Darwin, C. R. On the Origin of Species, or the Preservation of Favoured Races in the Struggle for Life (John Murray, 1859).Wallace, A. R. Natural Selection and Tropical Nature: Essays on Descriptive and Theoretical Biology 2nd edn (Macmillan, 1895).Darwin, C. R. A Naturalist’s Voyage Round the World (John Murray, 1913).Wallace, A. R. Colour in nature. Nature 19, 580–581 (1879).
    Google Scholar 
    Dalrymple, R. L. et al. Abiotic and biotic predictors of macroecological patterns in bird and butterfly coloration. Ecol. Monogr. 88, 204–224 (2018).
    Google Scholar 
    Adams, J. M., Kang, C. & June-Wells, M. Are tropical butterflies more colorful? Ecol. Res. 29, 685–691 (2014).
    Google Scholar 
    Bailey, S. F. Latitudinal gradients in colors and patterns of passerine birds. Condor 80, 372–381 (1978).
    Google Scholar 
    Wilson, M. F. & Von Neaumann, R. A. Why are neotropical birds more colourful than North American birds? Avicultural Mag. 78, 141–147 (1972).
    Google Scholar 
    Dalrymple, R. L. et al. Birds, butterflies and flowers in the tropics are not more colourful than those at higher latitudes. Glob. Ecol. Biogeogr. 24, 1424–1432 (2015).
    Google Scholar 
    Friedman, N. R. & Remeš, V. Ecogeographical gradients in plumage coloration among Australasian songbird clades. Glob. Ecol. Biogeogr. 26, 261–274 (2017).
    Google Scholar 
    Dale, J., Dey, C. J., Delhey, K., Kempenaers, B. & Valcu, M. The effects of life history and sexual selection on male and female plumage colouration. Nature 527, 367–370 (2015).CAS 

    Google Scholar 
    Dunn, P. O., Armenta, J. K. & Whittingham, L. A. Natural and sexual selection act on different axes of variation in avian plumage color. Sci. Adv. 1, e1400155 (2015).PubMed 
    PubMed Central 

    Google Scholar 
    Stoddard, M. C. & Prum, R. O. How colorful are birds? Evolution of the avian plumage color gamut. Behav. Ecol. 22, 1042–1052 (2011).
    Google Scholar 
    Renoult, J. P., Kelber, A. & Schaefer, H. M. Colour spaces in ecology and evolutionary biology. Biol. Rev. 92, 292–315 (2017).
    Google Scholar 
    Stoddard, M. C. & Prum, R. O. Evolution of avian plumage color in a tetrahedral color space: a phylogenetic analysis of New World buntings. Am. Nat. 171, 755–776 (2008).
    Google Scholar 
    Delhey, K. The colour of an avifauna: a quantitative analysis of the colour of Australian birds. Sci. Rep. 5, 18514 (2015).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Olson, D. M. et al. Terrestrial ecoregions of the world: a new map of life on Earth. Bioscience 51, 933–938 (2001).
    Google Scholar 
    Rabosky, D. L. et al. An inverse latitudinal gradient in speciation rate for marine fishes. Nature 559, 392–395 (2018).CAS 

    Google Scholar 
    Lynch, M. Methods for the analysis of comparative data in evolutionary biology. Evolution 45, 1065–1080 (1991).PubMed 
    PubMed Central 

    Google Scholar 
    Delhey, K. A review of Gloger’s rule, an ecogeographical rule of colour: definitions, interpretations and evidence. Biol. Rev. Camb. Phil. Soc. 94, 1294–1316 (2019).
    Google Scholar 
    Marchetti, K. Dark habitats and bright birds illustrate the role of the environment in species divergence. Nature 362, 149–152 (1993).
    Google Scholar 
    Endler, J. A. The color of light in forests and its implications. Ecol. Monogr. 63, 1–27 (1993).
    Google Scholar 
    Schemske, D. W. in Speciation and Patterns of Diversity Vol. 12 (eds Butlin, R. et al.) 219–239 (Cambridge Univ. Press, 2009).Schemske, D. W., Mittelbach, G. G., Cornell, H. V., Sobel, J. M. & Roy, K. Is there a latitudinal gradient in the importance of biotic interactions? Annu. Rev. Ecol. Evol. Syst. 40, 245–269 (2009).
    Google Scholar 
    MacArthur, R. H. Patterns of communities in the tropics. Biol. J. Linn. Soc. 1, 19–30 (1969).
    Google Scholar 
    Hadfield, J. D. & Nakagawa, S. General quantitative genetic methods for comparative biology: phylogenies, taxonomies and multi-trait models for continuous and categorical characters. J. Evol. Biol. 23, 494–508 (2010).CAS 

    Google Scholar 
    Cooney, C. R. et al. Sexual selection predicts the rate and direction of colour divergence in a large avian radiation. Nat. Commun. 10, 1773 (2019).PubMed 
    PubMed Central 

    Google Scholar 
    Cooney, C. R., MacGregor, H. E. A., Seddon, N. & Tobias, J. A. Multi-modal signal evolution in birds: re-assessing a standard proxy for sexual selection. Proc. R. Soc. B 285, 20181557 (2018).PubMed 
    PubMed Central 

    Google Scholar 
    van der Bijl, W. et al. Butterfly dichromatism primarily evolved via Darwin’s, not Wallace’s, model. Evol. Lett. 4, 545–555 (2020).PubMed 
    PubMed Central 

    Google Scholar 
    Darwin, C. R. The Descent of Man, and Selection in Relation to Sex (John Murray, 1871).Tobias, J. A., Montgomerie, R. & Lyon, B. E. The evolution of female ornaments and weaponry: social selection, sexual selection and ecological competition. Phil. Trans. R. Soc. B 367, 2274–2293 (2012).PubMed 
    PubMed Central 

    Google Scholar 
    Galván, I., Negro, J. J., Rodríguez, A. & Carrascal, L. M. On showy dwarfs and sober giants: body size as a constraint for the evolution of bird plumage colouration. Acta Ornithol. 48, 65–80 (2013).
    Google Scholar 
    Kiltie, R. A. Scaling of visual acuity with body size in mammals and birds. Funct. Ecol. 14, 226–234 (2000).
    Google Scholar 
    Zahavi, A. & Zahavi, A. The Handicap Principle (Oxford Univ. Press, 1997).Badyaev, A. V. & Hill, G. E. Avian sexual dichromatism in relation to phylogeny and ecology. Annu. Rev. Ecol. Evol. Syst. 34, 27–49 (2003).
    Google Scholar 
    Simpson, R. K., Johnson, M. A. & Murphy, T. G. Migration and the evolution of sexual dichromatism: evolutionary loss of female coloration with migration among wood-warblers. Proc. R. Soc. B 282, 20150375 (2015).PubMed 
    PubMed Central 

    Google Scholar 
    Helferich, G. Humboldt’s Cosmos (Tantor eBooks, 2011).Jetz, W., Thomas, G. H., Joy, J. B., Hartmann, K. & Mooers, A. O. The global diversity of birds in space and time. Nature 491, 444–448 (2012).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    He, Y. et al. Segmenting biological specimens from photos to understand the evolution of UV plumage in passerine birds. Preprint at bioRxiv https://doi.org/10.1101/2021.07.22.453339 (2021).Chen, L. C., Zhu, Y., Papandreou, G., Schroff, F. & Adam, H. Encoder–decoder with atrous separable convolution for semantic image segmentation. Preprint at arXiv https://doi.org/10.48550/arXiv.1802.02611 (2018).Hussein, B. R., Malik, O. A., Ong, W.-H. & Slik, J. W. F. in Computational Science and Technology Lecture Notes in Electrical Engineering (eds Alfred, R. et al.) 321–330 (Springer Singapore, 2020).Troscianko, J. & Stevens, M. Image calibration and analysis toolbox—a free software suite for objectively measuring reflectance, colour and pattern. Methods Ecol. Evol. 6, 1320–1331 (2015).PubMed 
    PubMed Central 

    Google Scholar 
    Hijmans, R. J. raster: Geographic Data Analysis and Modeling. R package version 3.5-15 https://CRAN.R-project.org/package=raster (2022).Maia, R., Gruson, H., Endler, J. A., White, T. E. & O’Hara, R. B. pavo 2: new tools for the spectral and spatial analysis of colour in R. Methods Ecol. Evol. 10, 1097–1107 (2019).
    Google Scholar 
    Stoddard, M. C. et al. Wild hummingbirds discriminate nonspectral colors. Proc. Natl Acad. Sci. USA 117, 15112–15122 (2020).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Gomez, D. & Théry, M. Simultaneous crypsis and conspicuousness in color patterns: comparative analysis of a neotropical rainforest bird community. Am. Nat. 169, S42–S61 (2007).
    Google Scholar 
    Blonder, B. Do hypervolumes have holes? Am. Nat. 187, E93–E105 (2016).
    Google Scholar 
    Schliep, K. P. phangorn: phylogenetic analysis in R. Bioinformatics 27, 592–593 (2011).CAS 

    Google Scholar 
    Fick, S. E. & Hijmans, R. J. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315 (2017).
    Google Scholar 
    Beckmann, M. et al. glUV: a global UV-B radiation data set for macroecological studies. Methods Ecol. Evol. 5, 372–383 (2014).
    Google Scholar 
    Running, S. W. et al. A continuous satellite-derived measure of global terrestrial primary production. Bioscience 54, 547–560 (2004).
    Google Scholar 
    Tobias, J. A. & Pigot, A. L. Integrating behaviour and ecology into global biodiversity conservation strategies. Phil. Trans. R. Soc. B 374, 20190012 (2019).PubMed 
    PubMed Central 

    Google Scholar 
    Dunn, P. O., Whittingham, L. A. & Pitcher, T. E. Mating systems, sperm competition, and the evolution of sexual dimorphism in birds. Evolution 55, 161–175 (2001).CAS 

    Google Scholar 
    Bivand, R. S. & Wong, D. W. S. Comparing implementations of global and local indicators of spatial association. TEST 27, 716–748 (2018).
    Google Scholar 
    Hawkins, B. A. et al. Structural bias in aggregated species-level variables driven by repeated species co-occurrences: a pervasive problem in community and assemblage data. J. Biogeogr. 44, 1199–1211 (2017).
    Google Scholar 
    Hadfield, J. D. MCMC methods for multi-response generalised linear mixed models: the MCMCglmm R package. J. Stat. Softw. 33, 1–22 (2010).
    Google Scholar 
    Healy, K. et al. Ecology and mode-of-life explain lifespan variation in birds and mammals. Proc. R. Soc. B 281, 20140298 (2014).PubMed 
    PubMed Central 

    Google Scholar 
    R Core Team R: A Language and Environment for Statistical Computing (R Foundation for Statistical Computing, 2021); https://www.R-project.org/ More

  • in

    Development and validation of an eDNA protocol for monitoring endemic Asian spiny frogs in the Himalayan region of Pakistan

    Lindenmayer, D. et al. A checklist of attributes for effective monitoring of threatened species and threatened ecosystems. J. Environ. Manage. 262, 110312 (2020).PubMed 

    Google Scholar 
    Reid, A. J. et al. Emerging threats and persistent conservation challenges for freshwater biodiversity. Biol. Rev. 94, 849–873 (2019).PubMed 

    Google Scholar 
    IUCN. The IUCN Red List of Threatened Species. Version 2019-3. http://www.iucnredlist.org (2021).Adams, M. J. et al. Trends in amphibian occupancy in the United States. PLoS ONE 8, e64347 (2013).ADS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Corn, P. S. Climate change and amphibians. Anim. Biodivers. Conserv. 28, 59–67 (2005).
    Google Scholar 
    Kiesecker, J. M., Blaustein, A. R. & Belden, L. K. Complex causes of amphibian population declines. Nature 410, 681–684 (2001).ADS 
    CAS 

    Google Scholar 
    Baldwin, R. F. & deMaynadier, P. G. Assessing threats to pool-breeding amphibian habitat in an urbanizing landscape. Biol. Conserv. 142, 1628–1638 (2009).Borzée, A., Kyong, C. N., Kil, H. K. & Jang, Y. Impact of water quality on the occurrence of two endangered Korean anurans: Dryophytes suweonensis and Pelophylax chosenicus. Herpetologica 74, 1–7 (2018).
    Google Scholar 
    Stuart, S. N. et al. Status and trends of amphibian declines and extinctions worldwide. Science 306, 1783–1786 (2004).ADS 
    CAS 

    Google Scholar 
    Caro, T., Rowe, Z., Berger, J., Wholey, P. & Dobson, A. An inconvenient misconception: climate change is not the principal driver of biodiversity loss. Conserv. Lett. e12868 (2022).Daszak, P. et al. Emerging infectious diseases and amphibian population declines. Emerg. Infect. 5, 735–748 (1999).CAS 

    Google Scholar 
    Fellers, G., Green, D. E. & Longcore, J. Oral chytridiomycosis in the mountain yellow-legged frog (Rana muscosa). Copeia 2001, 945–953Blaustein, A. R. et al. Effects of ultraviolet radiation on amphibians: field experiments. Am. Zool. 38, 799–812 (1999).
    Google Scholar 
    Langhelle, A., Lindell, M. J. & Nyström, P. Effects of ultraviolet radiation on amphibian embryonic and larval development. J. Herpetol. 33, 449–456 (1999).
    Google Scholar 
    Beebee, T. J. C. Amphibians breeding and climate. Nature 374, 219–220 (1995).ADS 
    CAS 

    Google Scholar 
    Donnelly, M. A. & Crump, M. L. Potential effects of climate change on two neotropical amphibian assemblages. In Potential Impacts of Climate Change on Tropical Forest Ecosystems (ed. Markham, A.) 401–421 (Springer Netherlands, 1998).Carey, C. & Alexander, M. A. Climate change and amphibian declines: is there a link? Divers. Distrib. 9, 111–121 (2003).
    Google Scholar 
    Fisher, R. N. & Shaffer, H. B. The decline of amphibians in California’s Great Central Valley. Conserv. Biol. 10, 1387–1397 (1996).
    Google Scholar 
    Sparling, D. W., Donald, W., Linder, G. & Bishop, C. A. Ecotoxicology of Amphibians and Reptiles. (SETAC Press, 2000).Rouse, M. J. & Daellenbach, U. S. Rethinking research methods for the resource-based perspective: isolating sources of sustainable competitive advantage. Strat. Manag. J. 20, 487–494 (1999).
    Google Scholar 
    Bridges, C. M. & Boone, M. D. The interactive effects of UV-B and insecticide exposure on tadpole survival, growth and development. Biol. Conserv. 113, 49–54 (2003).
    Google Scholar 
    Schmeller, D. S. et al. National responsibilities in European species conservation: a methodological review. Conserv. Biol. 22, 593–601 (2008).PubMed 

    Google Scholar 
    Anderson, S. Area and endemism. Q. Rev. Biol. 69, 451–471 (1994).
    Google Scholar 
    Strayer, D. L. & Dudgeon, D. Freshwater biodiversity conservation: recent progress and future challenges. J. N. Am. Benthol. Soc. 29, 344–358 (2010).
    Google Scholar 
    Gorman, C. E., Potts, B. M., Schweitzer, J. A. & Bailey, J. K. Shifts in species interactions due to the evolution of functional differences between endemics and non-endemics: an endemic syndrome hypothesis. PLoS ONE 9, e111190 (2014).Mace, G. M. et al. Quantification of extinction risk: IUCN’s system for classifying threatened species. Conserv. Biol. 22, 1424–1442 (2008).PubMed 

    Google Scholar 
    Fontaine, B. et al. The European Union’s 2010 target: putting rare species in focus. Biol. Conserv. 139, 167–185 (2007).
    Google Scholar 
    Saeed, M. et al. Rise in temperature causes decreased fitness and higher extinction risks in endemic frogs at high altitude forested wetlands in northern Pakistan. J. Therm. Biol. 95, 102809 (2021).McDonald, L. L. Sampling rare populations. In Sampling Rare or Elusive Species: Concepts, Designs, and Techniques for Estimating Population Parameters (ed. Thompson W. L.) 11–42 (Island Press, 2004).Dodd Jr. K. Monitoring Amphibians in Great Smoky Mountains National Park (USGS Survey Circular, 2003).Qu, C. & Stewart, K. A. Evaluating monitoring options for conservation : traditional and environmental DNA tools for a critically endangered mammal. Sci. Nat. 106, 9 (2019).
    Google Scholar 
    Deiner, K. et al. Environmental DNA metabarcoding: transforming how we survey animal and plant communities. Mol. Ecol. 26, 5872–5895 (2017).PubMed 

    Google Scholar 
    Schmidt, B. R., Kery, M., Ursenbacher, S., Hyman, O. J. & Collins, J. P. Site occupancy models in the analysis of environmental DNA presence/absence surveys: a case study of an emerging amphibian pathogen. Methods Ecol. Evol. 4, 646–653 (2013).
    Google Scholar 
    Iknayan, K. J., Tingley, M. W., Furnas, B. J. & Beissinger, S. R. Detecting diversity: emerging methods to estimate species diversity. Trends Ecol. Evol. 29, 97–106 (2014).PubMed 

    Google Scholar 
    Kéry, M. & Schmidt, B. R. Imperfect detection and its consequences for monitoring for conservation. Community Ecol. 9, 207–216 (2008).
    Google Scholar 
    Mackenzie, D. I. et al. Estimating site occupancy rates when detection probabilities are less than one. Ecology 83, 2248–2255 (2002).
    Google Scholar 
    Mackenzie, D. I. & Royle, J. A. Designing occupancy studies: general advice and allocating survey effort. J. Appl. Ecol. 42, 1105–1114 (2005).
    Google Scholar 
    Ficetola, G. F., Miaud, C., Pompanon, F. & Taberlet, P. Species detection using environmental DNA from water samples. Biol. Lett. 4, 423–425 (2008).PubMed 
    PubMed Central 

    Google Scholar 
    Darling, J. A. & Mahon, A. R. From molecules to management: adopting DNA-based methods for monitoring biological invasions in aquatic environments. Environ. Res. 111, 978–988 (2011).CAS 
    PubMed 

    Google Scholar 
    Goldberg, C. S., Pilliod, D. S., Arkle, R. S. & Waits, L. P. Molecular detection of vertebrates in stream water: a demonstration using Rocky Mountain tailed frogs and Idaho giant salamanders. PLoS ONE 6, e22746 (2011).Williams, M. R. et al. Isothermal amplification of environmental DNA (eDNA) for direct field-based monitoring and laboratory confirmation of Dreissena sp. PLoS ONE 12, e0186462 (2017).Agersnap, S. et al. Monitoring of noble, signal and narrow-clawed crayfish using environmental DNA from freshwater samples. PLoS ONE 12, e0179261 (2017).Barnes, M. A. & Turner, C. R. The ecology of environmental DNA and implications for conservation genetics. Conserv. Genet. 17, 1–17 (2016).CAS 

    Google Scholar 
    Bohmann, K. et al. Environmental DNA for wildlife biology and biodiversity monitoring. Trends Ecol. Evol. 29, 358–367 (2014).PubMed 

    Google Scholar 
    Sigsgaard, E. E., Carl, H., Møller, P. R. & Thomsen, P. F. Monitoring the near-extinct European weather loach in Denmark based on environmental DNA from water samples. Biol. Conserv. 183, 46–52 (2015).
    Google Scholar 
    Bedwell, M. E., Hopkins, K. V. S., Dillingham, C. & Goldberg, C. S. Evaluating Sierra Nevada yellow-legged frog distribution using environmental DNA. J. Wildl. Mangaement 85, 945–952 (2021).
    Google Scholar 
    Eiler, A., Löfgren, A., Hjerne, O., Nordén, S. & Saetre, P. Environmental DNA (eDNA) detects the pool frog (Pelophylax lessonae) at times when traditional monitoring methods are insensitive. Sci. Rep. 8, 5452 (2018).ADS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Brozio, S. et al. Development and application of an eDNA method to detect the critically endangered Trinidad golden tree frog (Phytotriades auratus) in bromeliad phytotelmata. PLoS ONE 12, e0170619 (2017).Pellet, J. & Schmidt, B. R. Monitoring distributions using call surveys: estimating site occupancy, detection probabilities and inferring absence. Biol. Conserv. 123, 27–35 (2005).
    Google Scholar 
    Weir, L. A., Royle, J. A., Nanjappa, P. & Jung, R. E. Modeling anuran detection and site occupancy on North American Amphibian Monitoring Program (NAAMP) routes in Maryland. J. Herpetol. 39, 627–639 (2005).
    Google Scholar 
    Fiske, I. J. & Chandler, R. B. Unmarked: an R package for fitting hierarchical models of wildlife occurrence and abundance. J. Stat. Softw. 43, 1–23 (2011).
    Google Scholar 
    Goldberg, C. S. et al. Critical considerations for the application of environmental DNA methods to detect aquatic species. Methods Ecol. Evol. 7, 1299–1307 (2016).
    Google Scholar 
    Holland, M. M. & Parsons, T. J. Mitochondrial DNA sequence analysis – validation and use for forensic casework. Forensic Sci. Rev. 11, 21–50 (1999).CAS 
    PubMed 

    Google Scholar 
    Willerslev, E. et al. Diverse plant and animal genetic records from Holocene and Pleistocene sediments. Science 300, 791–795 (2003).ADS 
    CAS 
    PubMed 

    Google Scholar 
    Waits, L. P. & Paetkau, D. Noninvasive genetic sampling tools for wildlife biologists: a review of applications and recommendations for accurate data collection. J. Wildl. Manage. 69, 1419–1433 (2006).
    Google Scholar 
    Shokralla, S. et al. Next-generation DNA barcoding: using next-generation sequencing to enhance and accelerate DNA barcode capture from single specimens. Mol. Ecol. Resour. 14, 892–901 (2014).CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Mills, L. S., Pilgrim, K. L., Schwartz, M. K. & McKelvey, K. Identifying lynx and other North American felids based on mtDNA analysis. Conserv. Genet. 1, 285–288 (2000).CAS 

    Google Scholar 
    Hajibabaei, M. et al. A minimalist barcode can identify a specimen whose DNA is degraded. Mol. Ecol. Notes 6, 959–964 (2006).CAS 

    Google Scholar 
    Kim, P., Kim, D., Yoon, T. J. & Shin, S. Early detection of marine invasive species, Bugula neritina (Bryozoa: Cheilostomatida), using species-specific primers and environmental DNA analysis in Korea. Mar. Environ. Res. 139, 1–10 (2018).CAS 
    PubMed 

    Google Scholar 
    Dejean, T. et al. Persistence of environmental DNA in freshwater ecosystems. PLoS ONE 6, e23398 (2011).Xia, Z. et al. Early detection of a highly invasive bivalve based on environmental DNA (eDNA). Biol. Invasions 20, 437–447 (2018).
    Google Scholar 
    Torresdal, J. D., Farrell, A. D. & Goldberg, C. S. Environmental DNA detection of the golden tree frog (Phytotriades auratus) in bromeliads. PLoS ONE 12, e0168787 (2017).Biggs, J. et al. Using eDNA to develop a national citizen science-based monitoring programme for the great crested newt (Triturus cristatus). Biol. Conserv. 183, 19–28 (2015).
    Google Scholar 
    Pilliod, D. S., Goldberg, C. S., Arkle, R. S. & Waits, L. P. Estimating occupancy and abundance of stream amphibians using environmental DNA from filtered water samples. Can. J. Fish. Aquat. Sci. 70, 1123–1130 (2013).CAS 

    Google Scholar 
    Smith, D. H. V., Jones, B., Randall, L. & Prescott, D. R. C. Difference in detection and occupancy between two anurans: the importance of species-specific monitoring. Herpetol. Conserv. Biol. 9, 267–277 (2014).
    Google Scholar 
    Bayley, P. B. & Peterson, J. T. An approach to estimate probability of presence and richness of fish species. Trans. Am. Fish. Soc. 130, 620–633 (2004).
    Google Scholar 
    Mehta, S. V., Haight, R. G., Homans, F. R., Polasky, S. & Venette, R. C. Optimal detection and control strategies for invasive species management. Ecol. Econ. 61, 237–245 (2007).
    Google Scholar 
    Scott, Jr., N. J. & Woodward, B. D. Surveys at breeding sites. In Measuring and Monitoring Biological Diversity: Standard Methods for Amphibians (eds. Heyer, W. R., Donnelly, M. A., McDiarmid, R. W., Hayek, L. C., & Foster, M. S.) 118–125 (Smithsonian Institution Press, 1994).Dejean, T. et al. Improved detection of an alien invasive species through environmental DNA barcoding: the example of the American bullfrog Lithobates catesbeianus. J. Appl. Ecol. 49, 953–959 (2012).
    Google Scholar 
    Goldberg, C. S., Sepulveda, A., Ray, A., Baumgardt, J. & Waits, L. P. Environmental DNA as a new method for early detection of New Zealand mudsnails (Potamopyrgus antipodarum). Freshw. Sci. 32, 792–800 (2013).
    Google Scholar 
    Mahon, A. R. et al. Validation of eDNA surveillance sensitivity for detection of Asian carps in controlled and field experiments. PLoS ONE 8, e58316 (2013).ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 
    Khan, M. S. Amphibians and Reptiles of Pakistan (Krieger Publishing Company, 2006).Ruppert, K. M., Davis, D. R., Rahman, M. S. & Kline, R. J. Development and assessment of an environmental DNA (eDNA) assay for a cryptic Siren (Amphibia: Sirenidae). Environ. Adv. 7, 100163 (2022).
    Google Scholar 
    Hobbs, J., Round, J. M., Allison, M. J. & Helbing, C. C. Expansion of the known distribution of the coastal tailed frog, Ascaphus truei, in British Columbia, Canada, using robust eDNA detection methods. PLoS ONE 14, e0213849 (2019).Barata, I. M., Griffiths, R. A., Fogell, D. J. & Buxton, A. S. Comparison of eDNA and visual surveys for rare and cryptic bromeliad-dwelling frogs. Herpetol. J. 31, 1–9 (2021).
    Google Scholar 
    Ahmed, W. et al. Site occupancy of two endemic stream frogs in different forest types in Pakistan. Herpetol. Conserv. Biol. 15, 506–511 (2020).
    Google Scholar 
    Richmond, O. M. W., Hines, J. E. & Beissinger, S. R. Two-species occupancy models: a new parameterization applied to co-occurrence of secretive rails. Ecol. Appl. 20, 2036–2046 (2010).PubMed 

    Google Scholar 
    Shea, C. P., Eaton, M. J. & MacKenzie, D. I. Implementation of an occupancy-based monitoring protocol for a widespread and cryptic species, the New England cottontail (Sylvilagus transitionalis). Wildl. Res. 46, 222–235 (2019).
    Google Scholar 
    Rota, C. T. et al. A multispecies occupancy model for two or more interacting species. Methods Ecol. Evol. 7, 1164–1173 (2016).
    Google Scholar 
    Ohler, A. & Dubois, A. Phylogenetic relationships and generic taxonomy of the tribe Paini (Amphibia, Anura, Ranidae, Dicroglossinae). Zoosystema 28, 769–784 (2006).
    Google Scholar 
    Jiang, J. et al. Phylogenetic relationships of the tribe Paini (Amphibia, Anura, Ranidae) based on partial sequences of mitochondrial 12s and 16s rRNA genes. Zool. Res. 362, 353–362 (2005).
    Google Scholar 
    Rais, M. et al. A note on recapture of Nanorana vicina (Anura: Amphibia) from Murree, Pakistan. J. Anim. Plant Sci. 24, 455–458 (2014).
    Google Scholar 
    Siddiqui, M. F., Ahmed, M., Khan, N. & Khan, I. A. A quantitative description of moist temperate conifer forests of Himalayan region of Pakistan and Azad Kashmir. Int. J. Biotechnol. 7, 175–185 (2010).
    Google Scholar 
    Beck, H. E. et al. Present and future köppen-geiger climate classification maps at 1-km resolution. Sci. Data 5, 180214 (2018).Sheikh, M. I. & Hafeez, S. M. Forest and Forestry in Pakistan (A-one Publishers, 2001).Lodhi, A. Conservation of leopards in Ayubia National Park, Pakistan (MS Thesis) (University of Montana, 2007).Palumbi, S. R. Nucleic acids II: the polymerase chain reaction. In Molecular Systematics, 2nd Edition (eds. Hillis, D. M. et al.) 205–247 (Sinauer, 1996).Vences, M., Thomas, M., Van Der Meijden, A., Chiari, Y. & Vieites, D. R. Comparative performance of the 16S rRNA gene in DNA barcoding of amphibians. Front. Zool. 2, 5 (2005).Pounds, J. A. & Crump, M. L. Amphibian declines and climate disturbance: the case of the golden toad and the harlequin frog. Conserv. Biol. 8, 72–85 (1994).
    Google Scholar 
    R Core Team. R: a language and environment for statistical computing. R foundation for statistical computing, Vienna, Austria. https://www.R-project.org/ (2021).Hutchinson, R. A., Valente, J. J., Emerson, S. C., Betts, M. G. & Dietterich, T. G. Penalized likelihood methods improve parameter estimates in occupancy models. Methods Ecol. Evol. 6, 949–959 (2015).
    Google Scholar 
    Clipp, H. L., Evans, A. L., Kessinger, B. E., Kellner, K., & Rota, C. T. A penalized likelihood for multispecies occupancy models improves predictions of species interactions. Ecology 102, e03520 (2021).PubMed 

    Google Scholar  More