- Data note
- Open Access
Using next generation sequencing of alpine plants to improve fecal metabarcoding diet analysis for Dall’s sheep
BMC Research Notes volume 14, Article number: 173 (2021)
Dall’s sheep (Ovis dalli dalli) are important herbivores in the mountainous ecosystems of northwestern North America, and recent declines in some populations have sparked concern. Our aim was to improve capabilities for fecal metabarcoding diet analysis of Dall’s sheep and other herbivores by contributing new sequence data for arctic and alpine plants. This expanded reference library will provide critical reference sequence data that will facilitate metabarcoding diet analysis of Dall’s sheep and thus improve understanding of plant-animal interactions in a region undergoing rapid climate change.
We provide sequences for the chloroplast rbcL gene of 16 arctic-alpine vascular plant species that are known to comprise the diet of Dall’s sheep. These sequences contribute to a growing reference library that can be used in diet studies of arctic herbivores.
Dall’s sheep (Ovis dalli dalli) are endemic to alpine ecosystems of northwestern North America, and their populations have been declining in recent decades [1,2,3,4]. Climate change may be altering alpine plant communities and contributing to these declines. Dall’s sheep have a generalist plant diet; they were observed eating 110 different plant species in the Yukon Territory, Canada through traditional observational methods . However, the diet of Dall’s sheep remains relatively poorly characterized and represents a gap in understanding how climate change is affecting plant-animal interactions in alpine ecosystems.
The level of taxonomic resolution of items consumed in a diet study greatly affects ecological analysis . DNA based tools can infer diet composition with higher resolution and reduces cost, time, and effort compared to observational, morphological, and microhistological methods [7, 8]. Specifically, DNA metabarcoding uses universal primers for multispecies identification to mass-amplify DNA barcodes using PCR that are then read using next generation sequencing and assigned to the appropriate taxon . DNA barcoding includes a reference database of potential diet components, providing the capability to identify diet items to a desirable taxonomic resolution, ensuring that all components will be detected and assigned . Next generation sequencing of DNA from fecal samples has been successfully used to characterize diets of a variety of species, including ungulates [11, 12]. However, metabarcoding has not yet been used to assess the diet of Dall’s sheep. Lack of sequence data for some arctic/alpine plants known to be grazed upon by Dall’s sheep currently limits the development and application of metabarcoding for alpine herbivore diet studies.
To improve capabilities for diet analysis of Dall’s sheep and other arctic herbivores, we used a python script  to identify gaps in archived nucleotide sequence data for species known to comprise the diet of Dall’s Sheep, then obtained specimens of 16 species of arctic/alpine vascular plants for which sequence information was missing or underrepresented in publicly archived databases. We then sequenced the rbcL gene of the plant chloroplast genome, which is one of the most commonly used barcoding regions for plants [9, 14].
Plant specimens were obtained from herbarium specimens collected from the various arctic or alpine sites across mainland Alaska (Additional file 1). Plant tissue was extracted at the U. S. Geological Survey Alaska Science Center, employing a CTAB-PVP protocol modified from Stewart and Via  as reported by Muñiz-Salazar et al. . Extracts were quantified and shipped to the School of Environmental and Forest Sciences Genetics Lab at the University of Washington for PCR amplification and NexteraXT library preparation for sequencing. The rbcL gene region of each specimen was amplified via a two-step PCR protocol  with a primary amplification with tailed primers (rbcLaf + adaptor, rbcLr506 + adaptor) followed by a second round of amplification to anneal NexteraXT indices. Amplicons were quantified using a Qubit 4 Fluorometer (ThermoFisher) and diluted with dH2O to the recommended starting concentration for library preparation, 0.2 ng/μL (Illumina). Tagmentation, library amplification, and clean-up steps were completed according to the NexteraXT library preparation protocol (Illumina) with a variation of using New England Biolabs AMPure XP beads for cleanup instead of Agentcourt AMPure beads. The libraries were normalized and pooled prior to sequencing on an Illumina Miseq platform. Samples were paired-end sequenced in a 2 × 300 bp format .
Illumina sequence reads were processed in Geneious Prime 2020.2.4. Forward and reverse read files (fastq) were paired upon import, then quality trimmed with BBDuk trimmer (minimum quality 20, minimum overlap 20, minimum length 20). Sequences were normalized, then aligned and assembled using the de novo assembly tool (Geneious Prime). Assembled contigs were uploaded and annotated using BankIt, then submitted to GenBank  (Table 1).
The following are limitations for these data files:
We sequenced one DNA extraction from each plant species.
The sequencing project was funded through a grant to train new users on Illumina Nextera sequencing.
Large subunit of ribulose 1, 5 bisphosphate carboxylase/oxygenase (RUBISCO or RuBPCase)
DNA extraction method using cetyltrimethylammonium bromide as a detergent-based extraction buffer and polyvinylpyrrolidone, which is added to remove phenolic compounds from plant DNA extracts [15, 16]
Polymerase chain reaction
NexteraXT DNA library preparation kit enables sequencing of small genomes, PCR amplicons, and plasmids (Illumina)
Illumina Miseq Next Generation Sequencer is an integrated instrument that performs clonal amplification, genomic DNA sequencing, and data analysis with base calling, alignment, variant calling, and reporting in a single run (Illumina)
Alaska Department of Fish and Game (ADFG). Division of Wildlife Conservation, Wildlife Management Report (Juneau, AK: Alaska Department of Fish and Game) Trends in Alaska sheep populations, hunting, and harvests. 2014.
Koizumi CL, Carey J, Branigan M, Callaghan, K. Status of Dall’s’s sheep (Ovis Dall’si Dall’si) in the Northern Richardson Mountains. Yukon Fish and Wildlife Branch Report. 2011. TRC-11–01. Whitehorse, Yukon, Canada.
Rattenbury KL, Schmidt JH. Declining sheep populations in Alaska's Arctic Parks. Alaska Park Science 2017; 16: 67– 69. https://www.nps.gov/articles/aps-16-1-15.htm.
Rattenbury KL, Schmidt JH, Swanson DK, Borg BL, Mangipane BA, Sousanes PJ. Delayed spring onset drives declines in abundance and recruitment in a mountain ungulate. Ecosphere. 2018;9(11):e02513. https://doi.org/10.1002/ecs2.2513.
Hoefs M, Cowan I. Ecological investigation of a population of Dall’s sheep (Ovis Dall’si Dall’si Nelson). Syesis. 1979;12:1–81.
Smith MA, Eveleigh ES, McCann KS, Merilo MT, McCarthy PC, Van Rooyen KI. Barcoding a quanitified food web: crypsis, concepts, ecology and hypotheses. PLoS ONE. 2011;6:e14424.
Alberdi A, Aizpurua O, Bohmann K, Gopalakrishnan S, Lynggaard C, Nielsen M, et al. Promises and pitfalls of using high-throughput sequencing for diet analysis. Mol Ecol Resour. 2019;19:327–48. https://doi.org/10.1111/1755-0998.12960.
Andriollo T, Gillet F, Michaux JR, Reudi M. The menu varies with metabarcoding practices: a case study with bat Plecotus auritus 2019. PLoS ONE. 2019;14(7):e0219135.
Pompanon F, Deagle BE, Symondson WOC, Brown DS, Jarman SN. Who is eating what: diet assessment using next generation sequencing. Mol Ecol. 2012;21:1931–50. https://doi.org/10.1111/j.1365-294X.2011.05403.x.
Nielsen JM, Clare EL, Hayden B, Brett MT, Kratina P. Diet tracing in ecology: Method comparison and selection. Methods Ecol Evol. 2018;9:278–91.
Kartzinel TR, Chen PA, Coverdale TC, Erickson DL, Kress WJ, Kuzmina ML, et al. DNA metabarcoding illuminates dietary niche partitioning by African large herbivores. Proc Natl Acad Sci USA. 2015;112:8019–24. https://doi.org/10.1073/pnas.1503283112.
McShea WJ, Sukmasuang R, Erickson DL, Herrmann V, Ngoprasert D, Bhumpakphan N, et al. Metabarcoding reveals diet diversity in an ungulate community in Thailand. Biotropica. 2019;51:923–37. https://doi.org/10.1111/btp.12720.
Menning DM, Talbot S. Python scripts for bioinformatics, 2017. U.S. Geological Survey data release. 2018. https://doi.org/https://doi.org/10.5066/F74F1NZ4..
Hollingsworth ML, Clark AA, Forrest LL, Richardson J, Pennington RT, Long DG, et al. Selecting barcoding loci for plants: evaluation of seven candidate loci with species-level sampling in three divergent groups of land plants. Mol Ecol Resour. 2009;9:439–57. https://doi.org/10.1111/j.1755-0998.2008.02439.x.
Stewart CN, Via LE. A rapid CTAB DNA isolation technique useful for RAPD fingerprinting and other PCR applications. Biotechniques. 1993;14(5):748–50 (PMID: 8512694).
Muñiz-Salazar RM, Talbot SL, Sage GK, Ward DH, Cabello-Pasini A. Population genetic structure of annual and perennial populations of Zostera marina L. along the Pacific coast of Baja California and the Gulf of California. Mol Ecol. 2005;14:711–22.
de Vere N, Jones LE, Gilmore T, Moscrop J, Lowe A, Smith D, et al. Using DNA metabarcoding to investigate honey bee foraging reveals limited flower use despite high floral availability. Sci Rep. 2017;7:42838. https://doi.org/10.1038/srep42838.
Benson DA, Cavanaugh M, Clark K, Karsch-Mizrachi I, Lipman DJ, Ostell J, et al. GenBank. Nucleic Acids Res. 2013;41(Database issue):D36–42. https://doi.org/10.1093/nar/gks1195 (PMID: 23193287).
Williams K, et al. Elymus borealis rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW538513. 2021.
Williams K, et al. Gentiana propinqua rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW538515. 2021.
Williams K, et al. Juncus mertensianus rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW548523. 2021.
Williams K, et al. Luzula arctica rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW548524. 2021.
Williams K, et al. Ranunculus kamchaticus rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW548525. 2021.
Williams K, et al. Oxytropsis scammaniana rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW548526. 2021.
Williams K, et al. Packera ogotorukensis rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW548527. 2021.
Williams K, et al. Penstemon gormanii rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW548528. 2021.
Williams K, et al. Saxifraga caespitosa rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW548529. 2021.
Williams K, et al. Silene tayloriae rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW548530. 2021.
Williams K, et al. Smelowskia integrifolia rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW548531. 2021.
Williams K, et al. Stellaria alaskana rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW548532. 2021.
Williams K, et al. Taraxacum lyratum rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW548533 (2021).
Williams K, et al. Anemone lithophilia rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW526257. 2021.
Williams K, et al. Carex pyrenaica rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW538514. 2021.
Williams K, et al. Elymus latiglumis rbcl partial GenBank https://identifiers.org/ncbi/insdc:MW537582. 2021.
The Illumina team at the University of Washington provided valuable training on Illumina sequencing technology. Plant specimens were obtained from the U.S. Geological Survey and the University of Alaska, Anchorage herbarium. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.
Illumina and NASA’s Arctic and Boreal Vulnerability Experiment program (Grant NNX15AU21A to LRP). The United States Geological Survey and National Park Service provided funding in terms of salary and field and laboratory processes.
Ethics approval and consent to participate
Consent for publication
There are no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Williams, K.E., Menning, D.M., Wald, E.J. et al. Using next generation sequencing of alpine plants to improve fecal metabarcoding diet analysis for Dall’s sheep. BMC Res Notes 14, 173 (2021). https://doi.org/10.1186/s13104-021-05590-z
- Dall’s sheep