Non-invasive genetic monitoring for the threatened valley elderberry longhorn beetle

Abstract

The valley elderberry longhorn beetle (VELB), Desmocerus californicus dimorphus (Coleoptera: Cerambycidae), is a federally threatened subspecies endemic to the Central Valley of California. The VELB range partially overlaps with that of its morphologically similar sister taxon, the California elderberry longhorn beetle (CELB), Desmocerus californicus californicus (Coleoptera: Cerambycidae). Current surveying methods are limited to visual identification of larval exit holes in the VELB/CELB host plant, elderberry (Sambucus spp.), into which larvae bore and excavate feeding galleries. Unbiased genetic approaches could provide a much-needed complementary approach that has more precision than relying on visual inspection of exit holes. In this study we developed a DNA sequencing-based method for indirect detection of VELB/CELB from frass (insect fecal matter), which can be easily and non-invasively collected from exit holes. Frass samples were collected from 37 locations and the 12S and 16S mitochondrial genes were partially sequenced using nested PCR amplification. Three frass-derived sequences showed 100% sequence identity to VELB/CELB barcode references from museum specimens sequenced for this study. Database queries of frass-derived sequences also revealed high similarity to common occupants of old VELB feeding galleries, including earwigs, flies, and other beetles. Overall, this non-invasive approach is a first step towards a genetic assay that could augment existing VELB monitoring and accurately discriminate between VELB, CELB, and other insects. Furthermore, a phylogenetic analysis of 12S and 16S data from museum specimens revealed evidence for the existence of a previously unrecognized, genetically distinct CELB subpopulation in southern California.

About the Authors:

Raman P. Nagarajan

Roles Investigation, Writing – original draft

* E-mail: [email protected]

Affiliation: Department of Animal Science, University of California Davis, Davis, CA, United States of America

ORCID logo http://orcid.org/0000-0003-4422-9934

Alisha Goodbla

Roles Investigation, Writing – review & editing

Affiliation: Department of Animal Science, University of California Davis, Davis, CA, United States of America

Emily Graves

Roles Investigation, Writing – review & editing

Affiliation: Department of Environmental Science and Policy, University of California Davis, Davis, CA, United States of America

Melinda Baerwald

Roles Conceptualization, Funding acquisition, Investigation, Writing – review & editing

Current address: California Department of Water Resources, Sacramento, CA, United States of America

Affiliation: Department of Animal Science, University of California Davis, Davis, CA, United States of America

Marcel Holyoak

Roles Conceptualization, Supervision, Writing – review & editing

Affiliation: Department of Environmental Science and Policy, University of California Davis, Davis, CA, United States of America

Andrea Schreier

Roles Conceptualization, Funding acquisition, Project administration, Supervision, Writing – review & editing

Affiliation: Department of Animal Science, University of California Davis, Davis, CA, United States of America

ORCID logo http://orcid.org/0000-0002-2241-3119

Introduction

The valley elderberry longhorn beetle (VELB) is a wood-boring subspecies of beetle of the Cerambycidae found only in the Central Valley of California, U.S.A. [1] (Fig 1). The subspecies is sexually dimorphic, with males displaying red-orange elytra (wing covers) with four dark elongate spots, and females showing dark, metallic green to black elytra with bright red-orange borders [2]. VELB inhabit elderberry shrubs (Sambucus spp.) for nearly their entire life cycle, and are commonly found in riparian forests and adjacent uplands near Central Valley waterways [3,4]. After hatching from eggs laid externally on leaves and stems, VELB larvae bore into, feed on the pithy center, and pupate within elderberry stems, culminating with the emergence of adults through distinctive exit holes [5]. During their 1–2 year inhabitation of the stems, VELB excavate feeding galleries and deposit fecal material mixed with wood shavings (frass) that remain after adults have exited the hole [5–7]. After emergence, adult males can live up to 5 days under laboratory conditions, whereas adult females live up to 3 weeks [8].

[Figure omitted. See PDF.]

Fig 1. Map of VELB and CELB ranges within California.

Individual points show frass collection locations that yielded sequencing data. One frass sample that was collected in the southern Central Valley outside of the VELB and CELB official ranges is labeled “VELB/CELB” to indicate uncertainty about which subspecies might be present. Arrows show locations where frass-derived sequences matched VELB or CELB museum specimens collected from nearby locations.

https://doi.org/10.1371/journal.pone.0227333.g001

Due to declining populations resulting from habitat loss, in 1980 VELB was listed as threatened under the U.S. Endangered Species Act [9]. The Central Valley of California has been extensively modified by farming, urbanization, and flood control measures, diminishing and altering the riparian habitats critical for VELB [10,11]. A proposal to delist the subspecies was made in 2012 [12] but ultimately withdrawn by the U.S. Fish and Wildlife Service after external review [11,13]. VELB continues to be threatened by habitat fragmentation and loss, climate change, predators, pesticides and invasive species [11,14].

