LIFEPLAN: A worldwide biodiversity sampling design

LAB PROTOCOL LIFEPLAN: A worldwide biodiversity sampling design Bess HardwickID1☯, Deirdre KerdraonID2☯*, Hanna M. K. Rogers2, Dimby Raharinjanahary3, Eric Tsiriniaina Rajoelison3, Tommi MononenID1, Petteri Lehikoinen4, Gaia Banelyte2, Arielle Farrell2, Brian L. Fisher3,5, Tomas Roslin1,2, Otso Ovaskainen1,6 1 Organismal and Evolutionary Biology Research Programme, University of Helsinki, Helsinki, Finland, 2 Department of Ecology, Swedish University of Agricultural Sciences, Uppsala, Sweden, 3 Madagascar Biodiversity Center, Antananarivo, Madagascar, 4 The Finnish Museum of Natural History, University of Helsinki, Helsinki, Finland, 5 California Academy of Sciences, San Francisco, CA, United States of America, 6 Department of Biological and Environmental Science, University of Jyva¨skyla¨, Jyva¨skyla¨ , Finland ☯ These authors contributed equally to this work. * deirdre.kerdraon@slu.se Abstract As the technology for mass identification of species is advancing rapidly, we developed a field sampling method that takes advantage of the emerging possibilities of combining sen- sor-based data with automated high-throughput data processing. This article describes the five field sampling methods used by the LIFEPLAN project to collect biodiversity data in a systematic manner, all over the world. These methods are designed for use by anyone with basic biology or ecology knowledge from the higher education or university level. We pres- ent the selection and characteristics of international sampling locations for urban and natural

sites, as well as the nested scale design in the Nordic countries and Madagascar. We describe the steps to collect sequences of animal images (.jpg) from infrared triggered cam- era traps, audio data (.WAV) of environment sounds from audio recorders, invertebrate samples in ethanol from Malaise traps for DNA metabarcoding, as well as both soil samples and 24-hour air samples obtained from cyclone samplers for fungal DNA metabarcoding. To ensure the usability and consistency of the data for future use, we pay particular attention to the metadata collected. In specifying the current sampling protocols, we note that technol- ogy will continue to improve and evolve. Hardware will also change within a short time period, with the advantage of improving the equipment used for collecting samples. Thus, we present examples of the samples collected by each current sampling method, to be used as a baseline or in comparison with different equipment models. Introduction We have so far only discovered an estimated 20% of all species on planet Earth [1]; at the same time, species are going extinct at an alarming rate [2]. If we are to understand current species diversity and predict future scenarios in a rapidly changing environment, monitoring PLOS ONE PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 1 / 15 a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 OPEN ACCESS Citation: Hardwick B, Kerdraon D, Rogers HMK, Raharinjanahary D, Rajoelison ET, Mononen T, et al. (2024) LIFEPLAN: A worldwide biodiversity sampling design. PLoS ONE 19(12): e0313353.

https://doi.org/10.1371/journal.pone.0313353 Editor: Bernd Schierwater, Tierarztliche Hochschule Hannover, GERMANY Received: March 11, 2024 Accepted: October 22, 2024 Published: December 31, 2024 Peer Review History: PLOS recognizes the benefits of transparency in the peer review process; therefore, we enable the publication of all of the content of peer review and author responses alongside final, published articles. The editorial history of this article is available here: https://doi.org/10.1371/journal.pone.0313353 Copyright: © 2024 Hardwick et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability Statement: No datasets were generated or analyzed during the current study. Unrestricted data from this study will be made available upon study completion.

biodiversity over extended time-periods, across a wide variety of different habitats, and in all major species groups, is essential. To establish a baseline against which to gauge future change, this monitoring needs to start now. Current biodiversity surveys are often limited in spatial coverage and timespan, biased towards regions with the most people and financial resources, and restricted to species that are easiest to see and identify [3–6]. As experts using morphological identification are commonly specialists in one or a few groups, there are simply not enough experts to produce biodiversity data with truly global coverage of the most diverse species groups [7]. In recent years, the development of molecular identification and machine learning methods, as well as increases in computing power, have enabled us to teach algorithms to identify species [8–10]. The sheer number of samples that can be identified by a computer is thus far greater than what human experts could ever process. To take full advantage of these new methods, we have devised sampling designs that gener- ate consistent data. In the LIFEPLAN project, sampling started in 2021 and is to continue until the end of 2025. A central aim was to create a standardised sampling scheme for monitoring terrestrial animal and fungal biodiversity under any climatic conditions. These methods require minimal technical or taxonomic expertise, allowing a wide range of participants to col- lect data. While many national long term monitoring endeavours already exist [11–13], LIFE-

