Philippa Kaur

Simpson, G. G. Tempo and Mode in Evolution (Columbia Univ. Press, 1944).Losos, J. B. Adaptive radiation, ecological opportunity, and evolutionary determinism. Am. Nat. 175, 623–639 (2010).PubMed 
PubMed Central 

Google Scholar 
Erwin, D. H. Novelty and innovation in the history of life. Curr. Biol. 25, 930–940 (2015).
Google Scholar 
Novack-Gottshall, P. M. General models of ecological diversification. I. Conceptual synthesis. Paleobiology 42, 185–208 (2016).
Google Scholar 
Marshall, C. R. Explaining the Cambrian “explosion” of animals. Annu. Rev. Earth Planet. Sci. 34, 355–384 (2006).CAS 

Google Scholar 
Hughes, M., Gerber, S. & Wills, M. A. Clades reach highest morphological disparity early in their evolution. Proc. Natl Acad. Sci. USA 110, 13875–13879 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
Hopkins, M. J. & Smith, A. B. Dynamic evolutionary change in post-Paleozoic echinoids and the importance of scale when interpreting changes in rates of evolution. Proc. Natl Acad. Sci. USA 112, 3758–3763 (2015).CAS 
PubMed 
PubMed Central 

Google Scholar 
Wright, D. F. Phenotypic innovation and adaptive constraints in the evolutionary radiation of Palaeozoic crinoids. Sci. Rep. 7, 13745 (2017).PubMed 
PubMed Central 

Google Scholar 
Cantalapiedra, J. L., Prado, J. L., Hernández Fernández, M. & Alberdi, M. T. Decoupled ecomorphological evolution and diversification in Neogene-Quaternary horses. Science 355, 627–630 (2017).CAS 
PubMed 
PubMed Central 

Google Scholar 
Crouch, N. M. A. & Ricklefs, R. E. Speciation rate is independent of the rate of evolution of morphological size, shape, and absolute morphological specialization in a large clade of birds. Am. Naturalist 193, E78–E91 (2019).
Google Scholar 
Slater, G. J. & Friscia, A. R. Hierarchy in adaptive radiation: a case study using the Carnivora (Mammalia). Evolution 73, 524–539 (2019).PubMed 
PubMed Central 

Google Scholar 
Ronco, F. et al. Drivers and dynamics of a massive adaptive radiation in cichlid fishes. Nature https://doi.org/10.1038/s41586-020-2930-4 (2020).Cole, S. R. & Hopkins, M. J. Selectivity and the effect of mass extinctions on disparity and functional ecology. Sci. Adv. https://doi.org/10.1126/sciadv.abf4072 (2021).Erwin, D. H. et al. The Cambrian conundrum: early divergence and later ecological success in the early history of animals. Science 334, 1091–1097 (2011).CAS 
PubMed 
PubMed Central 

Google Scholar 
Smith, A. B., Zamora, S. & Álvaro, J. J. The oldest echinoderm faunas from Gondwana show that echinoderm body plan diversification was rapid. Nat. Commun. 4, 1385 (2013).PubMed 
PubMed Central 

Google Scholar 
Deline, B. et al. Evolution and development at the origin of a phylum. Curr. Biol. 30, 1672–1679 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Pisani, D., Feuda, R., Peterson, K. J. & Smith, A. B. Resolving phylogenetic signal from noise when divergence is rapid: a new look at the old problem of echinoderm class relationships. Mol. Phylogenet. Evol. 62, 27–34 (2012).PubMed 
PubMed Central 

Google Scholar 
Smith, A. B. & Zamora, S. Cambrian spiral-plated echinoderms from Gondwana reveal the earliest pentaradial body plan. Proc. R. Soc. B 280, 20131197 (2013).PubMed 
PubMed Central 

Google Scholar 
Zamora, S. et al. in Early Palaeozoic Biogeography and Palaeogeography. Memoirs 38 (eds Harper, D. A. T. & Servais, T.) 157–171 (Geological Society, London, 2013).Lefebvre, B. et al. in Early Palaeozoic Biogeography and Palaeogeography. Memoirs 38 (eds Harper, D. A. T. & Servais, T.) 173–198 (Geological Society, London, 2013).Novack-Gottshall, P. M. Using a theoretical ecospace to quantify the ecological diversity of Paleozoic and modern marine biotas. Paleobiology 33, 273–294 (2007).
Google Scholar 
Wagner, P. J. On the probabilities of branch durations and stratigraphic gaps in phylogenies of fossil taxa when rates of diversification and sampling vary over time. Paleobiology 45, 30–55 (2019).
Google Scholar 
Sumrall, C. D. & Wray, G. A. Ontogeny in the fossil record: diversification of body plans and the evolution of “aberrant” symmetry in Paleozoic echinoderms. Paleobiology 33, 149–163 (2007).
Google Scholar 
Guensburg, T. E. & Sprinkle, J. in The Ecology of the Cambrian Radiation (eds Zhuravlev, A. I. U. & Riding, R.) 428–444 (Columbia Univ. Press, 2001).Zamora, S., Deline, B., Álvaro, J. J. & Rahman, I. A. The Cambrian Substrate Revolution and the early evolution of attachment in suspension-feeding echinoderms. Earth Sci. Rev. 171, 478–491 (2017).
Google Scholar 
Lloyd, G. T., Guillerme, T., Sherratt, E. & Wang, S. C. Claddis: measuring morphological diversity and evolutionary tempo. R package version 0.4.0 https://cran.r-project.org/package=Claddis (2020).Lloyd, G. T., Wang, S. C. & Brusatte, S. L. Identifying heterogeneity in rates of morphological evolution: discrete character change in the evolution of lungfish (Sarcopterygii; Dipnoi). Evolution 66, 330–348 (2012).PubMed 
PubMed Central 

Google Scholar 
Wagner, P. J. Early bursts of disparity and the reorganization of character integration. Proc. R. Soc. B 285, 20181604 (2018).PubMed 
PubMed Central 

Google Scholar 
Lewis, P. O. A likelihood approach to estimating phylogeny from discrete morphological character data. Syst. Biol. 50, 913–925 (2001).CAS 
PubMed 
PubMed Central 

Google Scholar 
Lloyd, G. T. Estimating morphological diversity and tempo with discrete character–taxon matrices: implementation, challenges, progress, and future directions. Biol. J. Linn. Soc. 118, 131–151 (2016).
Google Scholar 
Harmon, L. J. et al. Early bursts of body size and shape evolution are rare in comparative data. Evolution 64, 2385–2396 (2010).
Google Scholar 
Knope, M. L., Heim, N. A., Frishkoff, L. O. & Payne, J. L. Limited role of functional differentiation in early diversification of animals. Nat. Commun. 6, 6455–6461 (2015).CAS 

Google Scholar 
Cole, S. R., Wright, D. F. & Ausich, W. I. Phylogenetic community paleoecology of one of the earliest complex crinoid faunas (Brechin Lagerstätte, Ordovician). Palaeogeogr. Palaeoclimatol. Palaeoecol. 521, 82–98 (2019).
Google Scholar 
Revell, L. J. phytools: phylogenetic tools for comparative biology (and other things). R package version 0.7-47 https://cran.r-project.org/package=phytools (2020).Stayton, C. T. The definition, recognition, and interpretation of convergent evolution, and two new measures for quantifying and assessing the significance of convergence. Evolution 69, 2140–2153 (2015).PubMed 
PubMed Central 

Google Scholar 
Erwin, D. H. A conceptual framework of evolutionary novelty and innovation. Biol. Rev. 96, 1–15 (2021).
Google Scholar 
Lloyd, G. T. Journeys through discrete-character morphospace: synthesizing phylogeny, tempo, and disparity. Palaeontology 61, 637–645 (2018).
Google Scholar 
Anderson, P. S. L. & Friedman, M. Using cladistic characters to predict functional variety: experiments using early gnathostomes. J. Vertebr. Paleontol. 32, 1254–1270 (2012).
Google Scholar 
Deline, B. et al. Evolution of metazoan morphological disparity. Proc. Natl Acad. Sci. USA 15, E8909–E8918 (2018).
Google Scholar 
Matzke, N. J. & Irmis, R. B. Including autapomorphies is important for paleontological tip-dating with clocklike data, but not with non-clock data. PeerJ 6, e4553 (2018).PubMed 
PubMed Central 

Google Scholar 
Laing, A. M. et al. Giant taxon–character matrices: the future of morphological systematics. Cladistics 34, 333–335 (2018).PubMed 
PubMed Central 

Google Scholar 
Bush, A. M. & Bambach, R. K. Paleoecologic megatrends in marine Metazoa. Annu. Rev. Earth Planet. Sci. 39, 241–269 (2011).CAS 

Google Scholar 
Knope, M. L., Bush, A. M., Frishkoff, L. O., Heim, N. A. & Payne, J. L. Ecologically diverse clades dominate the oceans via extinction resistance. Science 367, 1035–1038 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
Pimiento, C. et al. Functional diversity of marine megafauna in the Anthropocene. Sci. Adv. 6, eaay7650 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
R Core Team. R: A Language and Environment for Statistical Computing v.4.0.2 (R Foundation for Statistical Computing, 2020).Zamora, S. & Rahman, I. A. Deciphering the early evolution of echinoderms with Cambrian fossils. Palaeontology 57, 1105–1119 (2014).
Google Scholar 
Cohen, K. M., Finney, S. C., Gibbard, P. L. & Fan, J.-X. The ICS International Chronostratigraphic Chart. Episodes 36, 199–204 (2013).
Google Scholar 
Bapst, D. W., Bullock, P. C., Melchin, M. J., Sheets, H. D. & Mitchell, C. E. Graptoloid diversity and disparity became decoupled during the Ordovician mass extinction. Proc. Natl Acad. Sci. USA 109, 3428–3433 (2012).CAS 
PubMed 
PubMed Central 

Google Scholar 
Bapst, D. W. A stochastic rate-calibrated method for time-scaling phylogenies of fossil taxa. Methods Ecol. Evol. 4, 724–733 (2013).
Google Scholar 
Bapst, D. W. Assessing the effect of time-scaling methods on phylogeny-based analyses in the fossil record. Paleobiology 40, 331–351 (2014).
Google Scholar 
Bapst, D. W. & Hopkins, M. J. Comparing cal3 and other a posteriori time-scaling approaches in a case study with the pterocephaliid trilobites. Paleobiology 43, 49–67 (2016).
Google Scholar 
Paradis, E. et al. ape: analyses of phylogenetics and evolution. R package version 5.4 https://cran.r-project.org/package=ape (2020).Yang, Z., Kumar, S. & Nei, M. A new method of inference of ancestral nucleotide and amino acid sequences. Genetics 141, 1641–1650 (1995).CAS 
PubMed 
PubMed Central 

Google Scholar 
Guillerme, T. et al. Disparities in the analysis of morphological disparity. Biol. Lett. 16, 20200199 (2020).PubMed 
PubMed Central 

Google Scholar 
Brusatte, S. L., Montanari, S., Yi, H.-Y. & Norell, M. A. Phylogenetic corrections for morphological disparity analysis: new methodology and case studies. Paleobiology 37, 1–22 (2011).
Google Scholar 
Guillerme, T. & Cooper, N. Time for a rethink: time sub-sampling methods in disparity-through-time analyses. Palaeontology 61, 481–493 (2018).
Google Scholar 
Wills, M. A. Crustacean disparity through the Phanerozoic: comparing morphological and stratigraphic data. Biol. J. Linn. Soc. 65, 455–500 (1998).
Google Scholar 
Lehmann, O. E. R., Ezcurra, M. D., Butler, R. J. & Lloyd, G. T. Biases with the Generalized Euclidean Distance measure in disparity analyses with high levels of missing data. Palaeontology 62, 837–849 (2019).
Google Scholar 
Hopkins, M. J. & John, K. S. A new family of dissimilarity metrics for discrete character matrices that include inapplicable characters and its importance for disparity studies. Proc. R. Soc. B 285, 20181784 (2018).PubMed 
PubMed Central 

Google Scholar 
Chessel, D., Dufour, A. B. & Thioulouse, J. ade4: analysis of ecological data–exploratory and Euclidean methods in environmental sciences. R package version 1.7-15 https://cran.r-project.org/package=ade4 (2020).Maire, E., Grenouillet, G., Brosse, S. & Villéger, S. How many dimensions are needed to accurately assess functional diversity? A pragmatic approach for assessing the quality of functional spaces. Glob. Ecol. Biogeogr. 24, 728–740 (2015).
Google Scholar 
Mouchet, M. A., Villéger, S., Mason, N. W. H. & Mouillot, D. Functional diversity measures: an overview of their redundancy and their ability to discriminate community assembly rules. Funct. Ecol. 24, 867–876 (2010).
Google Scholar 
Novack-Gottshall, P. M. General models of ecological diversification. II. Simulations and empirical applications. Paleobiology 42, 209–239 (2016).Villéger, S., Novack-Gottshall, P. M. & Mouillot, D. The multidimensionality of the niche reveals functional diversity changes in benthic marine biotas across geological time. Ecol. Lett. 14, 561–568 (2011).PubMed 
PubMed Central 

Google Scholar 
Mouillot, D. et al. The dimensionality and structure of species trait spaces. Ecol. Lett. https://doi.org/10.1111/ele.13778 (2021).Wills, M. A. in Fossils, Phylogeny, and Form: An Analytical Approach (eds Adrain, J. M. et al.) 55–143 (Kluwer Academic/Plenum, 2001).Foote, M. Discordance and concordance between morphological and taxonomic diversity. Paleobiology 19, 185–204 (1993).
Google Scholar 
Ciampaglio, C. N., Kemp, M. & McShea, D. W. Detecting changes in morphospace occupation patterns in the fossil record: characterization and analysis of measures of disparity. Paleobiology 27, 695–715 (2001).
Google Scholar 
Hopkins, M. J. & Gerber, S. in Evolutionary Developmental Biology: A Reference Guide (eds de la Rosa, L. N. & Müller, G.) 1–12 (Springer, 2017).Villéger, S., Mason, N. W. H. & Mouillot, D. New multidimensional functional diversity indices for a multifaceted framework in functional ecology. Ecology 89, 2290–2301 (2008).PubMed 
PubMed Central 

Google Scholar 
Anderson, M. J., Ellingsen, K. E. & McArdle, B. H. Multivariate dispersion as a measure of beta diversity. Ecol. Lett. 9, 683–693 (2006).PubMed 
PubMed Central 

Google Scholar 
Laliberté, E. & Legendre, P. A distance-based framework for measuring functional diversity from multiple traits. Ecology 91, 299–305 (2010).PubMed 
PubMed Central 

Google Scholar 
Novack-Gottshall, P. ecospace: simulating community assembly and ecological diversification using ecospace frameworks. R package version 1.0.1 https://cran.r-project.org/package=ecospace (2015).Laliberté, E. & Shipley, B. FD: measuring functional diversity from multiple traits, and other tools for functional ecology. R package version 1.0-12 https://cran.r-project.org/package=FD (2014).Dineen, A. A., Roopnarine, P. D. & Fraiser, M. L. Ecological continuity and transformation after the Permo-Triassic mass extinction in northeastern Panthalassa. Biol. Lett. 15, 20180902 (2019).PubMed 
PubMed Central 

Google Scholar 
Efron, B. Bootstrap methods: another look at the jackknife. Ann. Stat. 7, 1–26 (1979).
Google Scholar 
Kowalewski, M. & Novack-Gottshall, P. in Quantitative Methods in Paleobiology. Short Courses in Paleontology 16 (eds Alroy, J. & Hunt, G.) 19–54 (Paleontological Society and Paleontological Research Institute, 2010).Lingoes, J. C. Some boundary conditions for a monotone analysis of symmetric matrices. Psychometrika 36, 195–203 (1971).
Google Scholar 
Legendre, P. & Anderson, M. J. Distance-based redundancy analysis: testing multispecies responses in multifactorial ecological experiments. Ecol. Monogr. 69, 1–24 (1999).
Google Scholar 
Burnham, K. P. & Anderson, D. R. Model Selection and Multi-Model Inference: A Practical Information-Theoretic Approach 2nd edn (Springer, 2002). More