Despite its conservation significance, VELB has been historically difficult to survey and study. VELB population densities are low, and field surveys suffer from low probability of detection [4,15]. Live adult VELB are infrequently encountered in the wild [2,4,7], and occurrence records are mostly based on observation of exit holes in elderberry stems, which themselves occur at low density in elderberry stands. Exit hole surveys are non-invasive but can be inaccurate due to misidentification of holes produced by other xylophagous beetles or wasps [16]. The degree of misidentification of holes is not known. When adult VELB are (rarely) encountered, VELB females are indistinguishable from a closely related subspecies, the California elderberry longhorn beetle (CELB). CELB females are morphologically similar to VELB females, but CELB males have dark elytra bordered with a thin red margin and thus are distinguishable from VELB males [2]. In our experience in the field, the exit holes produced by VELB and CELB are indistinguishable, as might be expected with two closely related members of the same species that share many morphological features as well as life history and biology [17,18]. The CELB is not listed as threatened or endangered, although the subspecies is less studied than the VELB. CELB inhabit coastal California, parts of the Sierra Nevada, and the central and southern San Joaquin Valley and southwards, including the highly urbanized Los Angeles and San Diego areas [2,19,20] (Fig 1). VELB and CELB ranges overlap along the eastern edge of the Coast Range, and intermediate forms have been observed [2,8,16]. The exit holes of VELB cannot be distinguished from those of CELB, and in areas of overlap, exit hole survey methods are subject to ambiguity about which subspecies is present. For these reasons, there is a pressing need for an accurate survey tool that could improve the probability of detection of VELB and identify subspecies in the absence of adult specimens. A non-invasive survey technique that minimizes mortality and interference with natural VELB behaviors would be particularly useful.

Genetic barcoding of non-invasively collected DNA has been used to rapidly survey populations of many taxa including invertebrates (e.g. [21–24]). Species-specific DNA fragments can be used to document species presence in soil samples [25], gut contents [26,27], or fecal material [28]. Insect frass has been shown to be a viable source of non-invasive DNA sampling for butterfly caterpillars [29], scarab beetles [30], and bumblebees [31]. An indirect, genetic method could provide a much-needed complementary approach to existing exit hole survey methods for VELB and CELB. Frass would be an ideal source of material, since it is readily obtainable from unoccupied exit holes with minimal disturbance. In this study we generated VELB and CELB barcode reference sequences for multiple loci from a panel of museum voucher specimens. In parallel, we collected frass from exit holes at multiple locations and used nested PCR and direct Sanger sequencing to test the ability of indirect molecular techniques to identify VELB and CELB. In addition, phylogenetic analysis of VELB/CELB museum specimens identified a previously unknown, genetically distinct population of CELB in southern California.

Materials and methods

Non-destructive DNA extraction, PCR and sequencing of VELB and CELB museum specimens

Non-destructive DNA extraction involved full immersion of specimens in digestion buffer containing 3 mM CaCl2, 2% sodium dodecyl sulfate (SDS), 40 mM dithiotreitol (DTT), 250 mg/ml proteinase K, 100 mM Tris buffer pH 8 and 100 mM NaCl and incubated overnight at 55°C with gentle agitation. The specimens were then removed, immersed in 100% ethanol and air-dried. DNA was purified from the digestion buffer with the QIAquick PCR purification kit (Qiagen) following the manufacturer’s protocol. DNA concentrations were determined using a Qubit Fluorimeter along with the Broad Range dsDNA Assay Kit (ThermoFisher). Samples were normalized to a working concentration of 10 ng/μl and amplified in PCR reactions using the primers shown in S1 Table. PCR reactions were performed in a total volume of 25 μl with the following components: 1X TaKaRa Buffer, 1.4 μM MgCl2 (Roche), 0.4 μM dNTPs (TaKaRa), 0.4 μM of each primer, 0.2 μl TaKaRa Ex Taq, and 1 μl of template DNA. (For DNA concentrations less than 10 ng/μl, we used 2 μl of undiluted DNA). After initial denaturation at 94°C for 3 minutes, 35 cycles were performed of 94°C for 30s, 55°C for 60s, 72°C 60s. Final extension was 10 minutes at 72°C. PCR amplifications were visualized with 2% agarose electrophoresis, and positive PCR amplicons were purified with Ampure beads (or with gel extraction (Qiagen) if necessary) and submitted for Sanger sequencing to Quintara Biosciences. Sequencing data were analyzed with Sequencher 4.8 and alignments were made and trimmed in MEGA7 [32] using the Clustal W approach with default settings.

