Ecology

Anna Lena Bercht interviewed fishers in Lofoten, Norway, to assess how climate change was affecting their livelihoods.Credit: Anna Lena Bercht

I remember February 2011, when, in the Chinese megacity of Guangzhou, an older man finally overcame his scepticism about being interviewed and invited me to sit down next to him on a stone bench under a shady tree. I held my notebook on my lap, and we sat on either side of a translator and talked about his life and world for more than two hours. It was one of the most informative and revealing interviews that I had done during my fieldwork in the city.
Making it in the megacity
One of the most fundamental challenges in qualitative fieldwork is gaining access to research participants. This is often time-consuming and labour-intensive, particularly when the topic requires in-depth methods and addresses a sensitive subject.Advice that goes beyond the usual recommendations of establishing relationships with gatekeepers, ensuring anonymity for interviewees and relying on the snowball sampling technique (in which one research participant suggests further ones) is rare. In this light, I’m happy to share some simple, but often neglected, examples from my qualitative fieldwork in the lively Guangzhou (where I worked for 12 months)1 and on the remote, Arctic island chain of Lofoten, Norway (done over 4 months)2, that might offer some inspiration and encouragement.I have a background in human geography, and did my PhD on experiences of stress, coping and resilience among the Chinese population of Guangzhou in the face of the city’s rapid urbanization. I travelled there five times to help to establish research cooperation with Chinese scholars, make field observations, select a case-study site and interview locals. I, together with other PhD students, stayed in a typical Chinese high-rise apartment in a neighbourhood that wasn’t a common choice for expatriates. Living side-by-side with the locals gave us a perfect opportunity to experience genuine everyday life and Chinese culture.My first postdoctoral project after my PhD brought me to Lofoten, where I looked at psychological barriers to climate adaptation in small-scale coastal fisheries. I went to Lofoten twice. On my first visit, I travelled across the whole archipelago by bus for one month to get a profound overview of the fishing villages and local living conditions, and to conduct first interviews. During my second visit, I stayed for a total of three months in rental locations near fishing harbours, and conducted more extensive interviews.In both China and Norway, I used in-depth interviews to learn about the challenges that people face. I asked people about unemployment, about the possibility of being forced to move elsewhere and about how climate change might affect their livelihoods. This required a sensitive and thoughtful approach to ‘getting invited’ into people’s lives. In Guangzhou, German- and English-speaking Chinese students assisted me as translators (and interpreters, when needed). On Lofoten, I conducted the interviews myself in English.There are two ways to access research participants: physical access, which refers to the ability of the researcher to get in direct face-to-face contact with people, and mental access. Successful mental access means that interlocutors open up about why they think, feel and behave as they do. Physical access is a necessary condition for mental access; however, in my experience, both are equally valuable.

Chinese interviewees in Guangzou shared their feelings about the rapid urbanization of their city.Credit: Anna Lena Bercht

Compared with Lofoten, it took longer to get physical access to local inhabitants in China. Presumably, this was because of the language barrier and reliance on translators, as well as cultural differences. Trust is considered a central tenet in Chinese relationships, and time and effort are needed to let it grow. During my time in Guangzhou, I occasionally benefited from being a foreigner: people were touched that someone from abroad showed genuine interest in their well-being. In Lofoten, fishers appreciated talking to a social scientist instead of a natural scientist who would have mainly asked questions about fishing quotas and catch volume.My advice for other social scientists hoping to gain access to research participants falls into those two categories.How to get good physical accessUse local public transport. Using local public transport creates many unexpected opportunities to bump into people, get into conversations and gain relevant information. For example, while waiting at a bus stop in Lofoten, I came across an art-gallery owner from a fishing village. He wondered why I was travelling out of the peak tourism season. I ended up with an invitation to his gallery, where he introduced me to two retired fishers whom he had also invited. Without the gallerist and his proactive networking, I probably would not have been given the chance to interview these two very informative and engaging fishers.In a metro station in Guangzhou, a toddler kept staring at me and tried to touch my light hair. This small interaction led me to chat to the toddler’s father, who recommended that I talk to a local teacher to learn more about the area’s history. His advice opened up important insights into urban-restructuring processes that I would have missed otherwise.
Nine ‘brain food’ tips for researchers
Use local media. In Norway, a journalist was at the harbour to get first-hand information on the year’s cod catch, when he saw me interviewing fishers. He became curious and eager to learn more about my work. In the end, he wrote an article about my research, which was published a few days later across Lofoten. His article was a door-opener for me.People recognized me from my photo in the article and contacted me to tell me about their lives and the cod fisheries. They also invited me on their vessels and put me in touch with other key informants.Change your workplace. During fieldwork, a workplace is often needed for interview transcription, literature research and interim data analysis. Moving the workplace outside wherever you are staying during a field trip allows you to immerse yourself in the daily lives of local people and interact with them more easily. For me, such agile ‘mini-office’ locations were cafes, public libraries and picnic tables. In this way, I was able to recruit interview partners on the spot.How to create deeper mental accessWear appropriate outfits. First impressions count, always. Researchers are judged not only on what they say and how they say it, but also on how they look. Certain clothes, such as those with a political slogan or religious symbol, have certain meanings and connotations. Depending on the context and whom you talk to, your appearance could promote or impede making connections and building rapport. For instance, whereas my practical ‘outdoorsy’ get-dirty outfit was appropriate for interviews on fishing vessels, a modest appearance (non-branded clothes and a simple style) was useful in rural areas of Guangzhou.Show respect. Just like in any other relationship, respect and humility play a crucial part in building a trustworthy interviewer–interviewee relationship. Showing respect can be subtly embedded in conversations in many ways, including in the content of questions and the manner in which they are asked. When interviewees started to close down when asked about painful issues, such as underemployment or loss of identity, I upheld their privacy, comfort and security by not probing when given an evasive answer. Instead, I changed the interview focus and, when appropriate, cautiously reapproached the sensitive issue by using interview techniques such as roleplaying. Interviewees were asked to put themselves in the position of someone else, such as a spatial planner or politician, and assess the issue at hand from this perspective. Taking such an imaginary role can help to make the interviewees feel more secure and face pain more openly.Be humble. Having a modest view of yourself is essential to communicate at eye level with people. As a scientist, you can easily fall into the trap of thinking that your thoughts and concepts are somehow more valuable because you are well-educated and established. However, you are the one asking questions — and the interviewees, whether they are fishers, farmers or homeless people, often know more about many things than you do. Being aware of this is an expression of humility. I let the interviewees know that they were the local experts and I was the foreign learner.Use small talk. Small talk — including non-verbal communication, such as smiling, or connective gestures, for example handing out a handkerchief or offering some tea — has an essential bonding function. Talking about ‘safe’ topics can help the interviewee to overcome the feelings of otherness, newness and discomfort that can emerge in an interview, and fosters social cohesiveness. This can help to counteract the asymmetrical power relationship between the researcher (who asks) and the researched (who answers). For example, before substantive questioning, I created shared experiences by talking about last night’s storm or the world cod-fishing championship, which takes place every year in Lofoten. This took the relationship to a greater level of intimacy and togetherness — which small talk after finishing the interview can strengthen. I remember joking about my stamina for eating properly with chopsticks to one interviewee.Use self-disclosure. Revealing selected information about yourself and sharing your own thoughts with interlocutors can help to create and reaffirm a sphere of confidentiality and trust. Fishers in Norway would, for instance, often ask “What interested you in Lofoten coastal fisheries?” or “Why do you ask me and not the scientists from Tromsø University?” I answered such questions honestly, which assisted in creating a more balanced relationship, encouraging the interviewees to address sensitive subjects more openly and readily.Change interview sites. In several interviews, I found that the answers given tended to depend on where the interview was held and which identity that site evoked for the interviewee. For example, a fisher did not talk about climate-change concerns on his fishing vessel (any concern was masked by his existential fear of losing his livelihood as a coastal fisher), but he later that day freely discussed his worries in his home. Changing the interview site can be a helpful technique to access hidden thoughts and feelings.Above all, be realistic. You will probably make mistakes; I regretted not dressing warmly enough on a fishing vessel in Arctic weather. Locals will find you amusing, weird or impolite. They will keep out of your way, and you will never know why. And they will terminate interviews prematurely with no excuse. And that’s all right. In the end, fieldwork is a combination of planning, resources, time, skills, hard work, commitment, headache, joy — and luck. Learn from your mistakes, and accept the things you cannot change. More

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain
the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in
Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles
and JavaScript. More