Chemicals and reagentsEOs of basil (Ocimum basilicum L.), bergamot (Citrus bergamia Risso & Poit), camphor [Cinnamomum camphora (L.) J. Presl.], cinnamon (Cinnamomum zeylanicum Blume), citronella [Cymbopogon nardus (L.) Rendle], clove (Eugenia caryophyllus Wight), eucalyptus (Eucalyptus globulus Labill.), jasmine (Jasminum officinale L.), lavender (Lavandula angustifolia Mill.), lemon grass [Cymbopogan citratus (DC.) Stapf], mentha (Mentha piperita L.), rosemary (Rosmarinus officinalis L.), patchouli (Pogostemon patchouli Benth), and wild turmeric (Curcuma aromatica Salisb.) were procured from Talent Technologies (Talent Technologies, Kanpur, India). Acetylcholinesterase (AChE) activity assay kit, Anti-OBP2A antibody, ELISA kits, 1,1-diphenyl-2-picrylhydrazyl (DPPH), radioimmunoprecipitation (RIPA) buffer and phosphate buffer saline (PBS) were purchased from Sigma Aldrich (Sigma Aldrich Chemical Co., St. Luis, USA). TRPV1 antibody was purchased from Santa Cruz (Santa Cruz, California, USA). 1-chloro-2,4-dinitrobenzene (CDNB) was purchased from Cayman (Cayman Chemical Company, Michigan, USA). Human normal lung cell line (L-132) was obtained from the National Centre for Cell Sciences (NCCS), Pune, India. High performance liquid chromatography (HPLC) grade acetone was purchased from Merck (Merck Pvt. Ltd., Mumbai, India). All other chemicals used were of the highest analytical grade available.Test insects5–7 days old adult female Ae. albopictus mosquitoes were housed at the laboratory insectary, Division of Pharmaceutical Technology, Defence Research Laboratory, Tezpur, Assam, India. Mosquitoes were reared by maintaining temperature at 27 ± 2 °C, relative humidity: 75 ± 5% RH and 14L:10D h of light–dark alternative cycles in standard-sized wooden cages (75 cm × 60 cm × 60 cm) with a sleeve opening on one side as described previously63. 10% sucrose solution ad libitum were provided for nourishment. Before testing, the mosquitoes were starved for 24 h.Screening of EOsDose response study was performed to evaluate the best oils among the fourteen EOs. This study was approved (approval number: 032/2021TMCH, 28/08/2018) by the Institutional Human Ethical Committee (IHEC), of the Tezpur Medical College & Hospital (TMCH), Tezpur, Assam, India, and all experiments were performed in accordance with relevant guidelines and regulations. Five volunteers are chosen, not allergic to mosquito bite and all volunteers provided written informed consent. A volunteer’s thigh was marked according to the door opening hole of the K&D module as described by Klun and Debboun64. It is made of Plexiglas and the base of the rectangular cage (26 cm × 5 cm × 5 cm) has six holes, each with rectangular 3 × 4 cm holes that are opened and closed by a sliding door (Supplementary Fig. S8: Provide the photograph of K&D module). The flexor region of the forearms of a human volunteer was outlined with four rectangular (3 cm × 4 cm) test areas. A volume of 25 µL of each concentration of the EOs in soybean oil (40, 4 and 0.4 µg/cm2) and 25 µL of the soybean oil (diluent) as control was applied to the marked areas. After air drying for 5 min, a K&D module with matching cut outs in its floor was placed over the treated areas, containing five nulliparous 5–7 days old female mosquitoes in each hole. The doors of the cells were opened and the number of mosquitoes biting in each cell was recorded within a 2 min exposure, after which the doors were closed. After completion of each observation, mosquitoes were freed by opening cells of the K&D module in a sleeved screened cage. For each test, fresh sets of mosquitoes are used. Five replications for each test were carried out. The efficacy of EOs were determined by the percentage repellency against mosquitoes, using the formula or Eq. (2) described by WHO46.$$% ;{text{repellency}} = frac{C – T}{C} times 100$$
(2)
where, C is the number of mosquitoes landing, or biting at the control area; T is the number of mosquitoes landing or biting at the treated area.Fourier transform-infra red spectroscopy (FT-IR)Study of chemical compatibility for each formulation ingredients are necessary. All formulation ingredients possess specific value of vibrational frequency and have varied functional groups in their chemical structures. For compatibility study, each EOs, excipients to be used in cream formulation, and their physical mixture was placed one by one over the sample plate of the FT-IR instrument (Bruker, ALPHA, Billerica, MA, USA). The covering probe was placed over the sample and IR spectra was obtained over a wavelength of 2.5–25 μm at room temperature. Functional groups possessed by each individual ingredient should be identical in their physical mixture which confirms their compatibility37.Thermogravimetric analysis (TGA)The thermal behaviour of citronella oil, clove oil, lemon grass oil, their mixture and EO-MRC were evaluated using a thermal analyser (TG 209 F1 Libra®, NETZSCH-Gerätebau GmbH, 95100 Selb, Germany). Approximately about 10 mg sample weight was placed in the crucible each time. Nitrogen was used as a shielding gas. Heating program was fixed as 30–600 °C at a rate of 10 °C/min.Formulation development and optimizationFor optimization, a 17-run, 3-factor, 3-level Box-Behnken design (BBD) was utilized. A second order polynomial model was constructed by quadratic response surface methodology (RSM) using Design-Expert software (Version 6.0.8, Stat-Ease Inc., USA). Total seventeen formulations were obtained using EO concentrations as dependent variables against complete protection time (CPT) as independent variable or response variable. Analysis of variance (ANOVA) was performed using the same software to obtain the most effective formulation.Preparation of creamPhase inversion temperature method was applied for the preparation of EO-based mosquito repellent cream (EO-MRC). About 50 g cream sample was prepared in order to get enough for performing the various qualitative and quantitative assay. The oil phase (phase B) was prepared by dissolving the oil soluble excipients, except phase A (mosquito repellent active ingredients) under mild heating at 200 rpm in a hot magnetic plate stirrer (Magnetic Stirrer IKA RCT basic) and heated to 65 °C. The aqueous phase was prepared by mixing various aqueous soluble ingredients (phase C) under gentle heating and stirring. Temperature of the aqueous phase was raised to 65 °C. Phase A was gently added to the oil phase at a stirring speed of 200 rpm and 55 ± 2 °C. The mixture was then emulsified by adding phase C slowly and kept for 1 h at a stirring rate of 800 rpm and 60 ± 2 °C. The formulated EO-MRC was then kept for natural cooling.Efficacy assessmentCPT of the developed cream (EO-MRC) formulation was carried out by arm in cage bioassay. 1 mL EO-MRC was applied to ≈ 600 cm2 area of the forearm skin between the wrist and elbow and 1 mL of the 12% N, N-di ethyl benzamide (DEBA) based marketed cream (DBMC) was compared on the other arm. Two mosquito cages (size: 40 × 40 × 40 cm) each containing 200–250 non-blood-fed female Ae. Albopictus were used. One cage is designated for testing the EO-MRC and the other for the positive control (DBMC). During testing, hands were protected by surgical gloves for which the mosquitoes cannot bite while the volunteer avoids movement of the arm. EO-MRC and DBMC treated arms were exposed for 3 min at 30 min intervals to determine landing and/or probing activity. A single landing or probing of mosquito within a 3 min test interval concludes the test. CPT was calculated as the time (min) required for the first mosquito landing or probing after repellent application to the treated area. The median CPT and confidence intervals were estimated from the Kaplan–Meier Survival Function46.Efficacy was correlated with DEBA based marketed cream (DBMC). The inclusion of the specific commercial product DBMC is for comparison and does not constitute any recommendations.CharacterizationGas chromatography-mass spectroscopy (GC–MS)Qualitative studyDifferent chemical components in fourteen EOs and the selected blend were identified by a GC–MS system of Agilent Technologies (5301 Stevens Creek Blvd. Santa Clara, CA 95051, United States). Test sample concentration of 500 μg/mL was prepared in GC grade acetone. A sample volume of 1 μL was introduced into the injector held at 250 °C. Oven temperature of 40–300 °C was programmed at 20 °C/min. Helium was used as carrier gas at flow rate 1 mL/min. The injector and detector temperature were set at 250 °C and 230 °C (quad) and 150 °C (core) respectively37. Standard C7–C30 saturated alkanes were purchased from Sigma Aldrich Chemicals Co., St. Louis, USA. Retention indices (RI) of the identified components were determined for identification of the detected components.% Assay by GC–MS studyCalibration samples of eugenol and citronellol were prepared by dissolving an appropriate amount in GC grade acetone to get concentrations of 62.5 μg/mL, 125 μg/mL, 250 μg/mL and 500 μg/mL. Test samples of EO-MRC, clove oil and citronella oil were prepared by dissolving a required amount in acetone to quantify the EO components in the final formulation. A sample volume of 1 μL was introduced into the injector as described in ‘Qualitative study’ section.Physicochemical parametersPhysical parameters of the EO-MRC and placebo formulations were determined in order to establish aesthetic compliance and consumer acceptability. To determine the viscosity, a programmable viscometer was used (Model: DV2T, Ametek Brookfield, Middleboro, MA, USA); combined with software Rheo3000, version 1.2.2019.1 [R]. Sample volume was fixed at 30 g and viscosities were determined at 10 rpm for 40 s at room temperature using a T-Bar spindle (B-92) (Helipath spindle set, Brookfield Engineering Labs. Inc). Density was determined by using a pycnometer. pH of EO-MRC was checked by using digital pH meter (Labman Scientific instruments, Tamil Nadu, India).Spread ability of EO-MRC was determined as per the method reported earlier by Sabale65. In brief, 1 g of EO-MRC was placed on 1 cm2 pre-marked circular area on the glass slide (7.5 cm × 2.5 cm). EO-MRC was compressed using another glass slide placed from edge to centre of primary slide. 200 g of commercial weight was placed on the set up and allowed the gel to spread for the period of 1 min. The spread diameter was calculated with the aid of graph paper and spread ability was evaluated using formula expressed as Eq. (3):$$mathrm{Spread, ability}=mathrm{m}times frac{mathrm{l}}{mathrm{t}}$$
(3)
where, m is the commercial weight placed on the setup; l is the length of cream spread; and t is the time.Safety assessmentCytotoxicity by MTT assayThe reduction of tetrazolium salts is now widely accepted as a reliable way to examine cell proliferation. The yellow tetrazolium MTT (3-(4,5-dimethylthiazolyl-2)-2,5-diphenyltetrazolium bromide) is reduced by metabolically active cells, in part by the action of dehydrogenase enzymes, to generate reducing equivalents such as NADH and NADPH. With the help of spectrophotometric means, the resulting intracellular purple formazan can be quantified. The assay measures the cell proliferation rate and conversely, when metabolic events cause apoptosis or necrosis, the reduction in cell viability66.Cells cultured in T-25 flasks were trypsinized and aspirated into a 5 mL centrifuge tube. Cell pellet was obtained by centrifugation at 3000 rpm. The cell count was adjusted, using DMEM HG medium, such that 200 μL of suspension contained approximately 10,000 cells. To each well of the 96 well microtiter plate, 200 μL of the cell suspension was added and the plate was incubated at 37 ℃ and 5% CO2 atmosphere for 24 h. After 24 h, the spent medium was aspirated. 200 μL of different test concentrations viz. 62 µg/mL, 125 µg/mL, 250 µg/mL, 500 µg/mL, and 1000 µg/mL, of EO-MRC were added to the respective wells. The plate was then incubated at 37 °C and 5% CO2 atmosphere for 24 h. The plate was removed from the incubator and the drug containing media was aspirated. 200 μL of medium containing10% MTT reagent was then added to each well to get a final concentration of 0.5 mg/mL and the plate was incubated at 37 ℃ and 5% CO2 atmosphere for 3 h. Without disturbing the crystals formed in the wells, culture medium was completely removed. 100 μL of solubilisation solution (DMSO) was added to each well and the plate was then gently shake in a rocking shaker (ROCKYMAX™, Tarsons, Kolkata, India) to solubilize the formed formazan. The absorbance was measured at a wavelength of 570 nm and also at 630 nm using a microplate reader. The percentage growth inhibition was calculated and concentration of EO-MRC needed to inhibit cell growth by 50% (IC50) was generated from the dose–response curve for the cell line.Animals and ethics statementAll experimenting protocols using animal were performed according to the “Principles of Laboratory Animal care” (NIH publication 85–23, revised 1985) and approved by the Institutional Animal Ethical Committee (IAEC) of Defence Research Laboratory (DRL), Tezpur, Assam, India (approval no. CPCSEA/DRL/Protocol no. 3, 20/06/2018). All studies involving animals are reported in accordance with the ARRIVE guidelines for reporting experiments involving animals67. All efforts were made during the study period to minimize the suffering of animals and to reduce the number of animals used.5–8 weeks old, about 210–250 g of male healthy adult Wistar rats (Rattus norvegicus) and young and healthy New Zealand albino rabbits (Oryctolagus cuniculus) were obtained from the institutional animal housing facility and allowed to acclimatize for 7 days prior to the study. Standard food and purified water ad libitum were provided in clean and hygienic condition at 22–25 ℃, 40–70% RH with 12 h light–dark cycles.Acute dermal irritation studyAcute dermal irritation study was conducted on healthy New Zealand albino rabbits following the OECD test guidelines 40468. Approximately 24 h before the test, fur was removed from the dorsal area of the trunk. 0.5 g EO-MRC, was directly applied to the skin and after 4 h exposure period, residual EO-MRC was removed by using water without disturbing the integrity of the epidermis and examined for signs of erythema and oedema, at 60 min, and then at 24 h, 48 h and 72 h after EO-MRC removal. Dermal reactions are graded and recorded according to the grades in the Table 8. As per the method described by Banerjee et al.69; primary irritation index (PII) was calculated. Further, we have followed the Draize method of classification for PII scoring as non-irritant (if PII  More

1.Ceballos, G., Ehrlich, P. R. & Raven, P. H. Vertebrates on the brink as indicators of biological annihilation and the sixth mass extinction. Proc. Natl Acad. Sci. USA 117, 13596–13602 (2020).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
2.Committee, T. I. R. L. The IUCN Red List of Threatened Species – Strategic Plan 2017-2020. Tech. Rep., IUCN (2017).3.Witmer, G. W. Wildlife population monitoring: some practical considerations. Wild. Res. 32, 259–263 (2005).
Google Scholar 
4.McEvoy, J. F., Hall, G. P. & McDonald, P. G. Evaluation of unmanned aerial vehicle shape, flight path and camera type for waterfowl surveys: disturbance effects and species recognition. PeerJ 4, e1831 (2016).PubMed 
PubMed Central 

Google Scholar 
5.Burghardt, G. M. et al. Perspectives–minimizing observer bias in behavioral studies: a review and recommendations. Ethology 118, 511–517 (2012).
Google Scholar 
6.Giese, M. Effects of human activity on Adelie penguin Pygoscelis adeliae breeding success. Biol. Conserv. 75, 157–164 (1996).
Google Scholar 
7.Köndgen, S. et al. Pandemic human viruses cause decline of endangered great apes. Curr. Biol. 18, 260–264 (2008).PubMed 

Google Scholar 
8.Weissensteiner, M. H., Poelstra, J. W. & Wolf, J. B. W. Low-budget ready-to-fly unmanned aerial vehicles: an effective tool for evaluating the nesting status of canopy-breeding bird species. J. Avian Biol. 46, 425–430 (2015).
Google Scholar 
9.Sasse, D. B. Job-related mortality of wildlife workers in the united states, 1937–2000. Wildl. Soc. Bull. 31, 1015–1020 (2003).10.Kays, R., Crofoot, M. C., Jetz, W. & Wikelski, M. Terrestrial animal tracking as an eye on life and planet. Science 348, aaa2478 (2015).11.Altmann, J. Observational study of behavior: sampling methods. Behaviour 49, 227–266 (1974).CAS 
PubMed 

Google Scholar 
12.Hodgson, J. C. et al. Drones count wildlife more accurately and precisely than humans. Methods Ecol. Evolution 9, 1160–1167 (2018).
Google Scholar 
13.Betke, M. et al. Thermal imaging reveals significantly smaller Brazilian free-tailed bat colonies than previously estimated. J. Mammal. 89, 18–24 (2008).
Google Scholar 
14.Rollinson, C. R. et al. Working across space and time: nonstationarity in ecological research and application. Front. Ecol. Environ. 19, 66–72 (2021).
Google Scholar 
15.Junker, J. et al. A severe lack of evidence limits effective conservation of the world’s primates. BioScience 70, 794–803 (2020).PubMed 
PubMed Central 

Google Scholar 
16.Sherman, J., Ancrenaz, M. & Meijaard, E. Shifting apes: Conservation and welfare outcomes of Bornean orangutan rescue and release in Kalimantan, Indonesia. J. Nat. Conserv. 55, 125807 (2020).
Google Scholar 
17.O’Donoghue, P. & Rutz, C. Real-time anti-poaching tags could help prevent imminent species extinctions. J. Appl. Ecol. 53, 5–10 (2016).PubMed 

Google Scholar 
18.Lahoz-Monfort, J. J. & Magrath, M. J. L. A comprehensive overview of technologies for species and habitat monitoring and conservation. BioScience biab073. https://academic.oup.com/bioscience/advance-article/doi/10.1093/biosci/biab073/6322306 (2021).19.Gottschalk, T., Huettmann, F. & Ehlers, M. Thirty years of analysing and modelling avian habitat relationships using satellite imagery data: a review. Int. J. Remote Sens. 26, 2631–2656 (2005).
Google Scholar 
20.Steenweg, R. et al. Scaling-up camera traps: monitoring the planet’s biodiversity with networks of remote sensors. Front. Ecol. Environ. 15, 26–34 (2017).
Google Scholar 
21.Hausmann, A. et al. Social media data can be used to understand tourists’ preferences for nature-based experiences in protected areas. Conserv. Lett. 11, e12343 (2018).
Google Scholar 
22.Sugai, L. S. M., Silva, T. S. F., Ribeiro, J. W. & Llusia, D. Terrestrial passive acoustic monitoring: review and perspectives. BioScience 69, 15–25 (2018).
Google Scholar 
23.Wikelski, M. et al. Going wild: what a global small-animal tracking system could do for experimental biologists. J. Exp. Biol. 210, 181–186 (2007).PubMed 

Google Scholar 
24.Belyaev, M. Y. et al. Development of technology for monitoring animal migration on Earth using scientific equipment on the ISS RS. in 2020 27th Saint Petersburg International Conference on Integrated Navigation Systems (ICINS), 1–7 (IEEE, 2020).25.Harel, R., Loftus, J. C. & Crofoot, M. C. Locomotor compromises maintain group cohesion in baboon troops on the move. Proc. R. Soc. B 288, 20210839 (2021).PubMed 
PubMed Central 

Google Scholar 
26.Farley, S. S., Dawson, A., Goring, S. J. & Williams, J. W. Situating ecology as a big-data science: current advances, challenges, and solutions. BioScience 68, 563–576 (2018).
Google Scholar 
27.Lasky, M. et al. Candid critters: Challenges and solutions in a large-scale citizen science camera trap project. Citizen Science: Theory and Practice 6, https://doi.org/10.5334/cstp.343 (2021).28.Hastie, T., Tibshirani, R. & Friedman, J. The Elements of Statistical Learning: Data Mining, Inference, and Prediction (Springer, 2001).29.Christin, S., Hervet, É. & Lecomte, N. Applications for deep learning in ecology. Methods Ecol. Evolution 10, 1632–1644 (2019).
Google Scholar 
30.Kwok, R. Ai empowers conservation biology. Nature 567, 133–135 (2019).ADS 
CAS 
PubMed 

Google Scholar 
31.Kwok, R. Deep learning powers a motion-tracking revolution. Nature 574, 137–139 (2019).ADS 
CAS 
PubMed 

Google Scholar 
32.LeCun, Y., Bengio, Y. & Hinton, G. Deep learning. Nature 521, 436–444 (2015).ADS 
CAS 
PubMed 

Google Scholar 
33.Pichler, M., Boreux, V., Klein, A.-M., Schleuning, M. & Hartig, F. Machine learning algorithms to infer trait-matching and predict species interactions in ecological networks. Methods Ecol. Evolution 11, 281–293 (2020).
Google Scholar 
34.Knudby, A., LeDrew, E. & Brenning, A. Predictive mapping of reef fish species richness, diversity and biomass in Zanzibar using IKONOS imagery and machine-learning techniques. Remote Sens. Environ. 114, 1230–1241 (2010).ADS 

Google Scholar 
35.Rey, N., Volpi, M., Joost, S. & Tuia, D. Detecting animals in African savanna with UAVs and the crowds. Remote Sens. Environ. 200, 341–351 (2017).ADS 

Google Scholar 
36.Beery, S., Morris, D. & Yang, S. Efficient pipeline for camera trap image review. in Proceedings of the Workshop Data Mining and AI for Conservation, Conference for Knowledge Discovery and Data Mining (2019).37.Kellenberger, B., Marcos, D. & Tuia, D. When a few clicks make all the difference: improving weakly-supervised wildlife detection in UAV images. in Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition Workshops (2019).38.Schofield, D. et al. Chimpanzee face recognition from videos in the wild using deep learning. Sci. Adv. 5, eaaw0736 (2019).ADS 
PubMed 
PubMed Central 

Google Scholar 
39.Ditria, E. M. et al. Automating the analysis of fish abundance using object detection: optimizing animal ecology with deep learning. Front. Mar. Sci. 7, 429 (2020).
Google Scholar 
40.Kellenberger, B., Veen, T., Folmer, E. & Tuia, D. 21 000 birds in 4.5 h: efficient large-scale seabird detection with machine learning. Remote Sens. Ecol. Conserv. 7, 445–460 (2021).
Google Scholar 
41.Ahumada, J. A. et al. Wildlife insights: a platform to maximize the potential of camera trap and other passive sensor wildlife data for the planet. Environ. Conserv. 47, 1–6 (2020).MathSciNet 

Google Scholar 
42.Eikelboom, J. A. J. et al. Improving the precision and accuracy of animal population estimates with aerial image object detection. Methods Ecol. Evolution 10, 1875–1887 (2019).
Google Scholar 
43.Weinstein, B. G. A computer vision for animal ecology. J. Anim. Ecol. 87, 533–545 (2018).PubMed 

Google Scholar 
44.Valletta, J. J., Torney, C., Kings, M., Thornton, A. & Madden, J. Applications of machine learning in animal behaviour studies. Anim. Behav. 124, 203–220 (2017).
Google Scholar 
45.Peters, D. P. C. et al. Harnessing the power of big data: infusing the scientific method with machine learning to transform ecology. Ecosphere 5, art67 (2014).
Google Scholar 
46.Yu, Q. et al. Study becomes insight: ecological learning from machine learning. Methods Ecol. Evol. 12, 2117–2128 (2021).47.Lucas, T. C. D. A translucent box: interpretable machine learning in ecology. Ecol. Monogr. 90, https://doi.org/10.1002/ecm.1422 (2020).48.Reichstein, M. et al. Deep learning and process understanding for data-driven Earth system science. Nature 566, 195–204 (2019).ADS 
CAS 
PubMed 

Google Scholar 
49.Camps-Valls, G., Tuia, D., Zhu, X. X. & Reichstein, M. Deep Learning for the Earth Sciences: A Comprehensive Approach to Remote Sensing, Climate Science and Geosciences (Wiley & Sons, 2021).50.Karpatne, A. et al. Theory-guided data science: A new paradigm for scientific discovery from data. IEEE Trans. Knowl. Data Eng. 29, 2318–2331 (2017).
Google Scholar 
51.Oliver, R. Y., Meyer, C., Ranipeta, A., Winner, K. & Jetz, W. Global and national trends, gaps, and opportunities in documenting and monitoring species distributions. PLoS Biol 19, e3001336 https://doi.org/10.1371/journal.pbio.3001336 (2021).52.Beery, S., Wu, G., Rathod, V., Votel, R. & Huang, J. Context R-CNN: long term temporal context for per-camera object detection. in 2020 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), 13075–13085 (2020).53.Norouzzadeh, M. S. et al. Automatically identifying, counting, and describing wild animals in camera-trap images with deep learning. Proc. Natl Acad. Sci. USA 115, E5716–E5725 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
54.Schneider, S., Taylor, G. W., Linquist, S. & Kremer, S. C. Past, present and future approaches using computer vision for animal re-identification from camera trap data. Methods Ecol. Evolution 10, 461–470 (2019).
Google Scholar 
55.Beery, S., Van Horn, G. & Perona, P. Recognition in terra incognita. in 2018 European Conference on Computer Vision (ECCV), 456–473 (2018).56.Sugai, L. S. M., Silva, T. S. F., Ribeiro Jr, J. W. & Llusia, D. Terrestrial passive acoustic monitoring: review and perspectives. BioScience 69, 15–25 (2019).
Google Scholar 
57.Wrege, P. H., Rowland, E. D., Keen, S. & Shiu, Y. Acoustic monitoring for conservation in tropical forests: examples from forest elephants. Methods Ecol. Evolution 8, 1292–1301 (2017).
Google Scholar 
58.Desjonquères, C., Gifford, T. & Linke, S. Passive acoustic monitoring as a potential tool to survey animal and ecosystem processes in freshwater environments. Freshw. Biol. 65, 7–19 (2020).
Google Scholar 
59.Davis, G. E. et al. Long-term passive acoustic recordings track the changing distribution of North Atlantic right whales (eubalaena glacialis) from 2004 to 2014. Sci. Rep. 7, 1–12 (2017).
Google Scholar 
60.Wood, C. M. et al. Detecting small changes in populations at landscape scales: a bioacoustic site-occupancy framework. Ecol. Indic. 98, 492–507 (2019).
Google Scholar 
61.Kahl, S., Wood, C. M., Eibl, M. & Klinck, H. Birdnet: a deep learning solution for avian diversity monitoring. Ecol. Inform. 61, 101236 (2021).
Google Scholar 
62.Stowell, D., Wood, M. D., Pamuła, H., Stylianou, Y. & Glotin, H. Automatic acoustic detection of birds through deep learning: the first bird audio detection challenge. Methods Ecol. Evolution 10, 368–380 (2019).
Google Scholar 
63.Ford, J. K. B. in Encyclopedia of Marine Mammals 253–254 (Elsevier, 2018).64.Hughey, L. F., Hein, A. M., Strandburg-Peshkin, A. & Jensen, F. H. Challenges and solutions for studying collective animal behaviour in the wild. Philos. Trans. R. Soc. B: Biol. Sci. 373, 20170005 (2018).
Google Scholar 
65.Williams, H. J. et al. Optimizing the use of biologgers for movement ecology research. J. Anim. Ecol. 89, 186–206 (2020).PubMed 

Google Scholar 
66.Korpela, J. et al. Machine learning enables improved runtime and precision for bio-loggers on seabirds. Commun. Biol. 3, 1–9 (2020).
Google Scholar 
67.Yu, H. An evaluation of machine learning classifiers for next-generation, continuous-ethogram smart trackers. Mov. Ecol. 9, 14 (2021).
Google Scholar 
68.Browning, E. et al. Predicting animal behaviour using deep learning: GPS data alone accurately predict diving in seabirds. Methods Ecol. Evolution 9, 681–692 (2018).
Google Scholar 
69.Liu, Z. Y.-C. et al. Deep learning accurately predicts white shark locomotor activity from depth data. Anim. Biotelemetry 7, 1–13 (2019).
Google Scholar 
70.Wang, G. Machine learning for inferring animal behavior from location and movement data. Ecol. Inform. 49, 69–76 (2019).
Google Scholar 
71.Wijeyakulasuriya, D. A., Eisenhauer, E. W., Shaby, B. A. & Hanks, E. M. Machine learning for modeling animal movement. PLoS ONE 30, e0235750 (2020).72.Linchant, J., Lisein, J., Semeki, J., Lejeune, P. & Vermeulen, C. Are unmanned aircraft systems (UASs) the future of wildlife monitoring? A review of accomplishments and challenges. Mammal. Rev. 45, 239–252 (2015).
Google Scholar 
73.Hodgson, J. C., Baylis, S. M., Mott, R., Herrod, A. & Clarke, R. H. Precision wildlife monitoring using unmanned aerial vehicles. Sci. Rep. 6, 1–7 (2016).
Google Scholar 
74.Mathis, A. et al. DeepLabCut: markerless pose estimation of user-defined body parts with deep learning. Nat. Neurosci. 21, 1281–1289 (2018).CAS 
PubMed 

Google Scholar 
75.Graving, J. M. et al. DeepPoseKit, a software toolkit for fast and robust animal pose estimation using deep learning. Elife 8, e47994 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
76.Mathis, A., Schneider, S., Lauer, J. & Mathis, M. W. A primer on motion capture with deep learning: principles, pitfalls, and perspectives. Neuron 108, 44–65 (2020).CAS 
PubMed 

Google Scholar 
77.Kellenberger, B., Marcos, D. & Tuia, D. Detecting mammals in UAV images: best practices to address a substantially imbalanced dataset with deep learning. Remote Sens. Environ. 216, 139–153 (2018).ADS 

Google Scholar 
78.Kellenberger, B., Veen., T., Folmer, E. & Tuia, D. 21,000 birds in 4.5 hours: efficient large-scale seabird detection with machine learning. Remote Sens. Ecol. Conserv. https://doi.org/10.1002/rse2.200 (2021).79.Andrew, W., Greatwood, C. & Burghardt, T. Aerial animal biometrics: individual Friesian cattle recovery and visual identification via an autonomous UAV with onboard deep inference. in International Conference on Intelligent Robots and Systems (IROS) (2019).80.Schroeder, N. M., Panebianco, A., Gonzalez Musso, R. & Carmanchahi, P. An experimental approach to evaluate the potential of drones in terrestrial mammal research: a gregarious ungulate as a study model. R. Soc. open Sci. 7, 191482 (2020).ADS 
PubMed 
PubMed Central 

Google Scholar 
81.Bennitt, E., Bartlam-Brooks, H. L. A., Hubel, T. Y. & Wilson, A. M. Terrestrial mammalian wildlife responses to Unmanned Aerial Systems approaches. Sci. Rep. 9, 1–10 (2019).CAS 

Google Scholar 
82.Deneu, B., Servajean, M., Botella, C. & Joly, A. Evaluation of deep species distribution models using environment and co-occurrences. in International Conference of the Cross-Language Evaluation Forum for European Languages, 213–225 (Springer, 2019).83.Zhu, X. et al. Deep learning in remote sensing: A comprehensive review and list of resources. IEEE Geosci. Remote Sens. Mag. 5, 8–36 (2017).
Google Scholar 
84.Guirado, E., Tabik, S., Rivas, M. L., Alcaraz-Segura, D. & Herrera, F. Whale counting in satellite and aerial images with deep learning. Sci. Rep. 9, 1–12 (2019).CAS 

Google Scholar 
85.Duporge, I., Isupova, O., Reece, S., Macdonald, D. W. & Wang, T. Using very-high-resolution satellite imagery and deep learning to detect and count African elephants in heterogeneous landscapes. Remote Sens. Ecol. Conserv. https://doi.org/10.1002/rse2.195 (2020).86.Fretwell, P. T. & Trathan, P. N. Discovery of new colonies by Sentinel2 reveals good and bad news for emperor penguins. Remote Sens. Ecol. Conserv. https://doi.org/10.1002/rse2.176 (2020).87.Brodrick, P. G., Davies, A. B. & Asner, G. P. Uncovering ecological patterns with convolutional neural networks. Trends Ecol. Evolution 34, 734–745 (2019).
Google Scholar 
88.Audebert, N., Le Saux, B. & Lefèvre, S. Deep learning for classification of hyperspectral data: a comparative review. IEEE Geosci. Remote Sens. Mag. 7, 159–173 (2019).
Google Scholar 
89.McKinley, D. C. et al. Citizen science can improve conservation science, natural resource management, and environmental protection. Biol. Conserv. 208, 15–28 (2017).
Google Scholar 
90.Wäldchen, J. & Mäder, P. Machine learning for image based species identification. Methods Ecol. Evolution 9, 2216–2225 (2018).MATH 

Google Scholar 
91.Torney, C. J. et al. A comparison of deep learning and citizen science techniques for counting wildlife in aerial survey images. Methods Ecol. Evolution 10, 779–787 (2019).
Google Scholar 
92.Parham, J., Crall, J., Stewart, C., Berger-Wolf, T. & Rubenstein, D. I. Animal population censusing at scale with citizen science and photographic identification. in AAAI Spring Symposium-Technical Report (2017).93.Kühl, H. S. & Burghardt, T. Animal biometrics: quantifying and detecting phenotypic appearance. Trends Ecol. Evolution 28, 432–441 (2013).
Google Scholar 
94.Yu, X. et al. Automated identification of animal species in camera trap images. EURASIP J. Image Video Process. 2013, 1–10 (2013).ADS 

Google Scholar 
95.Mac Aodha, O. et al. Bat detective–deep learning tools for bat acoustic signal detection. PLoS Computat. Biol. 14, e1005995 (2018).
Google Scholar 
96.Schindler, F. & Steinhage, V. Identification of animals and recognition of their actions in wildlife videos using deep learning techniques. Ecol. Inform. 61, 101215 (2021).97.Avise, J. C. Molecular Markers, Natural History and Evolution (Springer Science & Business Media, 2012).98.Vidal, M., Wolf, N., Rosenberg, B., Harris, B. P. & Mathis, A. Perspectives on Individual Animal Identification from Biology and Computer Vision. Integr. Comp. Biol. 61, 900–916 https://doi.org/10.1093/icb/icab107 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
99.Berger-Wolf, T. Y. et al. Wildbook: crowdsourcing, computer vision, and data science for conservation. Preprint at https://arxiv.org/abs/1710.08880 (2017).100.Parham, J. et al. An animal detection pipeline for identification. in IEEE Winter Conference on Applications of Computer Vision (WACV), 1075–1083 (IEEE, 2018).101.Weideman, H. et al. Extracting identifying contours for African elephants and humpback whales using a learned appearance model. in Proceedings of the IEEE/CVF Winter Conference on Applications of Computer Vision (WACV) (2020).102.Brust, C.-A. et al. Towards automated visual monitoring of individual gorillas in the wild. in 2017 IEEE International Conference on Computer Vision Workshops (ICCVW), 2820–2830 (2017).103.Li, S., Li, J., Tang, H., Qian, R. & Lin, W. ATRW: a benchmark for Amur tiger re-identification in the wild. in 2020 ACM International Conference on Multimedia, 2590–2598 (2020).104.Bendale, A. & Boult, T. E. Towards open set deep networks. in 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), 1563–1572 (2016).105.Mathis, M. W. & Mathis, A. Deep learning tools for the measurement of animal behavior in neuroscience. Curr. Opin. Neurobiol. 60, 1–11 (2020).CAS 
PubMed 

Google Scholar 
106.Sanakoyeu, A., Khalidov, V., McCarthy, M. S., Vedaldi, A. & Neverova, N. Transferring dense pose to proximal animal classes. in 2020 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), 5233–5242 (2020).107.Zuffi, S., Kanazawa, A., Jacobs, D. W. & Black, M. J. 3D menagerie: modeling the 3D shape and pose of animals. in 2017 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), 6365–6373 (2017).108.Biggs, B., Roddick, T., Fitzgibbon, A. & Cipolla, R. Creatures great and smal: recovering the shape and motion of animals from video. in 2018 Asian Conference on Computer Vision (ACCV), 3–19 (Springer, 2018).109.Biggs, B., Boyne, O., Charles, J., Fitzgibbon, A. & Cipolla, R. Who left the dogs out? 3D animal reconstruction with expectation maximization in the loop. in 2020 European Conference on Computer Vision (ECCV), 195–211 (Springer, 2020).110.Zuffi, S., Kanazawa, A., Berger-Wolf, T. & Black, M. J. Three-D safari: learning to estimate zebra pose, shape, and texture from images” in the wild”. in 2019 IEEE International Conference on Computer Vision (ICCV), 5359–5368 (2019).111.Wang, Y., Kolotouros, N., Daniilidis, K. & Badger, M. Birds of a feather: capturing avian shape models from images. in 2021 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), 14739–14749 (2021).112.Haalck, L., Mangan, M., Webb, B. & Risse, B. Towards image-based animal tracking in natural environments using a freely moving camera. J. Neurosci. methods 330, 108455 (2020).PubMed 