Field survey and sampling protocol

Thirty-six sites were surveyed for VELB or CELB in 2016–2018. We surveyed known sites within the VELB range, based on the presence of VELB in past surveys [4,5] and accessibility to the public. We deliberately selected sites to include the core area of the VELB range, suspected hybrid zones with CELB (the southern Central Valley or south of Napa Valley), and some areas known to be inhabited by CELB. At each site, all elderberry shrubs that could be visually detected from marked trails and access points were inspected for larval exit holes. Exit holes were classified as either new (beetle emergence within the current year) or old (beetle emergence prior to the current year). New holes were those that appeared light-colored inside and had not yet begun to heal if in live stems. Surveying consisted of visual scanning of main stems and branches for exit holes, and recording the numbers of exit holes found. When an exit hole was found, tweezers were used to extract any material possible from inside the elderberry stem. Frass/feces and/or exuviae were collected from each exit hole and placed into a clean plastic Eppendorf tube. When multiple exit holes were detected on the same elderberry shrub, material from each exit hole was placed in a separate tube. When an exit hole was detected in a dead stem, the stem was collected when possible and brought back to the lab and broken apart to extract any exuviae and/or frass/feces. Live stems were not collected and were only sampled by removing frass/feces non-destructively. Upon returning from the field, samples were stored in a -20°C freezer until later processing. As our surveying did not result in any incidental take of either subspecies, since it was completely limited to collecting frass, no collecting permit was required. Nearly all of the locations that we sampled were public access and did not require any permissions. The only exceptions are sites listed as "private" for which we had explicit permission to access through our personal relationships with the landowners (contingent upon the locations remaining anonymous). We obtained permission to access the Russell Ranch site (UC Davis property) from the manager via email.

DNA extraction, PCR, and sequencing of frass samples

DNA extractions for the majority of frass samples were performed with the DNeasy PowerSoil Pro Kit (Qiagen) following the manufacturer’s protocol with the minor adjustment of eluting the DNA in 60 μl of the provided elution buffer pre-warmed to 56°C. DNA concentrations were determined using a Qubit Fluorimeter along with the Broad Range dsDNA Assay Kit (ThermoFisher). Samples were normalized to a working concentration of 10 ng/μl and 2 μl of each was used in PCR reactions. PCR primers used in this study are shown in S1 Table. PCR reaction setup and cycling were the same as described above, but with 2 μl template per reaction. For nested PCR, 2 μl of a 1:100 dilution of the first PCR product was used as the template for the second PCR. Amplification products were checked by 2% agarose gel electrophoresis, and PCR amplicons were purified with Ampure beads (or with Qiagen gel extraction if multiple bands were observed) and submitted for sequencing at Quintara Biosciences. Sequencing data were analyzed with Sequencher 4.8 and alignments were visualized with MEGA 7. To determine the species from which a frass specimen originated, frass sequences were queried with Basic Local Alignment Search Tool–Nucleotide (BLASTN) and compared to 12S and 16S rRNA alignments from VELB/CELB museum samples.

Phylogenetic tree methods

For phylogenetic tree construction, we selected 12S and 16S rRNA genes because they amplified consistently across samples. We calculated the number of transitions and transversions and nucleotide composition and used these metrics to parameterize construction of neighbor-joining trees in MEGA7 [33]. The 12S and 16S genes were concatenated before tree construction for all individuals that contained complete sequence at those genes (18 VELB, 18 CELB). Gaps and missing data were deleted and uniform evolutionary rates among sites and lineages was assumed. The only gap present was at 126–128 bp in 12S. (Analysis was performed with and without gap removal and the same topology was observed.) Trees were constructed using the Tamura 3-parameter model [34] and tested using 1000 bootstrap replicates [35].

Results

To generate reference barcode sequences for frass genetic identification, we first obtained VELB and CELB dry specimens from multiple museum collections. The specimens were collected across a wide geographic area, encompassing 22 California counties (S2 Table). To preserve these valuable specimens, we utilized a non-destructive DNA extraction protocol based on Thomsen and coworkers [36] to extract DNA from 44 specimens (19 VELB and 25 CELB), including those collected as far back as 1914 (S1 Fig). We tested 27 primer pairs targeting barcoding loci, several of which were developed for this study (S1 and S3 Tables). Fourteen primer pairs yielded successful PCR amplification for at least a subset of museum samples, and eight of these generated clean sequencing data (S3 Table). These data were from the following genes: cytochrome c oxidase 1 (COI, mitochondrial, 1 primer pair); 18S rRNA (nuclear, 2 primer pairs); 16S rRNA (mitochondrial, 2 primer pairs); and 12S rRNA (mitochondrial, 3 primer pairs). All sequences were deposited to GenBank (S4 Table). Sequences from all samples were aligned and inspected for the presence of variable sites. For COI and 18S, none were found. In contrast, multiple variable sites were identified for amplicons in the 12S and 16S genes (S5 Table). We searched these sequences for potential subspecies-specific SNPs that would enable discrimination between VELB and CELB, but none were found at these loci (S3 and S5 Tables).