PLAN is the first to create a worldwide sampling design to systematically collect data to target a broad range of taxonomic groups, which can be used for monitoring biodiversity throughout the year. By combining these sampling methods on a large scale with automated species identi- fication and novel statistical methods, LIFEPLAN aims to generate a predictive understanding of global biodiversity. These methods are replicable in time and space, which allows any future dataset that would be collected by identical methods to be compared to the original Lifeplan data. While the technology will improve, we will still be able to infer changes in biodiversity by matching the sampling design. The sampling design of LIFEPLAN is aimed at covering major terrestrial taxa excluding plants; thus, sampling is aimed at arthropods, fungi and vertebrates. Arthropods are collected with passive Malaise traps. To learn both what fungi are dispersing in the area as well as which species are actually growing in the soil, we collect fungal spores from the air with cyclone sam- plers as well as soil fungi from soil samples. For large or audible groups such as mammals and birds, we use camera traps and audio recorders, which are non-invasive and low-cost. Even these semi-automated methods do require regular personnel to collect the samples every week year-round. The personnel must be able to negotiate the common challenges of field work in their local conditions and have basic computer skills. Some financial resources

are also required to purchase the specialised equipment and consumables, and to fund the DNA sequencing. Having all these elements in place at a site that is also suitable for sampling has limited global LIFEPLAN coverage, as such sites are much easier to find in regions with more financial resources. The specific makes and models of equipment used have had to be cost-effective, but are challenged by some climatic conditions. Humidity is causing breakage in camera traps and AudioMoths, and cyclone samplers are struggling in freezing temperatures. International ship- ping of physical samples to a DNA sequencing facility requires special permits, and transfer- ring the large amounts of image and audio data collected requires a stable and reasonably fast internet connection. The methods for identifying species from these samples are still developing and improving. One limitation is the availability of annotated reference data [11, 14–21] that connects a DNA barcode sequence, image or sound to a known species. We are unable to identify all the species collected, and many of the species that we collect are novel and require taxonomical work to PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 2 / 15 Funding: OO, TR: grant number 856506. European Research Council (ERC), https://erc.europa.eu/ homepage. The funder did not play any role in the study design, data collection and analysis, decision to publish or preparation of the manuscript. The work of all the authors was supported by the same

grant awarded to OO and TR. Competing interests: The authors have declared that no competing interests exist.

even describe them. Yet, now is the time to collect an archive of samples, as the species we are looking at might be lost by the time we can identify them. This is not to say that species devoid of a traditional scientific binomial could not contribute to present-day analyses. Of the species detected but not yet assigned a traditional scientific name, we can form species hypotheses and thus efficiently estimate both community-level attributes such as species numbers, as well as spatial and temporal turnover in community composition [22–25]. We can also characterise phylogenetic and functional diversity at the community level [26–30]. And for the interim taxa delimited, we can characterise distributions and their drivers [31]. In terms of the focal ecosystems of LIFEPLAN, we have clearly targeted terrestrial biodiver- sity. Since any realistic sampling scheme will have to select some targets at the expense of oth- ers, we have chosen not to sample plants, or ground or tree crown-dwelling invertebrates. Furthermore, the LIFEPLAN project is not yet using the data collected to its full capacity due to funding and time constraints. More data can be acquired from the collected samples as tech- nology improves. Spore samples are likely also to include pollen, the DNA of which could be sequenced as well, and soil samples will also contain much more than fungi. More groups of vocalising taxa can be identified from audio recordings (such as amphibians and insects).The sampling sites can also be fitted with further sampling devices, such as sampling for aerosol-

borne DNA [32]. Current initiatives include adding site-specific measurements of ecosystem processes. Materials and methods The protocols described in this peer-reviewed article are published on protocols.io, updated March 6, 2024, and are included for printing as S1–S5 Files with this article. Camera. The protocol described in this peer-reviewed article is published on protocols.io, dx.doi.org/10.17504/protocols.io.q26g7pxp1gwz/v3 and is included for printing as S1 File with this article. Audio. The protocol described in this peer-reviewed article is published on protocols.io, dx.doi.org/10.17504/protocols.io.kqdg3xbp1g25/v2 and is included for printing as S2 File with this article. Cyclone sampler. The protocol described in this peer-reviewed article is published on pro- tocols.io, dx.doi.org/10.17504/protocols.io.6qpvr3zn2vmk/v1 and is included for printing as S3 File with this article. Malaise trap. The protocol described in this peer-reviewed article is published on proto- cols.io, dx.doi.org/10.17504/protocols.io.kqdg3xkdqg25/v2 and is included for printing as S4 File with this article. Soil sampling. The protocol described in this peer-reviewed article is published on proto- cols.io, dx.doi.org/10.17504/protocols.io.5jyl8pw9rg2w/v2 and is included for printing as S5 File with this article. Site selection Since biodiversity patterns may be differentially impacted by different community assembly processes at different spatial scales [33, 34], we opted for a hierarchical sampling design, span- ning six orders of magnitude in spatial scales (from 10−1 km to 104 km).