Google Scholar 
113.Pettorelli, N. et al. Satellite remote sensing for applied ecologists: opportunities and challenges. J. Appl. Ecol. 51, 839–848 (2014).
Google Scholar 
114.Davies, A. B., Tambling, C. J., Kerley, G. I. H. & Asner, G. P. Effects of vegetation structure on the location of lion kill sites in African thicket. PLoS ONE 11, e0149098 (2016).PubMed 
PubMed Central 

Google Scholar 
115.Froidevaux, J. S. P., Zellweger, F., Bollmann, K., Jones, G. & Obrist, M. K. From field surveys to LiDAR: shining a light on how bats respond to forest structure. Remote Sens. Environ. 175, 242–250 (2016).ADS 

Google Scholar 
116.Risse, B., Mangan, M., Stürzl, W. & Webb, B. Software to convert terrestrial LiDAR scans of natural environments into photorealistic meshes. Environ. Model. Softw. 99, 88–100 (2018).
Google Scholar 
117.Haalck, L. & Risse, B. Embedded dense camera trajectories in multi-video image mosaics by geodesic interpolation-based reintegration. in 2021 IEEE Winter Conference on Applications of Computer Vision (WACV), 1849–1858 (2021).118.Schonberger, J. L. & Frahm, J.-M. Structure-from-motion revisited. in 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), 4104–4113 (2016).119.Mur-Artal, R. & Tardós, J. D. ORB-SLAM2: an open-source SLAM system for monocular, stereo, and RGB-D cameras. IEEE Trans. Robot. 33, 1255–1262 (2017).
Google Scholar 
120.Kuppala, K., Banda, S. & Barige, T. R. An overview of deep learning methods for image registration with focus on feature-based approaches. Int. J. Image Data Fusion 11, 113–135 (2020).ADS 