To test the possibility of genetic detection from frass, samples were collected from elderberry exit holes in conjunction with VELB field survey work [37] (Fig 2A and 2B). Across these survey sites, 535 individual shrubs were examined for VELB/CELB exit holes (Table 1). Of the 212 total exit holes detected, 118 were old holes (beetle emergence occurred in a previous year) while 34 holes were determined to be new (beetle emergence occurred earlier in the current season). Additional CELB frass samples were collected separately in southern California by a consulting company. Altogether, 84 frass samples were collected.

[Figure omitted. See PDF.]

Fig 2. VELB exit holes in elderberry shrubs.

A. Archive photo from 1997 showing VELB pupae and frass in a live elderberry stem. B. New (current year) exit hole in elderberry stem from field survey/sampling performed during this study.

https://doi.org/10.1371/journal.pone.0227333.g002

[Figure omitted. See PDF.]

Table 1. Exit hole survey and frass collection.

https://doi.org/10.1371/journal.pone.0227333.t001

DNA extractions from frass samples yielded DNA concentrations ranging from 1–580 ng/μl, with an average concentration of 63 ng/μl. DNA was extracted from 84 samples, and a subset of these was checked for DNA quality by agarose gel electrophoresis. All samples examined showed mostly intact DNA, with relatively minor degradation, by agarose gel electrophoresis (S2 Fig).

Single band PCR products were successfully amplified from 23 frass samples and 26 sequences were queried using BLASTN (Table 2). The top BLASTN hits (all > 80% identity) for all 26 sequences corresponded to insects, including other beetles, earwigs, flies, and other insects. Since neither VELB nor CELB 12SS or 16S sequences were present in the GenBank database during the query, none of the top BLASTN hits corresponded to VELB or CELB. However, three 12S top hits were from other longhorn beetles (from the same taxonomic family, Cerambycidae, as VELB and CELB). We found that for these three samples, the 12S sequences exactly matched sequences from VELB or CELB museum specimens (Table 2 and S3 Fig). Furthermore, the frass sequences were identical to museum specimens collected at nearby locations. One frass specimen collected from the southern portion of the VELB range showed 100% sequence identity with VELB/CELB museum specimens from the VELB and northern CELB ranges (Fig 1, top arrow and S3 Fig). Two frass specimens collected approximately 300 meters from one another in San Bernardino county in southern California, in the southern part of the CELB range, showed 100% sequence identity with CELB museum specimens collected from southern California (Fig 1, bottom arrow and S3 Fig).

[Figure omitted. See PDF.]

Table 2. Top BLASTN hits for 12S and 16S rRNA frass-derived sequences.

https://doi.org/10.1371/journal.pone.0227333.t002

To further understand the genetic relationships across VELB and CELB populations, we used the concatenated 12S and 16S sequences from the museum specimens to construct a neighbor-joining phylogenetic tree (Fig 3). We identified 21 variable sites across 277 base pairs of the 12S rRNA sequence data, after removal of a three base pair gap in 18 VELB and 23 CELB samples. Nucleotide composition was A/T biased, with bases 41.9–43.0% T, 38.6–40.1% A, 11.2–12.6% C, and 6.1–6.9% G. Fifteen transitions and six transversions were identified. In the 16S rRNA sequence data, 13 variable sites were identified in 20 VELB and 21 CELB samples across 237 bases, with no gaps or missing data. Nucleotide composition at this locus was A/T biased, with bases 43.5–45.1% T, 39.2–40.9% A, 9.7–11.4% G, and 5.1–5.5% C. Eight transversions and five transitions were identified. Two sites were highly variable, with three (tri-allelic) and four (quad-allelic) different nucleotides detected.

[Figure omitted. See PDF.]

Fig 3. Neighbor-Joining tree constructed from 514 base pairs of the 12S and 16S rRNA genes.

The optimal tree with the sum of branch lengths = 0.071 is shown. Three monophyletic clades are labeled “A,” “B,” and “C.” Bootstrap values for highly supported nodes (>75%) are shown above the branches. Nodes without values indicate poor (<75%) bootstrap support. County from which each sample was collected is in parentheses, triangles indicate samples from areas of potential range overlap, and a star indicates a museum specimen initially misidentified as CELB.

https://doi.org/10.1371/journal.pone.0227333.g003