The Global (104 km scale) sampling consisted of 73 sampling stations spread across all con- tinents. Sampling at the National (103 km) scale was based on 43 sampling stations, and at the Nested (102–10−1 km) scale on another 43 sampling stations. Sampling at the National and Nested scales focused on two contrasting locations: the Nordic countries (Finland, Norway, Sweden) and Madagascar. Given their widely different patterns of climatic conditions and PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 3 / 15

regional and local endemism [35–37], this regional contrast yields insight into the context- dependency of the results achieved–as then validated by global patterns. We chose the Nordic countries and Madagascar for National and Nested scale sampling because they are of comparable area but have very different evolutionary histories. The Nordic countries were colonised after the most recent ice age, while in Madagascar, species have been evolving in isolation for millions of years [35]. Nordic species are very well known, while the enormous biodiversity of Madagascar is still poorly understood [35, 37–39]. The Swedish insect fauna consists of postglacial immigrants and is one of the best known in the world. Mad- agascar is the world’s most important biodiversity hotspot, with a long period of evolution in isolation–as combined with recent, massive deforestation [40]. Many species are critically endangered if not already extinct [41, 42]. By comparison, endemic species in the Nordic countries are next to absent. Given these contrasts, we aim to resolve–among other things— differences in range size between Nordic and Malagasy taxa, patterns of turnover in space and time, and differences in phylogenetic community structure between the respective regions. Beyond providing a general understanding of basic community assembly processes, each of these insights is needed for improved conservation prioritization—thus attesting to the value of large-scale biodiversity survey and biomonitoring programmes. To resolve patterns of community structure and turnover in space and time, we have

invented a hierarchically structured sampling design in Sweden and Madagascar. At the National scale in Madagascar, we have selected 25 sites spanning all major forest types. In the Nordic countries, we selected 18 sites based on their contribution to spatially even coverage of the region. For the Nested scale, we chose an area in southern Sweden and another in Mada- gascar that were accessible enough to be covered by a single sampling team, were covered by forest, and where landowners gave permission to sample. Overall, we will thus achieve site-by-time-by-species tables of fungi as airborne spores and mycelia in soil, insects, birds and mammals, across hierarchically Nested scales. Global design: Natural vs. urban The largest spatial scale in our study design is the Global scale. We initially invited researchers from all over the world to submit proposals for study sites. Of 184 sites proposed for the Global design, we selected 87, of which 14 eventually had to drop out (Fig 1). In the Global design, sites form natural/urban pairs, where the urban location is in or immediately adjacent to a built-up area, and the natural location is in an undisturbed setting with little direct human influence. In other respects, such as altitude and vegetation type, we have aimed to have the two locations be as similar as possible, and separated by 20–50 kilo- metres (Fig 2). We have tried to cover the world map as evenly as possible, and within the selected regions we have asked teams to sample a habitat that is typical for their region. As a

result, we are currently sampling in a wide range of habitats from various forest types to savan- nah, wetlands, desert and tundra. Teams alternate yearly between their natural and urban loca- tions, so that only one location of each pair is sampled at one time. To avoid confounding the effects of urbanisation and sampling year, we have randomised the design so that half of the teams started in the urban location and the other half started in the natural location. In practice, levels of urbanisation vary in different parts of the world so that a “natural” location in a densely populated area such as central Europe may be much more urban than an “urban” location in a sparsely populated area such as Canada. There are many different mea- sures of urbanisation, but using e.g. the Human Footprint Index 2009 [43], we see that for most, but not all, site pairs the “urban” location scores higher than the “natural” (Fig 3B). “Nat- ural” locations also often ended up at higher altitudes, as lower-altitude areas were fully built PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 4 / 15

up in many parts of the world (Fig 3A). Other practical considerations used to restrict site selection included accessibility, security, financial constraints and land ownership. National design The National scale was the intermediate between the Global and Nested scales. At this scale we do not alternate between a natural and urban location, and have instead chosen natural type Fig 1. Locations of LIFEPLAN sites around the world. Global sites are shown as black dots. The dense clusters of points without fill in Madagascar and the Nordic countries are the sites forming the National and Nested designs. https://doi.org/10.1371/journal.pone.0313353.g001 Fig 2. Histogram of actual distances within site pairs. In choosing natural / urban site pairs, we aimed for a distance of 20–50 km within pairs. Finding a natural / urban site pair 20–50 km apart that was also accessible and otherwise comparable was not always possible, and we accepted some site pairs that fell outside this range if the locations were otherwise valuable in terms of global coverage. The actual distribution of distances within the finally selected site pairs is shown here. https://doi.org/10.1371/journal.pone.0313353.g002 PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 5 / 15