Google Scholar 
121.Lisein, J., Linchant, J., Lejeune, P., Bouché, P. & Vermeulen, C. Aerial surveys using an unmanned aerial system (UAS): comparison of different methods for estimating the surface area of sampling strips. Tropical Conserv. Sci. 6, 506–520 (2013).
Google Scholar 
122.Wu, C. Critical configurations for radial distortion self-calibration. in 2014 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), 25–32 (2014).123.Ferrer, J., Elibol, A., Delaunoy, O., Gracias, N. & Garcia, R. Large-area photo-mosaics using global alignment and navigation data. in Mts/IEEE Oceans Conference, 1–9 (2007).124.Guisan, A. & Zimmermann, N. E. Predictive habitat distribution models in ecology. Ecol. Model. 135, 147–186 (2000).
Google Scholar 
125.Lehmann, A., Overton, J. M. & Austin, M. P. Regression models for spatial prediction: their role for biodiversity and conservation. Biodivers. Conserv. 11, 2085–2092 (2002).126.Breiman, L. Random forests. Mach. Learn. 45, 5–32 (2001).MATH 

Google Scholar 
127.Parravicini, V. et al. Global patterns and predictors of tropical reef fish species richness. Ecography 36, 1254–1262 (2013).
Google Scholar 
128.Smoliński, S. & Radtke, K. Spatial prediction of demersal fish diversity in the baltic sea: comparison of machine learning and regression-based techniques. ICES J. Mar. Sci. 74, 102–111 (2017).
Google Scholar 
129.Čandek, K., Čandek, U. P. & Kuntner, M. Machine learning approaches identify male body size as the most accurate predictor of species richness. BMC Biol. 18, 1–16 (2020).
Google Scholar 
130.Baltensperger, A. P. & Huettmann, F. Predictive spatial niche and biodiversity hotspot models for small mammal communities in Alaska: applying machine-learning to conservation planning. Landscape Ecol. 30, 681–697 (2015).131.Faisal, A., Dondelinger, F., Husmeier, D. & Beale, C. M. Inferring species interaction networks from species abundance data: a comparative evaluation of various statistical and machine learning methods. Ecol. Inform. 5, 451–464 (2010).
Google Scholar 
132.Van Horn, G. et al. The inaturalist species classification and detection dataset. in 2018 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), 8769–8778 (2018).133.Copas, K. et al. Training machines to improve species identification using GBIF-mediated datasets. in AGU Fall Meeting Abstracts, Vol. 2019, IN53C–0758 (2019).134.Lennox, R. J. et al. A novel framework to protect animal data in a world of ecosurveillance. BioScience 70, 468–476 (2020).
Google Scholar 
135.Strubell, E., Ganesh, A. & McCallum, A. Energy and policy considerations for deep learning in NLP. in Proceedings of the 57th Annual Meeting of the Association for Computational Linguistics, 3645–3650 (2019).136.Samek, W., Montavon, G., Vedaldi, A., Hansen, L. K. & Müller, K.-R. Explainable AI: Interpreting, Explaining and Visualizing Deep Learning, Vol. 11700 (Springer Nature, 2019).137.Swanson, A. et al. Snapshot Serengeti, high-frequency annotated camera trap images of 40 mammalian species in an African savanna. Sci. data 2, 1–14 (2015).
Google Scholar 
138.de Lutio, R. et al. Digital taxonomist: identifying plant species in community scientists’ photographs. ISPRS J. Photogramm. Remote Sens. 182, 112–121 (2021).139.Mac Aodha, O., Cole, E. & Perona, P. Presence-only geographical priors for fine-grained image classification. in Proceedings of the IEEE/CVF International Conference on Computer Vision, 9596–9606 (2019).140.Gurumurthy, S. et al. Exploiting Data and Human Knowledge for Predicting Wildlife Poaching. in Proceedings of the 1st ACM SIGCAS Conference on Computing and Sustainable Societies, 1–8, https://doi.org/10.1145/3209811.3209879 (ACM, 2018).141.Datta, S., Anderson, D., Branson, K., Perona, P. & Leifer, A. Computational neuroethology: a call to action. Neuron 104, 11–24 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
142.Joska, D. et al. AcinoSet: a 3D pose estimation dataset and baseline models for Cheetahs in the wild. 2021 IEEE International Conference on Robotics and Automation (ICRA) Preprint at https://arxiv.org/abs/2103.13282 (IEEE, Xi’an, China, 2021).143.Chen, Q. & Koltun, V. Photographic image synthesis with cascaded refinement networks. in 2019 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), 1511–1520 (2017).144.Lee, J., Hwangbo, J., Wellhausen, L., Koltun, V. & Hutter, M. Learning quadrupedal locomotion over challenging terrain. Sci. Robot. 5, eabc5986 (2020).145.Botella, C., Joly, A., Bonnet, P., Munoz, F. & Monestiez, P. Jointly estimating spatial sampling effort and habitat suitability for multiple species from opportunistic presence-only data. Methods Ecol. Evolution 12, 933–945 (2021).
Google Scholar 
146.Beery, S., Cole, E., Parker, J., Perona, P. & Winner, K. Species distribution modeling for machine learning practitioners: a review. in Proceedings of the 4th ACM SIGCAS Conference on Computing and Sustainable Societies (2021).147.Arzoumanian, Z., Holmberg, J. & Norman, B. An astronomical pattern-matching algorithm for computer-aided identification of whale sharks Rhincodon typus. J. Appl. Ecol. 42, 999–1011 (2005).
Google Scholar 
148.de Knegt, H. J., Eikelboom, J. A. J., van Langevelde, F., Spruyt, W. F. & Prins, H. H. T. Timely poacher detection and localization using sentinel animal movement. Sci. Rep. 11, 1–11 (2021).
Google Scholar 
149.Walter, T. & Couzin, I. D. TRex, a fast multi-animal tracking system with markerless identification, and 2D estimation of posture and visual fields. eLife 10, e64000 (2021).PubMed 
PubMed Central 

Google Scholar 
150.Kellenberger, B., Tuia, D. & Morris, D. AIDE: accelerating image-based ecological surveys with interactive machine learning. Methods Ecol. Evolution 11, 1716–1727 (2020).
Google Scholar 
151.Settles, B. Active learning. Synth. lectures Artif. Intell. Mach. Learn. 6, 1–114 (2012).MathSciNet 
MATH 

Google Scholar 
152.Ofli, F. et al. Combining human computing and machine learning to make sense of big (aerial) data for disaster response. Big Data 4, 47–59 (2016).PubMed 

Google Scholar 
153.Simpson, R., Page, K. R. & De Roure, D. Zooniverse: observing the world’s largest citizen science platform. in Proceedings of the 23rd International Conference on World Wide Web 1049–1054 (2014).154.Pocock, M. J. O., Roy, H. E., Preston, C. D. & Roy, D. B. The biological records centre: a pioneer of citizen science. Biol. J. Linn. Soc. 115, 475–493 (2015).
Google Scholar  More

The improved performance of the mangrove LUE model considering coastal environments in this study was mainly attributed to the determination of environmental scalars. Parameters determining environmental stressors (e.g., Topt, Tmin, Tmax, VPDmin, and VPDmax) were set based on the general characteristics of mangroves worldwide. It may not be as accurate for the mangroves in our study sites, but it generally reflects the response of mangroves to environmental changes. Furthermore, as can be seen in Fig. S1, it is applicable to our study sites. Despite the specific characteristics of each mangrove ecosystem at different sites being preferred, this study first offers the possibility to estimate mangrove productivity at a larger scale to track GPP, thus emphasizing the role of mangrove ecosystems nationally or worldwide.The validation results showed that the LUE values of the mangrove model agreed well with those estimated by EC method (Fig. 3) and indicated improved performance (slope = 0.8218–1.0108, intercept = -0.0006–0.0052, R2 = 0.54–0.64, RMSE = 0.0051–0.0068, Pearson’s r = 0.73–1), compared to the MOD17 model (slope = 0.4993–0.5566, intercept = 0.0311–0.0313, R2 = 0.24–0.45, RMSE = 0.0217–0.0220, Pearson’s r = 0.45–0.49). Firstly, the RS-based LUE model for terrestrial ecosystems (MOD17) considers only the environmental stressors of Tair and VPD. The photosynthesis in mangrove forests is influenced by other unique environmental factors caused by tidal inundation. According to Fig. S3, PAR caused the most significant effect on LUE, which is consistent with previous studies14,30,32. The impact of SST has not been quantitatively assessed, however, SST is a critical control that determines the upper latitudinal range of mangrove ecosystems12,33. In our study, the effects of SST and salinity on the mangrove LUE were quantified and helped improve LUE modeling.Secondly, LUEmax was typically defined for different land covers, however, there were no specific values for mangrove forests. In this study, the LUEmax of mangroves was first determined. It is worth noting that daytime NEE responses to PAR vary depending on the Tair23,30,34 so that LUEmax was determined separately at high, optimal, and low temperatures. The results showed that LUEmax reached a maximum when Tair was within the optimal range for mangroves, which represents the high productivity of mangrove ecosystems. Furthermore, the estimated LUEmax of mangrove forests (0.057) was larger than most terrestrial forests35,36,37, which could contribute to the high production and carbon sequestration in mangrove forests.Lastly, the relatively low stomatal conductance of mangroves leads to low LSP compared with terrestrial forests, which could result in the high-irradiance stress for photosynthesis38,39. Mangrove LSP ranges from about 0.2–1.2 mmol/m2/s, depending on the species and environments40,41,42. LUE was relatively low in April and May when seasonal PAR was high, as photosynthesis is more likely to reach saturation. Therefore, we assumed the LUE of mangroves decreased with increasing PAR. In addition, we found that the downscaling effect of PAR on LUE was not constant, but varied with increasing PAR. As follows, different PAR scalars were set for mangroves according to different PAR values. This is a first attempt at refining PARscalar considering different solar radiation, which represents a significant departure from the assumption of a constant downscaling effect of PAR in RS-driven models14,43. The accuracy of the LUE model was improved by refining the PARscalar with different downscaling slopes, especially in periods of high PAR values.Compared with the results obtained from flux-tower measurements, the modeled GPP was basically within the confidence interval of the measured results. The annual averages of GPP in Zhangjiang were 1729 g C/m2/year and 1924 g C/m2/year, in 2012 and 2016, and the annual mean value of GPP in Zhanjiang was 1434 g C/m2/year in 2015. The previous study showed that the GPP in Zhangjiang ranged from 1763 to 1919 g C/m2/year with a mean value of 1871 g C/m2/year32,44,45, which is in good agreement with the estimated values obtained in this study. Liu and Lai46 reported that the GPP of the Mai Po mangrove reserve was 2827 g C/m2/year. Rodda, et al.20 found a GPP value of 1271 g C/m2/year for Sunderbans mangroves in India. Gnanamoorthy, et al.47 estimated a GPP of 2305 g C/m2/year for Pichavaram mangroves. Variations in these estimates across sites were possibly caused by different climate-hydrological conditions, mangrove species, and ages. Differences in the same location may be due to different time scales and different methods of data gap filling and flux partitioning.In a similar way to the GPP model for terrestrial ecosystems48, the effect of the mangrove GPP model on the accuracy of GPP estimates can vary considerably under different environmental conditions. However, in comparison with the accuracy of models built for other vegetation types, the GPP model in this study performed substantially in two sites with RMSE of 2.54–3.41 g C/m2/day. Wang et al.49 adopted different models to estimate GPP for global vegetation and validation results showed the RMSE ranged from 1.79 to 2.33 g C/m2/day. Xiao, et al.50 demonstrated that the deviation between observed and predicted GPP was about 35–282 g C/m2 in an evergreen needleleaf forest. Also, the absolute GPP errors were 7.94–20.92% and 9.97–13.70% for maize cropland and degraded grassland36. Despite the discrepancy, our results were generally consistent with previous studies and were verified by field observations near the flux towers.The comparison of MODIS GPP and EC-estimated GPP showed that the MODIS GPP had a large fluctuation and weakly reflected productivity, being overestimated in 2012 and underestimated in 2015. Different meteorological inputs, different environmental scalars and fraction of absorbed photosynthetic active radiation (fAPAR) products in MODIS GPP and our mangrove GPP model can explain their different results. However, the improvements in our GPP model may help to obtain more accurate GPP estimates. The response of mangrove productivity to Tair has not been well-calibrated in the MODIS GPP product, which may partly account for the poor correlation between the MODIS GPP and EC estimates. Besides, MODIS GPP product was developed based on the International Geosphere-Biosphere Programme (IGBP) land cover map, which doesn’t include mangroves as a specific land cover37. Therefore, LUEmax and environmental parameters were not defined for mangroves, which varied with different environments. This may lead to uncertainty in MODIS GPP product for mangrove forests14. However, the GPP model generated in our study showed similar trends to the field measurements, capturing seasonal variations. The increase in the difference between MODIS GPP and EC estimates may be due to the assumption that the increase in GPP is linear with respect to PAR. In our model, the response of GPP to PAR was suppressed, resulting in seasonal changes in GPP that better match the observations. In addition, the GPP derived from this study was in higher agreement with measured values compared with GPP estimated from the vegetation photosynthesis model (VPM), as shown in Fig. S4. The improvement of this model was more obvious in winter (December to February), which may be due to the environmental stress of SST and PAR. The VPM without considering SSTscalar and PARscalar overestimated GPP in winter. It is indicated that the performance of the mangrove GPP model in this study varied with season. It is recommended to improve the estimation of GPP in the future by considering the seasonal variation of mangrove forests when determining environmental variables.Most studies provide EC-based estimates of GPP that are measurements from a limited footprint. It is possible to extrapolate results across similar vegetation types and geographic settings, but not to areas of heterogeneous vegetation. The RS-based GPP model offers spatial-scale estimates that can be directly incorporated into ecosystem-type models. PAR, SST, and salinity are the key environmental parameters of this RS-based mangrove GPP model. SST and salinity data were derived from the satellite images, while PAR was generated from the reconstructed PAR data, since it is more accurate than the existing RS data and has historical year data. However, PAR products from Hamawari-8, MERIS, and SeaWiFS are available now, which provide an opportunity to obtain large-scale PAR data using RS in the future. In addition to this, GPP of two mangrove forests was assessed and validated with three-year measurements. Validation at different sites and years showed similar results, which indicated the model has similar performance across mangrove forests. Nonetheless, these estimates need to be corroborated with EC databases, which are relatively accurate and provide many additional variables that are currently beyond the scope of higher spatial-resolution RS estimates. The proposed GPP model considering coastal environments was well suited to extend the study area by incorporating RS information and meteorological data. Currently, there are still few mangrove carbon flux towers worldwide. The LUE and GPP models proposed in this study are difficult to validate with measurements from flux towers in other countries. However, local measurements are available in many countries with large mangrove forests, such as Thailand, Vietnam, India, and Bangladesh. Therefore, it is expected that comparisons with measurements from previous studies can be conducted to show the consistency and applicability.The LUE model considering the effects of SST, salinity, and PAR performed well, however, the GPP estimated from the LUE, fAPAR, and PAR showed discrepancies and were generally lower than the measured values. Although the results are better than MODIS GPP products, limitations exist still.Firstly, the effects of salinity and SST on mangrove productivity were directly related to tidal activities. The soil pore water and surface water salinity could affect the osmotic pressure of mangroves especially for the submerged parts which would control the stomatal conductance. In the same way, SST could influence the temperature of mangrove root systems and soil sediments which has impacts on mangrove roots’ respiration and transpiration. Although, theoretically, salinity and SST should be considered as environmental variables affecting mangrove LUE, our results (Fig. S3) indicated that salinity and SST have little influence on mangrove productivity51. To date, the quantitative impact of SST has not been comprehensively unfolded, but it is a global control that determines the upper limit of the latitudinal range of mangroves12,33. The weak relationships between salinity, SST, and mangrove GPP could be due to the uncertainty caused by tidal inundation. Tide duration, tide height, and tide cycle would determine the effect of salinity and SST on the mangrove LUE and GPP. However, quantifying the influence from the tidal cycle remains a challenging task, which could result in the relatively poor performance of Salinityscalar and SSTscalar as shown in Fig. S3. Quantifying the soil temperature and surface water salinity considering the tidal cycle will contribute to model the LUE and GPP of mangrove forests.Secondly, mangroves of different species and ages exhibit diverse structural and physical conditions, resulting in different LUEmax, and optimal growing conditions such as Topt and VPDmin. The environmental settings would also vary from region to region. Liu and Lai46 found that LUE increased slightly with the increasing salinity below 15 ppt (R2 = 0.16). However, it was noted that photosynthetic activity of mangroves would be inhibited when the surface water salinity was high30,51,52,53. Probably, the mutual relationship between LUE and salinity depends on the salinity level and mangrove species. However, we have not specified the variables for different mangrove species, ages and locations which could be improved in the future. Besides, there are multicollinearities between different environmental variables. For example, Tair may have effects on SST and VPD, but as shown in Fig. S5, they are all important for mangrove photosynthesis. However, the correlations between them are not clear and need to be quantified in the future.Thirdly, the relatively low spatial and temporal resolution of the environmental data from RS would influence the accuracy of the model. The datasets have a relatively coarse resolution (usually 500 m–1 km and daily) and are thereby less suitable for smaller nature reserves, especially in the narrow patches of mangrove areas that are rapidly being exploited in coastal China. Moreover, the variability in LUE decreases with increasing temporal scale54. In our study, we determined the PARscalar based on the response of LUE to hourly-scale PAR and found the different down-regulation effects with increasing PAR. However, this phenomenon is not obvious in previous studies. Most RS-based LUE models were developed at a daily or 8-day temporal scale6,50,55,56,57. In terrestrial forests, the light saturated effect caused by increasing PAR was neglectable with coarse temporal scale because the average PAR was usually lower than the LSP. However, as the time scale increases, the effect of light saturation on LUE becomes more pronounced32,58,59. More importantly, this effect is more obvious in mangroves due to their lower LSP18,38, which makes it important in mangrove LUE modeling. The results in Fig. 3 show similar performances of LUE model on hourly and daily scale. Thus, we suggested that our model can be adopted in hourly and daily temporal resolution. However, the PARscalar developed in this study was based on the mangrove forests in one study site which may be influenced by the mangrove species with different LSP and light conditions. What’s more, VPD was on a monthly scale, which cannot reflect environmental dynamics. However, the hourly and daily VPD data are currently not available for coastal areas in China. Therefore, we used monthly averages to represent daily VPD, which may lead to uncertainty in the derived GPP estimates (Figs. 6 and 7). Besides, porewater salinity is controlled by sea surface salinity, precipitation, and river discharge. However, currently, pore water salinity was expressed in terms of sea surface salinity, which may lead to an underestimation of Salinityscalar. More systematic study is necessary to make it more applicable and accurate on a large scale, of which modeling the LUE for different mangrove species and locations is inevitable. However, serving as a fundamental and preliminary step, our study aims to provide a framework for RS-based mangrove GPP modeling. Recently, with the advancement of satellite imagery, hourly-scale RS data for PAR, temperature and SST are available. It can be expected that our current work could be further improved by investigating the light saturation effects in different mangrove forests and adopt higher temporal resolution RS products such as Himawari-8 and GCOM-C in the future.Lastly, the overall underestimation of GPP was mainly caused by the underestimation of fAPAR. Even though the fAPAR computed from Sentinel-2 had higher resolution and accuracy than MODIS fAPAR products, future improvements are still needed. Sentinel-2 fAPAR products (fAPAR-S2) was calculated as the instantaneous fAPAR obtained at 10:00 local solar time which only roughly represented the daily average but was not accurate. Besides, RS-derived fAPAR only considers the absorptions by living green vegetation elements, whereas the ground measured fAPAR refers to the contributions from all absorbing components60. The lower fAPAR-S2 values in mangrove forests may be due to the exposed-to-air root systems which absorb the radiation. Moreover, the spatial distribution of PAR was determined by Co-Kriging interpolation. The elevation was taken as the covariate to estimate spatial PAR. There are many other variables affecting the incoming PAR (e.g., slope and clearness)61. A more comprehensive set of variables needs to be included in the Co-kriging interpolation to improve the PAR estimation.The spatial and seasonal variations of the mangrove GPP were related to environmental changes along the shoreline. The low summer GPP was explained by the lower fAPAR in summer compared with other seasons, which was principally due to the underestimation of fAPAR in summer. Furthermore, PARscalar took a mean value of LSP as 1 mmol/m2/s, however, LSP varied with different species and environmental conditions. In summer, mangroves are more likely to obtain light saturation, and thus PARscalar may lead to an underestimation of LUE and thus GPP. On the contrary, PAR values in winter were relatively low but increased slightly with decreasing latitude. Thus, the inhibitory effect of PAR on LUE was not significant, and GPP increased with decreasing latitude. Salinity and VPD were more stable across years and locations and had no noticeable effect on the mangrove LUE and GPP. The seasonal latitudinal patterns and effects on mangrove productivity were similar for Tair and SST. Tair and SST were lower in winter, especially at high latitudes where mangroves were more sensitive to cold weather. Therefore, the GPP of mangroves at high latitudes in winter was the lowest throughout the year. However, hot weather in summer also limited the photosynthesis in mangroves, especially at low latitudes, where Tair and SST were higher. Nevertheless, there were some correlations among these environmental constants. For example, the Tair affects the vapor pressure and SST. There was a positive correlation between PAR and Tair. The multicollinearity among these variables and the various conditions of mangroves may affect the performance of the model and show variations along the coastline, which would be improved in future studies.Additionally, the GPP of mangroves increased from 2007 to 2018, which was mainly due to the expansion of mangrove forests in the coastal areas. As mangroves grow, canopy size and tree density increase, which may lead to higher LUE and less underestimation of fAPAR, thus contributing to high productivity. However, Zhejiang province (27° 02′ N–31° 11′ N) experienced extremely cold weather in January 2016 caused by the East Asia cold wave62,63, and large areas of mangrove forests died or became sick, leading to a decline in the mangrove GPP at high latitudes in 2018. More

