Environmental DNA Sampling Provides New Management Strategies for Vernal
Total Page:16
File Type:pdf, Size:1020Kb
bioRxiv preprint doi: https://doi.org/10.1101/2020.11.20.391052; this version posted November 20, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available for use under a CC0 license. 1 2 3 4 Environmental DNA sampling provides new management strategies for vernal 5 pool branchiopods in California 6 7 Shannon Rose Kieran1*, Joshua Hull2, Amanda Finger1 8 9 10 11 1Department of Animal Science, University of California, Davis, Davis, California, USA 12 2Sacramento Fish and Wildlife Office, United States Fish and Wildlife Service, Sacramento, 13 California, USA 14 15 16 * Corresponding author 17 E-mail: [email protected] 18 19 20 21 22 23 bioRxiv preprint doi: https://doi.org/10.1101/2020.11.20.391052; this version posted November 20, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available for use under a CC0 license. 24 ABSTRACT 25 California’s vernal pools are declining ecosystems that support valuable native plant and animal 26 diversity. Vernal pool branchiopods are particularly at risk from vernal pool habitat loss and 27 conservation efforts have targeted their long-term protection through the establishment of 28 preserves and conservation banks. These conservation strategies require repeated, perpetual 29 monitoring of preserved habitat, which is currently carried out through dip-net surveys and visual 30 identification of specimens. Dip-netting may be destructive and frequently requires some 31 sacrifice of protected species. Environmental DNA offers a new, modern method to monitor 32 many protected freshwater organisms. We designed qPCR-based species-specific assays for four 33 of California’s vernal pool branchiopods: The Vernal Pool Fairy Shrimp Branchinecta 34 lynchi (BRLY), the Midvalley Fairy Shrimp Branchinecta mesovallensis (BRME), and the 35 Conservancy Fairy Shrimp Branchinecta conservatio (BRCO), and the Vernal Pool Tadpole 36 Shrimp Lepidurus packardi (LEPA). We tested these assays using eDNA sampling protocols 37 alongside traditional dip-net surveys to assess their viability as an alternative method to monitor 38 vernal pool branchiopods. Based on occupancy modeling, each of our assays achieved a 95% or 39 higher detection rate when using optimized sampling protocols. bioRxiv preprint doi: https://doi.org/10.1101/2020.11.20.391052; this version posted November 20, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available for use under a CC0 license. 40 INTRODUCTION 41 42 Vernal pool wetlands support high levels of biodiversity and provide a wide range of 43 ecosystem services [1], but are facing decline. California’s vernal pool habitats support 44 ecologically and phylogenetically distinct biota, including a diverse assemblage of endemic 45 branchiopod crustaceans [2]. Since the mid-1800s, alteration of vernal pool wetlands and 46 conversion to agricultural and urban landscapes are estimated to have contributed to the 47 extinction of 15-33% of crustacean species in Central Valley vernal pools [3]. Of the remaining 48 vernal pool crustacean species, six are listed as threatened or endangered [4]. 49 50 A variety of conservation approaches have been implemented to help slow wetland loss 51 and conserve remaining vernal pool biodiversity, with particular emphasis on conservation of 52 vernal pool crustaceans . Conservation banks and habitat conservation plans are two important 53 tools that establish permanent, managed habitat protections for vernal pools. These land 54 ownership plans are designed to offset small-scale habitat loss from development and facilitate 55 the protection of strategically-located, large, unfragmented habitats. A central part of these 56 agreements are provisions for repeated, long-term monitoring of all listed species to help ensure 57 their persistence within protected lands. Repeated monitoring ensures that the habitat bank 58 perpetually supports the target, but comes at the cost of human impact on lands that are often 59 pristine and delicate. Monitoring tools that minimize repeated impact would be useful for long- 60 term monitoring of vernal pool conservation banks. 61 bioRxiv preprint doi: https://doi.org/10.1101/2020.11.20.391052; this version posted November 20, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available for use under a CC0 license. 62 Currently, vernal pool branchiopods (Arthropoda: Crustacea: Branchiopoda) are 63 primarily monitored during the wet season using dip-net surveys. Surveyors wade into vernal 64 pools and move a fine mesh net through the water column to disrupt the substrate and capture 65 benthic invertebrates. Once captured, adult vernal pool branchiopods are identified visually in- 66 field or in the lab under a microscope. Drawbacks to this labor-intensive survey method include 67 difficulty identifying morphologically-similar branchiopod species in-field, habitat disturbance, 68 and death of specimens used for laboratory identification. Additionally, no detection rates have 69 been established for dip-net surveys of adult vernal pool branchiopods, and juveniles in common 70 target genera cannot be identified to species visually. Long-term monitoring plans would benefit 71 from a survey method that was repeatable, with a known detection rate, that minimized the 72 sacrifice of protected species while providing highly accurate results. Environmental DNA 73 (eDNA) provides a viable alternative for sampling many habitats, including vernal pools, and 74 monitoring species of conservation concern [16-22]. 75 76 Here, we developed eDNA monitoring methods for three morphologically-similar vernal 77 pool species of conservation concern: The Vernal Pool Fairy Shrimp Branchinecta 78 lynchi (BRLY), the Midvalley Fairy Shrimp Branchinecta mesovallensis (BRME), and the 79 Conservancy Fairy Shrimp Branchinecta conservatio (BRCO), as well as for the endangered 80 Vernal Pool Tadpole Shrimp Lepidurus packardi (LEPA). We designed four species-specific 81 qPCR assays to detect the presence of each of these four large branchiopods, which are 82 commonly targeted and monitored in Central Valley dip-net surveys. We tested these protocols 83 on tissue-derived DNA and using eDNA samples collected from existing conservation lands 84 alongside dip-net surveys to determine their viability as an alternative monitoring method. bioRxiv preprint doi: https://doi.org/10.1101/2020.11.20.391052; this version posted November 20, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available for use under a CC0 license. 85 86 87 RESULTS 88 Laboratory Development and Validation 89 To develop qPCR assays able to identify our four species in the laboratory, we used 90 Sanger sequences from whole specimens of each target species sourced from at least two sites in 91 California. For each species we sequenced multiple genes to determine the best primer and probe 92 sequences (“assay”) for each species (see Materials and Methods). Using individuals from 93 multiple geographic sites allowed us to account for inter-population, intraspecific single- 94 nucleotide polymorphisms (SNPs) which might impact the assay efficiency. We determined the 95 best gene and assay independently for each species. We developed assays on the 12S Ribosomal 96 RNA Gene for BRLY and BRME, the 16S Ribosomal RNA Gene for LEPA and Cytochrome 97 Oxidase I for BRCO (see Table 1 for the sequences of each assay). We tested the assays on 98 tissue-derived target DNA and optimized our thermocycling protocols. Each assay amplifies 99 tissue-derived target DNA within 20 cycles of qPCR, except for BRLY, which amplifies within 100 30 cycles of PCR. Our optimized reaction recipe used 1X Taqman Environmental Mastermix 101 (ThermoFisher), 0.9 uM each forward and reverse primer, 0.15 uM probe, and 1X bovine serum 102 albumin. The thermocycling protocol consisted of a 10 minute initial denaturing at 95 °C 103 followed by 50 cycles of denaturing at 95 °C and a one minute elongation step that varied in 104 temperature for each assay, see Table 2 for full conditions. 105 Table 1. The complete DNA sequences of the primers and probes used in each species- 106 specific assay. Species Gene Oligo Assay Sequences bioRxiv preprint doi: https://doi.org/10.1101/2020.11.20.391052; this version posted November 20, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder. This article is a US Government work. It is not subject to copyright under 17 USC 105 and is also made available for use under a CC0 license. Primer F TGCAGAAAGGGGAGGATARACC Branchinecta conservatio (BRCO) Cytochrome Oxidase I Primer R TGCCTCCTGCCTTRACCTTRC Probe1,2 6-FAM/YCACCCAGT/ZEN/CCCAGCTCCACT/IBFQ/ Primer F GGATTTGGCGGTTCTTAAACTT Branchinecta lynchi (BRLY) 12S Ribosomal RNA Primer R TTTTCCTAGAAAAATGCATCCGT Probe1,2 6-FAM/TYAACAGCT/ZEN/TATATACCGTCGTTTAGAGGA/IBFQ/ Primer F CCGTCGCTTAGAGGATTACATT Branchinecta mesovallensis (BRME) 12S Ribosomal RNA Primer R ATGAGCTACGCCTTGATCTG Probe1,3 6-FAM/TTTAAATTCTTTTATTGGGAGTTCC/BHQ-1 Primer F CCGTGCGAAGGTAGCATAAT Lepidurus packardi (LEPA) 16S Ribosomal RNA Primer R AGGGTCTTATCGTCCCTCAA Probe1,3 6-FAM/TGAAGGCTGGTATGAATGGCTGGA /BHQ-1 107 1The probe is synthesized with a 5’ FAM fluorophore. 2The probe uses the proprietary IDT 108 internal ZEN™ quencher and 3’ Iowa Black® Quencher FQ. 3The probe uses the proprietary 109 IDT Black Hole Quencher®-1.