sites that are relatively undisturbed and representative of the local environment. At the National level, we sample in Madagascar as well as the Nordic countries of Sweden, Norway and Finland (Fig 4). Nested design To study hierarchical patterns of community structure in space, we have created a Nested sam- pling design with increasing distance between sites (Fig 5). We originally searched for sites in a Fig 3. Comparisons within site pairs. (A) Differences in elevation (metres above sea level) within natural/urban site pairs. (B) Differences in Human Footprint Index within natural/urban site pairs. https://doi.org/10.1371/journal.pone.0313353.g003 Fig 4. Locations of sites in the National design. Black points show National design sites in Madagascar (left) and the Nordic countries (right). https://doi.org/10.1371/journal.pone.0313353.g004 PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 6 / 15

Nested design that were on the order of 50 km, 5 km and 500 m apart. In Madagascar espe- cially, the search was constrained by the scarcity of remaining forest fragments and the high- est-level distances ended up being approximately 70–120 km. In Sweden we achieved highest- level distances of approximately 60–70 km. In Sweden, we implemented a triangular Nested design (n = 3) with 3 clusters each containing 5 sites, and in Madagascar, we had a square design (n = 4) with 4 clusters each containing 7 sites. The triangular design in used in Sweden was reduced from an originally square design due to the high costs of sampling in Sweden. We note that the original square design includes higher coverage of a wider range of pairwise dis- tances, simply by virtue of having more sampling points. Thus, this design was prioritized in Madagascar, as Madagascar has higher beta diversity [37]. If financially feasible, we recom- mend the square design especially in regions of high beta diversity. Site design The two site designs accommodate different budgets (Fig 6). The extensive design covers a one-hectare square: five cameras and audio recorders are placed in each corner and the centre. Soil samples are collected from each of the same five points. In the simplified design, three cameras and audio recorders are placed along a transect of 140 m, with soil samples collected at each of the three points. The reduced design corresponds to the extensive design with two diagonally opposite corners removed. We started with the extended design in October 2020 and switched to the reduced design in

May 2023. This switch was caused by financial constraints, due to equipment breaking sooner than expected. Fig 5. Locations of sites in the Nested design. (A) Schematic of the Nested scale field design in Madagascar (top) and Sweden (bottom). (B) Maps showing the actual but approximate site locations for the Nested design in Madagascar (top) and Sweden (bottom). https://doi.org/10.1371/journal.pone.0313353.g005 PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 7 / 15

Sampling methods Cameras and audio recorders collected data continuously year-round with weekly visits needed to change batteries and collect the data manually. Camera. The protocol described in this peer-reviewed article is published on protocols.io, dx.doi.org/10.17504/protocols.io.q26g7pxp1gwz/v3 and is included for printing as S1 File with this article. Audio. The protocol described in this peer-reviewed article is published on protocols.io, dx.doi.org/10.17504/protocols.io.kqdg3xbp1g25/v2 and is included for printing as S2 File with this article. Two 24-hour cyclone samples are collected while weather conditions allow the sampler to work, by visiting the site on three consecutive days each week or by using a timer for one of the samples and visiting the site two consecutive days per week. In a few cases, only one sample is collected using a timer due to the difficulty of reaching the site multiple days per week. The protocol described in this peer-reviewed article is published on protocols.io, dx.doi.org/10. 17504/protocols.io.6qpvr3zn2vmk/v1 and is included for printing as S3 File with this article. Flying arthropod samples are collected from Malaise traps once per week during the insect flying season producing a seven-day sample. The protocol described in this peer-reviewed arti- cle is published on protocols.io, dx.doi.org/10.17504/protocols.io.kqdg3xkdqg25/v2 and is included for printing as S4 File with this article. One soil sample from each of three or five points defined by the camera and audio recorder

locations is taken eight or four times per year when the temperatures are above freezing. The original frequency of eight times a year was reduced to four times a year in January 2023. The protocol described in this peer-reviewed article is published on protocols.io, dx.doi.org/10. 17504/protocols.io.5jyl8pw9rg2w/v2 and is included for printing as S5 File with this article. Fig 6. Two alternate LIFEPLAN site designs. (A) Design of LIFEPLAN one-hectare site with five sampling points. (B) Design of LIFEPLAN transect with three sampling points. Cameras and AudioMoths are placed within 10 m of each sampling point. Soil samples are collected in an area within approximately 5 m of the camera and audio recorder. The Malaise trap and cyclone sampler can be placed anywhere within the site, preferably close to the centre or middle but not in the camera’s field of view. https://doi.org/10.1371/journal.pone.0313353.g006 PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 8 / 15