1.Balmford, A. & Gaston, K. J. Why biodiversity surveys are good value How the ‘terror crocodile’ grew so big. Nature 398, 204–205 (1999).ADS 
CAS 
Article 

Google Scholar 
2.Rondinini, C., Marco, M. D., Visconti, P., Butchart, S. H. M. & Boitani, L. Update or outdate: Long-term viability of the IUCN Red List Making the Red List financially. Conserv. Lett. 7, 126–130 (2014).Article 

Google Scholar 
3.Braithwaite, M. E. & Walker, K. J. 50 Years of Mapping the British and Irish Flora 1962–2012 (Botanical Society of the British Isles, 2012).
Google Scholar 
4.Barrett, G., Silcocks, A., Barry, S., Cunningham, R. & Poulter, R. The New Atlas of Australian Birds (Royal Australasian Ornithologists Union, 2003).
Google Scholar 
5.Biodiversity Center of Japan. Animal Distribution Atlas of Japan (Ministry of Environment Japan, 2010).
Google Scholar 
6.McGowan, K. & Corwin, K. The Second Atlas of Breeding Birds in New York State (Comstock Publishing Associates, 2008).
Google Scholar 
7.Rhoads, A. F. & Klein, W. M. J. The vascular Flora of Pennsylvenia: Annoted Checklist and Atlas (American Philosophical Society, 1993).
Google Scholar 
8.Baker, H. et al. Population estimates of birds in Great Britain and the United Kingdom British Birds. Br. Birds 99, 25–44 (2006).
Google Scholar 
9.Bonn, A., Rodrigues, A. S. L. & Gaston, K. J. Threatened and endemic species: Are they good indicators of patterns of biodiversity on a national scale?. Ecol. Lett. 5, 733–741 (2002).Article 

Google Scholar 
10.Albuquerque, F. & Beier, P. Rarity-weighted richness: A simple and reliable alternative to integer programming and heuristic algorithms for minimum set and maximum coverage problems in conservation planning. PLoS ONE 10, e0119905 (2015).Article 

Google Scholar 
11.Williams, P. et al. A comparison of richness hotspots, rarity hotspots, and complementary areas for conserving diversity of British birds. Conserv. Biol. 10, 155–174 (1996).Article 

Google Scholar 
12.Platts, P. J. et al. Conservation implications of omitting narrow-ranging taxa from species distribution models, now and in the future. Divers. Distrib. 20, 1307-1320. (2014).Article 

Google Scholar 
13.Kujala, H., Moilanen, A. & Gordon, A. Spatial characteristics of species distributions as drivers in conservation prioritization. Methods Ecol. Evol. 9, 1121–1132 (2018).Article 

Google Scholar 
14.Kukkala, A. S. & Moilanen, A. Core concepts of spatial prioritisation in systematic conservation planning. Biol. Rev. Camb. Philos. Soc. 88, 443–464 (2013).Article 

Google Scholar 
15.Lawler, J. J., White, D., Sifneos, J. C. & Master, L. L. Rare species and the use of indicator groups for conservation planning. Conserv. Biol. 17, 875–882 (2003).Article 

Google Scholar 
16.Kujala, H., Lahoz-Monfort, J. J., Elith, J. & Moilanen, A. Not all data are equal: Influence of data type and amount in spatial conservation prioritisation. Methods Ecol. Evol. 9, 2249–2261 (2018).Article 

Google Scholar 
17.Grantham, H. S. et al. Diminishing return on investment for biodiversity data in conservation planning. Conserv. Lett. 1, 190–198 (2008).Article 

Google Scholar 
18.Margules, C. R. & Pressey, R. L. Systematic conservation planning. Nature 405, 243–253 (2000).CAS 
Article 

Google Scholar 
19.Moilanen, A., Wilson, K. A. & Possingham, H. Spatial Conservation Prioritization: Quantitative Methods and Computational Tools (Oxford University Press, 2009).
Google Scholar 
20.Akasaka, M., Kadoya, T., Ishihama, F., Fujita, T. & Fuller, R. A. Smart protected area placement decelerates biodiversity loss: A representation-extinction feedback leads rare species to extinction. Conserv. Lett. 10, 539–546 (2017).Article 

Google Scholar 
21.Boakes, E., McGowan, P. & Fuller, R. Distorted views of biodiversity: spatial and temporal bias in species occurrence data. PLoS Biol. 8, e1000385 (2010).Article 

Google Scholar 
22.Williams, P. H., Margules, C. R. & Hilbert, D. W. Data requirements and data sources for biodiversity priority area selection. J. Biosci. 27, 327–338 (2002).CAS 
Article 

Google Scholar 
23.Da Fonseca, G. A. B. et al. … following Africa’s lead Community groups could. Nature 405, 393–394 (2000).Article 

Google Scholar 
24.Possingham, H. P., Grantham, H. & Rondinini, C. How can you conserve species that haven’t been found? Commentary. J. Biogeogr. 34, 758–759 (2007).Article 

Google Scholar 
25.Ohlemuller, R. et al. The coincidence of climatic and species rarity: High risk to small-range species from climate change. Biol. Lett. 4, 568–572 (2008).Article 

Google Scholar 
26.Kier, G. et al. A global assessment of endemism and species richness across island and mainland regions. Proc. Natl. Acad. Sci. USA. 106, 9322–9327 (2009).ADS 
CAS 
Article 

Google Scholar 
27.Raes, N., Roos, M. C., Slik, J. W. F., Van Loon, E. E. & Steege, H. T. Botanical richness and endemicity patterns of Borneo derived from species distribution models. Ecography (Cop.) 32, 180–192 (2009).Article 

Google Scholar 
28.Tulloch, A. I. T., Mustin, K., Possingham, H. P., Szabo, J. K. & Wilson, K. A. To boldly go where no volunteer has gone before: predicting volunteer activity to prioritize surveys at the landscape scale. Divers. Distrib. 19, 465–480 (2013).Article 

Google Scholar 
29.Rodewald, A. D., Strimas-Mackey, M., Schuster, R. & Arcese, P. Tradeoffs in the value of biodiversity feature and cost data in conservation prioritization. Sci. Rep. 9, 5921 (2019).Article 

Google Scholar 
30.Rabinowitz, D. Seven forms of rarity. In The biological aspects of rare plant conservation. (ed. Synge, H.) 205–217 (Wiley, Chichester, 1981).
Google Scholar 
31.Ministry of Environment. Red Data Book 2014: Plants I (Gyousei, 2015).
Google Scholar 
32.Ministry of Environment. Red Data Book 2014: Plants II (Gyousei, 2015).
Google Scholar 
33.Yahara, T. et al. Red list of Japanese vascular plants: Summary of methods and results. Proc. Japan. Soc. Plant Taxon. 13, 89–96 (1998).
Google Scholar 
34.Ball, I. R., Possingham, H. P. & Watts, M. Marxan and relatives: Software for spatial conservation prioritisation. In Spatial Conservation Prioritisation: Quantitative Methods and Computational Tools (eds Moilanen, A. et al.) 185–195 (Oxford University Press, 2009).
Google Scholar 
35.Yoshioka, A., Akasaka, M. & Kadoya, T. Spatial prioritization for biodiversity restoration: A simple framework referencing past species distributions. Restor. Ecol. 22, 185–195 (2014).Article 

Google Scholar 
36.Naidoo, R. et al. Integrating economic costs into conservation planning. Trends Ecol. Evol. 21, 681–687 (2006).Article 

Google Scholar 
37.Japan Statistics. Population of Japan: Final Report of the 2005 Population Census. (Statistics Japan, 2010).38.R Core Team. R 4.0.0. (R Foundation for Statistical Computing, 2020).39.Zweig, M. H. & Campbell, G. Receiver-operating characteristic (ROC) plots: A fundamental evaluation tool in clinical medicine. Clin. Chem. 39, 561–577 (1993).CAS 
Article 

Google Scholar 
40.Hermoso, V., Ward, D. P. & Kennard, M. J. Using water residency time to enhance spatio-temporal connectivity for conservation planning in seasonally dynamic freshwater ecosystems. J. Appl. Ecol. 49, 1028–1035 (2012).Article 

Google Scholar 
41.Beyer, H. L., Dujardin, Y., Watts, M. E. & Possingham, H. P. Solving conservation planning problems with integer linear programming. Ecol. Modell. 328, 14–22 (2016).Article 

Google Scholar 
42.Schuster, R., Hanson, J. O., Strimas-Mackey, M. & Bennett, J. R. Integer linear programming outperforms simulated annealing for solving conservation planning problems. PeerJ 8, e9258. https://doi.org/10.7717/peerj.9258 (2020).Article 

Google Scholar  More