Three highly supported clades were identified in the neighbor-joining tree constructed using the Tamura 3-parameter model [34] (Fig 3). Clade A contained most VELB specimens and 33% of CELB specimens, all originating from northern California collection locations. All sub-branches within Clade A had relatively low bootstrap support, indicating statistically uncertain topology within this branch. Clade B consisted of the two remaining VELB that were collected from San Joaquin county, a location of potential range overlap for the two subspecies. Clade C included only CELB specimens originating from the southernmost part of that subspecies’ range.

Discussion

Insect frass has emerged as a non-invasive DNA source for multiple types of genetic studies, including those of beetles. This has included studies of pathogen infection of both beetles and their host plants [38] and dietary analysis [39]. Frass has also been utilized for molecular species identification, which can be indispensable in distinguishing closely related species [30]. In addition, frass genetic methods can reveal gender and identify transgenic individuals, even in visually similar early larval stage insects [40]. A PCR-based approach has been documented for invasive wood-boring insects [41] but there is also potential for non-lethal genetic sampling of rare and protected organisms. VELB and CELB frass could provide a source of DNA from these beetles, although the quality and quantity of DNA in situ is not currently known and might be variable. However, frass collection has the advantage that it does not interfere with or harm VELB adult beetles, which are protected by law and might be sensitive to disturbance.

Several factors may affect the efficiency of detecting VELB/CELB DNA in frass. First, some of the exit holes from which frass was collected could have been created by other xylophagous beetles or insects, and were never occupied by VELB/CELB. This limitation of the exit hole survey method has been previously documented [16] and underscores that a genetic detection method could improve the accuracy of visual surveying methods. Second, even if the exit hole was created by VELB/CELB, we hypothesize that some of these frass samples contain an admixture of DNA from multiple insects, potentially at varying proportions from sample to sample. It is possible that some of the older holes are created by VELB/CELB and leave behind feces but then other insects move into the holes and deposit fresher feces or exuviae (cast-off outer skin) with more intact DNA in the frass sample. We frequently observed other insects in or around the exit holes, and earwigs were particularly common around older holes. Earwigs are insect predators and their prey include beetles. Finally, any VELB/CELB DNA will degrade over time, likely decreasing its probability of detection in older samples. Additionally, DNA from microorganisms and from the elderberry plant itself might be present at much higher concentrations than VELB/CELB DNA or DNA from other insects, further hindering detection.

Our data support effects from all of these factors. Forty-three (of 84 total) frass samples could not be amplified by PCR, suggesting low VELB/CELB DNA concentrations or DNA degradation. PCR inhibitors could also contribute to this (see below). When PCR and sequencing were successful, all of the top BLAST hits were insects, consistent with exit hole creation and/or occupancy by multiple types of insects. These included earwigs, flies, midges, and other beetles. For future studies, cloning of individual PCR amplicons or deep sequencing of the admixed PCR product would allow a more granular and comprehensive picture of the organisms present (at least for sequences that can be PCR amplified). The frass samples might also contain PCR inhibitors from elderberry wood shavings, polysaccharides from feces, or other materials [42,43]. The extraction method we used does include a PCR inhibitor removal step, but this could potentially be further optimized. Finally, some mitochondrial loci might be more challenging to amplify from frass than others. Similar to another publication [40], we were unsuccessful in amplifying COI in our beetle frass samples, although these primers successfully amplified a subset of museum specimens (S3 Table).

Despite these challenges, sequences generated from three frass samples exactly matched references obtained from VELB/CELB museum specimens, strongly suggesting indirect genetic detection. One frass-derived sequence from a northern California collection location matched multiple museum specimens from clade A in the phylogenetic analysis, and frass collected from southern California matched museum specimens in clade C (Fig 1 and S3 Fig). Furthermore, the frass data, although limited, are consistent with the phylogenetic analysis from museum specimens pointing towards a genetically distinct CELB population in southern California. Both the museum and frass samples from this putative subpopulation were all collected from sites south of the Transverse Mountain Range in Southern California, which runs east to west and might present a geographic barrier promoting allopatric differentiation. Alternatively, the observed branching pattern could be due to decreased hybridization at the southernmost end of the CELB range, which is further from the VELB range. Additional genetic data, either from non-invasively collected frass, or high-quality DNA extracted from fresh specimens, would help to conclusively determine the true taxonomic relationships among these members of the Cerambycidae. Although there is scant published genetic data about either subspecies, a doctoral thesis that performed molecular taxonomic analysis of multiple Desmocerus species found a relatively large genetic distance between VELB and CELB clades, such that the data “might appear to suggest that the two populations are easily separable as distinct species” [44]. Notably, the three VELB specimens used in that study were from Ord Bend (Northern CA) and the two CELB were from San Bernardino (southern CA), which in our phylogeny would correspond to clades A and C, respectively. The author of that study was careful to emphasize the limitations of the data, and declined to suggest any changes to the existing VELB/CELB taxonomy. Similarly, although we have interrogated more samples across a broader geographic range, our phylogenetic data should also be considered preliminary and interpreted with caution. Additional samples combined with a genome-scale approach, such as RADseq [45] would likely allow more definitive taxonomic conclusions. Since a comprehensive genetic analysis of VELB and CELB has been hindered to date due to the lack of fresh specimens for DNA extraction, one potentially productive approach involves pheromone-based trapping of live beetles, using the recently discovered attractant sex pheromone (R)-Desmolactone [46,47]. Such an approach would enable efficient collection of fresh VELB and CELB samples from multiple, specific locations across their core ranges, and can be stopped at each site when sufficient sample numbers have been collected, minimizing the number of individuals removed from potentially sensitive populations. The application of new approaches, both for sample collection as well as genetic analysis, will afford a much deeper and complete understanding of these elusive and threatened beetles.