Metadata When collecting spatiotemporal data, it is essential to keep track of where and when each sam- ple was collected. At a minimum, each sample must be associated with information on the date and time of sample start and end, and the sampling location. We give each sample a unique identifier, which is a six-character code. A “sample” is either: a bottle containing a Malaise sample; a vial containing a cyclone sample; a bag containing a soil sample; or, for image and audio data, the contents of a memory card. For image and audio data, the uniquely identified memory cards are reused week after week and a sample is defined as one batch of data collected between memory card changes, so the unique identifier of an image or audio sample consisted of the six-character code of the memory card combined with the date on which the sample is collected. We also give unique six-character codes to the sampling equipment. All equipment and samples are labelled with durable stickers that contain the unique identifier both as a human readable code and a scannable QR code. When placing or collecting a sample, we use the LIFEPLAN mobile application [44] to scan the QR code of the equipment and the QR code of the sample. In the case of soil samples, we scanned only the sample code as there is no fixed equipment. This scanning creates an Activ- ity, which is a row in a metadata table (Table 1). Table 1. Metadata fields recorded for each activity. Field Description DarwinCore term Activity ID Running number N/A

Created at Date and time that the Activity was created. N/A Updated at Date and time that the Activity was last edited. N/A Deleted at Date and time that the Activity was deleted. All deletions are soft deletions, where the Activity remains in the database but is marked deleted. N/A Link Link to audio or image data: here users can paste a share link which gives access to the actual audio or image data stored in a cloud. After the data have been transferred to LIFEPLAN, the link becomes a link to the LIFEPLAN data repository. N/A Processing status For image and audio data, this field indicates whether or not the data have been transferred to LIFEPLAN. N/A https://doi.org/10.1371/journal.pone.0313353.t001 PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 9 / 15

For cases where the mobile application fails or metadata needs to be corrected, we also have a browser version: LIFEPLAN Web Admin [45]. Because metadata can be edited after collect- ing or created after the fact, we also record the fields described in Table 1. For each sample, we thus hold spatiotemporal metadata of where and when it was placed in the field and collected. Our database has additional tables with information on the Sites, Teams, Equipment, Sample codes and Users, which we can associate with the Activity by the relevant IDs. Ethics declarations The acquisition and transport of biological samples is strictly regulated by international agree- ments. All sampling and processing of samples within LIFEPLAN adheres to both national and international legislation and agreements. To comply with the Convention on Biological Diversity, with each team contributing to LIFEPLAN we have a Material Transfer Agreement (MTA). With each nation, we have estab- lished via the ABS clearing house whether Prior Informed Consent (PIC) is required, and acquired it where necessary. The MTA and PIC define Access and Benefit Sharing clauses to provide non-monetary benefits to the data providers. This study was subject to review by European Research Council Executive Agency (ERCEA) Unit B1 –Ethics Review and Expert Management. Expected results Samples for DNA based identification Physical samples can be identified with metabarcoding, using different primers depending on the species group of interest. Malaise samples can also be individually sequenced or barcoded

and photographed. In LIFEPLAN, we metabarcode all physical samples for the COI (cyto- chrome oxidase I) region for Malaise samples or ITS2 (internal transcribed spacer 2) region for cyclone and soil samples, as well as individually COI-barcoding and bulk imaging a subset of Malaise samples. Cyclone samples are collected in 1.5 ml sterile microtubes, then freeze dried and stored at -20˚C. They usually contain a small amount of fungal spores, pollen and other particles under 1 μm in diameter. Samples may contain small insects or water, which are removed and their presence noted before analysis (Fig 7A; Table 2). Malaise samples are collected in Nalgene bottles filled with pure ethanol to euthanise and conserve the arthropods, and stored at -20 C˚ (Fig 7B). Fig 7. Physical samples collected in project LIFEPLAN. (A) Two cleaned and freeze-dried cyclone samples; an insect was removed from the tube on the left. (B) A Malaise sample collected after one week in the field. (C) A soil sample inside a glassine bag. https://doi.org/10.1371/journal.pone.0313353.g007 PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 10 / 15

Soil samples are collected in 75 x 102 mm glassine bags (Fig 7C). They contain a compound sample of about three tablespoons of soil. Samples are dried immediately after collection, then freeze dried and stored at -20˚ C until analysis. Table 2. Example of a processing metadata table for cyclone samples. Sample Team Site Placement Collection Duration Notes Equipment condition Sample field condition Insect Water Freeze- dried CWCBAW SWE SH11 7.10.2021 9:00 8.10.2021 8:58 0d 23h 58 min battery died ok X X C4LY21 SWE UA21 11.10.2021 8:55 12.10.2021 08:56 1d 0h 1min ok insect X X C6T1SX AUS AU4_Urban 11.10.2021 3:10 12.10.2021 3:11 1d 0h 1min ok ok X https://doi.org/10.1371/journal.pone.0313353.t002 Fig 8. Example of photo series from LIFEPLAN cameras. The cameras record a five-image burst when triggered. (A) In daylight. (B) In darkness. https://doi.org/10.1371/journal.pone.0313353.g008 PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 11 / 15