Evidence for biological and chemical isoprene consumption in coastal seawaterThe time course of isoprene concentration in coastal seawater samples incubated in closed glass bottles at the in situ temperature and in the dark demonstrated sustained loss for at least 45 h (Fig. 1a). Enclosure without headspace prevented isoprene loss by ventilation, and darkness was assumed to arrest all or most of the biological production25 and any photochemical production15 or degradation. Thus, the measured loss was considered the result of microbial degradation and chemical oxidation. In most cases an exponential function fitted better the decay than a linear function (Supplementary Table 1), indicating first-order (concentration-dependent) kinetics for isoprene loss.Fig. 1: Isoprene loss in dark incubations of coastal seawater.a Time course of isoprene concentration in 2 L dark incubations of non-filtered seawater samples from the back-reef lagoon of Mo’orea in April (blue) and the coastal Mediterranean in March (red) and May (green). Filled and open symbols correspond to duplicate incubations. Exponential fits to the data are shown by lines. See Supplementary Table 1 for fit equations and metrics, water temperatures and chlorophyll a concentrations. b Time course of isoprene concentration in series of 30 mL dark incubations of coastal Mediterranean seawater. Dark blue: non-filtered; red: filtered through 0.2 µm; green: filtered + 10 µmol L−1 H2O2; purple: filtered + 0.0025 units mL−1 bromoperoxidase (BrPO); light blue: filtered + H2O2 + BrPO. Exponential fit results in Supplementary Table 2.Full size imageIncubation of microorganism-devoid (filtered through 0.2 µm) coastal seawater sampled next to seaweeds showed an isoprene loss (0.12 d−1) that was half the loss in non-filtered water (0.20 d−1; Fig. 1b and Supplementary Table 2), implying that chemical oxidation accounted for half the total loss. Oxidation by OH·, the fastest amongst isoprene reactions with oxidative transients for which reaction rate data exist19, could account for the observed chemical loss. However, the possibility of oxidation by hitherto overlooked, pervasive oxidants like H2O2 deserved consideration. The addition of unrealistically high concentrations of either H2O2 or the enzyme bromoperoxidase (BrPO), substantially speeded up the chemical loss (0.91 d−1 with 10 µmol H2O2 L−1, 0.31 d−1 with 0.0025 units BrPO mL−1; Fig. 1b and Supplementary Table 2). Isoprene could have reacted with H2O2 in seawater as it does in acidic aerosols26. Besides, should dissolved27 BrPOs from seaweeds or outer-membrane-bound28 BrPOs from phytoplankton occur, they would have reacted with added H2O2 to produce hypobromous acid (HOBr), a strong oxidant29 that would further remove isoprene. Indeed, the addition of BrPO consumed isoprene because it produced HOBr by reaction with the naturally occurring H2O2. Confirming this interpretation, large HOBr production by simultaneous addition of BrPO and H2O2 caused complete isoprene removal in less than 4 h (Fig. 1b). Therefore, the results shown in Fig. 1b indicate that isoprene is reactive to pervasive H2O2 either directly or through the formation of enzymatically derived HOBr. All in all, first-order total isoprene loss (Fig. 1a) is expected to depend on photochemically-produced oxidants30 like H2O2, OH· and 1O2 as well as on microbiota through (a) microbial uptake and catabolism11 and (b) reaction with biologically produced oxidants26,31,32 like HOBr, H2O2 or superoxide.Variability of isoprene loss rate constants in the open oceanTen of the eleven offshore experimental sites were located in the open ocean, and one was located on the Southwestern Atlantic Shelf. Altogether they covered wide ranges of latitude (40°N–61°S), sea surface temperature (−0.8–28.6 °C), daily-averaged wind speed (3–12 m s−1), fluorometric chlorophyll-a (chla) concentration (0.1–5.8 mg m−3), and isoprene concentration (4–104 nmol m−3) (Fig. 2, Table 1 and Supplementary Table 3). Unfiltered seawater samples from the surface ocean were incubated in glass bottles for 24 h, at the in situ temperature and in the dark, and first-order loss rate constants were determined from initial and final isoprene concentrations (see Methods). Note that loss was determined under the assumption that isoprene production was arrested in the dark25. There is published evidence that residual isoprene production may occur in the dark33, but in our incubations, it was insufficient to counteract loss. Thus, isoprene losses caused by processes other than ventilation may have been underestimated.Fig. 2: Geographical distribution of the offshore experiments.Location of the sampling and incubation sites are shown by circles, coloured for isoprene concentration.Full size imageTable 1 Measured biological variables and isoprene process rate constants.Full size tableLoss rate constants (kloss = kbio + kchem) varied over an order of magnitude, ranging 0.03–0.64 d−1 with a median of 0.08 d−1 (Table 1). They did not show any significant relationship to sea surface temperature (SST) (Supplementary Fig. 1) but showed proportionality to the chla concentration (Fig. 3a) that was best described by the following linear regression equation:$${k}_{{{{{{rm{loss}}}}}}}=0.10; (pm 0.01),{{{{{rm{x}}}}}}, [{{{{{rm{chl}}}}}}a]+0.05; (pm 0.01)$$
(1)
Fig. 3: Isoprene processes and their main drivers.a Rate constant of isoprene loss in dark incubations (kloss, considered to be microbial and chemical consumption) vs. chlorophyll-a concentration. The linear regression equation is kloss = 0.10 × [chla] + 0.05 (R2 = 0.96, p = 10−7, n = 11). The standard error of the slope is 0.01 L mg−1 d−1, and the standard error of the intercept is 0.01 d−1. Error bars represent the experimentally determined standard error of kloss. The colour scale of the circles indicates bacterial abundances. b Specific (chla-normalised) rate of isoprene production vs seawater temperature (SST) across the sample series. The dashed line is the general smoothed trend. The blue line is the exponential adjustment at SST , 1000)$$
(2)
Substitution in Eq. (1) results in:$${k}_{{{{{{rm{loss}}}}}}}=0.14,{{{{{rm{x}}}}}}, {[{{{{{rm{chl}}}}}}{a}_{{{{{{rm{sat}}}}}}}]}^{1.28}+0.05$$
(3)
which is our recommended equation for kloss prediction from satellite chla. Note that only the variable term (kbio) changes from Eq. (1), while the intercept (kchem) is maintained at 0.05 d−1.Comparison of isoprene sinks and total turnover timeThe change of isoprene concentration ([iso]) in the surface mixed layer over time can be described as the budget of sources and sinks:$$varDelta [{{{{{rm{iso}}}}}}]/varDelta {{{{{rm{t}}}}}}=[{{{{{rm{iso}}}}}}]cdot ({k}_{{{{{{rm{prod}}}}}}} – {k}_{{{{{{rm{loss}}}}}}} – {k}_{{{{{{rm{vent}}}}}}} – {k}_{{{{{{rm{mix}}}}}}})$$
(4)
where kprod, kvent and kmix are the rate constants of isoprene production, ventilation to the atmosphere and vertical downward mixing by turbulent diffusion, respectively.We calculated kvent from our sampling sites over a period of 24 h (Table 1). Ventilation has been considered the main isoprene sink from the upper mixed layer of the ocean18. In our sampling sites, kloss was 0.4 to 10 times the kvent (median factor: 1.2). That is, loss through microbial + chemical consumption was of the same order as ventilation, sometimes considerably faster. Vertical mixing, kmix, was estimated to be one order of magnitude lower than the other process rates (Table 1), and in all cases but one it was calculated or assumed not to be a loss term but an import term into the mixed layer, because vertical profiles generally show maximum isoprene concentrations below the mixed layer and turbulent diffusion causes upward transport14,17. Altogether, the microbial, chemical, ventilation, and, where relevant, mixing losses resulted in total turnover times (1/(kloss + kvent + kmix)) of isoprene between 1.4 and 16 days, median 5 days (Table 1).Isoprene productionAssuming steady-state for isoprene concentrations over 24 h (Supplementary Fig. 2), i.e. Δ[iso]/Δt = 0 in Eq. (4), the sum of the daily rate constants of all sinks (kloss + kvent) equals the rate constant of isoprene production (kprod), with kmix adding to either side depending on whether it is an import to or an export from the mixed layer (Table 1). Note that kprod was the highest coinciding with higher [chla]. This is consistent with a recent study44 where measurement of the net biological isoprene production (i.e. production — consumption rates) across seasons in the open ocean was attempted; net production rates increased in May, coinciding with a large increase in [chla] and phytoplankton cell abundance.The product of kprod by the isoprene concentration gives the daily isoprene production rate, which can be normalised by dividing it by the chla concentration. In our study, this specific isoprene production rate varied between 1 and 38 nmol (mg chla)−1 d−1 (Table 1), median 8 nmol (mg chla)−1 d−1. These values are within the broad range reported across phytoplankton taxa from laboratory studies with monocultures41,45 (0.3–32, median 3 nmol (mg chla)−1 d−1, n = 124). Five of the eleven sites gave values >13 nmol (mg chla)−1 d−1, i.e. in the higher end of the laboratory data range. This is not unexpected, since measurements in monoculture experiments are typically conducted before reaching nutrient limitation, below light saturation and in the absence of UV radiation, to mention three stressors commonly occurring in the surface open ocean. If isoprene biosynthesis and release is enhanced by any of these stressors, as is the case in vascular plants7,10, then monoculture-derived results will easily render underestimates of isoprene production in the open ocean. Production by heterotrophic bacteria46 could have also contributed to increase apparent specific isoprene production rates, but the occurrence and importance of this process in the marine environment is unknown.When plotted against the SST, which was also the temperature of the incubations, specific isoprene production rates increased exponentially between −0.8 and 23 °C and dropped drastically at higher SST (Fig. 3b). Several studies with phytoplankton monocultures have reported positive dependence of specific isoprene production rates on temperature45,47,48,49,50. One of these studies45 described that the increase with temperature reaches an optimum for production that varies among phytoplankton strains and with light intensity, but falls around 23–26 °C. The most detailed study47 was conducted with a Prochlorococcus strain; remarkably, the shape of the specific production rate vs. temperature curve for this cyanobacterium strain was almost identical to that of Fig. 3b, with an exponential increase until 23 °C and a drop thereafter. This is the canonical curve type of enzymatic activities, but the thermal behaviour of the enzymes for isoprene synthesis in marine unicellular algae has not yet been characterised12.Revising the magnitude and players of the marine isoprene cycleOur results allow redrawing the isoprene cycle in the surface mixed layer of the ocean. Figure 4 sketches the magnitude of the rate constants for production and sinks presented in Table 1, averaged according to a chla concentration threshold: the blue and green arrows correspond to the experiments in waters with [chla] lower and higher than 0.4 mg m−3, respectively. Isoprene production in productive (chla-richer) waters is faster than in oligotrophic (chla-poorer) waters. Vertical mixing is assumed to majorly constitute an input into the mixed layer, yet very small. Photochemical production and emission from surfactants15 in the surface microlayer of productive waters is depicted as uncertain. Among sinks, the microbiota-dependent consumption is much faster in productive waters; actually, the statistical uncertainty of Eq. (1) and the uneven distribution of incubation results along the [chla] axis hamper resolving kbio in phytoplankton-poor waters ( More

1.Cowen, R. K. & Sponaugle, S. Larval dispersal and marine population connectivity. Ann. Rev. Mar. Sci. 1, 443–466 (2009).PubMed 

Google Scholar 
2.Hellberg, M. E. Gene flow and isolation among populations of marine animals. Annu. Rev. Ecol. Evol. Syst. 40, 291–310 (2009).
Google Scholar 
3.Selkoe, K. A. et al. Taking the chaos out of genetic patchiness: Seascape genetics reveals ecological and oceanographic drivers of genetic patterns in three temperate reef species. Mol. Ecol. 19, 3708–3726 (2010).PubMed 

Google Scholar 
4.Guo, X. et al. Phylogeography of the rock shell Thais clavigera (Mollusca): Evidence for long-distance dispersal in the Northwestern Pacific. PLoS ONE 10, e0129715 (2015).PubMed 
PubMed Central 

Google Scholar 
5.Hoffman, J. I., Clarke, A., Linse, K. & Peck, L. S. Effects of brooding and broadcasting reproductive modes on the population genetic structure of two Antarctic gastropod molluscs. Mar. Biol. 158, 287–296 (2011).
Google Scholar 
6.Modica, M. V., Russini, V., Fassio, G. & Oliverio, M. Do larval types affect genetic connectivity at sea? Testing hypothesis in two sibling marine gastropods with contrasting larval development. Mar. Environ. Res. 127, 92–101 (2017).CAS 
PubMed 

Google Scholar 
7.Je Lee, H. & Boulding, E. G. Spatial and temporal population genetic structure of four northeastern Pacific littorinid gastropods: The effect of mode of larval development on variation at one mitochondrial and two nuclear DNA markers. Mol. Ecol. 18, 2165–2184 (2009).
Google Scholar 
8.Barbosa, S. S., Klanten, S. O., Puritz, J. B., Toonen, R. J. & Byrne, M. Very fine-scale population genetic structure of sympatric asterinid sea stars with benthic and pelagic larvae: Influence of mating system and dispersal potential. Biol. J. Linn. Soc. 108, 821–833 (2013).
Google Scholar 
9.Shanks, A. L. Pelagic larval duration and dispersal distance revisited. Biol. Bull. 216, 373–385 (2009).PubMed 

Google Scholar 
10.Riginos, C., Buckley, Y. M., Blomberg, S. P. & Treml, E. A. Dispersal capacity predicts both population genetic structure and species richness in reef fishes. Am. Nat. 184, 52–64 (2014).PubMed 

Google Scholar 
11.Wort, E. J. G. et al. Contrasting genetic structure of sympatric congeneric gastropods: Do differences in habitat preference, abundance and distribution matter?. J. Biogeogr. 46, 369–380 (2019).
Google Scholar 
12.Ayre, D. J., Minchinton, T. E. & Perrin, C. Does life history predict past and current connectivity for rocky intertidal invertebrates across a marine biogeographic barrier?. Mol. Ecol. 18, 1887–1903 (2009).CAS 
PubMed 

Google Scholar 
13.Meyer, C. P., Geller, J. B. & Paulay, G. Fine scale endemism on coral reefs: Archipelagic differentiation in turbinid gastropods. Evolution (N. Y.) 59, 113–125 (2005).
Google Scholar 
14.White, C. et al. Ocean currents help explain population genetic structure. Proc. R. Soc. B Biol. Sci. 277, 1685–1694 (2010).
Google Scholar 
15.Marko, P. B. ‘What’s larvae got to do with it?’ Disparate patterns of post-glacial population structure in two benthic marine gastropods with identical dispersal potential. Mol. Ecol. 13, 597–611 (2004).CAS 
PubMed 

Google Scholar 
16.Edmands, S. Phylogeography of the intertidal copepod Tigriopus californicus reveals substantially reduced population differentiation at northern latitudes. Mol. Ecol. 10, 1743–1750 (2001).CAS 
PubMed 

Google Scholar 
17.Ni, G., Li, Q., Kong, L. & Yu, H. Comparative phylogeography in marginal seas of the northwestern Pacific. Mol. Ecol. 23, 534–548 (2014).PubMed 

Google Scholar 
18.Vendrami, D. L. J. et al. RAD sequencing sheds new light on the genetic structure and local adaptation of European scallops and resolves their demographic histories. Sci. Rep. 9, 1–13 (2019).CAS 

Google Scholar 
19.Sandoval-Castillo, J., Robinson, N. A., Hart, A. M., Strain, L. W. S. & Beheregaray, L. B. Seascape genomics reveals adaptive divergence in a connected and commercially important mollusc, the greenlip abalone (Haliotis laevigata), along a longitudinal environmental gradient. Mol. Ecol. 27, 1603–1620 (2018).PubMed 

Google Scholar 
20.Hirai, J. Insights into reproductive isolation within the pelagic copepod Pleuromamma abdominalis with high genetic diversity using genome-wide SNP data. Mar. Biol. 167, 1–6 (2020).CAS 

Google Scholar 
21.Hosoya, S. et al. Random PCR-based genotyping by sequencing technology GRAS-Di (genotyping by random amplicon sequencing, direct) reveals genetic structure of mangrove fishes. Mol. Ecol. Resour. 19, 1153–1163 (2019).CAS 
PubMed 

Google Scholar 
22.Losos, J. B. & Ricklefs, R. E. Adaptation and diversification on islands. Nature https://doi.org/10.1038/nature07893 (2009).Article 
PubMed 

Google Scholar 
23.Savolainen, V. et al. Sympatric speciation in palms on an oceanic island. Nature https://doi.org/10.1038/nature04566 (2006).Article 
PubMed 

Google Scholar 
24.Parent, C. E. & Crespi, B. J. Ecological opportunity in adaptive radiation of Galápagos endemic land snails. Am. Nat. https://doi.org/10.1086/646604 (2009).Article 
PubMed 

Google Scholar 
25.Chiba, S. & Cowie, R. H. Evolution and extinction of land snails on oceanic islands. Annu. Rev. Ecol. Evol. Syst. 47, 123–141 (2016).
Google Scholar 
26.Grant, P. R. & Grant, B. R. Unpredictable evolution in a 30-year study of Darwin’s finches. Science (80-.) 296, 707–711 (2002).CAS 
ADS 

Google Scholar 
27.Scheltema, R. The relevance of passive dispersal for the biogeography of Caribbean mollusks. Am. Malacol. Bull. 11, 95–115 (1995).
Google Scholar 
28.Bernardi, G. et al. Darwin’s fishes: Phylogeography of Galápagos Islands reef fishes. Bull. Mar. Sci. 90, 533–549 (2014).
Google Scholar 
29.Eble, J. A., Toonen, R. J. & Bowen, B. W. Endemism and dispersal: Comparative phylogeography of three surgeonfishes across the Hawaiian Archipelago. Mar. Biol. 156, 689–698 (2009).
Google Scholar 
30.Tomokuni, M. M. Aquatic and Semiaquatic Insects of the Bonin Islands (including the Volcano Islands). Mem. Natl. Sci. Museum (1978).31.Sugawara, T., Watanabe, K., Kato, H. & Yasuda, K. Dioecy in Wikstroemia pseudoretusa (Thymelaeaceae) endemic to the Bonin (Ogasawara) islands. APG Acta Phytotaxon. Geobot. https://doi.org/10.18942/apg.KJ00004622804 (2004).Article 

Google Scholar 
32.Chiba, S. Species diversity and conservation of Mandarina, an endemic land snail of the Ogasawara Islands. In Restoring the Oceanic Island Ecosystem: Impact and Management of Invasive Alien Species in the Bonin Islands (eds Kawakami, K. & Okochi, I.) 117–125 (Springer, 2010). https://doi.org/10.1007/978-4-431-53859-2_18.Chapter 

Google Scholar 
33.Mukai, T., Nakamura, S., Suzuki, T. & Nishida, M. Mitochondrial DNA divergence in yoshinobori gobies (Rhinogobius species complex) between the Bonin Islands and the Japan-Ryukyu Archipelago. Ichthyol. Res. 52, 410–413 (2005).
Google Scholar 
34.Shih, H. T., Komai, T. & Liu, M. Y. A new species of fiddler crab from the Ogasawara (Bonin) Islands, Japan, separated from the widely-distributed sister species Uca (Paraleptuca) crassipes (White, 1847) (Crustacea: Decapoda: Brachyura: Ocypodidae). Zootaxa 3746, 175–193 (2013).PubMed 

Google Scholar 
35.Yamazaki, D. et al. Genetic diversification of intertidal gastropoda in an archipelago: The effects of islands, oceanic currents, and ecology. Mar. Biol. https://doi.org/10.1007/s00227-017-3207-9 (2017).Article 

Google Scholar 
36.Nakano, T., Takahashi, K. & Ozawa, T. Description of an endangered new species of Lunella (Gastropoda:Turbinidae) from the Ogasawara Islands, Japan. Venus J. Malacol. Soc. Japan 66, 1–10 (2007).
Google Scholar 
37.Nakano, T., Yazaki, I., Kurokawa, M., Yamaguchi, K. & Kuwasawa, K. The origin of the endemic patellogastropod limpets of the Ogasawara Islands in the northwestern Pacific. J. Molluscan Stud. 75, 87–90 (2009).
Google Scholar 
38.González-Wevar, C. A., Nakano, T., Palma, A. & Poulin, E. Biogeography in cellana (patellogastropoda, nacellidae) with special emphasis on the relationships of southern hemisphere oceanic island species. PLoS ONE 12, 1–16 (2017).
Google Scholar 
39.Tenggardjaja, K. A., Bowen, B. W. & Bernardi, G. Reef fish dispersal in the Hawaiian Archipelago: Comparative phylogeography of three endemic damselfishes. J. Mar. Sci. https://doi.org/10.1155/2016/3251814 (2016).Article 

Google Scholar 
40.Tenggardjaja, K. A., Bowen, B. W. & Bernardi, G. Comparative phylogeography of widespread and endemic damselfishes in the Hawaiian Archipelago. Mar. Biol. 165, 1–21 (2018).
Google Scholar 
41.Kurozumi, T. & Asakura, A. Marine molluscs from the northern Mariana Islands, Micronesia. Nat. Hist. Res. Spec. Issue 1, 121–168 (1994).
Google Scholar 
42.Nakano, D. & Makoto, N. Age structure and growth in a population of Monodonta labio (Linnaeus) at Shima Peninsula, Japan. Venus J. Malacol. Soc. Japan 40, 34–40 (1981).
Google Scholar 
43.Hashino, T. & Tomiyama, K. Life history of Monodonta labio confusa Tapprone-Canefri, 1874 in Kagoshima Bay, Kyushu, Japan and age estimation based on annual ring analysis of shell. Nat. Kagoshima 39, 143–155 (2013).
Google Scholar 
44.Yoh, A. & Sakurai, I. Reproductive cycle and food habits of the herbivorous snail Monodonta confusa off the coast of Suttsu Bay in southwestern Hokkaido, Japan. Proc. Sch. Biol. Sci. Tokai Univ. 6, 17–23 (2017).
Google Scholar 
45.Sasaki, R. Larval identification and occurrence of ezo abalone, Haliotis discus hannai, in the adjacent waters of Kesennuma Bay, Miyagi Prefecture. Suisan Zoushoku 32, 199–206 (1985).
Google Scholar 
46.Yamazaki, D., Miura, O., Uchida, S., Ikeda, M. & Chiba, S. Comparative seascape genetics of co-distributed intertidal snails Monodonta spp. in the Japanese and Ryukyu archipelagoes. Mar. Ecol. Prog. Ser. 657, 135–146 (2020).ADS 

Google Scholar 
47.Ballard, J. W. O. & Whitlock, M. C. The incomplete natural history of mitochondria. Mol. Ecol. 13, 729–744 (2004).PubMed 

Google Scholar 
48.Parham, J. F. et al. Genetic introgression and hybridization in Antillean freshwater turtles (Trachemys) revealed by coalescent analyses of mitochondrial and cloned nuclear markers. Mol. Phylogenet. Evol. 67, 176–187 (2013).CAS 
PubMed 

Google Scholar 
49.Hirano, T. et al. Enigmatic incongruence between mtDNA and nDNA revealed by multi-locus phylogenomic analyses in freshwater snails. Sci. Rep. 9, 6223 (2019).PubMed 
PubMed Central 
ADS 

Google Scholar 
50.Funk, D. J. & Omland, K. E. Species-level paraphyly and polyphyly: Frequency, causes, and consequences, with insights from animal mitochondrial DNA. Annu. Rev. Ecol. Evol. Syst. 34, 397–423 (2003).
Google Scholar 
51.Toews, D. P. L. & Brelsford, A. The biogeography of mitochondrial and nuclear discordance in animals. Mol. Ecol. 21, 3907–3930 (2012).CAS 
PubMed 

Google Scholar 
52.Hirase, S. et al. Integrative genomic phylogeography reveals signs of mitonuclear incompatibility in a natural hybrid goby population. Evolution (N.Y.) 75, 176–194 (2021).
Google Scholar 
53.Zhao, D., Li, Q., Kong, L. & Yu, H. Cryptic diversity of marine gastropod Monodonta labio (Trochidae): Did the early Pleistocene glacial isolation and sea surface temperature gradient jointly drive diversification of sister species and/or subspecies in the Northwestern Pacific?. Mar. Ecol. https://doi.org/10.1111/maec.12443 (2017).Article 

Google Scholar 
54.Mukai, T., Nakamura, S. & Nishida, M. Genetic population structure of a reef goby, Bathygobius cocosensis, in the northwestern Pacific. Ichthyol. Res. 56, 380–387 (2009).
Google Scholar 
55.Keith, S. A., Herbert, R. J. H., Norton, P. A., Hawkins, S. J. & Newton, A. C. Individualistic species limitations of climate-induced range expansions generated by meso-scale dispersal barriers. Divers. Distrib. 17, 275–286 (2011).
Google Scholar 
56.Faurby, S. & Barber, P. H. Theoretical limits to the correlation between pelagic larval duration and population genetic structure. Mol. Ecol. 21, 3419–3432 (2012).PubMed 

Google Scholar 
57.Funk, W. C. et al. Adaptive divergence despite strong genetic drift: Genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis). Mol. Ecol. 25, 2176–2194 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
58.Whiteley, A. R. et al. Genetic variation and effective population size in isolated populations of coastal cutthroat trout. Conserv. Genet. 11, 1929–1943 (2010).
Google Scholar 
59.Riginos, C., Douglas, K. E., Jin, Y., Shanahan, D. F. & Treml, E. A. Effects of geography and life history traits on genetic differentiation in benthic marine fishes. Ecography (Cop.) 34, 566–575 (2011).
Google Scholar 
60.Kuriiwa, K., Chiba, S. N., Motomura, H. & Matsuura, K. Phylogeography of Blacktip Grouper, Epinephelus fasciatus (Perciformes: Serranidae), and influence of the Kuroshio Current on cryptic lineages and genetic population structure. Ichthyol. Res. 61, 361–374 (2014).
Google Scholar 
61.Tachikawa, H. Nature profile of the isolated oceanic island, the Bonin Islands. Midoriishi 5, 27–29 (1994).
Google Scholar 
62.Setsuko, S. et al. Genetic diversity, structure, and demography of Pandanus boninensis (Pandanaceae) with sea drifted seeds, endemic to the Ogasawara Islands of Japan: Comparison between young and old islands. Mol. Ecol. 29, 1050–1068 (2020).PubMed 

