Biodiversity Data Journal :
Research Article
|
Corresponding author: Leighton J Thomas (l.thomas@leibniz-lib.de), Sarah J Bourlat (S.Bourlat@leibniz-zfmk.de)
Academic editor: Paolo Biella
Received: 16 May 2024 | Accepted: 04 Jul 2024 | Published: 30 Jul 2024
© 2024 Leighton J Thomas, Ameli Kirse, Hanna Raus, Kathrin Langen, Björn Nümann, Georg Tschan, Birgit Gemeinholzer, J. Wolfgang Wägele, Sarah Bourlat
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Thomas LJ, Kirse A, Raus H, Langen K, Nümann B, Tschan GF, Gemeinholzer B, Wägele JW, Bourlat S (2024) Synchronised monitoring of plant and insect diversity: a case study using automated Malaise traps and DNA-based methods. Biodiversity Data Journal 12: e127669. https://doi.org/10.3897/BDJ.12.e127669
|
The occurrence and distribution of insects and their possible associations with plant species are largely unknown in Germany and baseline data to monitor future trends are urgently needed. Using newly-designed automated Malaise trap multi-samplers, the occurrence of insect species and their potential associations with plants was monitored synchronously at two contrasting field sites in Germany: an urban botanical garden and a forest research station. Taxa were identified by metabarcoding of the insects and the plant traces present in the preservative ethanol of the Malaise trap samples. For comparison, a botanical survey was conducted in the vicinity of the traps. Across both sites, we identified a total of 1290 exact sequence variants (ESVs) assigned to Insecta, of which 205 are known to be pollinators. In the botanical garden, we detected the occurrence of 128 plant taxa, of which 41 also had one of their known insect pollinator species detected. Insect species richness was highest in May, mainly attributed to an increase in Diptera. These results present a case study of the applicability of automated sampling and DNA-based methods to monitor the timings of flowering and corresponding activity of plant-visiting insects.
metabarcoding, plant traces, plant-insect interaction, automated Sampling, Malaise trap, biomonitoring, phenology
A comprehensive understanding of insect diversity and distributions requires consistent and standardised monitoring. However, monitoring of insects is currently hindered by the complexity and time-consuming nature of species identification, which is further impeded by a decline in trained taxonomists (
Despite their well-established importance, concerns still persist regarding the decline of pollinator and plant populations (
Here, we use automated Malaise trap samplers to monitor insect and plant diversity over time, as described in
In this study, we focus on the detection of insects and their associated plants, greatly increasing the knowledge of local species occurrences. Since the data have a temporal component, information on insect flight times can also be obtained. The co-occurrences of plant traces and insects can be traced across a temporal scale, making it possible to establish a baseline dataset of flower-visiting insect flight times for future comparisons. Using these newly-developed tools, we sought to answer the following questions: which plant species can be detected from pollen and other plant traces accumulated in the insect samples? Do the detected plants match existing species lists from the vicinity of the traps? Do temporal patterns of insect occurrence correspond to the detection of individual plant species as recorded in the database of pollinator interactions (DoPI,
The insect samplers are Malaise traps complemented with an automated multisampler, containing 12 sample bottles and a control unit. The sample vials can be changed automatically following predetermined time settings (see
The samples collected using the insect multisampler. Due to a malfunction of the Malaise trap at the Melbgarten collection site, only samples B01 to B08 were collected. Of these samples, B01 to B05 were used to identify the insect species and samples B01 to B08 were used to identify the plant species (see Results section for details).
Sampling Interval |
Collection time per bottle |
Melbgarten (Bonn) |
Britz (Eberswalde) |
1 |
31 Mar – 14 Apr 2022 |
B01 |
B01 |
2 |
14 Apr – 28 Apr 2022 |
B02 |
B02 |
3 |
28 Apr – 12 May 2022 |
B03 |
B03 |
4 |
12 May – 26 May 2022 |
B04 |
B04 |
5 |
26 May – 09 June 2022 |
B05 |
B05 |
6 |
09 June – 23 June 2022 |
B06 |
B06 |
7 |
23 June – 07 July 2022 |
B07 |
B07 |
8 |
07 July – 21 July 2022 |
B08 |
B08 |
9 |
21 July – 04 Aug 2022 |
- |
B09 |
10 |
04 Aug – 18 Aug 2022 |
- |
B10 |
11 |
18 Aug – 01 Sep 2022 |
- |
B11 |
12 |
01 Sep – 14 Sep 2022 |
- |
B12 |
Total number of samples |
5 insect / 8 plant |
12 insect / 12 plant |
The Melbgarten contains both planted and naturally occurring species, which are well documented. The planted material originates predominantly from Georgia and China, but large parts of the garden are not actively cultivated. A comprehensive list of all of the planted species in the Melbgarten was provided by the curators (Supplementary Table 2). Naturally occurring species within a five-metre radius around the traps were identified using standard botanical identification literature, supported by the identification apps ‘ObsIdentify’ (
Following removal of insects by sieving, the preservative ethanol from the Malaise trap samples was vacuum filtered using a cellulose nitrate membrane (GVS Filter Technology, Sanford, USA; diameter 47 mm and 0.22 μl pore size). After filtration, the cellulose nitrate filter was cut in half and each half was placed in a separate 2 ml SafeSeal micro tube (Sarstedt AG & Co. KG, Nümbrecht, Germany). One half of the filter paper was used for DNA isolation, the other half was kept as a voucher.
DNA extraction from the plant parts consumed or pollen externally attached to insects found in the preservative ethanol of the Malaise trap was performed with the NucleoMag Plant Kit (Macherey-Nagel, Düren, Germany). Two DNA extraction negative controls were processed for each batch of up to 24 samples. The detailed protocol for DNA extraction from plant material, as well as clean lab methods used for quality assurance, are described in
The ITS2 region was used to generate amplicon libraries using a dual PCR protocol in which the barcode region was amplified in the first round of PCR and the flow cell binding adapter, a sequencing primer binding site and an index were attached in the second round of PCR (
Demultiplexing and primer removal was performed using Cutadapt v.1.9.1 (
Preservative ethanol was removed from the Malaise trap samples and the ethanol stored for the subsequent extraction of plant material. Insects were sieved using a 4 x 4 mm mesh (wire diameter 0.5 mm, untreated stainless steel), resulting in two size fractions: S (small, ≤ 4 mm) and L (large, > 4 mm). Individuals of both size fractions were transferred to either disposable grinding chambers (IKA, 40 or 100 ml) or 30 ml Nalgene tubes with metal beads (5 mm in diameter) and dried in an incubator at 50°C for up to 5 days until complete ethanol evaporation. Dried insect tissue was homogenised for 3 minutes either with a batch mill (Tube Mill 100 Control, IKA) at 25,000 rpm or a mixer mill (MM400, Retsch) at 30 Hz for 5 minutes. Approximately 20 mg of finely homogenised tissue were transferred to a 1.5 ml Eppendorf tube and 190 µl of ATL buffer (Qiagen, Hilden, Germany) and 10 µl of Proteinase K (Qiagen, Hilden, Germany) were added. The samples were incubated overnight at 56°C using a shaking incubator to allow for tissue lysis (110 rpm, INCU Line ILS 6, VWR, Radnor, PA, USA). Twelve negative controls containing 200 μl of ATL buffer only were added for each batch of 84 samples during processing in 96 well plate format. DNA was extracted from the Melbgarten samples for each size fraction (S and L) separately with the DNeasy 96 Blood and Tissue Kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. For the Britz samples, 135 µl of the S fraction lysate and 15 µl of the L fraction lysate of each sample were merged according to
The dual PCR protocol of
PCR success was checked on a 1% agarose gel before PCR products were normalised using a SequalPrep normalisation plate (Thermo Fisher Scientific, MA, USA) following the manufacturer’s instructions, resulting in a final DNA yield of 25 ng per sample (20 μl volume). For each sample, 10 μl aliquots were pooled and two rounds of left-sided size selection were carried out on the sample pool with magnetic beads to remove primer dimers (ratio 1:0.7, SPRIselect Beckman Coulter). Library concentration was measured with a Quantus fluorometer (Promega, Madison, USA) and on a Fragment Analyzer (Agilent Technologies, Santa Clara, CA, USA). The pool was sent for sequencing on two Hiseq 2500 runs (2 x 250 bp) (Macrogen Europe B.V., Netherlands). Raw data were uploaded to the GenBank SRA archive under bioproject accession number PRJNA1068928.
Following the APSCALE pipeline (
A) Total pollinator species richness (orange) across time for both Britz (sampling points 1-12) and the Melbgarten (sampling points 1-5). Species richness across time per order is presented for Coleoptera (red), Diptera (blue), Hymenoptera (green) and Lepidotera (purple); B) Community composition bar plot showing the number of pollinator species per order at both sampling sites; C) Taxon richness across time for plants detected in the insect multisampler.
UpsetR plot showing the species overlap between sampling points for the insect pollinators. Horizontal bars on the left indicate the total number of detected insect pollinator species per sampling point at Britz and the Melbgarten (Melb). Vertical bars indicate the number of shared and unique species within and between sampling points, as well as their taxonomic composition.
UpsetR plot showing the overlap between sampling points for the plants found in the Malaise trap samples. Horizontal bars on the left indicate the total number of detected plant taxa per sampling point at Britz and the Melbgarten (Melb). Vertical bars indicate the number of shared and unique taxa within and between sampling points, as well as their taxonomic composition.
Due to a malfunction of the insect multisampler at the Melbgarten sampling site, only eight pollen and five insect samples were reliably collected. From sampling interval 9 onwards (see Table
A total of 163 ESVs detected in the Britz site were assigned to 41 families, of which 45 could only be determined to genus level and one could only be assigned to family level, as polyploidy and hybridisation in Hieracium and Pilosella hindered unambiguous genus assignment. A total of 128 ESVs representing 37 different plant families were recorded at Melbgarten, of which 31 could only be determined to genus level and one could only be assigned to family level (Suppl. material
Across the two sampling sites, a total of 1290 ESVs assigned to Insecta were detected (Suppl. material
Across both sites, a total of 205 insect species recognised as pollinators were identified (Suppl. material
Taxonomic information for the insect species shown in Fig.
Image | Order | Family | Species | Author | Sex | Month | Potential plants pollinated according to DoPI ( |
A | Hymenoptera | Vespidae | Vespula vulgaris | (Linnaeus, 1758) | female | July | Borago officinalis, Calluna vulgaris, Centaurea cyanus, Crepis capillaris, Daucus carota, Hedera helix, Heracleum sphondylium, Pastinaca sativa, Plantago lanceolata, Potentilla reptans, Trifolium repens, Tripleurospermum inodorum |
B | Hymenoptera | Apidae | Bombus pascuorum | (Scopoli, 1793) | female | August | Ajuga reptans, Alliaria petiolata, Allium ursinum, Bellis perennis, Borago officinalis, Clematis vitalba, Crepis capillaris, Daucus carota, Diplotaxis tenuifolia, Dipsacus fullonum, Glechoma hederacea, Heracleum sphondylium, Hypericum perforatum, Hypochaeris radicata, Lamium galeobdolon, Leontodon hispidus, Lotus pedunculatus, Ononis spinosa, Papaver rhoeas, Prunella vulgaris, Ranunculus repens, Stachys sylvatica, Torilis japonica, Trifolium pratense, Trifolium repens, Vicia sepium |
C | Coleoptera | Cerambycidae | Stenurella melanura | (Linnaeus, 1758) | male | May | NA |
D | Lepidoptera | Pieridae | Pieris napi | (Linnaeus, 1758) | unknown | August | Bellis perennis, Brassica oleracea, Cardamine pratensis, Dipsacus fullonum, Eupatorium cannabinum, Geranium robertianum, Glechoma hederacea, Jasione montana, Prunella vulgaris, Ranunculus repens, Syringa vulgaris, Tripleurospermum inodorum |
E | Lepidoptera | Pieridae | Pieris rapae | (Linnaeus, 1758) | unknown | July | Bellis perennis, Brassica napus, Centaurea cyanus, Crepis capillaris, Eupatorium cannabinum, Prunus avium, Tripleurospermum inodorum |
F | Lepidoptera | Nymphalidae | Pararge aegeria | (Linnaeus, 1758) | female | July | NA |
G | Lepidoptera | Nymphalidae | Vanessa atalanta | (Linnaeus, 1758) | female | July | NA |
H | Lepidoptera | Lycaenidae | Polyommatus icarus | (Rottemburg, 1775) | male | July | Bellis perennis, Calluna vulgaris, Eupatorium cannabinum, Hypochaeris radicata |
I | Lepidoptera | Noctuidae | Noctua pronuba | (Linnaeus, 1758) | unknown | July | Hedera helix |
J | Diptera | Bibionidae | Bibio marci | (Linnaeus, 1758) | female | April | Brassica napus |
K | Diptera | Syrphidae | Myathropa florea | (Linnaeus, 1758) | male | May | Calluna vulgaris, Eupatorium cannabinum, Hedera helix |
L | Mecoptera | Panorpidae | Panorpa germanica | Linnaeus, 1758 | female | April | Heracleum sphondylium |
Examples of the insect species found in the Malaise trap sampling of this study that are also pollinators according to the Database of Pollinator Interactions (DoPI) (
Peak pollinator species richness occurred at sampling round 4 (12 May 2022 – 26 May 2022) in both locations with species richness overall higher in the Melbgarten (Fig.
In temperate regions, 78% of all flowering plants are animal-pollinated, of which most are insects (
Recording the interactions between plants and insects in the field can be tedious. In addition, the quality of the data collected depends on the training and working accuracy of each individual involved. The use of the automated insect multisamplers combined with metabarcoding techniques is a contribution to reducing sampling effort and standardising the acquisition of temporal monitoring data, in a user-independent way. In the AMMOD project, plant species were detected using: (i) airborne pollen traps (
The unveiling of plant-insect interactions using DNA-based methods is an emerging field (
Some of the co-occurrences can be revealed indirectly using statistical methods, for example, with the R package cooccur (
Some interactions between plants and insects are already well known (DoPI,
These examples emphasise the complexity of the interaction between plants and insects, some of which are known, but many of which have yet to be investigated. Until recently, it was not even clear how many of the more than 350,000 angiosperm species interacted with pollinators, with figures ranging from less than 70 and close to 100 percent; for temperate-zone communities, the proportion is – on average – a bit more than three quarters (
In our samples, there was a low degree of overlap between insects detected in each time-frame, possibly due to short species-specific flying time-frames of plant-visting insects. The low degree of overlap is a promising indication that the new methodologies can be used to determine flight and flowering times with a higher degree of precision. However, this phenomenon is also a characteristic of Malaise traps: it has been demonstrated that – even when placed in close vicinity – there is still a low degree of overlap between the species caught in adjacent traps (
The methodologies described here provide an example of how trends could be monitored in the future. For example, the first and last seasonal detection of insect and plant taxa could be used to estimate flight onset and duration for floral visitors, as well as peak flower visiting periods. Additionally, a greater understanding of the biology of insects needs to be considered and, here, the data can be used to elucidate different phenological patterns. For example, Andrena cineraria (Hymenoptera) has multiple generations per year, but in our study, the species was detected in only one time interval at both sites; it might perhaps be detected for each of its generation in a given year if monitoring occurred over longer temporal scales. The example also shows that the methods must be standardised and the instruments calibrated. Nevertheless, long term monitoring could, indeed, be used to investigate phenological changes triggered by climate change.
Consistent and standardised methods for monitoring populations are needed for the assessment of extinction risk faced by invertebrate species (
The Automated Multisensor stations for Monitoring of bioDiversity (AMMOD) project was funded by the German Federal Ministry of Education and Research (Bundesministerium für Bildung und Forschung), project number 16LC1903. We thank Marco Nathkin for technical support at the Britz site and the staff of the University of Bonn Botanic Gardens for their help during many visits to the Melbgarten.
Presence/absence of insect ESVs detected in the insect malaise trap multisampler in samples from the Melbgarten and Britz sites. Percentage similarity to BOLD records, Red List category and invasive or pollinator status of the species are also shown.
Presence/absence of plant ESVs detected in the insect Malaise trap multisampler in samples from the Melbgarten and Britz sites.
This table shows the plant-pollinator interactions found in the Melbgarten. In order to review the presence of known pollinator interactions, the insect species detected in our study were compared to the database of pollinator interactions (DoPI) (https://www.sussex.ac.uk/lifesci/ebe/dopi/) (Balfour et al. 2022).
This table shows the plant-pollinator interactions found at the Britz site. In order to review the presence of known pollinator interactions, the insect species detected in our study were compared to the database of pollinator interactions (DoPI) (https://www.sussex.ac.uk/lifesci/ebe/dopi/) (Balfour et al. 2022).