- Research Note
- Open access
- Published:
GoEnrich: creating high quality genomic DNA resources from limited voucher specimen tissues or museum specimens of at-risk species for conservation-friendly use in the validation of environmental DNA assays
BMC Research Notes volume 17, Article number: 255 (2024)
Abstract
Objective
Environmental DNA (eDNA) methods are crucial for monitoring populations, particularly rare and cryptic species. For confident eDNA application, rigorous assay validation is required including specificity testing with genomic DNA (gDNA). However, this critical step is often difficult to achieve as obtaining fresh tissue samples from at-risk species can be difficult, highly limited, or impossible. Natural history museum collections could serve as a valuable and ethical voucher specimen resource for eDNA assay validation. The present study demonstrates the effectiveness of whole genome amplification (WGA) in providing enough gDNA to assemble high quality mitogenomes from which robust targeted eDNA assays can be designed.
Results
Using fresh and historical museum tissue samples from six species spanning fish, birds, and mammals, we successfully developed a WGA method with an average yield of 380 to 1,268 ng gDNA per 20 µL reaction. This gDNA was used for whole genome shotgun sequencing and subsequent assembly of high quality mitogenomes using mtGrasp. These mitogenomes were then used to develop six new robust, targeted quantitative real time polymerase chain reaction-based eDNA assays and 200 ng WGA-enriched yielded satisfactory Cq values and near 100% detection frequencies for all assays tested. This approach offers a cost-effective and non-invasive alternative, streamlining eDNA research processes and aiding in conservation efforts.
Introduction
Environmental DNA (eDNA), genetic material obtained from environmental samples, has enhanced biodiversity monitoring offering a sensitive and cost-effective means of detecting cryptic species compared to traditional methods [1, 2]. However, data quality and reliability issues have undermined trust in eDNA applications [3, 4]. Validating targeted qPCR-based eDNA assays requires ensuring assay specificity using genomic DNA (gDNA) from target and non-target species; however, obtaining tissue samples for gDNA isolation is challenging for elusive or endangered species due to permit processes [3,4,5]. Researchers have explored natural history museums as potential sources of DNA [6], extracting gDNA from voucher specimens, though optimal methods for gDNA isolation and processing remain underexplored due to variations in specimen preparation and archival conditions [7].
Museum specimens are commonly stored under ambient conditions, either desiccated or in liquid-preserved collections, with some housed in refrigerated cold rooms. Over time, prolonged exposure to these conditions causes long strands of high molecular weight DNA to degrade into shorter fragments, diminishing both the size and quantity of DNA [8]. A substantial portion of dried specimens, often taxidermized hides, pose challenges in preserving intact DNA sequences suitable for downstream molecular analysis [9]. Furthermore, DNA degradation can occur when specimens are preserved in alcohol, typically 70–75% ethanol, where decreasing alcohol concentrations due to evaporation may lead to DNA hydrolysis [10]. Fragmented DNA and DNA that is chemically crosslinked by formaldehyde-based preservatives complicate downstream analysis such as qPCR due to the limited availability of intact gDNA [11].
Recent advances in museum sample genomic studies (museomics) have employed whole genome amplification (WGA) to amplify genomic DNA in an unbiased way [12]. However, these techniques were not tested for enhancing genomic resources for validating qPCR-based eDNA assays. The present study presents a workflow utilizing WGA-enriched gDNA from fresh and historical museum samples for genome skimming and mitochondrial genome assembly. Additionally, we demonstrate the use of enriched gDNA to validate six targeted qPCR-based assays.
Main text
Methods
Tissue collection and extraction
The present study includes freshly collected and historical museum tissue samples. The goldeye (Actinopterygii: Hiodon alosoides, te-HIAL) and cisco (Actinopterygii: Coregonus artedi, te-COAR) fresh tissues were obtained from the Alberta Ministry of Environment and Protected Areas. Fresh tissue samples of Northern goshawk (Aves: Accipiter gentilis, av-ACGE) were sourced from the Beatty Biodiversity Museum. Fresh tissue samples were preserved in 70% ethanol and stored at -20 °C until DNA extraction. Museum voucher specimens, comprising dried tissue samples of snowshoe hare (Mammalia: Lepus americanus, ma-LEAM), bison (Mammalia: Bison bison, ma-BIBI), and bald eagle (Aves: Haliaeetus leucocephalus, av-HALE), were obtained from the Royal British Columbia Museum in Victoria, Canada. Possession and disposal of these wildlife tissue samples for genetic analysis are covered under British Columbia Wildlife Act permit #NA22-782623. DNA extraction was performed using the DNeasy Blood and Tissue kit (Qiagen, Ontario, Canada), and quantification was conducted using the Qubit™ 1X dsDNA High Sensitivity Assay Kit (Invitrogen, Massachusetts, USA). Further details on the samples used in the present study are provided in Table 1.
Whole genome amplification (WGA) of genomic DNA
A modified method by Hutchison et al. [13] was adapted in the present study. A master mix was created by combining 20 ng gDNA, 1 µL 100 µM random hexamer primers (Cat. R016 and R106, Fidelity Systems, Maryland, USA), and 0.5 µL Phi29 10X buffer (Cat. M0269S, New England Biolabs, Ontario, Canada), adjusted to 5 µL with UltraPure-dH2O (Invitrogen, Massachusetts, USA). The mixture underwent heating at 95 °C for 5 min, then quick cooling by plunging on ice to 4 °C. This mix was added to a mixture of 1.5 µL Phi29 10X buffer, 1 µL Phi29 DNA Polymerase (10,000 U/mL, New England Biolabs), 8 µL 2.5 mM dNTPs (FroggaBio, Ontario, Canada), and 4.5 µL of UltraPure-dH2O, resulting in a 20 µL final reaction volume. Incubation at 30 °C for 18 h followed, with subsequent deactivation of Phi29 polymerase at 65 °C for 10 min in a BioRad thermocycler. WGA products were purified using the QIAquick PCR purification kit (Cat. 28104, Qiagen, Ontario, Canada), and a portion was subjected to electrophoresis on a 0.8% agarose gel and visualized with Gel Red (Gold Biotechnology, Missouri, USA) for DNA size and integrity confirmation. DNA sizes were compared with a 1Kb DNA Ladder RTU (GeneDireX, Taipei, Taiwan).
Library preparation, whole genome shotgun sequencing, and mitogenome assembly
The WGA-enriched gDNA samples were sent to the Michael Smith Genome Sciences Centre, BC Cancer, Vancouver, Canada for library preparation and whole genome sequencing. The TruSeq DNA PCR-free library prep kit (Cat. 20015963, Illumina Inc., California, US) was used for library preparation using the Microlab NIMBUS liquid handling robot (Hamilton Robotics, Nevada, USA). A total of 500 ng gDNA was fragmented by sonication with Covaris LE220 (Covaris, Massachusetts, USA) and underwent end repair and size selection of 300–400 bp fragments using PCRClean DX paramagnetic beads (Cat. C-1003-50, Aline Biosciences, Massachusetts, USA). After 3’ A-tailing, TruSeq adapters (Cat. 20015963, Illumina Inc., California, USA) were added and purified with PCRClean DX paramagnetic beads. Quantification of PCR-free genome libraries was done using a qPCR Library Quantification kit (Cat. 07960140001, KAPA, Massachusetts, USA) before sequencing on an Illumina NovaSeq 6000 instrument with S4 reagents in a single run, generating 150 bp paired end reads (Illumina Inc., California, USA). Mitogenomes were assembled, circularized, and standardized using mtGrasp (Mitochondrial Reference-Grade Genome Assembly and Standardization Pipeline, v. 1.1.0), aligning sequences from 5’ to 3’ end with tRNA-Phe as the initial sequence [14, 15].
eDNA assay design and validation
The eDNA assays for detecting L. americanus (eLEAM2), H. leucocephalus (eHALE), and C. artedii (eCOAR7) were fully validated through prior studies [15, 19], while the rest of the assays were validated as part of the current work (Table S2). All eDNA assays were designed and validated based on our established workflow [5], with performance characteristics meeting or exceeding Canadian standards [3, 4]. Mitogenomes of target and confounding species (from the National Center for Biotechnology Information (NCBI); Table S1) were aligned using MAFFT v7.490 [16], and phylogenetic trees were generated with RAXml-NG v. 1.2.2 [17]. Unikseq was used to identify unique regions in the target species’ whole mitogenome sequences to increase the likelihood of finding suitable regions for robust qPCR assay design [18]. Primers and probes were selected using Beacon Designer 8.21 (PREMIER Biosoft, California, USA) based on identified regions. Successful primer and probe sequences for each eDNA assay are listed in Tables S2 and S3.
In vitro specificity validation involved testing 10 ng WGA-enriched target species’ gDNA (2 µL of 5 ng/µL gDNA per 20 µL reaction) and non-target species’ gDNA following methodologies outlined in previous studies [5, 18, 19]. qPCR validation of multiple primer pairs tested on gDNA from both target and sympatric species was done using QIAcuity EG PCR Kit (Cat. 250111, Qiagen). The thermocycler profile included an initial denaturation at 95 °C for 2 min, followed by 50 cycles of 15 s at 95 °C, 30 s at 64 °C, and 45 s at 72 °C. High-end specificity validation of primer-probe assay was carried out using QIAcuity Probe PCR Kit (Cat. 250101, Qiagen) with 25 technical replicates per sample. Amplification was observed in only a few replicates with the 10 ng WGA-enriched gDNA from older museum samples as input. To enhance amplification frequency, subsequent tests utilized higher input amounts of 100, 200, and 1,000 ng per reaction. Lastly, serial dilutions of synthetic double-stranded DNA (gBlocks®, IDT, Iowa, USA) were prepared to construct standard curves (Figures S1 – S6) to assess assay sensitivity [20].
Statistical analysis
To minimize the impact of outliers, the data is presented as median Cq values. Friedman repeated measures test was employed to identify significant differences (p ≤ 0.05) among the test groups. Within each group, pairwise comparisons between treatments were conducted using the Wilcoxon Signed Rank test to determine statistical significance (p ≤ 0.05).
Results
WGA enrichment of fresh and museum gDNA samples
After WGA-enrichment, the sample concentrations ranged from 19 to 63.4 ng/µL with overall yields between 380 and 1,268 ng per 20 µL reaction. Changes in the fragment size of gDNA samples before and after WGA enrichment were observed (Fig. 1). The gDNA samples extracted from fresh tissues showed less smearing, whereas museum voucher specimens exhibited significant fragmentation, ranging from 100 to 1,000 bp. The WGA-enriched gDNA displayed enhanced fragment size, with a smear observed at 3,000 bp and larger.
Assembled mitogenomes
Libraries were constructed from WGA-enriched gDNA and resulted in 1.92 to 2.3 million reads (1.9 ± 0.4 million) per library. We successfully constructed mitogenomes for all six species examined in the present study (Table 1). With mtGrasp, we successfully assembled a single contig for each sample, suggesting a high sequencing depth of the mitogenome. The assembled mitogenomes ranges from 16,197 to 18,645 base pairs.
WGA-enriched gDNA samples for eDNA assay specificity validation
The assays developed in the present study were thoroughly validated for their specificity by testing them against the gDNA of the target species, sympatric species in the same environment, and common contaminating DNA such as human, cat, and dog. There was no cross-amplification of non-target species for any of the six eDNA assays (Table S4).
With limited access to voucher specimens, being able to confidently use WGA-enriched gDNA for assay validation enables individual labs to confidently demonstrate that an eDNA assay is performing as expected. We compared the abilities of the non-enriched and WGA-enriched gDNA to effectively act as templates to amplify amplicons. Utilizing 10 ng/reaction non-enriched and WGA-enriched gDNA from fresh tissue samples (Fig. 2A) resulted in comparable Cq values and detection frequencies for C. artedi (te-COAR) and H. alosoides (te-HIAL), but a 14 Cq increase for A. gentilis (av-ACGE) WGA-enriched DNA compared to non-enriched gDNA was observed.
In the case of museum samples (Fig. 2B), increasing the input amounts of WGA-enriched gDNA improved Cq values and detection frequencies. Specifically, 10 ng/reaction WGA-enriched gDNA had very high Cq values (43 to 49) and reduced detection frequencies (12 to 64%) compared to 10 ng/reaction non-enriched gDNA. Increasing the concentration of WGA-enriched gDNA to 200–2000 ng/reaction lowered the Cq values while greatly improving detection frequency to near-100% among replicates for all assays tested (Fig. 2B). All eDNA assays in the present study meet the sensitivity criteria set by the Canadian Standard Association (Table S5).
Discussion
The present study shows the efficacy of WGA enrichment in generating complete mitogenomes from gDNA samples extracted from limited fresh or museum voucher specimens. Moreover, the enriched gDNA can serve as a replenishable source of genetic material for eDNA assay validation. Providing an accessible way to grow the mitogenome sequence repositories will enhance the development of more reliable qPCR-based eDNA tools [5]. Unbiased selection of unique sequences of the target taxon rather than relying on short, barcoded regions of the mitogenome increases the likelihood of designing a taxon-specific assay while satisfying the rules of good primer and probe design for qPCR [15].
In vitro specificity assessments necessitate the collection of gDNA from target and sympatric taxa, sourced from voucher specimens [5, 21]. Acquiring tissue samples for gDNA isolation, especially from elusive, endangered, and regulated species, often entails cumbersome permit procedures and can negatively impact animal health. Leveraging natural history museums as repositories of gDNA resources can circumvent these challenges. WGA-enriched gDNA facilitates the preservation of requisite samples, bolstering their concentration and volume through unbiased amplification [7].
Limitations
We used six separate targeted qPCR assays and evaluated the specificity of the primers and probes using WGA-enriched gDNA. More assays on additional specimens, including even older museum specimens, is desired.
Conclusion
WGA is an innovative approach that effectively amplifies the original template DNA, enabling the production of a substantial quantity of genetic information from small quantities of gDNA. The present study demonstrates the successful application of WGA enrichment to generate complete mitogenomes using small amounts of gDNA from endangered species and museum voucher specimens. WGA-enriched gDNA samples were also useful in the comprehensive validation of eDNA assay specificity by increasing the concentration and volume of the original gDNA samples needed for the validation. This is particularly advantageous in minimizing the need for continuous destructive sampling of museum samples, arduous permit applications, and expensive and intrusive collection of rare and endangered species.
Data availability
The datasets supporting the conclusions of this article are included within the article and its additional files. Mitochondrial genomes assembled in the present study were uploaded to National Center for Biotechnology Information’s (NCBI) nucleotide database (see Table 1 for accession numbers).
Abbreviations
- av-ACGE:
-
Northern goshawk (Aves: Accipiter gentilis)
- av-HALE:
-
Bald eagle (Aves: Haliaeetus leucocephalus)
- eDNA:
-
Environmental deoxyribonucleic acid
- gDNA:
-
Genomic deoxyribonucleic acid
- ma-BIBI:
-
Bison (Mammalia: Bison bison)
- ma-LEAM:
-
Snowshoe hare (Mammalia: Lepus americanus)
- qPCR:
-
Quantitative real-time polymerase chain reaction
- te-HIAL:
-
Goldeye (Actinopterygii: Hiodon alosoides)
- te-COAR:
-
Cisco (Actinopterygii: Coregonus artedi)
- WGA:
-
Whole genome amplification
References
Deiner K, Bik HM, Mächler E, Seymour M, Lacoursière-Roussel A, Altermatt F, Creer S, Bista I, Lodge DM, Pfrender ME, Bernatchez L. Environmental DNA metabarcoding: transforming how we survey animal and plant communities. Mol Ecol. 2017;26(21):5872–95. https://doi.org/10.1111/mec.14350.
Goldberg CS, Turner CR, Deiner K, Klymus KE, Thomsen PF, Murphy MA, Spear SF, McKee A, Oyler-McCance SJ, Cornman RS, Laramie MB, Mahon AR, Lance RF, Pilliod DS, Strickler KM, Waits LP, Fremier AK, Takahara T, Herder JE, Taberlet P. Critical considerations for the application of environmental DNA methods to detect aquatic species. Methods Ecol Evol. 2016;7(11):1299–307. https://doi.org/10.1111/2041-210X.12595.
Gagné N, Bernatchez L, Bright D, Côté G, Coulson M, Gurney K, Hanner R, Helbing C, Hobbs J, Hocking M, Khan I, Naumann C, Parent G, Richter C, Silverio C, Skinner M, Weir A, Wilcox T, Wilson C, Clogg-Wright K. (2021). Environmental DNA (eDNA) reporting requirements and terminology (p. 31). National standard of Canada, CSA W214:21, Canadian Standards Association. https://www.csagroup.org/store/product/CSA%20W214%3A21/
Abbott C, Bright D, Bryant H, Côté G, Crookes S, Gurney K, Hanner R, Helbing C, Hocking M, Khan I, Langlois VS, Lemay M, Marshall N, Miliano R, Mirabzadeh-Ardakani A, Parent G, Richter C, Wagener A, Wilson C, Clogg-Wright. K. Performance criteria for the analyses of environmental DNA by targeted quantitative polymerase chain reaction. National standard of Canada, CSA W219:23. Canadian Standards Association,. 2023. 29 pp. https://www.csagroup.org/store/product/CSA_W219%3A23/
Langlois VS, Allison MJ, Bergman LC, To TA, Helbing CC. The need for robust qPCR-based eDNA detection assays in environmental monitoring and species inventories. Environ DNA. 2021;3(3):519–27. https://doi.org/10.1002/edn3.164.
Fortes GG, Paijmans JLA. Analysis of whole mitogenomes from ancient samples. In: Kroneis T, editor. Whole genome amplification. Methods in Molecular Biology. Volume 1347. New York, NY: Humana; 2015. https://doi.org/10.1007/978-1-4939-2990-0_13.
Raxworthy CJ, Smith BT. Mining museums for historical DNA: advances and challenges in museomics. Trends Ecol Evol. 2021;36(11):1049–60. https://doi.org/10.1016/j.tree.2021.07.009.
Camacho-Sanchez M, Burraco P, Gomez-Mestre I, Leonard JA. Preservation of RNA and DNA from mammal samples under field conditions. Mol Ecol Resour. 2013;13(4):663–73. https://doi.org/10.1111/1755-0998.12108.
Mundy N, Unitt P, Woodruff DS. Skin from feet of museum specimens as a non-destructive source of DNA for avian genotyping. Auk. 1997;114(1):126–9. https://doi.org/10.2307/4089075.
McGuire JA, Cotoras DD, O’Connell B, Lawalata SZS, Wang-Claypool CY, Stubbs A, Huang X, Wogan GOU, Hykin SM, Reilly SB, Bi K, Riyanto A, Arida E, Smith LL, Milne H, Streicher JW, Iskandar DT. Squeezing water from a stone: high-throughput sequencing from a 145-year old holotype resolves (barely) a cryptic species problem in flying lizards. PeerJ. 2018;6:e4470. https://doi.org/10.7717/peerj.4470.
Tsai WL, Schedl ME, Maley JM, McCormack JE. More than skin and bones: comparing extraction methods and alternative sources of DNA from avian museum specimens. Mol Ecol Resour. 2020;20(5):1220–7. https://doi.org/10.1111/1755-0998.13077.
Wang X, Liu Y, Liu H, Pan W, Ren J, Zheng X, Tan Y, Chen Z, Deng Y, He N, Chen H, Li S. Recent advances and application of whole genome amplification in molecular diagnosis and medicine. MedComm. 2022;3(1):e116. https://doi.org/10.1002/mco2.116.
Hutchison CA, Smith HO, Pfannkoch C, Venter JC. (2005). Cell-free cloning using φ29 DNA polymerase. Proceedings of the National Academy of Sciences, 102(48), 17332–17336. https://doi.org/10.1073/pnas.0508809102
Yang C, Coombe L. bcgsc/mtGrasp: mtGrasp v1.1.0 (v1.1.0). Zenodo. 2023. https://doi.org/10.5281/zenodo.10073719.
Lopez MLD, Yang CL, Coombe L, Warren RL, Allison MJ, Imbery JJ, Birol I, Helbing CC. (2024). Supporting data for: mtGrasp: Streamlined mitochondrial genome reference-grade assembly and standardization to enhance mitogenome resources and improve the development of environmental DNA assays [Data set]. Zenodo. https://doi.org/10.5281/zenodo.11375303.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30(4):772–80. https://doi.org/10.1093/molbev/mst010.
Stamatakis A. Bioinformatics. 2014;30(9):1312–3. https://doi.org/10.1093/bioinformatics/btu033. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies.
Allison MJ, Warren RL, Lopez ML, Acharya-Patel N, Imbery JJ, Coombe L, Yang CL, Birol I, Helbing CC. Enabling robust environmental DNA assay design with unikseq for the identification of taxon-specific regions within whole mitochondrial genomes. Environ DNA. 2023;5(5):1032–47. https://doi.org/10.1002/edn3.438.
Lopez MLD, Bonderud M, Allison MJ, MacDermid F, Ussery EJ, McMaster ME, Dersch A, Staniszewska KJ, Cooke CA, Drevnick P, Helbing CC. qPCR-based eDNA Workflow for humic-rich lake sediments: combined use of sedimentary DNA (sedDNA) and indigenous knowledge in reconstructing historical fish records. Ecol Ind. 2023;155:111014. https://doi.org/10.1016/j.ecolind.2023.111014.
Lesperance ML, Allison MJ, Bergman LC, Hocking MD, Helbing CC. A statistical model for calibration and computation of detection and quantification limits for low copy number environmental DNA samples. Environ DNA. 2021;3(5):970–81. https://doi.org/10.1002/edn3.220.
Thalinger B, Deiner K, Harper LR, Rees HC, Blackman RC, Sint D, Traugott M, Goldberg CS, Bruce K. A validation scale to determine the readiness of environmental DNA assays for routine species monitoring. Environ DNA. 2021;3(4):823–36. https://doi.org/10.1002/edn3.189.
Acknowledgements
Not applicable.
Funding
This work was funded under the Genome Canada, Genome British Columbia, Genome Québec large-scale applied research project #312ITD awarded to CCH. MLDL is supported by a Liber Ero postdoctoral fellowship. The funders had no role in study design; in the collection, analysis, and interpretation of data; in the writing of the report; and in the decision to submit the article for publication.
Author information
Authors and Affiliations
Contributions
MLDL, MTB, CCH designed the research; MLDL, MTB, IGM performed the experiments; MLDL, VCT, CCH analyzed the data; MLDL, MTB, IGM, VCT, CCH prepared the manuscript. MLDL, MTB, IGM, VCT, CCH read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Ethics approval from the University of Victoria’s Animal Care Committee was not required for the present study according to the Canadian Council on Animal Care as the obtained samples were opportunistically subsampled from existing provincial government and museum collections. Possession and disposal of these wildlife tissue samples for genetic analysis are covered under British Columbia Wildlife Act permit #NA22-782623.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
About this article
Cite this article
Lopez, M.L.D., Bonderud, M.T., Ma, I.G. et al. GoEnrich: creating high quality genomic DNA resources from limited voucher specimen tissues or museum specimens of at-risk species for conservation-friendly use in the validation of environmental DNA assays. BMC Res Notes 17, 255 (2024). https://doi.org/10.1186/s13104-024-06936-z
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13104-024-06936-z