Google Scholar 
63.Asakura, A. & Nishihama, S. Studies on the biology and ecology of the intertidal animals of Chichijima Island in the Ogasawara (Bonin) Islands III: Description, form and habitat of the trochid snail, Monodonta perplexa boninensis n. subsp. in comparison with those in Monodonta perpl. Venus J. Malacol. Soc. Japan 46, 194–201 (1987).
Google Scholar 
64.Nakano, T. & Minato, R. Marine organisms in the intertidal zone of Nishinoshima Island. Ogasawara Res. 46, 109–121 (2019).
Google Scholar 
65.Sasaki, T. & Horikoshi, K. Marine animals of Minami-Iw-To lsland. Ogasawara Res. 33, 155–171 (2008).
Google Scholar 
66.Williams, S., Apte, D., Ozawa, T., Kaligis, F. & Nakano, T. Speciation and dispersal along continental coastlines and island arcs in the indo-west pacific turbinid gastropod genus lunella. Evolution (N. Y.) 65, 1752–1771 (2011).
Google Scholar 
67.Siddall, M. et al. Sea-level fluctuations during the last glacial cycle. Nature 423, 853–858 (2003).CAS 
PubMed 
ADS 

Google Scholar 
68.Setsuko, S. et al. Genetic variation of pantropical Terminalia catappa plants with sea-drifted seeds in the Bonin Islands: Suggestions for transplantation guidelines. Plant Species Biol. 32, 13–24 (2017).
Google Scholar 
69.Hedgecock, D. Is gene flow from pelagic larval dispersal important in the adaptation and evolution of marine invertebrates?. Bull. Mar. Sci. 39, 550–564 (1986).
Google Scholar 
70.Parsons, K. E. The genetic effects of larval dispersal depend on spatial scale and habitat characteristics. Mar. Biol. 126, 403–414 (1996).CAS 

Google Scholar 
71.Pechenik, J. A. On the advantages and disadvantages of larval stages in benthic marine invertebrate life cycles. Mar. Ecol. Prog. Ser. 177, 269–297 (1999).ADS 

Google Scholar 
72.Scheltema, R. S. Larval dispersal as a means of genetic exchange between geographically separated populations of shallow-water benthic marine gastropods. Biol. Bull. 140, 284–322 (1971).
Google Scholar 
73.Wright, L. I., Tregenza, T. & Hosken, D. J. Inbreeding, inbreeding depression and extinction. Conserv. Genet. 9, 833–843 (2008).
Google Scholar 
74.Caley, M. J. et al. Recruitment and the local dynamics of open marine populations. Annu. Rev. Ecol. Syst. 27, 477–500 (1996).
Google Scholar 
75.Johannesson, K. The paradox of Rockall: Why is a brooding gastropod (Littorina saxatilis) more widespread than one having a planktonic larval dispersal stage (L. littorea)?. Mar. Biol. 99, 507–513 (1988).
Google Scholar 
76.Nakajima, Y., Nishikawa, A., Iguchi, A. & Sakai, K. Regional genetic differentiation among northern high-latitude island populations of a broadcast-spawning coral. Coral Reefs 31, 1125–1133 (2012).ADS 

Google Scholar 
77.Bowen, B. W. et al. Comparative phylogeography of the ocean planet. Proc. Natl. Acad. Sci. U. S. A. 113, 7962–7969 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
78.Funk, W. C., McKay, J. K., Hohenlohe, P. A. & Allendorf, F. W. Harnessing genomics for delineating conservation units. Trends Ecol. Evol. 27, 489–496 (2012).PubMed 
PubMed Central 

Google Scholar 
79.Palumbi, S. R. Population genetics, demographic connectivity, and the design of marine reserves. Ecol. Appl. 13, 146–158 (2003).
Google Scholar 
80.Jones, G., Srinivasan, M. & Almany, G. Population connectivity and conservation of marine biodiversity. Oceanography 20, 100–111 (2007).
Google Scholar 
81.Colgan, D. J., Ponder, W. F., Beacham, E. & Macaranas, J. M. Gastropod phylogeny based on six segments from four genes representing coding or non-coding and mitochondrial or nuclear DNA. Molluscan Res. https://doi.org/10.1071/MR03002 (2003).Article 

Google Scholar 
82.Griekspoor, A. & Groothuis, T. 4peaks. Ver. 1.7.1. http://nucleobytes.com/4peaks/ (2005).83.Thompson, J. D., Higgins, D. G. & Gibson, T. J. CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. https://doi.org/10.1093/nar/22.22.4673 (1994).Article 
PubMed 
PubMed Central 

Google Scholar 
84.Edgar, R. C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. https://doi.org/10.1093/nar/gkh340 (2004).Article 
PubMed 
PubMed Central 

Google Scholar 
85.Excoffier, L. & Lischer, H. E. L. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567 (2010).PubMed 
PubMed Central 

Google Scholar 
86.Excoffier, L., Smouse, P. E. & Quattro, J. M. Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics 131, 479–491 (1992).CAS 
PubMed 
PubMed Central 

Google Scholar 
87.Bandelt, H., Forster, P. & Röhl, A. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16, 37–48 (1999).CAS 
PubMed 
PubMed Central 

Google Scholar 
88.Leigh, J. W. & Bryant, D. POPART: Full-feature software for haplotype network construction. Methods Ecol. Evol. 6, 1110–1116 (2015).
Google Scholar 
89.Peterson, B. K., Weber, J. N., Kay, E. H., Fisher, H. S. & Hoekstra, H. E. Double digest RADseq: An inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS ONE 7, e37135 (2012).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
90.Eaton, D. A. R. & Overcast, I. ipyrad: Interactive assembly and analysis of RADseq datasets. Bioinformatics https://doi.org/10.1093/bioinformatics/btz966 (2020).Article 
PubMed 

Google Scholar 
91.Meirmans, P. G. & Van Tienderen, P. H. genotype and genodive: Two programs for the analysis of genetic diversity of asexual organisms. Mol. Ecol. Notes 4, 792–794 (2004).
Google Scholar 
92.Hammer, Ø., Harper, D. A. T. & Ryan, P. D. PAST: Paleontological statistics software package for education and data analysis. Palaeontol. Electron. 4, 9 (2001).
Google Scholar 
93.Mussmann, S. M., Douglas, M. R., Chafin, T. K. & Douglas, M. E. BA3-SNPs: Contemporary migration reconfigured in BayesAss for next-generation sequence data. Methods Ecol. Evol. 10, 1808–1813 (2019).
Google Scholar 
94.Rambaut, A. & Drummond, A. J. Tracer v1.6. http://tree.bio.ed.ac.uk/software/tracer/ (2013).95.Excoffier, L., Dupanloup, I., Huerta-Sánchez, E., Sousa, V. C. & Foll, M. Robust demographic inference from genomic and SNP data. PLoS Genet. 9, e1003905 (2013).PubMed 
PubMed Central 

Google Scholar  More

1.Falkowski, P. G., Fenchel, T. & Delong, E. F. The microbial engines that drive earth’s biogeochemical cycles. Science 320, 1034–1039 (2008).CAS 
ADS 

Google Scholar 
2.Fuhrman, J. A., Cram, J. A. & Needham, D. M. Marine microbial community dynamics and their ecological interpretation. Nat. Rev. Microbiol. 13, 133–146 (2015).CAS 
PubMed 

Google Scholar 
3.Giovannoni, S. J. & Stingl, U. Molecular diversity and ecology of microbial plankton. Nature 437, 343–348 (2005).CAS 
PubMed 
ADS 

Google Scholar 
4.Kujawinski, E. B. The impact of microbial metabolism on marine dissolved organic matter. Ann. Rev. Mar. Sci. 3, 567–599 (2011).PubMed 

Google Scholar 
5.Moran, M. A. The global ocean microbiome. Science 350, aac8455 (2015).PubMed 

Google Scholar 
6.Pedrós-Alió, C. The rare bacterial biosphere. Ann. Rev. Mar. Sci. 4, 15.1-15.18 (2012).
Google Scholar 
7.Salazar, G. et al. Global diversity and biogeography of deep-sea pelagic prokaryotes. ISME J. 10, 596–608 (2015).PubMed 
PubMed Central 

Google Scholar 
8.Sogin, M. L. et al. Microbial diversity in the deep sea and the underexplored ‘rare biosphere’. Proc. Natl. Acad. Sci. 103, 12115–12120 (2006).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
9.Kirchman, D. L. Growth rates of microbes in the oceans. Ann. Rev. Mar. Sci. 8, 285–309 (2016).PubMed 

Google Scholar 
10.Campbell, B. J., Yu, L., Heidelberg, J. F. & Kirchman, D. L. Activity of abundant and rare bacteria in a coastal ocean. Proc. Natl. Acad. Sci. 108, 12776–12781 (2011).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
11.Salter, I. et al. Seasonal dynamics of active SAR11 ecotypes in the oligotrophic Northwest Mediterranean Sea. ISME J. 9, 347–360 (2015).CAS 
PubMed 

Google Scholar 
12.Giovannoni, S. J., Cameron Thrash, J. & Temperton, B. Implications of streamlining theory for microbial ecology. ISME J. 8, 1553–1565 (2014).PubMed 
PubMed Central 

Google Scholar 
13.Våge, S., Storesund, J. E. & Thingstad, T. F. Adding a cost of resistance description extends the ability of virus-host model to explain observed patterns in structure and function of pelagic microbial communities. Environ. Microbiol. 15, 1842–1852 (2012).
Google Scholar 
14.Våge, S., Storesund, J. E. & Thingstad, T. F. SAR11 viruses and defensive host strains. Nature 499, 9–11 (2013).
Google Scholar 
15.Giovannoni, S., Temperton, B. & Zhao, Y. Giovannoni et al. reply. Nature 499, 9–11 (2013).
Google Scholar 
16.Zhao, Y. et al. Abundant SAR11 viruses in the ocean. Nature 494, 357–360 (2013).CAS 
PubMed 
ADS 