Supporting information

[Figure omitted. See PDF.]

S1 Table. PCR primers used in this study.

Mitochondrial and nuclear DNA regions were screened in VELB and CELB museum samples to identify diagnostic loci, and nested primers for mitochondrial 12S and 16S genes were used to amplify DNA extracted from frass material.

https://doi.org/10.1371/journal.pone.0227333.s001

(XLSX)

S2 Table. Collection information for VELB and CELB beetle samples used in genetic analyses.

https://doi.org/10.1371/journal.pone.0227333.s002

(XLSX)

S3 Table. PCR and sequencing results for 27 tested primer sets on museum specimens.

https://doi.org/10.1371/journal.pone.0227333.s003

(XLSX)

S4 Table. GenBank Accession numbers for VELB and CELB museum specimens.

https://doi.org/10.1371/journal.pone.0227333.s004

(XLSX)

S5 Table. Highly variable sites (SNPs) in the mitochondrial 12S and 16S rRNA genes across VELB and CELB museum specimens.

Based on the original geographical collection location, all specimens were annotated either high or low confidence (column 3). High confidence samples are those collected from core regions of the known range for each subspecies and are less likely to be subspecies or hybrids.

https://doi.org/10.1371/journal.pone.0227333.s005

(XLSX)

S1 Fig. Two VELB and two CELB beetle specimens from the University of California Riverside's Entomology Research Museum.

Left, pre-extraction; right, post-extraction. Use of a non-destructive DNA isolation protocol allowed for successful DNA extraction without morphological damage.

https://doi.org/10.1371/journal.pone.0227333.s006

(PDF)

S2 Fig. Frass DNA extractions (samples numbered 1–24) from the Qiagen PowerSoil Pro kit.

Purified DNAs were run on a 1% agarose gel to examine DNA quality.

https://doi.org/10.1371/journal.pone.0227333.s007

(PDF)

S3 Fig. Multiple sequence alignment showing the exact DNA match between frass samples and multiple museum specimens.

Frass sample 4 showed a 100% match to multiple museums specimens that had been assigned to clade A (from the 12S/16S phylogenetic tree in Fig 2). In addition, Frass sample 4 showed a 100% match to the two museum specimens comprising clade B (not shown). Frass samples F02A and F04A both showed 100% sequence identity to specimens assigned to clade C.

https://doi.org/10.1371/journal.pone.0227333.s008

(PDF)

Acknowledgments

We thank the California Academy of Sciences, the Natural History Museum of Los Angeles County, the University of California Davis Bohart Museum of Entomology, and the University of California Riverside Entomology Research Museum for providing specimens. Robyn Powers (H.T. Harvey & Associates Ecological Consultants) collected frass samples in southern California. We thank Philip A. Schapker for sharing his thesis, Amanda Coen for assistance with map creation, Dr. Amanda Finger for advice regarding NJ tree construction, and Dr. Richard Arnold for general advice and entomological expertise.

Citation: Nagarajan RP, Goodbla A, Graves E, Baerwald M, Holyoak M, Schreier A (2020) Non-invasive genetic monitoring for the threatened valley elderberry longhorn beetle. PLoS ONE 15(1): e0227333. https://doi.org/10.1371/journal.pone.0227333

References

1. Fisher WS. A New Cerambycid Beetle from California. Proc Entomol Soc Washingt. 1921;23: 206–208.

2. Talley TS, Wright D, Holyoak M. Assistance with the 5-Year Review of the Valley elderberry longhorn beetle (Desmocerus californicus dimorphus). Rep to US Fish Wildl Serv Sacramento Fish Wildl Off Sacramento, Calif. 2006.

3. Holyoak M, Talley TS, Hogle SE. The effectiveness of US mitigation and monitoring practices for the threatened Valley elderberry longhorn beetle. J Insect Conserv. 2010.