Methods for data analysis in figuresAll analyses in figures were performed using Mathematica 12.3 (Wolfram Research, Inc., Champaign, IL, USA).Analysis in Fig. 1
C&D. To calculate integrated abundances of E. huxleyi cells and EhV, we first selected days for which all the bags had a non-null value. Values were then summed up to obtain the integrated abundance.E&J. We computed a standard linear fit between the E. huxleyi total abundances and total EhV abundances for covered and uncovered bags separately. We followed the same procedure for the correlations in panel J and provide a comparison between different models in Supplementary Fig. 5.Analysis in Fig. 2
A. The ASVs that were selected appeared at a relative abundance of at least 2% in at least 4 samples for the 0.2–2 µm 16S sequences and at least in 8 samples for the 2–20 µm 18S sequences. Abundances were concatenated for each time point and normalized by row, to have maximum relative abundance of 1 across all samples. ASVs were sorted by the position of their individual center of mass ({t}_{{CM}}) defined by$${t}_{{CM}}=,frac{mathop{sum}limits_{i}{t}_{i}f({t}_{i})}{mathop{sum}limits_{i}f({t}_{i})}$$
(1)
with i representing the different time points and f(({t}_{i})) the relative abundance of the ASV. The same figure for the individual bags in shown in Supplementary Fig. 14 and Supplementary Fig. 15.B. We selected 18S ASVs with a maximum relative abundance of at least 2% and observed in at least five samples. We averaged relative abundance across bags and then smoothed the time series with a moving average filter (width 2). Then, we grouped all ASVs into clusters based on their cosine distance using Mathematica’s FindClusters function and the KMeans method. The number of possible clusters ranged from 2 to 12, and the final number of clusters was decided using the silhouette method71. Only silhouette scores for 2 and 6 clusters were positive (between-cluster distance minus within-cluster distance).D. We subset reads that map to either Flavobacteriales or Fhodobacterales, then renormalized within each class, taking the mean over bags. Results per bag are shown in Supplementary Fig. 9.F. The turnover time was defined by the exponential rate k at which the Bray-Curtis similarity ({BC}(t)) declined over time. To this end, for a given bag, we computed the Bray–Curtis similarity between the composition vector at a starting day t’ with all following days t, giving a curve that declined roughly exponentially. For earlier starting days (for which the similarity curves declined the furthest), we found that the Bray–Curtis similarity never reached 0 but instead leveled out around ({{BC}}_{infty }=0.05) (due to ASVs that are constantly present in all the samples and maintain a minimal level of similarity between bags). Thus, we imposed an offset at(,{{BC}}_{infty }) for all fits (using Mathematica’s FindFit function) with the function:$${BC}(t)=(1-{{BC}}_{{{infty }}}) times {e}^{-kleft({t}^{{prime} }-tright)}+{{BC}}_{{{infty }}}$$
(2)
The turnover is averaged over bags, showing the standard deviation as error bars in the figure.G. To find differentially abundant ASVs, we first selected a subset of ASVs that had a maximum abundance of at least 10%, and performed Mann–Whitney U-Tests between the relative abundance values of a given ASV in the focal bag and all the other bags over all timepoints of the bloom’s demise. Correcting for multiple testing, we found four 16S ASVs that were differentially abundant in any of the bags, three of which were specific to bag 7, shown in Fig. 2g; and five 18S ASVs, two specific to bags 5 and 6 (Rhizosolenia delicatula and Aplanochytrium), one specific to bag 4 (Pterosperma), and two specific to bag 7 (MAST-1C and Woloszynskia halophila, shown in Fig. 2g).H. The divergence between bags was calculated as follows: we first measured, for each bag, the Bray–Curtis distance between this given bag and all the other bags at the end of the experiment (Supplementary Fig. 13). In order to control for the existing differences between bags at the beginning of the bloom, Bray–Curtis distances were normalized according to the differences between bags at the starting day of the E. huxleyi bloom. As the exact starting days of the bloom is not clear, we normalized for starting days 11, 12, or 13. The plot shows averages with the standard deviation as error bars. For the 18S microbiome, we first removed reads that map to E. huxleyi to reduce bias toward bag 7 (which had by far the lowest E. huxleyi abundance, Fig. 1c).Analysis in Fig. 3
A. Functional annotation of dominant 18S ASVs was based on manual literature search for the 100 most abundant 18S ASVs. Automatic annotation using the functional database created by72 gave qualitatively identical results but contained fewer organisms (covering about 50% of reads). The relative abundance of each trait was obtained by summing up the relative abundance of all the species harboring a specific trait. We used the annotations from72 to further subdivide heterotrophs into osmotrophs, saprotrophs, and other types of heterotrophy (e.g., grazing), ignoring ASVs with missing annotations.D. Growth rates were computed by fitting a linear model to the log-transformed absolute abundances. For thraustochytrids, we measured growth rates until the abundances reached their maximum, i.e., for days indicated by solid lines in Fig. 3b. For bacteria in the 0.2–2 micron fraction, we measured growth rates during the bloom and demise of E. huxleyi, i.e., for the time period after day 15 until the final day, except for bag 4 (until day 22) and bag 7 (until day 18) to account for their different bloom and demise dynamics. For bacteria in the 2–20 micron fraction, we measured growth rates similarly, starting after day 10 until the final day, except bags 4 and 7 (until day 22).E. To quantify the rate of change k of the biomass ratio of thraustochytrids to bacteria we fit a linear function to the log of biomass ratio from day 10 to the time point t where the ratio was maximal; for bag 7, this was day 18, for all others, day 23. We thus have:$$,{{log }},{BR},(t)={kt},+,{{log }},{BR},(0)$$
(3)
Analysis in Fig. 4
C&D. Since TEP accumulates over time, it cannot be expressed as a weighted sum of phytoplankton abundances. Instead, we formulate the model as a recursive relation where TEP can be produced by E. huxleyi, naked nanophytoplankton, and picophytoplankton, and degraded or lost through sinking:$${TEP}left(tright)=left(1-dright){TEP}left(t-1right)+{a}_{E}Eleft(tright)+{a}_{N}Nleft(tright)+{a}_{P}Pleft(tright),$$
(4)
The amount of TEP at time t is given by the fraction (1-d) of TEP at time t-1, where d corresponds to the fraction of TEP that is degraded between time points, plus the amount of TEP produced by the phytoplankton cells present at time t (or time t-1, which gives equivalent results). E, N, and P correspond to E. huxleyi, naked nanophytoplankton, and picophytoplankton, respectively. The parameter ({a}_{E}) corresponds to the amount of TEP produced per E. huxleyi cell, reported in panel D. ({a}_{E}) is set to be fixed through time, and different for each bag. This recursion can be solved to give an explicit expression for TEP(t):$${TEP}left(tright)=mathop{sum }limits_{{t}^{{prime} }=0}^{t}{left(1-dright)}^{t-{t}^{{prime} }}[{a}_{E}Eleft({t}^{{prime} }right)+{a}_{N}Nleft({t}^{{prime} }right)+{a}_{P}Pleft({t}^{{prime} }right)].$$
(5)
This functional form was then used to perform a linear model fitting with the constraint ({a}_{i}ge 0) for various values of the parameter d. The best fit, defined by maximum ({R}^{2}) over the resulting linear model, was used to fix d = 0.12. Our model considers that the fraction of non-calcified E. huxleyi cells in the nanophytoplankton counts is small.Larger phytoplankton cells ( >40 μm) filtered out from flow-cytometry measurements can also be a major source of TEP, despite low cell density. In order to verify this, FlowCam data was analyzed. None of the identified classes of larger phytoplankton (such as Phaeocystis or Dinobryon) increased in a systematic manner toward later stages of the bloom, explaining why larger phytoplankton were not included in the TEP model (Supplementary Fig. 24 and Supplementary Fig. 25).E. Using the smFISH method that reports the proportion of infected E. huxleyi cells, we estimated the amount of TEP produced from infected cells. We first used the least infected uncovered bags (bags 1 and 3) as a baseline to fix model parameters such as how much TEP does a non-infected cell produce. We then split the E. huxleyi abundance into an uninfected subpopulation producing T TEP/cell as in the uninfected bags, and an infected subpopulation producing I×T TEP/cells. To define I, we combined the fixed model parameters (i.e., amount of TEP produced per cell from Fig. 4d for bags 1 and 3) with the measured fraction of infected cells. We adjusted the factor I = 4 to minimize deviation of the measure total TEP concentration from the model prediction including the two subpopulations. The same procedure was used for panel H, using the corresponding model for PIC.F&G. To model the amount of PIC produced per cell we assume that the measured PIC only increases via new E. huxleyi coccoliths. The equivalent model for PIC reads$${PIC}left(tright)=left(1-dright){PIC}left(t-1right)+{a}_{E}{{max }}left(Eleft(tright)-Eleft(t-1right)right).$$
(6)
Where ({a}_{E}) is the amount of PIC produced per cell, and displayed in panel G. Using the same procedure as for TEP, we obtain the best fit for d = 0.0075. Our PIC model assumes that all PIC production comes from E. huxleyi, supported by large occurrence of E. huxleyi cells observed in scanning electron microscopy (Supplementary Fig. 1).Methods for data collectionMesocosm core setupThe mesocosm experiment AQUACOSM VIMS-Ehux was carried out for 24 days between 24th May (day 0) and 16th June (day 23) 2018 in Raunefjorden at the University of Bergen’s Marine Biological Station Espegrend, Norway (60°16′11 N; 5°13′07E). The experiment consisted of seven enclosure bags made of transparent polyethylene (11 m3, 4 m deep and 2 m wide, permeable to 90% photosynthetically active radiation) mounted on floating frames and moored to a raft in the middle of the fjord. The bags were filled with surrounding fjord water (day −1; pumped from 5 m depth) and continuously mixed by aeration (from day 0 onwards). Each bag was supplemented with nutrients at a nitrogen to phosphorus ratio of 16:1 according to the optimal Redfield Ratio (1.6 µM NaNO3 and 0.1 µM KH2PO4 final concentration) on days 0–5 and 14–17, whereas on days 6, 7 and 13 only nitrogen was added to limit the growth of pico-eukaryotes and favor the growth of E. huxleyi that is more resistant to phosphate limited conditions. Silica was not added as a nutrient source in order to suppress the growth of diatoms and to enhance E. huxleyi proliferation. Bags 5, 6, 7 were covered to collect aerosols and guarantee minimal contamination while sampling for core variables. Bags 1, 2, 3, 4 were sampled for additional assays such as metabolomics, polysaccharides profiling, and vesicles, which increase sampling time and potential for contamination.Measurement of dissolved inorganic nutrientsUnfiltered seawater aliquots (10 mL) were collected from each bag and the surrounding fjord water in 12 mL polypropylene tubes and stored frozen at −20 °C. Dissolved inorganic nutrients were measured with standard segmented flow analysis with colorimetric detection73, using a Bran & Luebe autoanalyser. Data are available in ref. 74 and values for individual bags are plotted in Supplementary Fig. 26.Measurement of water temperature and salinityWater temperature and salinity were measured in each bag and the surrounding fjord water using a SD204 CTD/STD (SAIV A/S, Laksevag, Norway). Data points were averaged for 1–3 m depth (descending only). When this depth was not available, the available data points were taken. Data are missing for the fjord in days 0–1. Outliers were removed for the following samples: bag 1 at days 0, 4, 15; bag 7 at day 15. Data are available in ref. 74.Flow cytometry measurementsSamples for flow cytometric counts were collected twice a day, in the morning (7:00 a.m.) and evening (8:00–9:00 p.m.) from each bag and the surrounding fjord, which served as an environmental reference. Water samples were collected in 50 mL centrifugal tubes from 1 m depth, pre-filtered using 40 µm cell strainers, and immediately analyzed with an Eclipse iCyt (Sony Biotechology, Champaign, IL, USA) flow cytometer. A total volume of 300 µL with a flow rate of 150 µL/min was analyzed with the machine’s software ec800 v1.3.7. A threshold was applied based on the forward scatter signal to reduce the background noise.Phytoplankton populations were identified by plotting the autofluorescence of chlorophyll versus phycoerythrin and side scatter: calcified E. huxleyi (high side scatter and high chlorophyll), Synechococcus (high phycoerythrin and low chlorophyll), nano- and picophytoplankton (high and low chlorophyll, respectively). Chlorophyll fluorescence was detected by FL4 (excitation (ex): 488 nm and emission (em): 663–737 nm). Phycoerythrin was detected by FL3 (ex: 488 nm and em: 570–620 nm). Raw.fcs files were extracted and analyzed in R using ‘flowCore’ and ‘ggcyto’ packages and all data are available on Dryad74. In particular, the gating strategy was adapted to each day and each bag and individual plots for each days and each bag can be found in the Dryad link.For bacteria and viral counts, 200 µL of sample were fixed with 4 µL of 20% glutaraldehyde (final concentration of 0.5%) for 1 h at 4 °C and flash frozen. They were thawed and stained with SYBR gold (Invitrogen) that was diluted 1:10,000 in Tris-EDTA buffer, incubated for 20 min at 80 °C and cooled to room temperature. Bacteria and viruses were counted and analyzed using a Cytoflex and identified based on the Violet SSC-A versus FITC-A by comparing to reference samples containing fixed bacteria and viruses from lab cultures. A total volume of 60 µL with a flow rate of 10 µL/min was analyzed. A threshold was applied based on the forward scatter signal to reduce the background noise. For plotting bacteria (Fig. 1h), a moving average of three successive days was used.Enumeration of extracellular EhV abundance by qPCRDNA extracts from filters from the core sampling (see above) were diluted 100 times, and 1 µL was then used for qPCR analysis. EhV abundance was determined by qPCR for the major capsid protein (mcp) gene: 5′-acgcaccctcaatgtatggaagg-3′ (mcp1F) and 5′-rtscrgccaactcagcagtcgt -3′ (mcp94Rv). All reactions were carried out in technical triplicates using water as a negative control. For all reactions, Platinum SYBER Green qPCR SuperMix-UDG with ROX (Invitrogen, Carlsbad, CA, USA) was used as described by the manufacturer. Reactions were performed on a QuantStudio 5 Real-Time PCR System equipped with the QuantStudio Design and Analysis Software version 1.5.1 (Applied Biosystems, Foster City, CA, USA) as follows: 50 °C for 2 min, 95 °C for 5 min, 40 cycles of 95 °C for 15 s, and 60 °C for 30 s. Results were calibrated against serial dilutions of EhV201 DNA at known concentrations, enabling exact enumeration of viruses. Samples showing multiple peaks in melting curve analysis or peaks that were not corresponding to the standard curves were omitted. Data are available in ref. 74. A comparison of viral counts based on flow-cytometry and qPCR is shown in Supplementary Fig. 2.FlowCam analysisSamples for automated flow imaging microcopy were collected once a day in the morning (7:00 a.m.) from each bag and the surrounding fjord, which served as an environmental reference. Water samples were collected in 50 mL centrifugal tubes from 1 m depth, kept at 12 °C in darkness, and analyzed within 2 h of sampling, using a FlowCAM II (Fluid Imaging Technologies Inc., Scarborough, ME, USA) fitted with a 300 µm path length flow cell and a 4× microscope objective. Images were collected using auto-image mode at a rate of 7 frames/second. A sample volume of 10 mL was processed at a flow rate of 0.7 mL/min. Individual objects within each sample were clustered and annotated using the Ecotaxa platform75. Absolute counts for major groups, including the most abundant ciliate category Ciliophora U04, were then exported and normalized by the individual amount of water volume processed for each sample.Data are available under “Flowcam Composite Aquacosm_2018_VIMS-Ehux” project on Ecotaxa.Scanning electron microscopy50 ml of water samples from bags or fjord were collected on polycarbonate filters (0.2 µm pore size, 47 mm diameter, Millipore). The filters were air dried and stored on petri-slides (Millipore) at room temperature. Prior to observation, a small fraction of the filter was cut and coated with 2 nm of iridium using a Safematic CCU-010 coater (Safematic GMBH, Switzerland). Samples were observed on a Zeiss Ultra SEM that was set at a working distance of 6.2 ± 0.1 mm, an acceleration voltage of 3.0 kV and an aperture size of 30 mm. The secondary electron detector was used for image acquisition.Paired dilution experimentPhytoplankton growth and microzooplankton grazing rates were estimated using the dilution method76,77. A slightly modified version of the method was used with only one low dilution level (20%) and an undiluted treatment used78. Rates calculated using this method are considered conservative but accurate when compared with those using multiple dilution levels and a linear regression. Water from bags 1–4 was collected using a peristaltic pump at ~1 m depth and mixed into a 20 L clean carboy. Water was screened through a 200 µm mesh to remove larger mesozooplankton. The collected water was shaded with black plastic and returned to shore. Dilution experiments were set-up in a temperature-controlled room, set to ambient water temperature (±2 °C). Particle-free diluent (FSW) was prepared by gravity filtering whole seawater (WSW) through a 0.45 µm inline filter (PALL Acropak™ Membrane capsule) into a clean carboy. To the FSW, WSW was gently siphoned at a proportion of 20%. The 20% dilution and 100% WSW treatments were prepared in single carboys and then siphoned into triplicate 1.2 L Nalgene™ incubation bottles. To control for nutrient limitation, additional triplicate bottles of 100% WSW were incubated without added nutrients (10 µM nitrate and 1 µM phosphate). The incubation bottles were incubated for 24 h in an outdoor tank maintained at in-situ water temperatures by a flow-through system of ambient seawater. Bottles could float freely, and the seawater inflow caused gentle agitation throughout the 24 h period. A screen was used to mimic light conditions experienced within the mesocosm bags.To quantify viral mortality, we used the paired dilution method79 which involves setting up an extra low dilution level (20%) containing water filtered through a tangential flow filter (TFF) of 100 kDå to remove viral particles. During this experiment, TFF water was produced 1–2 days prior to the dilution experiment, to ensure the chemical composition of the water was as similar as possible, and experiments could be set up in a timely manner.At T0 hours and T24 hours from all dilution experiments, sub-samples were taken for the determination of chlorophyll-a and flow cytometry. For chlorophyll-a, 100–150 mL of seawater was filtered under low vacuum pressure through a 47 mm Whatman GF/F filters (effective pore size 0.7 µm), and then extracted in 7 mL of 97% methanol at 4 °C in the dark for 12 h. All chlorophyll readings were conducted on a Turner TD700 fluorometer80. Methanol blanks were included, and all samples were corrected for phaeophytin using a drop of 10% hydrochloric acid and then reading the sample again81.Water samples (2 × 1 mL) for flow cytometry were taken at T0 and T24 of dilution experiments for the determination of phytoplankton abundances. Water samples were taken in triplicate from T0, and from each bottle at T24. Samples were immediately fixed in 20 µL of glutaraldehyde (final concentration More

Butchart, S. H. M. et al. Global biodiversity: Indicators of recent declines. Science 328, 1164–1168 (2010).ADS 
CAS 

Google Scholar 
McGill, B. J., Dornelas, M., Gotelli, N. J. & Magurran, A. E. Fifteen forms of biodiversity trend in the Anthropogene. Trends Ecol. Evol. 30, 104–113 (2015).
Google Scholar 
Bradshaw, C. J. A., Sodhi, N. S. & Brook, B. W. Tropical turmoil: A biodiversity tragedy in progress. Front. Ecol. Environ. 7, 79–87 (2009).
Google Scholar 
Newbold, T. et al. Global effects of land use on local terrestrial biodiversity. Nature 520, 45–50 (2015).ADS 
CAS 

Google Scholar 
Loreau, M. et al. Biodiversity and ecosystem functioning: Current knowledge and future challenges. Science 294, 804–808 (2001).ADS 
CAS 

Google Scholar 
Hooper, D. U. et al. Effects of biodiversity on ecosystem functioning: A consensus of current knowledge. Ecol. Monogr. 75, 3–35 (2005).
Google Scholar 
Hooper, D. U. et al. A global synthesis reveals biodiversity loss as a major driver of ecosystem change. Nature 486, 105–108 (2012).ADS 
CAS 

Google Scholar 
Balvanera, P. et al. Quantifying the evidence for biodiversity effects on ecosystem functioning and services. Ecol. Lett. 9, 1146–1156 (2006).
Google Scholar 
Cardinale, B. J. et al. Biodiversity loss and its impact on humanity. Nature 486, 59–67 (2012).ADS 
CAS 

Google Scholar 
Pasari, J. R., Levi, T., Zavaleta, E. S. & Tilman, D. Several scales of biodiversity affect ecosystem multifunctionality. Proc. Natl. Acad. Sci. U.S.A. 110, 10219–10222 (2013).ADS 
CAS 

Google Scholar 
Tilman, D., Isbell, F. & Cowles, J. M. Biodiversity and ecosystem functioning. Annu. Rev. Ecol. Evol. Syst. 45, 471–493 (2014).
Google Scholar 
Murphy, G. E. P. & Romanuk, T. N. A meta-analysis of declines in local species richness from human disturbances. Ecol. Evol. 4, 91–103 (2014).
Google Scholar 
Johnson, C. N. et al. Biodiversity losses and conservation responses in the Anthropocene. Science 356, 270–275 (2017).ADS 
CAS 

Google Scholar 
de Coster, G., Banks-Leite, C. & Metzger, J. P. Atlantic forest bird communities provide different but not fewer functions after habitat loss. Proc. R. Soc. B 282, 20142844 (2015).
Google Scholar 
Riemann, J. C., Ndriantsoa, S. H., Rödel, M.-O. & Glos, J. Functional diversity in a fragmented landscape—habitat alterations affect functional trait composition of frog assemblages in Madagascar. Global Ecol. Conserv. 10, 173–183 (2017).
Google Scholar 
McKinney, M. L. & Lockwood, J. L. Biotic homogenization: A few winners replacing many losers in the next mass extinction. Trends Ecol. Evol. 14, 450–453 (1999).CAS 

Google Scholar 
Socolar, J. B., Gilroy, J. J., Kunin, W. E. & Edwards, D. P. How should beta-diversity inform biodiversity conservation?. Trends Ecol. Evol. 31, 67–80 (2016).
Google Scholar 
van der Plas, F. et al. Biotic homogenization can decrease landscape-scale forest multi-functionality. Proc. Natl. Acad. Sci. U.S.A. 113, 3557–3562 (2016).ADS 

Google Scholar 
Mori, A. S., Isbell, F. & Seidl, R. β-diversity, community assembly, and ecosystem functioning. Trends Ecol. Evol. 33, 549–564 (2018).
Google Scholar 
Dehling, J. M. & Dehling, D. M. Conserving ecological functions of frog communities in Borneo requires diverse forest landscapes. Global Ecol. Conserv. 26, e01481 (2021).
Google Scholar 
Hector, A. & Bagchi, R. Biodiversity and ecosystem multifunctionality. Nature 448, 188–190 (2007).ADS 
CAS 

Google Scholar 
Isbell, F. et al. High plant diversity is needed to maintain ecosystem services. Nature 477, 199–202 (2011).ADS 
CAS 

Google Scholar 
Loreau, M., Mouquet, N. & Gonzalez, A. Biodiversity as spatial insurance in heterogeneous landscapes. Proc. Natl. Acad. Sci. U.S.A. 100, 12765–12770 (2003).ADS 
CAS 

Google Scholar 
Seibold, S. et al. Arthropod decline in grasslands and forests is associated with landscape-level drivers. Nature 574, 671–674 (2019).ADS 
CAS 

Google Scholar 
Felipe-Lucia, M. R. et al. Land-use intensity alters networks between biodiversity, ecosystem functions, and services. Proc. Natl. Acad. Sci. U.S.A. 117, 28140–28149 (2020).ADS 
CAS 

Google Scholar 
Tilman, D. Functional diversity in Encyclopedia of biodiversity, Vol. 3. (ed. Levin S. A.) 109–120 (Academic Press, 2001)Cadotte, M. W., Carscadden, K. & Mirotchnick, N. Beyond species: functional diversity and the maintenance of ecological processes and services. J. Appl. Ecol. 48, 1079–1087 (2011).
Google Scholar 
Flynn, D. F. B., Mirotchnick, N., Jain, M., Palmer, M. I. & Naeem, S. Functional and phylogenetic diversity as predictors of biodiversity-ecosystem function relationships. Ecology 92, 1573–1581 (2011).
Google Scholar 
Lean, C. & Maclaurin, J. The value of phylogenetic diversity in Biodiversity conservation and phylogenetic systematics. Topics in Biodiversity and Conservation 14. (eds. Pellens, R., Grandcolas, P.) 19–38 (Springer, 2016).Owen, N. R., Gumbs, R., Gray, C. L. & Faith, D. P. Global conservation of phylogenetic diversity captures more than just functional diversity. Nat. Commun. 10, 859 (2019).ADS 

Google Scholar 
Gumbs, R., Williams, R. C., Lowney, A. M. & Smith, D. Spatial and species-level metrics reveal global patterns of irreplaceable and imperiled gecko phylogenetic diversity. Israel J. Ecol. Evolut. 66, 239–252 (2020).
Google Scholar 
Brooks, D. R., Mayden, R. L. & McLennan, D. A. Phylogeny and biodiversity: Conserving our evolutionary legacy. Trends Ecol. Evol. 7, 55–59 (1992).CAS 

Google Scholar 
Phillimore, A. B. et al. Biogeographical basis of recent phenotypic divergence among birds: a global study of subspecies richness. Evolution 61, 942–957 (2007).
Google Scholar 
Miraldo, A. et al. An Anthropocene map of genetic diversity. Science 353, 1532–1535 (2016).ADS 
CAS 

Google Scholar 
Smith, B. T., Seeholzer, G. F., Harvey, M. G., Cuervo, A. M. & Brumfield, R. T. A latitudinal phylogeographic diversity gradient in birds. PLoS Biol. 15, e2001073 (2017).
Google Scholar 
Tucker, C. M. et al. Assessing the utility of conserving evolutionary history. Biol. Rev. 94, 1740–1760 (2019).
Google Scholar 
Flynn, D. F. B. et al. Loss of functional diversity under land use intensification across multiple taxa. Ecol. Lett. 12, 22–33 (2009).
Google Scholar 
Villéger, S., Miranda, J. R., Hernández, D. F. & Mouillot, D. Contrasting changes in taxonomic vs. functional diversity of tropical fish communities after habitat degradation. Ecological Applications 20, 1512–1522 (2010).Gibbons, J. W. et al. Remarkable amphibian biomass and abundance in an isolated wetland: Implications for wetland conservation. Conserv. Biol. 20, 1457–1465 (2006).
Google Scholar 
Hocking, D. J. & Babbitt, K. J. Amphibian contributions to ecosystem services. Herpetol. Conserv. Biol. 9, 1–17 (2014).
Google Scholar 
Beebee, T. J. C. Amphibian breeding and climate change. Nature 374, 219–220 (1995).ADS 
CAS 

Google Scholar 
Kiesecker, J. M., Blaustein, A. R. & Belden, L. K. Complex causes of amphibian population declines. Nature 410, 681–684 (2001).ADS 
CAS 

Google Scholar 
Cheng, T. L., Rovito, S. M., Wake, D. B. & Vredenburg, V. T. Coincident mass extirpation of neotropical amphibians with the emergence of the infection fungal pathogen Batrachochytrium dendrobatidis. Proc. Natl. Acad. Sci. U.S.A. 108, 9502–9507 (2011).ADS 
CAS 

Google Scholar 
Wake, D. B. & Vredenburg, V. T. Are we in the midst of the sixth mass extinction? A view from the world of amphibians. Proc. Natl. Acad. Sci. U.S.A. 105, 11466–11473 (2008).ADS 
CAS 

Google Scholar 
Ernst, R. & Rödel, M.-O. Patterns of community composition in two tropical tree frog assemblages: Separating spatial structure and environmental effects in disturbed and undisturbed forests. J. Trop. Ecol. 24, 111–120 (2008).
Google Scholar 
Gardner, T. A. et al. The value of primary, secondary, and plantation forests for a Neotropical Herpetofauna. Conserv. Biol. 21, 775–787 (2007).
Google Scholar 
Gardner, T. A., Fitzherbert, E. B., Drewes, R. C., Howell, K. M. & Caro, T. Spatial and temporal patterns of abundance and diversity of an East African leaf litter amphibian fauna. Biotropica 39, 105–113 (2007).
Google Scholar 
Gillespie, G. R. et al. Conservation of amphibians in Borneo: relative value of secondary tropical forest and non-forest habitats. Biol. Cons. 152, 136–144 (2012).
Google Scholar 
Angarita-M., O., Montes-Correa, A. C. & Renjifo, J. M. Amphibians and reptiles of an agroforestry system in the Colombian Caribbean. Amphibian & Reptile Conservation 8, 33–52 (2015).Jiménez-Robles, O., Guayasamin, J. M., Ron, S. R. & De la Riva, I. Reproductive traits associated with species turnover of amphibians in Amazonia and its Andean slopes. Ecol. Evol. 7, 2489–2500 (2017).
Google Scholar 
Ernst, R., Linsenmair, K. E. & Rödel, M.-O. Diversity erosion beyond the species level: dramatic loss of functional diversity after selective logging in two tropical amphibian communities. Biol. Cons. 133, 143–155 (2006).
Google Scholar 
Oda, F. H. et al. Anuran species richness, composition, and breeding habitat preferences: a comparison between forest remnants and agricultural landscapes in Southern Brazil. Zool. Stud. 55, 34 (2016).
Google Scholar 
Sinsch, U., Lümkemann, K., Rosar, K., Schwarz, C. & Dehling, J. M. Acoustic niche partitioning in an anuran community inhabiting an Afromontane wetland (Butare, Rwanda). African Zool. 47, 60–73 (2012).
Google Scholar 
Tumushimire, L., Mindje, M., Sinsch, U. & Dehling, J. M. The anuran diversity of cultivated wetlands in Rwanda: Melting pot of generalists?. Salamandra 56, 99–112 (2020).
Google Scholar 
REMA. Rwanda State of Environment and Outlook Report 2017 – Achieving Sustainable Urbanization. (Rwanda Environment Management Authority, Government of Rwanda, 2017).Su, J. C., Debinski, D. M., Jakubauskas, M. E. & Kindscher, K. Beyond species richness: Community similarity as a measure of cross-taxon congruence for coarse-filter conservation. Conserv. Biol. 18, 167–173 (2004).
Google Scholar 
Gibson, L. et al. Primary forests are irreplaceable for sustaining tropical biodiversity. Nature 478, 378–381 (2011).ADS 
CAS 

Google Scholar 
Zimkus, B. M., Rödel, M.-O. & Hillers, A. Complex patterns of continental speciation: Molecular phylogenetics and biogeography of sub-Saharan puddle frogs (Phrynobatrachus). Mol. Phylogenet. Evol. 55, 883–900 (2010).
Google Scholar 
Dehling, J. M. & Sinsch, U. Partitioning of morphospace in larval and adult reed frogs (Anura: Hyperoliidae: Hyperolius) of the Central African Albertine Rift. Zool. Anz. 280, 65–77 (2019).
Google Scholar 
Mazel, F. et al. Prioritizing phylogenetic diversity captures functional diversity unreliably. Nat. Commun. 9, 2888 (2018).ADS 

Google Scholar 
Haddad, C. F. B. & Prado, C. P. A. Reproductive modes and their unexpected diversity in the Atlantic forest of Brazil. Bioscience 55, 207–217 (2005).
Google Scholar 
Capinha, C., Essl, F., Seebens, H., Moser, D. & Pereira, H. M. The dispersal of alien species redefines biogeography in the Anthropocene. Science 348, 1248–1251 (2015).ADS 
CAS 

Google Scholar 
Alroy, J. Effects of habitat disturbance on tropical forest biodiversity. Proc. Natl. Acad. Sci. U.S.A. 114, 6056–6061 (2017).ADS 
CAS 

Google Scholar 
Dehling, J. M. & Sinsch, U. Diversity of Ptychadena in Rwanda and taxonomic status of P. chrysogaster Laurent, 1954 (Amphibia, Anura, Ptychadenidae). ZooKeys 356, 69–102 (2013).IUCN. The IUCN Red List of Threatened Species. Version 2020–1. https://www.iucnredlist.org (2020).Portillo, F., Greenbaum, E., Menegon, M., Kusamba, C. & Dehling, J. M. Phylogeography and species boundaries of Leptopelis (Anura: Arthroleptidae) from the Albertine Rift. Mol. Phylogenet. Evol. 82, 75–86 (2015).
Google Scholar 
Channing, A., Dehling, J. M., Lötters, S. & Ernst, R. Species boundaries and taxonomy of the African River Frogs (Anura: Pyxicephalidae: Amietia). Zootaxa 4155, 1–76 (2016).CAS 

Google Scholar 
Rödel, M.-O. & Ernst, R. Measuring and monitoring amphibian diversity in tropical forests. I. An evaluation of methods with recommendations for standardization. Ecotropica 10, 1–14 (2004).Channing, A. & Howell, K. M. Amphibians of East Africa. (Chimaira, 2006).Jetz, W. & Pyron, R. A. The interplay of past diversification and evolutionary isolation with present imperilment across the amphibian tree of life. Nat. Ecol. Evolut. 2, 850–858 (2018).
Google Scholar 
Villéger, S., Mason, N. W. & Mouillot, D. New multidimensional functional diversity indices for a multifaceted framework in functional ecology. Ecology 89, 2290–2301 (2008).
Google Scholar 
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 
Faith, D. P. Conservation evaluation and phylogenetic diversity. Biol. Cons. 61, 1–10 (1992).
Google Scholar 
Dehling, D. M. et al. Functional and phylogenetic diversity and assemblage structure of frugivorous birds along an elevational gradient in the tropical Andes. Ecography 37, 1047–1055 (2014).
Google Scholar 
Baselga, A. et al. betapart: partitioning beta diversity into turnover and nestedness components. R package version 1.5.6. https://CRAN.R-project.org/package=betapart (2022).Dehling, D. M. et al. Specialists and generalists fulfil important and complementary functional roles in ecological processes. Funct. Ecol. 35, 1810–1821 (2021).CAS 

Google Scholar 
Dehling, D. M., Barreto, E. & Graham, C. H. The contribution of mutualistic interactions to functional and phylogenetic diversity. Trends Ecol. Evol. https://doi.org/10.1016/j.tree.2022.05.006 (2022).Article 

Google Scholar 
R Core Team. R: a language and environment for statistical computing. (R Foundation for Statistical Computing, 2021). More

German CR, Von Damm KL. Hydrothermal processes. In: Holland HD, Turekian KK and Elderfield H, editors. Treatise geochem, Vol. 6. The oceans and marine geochemistry. Oxford, UK:Elsevier-Pergamon, 2004;181–222.Bell JB, Woulds C, Oevelen DV. Hydrothermal activity, functional diversity and chemoautotrophy are major drivers of seafloor carbon cycling. Sci Rep. 2017;7:1–3.
Google Scholar 
McCollom TM. Geochemical constraints on primary productivity in submarine hydrothermal vent plumes. Deep Res Part I Oceanogr Res Pap. 2000;47:85–101.CAS 

Google Scholar 
Tunnicliffe V, Baross JA, Gebruk AV, Giere O, Holland ME, Koschinsky A, et al. Group report: what are the interactions between biotic processes at vents and physical, chemical, and geological conditions. In: Halbach PE, Tunnicliffe V, and Hein JR, editors. Energy and Mass Transfer in Marine Hydrothermal Systems. Berlin-Dahlem:University Press; 2003;251–70.Nakamura K, Takai K. Theoretical constraints of physical and chemical properties of hydrothermal fluids on variations in chemolithotrophic microbial communities in seafloor hydrothermal systems. Prog Earth Planet Sci. 2014;1:1–24.
Google Scholar 
Wang W, Li Z, Zeng L, Dong C, Shao Z. The oxidation of hydrocarbons by diverse heterotrophic and mixotrophic bacteria that inhabit deep-sea hydrothermal ecosystems. ISME J. 2020;14:1994–2006.CAS 

Google Scholar 
Sinha RK, Krishnan KP, Kurian PJ. Complete genome sequence and comparative genome analysis of Alcanivorax sp. IO_7, a marine alkane-degrading bacterium isolated from hydrothermally-influenced deep seawater of southwest Indian ridge. Genomics 2021;113:884–91.CAS 

Google Scholar 
Li J, Yang J, Sun M, Su L, Wang H, Gao J, et al. Distribution and succession of microbial communities along the dispersal pathway of hydrothermal plumes on the Southwest Indian Ridge. Front Mar Sci. 2020;7:581381.
Google Scholar 
Meier DV, Bach W, Girguis PR, Gruber-Vodicka HR, Reeves EP, Richter M, et al. Heterotrophic Proteobacteria in the vicinity of diffuse hydrothermal venting. Environ Microbiol. 2016;18:4348–68.
Google Scholar 
Li WL, Huang JM, Zhang PW, Cui GJ, Wei ZF, Wu YZ, et al. Periodic and spatial spreading of alkanes and Alcanivorax bacteria in deep waters of the Mariana Trench. Appl Environ Microbiol. 2019;85:e02089–18.CAS 

Google Scholar 
Brooijmans RJW, Pastink MI, Siezen RJ. Hydrocarbon-degrading bacteria: The oil-spill clean-up crew. Micro Biotechnol. 2009;2:587.CAS 

Google Scholar 
Scoma A, Barbato M, Borin S, Daffonchio D, Boon N. An impaired metabolic response to hydrostatic pressure explains Alcanivorax borkumensis recorded distribution in the deep marine water column. Sci Rep. 2016;6:1–3.
Google Scholar 
Lai Q, Wang L, Liu Y, Fu Y, Zhong H, Wang B, et al. Alcanivorax pacificus sp. nov., isolated from a deep-sea pyrene-degrading consortium. Int J Syst Evol Microbiol. 2011;61:1370–4.CAS 

Google Scholar 
Wu Y, Lai Q, Zhou Z, Qiao N, Liu C, Shao Z. Alcanivorax hongdengensis sp. nov., an alkane-degrading bacterium isolated from surface seawater of the straits of Malacca and Singapore, producing a lipopeptide as its biosurfactant. Int J Syst Evol Microbiol. 2009;59:1474–9.CAS 

Google Scholar 
Fernández-Martínez J, Pujalte MJ, García-Martínez J, Mata M, Garay E, Rodríguez-Valera F. Description of Alcanivorax venustensis sp. nov. and reclassification of Fundibacter jadensis DSM 12178T (Bruns and Berthe-Corti 1999) as Alcanivorax jadensis comb. nov., members of the emended genus Alcanivorax. Int J Syst Evol Microbiol. 2003;53:331–8.
Google Scholar 
Radwan SS, Khanafer MM, Al-Awadhi HA. Ability of the so-called obligate hydrocarbonoclastic bacteria to utilize nonhydrocarbon substrates thus enhancing their activities despite their misleading name. BMC Microbiol. 2019;19:1–2.
Google Scholar 
Kalscheuer R, Stöveken T, Malkus U, Reichelt R, Golyshin PN, Sabirova JS, et al. Analysis of storage lipid accumulation in Alcanivorax borkumensis: Evidence for alternative triacylglycerol biosynthesis routes in bacteria. J Bacteriol. 2007;189:918–28.CAS 

Google Scholar 
Timm C, Davy B, Haase K, Hoernle KA, Graham IJ, De Ronde CEJ, et al. Subduction of the oceanic Hikurangi Plateau and its impact on the Kermadec arc. Nat Commun. 2014;5:1–9.
Google Scholar 
Haase KM, Beier C, Bach W, Kleint C, Anderson MO, Rubin K, et al. SO-263 Cruise Report: Tonga Rift. 2018. https://doi.org/10.13140/RG.2.2.23035.16169.Gartman A, Hannington M, Jamieson JW, Peterkin B, Garbe-Schönberg D, Findlay AJ, et al. Boiling-induced formation of colloidal gold in black smoker hydrothermal fluids. Geology 2018;46:39–42.CAS 

Google Scholar 
Falkenberg JJ, Keith M, Haase KM, Bach W, Klemd R, Strauss H, et al. Effects of fluid boiling on Au and volatile element enrichment in submarine arc-related hydrothermal systems. Geochim Cosmochim Acta. 2021;307:105–32.CAS 

Google Scholar 
Peters C, Strauss H, Haase K, Bach W, de Ronde CEJ, Kleint C, et al. SO2 disproportionation impacting hydrothermal sulfur cycling: Insights from multiple sulfur isotopes for hydrothermal fluids from the Tonga-Kermadec intraoceanic arc and the NE Lau Basin. Chem Geol. 2021;586:120586.CAS 

Google Scholar 
Baker ET, Walker SL, Massoth GJ, Resing JA. The NE Lau Basin: Widespread and abundant hydrothermal venting in the back-arc region behind a superfast subduction zone. Front Mar Sci. 2019;6:382.
Google Scholar 
Kim J, Lee KY, Kim JH. Metal-bearing molten sulfur collected from a submarine volcano: Implications for vapor transport of metals in seafloor hydrothermal systems. Geology 2011;39:351–4.CAS 

Google Scholar 
Klose L, Keith M, Hafermaas D, Kleint C, Bach W, Diehl A, et al. Trace element and isotope systematics in vent fluids and sulphides from Maka volcano, North Eastern Lau Spreading Centre: Insights into three-component fluid mixing. Front Earth Sci. 2021;9:1–26.
Google Scholar 
Herlemann DPR, Labrenz M, Jürgens K, Bertilsson S, Waniek JJ, Andersson AF. Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. ISME J. 2011;5:1571–9.CAS 

Google Scholar 
Dede B, Hansen CT, Neuholz R, Schnetger B, Kleint C, Walker S, et al. Niche differentiation of sulfur-oxidizing bacteria (SUP05) in submarine hydrothermal plumes. ISME J. 2022;16:1479–90.CAS 

Google Scholar 
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17:10–2.
Google Scholar 
R Core Team. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2013.Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.CAS 

Google Scholar 
McMurdie PJ, Holmes S. Phyloseq: An R Package for reproducible interactive analysis and graphics of microbiome census data. PLoS ONE. 2013;8:e61217.CAS 

Google Scholar 
Diehl A, Bach W. MARHYS (MARine HYdrothermal Solutions) Database: A global compilation of marine hydrothermal vent fluid, end member, and seawater compositions. Geochem Geophys Geosystems. 2020;21:e2020GC009385.
Google Scholar 
Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.CAS 

Google Scholar 
Pruesse E, Peplies J, Glöckner FO. SINA: Accurate high-throughput multiple sequence alignment of ribosomal RNA genes. Bioinformatics 2012;28:1823–9.CAS 

Google Scholar 
Ludwig W, Strunk O, Westram R, Richter L, Meier H, Yadhukumar A, et al. ARB: A software environment for sequence data. Nucleic Acids Res. 2004;32:1363–71.CAS 

Google Scholar 
Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0. Syst Biol. 2010;59:307–21.CAS 

Google Scholar 
Stamatakis A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 2014;30:1312–3.CAS 

Google Scholar 
Pernthaler A, Pernthaler J, Amann R. Fluorescence in situ hybridization and catalyzed reporter deposition for the identification of marine bacteria. Appl Environ Microbiol. 2002;68:3094–101.CAS 

Google Scholar 
Amann RI, Binder BJ, Olson RJ, Chisholm SW, Devereux R, Stahl DA. Combination of 16S rRNA-targeted oligonucleotide probes with flow cytometry for analyzing mixed microbial populations. Appl Environ Microbiol. 1990;56:1919–25.CAS 

Google Scholar 
Daims H, Brühl A, Amann R, Schleifer KH, Wagner M. The domain-specific probe EUB338 is insufficient for the detection of all Bacteria: Development and evaluation of a more comprehensive probe set. Syst Appl Microbiol. 1999;22:434–44.CAS 

Google Scholar 
Wallner G, Amann R, Beisker W. Optimizing fluorescent in situ hybridization with rRNA‐targeted oligonucleotide probes for flow cytometric identification of microorganisms. Cytometry 1993;14:136–43.CAS 

Google Scholar 
Stahl DA, Amann R. Development and application of nucleic acid probes in bacterial systematics. In: Nucleic acid techniques in bacterial systematics. Stackebrandt, E, Goodfellow M, editors. Chichester, UK: John Wiley & Sons Ltd; 1991. pp. 205–48.Manz W, Amann R, Ludwig W, Wagner M, Schleifer KH. Phylogenetic oligodeoxynucleotide probes for the major subclasses of Proteobacteria: Problems and solutions. Syst Appl Microbiol. 1992;15:593–600.
Google Scholar 
Eilers H, Pernthaler J, Glöckner FO, Amann R. Culturability and in situ abundance of pelagic Bacteria from the North Sea. Appl Environ Microbiol. 2000;66:3044–51.CAS 

Google Scholar 
Syutsubo K, Kishira H, Harayama S. Development of specific oligonucleotide probes for the identification and in situ detection of hydrocarbon-degrading Alcanivorax strains. Environ Microbiol. 2001;3:371–9.CAS 

Google Scholar 
Morris RM, Rappé MS, Urbach E, Connon SA, Giovannoni SJ. Prevalence of the Chloroflexi-related SAR202 bacterioplankton cluster throughout the mesopelagic zone and deep ocean. Appl Environ Microbiol. 2004;70:2836–42.CAS 

Google Scholar 
Bushnell B BBMap (version 35.14). 2015. https://sourceforge.net/projects/bbmap/.Andrews S. FastQC: A quality control tool for high throughput sequence data. Babraham Bioinforma. 2010; http://www.bioinformatics.babraham.ac.uk/projects/.Rodriguez-R LM, Gunturu S, Tiedje JM, Cole JR, Konstantinidis KT. Nonpareil 3: Fast estimation of metagenomic coverage and sequence diversity. mSystems 2018;3:e00039–18.
Google Scholar 
Menzel P, Ng KL, Krogh A. Fast and sensitive taxonomic classification for metagenomics with Kaiju. Nat Commun. 2016;7:1–9.
Google Scholar 
Kopylova E, Noé L, Touzet H. SortMeRNA: Fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics 2012;28:3211–7.CAS 

Google Scholar 
Li D, Liu CM, Luo R, Sadakane K, Lam TW. MEGAHIT: An ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 2015;31:1674–6.CAS 

Google Scholar 
Gurevich A, Saveliev V, Vyahhi N, Tesler G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 2013;29:1072–5.CAS 

Google Scholar 
Alneberg J, Bjarnason BS, De Bruijn I, Schirmer M, Quick J, Ijaz UZ, et al. Binning metagenomic contigs by coverage and composition. Nat Methods. 2014;11:1144–6.CAS 

Google Scholar 
Eren AM, Kiefl E, Shaiber A, Veseli I, Miller SE, Schechter MS, et al. Community-led, integrated, reproducible multi-omics with anvi’o. Nat Microbiol. 2021;6:3–6.CAS 

Google Scholar 
Bankevich A, Nurk S, Antipov D, Gurevich AA, Dvorkin M, Kulikov AS, et al. SPAdes: A new genome assembly algorithm and its applications to single-cell sequencing. J Comput Biol. 2012;19:455–77.CAS 

Google Scholar 
Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. CheckM: Assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 2015;25:1043–55.CAS 

Google Scholar 
Bowers RM, Kyrpides NC, Stepanauskas R, Harmon-Smith M, Doud D, Reddy TBK, et al. Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea. Nat Biotech. 2017;35:725–31.CAS 

Google Scholar 
Chaumeil P-A, Mussig AJ, Hugenholtz P, Parks DH. GTDB-Tk: a toolkit to classify genomes with the Genome Taxonomy Database. Bioinformatics 2019;36:1925–7.
Google Scholar 
Seemann T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 2014;30:2068–9.CAS 

Google Scholar 
Priest T, Heins A, Harder J, Amann R, Fuchs BM. Niche partitioning of the ubiquitous and ecologically relevant NS5 marine group. ISME J. 2022;16:1570–82.CAS 

Google Scholar 
Eddy SR. Accelerated profile HMM searches. PLoS Comput Biol. 2011;7:e1002195.CAS 

Google Scholar 
Karthikeyan S, Rodriguez‐R LM, Heritier‐Robbins P, Hatt JK, Huettel M, Kostka JE, et al. Genome repository of oil systems: An interactive and searchable database that expands the catalogued diversity of crude oil‐associated microbes. Environ Microbiol. 2020;22:2094–106.CAS 

Google Scholar 
Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49:W293–6.CAS 

Google Scholar 
Arndt D, Grant JR, Marcu A, Sajed T, Pon A, Liang Y, et al. PHASTER: a better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016;44:W16–21.CAS 

Google Scholar 
Bolger AM, Lohse M, Usadel B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014;30:2114–20.CAS 

Google Scholar 
Gomes AÉI, Stuchi LP, Siqueira NMG, Henrique JB, Vicentini R, Ribeiro ML, et al. Selection and validation of reference genes for gene expression studies in Klebsiella pneumoniae using Reverse Transcription Quantitative real-time PCR. Sci Rep. 2018;8:1–4.
Google Scholar 
Guidi L, Chaffron S, Bittner L, Eveillard D, Larhlimi A, Roux S, et al. Plankton networks driving carbon export in the oligotrophic ocean. Nature 2016;532:465–70.CAS 

Google Scholar 
Duarte CM. Seafaring in the 21st century: the Malaspina 2010 circumnavigation expedition. Limnol Oceanogr Bull. 2015;24:11–4.
Google Scholar 
Anantharaman K, Breier JA, Dick GJ. Metagenomic resolution of microbial functions in deep-sea hydrothermal plumes across the Eastern Lau Spreading Center. ISME J. 2016;10:225–39.CAS 

Google Scholar 
Waite DW, Vanwonterghem I, Rinke C, Parks DH, Zhang Y, Takai K, et al. Comparative genomic analysis of the class Epsilonproteobacteria and proposed reclassification to Epsilonbacteraeota (phyl. nov.). Front Microbiol. 2017;8:682.
Google Scholar 
Waite DW, Vanwonterghem I, Rinke C, Parks DH, Zhang Y, Takai K, et al. Addendum: Comparative genomic analysis of the class Epsilonproteobacteria and proposed reclassification to Epsilonbacteraeota (phyl.nov.). Front Microbiol. 2017;9:772.
Google Scholar 
Green DH, Llewellyn LE, Negri AP, Blackburn SI, Bolch CJS. Phylogenetic and functional diversity of the cultivable bacterial community associated with the paralytic shellfish poisoning dinoflagellate Gymnodinium catenatum. FEMS Microbiol Ecol. 2004;47:345–57.CAS 

Google Scholar 
Ramasamy KP, Rajasabapathy R, Lips I, Mohandass C, James RA. Genomic features and copper biosorption potential of a new Alcanivorax sp. VBW004 isolated from the shallow hydrothermal vent (Azores, Portugal). Genomics 2020;112:3268–73.CAS 

Google Scholar 
Barbato M, Scoma A, Mapelli F, De Smet R, Banat IM, Daffonchio D, et al. Hydrocarbonoclastic Alcanivorax isolates exhibit different physiological and expression responses to N-dodecane. Front Microbiol. 2016;7:2056.
Google Scholar 
Sevilla E, Yuste L, Rojo F. Marine hydrocarbonoclastic bacteria as whole-cell biosensors for n-alkanes. Micro Biotechnol. 2015;8:693–706.CAS 

Google Scholar 
Tivey MK. Black and white smokers. In: Harff J, Meschede M, Petersen S, Thiede Jö, editors. Encyclopedia of Marine Geosciences. Dordrecht: Springer Netherlands; 2016. p. 58–62.Djurhuus A, Mikalsen SO, Giebel HA, Rogers AD. Cutting through the smoke: The diversity of microorganisms in deep-sea hydrothermal plumes. R Soc Open Sci. 2017;4:160829.
Google Scholar 
Leahy JG, Colwell RR. Microbial degradation of hydrocarbons in the environment. Microbiol Rev. 1990;54:305–15.CAS 

Google Scholar 
Atlas R, Bragg J. Bioremediation of marine oil spills: When and when not – The Exxon Valdez experience. Micro Biotechnol. 2009;2:213–21.CAS 

Google Scholar 
Reva ON, Hallin PF, Willenbrock H, Sicheritz-Ponten T, Tümmler B, Ussery DW. Global features of the Alcanivorax borkumensis SK2 genome. Environ Microbiol. 2008;10:614–25.CAS 

Google Scholar 
Gregory GJ, Morreale DP, Carpenter MR, Kalburge SS, Boyd EF. Quorum sensing regulators AphA and OpaR control expression of the operon responsible for biosynthesis of the compatible solute ectoine. Appl Environ Microbiol. 2019;85:e01543–19.CAS 

Google Scholar 
Richter AA, Mais CN, Czech L, Geyer K, Hoeppner A, Smits SHJ, et al. Biosynthesis of the stress-protectant and chemical chaperon ectoine: biochemistry of the transaminase EctB. Front Microbiol. 2019;10:2811.
Google Scholar 
Schneiker S, Dos Santos VAPM, Bartels D, Bekel T, Brecht M, Buhrmester J, et al. Genome sequence of the ubiquitous hydrocarbon-degrading marine bacterium Alcanivorax borkumensis. Nat Biotechnol. 2006;24:997–1004.CAS 

Google Scholar 
Wang W, Shao Z. Enzymes and genes involved in aerobic alkane degradation. Front Microbiol. 2013;4:116.
Google Scholar 
Barclay W, Rodd JA, Pflueger JC, Havard KR, Helu SP. Oil plays in the kingdom of Tonga, Southwest Pacific. PESA J. 1993;21:79–92.
Google Scholar 
Chadwick WW, Rubin KH, Merle SG, Bobbitt AM, Kwasnitschka T, Embley RW. Recent eruptions between 2012-2018 discovered at West Mata submarine volcano (NE Lau Basin, SW Pacific) and characterized by new ship, AUV, and ROV data. Front Mar Sci. 2019;6:495.
Google Scholar 
Baumberger T, Lilley MD, Lupton JE, Baker ET, Resing JA, Buck NJ, et al. Dissolved gas and metal composition of hydrothermal plumes from a 2008 submarine eruption on the Northeast Lau Spreading Center. Front Mar Sci. 2020;7:171.
Google Scholar 
Lupton J, Rubin KH, Arculus R, Lilley M, Butterfield D, Resing J, et al. Helium isotope, C/3 He, and Ba‐Nb‐Ti signatures in the northern Lau Basin: Distinguishing arc, back‐arc, and hotspot affinities. Geochem Geophys. 2015;16:1133–55.CAS 

Google Scholar 
Graham DW. Noble gas isotope geochemistry of mid-ocean ridge and ocean island basalts: Characterization of mantle source reservoirs. In: Porcelli D, Wieler R, Ballentine C, editors. Noble gases in Geochemistry and cosmochemistry, Rev Mineral Geochem. Vol 47. Washington D.C.: Mineral Soc. Of Am; 2002. p. 247–318.Lupton JE, Arculus RJ, Greene RR, Evans LJ, Goddard CI. Helium isotope variations in seafloor basalts from the Northwest Lau Backarc Basin: Mapping the influence of the Samoan hotspot. Geophys Res Lett. 2009;36:L17313.
Google Scholar 
Gordon GW. Naturally occurring organohalogen compounds – A comprehensive survey. Prog Chem Org Nat Prod. 1996;68:1–423.
Google Scholar 
Spietz RL, Butterfield DA, Buck NJ, Larson BI, Chadwick WW, Walker SL, et al. Deep-sea volcanic eruptions create unique chemical and biological linkages between the subsurface lithosphere and the oceanic hydrosphere. Oceanography. 2018;31:128–35.
Google Scholar 
Huber JA, Butterfield DA, Baross JA. Bacterial diversity in a subseafloor habitat following a deep-sea volcanic eruption. FEMS Microbiol Ecol. 2003;43:393–409.CAS 

Google Scholar  More

Brown, J. S., Laundre, J. W. & Gurung, M. The ecology of fear: optimal foraging, game theory, and trophic interactions. J. Mammal. 80, 385–399 (1999).
Google Scholar 
Laundré, J. W., Hernández, L. & Altendorf, K. B. Wolves, elk, and bison: reestablishing the ‘landscape of fear’ in Yellowstone National Park, US.A. Can. J. Zool. 79, 1401–1409 (2001).
Google Scholar 
Atkins, J. L. et al. Cascading impacts of large-carnivore extirpation in an African ecosystem. Science 364, 173–177 (2019).CAS 

Google Scholar 
Laundre, J. W., Hernandez, L. & Ripple, W. J. The landscape of fear: ecological implications of being afraid. Open Ecol. J. 3, 1–7 (2010).
Google Scholar 
Loggins, A. A., Shrader, A. M., Monadjem, A. & McCleery, R. A. Shrub cover homogenizes small mammals’ activity and perceived predation risk. Sci. Rep. 9, 16857 (2019).
Google Scholar 
Whittingham, M. J. & Evans, K. L. The effects of habitat structure on predation risk of birds in agricultural landscapes. Ibis 146, 210–220 (2004).
Google Scholar 
Marshall, K. L. A., Philpot, K. E. & Stevens, M. Microhabitat choice in island lizards enhances camouflage against avian predators. Sci. Rep. 6, 19815 (2016).CAS 

Google Scholar 
Stevens, M., Troscianko, J., Wilson-Aggarwal, J. K. & Spottiswoode, C. N. Improvement of individual camouflage through background choice in ground-nesting birds. Nat. Ecol. Evol. 1, 1325–1333 (2017).
Google Scholar 
Wilson-Aggarwal, J. K., Troscianko, J. T., Stevens, M. & Spottiswoode, C. N. Escape distance in ground-nesting birds differs with individual level of camouflage. Am. Nat. 188, 231–239 (2016).
Google Scholar 
Troscianko, J., Wilson-Aggarwal, J., Stevens, M. & Spottiswoode, C. N. Camouflage predicts survival in ground-nesting birds. Sci. Rep. 6, 19966 (2016).CAS 

Google Scholar 
Gaston, K. J., Duffy, J. P., Gaston, S., Bennie, J. & Davies, T. W. Human alteration of natural light cycles: causes and ecological consequences. Oecologia 176, 917–931 (2014).
Google Scholar 
Gaston, K. J., Davies, T. W., Nedelec, S. L. & Holt, L. A. Impacts of artificial light at night on biological timings. Annu. Rev. Ecol. Evol. Syst. 48, 49–68 (2017).
Google Scholar 
Falchi, F. et al. The new world atlas of artificial night sky brightness. Sci. Adv. 2, e1600377 (2016).
Google Scholar 
Gaston, K. J. et al. Pervasiveness of biological impacts of artificial light at night. Integr. Comp. Biol. 61, 1098–1110 (2021).
Google Scholar 
Sanders, D., Frago, E., Kehoe, R., Patterson, C. & Gaston, K. J. A meta-analysis of biological impacts of artificial light at night. Nat. Ecol. Evol. 5, 74–81 (2021).
Google Scholar 
Kronfeld-Schor, N., Visser, M. E., Salis, L. & van Gils, J. A. Chronobiology of interspecific interactions in a changing world. Philos. Trans. R. Soc. B Biol. Sci. 372, 20160248 (2017).
Google Scholar 
Underwood, C. N., Davies, T. W. & Queir Os, A. M. Artificial light at night alters trophic interactions of intertidal invertebrates. J. Anim. Ecol. 86, 781–789 (2017).
Google Scholar 
Burger, J., Howe, M. A., Hahn, D. C. & Chase, J. Effects of tide cycles on habitat selection and habitat partitioning by migrating shorebirds. Auk 94, 743–758 (1977).
Google Scholar 
Granadeiro, J. P., Dias, M. P., Martins, R. C. & Palmeirim, J. M. Variation in numbers and behaviour of waders during the tidal cycle: implications for the use of estuarine sediment flats. Acta Oecologica 29, 293–300 (2006).
Google Scholar 
Lourenço, P. M. et al. The energetic importance of night foraging for waders wintering in a temperate estuary. Acta Oecologica 34, 122–129 (2008).
Google Scholar 
McNeil, R., Drapeau, P. & Goss-Custard, J. D. The occurrence and adaptive significance of nocturnal habits in waterfowl. Biol. Rev. 67, 381–419 (1992).
Google Scholar 
Martin, G. R. Visual fields and their functions in birds. J. Ornithol. 148, 547–562 (2007).
Google Scholar 
Martin, G. R. What is binocular vision for? A birds’ eye view. J. Vis. 9, 1–19 (2009).
Google Scholar 
Davies, T. W., Duffy, J. P., Bennie, J. & Gaston, K. J. The nature, extent, and ecological implications of marine light pollution. Front. Ecol. Environ. 12, 347–355 (2014).
Google Scholar 
Leopold, M. F., Philippart, C. J. M. & Yorio, P. Nocturnal feeding under artificial light conditions by Brown-hooded Gull (Larus maculipennis) in Puerto Madryn harbour (Chubut Province, Argentina). Hornero 25, 55–60 (2010).
Google Scholar 
Pugh, A. R. & Pawson, S. M. Artificial light at night potentially alters feeding behaviour of the native southern black-backed gull (Larus dominicanus). Notornis 63, 37–39 (2016).
Google Scholar 
Santos, C. D. et al. Effects of artificial illumination on the nocturnal foraging of waders. Acta Oecologica 36, 166–172 (2010).
Google Scholar 
Montevecchi, W. A. Influences of Artificial Light on Marine Birds. in Ecological Consequences of Artificial Night Lighting (eds. Rich, C. & Longcore, T.) 94–113 (Island Press, 2006).Dwyer, R. G., Bearhop, S., Campbell, H. A. & Bryant, D. M. Shedding light on light: benefits of anthropogenic illumination to a nocturnally foraging shorebird. J. Anim. Ecol. 82, 478–485 (2013).
Google Scholar 
Blumstein, D. T. Developing an evolutionary ecology of fear: how life history and natural history traits affect disturbance tolerance in birds. Anim. Behav. 71, 389–399 (2006).
Google Scholar 
Stankowich, T. & Blumstein, D. T. Fear in animals: a meta-analysis and review of risk assessment. Proc. R. Soc. B Biol. Sci. 272, 2627–2634 (2005).
Google Scholar 
Caro, T. Antipredator Defenses in Birds and Mammals. (University of Chicago Press, 2005).Tillmann, J. E. Fear of the dark: night-time roosting and anti-predation behaviour in the grey partridge (Perdix perdix L.). Behaviour 146, 999–1023 (2009).
Google Scholar 
IUCN. The IUCN Red List of Threatened Species. Version 2022-1. https://www.iucnredlist.org/species/22693190/117917038 (2022).Brown, D. et al. The Eurasian Curlew—the most pressing bird conservation priority in the UK? Br. Birds 108, 660–668 (2015).
Google Scholar 
Franks, S. E., Douglas, D. J. T., Gillings, S. & Pearce-Higgins, J. W. Environmental correlates of breeding abundance and population change of Eurasian Curlew Numenius arquata in Britain. Bird. Study 64, 393–409 (2017).
Google Scholar 
Desholm, M. & Kahlert, J. Avian collision risk at an offshore wind farm. Biol. Lett. 1, 296–298 (2005).
Google Scholar 
Clarke, J. A. Moonlight’s influence on predator/prey interactions between short-eared owls (Asio flammeus) and Deermice (Peromyscus maniculatus). Behav. Ecol. Sociobiol. 13, 205–209 (1983).
Google Scholar 
Mandelik, Y., Jones, M. & Dayan, T. Structurally complex habitat and sensory adaptations mediate the behavioural responses of a desert rodent to an indirect cue for increased predation risk. Evol. Ecol. Res. 5, 501–515 (2003).
Google Scholar 
Alexander, R. D. The Evolution of Social Behavior | Annual Review of Ecology, Evolution, and Systematics. Annu. Rev. Ecol. Syst. 5, 325–383 (1974).
Google Scholar 
Pulliam, H. R. On the advantages of flocking. J. Theor. Biol. 38, 419–422 (1973).CAS 

Google Scholar 
Barnard, C. J. Flock feeding and time budgets in the house sparrow (Passer domesticus L.). Anim. Behav. 28, 295–309 (1980).
Google Scholar 
Cooper, W. E. Jr. et al. Effects of risk, cost, and their interaction on optimal escape by nonrefuging Bonaire whiptail lizards, Cnemidophorus murinus. Behav. Ecol. 14, 288–293 (2003).
Google Scholar 
Lagos, P. A. et al. Flight initiation distance is differentially sensitive to the costs of staying and leaving food patches in a small-mammal prey. Can. J. Zool. 87, 1016–1023 (2009).
Google Scholar 
Ydenberg, R. C. & Dill, L. M. The economics of fleeing from predators. Adv. Study Behav. 16, 229–249 (1986).
Google Scholar 
Tucker, V. A., Tucker, A. E., Akers, K. & Enderson, J. H. Curved flight paths and sideways vision in peregrine falcons (Falco peregrinus). J. Exp. Biol. 203, 3755–3763 (2000).CAS 

Google Scholar 
Carr, J. M. & Lima, S. L. Wintering birds avoid warm sunshine: predation and the costs of foraging in sunlight. Oecologia 174, 713–721 (2014).
Google Scholar 
van den Hout, P. J. & Martin, G. R. Extreme head-tilting in shorebirds: predator detection and sun avoidance. Wader Study Group Bull. 118, 18–21 (2011).
Google Scholar 
Ferguson, J. W. H., Galpin, J. S. & de Wet, M. J. Factors affecting the activity patterns of black-backed jackals Canis mesomelas. J. Zool. 214, 55–69 (1988).
Google Scholar 
Pyke, G. H. Optimal foraging theory: a critical review. Annu. Rev. Ecol. Syst. 15, 523–575 (1984).
Google Scholar 
Stephens, D. W. & Krebs, J. R. Foraging Theory. (Princeton University Press, 1986).Mouritsen, K. N. Predator avoidance in night-feeding dunlins calidris alpina: a matter of concealment. Ornis Scand. 23, 195–198 (1992).
Google Scholar 
Blumstein, D. T. Flight-initiation distance in birds is dependent on intruder starting distance. J. Wildl. Manag. 67, 852–857 (2003).
Google Scholar 
Troscianko, J. OSpRad; an open-source, low-cost, high-sensitivity spectroradiometer (p. 2022.12.09.519768). bioRxiv https://doi.org/10.1101/2022.12.09.519768 (2022).Article 

Google Scholar 
Hartig, F. DHARMa: residual diagnostics for hierarchical (multi-level/mixed) regression models. R package version 0.4.4. http://florianhartig.github.io/DHARMa/ (2022).Core Team, R. R: a Language and Environment for Statistical Computing. (R Foundation for Statistical Computing, Vienna, 2022).
Google Scholar  More

Lovegrove, B. G. A phenology of the evolution of endothermy in birds and mammals. Biol. Rev. 92, 1213–1240 (2017).
Google Scholar 
Cuthill, I. C. et al. The biology of color. Science 357, 1–7 (2017).
Google Scholar 
Stuart-Fox, D., Newton, E. & Clusella-Trullas, S. Thermal consequences of colour and near-infrared reflectance. Philos. Trans. R. Soc. B Biol. Sci. 372, 20160345 (2017).
Google Scholar 
Smith, K. R. et al. Color change for thermoregulation versus camouflage in free-ranging lizards. Am. Nat. 188, 668–678 (2016).
Google Scholar 
Rudh, A. & Qvarnström, A. Adaptive colouration in amphibians. Semin. Cell Dev. Biol. 24, 553–561 (2013).
Google Scholar 
Geen, M. R. S. & Johnston, G. R. Coloration affects heating and cooling in three color morphs of the Australian bluetongue lizard, Tiliqua scincoides. J. Therm. Biol. 43, 54–60 (2014).
Google Scholar 
Tattersall, G. J., Eterovick, P. C. & de Andrade, D. V. Tribute to R. G. Boutilier: skin colour and body temperature changes in basking Bokermannohyla alvarengai (Bokermann 1956). J. Exp. Biol. 209, 1185–1196 (2006).
Google Scholar 
Tattersall, G. J., Hillman, S. S., Drewes, R. C. & Sokol, O. M. The thermogenesis of digestion in rattlesnakes. J. Exp. Biol. 207, 579–585 (2004).
Google Scholar 
Seebacher, F. & Murray, S. A. Transient receptor potential ion channels control thermoregulatory behaviour in reptiles. PLoS One 2, e281, 1–7 (2007).Forget-Klein, É. & Green, D. M. Toads use the subsurface thermal gradient for temperature regulation underground. J. Therm. Biol. 99, 1–9 (2021).
Google Scholar 
Kiefer, M. C., Van Sluys, M. & Rocha, C. F. D. Thermoregulatory behaviour in Tropidurus torquatus (Squamata, Tropiduridae) from Brazilian coastal populations: an estimate of passive and active thermoregulation in lizards. Acta Zool. 88, 81–87 (2007).
Google Scholar 
Spencer, K. et al. Growth at cold temperature increases the number of motor neurons to optimize locomotor function. Curr. Biol. 29, 1787–1799.e5 (2019).CAS 

Google Scholar 
Herrel, A. & Bonneaud, C. Temperature dependence of locomotor performance in the tropical clawed frog, Xenopus tropicalis. J. Exp. Biol. 215, 2465–2470 (2012).
Google Scholar 
Casterlin, M. E. & Reynolds, W. W. Diel activity and thermoregulatory behavior of a fully aquatic frog: Xenopus laevis. Hydrobiologia 75, 189–191 (1980).
Google Scholar 
Guo, K. et al. The thermal dependence and molecular basis of physiological color change in Takydromus septentrionalis (Lacertidae). Biol. Open 10, 1–9 (2021).
Google Scholar 
De Velasco, J. B. & Tattersall, G. J. The influence of hypoxia on the thermal sensitivity of skin colouration in the bearded dragon, Pogona vitticeps. J. Comp. Physiol. B. 178, 867–875 (2008).CAS 

Google Scholar 
Stuart-Fox, D. & Moussalli, A. Camouflage, communication and thermoregulation: lessons from colour changing organisms. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 364, 463–470 (2009).
Google Scholar 
Sanabria, E. A., Vaira, M., Quiroga, L. B., Akmentins, M. S. & Pereyra, L. C. Variation of thermal parameters in two different color morphs of a diurnal poison toad, Melanophryniscus rubriventris (Anura: Bufonidae). J. Therm. Biol. 41, 1–5 (2014).
Google Scholar 
Clusella-Trullas, S., van Wyk, J. H. & Spotila, J. R. Thermal benefits of melanism in cordylid lizards: a theoretical and field test. Ecology 90, 2297–2312 (2009).
Google Scholar 
Duarte, R. C., Flores, A. A. V. & Stevens, M. Camouflage through colour change: mechanisms, adaptive value and ecological significance. Philos. Trans. R. Soc. B: Biol. Sci. 372, 1–7 (2017).Bertolesi, G. E. & McFarlane, S. Seeing the light to change colour: an evolutionary perspective on the role of melanopsin in neuroendocrine circuits regulating light-mediated skin pigmentation. Pigment Cell Melanoma Res. 31, 354–373 (2018).CAS 

Google Scholar 
Bertolesi, G. E. et al. The regulation of skin pigmentation in response to environmental light by pineal type II opsins and skin melanophore melatonin receptors. J. Photochem. Photobiol. B Biol. 212, 112024 (2020).CAS 

Google Scholar 
Bagnara, J. T. Pineal regulation of the body lightening reaction in amphibian larvae. Sci. (80-.). 132, 1481–1483 (1960).CAS 

Google Scholar 
Bertolesi, G. E., Song, Y. N., Atkinson-Leadbeater, K., Yang, J.-L. J. & McFarlane, S. Interaction and developmental activation of two neuroendocrine systems that regulate light-mediated skin pigmentation. Pigment Cell Melanoma Res. 30, 413–423 (2017).CAS 

Google Scholar 
Wang, H. & Siemens, J. TRP ion channels in thermosensation, thermoregulation and metabolism. Temp. (Austin, Tex.) 2, 178–187 (2015).
Google Scholar 
Hoffstaetter, L. J., Bagriantsev, S. N. & Gracheva, E. O. TRPs et al.: a molecular toolkit for thermosensory adaptations. Pflug. Arch. Eur. J. Physiol. 470, 745–759 (2018).CAS 

Google Scholar 
Kashio, M. Thermosensation involving thermo-TRPs. Mol. Cell. Endocrinol. 520, 1–8 (2021).
Google Scholar 
Señarís, R., Ordás, P., Reimúndez, A. & Viana, F. Mammalian cold TRP channels: impact on thermoregulation and energy homeostasis. Pflug. Arch. 470, 761–777 (2018).
Google Scholar 
Guo, H., Carlson, J. A. & Slominski, A. Role of TRPM in melanocytes and melanoma. Exp. Dermatol. 21, 650–654 (2012).CAS 

Google Scholar 
Kadowaki, T. Evolutionary dynamics of metazoan TRP channels. Pflug. Arch. 467, 2043–2053 (2015).CAS 

Google Scholar 
Saito, S. & Tominaga, M. Evolutionary tuning of TRPA1 and TRPV1 thermal and chemical sensitivity in vertebrates. Temp. (Austin, Tex.) 4, 141–152 (2017).
Google Scholar 
Saito, S. et al. Analysis of transient receptor potential ankyrin 1 (TRPA1) in frogs and lizards illuminates both nociceptive heat and chemical sensitivities and coexpression with TRP vanilloid 1 (TRPV1) in ancestral vertebrates. J. Biol. Chem. 287, 30743–30754 (2012).CAS 

Google Scholar 
Saito, S. et al. Evolution of heat sensors drove shifts in thermosensation between xenopus species adapted to different thermal niches. J. Biol. Chem. 291, 11446–11459 (2016).CAS 

Google Scholar 
Gracheva, E. O. et al. Molecular basis of infrared detection by snakes. Nature 464, 1006–1011 (2010).CAS 

Google Scholar 
Laursen, W. J., Anderson, E. O., Hoffstaetter, L. J., Bagriantsev, S. N. & Gracheva, E. O. Species-specific temperature sensitivity of TRPA1. Temp. (Austin, Tex.) 2, 214–226 (2015).
Google Scholar 
Bertolesi, G. E., Hehr, C. L. & McFarlane, S. Melanopsin photoreception in the eye regulates light-induced skin colour changes through the production of α-MSH in the pituitary gland. Pigment Cell Melanoma Res. 28, 559–571 (2015).CAS 

Google Scholar 
Bagnara, J. T. The pineal and the body lightening reaction of larval amphibians. Gen. Comp. Endocrinol. 3, 86–100 (1963).CAS 

Google Scholar 
Nisembaum, L. et al. In the heat of the night: thermo-TRPV channels in the salmonid pineal photoreceptors and modulation of melatonin secretion. Endocrinology 156, 4629–4638 (2015).CAS 

Google Scholar 
Schartl, M. et al. What is a vertebrate pigment cell? Pigment Cell Melanoma Res. 29, 8–14 (2016).
Google Scholar 
Slominski, A. Cooling skin cancer: menthol inhibits melanoma growth. Focus on ‘TRPM8 activation suppresses cellular viability in human melanoma’. Am. J. Physiol. – Cell Physiol. 295, C293–C295 (2008).CAS 

Google Scholar 
Yamamura, H., Ugawa, S., Ueda, T., Morita, A. & Shimada, S. TRPM8 activation suppresses cellular viability in human melanoma. Am. J. Physiol. Cell Physiol. 295, C296–C301 (2008).CAS 

Google Scholar 
Knowlton, W. M. et al. A sensory-labeled line for cold: TRPM8-expressing sensory neurons define the cellular basis for cold, cold pain, and cooling-mediated analgesia. J. Neurosci. 33, 2837–2848 (2013).CAS 

Google Scholar 
Weyer-Menkhoff, I., Pinter, A., Schlierbach, H., Schänzer, A. & Lötsch, J. Epidermal expression of human TRPM8, but not of TRPA1 ion channels, is associated with sensory responses to local skin cooling. Pain 160, 2699–2709 (2019).Kumasaka, M., Sato, S., Yajima, I. & Yamamoto, H. Isolation and developmental expression of tyrosinase family genes in Xenopus laevis. Pigment Cell Res. 16, 455–462 (2003).CAS 

Google Scholar 
Rodionov, V. I., Hope, A. J., Svitkina, T. M. & Borisy, G. G. Functional coordination of microtubule-based and actin-based motility in melanophores. Curr. Biol. 8, 165–169 (1998).CAS 

Google Scholar 
Session, A. M. et al. Genome evolution in the allotetraploid frog Xenopus laevis. Nature 538, 336–343 (2016).CAS 

Google Scholar 
Gosset, J. R. et al. A cross-species translational pharmacokinetic-pharmacodynamic evaluation of core body temperature reduction by the TRPM8 blocker PF-05105679. Eur. J. Pharm. Sci. 109S, S161–S167 (2017).
Google Scholar 
Winchester, W. J. et al. Inhibition of TRPM8 channels reduces pain in the cold pressor test in humans. J. Pharmacol. Exp. Ther. 351, 259–269 (2014).
Google Scholar 
Bianchi, B., Smith, P. A. & Abriel, H. The ion channel TRPM4 in murine experimental autoimmune encephalomyelitis and in a model of glutamate-induced neuronal degeneration. Mol. Brain 11, 1–10 (2018).
Google Scholar 
Li, K., Shi, Y., Gonye, E. C. & Bayliss, D. A. TRPM4 contributes to subthreshold membrane potential oscillations in multiple mouse pacemaker neurons. eNeuro 8, 1–13 (2021).
Google Scholar 
Dong, W. et al. Visual avoidance in Xenopus tadpoles is correlated with the maturation of visual responses in the optic tectum. J. Neurophysiol. 101, 803–815 (2009).
Google Scholar 
Bertolesi, G. E., Debnath, N., Atkinson-Leadbeater, K., Niedzwiecka, A. & McFarlane, S. Distinct type II opsins in the eye decode light properties for background adaptation and behavioural background preference. Mol. Ecol. 30, 6659–6676 (2021).CAS 

Google Scholar 
Viczian, A. S. & Zuber, M. E. A simple behavioral assay for testing visual function in xenopus laevis. J. Vis. Exp. 12, 51726 (2014).
Google Scholar 
Myers, B. R., Sigal, Y. M. & Julius, D. Evolution of thermal response properties in a cold-activated TRP channel. PLoS One 4, e5741 (2009).
Google Scholar 
Furman, B. L. S. et al. Pan-African phylogeography of a model organism, the African clawed frog ‘Xenopus laevis’. Mol. Ecol. 24, 909–925 (2015).CAS 

Google Scholar 
Wilson, R. S., James, R. S. & Johnston, I. A. Thermal acclimation of locomotor performance in tadpoles and adults of the aquatic frog Xenopus laevis. J. Comp. Physiol. B. 170, 117–124 (2000).CAS 

Google Scholar 
Kashiwagi, K. et al. Xenopus tropicalis: an ideal experimental animal in amphibia. Exp. Anim. 59, 395–405 (2010).CAS 

Google Scholar 
Martínez-Freiría, F., Toyama, K. S., Freitas, I. & Kaliontzopoulou, A. Thermal melanism explains macroevolutionary variation of dorsal pigmentation in Eurasian vipers. Sci. Rep. 10, 72871–1 (2020).Tanaka, K. Does the thermal advantage of melanism produce size differences in color-dimorphic snakes? Zool. Sci. 26, 698–703 (2009).
Google Scholar 
Moreno Azócar, D. L., Nayan, A. A., Perotti, M. G. & Cruz, F. B. How and when melanic coloration is an advantage for lizards: the case of three closely-related species of Liolaemus. Zool. (Jena.) 141, 125774 (2020).
Google Scholar 
Azócar, D. L. M. et al. Effect of body mass and melanism on heat balance in Liolaemus lizards of the goetschi clade. J. Exp. Biol. 219, 1162–1171 (2016).
Google Scholar 
Smith, K. R. et al. Colour change on different body regions provides thermal and signalling advantages in bearded dragon lizards. Proc. R. Soc. B Biol. Sci. 283, 20160626 (2016).
Google Scholar 
Rowe, J. W. et al. Thermal and substrate color-induced melanization in laboratory reared red-eared sliders (Trachemys scripta elegans). J. Therm. Biol. 61, 125–132 (2016).
Google Scholar 
Larsen, E. H. Dual skin functions in amphibian osmoregulation. Comp. Biochem. Physiol. A. Mol. Integr. Physiol. 253, 110869 (2021).CAS 

Google Scholar 
Franco-Belussi, L., Sköld, H. N. & De Oliveira, C. Internal pigment cells respond to external UV radiation in frogs. J. Exp. Biol. 219, 1378–1383 (2016).
Google Scholar 
Langhelle, A., Lindell, M. J. & Nyström, P. Effects of ultraviolet radiation on amphibian embryonic and larval development. J. Herpetol. 33, 449–456 (1999).
Google Scholar 
Mueller, K. P. & Neuhauss, S. C. F. Sunscreen for fish: co-option of UV light protection for camouflage. PLoS One 9, e87372 (2014).
Google Scholar 
Perotti, M. G., Diéguez, M. & Del, C. Effect of UV-B exposure on eggs and embryos of patagonian anurans and evidence of photoprotection. Chemosphere 65, 2063–2070 (2006).CAS 

Google Scholar 
Nilsson Sköld, H., Aspengren, S. & Wallin, M. Rapid color change in fish and amphibians – function, regulation, and emerging applications. Pigment Cell Melanoma Res. 26, 29–38 (2013).
Google Scholar 
Vences, M. et al. Field body temperatures and heating rates in a montane frog population: the importance of black dorsal pattern for thermoregulation on JSTOR. Ann. Zool. Fennici 39, 209–220 (2002).
Google Scholar 
Lindgren, J. et al. Skin pigmentation provides evidence of convergent melanism in extinct marine reptiles. Nature 506, 484–488 (2014).CAS 

Google Scholar 
Bonino, M. F., Cruz, F. B. & Perotti, M. G. Does temperature at local scale explain thermal biology patterns of temperate tadpoles? J. Therm. Biol. 94, 102744 (2020).
Google Scholar 
Kumar, S., Stecher, G., Li, M., Knyaz, C. & Tamura, K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549 (2018).CAS 

Google Scholar 
Liu, T. et al. RNA interference-mediated depletion of TRPM8 enhances the efficacy of epirubicin chemotherapy in prostate cancer LNCaP and PC3 cells. Oncol. Lett. 15, 4129–4136 (2018).
Google Scholar 
Kashina, A. S. et al. Protein Kinase A, which regulates intracellular transport, forms complexes with molecular motors on organelles. Curr. Biol. 14, 1877–1881 (2004).CAS 

Google Scholar  More

Daskalova, G. N. et al. Landscape-scale forest loss as a catalyst of population and biodiversity change. Science 368(6497), 1341–1347 (2020).ADS 
CAS 

Google Scholar 
Betts, M. G. et al. Extinction filters mediate the global effects of habitat fragmentation on animals. Science 366(6470), 1236–1239 (2019).ADS 
CAS 

Google Scholar 
Siddig, A. A., Ellison, A. M., Ochs, A., Villar-Leeman, C. & Lau, M. K. How do ecologists select and use indicator species to monitor ecological change? Insights from 14 years of publication in Ecological Indicators. Ecol. Ind. 60, 223–230 (2016).
Google Scholar 
Thancharoen, A. Well managed firefly tourism: A good tool for firefly conservation in Thailand. Lampyrid. 2, 142–148 (2012).
Google Scholar 
Hwang, Y. T., Moon, J., Lee, W. S., Kim, S. A. & Kim, J. Evaluation of firefly as a tourist attraction and resource using contingent valuation method based on a new environmental paradigm. J. Qual. Assur. Hosp. Tour. 21(3), 320–336 (2019).Carlson, A. D. & Copeland, J. Flash communication in fireflies. Q. Rev. Biol. 60(4), 415–436 (1985).
Google Scholar 
Evans, T. R., Salvatore, D., van de Pol, M. & Musters, C. J. M. Adult firefly abundance is linked to weather during the larval stage in the previous year. Ecol. Entomol. 44(2), 265–273 (2018).
Google Scholar 
Lewis, S. M. et al. A global perspective on firefly extinction threats. Bioscience 70(2), 157–167 (2020).
Google Scholar 
Cao, C. Q., Zhang, Y., Wang, Y. Z. & He, H. Progress in the research, protection, development and utilization of fireflies. J. Environ. Entomol.1–36 (2022).Myers, N., Mittermeier, R. A., Mittermeier, C. G., da Fonseca, G. A. B. & Kent, J. Biodiversity hotspots for conservation priorities. Nature 403(6772), 853–858 (2000).ADS 
CAS 

Google Scholar 
Thorn, J. S., Nijman, V., Smith, D. & Nekaris, K. A. I. Ecological niche modelling as a technique for assessing threats and setting conservation priorities for Asian slow lorises (Primates:Nycticebus). Divers. Distrib. 15(2), 289–298 (2009).
Google Scholar 
Elith, J. & Leathwick, J. R. Species distribution models: Ecological explanation and prediction across space and time. Annu. Rev. Ecol. Evol. Syst. 40(1), 677–697 (2009).
Google Scholar 
Phillips, S. J., Anderson, R. P. & Schapire, R. E. Maximum entropy modeling of species geographic distributions. Ecol. Model. 190(3–4), 231–259 (2006).
Google Scholar 
Hirzel, A. H., Hausser, J., Chessel, D. & Perrin, N. Ecological-Niche Factor Analysis: How to compute habitat-suitability maps without absence data?. Ecology 83(7), 2027–2036 (2002).
Google Scholar 
Nelder, J. A. & Wedderburn, R. W. Generalized linear models. J. R. Stat. Soc. Ser. A (General). 135(3), 370–384 (1972).
Google Scholar 
Hastie, T. J. Generalized additive models. Statistical models in S. Routledge. 249–307 (2017).Stockwell, D. R. & Noble, I. R. Induction of sets of rules from animal distribution data: A robust and informative method of data analysis. Math. Comput. Simul. 33(5–6), 385–390 (1992).
Google Scholar 
Beaumont, L. J., Hughes, L. & Poulsen, M. Predicting species distributions: use of climatic parameters in BIOCLIM and its impact on predictions of species’ current and future distributions. Ecol. Model. 186(2), 251–270 (2005).
Google Scholar 
Jung, J. M., Lee, W. H. & Jung, S. Insect distribution in response to climate change based on a model: Review of function and use of CLIMEX. Entomol. Res. 46(4), 223–235 (2016).
Google Scholar 
Phillips, S. J. & Dudík, M. Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation. Ecography 31(2), 161–175 (2008).
Google Scholar 
Moreno, R., Zamora, R., Molina, J. R., Vasquez, A. & Herrera, M. Á. Predictive modeling of microhabitats for endemic birds in South Chilean temperate forests using Maximum entropy (Maxent). Eco. Inform. 6(6), 364–370 (2011).
Google Scholar 
Wang, Z. et al. Prediction of potential distribution of the invasive Chrysanthemum Lace Bug, Corythucha marmorata in China based on Maxent. J. Environ. Entomol. 41(3), 626–633 (2019).
Google Scholar 
Li, A. et al. MaxEnt modeling to predict current and future distributions of Batocera lineolata (Coleoptera: Cerambycidae) under climate change in China. Ecoscience 27(1), 23–31 (2020).
Google Scholar 
Sutherland, L. N., Powell, G. S. & Bybee, S. M. Validating species distribution models to illuminate coastal fireflies in the South Pacific (Coleoptera: Lampyridae). Sci. Rep. 11(1), 1–12 (2021).ADS 

Google Scholar 
Fu, X. H., Ballantyne, L. A. & Lambkin, C. Emeia gen. nov., a new genus of Luciolinae fireflies from China (Coleoptera: Lampyridae) with an unusual trilobite-like larva, and a redescription of the genus Curtos Motschulsky. Zootaxa. 3403(1), 1–53 (2012).Idris, N. S. et al. The dynamics of landscape changes surrounding a firefly ecotourism area. Glob. Ecol. Conserv. 29, e01741 (2021).
Google Scholar 
Santiago-Blay, J. A. Silent Sparks: The Wondrous World of Fireflies. Life: The Excitement of Biology. (2016).Picchi, M. S., Avolio, L., Azzani, L., Brombin, O. & Camerini, G. Fireflies and land use in an urban landscape: the case of Luciola italica L.(Coleoptera: Lampyridae) in the city of Turin. J. Insect Conserv. 17(4), 797–805 (2013).Pearsons, K. A., Lower, S. E. & Tooker, J. F. Toxicity of clothianidin to common Eastern North American fireflies. PeerJ 9, e12495 (2021).
Google Scholar 
Madruga Rios, O. & Hernández Quinta, M. Larval Feeding Habits of the Cuban Endemic FireflyAlecton discoidalisLaporte (Coleoptera: Lampyridae). Psyche J. Entomol. 2010, 1–5 (2010).Roberge, J. M. & Angelstam, P. E. R. Usefulness of the umbrella species concept as a conservation tool. Conserv. Biol. 18(1), 76–85 (2004).
Google Scholar 
Bowen-Jones, E. & Entwistle, A. Identifying appropriate flagship species: The importance of culture and local contexts. Oryx 36(2), 189–195 (2002).
Google Scholar 
Walpole, M. J. & Leader-Williams, N. Tourism and flagship species in conservation. Biodivers. Conserv. 11(3), 543–547 (2002).Zhejiang Provincial Bureau of Statistics. Zhejiang physical geography profile, http://tjj.zj.gov.cn/col/col1525489/index.html (2022).Zhejiang Provincial Forestry Department. Announcement of Forest Resources and Their Ecological Function Value in Zhejiang Province. Zhejiang Daily. https://doi.org/10.38328/n.cnki.nzjrb.2016.002829 (2016).Boria, R. A., Olson, L. E., Goodman, S. M. & Anderson, R. P. Spatial filtering to reduce sampling bias can improve the performance of ecological niche models. Ecol. Model. 275, 73–77 (2014).
Google Scholar 
Brown, J. L. SDM toolbox: A python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. Methods Ecol. Evol. 5(7), 694–700 (2014).
Google Scholar 
Fick, S. E. & Hijmans, R. J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37(12), 4302–4315 (2017).
Google Scholar 
Elvidge, C. D., Zhizhin, M., Ghosh, T., Hsu, F.-C. & Taneja, J. Annual time series of global VIIRS nighttime lights derived from monthly averages: 2012 to 2019. Remote Sens. 13(5), 922 (2021).ADS 

Google Scholar 
WAN, J. et al. Predicting the potential geographic distribution of Bactrocera bryoniae and Bactrocera neohumeralis (Diptera: Tephritidae) in China using MaxEnt ecological niche modeling. J. Integr. Agric. 19(8), 2072–2082 (2020).Zhou, R. et al. Projecting the potential distribution of glossina morsitans (Diptera: Glossinidae) under climate change using the MaxEnt model. Biology. 10(11), 1150 (2021).
Google Scholar 
Hill, M. P., Hoffmann, A. A., McColl, S. A. & Umina, P. A. Distribution of cryptic blue oat mite species in Australia: current and future climate conditions. Agric. For. Entomol. 14(2), 127–137 (2011).
Google Scholar 
Su, H., Bista, M. & Li, M. Mapping habitat suitability for Asiatic black bear and red panda in Makalu Barun National Park of Nepal from Maxent and GARP models. Sci. Rep. 11(1), 1 (2021).ADS 
CAS 

Google Scholar 
Proosdij, A. J., Sosef, M., Wieringa, J. & Raes, N. Minimum required number of specimen records to develop accurate species distribution models. Ecography 39, 542–552 (2016).
Google Scholar 
Lobo, J. M., Jiménez-Valverde, A. & Real, R. AUC: A misleading measure of the performance of predictive distribution models. Glob. Ecol. Biogeogr. 17(2), 145–151 (2008).
Google Scholar 
Allouche, O., Tsoar, A. & Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 43(6), 1223–1232 (2006).
Google Scholar 
Liu, C., Newell, G. & White, M. On the selection of thresholds for predicting species occurrence with presence-only data. Ecol. Evol. 6(1), 337–348 (2016).
Google Scholar 
Swets, J. A. Measuring the accuracy of diagnostic systems. Science 240(4857), 1285–1293 (1988).ADS 
CAS 
MATH 

Google Scholar 
Pearce, J. & Ferrier, S. Evaluating the predictive performance of habitat models developed using logistic regression. Ecol. Model. 133(3), 225–245 (2000).
Google Scholar 
Gama, M., Crespo, D., Dolbeth, M. & Anastácio, P. M. Ensemble forecasting of Corbicula fluminea worldwide distribution: projections of the impact of climate change. Aquat. Conserv. Mar. Freshwat. Ecosyst. 27(3), 675–684 (2017).
Google Scholar 
Zhao, Y., Deng, X., Xiang, W., Chen, L. & Ouyang, S. Predicting potential suitable habitats of Chinese fir under current and future climatic scenarios based on Maxent model. Eco. Inform. 64, 101393 (2021).
Google Scholar 
Evans, T. R., Salvatore, D., van de Pol, M. & Musters, C. J. M. Adult firefly abundance is linked to weather during the larval stage in the previous year. Ecol. Entomol. 44(2), 265–273 (2018).Chettri, B., Bhupathy, S. & Acharya, B. K. Distribution pattern of reptiles along an eastern Himalayan elevation gradient India. Acta Oecol. 36(1), 16–22 (2010).ADS 

Google Scholar 
Brown, J. H. Mammals on mountainsides: elevational patterns of diversity. Global Ecol. Biogeogr. 10(1), 101–109 (2001).Gairola, S., Sharma, C. M., Ghildiyal, S. K. & Suyal, S. Tree species composition and diversity along an altitudinal gradient in moist tropical montane valley slopes of the Garhwal Himalaya India. For. Sci. Technol. 7(3), 91–102 (2011).
Google Scholar 
Pearson, R. G., Raxworthy, C. J., Nakamura, M. & Townsend Peterson, A. Predicting species distributions from small numbers of occurrence records: A test case using cryptic geckos in Madagascar. J. Biogeogr. 34(1), 102–117 (2007).Hernandez, P. A., Graham, C. H., Master, L. L. & Albert, D. L. The effect of sample size and species characteristics on performance of different species distribution modeling methods. Ecography 29(5), 773–785 (2006).
Google Scholar 
Abe, N. Kansei estimation on luminescence of Firefly-Kansei information measurement and welfare utilization. J. Japan Soc. Kansei Eng. 3(2), 41–50 (2004).
Google Scholar 
Buckley, R. et al. Economic value of protected areas via visitor mental health. Nat. Commun. 10(1), 1 (2019).
Google Scholar 
Lewis, S. M. et al. Firefly tourism: Advancing a global phenomenon toward a brighter future. Conserv. Sci. Pract. 3(5), 1 (2021).
Google Scholar  More