Google Scholar 
17.Herndl, G. J. et al. Contribution of Archaea to total prokaryotic production in the deep Atlantic Ocean. Appl. Environ. Microbiol. 71, 2303–2309 (2005).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
18.Teira, E., Lebaron, P., Van Aken, H. & Herndl, G. J. Distribution and activity of Bacteria and Archaea in the deep water masses of the North Atlantic. Limnol. Oceanogr. 51, 2131–2144 (2006).CAS 
ADS 

Google Scholar 
19.Newton, R. J. & Shade, A. Lifestyles of rarity: Understanding heterotrophic strategies to inform the ecology of the microbial rare biosphere. Aquat. Microb. Ecol. 78, 51–63 (2016).
Google Scholar 
20.Hamasaki, K., Taniguchi, A., Tada, Y., Long, R. A. & Azam, F. Actively growing bacteria in the Inland Sea of Japan, identified by combined bromodeoxyuridine immunocapture and denaturing gradient gel electrophoresis. Appl. Environ. Microbiol. 73, 2787–2798 (2007).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
21.Tada, Y., Makabe, R., Kasamatsu-Takazawa, N., Taniguchi, A. & Hamasaki, K. Growth and distribution patterns of Roseobacter/Rhodobacter, SAR11, and Bacteroidetes lineages in the Southern Ocean. Polar Biol. 36, 691–704 (2013).
Google Scholar 
22.Suttle, C. A. Marine viruses—Major players in the global ecosystem. Nat. Rev. Microbiol. 5, 801–812 (2007).CAS 
PubMed 

Google Scholar 
23.Pernthaler, J. Predation on prokaryotes in the water column and its ecological implications. Nat. Rev. Microbiol. 3, 537–546 (2005).CAS 
PubMed 

Google Scholar 
24.Mena, C. et al. Seasonal niche partitioning of surface temperate open ocean prokaryotic communities. Front. Microbiol. 11, 1749 (2020).PubMed 
PubMed Central 

Google Scholar 
25.Mena, C. et al. Dynamic prokaryotic communities in the dark western Mediterranean Sea. Sci. Rep. 11, 17859 (2021).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
26.Urbach, E., Vergin, K. L. & Giovannoni, S. J. Immunochemical detection and isolation of DNA from metabolically active bacteria. Appl. Environ. Microbiol. 65, 1207–1213 (1999).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
27.Hatzenpichler, R. et al. In situ visualization of newly synthesized proteins in environmental microbes using amino acid tagging and click chemistry. Environ. Microbiol. 16, 2568–2590 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
28.Emerson, J. et al. Schrödinger’s microbes: Tools for distinguishing the living from the dead in microbial ecosystems. Micobiome 5, 86 (2017).
Google Scholar 
29.Smriga, S., Samo, T., Malfatti, F., Villareal, J. & Azam, F. Individual cell DNA synthesis within natural marine bacterial assemblages as detected by ‘click’ chemistry. Aquat. Microb. Ecol. 72, 269–280 (2014).
Google Scholar 
30.Reichart, N. et al. Activity-based cell sorting reveals responses of uncultured archaea and bacteria to substrate amendment. ISME J. 14, 2851–2861 (2020).CAS 
PubMed 
PubMed Central 

Google Scholar 
31.Bakenhus, I. et al. Composition of total and cell-proliferating bacterioplankton community in early summer in the North Sea—Roseobacters are the most active component. Front. Microbiol. 8, 1–14 (2017).
Google Scholar 
32.Morris, R. M. et al. SAR11 clade dominates ocean surface bacterioplankton communities. Nature 420, 806–810 (2002).CAS 
PubMed 
ADS 

Google Scholar 
33.Giovannoni, S. J. SAR11 Bacteria: The most abundant plankton in the oceans. Ann. Rev. Mar. Sci. 9, 231–255 (2017).PubMed 

Google Scholar 
34.Clifford, E. L. et al. Taurine is a major carbon and energy source for marine prokaryotes in the North Atlantic Ocean off the Iberian Peninsula. Microb. Ecol. 78, 299–312 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
35.Tripp, H. J. et al. SAR11 marine bacteria require exogenous reduced sulphur for growth. Nature 452, 741–744 (2008).CAS 
PubMed 
ADS 

Google Scholar 
36.Carlson, C. A. et al. Seasonal dynamics of SAR11 populations in the euphotic and mesopelagic zones of the northwestern Sargasso Sea. ISME J. 3, 283–295 (2009).CAS 
PubMed 

Google Scholar 
37.Winter, C., Bouvier, T., Weinbauer, M. G. & Thingstad, T. F. Trade-offs between competition and defense specialists among unicellular planktonic organisms: The ‘Killing the Winner’ hypothesis revisited. Microbiol. Mol. Biol. Rev. 74, 42–57 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
38.Vergin, K. L. et al. High-resolution SAR11 ecotype dynamics at the Bermuda Atlantic Time-series Study site by phylogenetic placement of pyrosequences. ISME J. 7, 1322–1332 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
39.Hugoni, M. et al. Structure of the rare archaeal biosphere and seasonal dynamics of active ecotypes in surface coastal waters. Proc. Natl. Acad. Sci. 110, 1–6 (2013).
Google Scholar 
40.Qin, W. et al. Marine ammonia-oxidizing archaeal isolates display obligate mixotrophy and wide ecotypic variation. Proc. Natl. Acad. Sci. 111, 12504–12509 (2014).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
41.Pollard, P. C. & Moriarty, D. J. W. Validity of the tritiated thymidine method for estimating bacterial growth rates: Measurement of isotope dilution during DNA synthesis. Appl. Environ. Microbiol. 48, 1076–1083 (1984).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
42.Wellsbury, P., Herbert, R. A. & John Parkes, R. Incorporation of [methyl-3H]thymidine by obligate and facultative anaerobic bacteria when grown under defined culture conditions. FEMS Microbiol. Ecol. 12, 87–95 (1993).CAS 

Google Scholar 
43.Clausen, A., Matakos, A., Sandrini, M. & Piskur, J. Thymidine kinases in Archaea. Nucleosides Nucleotides Nucleic Acids 25, 1159–1163 (2006).CAS 
PubMed 

Google Scholar 
44.Hamasaki, K., Long, R. A. & Azam, F. Individual cell growth rates of marine bacteria, measured by bromodeoxyuridine incorporation. Aquat. Microb. Ecol. 35, 217–227 (2004).
Google Scholar 
45.Qin, W. et al. Influence of oxygen availability on the activities of ammonia-oxidizing Archaea. Environ. Microbiol. Rep. 9, 250–256 (2017).CAS 
PubMed 

Google Scholar 
46.Reji, L., Tolar, B. B., Smith, J. M., Chavez, F. P. & Francis, C. A. Differential co-occurrence relationships shaping ecotype diversification within Thaumarchaeota populations in the coastal ocean water column. ISME J. 13, 1144–1158 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
47.Sebastián, M. et al. Deep ocean prokaryotic communities are remarkably malleable when facing long-term starvation. Environ. Microbiol. 20, 713–723 (2018).PubMed 

Google Scholar 
48.Vergin, K. L., Done, B., Carlson, C. A. & Giovannoni, S. J. Spatiotemporal distributions of rare bacterioplankton populations indicate adaptive strategies in the oligotrophic ocean. Aquat. Microb. Ecol. 71, 1–13 (2013).
Google Scholar 
49.Tada, Y., Taniguchi, A., Sato-Takabe, Y. & Hamasaki, K. Growth and succession patterns of major phylogenetic groups of marine bacteria during a mesocosm diatom bloom. J. Oceanogr. 68, 509–519 (2012).
Google Scholar 
50.Mestre, M. et al. Sinking particles promote vertical connectivity in the ocean microbiome. Proc. Natl. Acad. Sci. 115, E6799–E6807 (2018).CAS 
PubMed 
PubMed Central 

Google Scholar 
51.Ruiz-González, C. et al. Major imprint of surface plankton on deep ocean prokaryotic structure and activity. Mol. Ecol. 29, 1820–1838 (2020).PubMed 

Google Scholar 
52.Chen, X., Ma, R., Yang, Y., Jiao, N. & Zhang, R. Viral regulation on bacterial community impacted by lysis-lysogeny switch: A microcosm experiment in eutrophic coastal waters. Front. Microbiol. 10, 1763 (2019).PubMed 
PubMed Central 

Google Scholar 
53.McCarren, J. et al. Microbial community transcriptomes reveal microbes and metabolic pathways associated with dissolved organic matter turnover in the sea. Proc. Natl. Acad. Sci. 107, 16420–16427 (2010).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
54.Reintjes, G., Arnosti, C., Fuchs, B. & Amann, R. Selfish, sharing and scavenging bacteria in the Atlantic Ocean: A biogeographical study of bacterial substrate utilisation. ISME J. 13, 1119–1132 (2019).CAS 
PubMed 

Google Scholar 
55.Middelboe, M. Bacterial growth rate and marine virus-host dynamics. Microb. Ecol. 40, 114–124 (2000).CAS 
PubMed 

Google Scholar 
56.Buchan, A., Lecleir, G. R., Gulvik, C. A. & González, J. M. Master recyclers: Features and functions of bacteria associated with phytoplankton blooms. Nat. Rev. Microbiol. 12, 686–698 (2014).CAS 
PubMed 

Google Scholar 
57.Mou, X. et al. Bromodeoxyuridine labelling and fluorescence-activated cell sorting of polyamine-transforming bacterioplankton in coastal seawater. Environ. Microbiol. 17, 876–888 (2014).PubMed 

Google Scholar 
58.Azam, F. & Malfatti, F. Microbial structuring of marine ecosystems. Nat. Rev. Microbiol. 5, 782–791 (2007).CAS 
PubMed 

Google Scholar 
59.Yilmaz, P., Yarza, P., Rapp, J. Z. & Glöckner, F. O. Expanding the world of marine bacterial and archaeal clades. Front. Microbiol. 6, 1524 (2016).PubMed 
PubMed Central 

Google Scholar 
60.Coe, A. et al. Survival of Prochlorococcus in extended darkness. Limnol. Oceanogr. 61, 1375–1388 (2016).ADS 

Google Scholar 
61.Cottrell, M. T. & Kirchman, D. L. Photoheterotrophic microbes in the arctic ocean in summer and winter. Appl. Environ. Microbiol. 75, 4958–4966 (2009).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
62.Zeder, M., Peter, S., Shabarova, T. & Pernthaler, J. A small population of planktonic Flavobacteria with disproportionally high growth during the spring phytoplankton bloom in a prealpine lake. Environ. Microbiol. 11, 2676–2686 (2009).PubMed 

Google Scholar 
63.Cottrell, M. T. & Kirchman, D. L. Natural assemblages of marine proteobacteria and members of the Cytophaga-flavobacter cluster consuming low- and high-molecular-weight dissolved organic matter. Appl. Environ. Microbiol. 66, 1692–1697 (2000).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
64.Banning, E. C., Casciotti, K. L. & Kujawinski, E. B. Novel strains isolated from a coastal aquifer suggest a predatory role for flavobacteria. FEMS Microbiol. Ecol. 73, 254–270 (2010).CAS 
PubMed 

Google Scholar 
65.Uchimiya, M. et al. Coupled response of bacterial production to a wind-induced fall phytoplankton bloom and sediment resuspension in the chukchi sea shelf, Western Arctic Ocean. Front. Mar. Sci. 3, 1–12 (2016).
Google Scholar 
66.Ivancic, I. et al. Seasonal variations in extracellular enzymatic activity in marine snow-associated microbial communities and their impact on the surrounding water. FEMS Microbiol. Ecol. 94, fiy198 (2018).CAS 

Google Scholar 
67.Cram, J. A. et al. Seasonal and interannual variability of the marine bacterioplankton community throughout the water column over ten years. ISME J. 9, 563–580 (2015).PubMed 

Google Scholar 
68.Manca, B. et al. Physical and biochemical averaged vertical profiles in the Mediterranean regions: An important tool to trace the climatology of water masses and to validate incoming data from operational oceanography. J. Mar. Syst. 48, 83–116 (2004).
Google Scholar 
69.Puig, P. & Palanques, A. Temporal variability and composition of settling particle fluxes on the Barcelona continental margin (Northwestern Mediterranean). J. Mar. Res. 56, 639–654 (1998).
Google Scholar 
70.Buesseler, K. O. & Boyd, P. W. Shedding light on processes that control particle export and flux attenuation in the twilight zone of the open ocean. Limnol. Oceanogr. 54, 1210–1232 (2009).CAS 
ADS 

Google Scholar 
71.Alonso-González, I. J., Arístegui, J., Lee, C. & Calafat, A. Regional and temporal variability of sinking organic matter in the subtropical northeast Atlantic Ocean: A biomarker diagnosis. Biogeosciences 7, 2101–2115 (2010).ADS 

Google Scholar 
72.Hunt, D. E. et al. Relationship between abundance and specific activity of bacterioplankton in open ocean surface waters. Appl. Environ. Microbiol. 79, 177–184 (2013).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
73.Campbell, B. J. & Kirchman, D. L. Bacterial diversity, community structure and potential growth rates along an estuarine salinity gradient. ISME J. 7, 210–220 (2013).CAS 
PubMed 

Google Scholar 
74.Alderkamp, A. C., Sintes, E. & Herndl, G. J. Abundance and activity of major groups of prokaryotic plankton in the coastal North Sea during spring and summer. Aquat. Microb. Ecol. 45, 237–246 (2006).
Google Scholar 
75.De Corte, D., Sintes, E., Yokokawa, T. & Herndl, G. J. Comparison between MICRO-CARD-FISH and 16S rRNA gene clone libraries to assess the active versus total bacterial community in the coastal Arctic. Environ. Microbiol. Rep. 5, 272–281 (2013).PubMed 

Google Scholar 
76.Bergauer, K. et al. Organic matter processing by microbial communities throughout the Atlantic water column as revealed by metaproteomics. Proc. Natl. Acad. Sci. 115, E400–E408 (2018).CAS 
PubMed 

Google Scholar 
77.Georges, A. A., El-Swais, H., Craig, S. E., Li, W. K. W. & Walsh, D. A. Metaproteomic analysis of a winter to spring succession in coastal northwest Atlantic Ocean microbial plankton. ISME J. 8, 1301–1313 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
78.Couradeau, E. et al. Probing the active fraction of soil microbiomes using BONCAT-FACS. Nat. Commun. 10, 2770 (2019).PubMed 
PubMed Central 
ADS 

Google Scholar 
79.Wu, X. et al. Culturing of ‘unculturable’ subsurface microbes: Natural organic carbon source fuels the growth of diverse and distinct bacteria from groundwater. Front. Microbiol. 11, 610001 (2020).PubMed 
PubMed Central 

Google Scholar 
80.Alonso-Sáez, L., Díaz-Pérez, L. & Morán, X. A. G. The hidden seasonality of the rare biosphere in coastal marine bacterioplankton. Environ. Microbiol. 17, 3766–3780 (2015).PubMed 

Google Scholar 
81.Liu, J., Meng, Z., Liu, X. & Zhang, X.-H. Microbial assembly, interaction, functioning, activity and diversification: A review derived from community compositional data. Mar. Life Sci. Technol. 1, 112–128 (2019).ADS 

Google Scholar 
82.Long, R. A. & Azam, F. Antagonistic interactions among marine pelagic bacteria. Appl. Environ. Microbiol. 67, 4975–4983 (2001).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
83.López-Jurado, J. L. et al. The RADMED monitoring programme as a tool for MSFD implementation: Towards an ecosystem-based approach. Ocean Sci. 11, 897–908 (2015).ADS 

Google Scholar 
84.Strickland, J. D. H. & Parsons, T. R. A Practical Handbook of Seawater Analysis (Fisheries Research Board of Canada, 1968).
Google Scholar 
85.Grasshoff, K., Ehrhardt, M. & Kremling, K. Methods of seawater analysis (Verlag Chemie GmbH, 1983). https://doi.org/10.1002/iroh.19850700232.Book 

Google Scholar 
86.Murphy, J. & Riley, J. P. A modified single solution method for the determination of phosphate in natural waters. Anal. Chim. Acta 27, 31–36 (1962).CAS 

Google Scholar 
87.Brussaard, C. P. D. Optimization of procedures for counting viruses by flow cytometry. Appl. Environ. Microbiol. 70, 1506–1513 (2004).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
88.Gasol, J. M. & del Giorgio, P. A. Using flow cytometry for counting natural planktonic bacteria and understanding the structure of planktonic bacterial communities. Sci. Mar. 64, 197–224 (2000).
Google Scholar 
89.Parada, A. E., Needham, D. M. & Fuhrman, J. A. Every base matters: Assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples. Environ. Microbiol. 18, 1403–1414 (2016).CAS 
PubMed 

Google Scholar 
90.Bolyen, E. et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME2. Nat. Biotechnol. 37, 852–857 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
91.Callahan, B. J. et al. Dada2: High-resolution sample inference from illumina amplicon data. Nat. Methods 13, 581–583 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
92.Katoh, K. & Standley, D. M. MAFFT Multiple sequence aligment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 30, 772–780 (2013).CAS 
PubMed 
PubMed Central 

Google Scholar 
93.Eren, A. M. et al. Oligotyping: Differentiating between closely related microbial taxa using 16S rRNA gene data. Methods Ecol. Evol. 4, 1111–1119 (2013).PubMed Central 

Google Scholar  More