4. Collinge SK, Holyoak M, Barr CB, Marty JT. Riparian habitat fragmentation and population persistence of the threatened valley elderberry longhorn beetle in central California. Biol Conserv. 2001. (00)00211–1

5. Barr C. The Distribution, Habitat, and Status of the Valley Elderberry Longhorn Beetle Desmocerus californicus dimorphus Fisher (Insecta: Coleoptera: Cerambycidae). Sacramento; 1991.

6. Burke HE. Biological notes on Desmocerus, a genus of roundhead borers, the species of which infest various elders. J Econ Entomol. 1921;14: 450–452.

7. Lang FJ, Jokerst JD, Sutter GE. HABITAT AND POPULATIONS OF THE VALLEY ELDERBERRY LONGHORN BEETLE ALONG THE SACRAMENTO RIVER 1.

8. Fish US and Service Wildlife. Framework for Assessing Impacts to the Valley Elderberry Longhorn Beetle (Desmocerus californicus dimorphus). 2017.

9. Fish US and Service Wildlife. Listing the Valley Elderberry Longhorn Beetle as a Threatened Species with Critical Habitat. Fed Regist. 1980;45.

10. Katibah EF, Dummer KJ, Nedeff NE. Current Condition of Riparian Resources in the Central Valley of California. California Riparian Systems: Ecology, Conservation, and Productive Management. 1984.

11. Fish US and Service Wildlife. Endangered and Threatened Wildlife and Plants; Withdrawal of the Proposed Rule To Remove the Valley Elderberry Longhorn Beetle From the Federal List of Endangered and Threatened Wildlife; Proposed Rule. Fed Regist. 2014;79.

12. Fish US and Service Wildlife. Endangered and Threatened Wildlife and Plants; Removal of the Valley Elderberry Longhorn Beetle From the Federal List of Endangered and Threatened Wildlife; Proposed Rule E:\FR\FM\02OCP4. Fed Regist. 2012;77.

13. USFWS. Withdrawal of the Proposal to Delist the Valley Elderberry Longhorn Beetle (VELB) Questions and Answers. 2014.

14. USFWS. Framework for Assessing Impacts to the Valley Elderberry Longhorn Beetle. 2017.

15. Talley T, Fleishman E, Holyoak M, Murphy D, Ballard a. Rethinking a rare-species conservation strategy in an urban landscape: The case of the valley elderberry longhorn beetle. Biol Conserv. 2007;135: 21–32.

16. Fish US and Service Wildlife. Valley Elderberry Longhorn Beetle (Desmocerus californicus dimorphus) 5-Year Review: Summary and Evaluation. 2006.

17. Linsley EG, Chemsak JA. Cerambycidae of North America. Part VI, No. 1. Taxonomy and classification of the subfamily Lepturinae. Univ Calif Publ Entomol. 1972;69: 1–138.

18. Chemsak JA. Illustrated revision of the Cerambycidae of North America (Vol II. Lepturinae). Wolfsgarden Books, Chino.; 2005.

19. US Fish and Wildlife Service. Recovery plan for the valley elderberry longhorn beetle. Portland, OR; 1984.

20. Halstead JA, Oldham J. Revision of the nearctic Desmocerus Audinet-Serville with emphasis on the federally threatened valley elderberry longhorn beetle (Coleoptera: Cerambycidae). Fresno, CA; 1990.

21. Ahrens D, Monaghan MT, Vogler AP. DNA-based taxonomy for associating adults and larvae in multi-species assemblages of chafers (Coleoptera: Scarabaeidae). Mol Phylogenet Evol. 2007. pmid:17420144

22. Raupach MJ, Astrin JJ, Hannig K, Peters MK, Stoeckle MY, Wägele JW. Molecular species identification of Central European ground beetles (Coleoptera: Carabidae) using nuclear rDNA expansion segments and DNA barcodes. Front Zool. 2010. pmid:20836845

23. Briski E, Cristescu ME, Bailey SA, MacIsaac HJ. Use of DNA barcoding to detect invertebrate invasive species from diapausing eggs. Biol Invasions. 2011.

24. Jackson JK, Battle JM, White BP, Pilgrim EM, Stein ED, Miller PE, et al. Cryptic biodiversity in streams: a comparison of macroinvertebrate communities based on morphological and DNA barcode identifications. Freshw Sci. 2014.

25. Taberlet P, Prud’Homme SM, Campione E, Roy J, Miquel C, Shehzad W, et al. Soil sampling and isolation of extracellular DNA from large amount of starting material suitable for metabarcoding studies. Mol Ecol. 2012. pmid:22300434