Samples for automated species identification Image data is in JPG format, in sequences of five images taken in a burst, approximately 1 sec- ond apart. The camera is triggered by movement, and images are taken both during the day and at night (Fig 8). Audio recording data consist of regular one-minute segments, 48-hour continuous record- ings, and opportunistic recordings in between regular segments that are triggered by bat calls. Audio recordings are converted to WAV format for analysis. Animal species can be recognised from spectrogram images of sounds by machine learning algorithms (Fig 9; [9]), as can camera trap images The LIFEPLAN automatic species identification software will be published else- where and made openly accessible. Supporting information S1 File. Camera trapping protocol. (PDF) S2 File. Audio recording protocol. (PDF) S3 File. Cyclone sampling protocol. (PDF) S4 File. Malaise trapping protocol. (PDF) S5 File. Soil sampling protocol. (PDF) Acknowledgments Open Acoustic Devices provided the published version of the Lifeplan AudioMoth firmware. Pekka Niittynen calculated HFI and elevations for sites. Author Contributions Conceptualization: Tomas Roslin, Otso Ovaskainen. Data curation: Bess Hardwick, Deirdre Kerdraon. Funding acquisition: Tomas Roslin, Otso Ovaskainen. Fig 9. Audio recording converted into a spectrogram for analysis. https://doi.org/10.1371/journal.pone.0313353.g009 PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 12 / 15

Methodology: Bess Hardwick, Deirdre Kerdraon, Hanna M. K. Rogers, Dimby Raharinjanah- ary, Eric Tsiriniaina Rajoelison, Tommi Mononen, Petteri Lehikoinen, Gaia Banelyte, Arielle Farrell. Project administration: Bess Hardwick, Deirdre Kerdraon, Brian L. Fisher. Supervision: Brian L. Fisher, Tomas Roslin, Otso Ovaskainen. Writing – original draft: Bess Hardwick, Deirdre Kerdraon. Writing – review & editing: Hanna M. K. Rogers, Dimby Raharinjanahary, Eric Tsiriniaina Rajoelison, Tommi Mononen, Petteri Lehikoinen, Gaia Banelyte, Arielle Farrell, Brian L. Fisher, Tomas Roslin, Otso Ovaskainen. References 1. Mora C, Tittensor DP, Adl S, Simpson AGB, Worm B. How Many Species Are There on Earth and in the Ocean? Mace GM, editor. PLoS Biol. 2011; 9: e1001127. https://doi.org/10.1371/journal.pbio.1001127 PMID: 21886479 2. Ceballos G, Ehrlich PR, Barnosky AD, Garcı´a A, Pringle RM, Palmer TM. Accelerated modern human —induced species losses: Entering the sixth mass extinction. Sci Adv. 2015; 1: e1400253. https://doi. org/10.1126/sciadv.1400253 PMID: 26601195 3. Beck J, Bo¨ ller M, Erhardt A, Schwanghart W. Spatial bias in the GBIF database and its effect on model- ing species’ geographic distributions. Ecol Inform. 2014; 19: 10–15. https://doi.org/10.1016/j.ecoinf. 2013.11.002 4. Boakes EH, McGowan PJK, Fuller RA, Chang-qing D, Clark NE, O’Connor K, et al. Distorted Views of Biodiversity: Spatial and Temporal Bias in Species Occurrence Data. PLoS Biol. 2010; 8: e1000385. https://doi.org/10.1371/journal.pbio.1000385 PMID: 20532234

5. Troudet J, Grandcolas P, Blin A, Vignes-Lebbe R, Legendre F. Taxonomic bias in biodiversity data and societal preferences. Sci Rep. 2017; 7: 9132. https://doi.org/10.1038/s41598-017-09084-6 PMID: 28831097 6. Yesson C, Brewer PW, Sutton T, Caithness N, Pahwa JS, Burgess M, et al. How Global Is the Global Biodiversity Information Facility? Beach J, editor. PLoS ONE. 2007; 2: e1124. https://doi.org/10.1371/ journal.pone.0001124 PMID: 17987112 7. Engel MS, Cerı´aco LMP, Daniel GM, Dellape´ PM, Lo¨ bl I, Marinov M, et al. The taxonomic impediment: a shortage of taxonomists, not the lack of technical approaches. Zool J Linn Soc. 2021; 193: 381–387. https://doi.org/10.1093/zoolinnean/zlab072 8. Abarenkov K, Somervuo P, Nilsson RH, Kirk PM, Huotari T, Abrego N, et al. PROTAX -fungi: a web- based tool for probabilistic taxonomic placement of fungal internal transcribed spacer sequences. New Phytol. 2018; 220: 517–525. https://doi.org/10.1111/nph.15301 PMID: 30035303 9. Lauha P, Somervuo P, Lehikoinen P, Geres L, Richter T, Seibold S, et al. Domain-specific neural net- works improve automated bird sound recognition already with small amount of local data. Methods Ecol Evol. 2022; 13: 2799–2810. https://doi.org/10.1111/2041-210X.14003 10. Roslin T, Somervuo P, Pentinsaari M, Hebert PDN, Agda J, Ahlroth P, et al. A molecular-based identifi- cation resource for the arthropods of Finland. Mol Ecol Resour. 2022; 22: 803–822. https://doi.org/10. 1111/1755-0998.13510 PMID: 34562055

