- Open Access
Analysing diet of small herbivores: the efficiency of DNA barcoding coupled with high-throughput pyrosequencing for deciphering the composition of complex plant mixtures
Frontiers in Zoology volume 6, Article number: 16 (2009)
In order to understand the role of herbivores in trophic webs, it is essential to know what they feed on. Diet analysis is, however, a challenge in many small herbivores with a secretive life style. In this paper, we compare novel (high-throughput pyrosequencing) DNA barcoding technology for plant mixture with traditional microhistological method. We analysed stomach contents of two ecologically important subarctic vole species, Microtus oeconomus and Myodes rufocanus, with the two methods. DNA barcoding was conducted using the P6-loop of the chloroplast trn L (UAA) intron.
Although the identified plant taxa in the diets matched relatively well between the two methods, DNA barcoding gave by far taxonomically more detailed results. Quantitative comparison of results was difficult, mainly due to low taxonomic resolution of the microhistological method, which also in part explained discrepancies between the methods. Other discrepancies were likely due to biases mostly in the microhistological analysis.
We conclude that DNA barcoding opens up for new possibilities in the study of plant-herbivore interactions, giving a detailed and relatively unbiased picture of food utilization of herbivores.
Small mammalian herbivores, such as voles and lemmings, play a key role in many boreal ecosystems where they function as the main link between the vegetation and predators [1, 2]. Their dramatic population fluctuations and consequent ecosystem implications [3, 4] have been subject to intensive studies for decades [5–7]. Trophic interactions have been emphasized as a main determinant of these fluctuations . While predator-prey interactions have been much emphasized , the interactions between voles and plants have gained new attention lately (see e.g. [9, 10]), as new interaction pathways have been identified [11–13]. However, plant-herbivore interactions cannot be understood without knowing, among other aspects, herbivore diets in natural settings. Yet most studies on plant-vole interactions have used only limited and indirect diet data [14–16], such as recording signs of rodent feeding on vegetation [17, 18] or cafeteria experiments in artificial settings [19, 20]. The most direct information comes from microhistological analysis of stomach content [21–24]. However, this is a very time-consuming method giving unspecific and context-dependent results , due to the small particle size as well as the composition of encountered plant material in rodent stomachs.
DNA barcoding, i.e. taxon identification using a standardized DNA region , is now increasingly used in ecological studies (see Valentini et al.  for a review), including diet analysis. DNA barcoding is particularly useful in diet determination when the food is not identifiable by morphological criteria, such as in the case of liquid feeders , or when the diet cannot be deduced by observing the feeding behavior, (e.g. diatom-feeding krill ). Diet analyses based on DNA markers concern so far mostly carnivorous animals (e.g. ). A universal approach for herbivorous diet analysis has been developed recently . This approach combines the new highly parallel sequencing systems , with the amplification of the P6 loop of the chloroplast trn L (UAA) intron .
In this paper, we present the first results on species-level food selection of two ecologically important vole species (Microtus oeconomus and Myodes (formerly Clethrionomys) rufocanus) in sub-arctic ecosystems based on DNA barcoding of ingested plants. We compare the results with the traditional microhistological method, which, for small rodents, has previously not been assessed against other methods.
The voles were sampled in low arctic tundra at Varanger peninsula (Finnmark, Norway: 70° 20' N, 30° 00' E) in July and September 2007. The sampling was conducted by snap trapping according to standard methods , after which the voles were dissected and their stomachs were stored in 70% ethanol for approximately 1/2 year prior to analysis. A sample of 48 individuals was chosen to be analyzed with both methods, based on the available space in DNA analysis batch and on stratifying according to species and season. The stomachs were dissected and the contents homogenized. Two subsamples were taken from each stomach for DNA analysis. They were stored in paper filter bags and submerged in silicagel to dry them. Remaining stomach contents were stored in 70% ethanol until the microhistological analysis.
Microhistological recognition of food particles in vole stomachs was based on leaf epiderm morphology. The shape of epidermal cells is taxon-specific, and several additional features, such as trichomes, hairs and characteristics of cells surrounding stomata, can be used for species identification [25, 35]. A photography-guide (Soininen & Nielsen, unpublished data) of epidermis of all vascular plants recorded at the sampling area (Ravolainen et al. unpublished data) was prepared using a method modified from Carrière . Dry plant samples were soaked overnight, scraped to reveal the epidermis and bleached with household bleach to clear the tissue of chlorophyll. Hard leaves were first boiled in table vinegar to soften the mesophyll tissue. Microphotographs (40×) were taken of abaxial and adaxial leaf side and leaf edge of all plants. Additional photographs were taken of stems and seeds of certain species of special interest. In addition to the specifically prepared epidermis photographs, photographs and microscopy slides of arctic plant epidermis, received from C. Hübner and E. Bjørkevoll, were used to aid identification.
A method modified from Hansson  was used for microscopy analysis. After taking subsamples for DNA analysis, stomach contents were filtered to > 0.16 mm and > 0.56 mm fractions. These were bleached with approximately 2 mL of household bleach for approximately 1/2 hour. One sample per fraction was analyzed, mounting a droplet of it on a microscopy slide. The frequency of occurrence of food items was recorded by light microscope (40×), by counting 25 hits on identifiable material along a measure grid. When approximately 95% of fragments were unidentifiable, the slide was discarded. For four individuals, no slide with adequate amount of identifiable material could be made. Therefore, they were discarded from the microhistological analysis.
Two samples were analyzed for each individual. Total DNA was extracted from about 10 mg of sample with the DNeasy Tissue Kit (Qiagen GmbH, Hilden, Germany), following the manufacturer's instructions. The DNA extracts were recovered in a total volume of 300 μ L. Mock extractions without samples were systematically performed to monitor possible contaminations.
DNA amplifications were carried out in a final volume of 25 μ L, using 2.5 μ L of DNA extract as template. The amplification mixture contained 1 U of AmpliTaq® Gold DNA Polymerase (Applied Biosystems, Foster City, CA), 10 mM Tris-HCl, 50 mM KCl, 2 mM of MgCl2, 0.2 mM of each dNTP, 0.1 μ M of each primer, and 0.005 mg of bovine serum albumin (BSA, Roche Diagnostic, Basel, Switzerland). The mixture was denatured at 95°C for 10 min, followed by 35 cycles of 30 s at 95°C, and 30 s at 55°C; the elongation was removed in order to reduce the +A artifact [37, 38]. Samples were amplified with using the universal primers g and h described by Taberlet et al. . The addition of a specific tag on the 5' end allowed an assignement of sequences to the respective samples. After amplification all samples were pooled for the pyrosequencing run. Each sample was recognized by a specific five bases long tag with at least two differences between tags for a better assignation of sequences to samples during bioinformatic segregation of sequences.
PCR products were purified using the MinElute PCR purification kit (Qiagen GmbH, Hilden, Germany). DNA quantification was carried out using the BioAnalyzer (Agilent Technologies, Inc., Santa Clara, CA). Taking these concentrations into account PCR products were pooled leading to equal amounts per sample. Large-scale pyrosequencing was carried out using GS FLX sequencer (Roche, Basel, Switzerland) following the manufacturer's instructions.
The first step of analyzing the output of the pyrosequencing consisted of sorting the different sequences according to the tag present on the 5' end of the primers. Thus, for each sample (each stomach content), a new file was generated, containing all the sequences having the relevant tag. Then, these sequences were analyzed to determine the diet. To limit the influence of sequence errors , only sequences that were present more than three times were considered in the subsequent analyses.
The sequences were compared to a database of 842 species representing all widespread and/or ecologically important taxa of the arctic flora (GenBank accession number GQ244527 to GQ245667) (Sønstebø et al: A minimalist DNA barcoding approach for reconstructing past Arctic vegetation and climate, submitted). It was developed by sequencing the whole chloroplast trn L (UAA) intron of these species using primer pair designed by Taberlet et al. , and following the protocol described and evaluated in Taberlet et al. . In the database a total of 33,5% of species and 77,1% of genera could be identified by the P6 loop. All families were unambiguously identified (Sønstebø et al: A minimalist DNA barcoding approach for reconstructing past Arctic vegetation and climate, submitted). When sequences were not fully identified using the arctic plant database, they were compared with sequences retrieved from GenBank, using ecoPCR ; http://www.grenoble.prabi.fr/trac/ecoPCR. The taxon was assigned to each sequence in a dataset by similarity assessment with a reference database using FASTA  algorithm. A FASTA alignment was retrieved if there was at least 98% of identity between query and database sequences and 100% of query coverage. If two or more taxa could be assigned with the same score for a given sequence, we assigned this sequence to the higher taxonomic level that included both taxa. This method resulted in some sequenced taxa being assigned to the rank of genus or family.
Chimeric sequences are a well know problem when amplifying a mixture of homologous genes, and it is impossible to avoid their formation . But if two unrelated taxa compose the chimeric sequence the resulting sequence is not taken into account because for taxon identification the FASTA alignments is retrieved only if the sequence have 100% of query coverage with the reference sequence. If two related taxa compose the chimeric sequence, this sequence is assigned to the higher taxonomic level that included both taxa (e.g. genus, family, order, etc.).
Results of the two methods were compared in two ways. First, the taxonomic resolution obtained by the two methods was examined by comparing the relative frequencies of hits (microscopy) and sequences (DNA) at different taxonomic levels. For this comparison, the four individuals with no microscopy data were excluded also from DNA dataset.
Second, the relative frequencies of food items in diets were compared between the methods. Taxonomic adjustments were first made to make the results from the two methods comparable. Nomenclature from species to family level follows Lid & Lid  and Elven , and for higher taxonomy Judd et al. . Several genera are represented only by one species in the study area and sequences assigned to these were therefore attributed to the respective species (e.g. Arctous alpinus, Bistorta vivipara, Chamaepericlymenum suecicum, Rumex acetosa). Similar adjustments were done at other taxonomic levels (e.g. Salicaceae to Salix and Ranunculales to Ranunculaceae). Then, proportions of food items estimated at the level of individual voles were averaged (mean ± standard deviation) across species and sampling season for both methods.
Using DNA barcoding, 75% of all sequences were identified at least to the genus level (Figure 1), whereas with the microhistological method, less than 20% of the identified fragments could be specified at this level. Consequently, more plant species and genera were identified in vole diets with the DNA barcoding than microhistology (Table 1 and 2). For the M. oeconomus diet, DNA barcoding identified 13 species and 9 genera compared with 9 and 5 with microscopy. Corresponding numbers for M. rufocanus were 17 and 8 (DNA barcoding) compared with 11 and 7 (microscopy).
Differences in taxonomic resolution between the methods made diet comparison complicated. For instance using microhistology, a large number of fragments could only be assigned to high taxonomic levels such as eudicotyledons and monocotyledons. Still, most of the species or genera that according to DNA barcoding were found to be quantitatively important in the diet, were also found to be so in the microhistological data, after differences in taxonomic accuracy were taken into account (Table 1 and 2). For example, for M. rufocanus both methods suggested that Vaccinium spp. are important food resources during both seasons, but especially in autumn (Table 2). Similarly, there was an agreement between the methods with respect to the importance of the family Polygonaceae (including Rumex acetosa and other species within this family) for both vole species in summer and that family Caryophyllaceae (with the species Stellaria nemorum and Cerastium spp.) is prevalent in M. oeconomus diet in the same season (Table 1 and 2). On the contrary, the substantial amounts of graminoids and horsetails (Equisetum spp.) found in the diet of M. oeconomus with microscopy were not evident from DNA analysis. Moreover, the prevalence of non-plant food items (fungi) and various plant structures (bark, root, seed) identified with microhistology could not be identified by DNA barcoding.
Although there was an agreement between the methods with respect to the importance of the main plant groups, DNA barcoding gave by far a taxonomically more detailed picture of the diet of the two vole species than did the microhistological analysis. Indeed, much of the discrepancy in the results derived from the two methods could be explained by differences in taxonomic resolution. Other discrepancies must be attributed to particular biases in the methods. Because a substantial fraction of the stomach content is left unidentified by the microhistological diet analysis, the proportion of the different plant groups in the diets will be biased towards the most easily identified groups by this method. One example might be the fairly large discrepancy in the abundance of monocotyledons between the methods. Another bias is due to varying epiderm/mesophyll ratios between taxa. Finally, a general problem is introduced by a great deal of subjectivity in the microhistological identification processes so that experience of the observer will have a major impact on the outcome of the analysis.
However, quantitative comparison between the methods should be done cautiously. Quantitative interpretation of DNA barcoding results is not straightforward; partly because it is based on chloroplast DNA. Thus, species from which chloroplast-rich tissues are eaten, are likely to be overrepresented compared to species mostly represented by e.g. seed or root in the diet [31, 46]. In addition, even if the trn L approach using g-h primers is universal for angiosperms and gymnosperms , the horsetails, mosses and fungi identified with microhistological method are mostly omitted by the DNA analysis. Together with easy microhistological identification this explains discrepancy between methods in the amount of Equisetum. Although DNA barcoding using the trn L approach also has its shortcomings, we conclude that it is superior to traditional methods for establishing diet analysis based on stomach content in small rodents. This novel method is less prone to the biases and context-dependencies that hamper microhistological analyses and yields a vastly improved taxonomic resolution. Information about diets at the level of individual plant species opens for testing more precise hypotheses on interactions between voles and plants. Furthermore, the new DNA-based technology makes it possible to study vole-plant interaction by non-destructive sampling of faeces in the natural habitats of voles. In faeces, DNA is highly degraded and only small fragments remain . The short sequences obtained using the trn L approach (10–143 base pairs) allow applying it on faecal samples . In this case, the first step will be to identify the rodent species using a mitochondrial DNA marker, and the second step will be the diet analysis. The analysis can even be more specific, by performing individual and sex identification using microsatellite polymorphism and Y-chromosome amplification [48, 49]. Thus, diet comparisons among species, individuals and sexes can be carried out, even without observing the animals (e.g. ).
For increasing the resolution in genera where sequences do not vary among the species (e.g. Carex and Salix), the trn L approach using the g-h primers can be complemented by one or several additional systems, specially designed for amplifying a short and variable region in these genera (as suggested in Valentini et al. ). Another strategy would be to amplify the whole trn L intron using universal primers such as c and d designed by Taberlet et al.  (254–767 bp), which strongly increases the resolution . However, even if this second strategy can be easily implemented for the analysis of stomach content (using the Titanium upgrade of the 454 FLX), it will not be suitable for faeces analysis because of the shortness of the degraded plant DNA fragments.
Krebs C, Danell K, Angerbjorn A, Agrell J, Berteaux D, Brathen K, Danell O, Erlinge S, Fedorov V, Fredga K, Hjalten J, Hogstedt G, Jonsdottir I, Kenney A, Kjellen N, Nordin T, Roininen H, Svensson M, Tannerfeldt M, Wiklund C: Terrestrial trophic dynamics in the Canadian Arctic. Can J Zool. 2003, 81: 827-843. 10.1139/z03-061.
Ims RA, Fuglei E: Trophic Interaction Cycles in Tundra Ecosystems and the Impact of Climate Change. Bioscience. 2005, 55: 311-322. 10.1641/0006-3568(2005)055[0311:TICITE]2.0.CO;2.
Gilg O, Hanski I, Sittler B: Cyclic Dynamics in a Simple Vertebrate Predator-Prey Community. Science. 2003, 302: 866-868. 10.1126/science.1087509.
Henden JA, Bårdsen BJ, Yoccoz RA, Nigel Gand Ims: Impacts of differential prey dynamics on the potential recovery of endangered arctic fox populations. J Appl Ecol. 2008, 45: 1086-1093.
Elton C: Periodic fluctuations in the numbers of animals: their causes and effects. J Exp Biol. 1924, 2: 119-163.
Hanski I, Henttonen H, Korpimaki E, Oksanen L, Turchin P: Small-Rodent Dynamics and Predation. Ecology. 2001, 82: 1505-1520.
Kausrud KL, Mysterud A, Steen H, Vik JO, Ostbye E, Cazelles B, Framstad E, Eikeset AM, Mysterud I, Solhoy T, Stenseth NC: Linking climate change to lemming cycles. Nature. 2008, 456 (7218): 93-U3. 10.1038/nature07442.
Berryman AA: Population cycles: the case of trophic interactions. 2002, Oxford: Oxford University Press
Rammul U, Oksanen T, Oksanen L, Lehtelä J, Virtanen R, Olofsson J, Strengbom J, Rammul I, Ericson L: Vole-vegetation interactions in an experimental, enemy free taiga floor system. Oikos. 2007, 116: 1501-1513.
Massey FP, Hartley SE: Experimental demonstration of the antiherbivore effects of silica in grasses: impacts on foliage digestibility and vole growth rates. Proc R Soc Lond B Biol Sci. 2006, 273: 2299-2304. 10.1098/rspb.2006.3586.
Conover M: Impact of the consumption of endophyte-infected perennial ryegrass by meadow voles. Agric Ecosyst Environ. 2003, 97: 199-203. 10.1016/S0167-8809(03)00037-9.
Huitu O, Helander M, Lehtonen P, Saikkonen K: Consumption of grass endophytes alters the ultraviolet spectrum of vole urine. Oecologia. 2008, 156: 333-340. 10.1007/s00442-008-0984-3.
Massey FP, Smith MJ, Lambin X, Hartley SE: Are silica defences in grasses driving vole population cycles?. Biol Lett. 2008, 4: 419-422. 10.1098/rsbl.2008.0106.
Laine KM, Henttonen H: Phenolics/nitrogen ratios in the blueberry Vaccinium myrtillus in relation to temperature and microtine density in Finnish Lapland. Oikos. 1987, 50: 389-395. 10.2307/3565500.
Hörnfeldt B: Long-term decline in numbers of cyclic voles in boreal Sweden: analysis and presentation of hypotheses. Oikos. 2004, 107: 376-392. 10.1111/j.0030-1299.2004.13348.x.
Högstedt G, Seldal T, Breistøl A: Period length in cyclic animal populations. Ecology. 2005, 86: 373-378. 10.1890/02-0561.
Kalela O: Regulation of reproduction rate in subarctic populations of the vole Clethrionomys rufocanus (Sund.). Ann Acad Sci Fenn. 1957, 34: 1-60.
Hambäck PA, Grellmann D, Hjältén J: Winter herbivory by voles during a population peak: the importance of plant quality. Ecography. 2002, 25: 74-80. 10.1034/j.1600-0587.2002.250108.x.
Myllymäki A: Interactions between the Field Vole Microtus agrestis and Its Microtine Competitors in Central-Scandinavian Populations. Oikos. 1977, 29: 570-580. 10.2307/3543595.
Hjalten , Joakim , Danell , Kjell , Ericson , Lars : Food Selection by Two Vole Species in Relation to Plant Growth Strategies and Plant Chemistry. Oikos. 1996, 76: 181-190. 10.2307/3545760.
Tast J: The food and feeding habits of the root vole, Microtus oeconomus, in Finnish Lapland. Aquilo Ser Zool. 1974, 15: 25-32.
Hansson L, Larsson TB: Vole diet on experimentally managed reforestation areas in northern Sweden. Holarct Ecol. 1978, 1: 16-26.
Batzli G, Jung HJG: Nutritional ecology of microtine rodents: resource utilization near Atkasook, Alaska. Arct Antarct Alp Res. 1980, 12 (4): 483-499.
Lindroth RL, Batzli GO: Food Habits of the Meadow Vole (Microtus pennsylvanicus) in Bluegrass and Prairie Habitats. J Mammal. 1984, 65: 600-606. 10.2307/1380843.
Carriére S: Photographic Key for the Microhistological Identification of some Arctic Vascular Plants. Arctic. 2002, 55: 247-268.
Hebert P, Ratnasingham S, deWaard J: Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc R Soc Lond B Biol Sci. 2003, 270 (Suppl 1): S96-S99. 10.1098/rsbl.2003.0025.
Valentini A, Pompanion F, Taberlet P: DNA barcoding for ecologist. Trends Ecol Evol. 2009, 9: 51-60.
Agusti N, Shayler S, Harwood J, Vaughan I, Sunderland K, Symondson W: Collembola as alternative prey sustaining spiders in arable ecosystems: prey detection within predators using molecular markers. Mol Ecol. 2003, 12 (12): 3467-3475. 10.1046/j.1365-294X.2003.02014.x.
Passmore AJ, Jarman SN, Swadling KM, Kawaguchi S, McMinn A, Nicol S: DNA as a dietary biomarker in antarctic krill, Euphausia superba. Mar Biotechnol. 2006, 8 (6): 686-696. 10.1007/s10126-005-6088-8.
Deagle BE, Gales NJ, Evans K, Jarman SN, Robinson S, Trebilco R, Hindel MA: Studying Seabird Diet through Genetic Analysis of Faeces: A Case Study on Macaroni Penguins (Eudyptes chrysolophus). PLoS ONE. 2007, 2: e831-10.1371/journal.pone.0000831.
Valentini A, Miquel C, Nawaz MA, Bellemain E, Coissac E, Pompanon F, Gielly L, Cruaud C, Nascetti G, Wincker P, Swenson JE, Taberlet P: New perspectives in diet analysis based on DNA barcoding and parallel pyrosequencing: the trnL approach. Mol Ecol Resour. 2009, 24 (2): 110-117.
Margulies M, Egholm M, Altman W, Attiya S, Bader J, Bemben L, Berka J, Braverman M, Chen Y, Chen Z, Dewell S, Du L, Fierro J, Gomes X, Godwin B, He W, Helgesen S, Ho C, Irzyk G, Jando S, Alenquer M, Jarvie T, Jirage K, Kim J, Knight J, Lanza J, Leamon J, Lefkowitz S, Lei M, Li J, Lohman K, Lu H, Makhijani V, McDade K, McKenna M, Myers E, Nickerson E, Nobile J, Plant R, Puc B, Ronan M, Roth G, Sarkis G, Simons J, Simpson J, Srinivasan M, Tartaro K, Tomasz A, Vogt K, Volkmer G, Wang S, Wang Y, Weiner M, Yu P, Begley R, Rothberg J: Genome sequencing in microfabricated high-density picolitre reactors. Nature. 2005, 437 (7057): 376-380.
Taberlet P, Coissac E, Pompanon F, Gielly L, Miquel C, Valentini A, Vermat T, Corthier G, Brochmann C, Willerslev E: Power and limitations of the chloroplast trnL (UAA) intron for plant DNA barcoding. Nucleic Acids Res. 2007, 35 (3): e14-10.1093/nar/gkl938.
Myllymäki A, Paasikallio A, Pankakoski E, Kanevo V: Removal experiments on small quadrats as a means of rapid assessment of the abundance of small mammals. Ann Zool Fenni. 1971, 8: 177-185.
Johnson MK, Wofford H, Pearson HA: Microhistological techniques for food habits analyses Volume 199 of Research paper SO. 1983, U.S. Department of Agriculture, Forest Service, Southern Forest Experimental Station
Hansson L: Methods of Morphological Diet Micreo-Analysis in Rodents. Oikos. 1970, 21: 255-266. 10.2307/3543682.
Brownstein M, Carpten J, Smith J: Modulation of non-templated nucleotide addition by tag DNA polymerase: Primer modifications that facilitate genotyping. Biotechniques. 1996, 20 (6): 1004-6.
Magnuson VL, Ally DS, Nylund SJ, Karanjawala ZE, Rayman JB, Knapp JI, Lowe AL, Ghosh S, Collins FS: Substrate nucleotide-determined non-templated addition of adenine by Taq DNA polymerase: implications for PCR-based genotyping and cloning. Biotechniques. 1996, 21: 700-9.
Huse SM, Huber JA, Morrison HG, Sogin ML, Mark Welch D: Accuracy and quality of massively parallel DNA pyrosequencing. Genome Biol. 2007, 8 (7): 10.1186/gb-2007-8-7-r143.
Taberlet P, Gielly L, Pautou G, Bouvet J: Universal primers for amplification of 3 noncoding regions of chloroplast DNA. Plant Mol Biol. 1991, 17 (5): 1105-1109. 10.1007/BF00037152.
Pearson WR, Lipman DJ: Improved tools for biological sequence comparison. Proc Natl Acad Sci USA. 1988, 85 (8): 2444-2448. 10.1073/pnas.85.8.2444.
Qiu XY, Wu LY, Huang HS, McDonel PE, Palumbo AV, Tiedje J, Zhou JZ: Evaluation of PCR-generated chimeras: Mutations, and heteroduplexes with 16S rRNA gene-based cloning. App Env Microbiol. 2001, 67: 880-887. 10.1128/AEM.67.2.880-887.2001.
Lid J, Lid D: Norsk flora. 2005, Oslo: Det norske samlaget
Elven RE: Checklist of the Panarctic Flora (PAF) Vascular Plants. Version: May 2007. 2007, [Http://www.binran.ru/infsys/paflist/index.htm]
Judd WS, Campbell CS, Kellogg EA, Stevens PF, Donoghue MJ: Plant systematics: a phylogenetic approach. 2008, Sunderland, Massachusetts U.S.A.: Sinauer Associates, inc, 3
Bradley BJ, Stiller M, Doran-Sheehy DM, Harris T, Chapman CA, Vigilant L, Poinar H: Plant DNA sequences from feces: Potential means for assessing diets of wild primates. Am J Primatol. 2007, 69 (6): 699-705. 10.1002/ajp.20384.
Deagle BE, Eveson JP, Jarman SN: Quantification of damage in DNA recovered from highly degraded samples – a case study on DNA in faeces. Front Zool. 2006, 3: 11-10.1186/1742-9994-3-11.
Griffits R, Tiwari B: Primers for the differential amplification of the sex-determining region y-gene in a range of mammal species. Mol ecol. 1993, 2 (6): 405-406. 10.1111/j.1365-294X.1993.tb00034.x.
Taberlet P, Luikart G: Non-invasive genetic sampling and individual identification. Biol J Linn Soc Lond. 1999, 68 (1–2): 41-55. 10.1111/j.1095-8312.1999.tb01157.x.
This study was supported by the Norwegian Research Council (project Ecosystem Finnmark), the EC FP6 ECOCHANGE project (Contract No 066866 GOCE) and Oskar Huttunen Foundation (personal scholarship for EMS). We thank Eirin Bjørkevoll, John-Andre Henden, Anne Grethe Hestnes, Christiane E. Hübner, Sissel Kaino, Siw T. Killengreen, Mette Nielsen and Virve Ravolainen for either supplying samples and data or assistance and expertice in microhistological analysis.
The authors declare that they have no competing interests.
EMS carried out the microhistological analysis, compared the results with DNA analysis and took the lead in writing the manuscript. AV, EC, CM and LG did the DNA analyses. CB, AKB and JHS provided the plant DNA database. RAI, NGY and PT developed the framework for the analyses. AV, LG, CB, AKB, RAI, NGY and PT assisted with the manuscript writing.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.