26. Rougerie R, Smith MA, Fernandez-Triana J, Lopez-Vaamonde C, Ratnasingham S, Hebert PDN. Molecular analysis of parasitoid linkages (MAPL): Gut contents of adult parasitoid wasps reveal larval host. Mol Ecol. 2011. pmid:21083857

27. Brandl S, Schumer G, Schreier BM, Conrad JL, May B, Baerwald MR. Ten real-time PCR assays for detection of fish predation at the community level in the San Francisco Estuary-Delta. Mol Ecol Resour. 2015. pmid:25042458

28. Zeale MRK, Butlin RK, Barker GLA, Lees DC, Jones G. Taxon-specific PCR for DNA barcoding arthropod prey in bat faeces. Mol Ecol Resour. 2011. pmid:21429129

29. Feinstein J. DNA sequence from butterfly frass and exuviae. Conserv Genet. 2004.

30. Lefort MC, Boyer S, Worner SP, Armstrong K. Noninvasive molecular methods to identify live scarab larvae: An example of sympatric pest and nonpest species in New Zealand. Mol Ecol Resour. 2012. pmid:22189059

31. Scriven JJ, Woodall LC, Goulson D. Nondestructive DNA sampling from bumblebee faeces. Mol Ecol Resour. 2013. pmid:23190789

32. Kumar S, Stecher G, Tamura K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol Biol Evol. 2016;33: 1870–1874. pmid:27004904

33. Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4: 406–25. pmid:3447015

34. Tamura K. Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G+C-content biases. Mol Biol Evol. 1992;9: 678–87. pmid:1630306

35. Felsenstein J. CONFIDENCE LIMITS ON PHYLOGENIES: AN APPROACH USING THE BOOTSTRAP. Evolution (N Y). 1985;39: 783–791. pmid:28561359

36. Thomsen PF, Elias S, Gilbert MTP, Haile J, Munch K, Kuzmina S, et al. Non-destructive sampling of ancient insect DNA. PLoS One. 2009. pmid:19337382

37. Nagarajan R, Schreier A, Baerwald M, Goodbla A, Graves E, Holyoak M. Using fecal DNA survey protocols to optimize traditional exit hole surveys and increase detection rates for VELB. Final Report to USFWS. 2018.

38. Roy K, Ewing CP, Hughes MA, Keith L, Bennett GM. Presence and viability of Ceratocystis lukuohia in ambrosia beetle frass from Rapid ʻŌhiʻa Death-affected Metrosideros polymorpha trees on Hawaiʻi Island. For Pathol. 2019.

39. Kamenova S, Mayer R, Rubbmark OR, Coissac E, Plantegenest M, Traugott M. Comparing three types of dietary samples for prey DNA decay in an insect generalist predator. Mol Ecol Resour. 2018. pmid:29509995

40. Ali B, Zhou Y, Zhang Q, Niu C, Zhu Z. Development of an easy and cost-effective method for non-invasive genotyping of insects. Kalendar R, editor. PLoS One. 2019;14: e0216998. pmid:31158235

41. Ide T, Kanzaki N, Ohmura W, Okabe K. Molecular identification of an invasive wood-boring insect Lyctus brunneus (Coleoptera: Bostrichidae: Lyctinae) using Frass by loop-mediated isothermal amplification and nested PCR assays. J Econ Entomol. 2016. pmid:27018442

42. Schrader C, Schielke A, Ellerbroek L, Johne R. PCR inhibitors—occurrence, properties and removal. Journal of Applied Microbiology. 2012. pmid:22747964

43. Lee AB, Cooper TA. Improved direct PCR screen for bacterial colonies: Wooden toothpicks inhibit PCR amplification. Biotechniques. 1995.

44. Schapker PA. A Review of the Lepturinae (Coleoptera: Cerambycidae) of the Pacific Northwest, with Special Attention to the Elderberry Longhorn Beetle, Desmocerus Dejean. Oregon State. 2014.

45. Baird NA, Etter PD, Atwood TS, Currey MC, Shiver AL, Lewis ZA, et al. Rapid SNP Discovery and Genetic Mapping Using Sequenced RAD Markers. Fay JC, editor. PLoS One. 2008;3: e3376. pmid:18852878

46. Ray AM, Swift IP, McElfresh JS, Alten RL, Millar JG. (R)-Desmolactone, A Female-produced Sex Pheromone Component of the Cerambycid Beetle Desmocerus californicus californicus (subfamily Lepturinae). J Chem Ecol. 2012. pmid:22350519

47. Ray AM, Arnold RA, Swift I, Schapker PA, McCann S, Marshall CJ, et al. (R)-desmolactone is a sex pheromone or sex attractant for the endangered valley elderberry longhorn beetle Desmocerus californicus dimorphus and several congeners (Cerambycidae: Lepturinae). PLoS One. 2014;9: 1–18. pmid:25521293