11. Geiger M, Moriniere J, Hausmann A, Haszprunar G, Wa¨ gele W, Hebert P, et al. Testing the Global Mal- aise Trap Program—How well does the current barcode reference library identify flying insects in Ger- many? Biodivers Data J. 2016; 4: e10671. https://doi.org/10.3897/BDJ.4.e10671 PMID: 27932930 12. Karlsson D, Forshage M, Holston K, Ronquist F. The data of the Swedish Malaise Trap Project, a coun- trywide inventory of Sweden’s insect fauna. Biodivers Data J. 2020; 8: e56286. https://doi.org/10.3897/ BDJ.8.e56286 PMID: 33177946 13. Vaughan H, Brydges T, Fenech A, Lumb A. Monitoring Long-Term Ecological Changes Through the Ecological Monitoring and Assessment Network: Science-Based and Policy Relevant. Environ Monit Assess. 2001; 67: 3–28. https://doi.org/10.1023/a:1006423432114 PMID: 11339704 14. A Community for Naturalists � iNaturalist. [cited 5 Feb 2024]. https://www.inaturalist.org/. 15. Home | Wildlife Insights. [cited 5 Feb 2024]. https://www.wildlifeinsights.org/. 16. Home—LILA BC. [cited 5 Feb 2024]. https://lila.science/. PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 13 / 15

17. Cornell Lab of Ornithology. Macaulay Library at the Cornell Lab of Ornithology. https://www. macaulaylibrary.org/. 18. deWaard JR, Ratnasingham S, Zakharov EV, Borisenko AV, Steinke D, Telfer AC, et al. A reference library for Canadian invertebrates with 1.5 million barcodes, voucher specimens, and DNA samples. Sci Data. 2019; 6: 308. https://doi.org/10.1038/s41597-019-0320-2 PMID: 31811161 19. Kwong S, Srivathsan A, Meier R. An update on DNA barcoding: low species coverage and numerous unidentified sequences. Cladistics. 2012; 28: 639–644. https://doi.org/10.1111/j.1096-0031.2012. 00408.x PMID: 34856738 20. Museum fu¨r Naturkunde Berlin. Tierstimmenarchiv—Museum fu¨ r Naturkunde Berlin. https://www. tierstimmenarchiv.de/. 21. Xeno-canto foundation. Xeno-canto. https://www.xeno-canto.org. 22. Hartop E, Srivathsan A, Ronquist F, Meier R. Towards Large-Scale Integrative Taxonomy (LIT): Resolving the Data Conundrum for Dark Taxa. Dayrat B, editor. Syst Biol. 2022; 71: 1404–1422. https:// doi.org/10.1093/sysbio/syac033 PMID: 35556139 23. Hebert PDN, Ratnasingham S, Zakharov EV, Telfer AC, Levesque-Beaudin V, Milton MA, et al. Count- ing animal species with DNA barcodes: Canadian insects. Philos Trans R Soc B Biol Sci. 2016; 371: 20150333. https://doi.org/10.1098/rstb.2015.0333 PMID: 27481785 24. Srivathsan A, Ang Y, Heraty JM, Hwang WS, Jusoh WFA, Kutty SN, et al. Convergence of dominance and neglect in flying insect diversity. Nat Ecol Evol. 2023; 7: 1012–1021. https://doi.org/10.1038/ s41559-023-02066-0 PMID: 37202502

25. Srivathsan A, Hartop E, Puniamoorthy J, Lee WT, Kutty SN, Kurina O, et al. Rapid, large-scale species discovery in hyperdiverse taxa using 1D MinION sequencing. BMC Biol. 2019; 17: 96. https://doi.org/ 10.1186/s12915-019-0706-9 PMID: 31783752 26. Chesson P. Mechanisms of Maintenance of Species Diversity. Annu Rev Ecol Syst. 2000; 31: 343–366. https://doi.org/10.1146/annurev.ecolsys.31.1.343 27. Ramm T, Cantalapiedra JL, Wagner P, Penner J, Ro¨del M-O, Mu¨ ller J. Divergent trends in functional and phylogenetic structure in reptile communities across Africa. Nat Commun. 2018; 9: 4697. https:// doi.org/10.1038/s41467-018-07107-y PMID: 30409973 28. Tiusanen M, Huotari T, Hebert PDN, Andersson T, Asmus A, Bêty J, et al. Flower-visitor communities of an arcto-alpine plant—Global patterns in species richness, phylogenetic diversity and ecological functioning. Mol Ecol. 2019; 28: 318–335. https://doi.org/10.1111/mec.14932 PMID: 30418699 29. Webb CO. Exploring the Phylogenetic Structure of Ecological Communities: An Example for Rain For- est Trees. Am Nat. 2000; 156: 145–155. https://doi.org/10.1086/303378 PMID: 10856198 30. Webb CO, Ackerly DD, McPeek MA, Donoghue MJ. Phylogenies and Community Ecology. Annu Rev Ecol Syst. 2002; 33: 475–505. https://doi.org/10.1146/annurev.ecolsys.33.010802.150448 31. Vasar M, Davison J, Sepp S-K, Oja J, Al-Quraishy S, Bueno CG, et al. Global taxonomic and phyloge- netic assembly of AM fungi. Mycorrhiza. 2022; 32: 135–144. https://doi.org/10.1007/s00572-022- 01072-7 PMID: 35138435

32. Bohmann K, Lynggaard C. Transforming terrestrial biodiversity surveys using airborne eDNA. Trends Ecol Evol. 2023; 38: 119–121. https://doi.org/10.1016/j.tree.2022.11.006 PMID: 36513529 33. Early R, Keith SA. Geographically variable biotic interactions and implications for species ranges. Dor- nelas M, editor. Glob Ecol Biogeogr. 2019; 28: 42–53. https://doi.org/10.1111/geb.12861 34. Mod HK, Heikkinen RK, Le Roux PC, Wisz MS, Luoto M. Impact of biotic interactions on biodiversity varies across a landscape. J Biogeogr. 2016; 43: 2412–2423. https://doi.org/10.1111/jbi.12794 35. Antonelli A, Smith RJ, Perrigo AL, Crottini A, Hackel J, Testo W, et al. Madagascar’s extraordinary biodi- versity: Evolution, distribution, and use. Science. 2022; 378: eabf0869. https://doi.org/10.1126/science. abf0869 PMID: 36454829 36. Kremen C, Cameron A, Moilanen A, Phillips SJ, Thomas CD, Beentje H, et al. Aligning Conservation Priorities Across Taxa in Madagascar with High-Resolution Planning Tools. Science. 2008; 320: 222– 226. https://doi.org/10.1126/science.1155193 PMID: 18403708 37. Wilme´ L, Goodman SM, Ganzhorn JU. Biogeographic Evolution of Madagascar’s Microendemic Biota. Science. 2006; 312: 1063–1065. https://doi.org/10.1126/science.1122806 PMID: 16709785 38. Goodman SM, editor. The New Natural History of Madagascar. Princeton University Press; 2022. 39. Vences M, Wollenberg KC, Vieites DR, Lees DC. Madagascar as a model region of species diversifica- tion. Trends Ecol Evol. 2009; 24: 456–465. https://doi.org/10.1016/j.tree.2009.03.011 PMID: 19500874

40. Suzzi-Simmons A. Status of deforestation of Madagascar. Glob Ecol Conserv. 2023; 42: e02389. https://doi.org/10.1016/j.gecco.2023.e02389 PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 14 / 15

41. Myers N, Mittermeier RA, Mittermeier CG, Da Fonseca GAB, Kent J. Biodiversity hotspots for conserva- tion priorities. Nature. 2000; 403: 853–858. https://doi.org/10.1038/35002501 PMID: 10706275 42. Hanski I, Koivulehto H, Cameron A, Rahagalala P. Deforestation and apparent extinctions of endemic forest beetles in Madagascar. Biol Lett. 2007; 3: 344–347. https://doi.org/10.1098/rsbl.2007.0043 PMID: 17341451 43. Venter O, Sanderson EW, Magrach A, Allan JR, Beher J, Jones KR, et al. Sixteen years of change in the global terrestrial human footprint and implications for biodiversity conservation. Nat Commun. 2016; 7: 12558. https://doi.org/10.1038/ncomms12558 PMID: 27552116 44. Slagkryssaren. LIFEPLAN mobile application. 2020. https://play.google.com/store/apps/details?id= com.lifeplanapp. 45. Slagkryssaren. LIFEPLAN web admin. 2020. https://lifeplan.luomus.fi/. PLOS ONE LIFEPLAN: A worldwide biodiversity sampling design PLOS ONE | https://doi.org/10.1371/journal.pone.0313353 December 31, 2024 15 / 15