DNA metabarcoding-based diet survey for the Eurasian otter (Lutra lutra): Development of a Eurasian otter-specific blocking oligonucleotide for 12S rRNA gene sequencing for vertebrates
Authors:
Priyanka Kumari aff001; Ke Dong aff001; Kyung Yeon Eo aff003; Woo-Shin Lee aff004; Junpei Kimura aff005; Naomichi Yamamoto aff001
Authors place of work:
Department of Environmental Health Sciences, Graduate School of Public Health, Seoul National University, Seoul, South Korea
aff001; Institute of Health and Environment, Graduate School of Public Health, Seoul National University, Seoul, South Korea
aff002; Conservation and Research Center, Seoul Zoo, Gwacheon, South Korea
aff003; Department of Forest Sciences, College of Agriculture and Life Science, Seoul National University, Seoul, South Korea
aff004; College of Veterinary Medicine, Seoul National University, Seoul, South Korea
aff005
Published in the journal:
PLoS ONE 14(12)
Category:
Research Article
doi:
https://doi.org/10.1371/journal.pone.0226253
Summary
The Eurasian otter (Lutra lutra) is an endangered species for which diet analyses are needed as part of its conservation efforts. Eurasian otters feed on vertebrates, such as fishes, and invertebrates, such as crustaceans, but their detailed taxonomies are not fully understood in part due to limited resolving power of traditional morphological identification methods. Here, we used high-throughput sequencing (HTS)-based DNA metabarcoding approaches to analyze diet profiles of Eurasian otters inhabiting a marshy estuary area in Korea. We investigated their diet profiles based on spraint sampling followed by DNA metabarcoding analyses targeting 12S rRNA gene region for vertebrates, 16S rRNA gene region for invertebrates, and cytochrome c oxidase 1 (COI) gene region for fishes. For the vertebrate analysis, a blocking oligonucleotide (OBS1) was designed to suppress amplification of DNA fragments derived from the otters. The 12S rRNA gene sequencing assay detected species belonging to fishes (95%) and amphibians (3.3%). Fishes detected by 12S rRNA gene sequencing included crucian carp (Carassius auratus), mullets (Mugil spp.), bluegill (Lepomis macrochirus), and northern snakehead (Channa argus), which were also detected by COI gene sequencing. Among invertebrates, mud flat crabs (Helicana spp.) and shrimps (Palaemon spp.) were abundant. The designed blocking oligonucleotide OBS1 effectively inhibited amplification of the otter’s DNA, with only up to 0.21% of vertebrate sequence reads assigned to the otter. This study demonstrated that HTS-based DNA metabarcoding methods were useful to provide in-depth information regarding diet profiles of the otters at our sampling site. By using HTS-based DNA metabarcoding approaches, future research will explore detailed taxonomies of their diets across locations and seasons.
Keywords:
DNA – Invertebrates – Vertebrates – Sequence databases – Ribosomal RNA – Otters – Gene sequencing – Oligonucleotides
Introduction
Eurasian otters (Lutra lutra) are widely distributed from Europe to Asia, including South Korea. Due to its declining population, however, the Eurasian otter is listed as a near threatened species by the International Union for Conservation of Nature Red List [1] and as a class I endangered species by the Ministry of Environment of Korea [2]. Eurasian otters feed on fishes, amphibians, crabs, birds, and insects [3, 4], and are at risk of intake of polluted diets, which may lead to a further decline in their population [5]. Therefore, it is critical to survey their diets, including potentially contaminated preys, as an effort for their conservation.
The otters’ feces (spraints) are useful for their diet analyses [3, 4, 6] and individual identification [7, 8]. Morphological observation of undigested food remains, such as bones, shells, and hair, in spraints provides information regarding the otters’ diets [3, 6]. Previous studies reported fishes as the otters’ main preys followed by amphibians, crabs, birds, and insects [3, 4]. Nonetheless, knowledge is limited regarding their detailed taxonomies in part due to limited resolving power of traditional morphological methods [9]. Meanwhile, Hong et al. [10] used a Sanger sequencing-based approach to identify vertebrate species for each individual bone remain isolated from spraints. However, this approach is laborious and requires technical expertise that limits the capacity of information generation.
Emerging high-throughput sequencing (HTS)-based DNA metabarcoding methods can circumvent such difficulties, and are used for diet analyses of wild animals, such as leopard cat in Pakistan [11], puma, jaguar, ocelot, and crab-eating fox in Venezuela [12], and Iberian lynx in Spain [13]. Due to its large capacity of information generation, the HTS-based methods can improve taxonomic resolution in greater details with limited quantities of fecal analytes [11, 14, 15]. However, the DNA metabarcoding-based diet analyses can potentially suffer from amplification of DNA fragments derived from a predator rather than preys, resulting in the reduced sensitivity in rare prey detection [11]. The most effective solution to overcome this problem is to use blocking oligonucleotides that suppress amplification of predator DNA [16]. We expect that a similar approach would be useful for diet analyses of the Eurasian otters inhabiting Korea.
Here, we investigated vertebrate and invertebrate diet profiles of the Eurasian otters inhabiting a marshy estuary area of South Korea based on spraint sampling followed by HTS-based DNA metabarcoding analyses targeting 12S rRNA gene region for vertebrates, 16S rRNA gene region for invertebrates, and cytochrome c oxidase 1 (COI) gene for fishes. We used previously reported universal primers for identification of vertebrates [17], invertebrates [15], and fishes [18]. In particular, we aimed to develop a blocking oligonucleotide that inhibits amplification of the otter’s DNA when preparing libraries for the vertebrate-specific 12S rRNA gene sequencing assay. The performance of the designed blocking oligonucleotide was assessed in this study.
Materials and methods
Study area and sample collection
A total of seven spraints (S1 Fig) were collected on June 6, 2017 in and with permission from the Ansan Reed Marshy Park (37°16'22.6"N 126°50'24.7"E) in an estuary area in Ansan-si in Gyeonggi-do in South Korea (Fig 1). The park was situated along a branch stream flowing into Sihwaho Lake, a regulating brackish lake separated from the Yellow Sea by a seawall. About 10 g of each spraint sample was collected using a sterile wooden spatula into a 50 ml tube. The collected samples were kept on ice packs and transported to the laboratory on the same day of sample collection. The samples were kept at –80°C until DNA extraction.
DNA extraction and sample confirmation
DNA was extracted from spraints using the PowerMax® Soil DNA Isolation Kit (Mobio Laboratory, Inc., Carlsbad, CA, USA) with modifications [19]. About 5 ml of ultra-pure water was added to each sample, and manually homogenized in the 50 ml tube by a sterile wooden spatula. About 0.2 g of each homogenized sample, including visible bone remains and scales, was transferred into a 2 ml tube containing the Mobio’s Power Beads and solution (750 ml) with additional 0.1 mm diameter glass beads (300 mg) and 0.5 mm diameter glass beads (100 mg) [20]. The samples were homogenized for 3 min by a bead beater (BioSpec Products, OK, USA). After homogenization, DNA was extracted and eluted into 50 μl of TE (10 mM Tris-HCl, 1 mM EDTA, pH = 8.0) according to the kit’s protocol. For the 12S and 16S rRNA gene analyses, three different subsamples (each with 0.2 g) of each homogenized sample in the 50 ml tube were transferred into three different 2 ml tubes for DNA extraction and purification, and recombined for elution into a tube by 50 μl of TE. For the COI analysis, only one part was used for DNA extraction and elution into a tube by 50 μl of TE. The extracted DNA concentrations ranged from 16.6 to 38.9 ng μl-1 for the 12S and 16S rRNA gene analyses, and from 2.5 to 10.6 ng μl-1 for the COI analysis (S1 Table). Each collected sample was analyzed for its identity by a Eurasian otter-specific PCR assay with primers LutcytF and LutcytR targeting the partial cytochrome b sequence [7]. All of the samples were confirmed to be of spraints (S2 Fig).
Designing the blocking oligonucleotide
To block amplification of DNA of Eurasian otters by a universal PCR assay targeting the 12S rRNA gene of vertebrates [17], we designed a Eurasian otter-specific blocking oligonucleotide according to the method reported elsewhere [16]. The designed oligonucleotide OBS1 with a 3-carbon spacer at the 3’-end specifically binds to and blocks amplification of the12S rRNA gene region of the Eurasian otter (Table 1). It was designed not to block amplification of 12S rRNA gene of its potential preys (Table 1).
To confirm the blocking efficacy, PCR assays with or without OBS1 targeting the 12S rRNA gene region of vertebrates with primers 12SV5F and 12SV5R [17] were performed using synthesized Lutra lutra’s DNA as a template. The DNA template containing priming sites of 12SV5F, 12SV5R, and OBS1 was synthesized from 499th to 634th nucleotide positions of the mitochondrial genome of Lutra lutra (FJ236015.1) and cloned into the pTOP Blunt V2 vector (Enzynomics, Seoul, Korea) at Macrogen Incorporation (Seoul, Korea). Each PCR reaction mixture (30 μl) contained 0.26 nM of the synthesized DNA template of Lutra lutra and 0.08 μM of each primer in Premix TaqTM (Takara Bio Inc., Shiga, Japan) with or without 0.8 μM of the blocking oligonucleotide OBS1. The thermal condition was at 95°C for 15 min for initial denaturation, followed by 55 cycles of 30 s at 95°C for denaturation and 30 s at 50°C for annealing. The PCR products were electrophoresed and visualized on a 2.0% (w/w) agarose gel (OmniPur Agarose, Merck, NJ, USA) prepared in 1×TAE buffer by staining with SYBR® Green I Nucleic Acid Gel Stain (Invitrogen, MA, USA).
DNA sequencing
PCR was performed with three sets of taxon-specific primers (Table 2), which are 12SV5F and 12SV5R targeting the 12S rRNA gene region of vertebrates [17], 16SMAV‐F and 16SMAV‐R targeting the 16S rRNA gene region of invertebrates [15], and VF2_t1, FishF2_t1, FishR2_t1 and FR1d_t1 targeting the cytochrome c oxidase 1 (COI) gene for fishes [18]. The expected sizes of amplicons, excluding the lengths of primers, are 98 bp for 12S rRNA gene [17], 36 bp for 16S rRNA gene [15], and 631 bp for COI [18]. These primers were attached with the adapter sequences for Illumina MiSeq (Illumina, Inc., CA, USA). Each PCR was performed in a 30 μl mixture containing 2 μl of DNA extract in Premix TaqTM (Takara Bio Inc., Shiga, Japan). For the vertebrate-specific 12S rRNA gene sequencing assay, 0.08 μM of each primer and 0.8 μM of each of blocking oligonucleotides OBS1 and HomoB [15] were added. The additional assay was conducted without OBS1 and HomoB to check for unintended inhibition of amplification of prey DNA by these oligonucleotides. The thermal condition was at 95°C for 15 min for initial denaturation, followed by 55 cycles of 30 s at 95°C for denaturation and 30 s at 50°C for annealing. For the invertebrate-specific assay, 0.2 μM of each primer and 2 μM of a blocking oligonucleotide for mammals (MamMAVB1) [15] were added. The thermal condition was at 95°C for 15 min for initial denaturation, followed by 55 cycles of 30 s at 95°C for denaturation and 30 s at 55°C for annealing. For the COI assay, 0.167 μM of each primer were added. The thermal condition was at 95°C for 15 min for initial denaturation, followed by 35 cycles of 30 s at 95°C, 40 s at 52°C and 1 min at 72°C, followed by a final extension step at 72°C for 10 min. The amplicons were purified and indexed using the Nextera XT Index kit (Illumina). Each of the purified indexed amplicons was normalized, pooled, thermally-denatured, and loaded onto a v3 600 cycle-kit reagent cartridge (Illumina) for 2 × 300 bp paired-end sequencing by Illumina MiSeq. Raw sequence data are available at NCBI under the BioProject accession numbers PRJNA497180, PRJNA548047, PRJNA548048, and PRJNA565647.
DNA sequence processing and analyses
The raw sequence reads were quality trimmed, and the adapter sequences were removed using Trimmomatic v0.33 [21]. The remaining high-quality sequence reads of 12S and 16S rRNA genes were processed in OBITools [22]. The illuminapairedend command was used to assemble the paired-end reads followed by the obigrep command to remove unaligned sequences. As the barcodes had been already removed and sequences were already demultiplexed, we run the ngsfilter command with -:- option in place of the barcode tag sequence. The sample information was added as attributes in sequence headers of each sample, and the resultant sequences were concatenated and dereplicated using the obiuniq command. The reads shorter than 80 bp for 12S rRNA gene and 20 bp for 16S rRNA gene were removed using the obigrep command. The obiclean command was used to remove the erroneous reads. The remaining sequences were taxonomically assigned using the ecotag command against custom-made 12S and 16S rRNA gene sequences reference databases. These databases were built by extracting relevant regions of 12S and 16S rRNA gene from EMBL nucleotide database (release 138, February 21, 2019) using the ecoPCR program [23]. The ecotag command assigned the query sequence to the least common ancestor in case identification was ambiguous among sibling taxa. For the COI gene sequence analyses, the forward and reverse reads were not joined since the amplicon length (631 bp) was too long to be joined by 2×300 bp paired-end sequencing of Illumina Miseq. They were analyzed separately in mothur v1.40.5 [24]. First, identical reads were binned to generate a set of unique sequences, each with an identical length and nucleotide sequence. Next, the chimeric reads were removed by VSEARCH [25]. Finally, the remaining sequences were classified using the method reported by Wang et al. [26] in classify.seq command with a kmer search (kmer length of 8) against MIDORI database for COI gene [27].
Statistical analyses
R version 3.6.0 was used for statistical analyses. The vegan package version 2.5–6 [28] was used to check for possible changes in the sequence structures by the blocking oligonucleotides OBS1 and for the difference in fish memberships detected by sequencing of two different loci of 12S rRNA and COI genes. For the first purpose, the Bray-Curtis dissimilarities of unique sequence structures were characterized and compared across the 12S rRNA gene libraries prepared with vs. without blocking oligonucleotides OBS1 and HomoB after excluding the sequence reads assigned to Lutra lutra and Homo sapiens. For the second purpose, the Jaccard distances were calculated to characterize the differences of memberships of fish genera detected by 12S rRNA vs. COI gene sequencing analyses. Principal coordinate analysis plot and cluster dendrogram were created based on the Bray-Curtis and Jaccard distances calculated for the first and second purposes, respectively. Permutational multivariate analysis of variance (PERMANOVA) tests were conducted and followed by Kruskal-Wallis rank sum tests with post hoc Wilcoxon rank-sum tests with Bonferroni correction to compare pair-wise intra- and inter-sample distances to assess the degrees of variabilities across the samples and across the library preparation methods. P-values less than 0.05 were considered statistically significant.
Results
Sequencing statistics
The sample O4 was not amplified for 16S rRNA gene sequencing, resulting in n = 6 for the invertebrate analysis. After excluding the reads shorter than 80 bp, 727,245 and 417,605 joined paired-end reads were obtained for the vertebrate 12S rRNA gene sequencing with and without blocking oligonucleotides OBS1 and HomoB, respectively (S2 Table and S3 Table). After excluding the reads shorter than 20 bp, a total of 355,759 sequence reads were obtained by the invertebrate 16S rRNA gene sequencing (S4 Table). In total, 53,468 forward and 53,423 reverse reads of COI gene sequences were obtained (S5 Table). The rarefaction curves appear to reach asymptotes for most libraries, except for those with the small number of sequence reads (e.g., O1 for 16S rRNA gene) (S3 Fig).
Performance of blocking oligonucleotide OBS1
The blocking oligonucleotide OBS1 statistically significantly inhibited the amplification of the predator DNA (p<0.05; Wilcoxon signed rank test), as we found only up to 0.21% of the 12S rRNA gene sequence reads assigned to Lutra lutra by the assay with OBS1 while 10–80% of the reads assigned without OBS1 (Table 3). We confirmed the absence of amplification of the otter’s DNA by the conventional vertebrate-specific PCR assay with OBS1 (Fig 2). Additionally, we confirmed that the variabilities in unique sequence structures were smaller across the library preparation methods (i.e., with vs. without OBS1 and HomoB) than across the samples (Fig 3), indicating little inhibitory effects on prey DNA by OBS1 and HomoB. We found few human DNA sequences in our dataset (Table 3).
Vertebrate 12S rRNA gene
The vertebrate families identified in spraints are shown in Fig 4A and 4B. Based on the assay with OBS1 and HomoB, the most abundant sequences were from the family of Cyprinidae (55%), followed by Channidae (17%), Centrarchidae (9.5%), Mugilidae (7.6%), Ranidae (3.3%), Gobiidae (2.7%), Anguillidae (1.8%), Acheilognathidae (1.0%), and Osphronemidae (0.8%). Except the family Ranidae (amphibian), all other families belong to fishes. At the species or genus levels, 11 taxa were identified with more than 1% of relative abundances from at least one of the analyzed libraries (Fig 5A and 5B). Based on the assay with OBS1 and HomoB, the detected taxa and their mean relative abundances were Carassius auratus (crucian carp) (47%), Channa argus (northern snakehead) (17%), Lepomis macrochirus (bluegill) (8.5%), Mugil cephalus (flathead grey mullet) (7.2%), Anguilla japonica (glass eel) (1.8%), Cyprinus capio (common carp) (1.5%), Acheilognathus barbatus (bitterling fish) (1.0%), Macropodus spp. (paradise fishes) (0.8%), and Micropterus salmoides (largemouth bass) (0.8%). The non-fish taxa detected were Pelophylax nigromaculatus (water frog) (2.7%) and Rana spp. (brown frogs) (0.5%). The vertebrate taxa identified at the species level and their sibling species included in the reference database are listed in S6 Table.
Invertebrate 16S rRNA gene
The invertebrate families detected from spraints are shown in Fig 6. Varunidae was the most abundant (46%), followed by Palaemonidae (18%), Reduviidae (11%), Chironomidae (3.4%), Plumatellidae (2.6%), Atyidae (1.2%), Cicadellidae (1.2%), and Gyrinidae (1.0%). At the species or genus levels, Helicana spp. (mud flat crabs) was the most abundant (46%) and detected from 5 out of 6 spraints (Fig 7). Other species detected included shrimps such as Palaemon paucidens (14%) and Palaemon modestus (3.9%), Paratanytarsus grimmii (chironomid) (3.4%), Macrogyrus oblongus (whirligig beetle) (1.0%), and the ostracod species such as Physocypria sp. X IK-2017 (0.7%) and Cypridopsis vidua (0.2%). The invertebrate taxa identified at the species level and their sibling species included in the reference database are listed in S7 Table.
Fish COI gene
The large numbers of non-fish sequence reads were detected, i.e., 45,613 forward reads (85%) and 44,817 reverse reads (84%) (S5 Table), due to lack of selectivity of the primers for fish detection from environmental samples [18]. Nonetheless, the remaining 7,855 forward reads (15%) and 8,606 reverse reads (16%) were assigned to Actinopteri, a clade encompassing most species of fishes. The identified fish families were Cyprinidae, Channidae, Osphronemidae, Mugilidae, Centrarchidae, and Gobiidae (Fig 4C and 4D), which included Carassius spp. including Carassius auratus and Carassius cuvieri (white crucian carp), Channa spp. including Channa argus, Cyprinus spp., Macropodus chinensis, and Mugil spp. (Fig 5C and 5D). The fish taxa identified at the species level and their sibling species included in the database are listed in S8 Table.
The fish genera detected by the COI sequencing appear to be similar to those by the 12S rRNA gene sequencing (Fig 8), with examples of the identical genera detected in the samples O2 and O5 between the two different sequencing loci (S9 Table). Fig 8 shows that the variabilities in the memberships of detected fish genera were smaller across the sequencing loci (i.e., 12S rRNA gene vs. COI) than across the samples, with a statistical significance observed between the intra-sample distances based on the two different loci and the inter-sample distances based on 12S rRNA gene sequencing, indicating that the detected fish genera are in agreement between the two locus analyses in terms of detected fish memberships.
Discussion
In this study, we designed a Eurasian otter-specific blocking oligonucleotide (OBS1) for their diet analyses based on 12S rRNA gene sequencing for vertebrates. The use of blocking oligonucleotides is an effective solution to suppress amplification of DNA fragments derived from the predators such as brown bears [15], leopard cats [11], penguins [29], and seals [30]. It has been shown that blocking oligonucleotides can increase sensitivity of detecting rare preys by suppressing amplification of predator DNA [11]. The blocking oligonucleotide OBS1 designed to specifically bind to the 12S rRNA gene region of the Eurasian otter (Table 1) successfully inhibited amplification of the synthesized DNA fragments of the Eurasian otter (Fig 2), with only up to 0.21% of the total reads assigned to the otter by the vertebrate-specific 12S rRNA gene sequencing (Table 3). We confirmed that the types of fishes detected by the 12S rRNA gene sequencing were congruent with those detected by the COI gene sequencing in terms of their memberships (Fig 8 and S9 Table). Moreover, we confirmed that the variabilities in unique sequence structures were smaller across the library preparation methods (i.e., with vs. without OBS1) than across the samples (Fig 3), suggesting little inhibitory effects on prey DNA by the designed blocking oligonucleotide OBS1.
Our analyses of 12S rRNA gene sequences of vertebrates showed that fishes (95%) were the primary diet, followed by amphibians (3.3%) for the Eurasian otters inhabiting a marshy estuary area of South Korea in the early summer (Fig 4A). Previous studies [31–33] demonstrated that 12S rRNA gene sequences provide accurate identification of most of fishes including teleosts at species level, highlighting the 12S rRNA gene as an ideal marker for diet analyses of Eurasian otters. Additionally, we confirmed that fish species identified in this study are in agreements with fish species detected from spraints by a previous study [10] and spatial distributions of their habitats in Korea [34].
Our observation of fishes as dominant preys of the Eurasian otters in South Korea is in congruence with previous studies showing dominance of fish preys of the Eurasian otters in Europe based on conventional morphological identification methods [3, 4]. Similar to the previous findings [35–37], the majority of the detected fishes belonged to the family Cyprinidae (Fig 4) among which Carassius auratus (crucian carp) was the most dominant species detected across all spraints (Fig 5). A previous study also reported that Cyprinidae was the most dominant diet by the Eurasian otters along a river basin in South Korea [10]. C. auratus dominates fresh water ecosystems in South Korea [38], and therefore these might represent an abundant prey in their local habitats. The other freshwater species detected in this study included Channa argus (northern snakehead), Lepomis macrochirus (bluegill) and Cyprinus capio (common carp). These freshwater species are known to inhabit Korea [39].
Euryhaline species, such as Mugil cephalus (flathead grey mullet), were also detected, and they are known to inhabit freshwater and seawater ecosystems in Korea [40]. The detection of these species seems reasonable since the spraints were collected in an estuary area with a freshwater stream flowing into a regulating brackish lake neighboring the Yellow Sea (Fig 1). M. cephalus possibly migrates into the freshwater stream in the area where the Eurasian otters inhabit, and/or the Eurasian otters possibly feed in the areas of the brackish lake where M. cephalus inhabits.
Introduced species such as Lepomis macrochirus (bluegill) and Micropterus salmoides (largemouth bass) were detected (Fig 5A and 5B), with their respective mean relative abundances of 8.5% and 0.8%. L. macrochirus is native to North America, and was first imported to Japan in 1960 [41]. Later, a group of L. macrochirus was imported from Japan and released into a reservoir in Korea from which it was dispersed and became a dominant species in freshwater systems in Korea [41]. The Eurasian otters fed on these exotic species, but their quantities, i.e., 8.5% for L. macrochirus and 0.8% for M. salmoides, were smaller than those of the native species such as Carassius auratus (crucian carp) (47%) and Channa argus (northern snakehead) (17%), which was similar to the previous finding that the exotic fishes, including Lepomis gibbosus (pumpkinseed), were preyed at disproportionately lower levels than the native fishes, including Anguilla anguilla (European eel), by the Eurasian otters in England [42].
Non-fish vertebrates were also detected, which included frogs such as Pelophylax nigromaculatus (water frog) and Rana spp. (brown frogs) (Fig 5A and 5B). Eurasian otters prey on amphibians when fishes become scarce [37]. Both of these frog species inhabit Korea [43], and a study reported Rana spp. in spraints collected along a river basin in Korea [10]. However, they were not abundant, with mean relative abundances of 2.7% for Pelophylax nigromaculatus and 0.5% for Rana spp., suggesting that frogs are not primary diets of the Eurasian otters during this time of the year at our sampling site possibly because of the sufficient availability of fishes.
Among invertebrates, Helicana spp. (mud flat crabs) were the most abundant (46%) (Fig 7), followed by shrimps such as Palaemon paucidens and Palaemon modestus. Helicana spp. are widely distributed in coastal wetlands in Korea [44, 45]. Similarly, P. modestus and P. paucidens are endemic in freshwater ecosystems of Korea [46]. Eurasian otters are known to feed on crabs and other crustaceans [37, 47, 48]. Additionally, the insects were detected from some spraints with relatively high abundances. Insects are known to be fed by Eurasian otters [3, 37]. We detected Macrogyrus oblongus (whirligig beetle) with 6.2% of relative abundance from the sample O2. This water beetle is known to inhabit parts of Oceania and South America, but not East Asia [49]. This suggests potential misidentification of a closely related genus such as Dineutus that is known to inhabit Korea [49]. The family Reduviidae (assassin bugs) was also detected (Fig 6), with a high relative abundance (65%) in the sample O5, but the sequences were not identified down to the genus level. Reduviidae is unlikely a common prey for the Eurasian otters since no such observation has been reported by previous research. The otter might opportunistically prey Reduviidae that might be present in its local habitat.
There are several limitations in our DNA metabarcoding-based diet analyses. First, caution must be taken when drawing conclusions of relative biomass consumed based on relative sequence abundance in environmental DNA since the relationship between sequence abundance and biomass can often vary and be weakly correlated due to differences in amplification efficiency of primers particularly in COI gene [50–52]. A semiquantitative estimate of diet composition could be obtained in future studies after applying the robust filtering procedure described by Corse et al. [53]. It is also impossible to compare relative importance between invertebrates and vertebrates based on our multiple locus approach. Second, we detected millimeter-size organisms such as Paratanytarsus grimmii (chironomid) (3.4%), and the species of ostracods Physocypria sp. (0.7%) and Cypridopsis vidua (0.2%) (Fig 7). Due to their small sizes, these organisms were unlikely predated by the Eurasian otters. Instead, these organisms were likely incidentally ingested by and/or attached to larger organisms predated by the otters. The secondary predation is also possible, for instance, that the Eurasian otters feed on piscivorous fishes (e.g., northern snakehead) that previously fed on smaller fishes (e.g., crucial carp) and contained them in their guts. We do not know such relationships. Other potential limitations include the inability of detecting possible cannibalism [54] because, unless special approaches, such as the so-called the variant-centered approach for haplotype detection [53, 55], are employed, DNA metabarcoding cannot distinguish DNA fragments from cannibalized and cannibalizing individuals. Finally, a limitation of this study was small sample size as the accumulation curves of the numbers of detected prey genera do not appear to reach to the asymptotes (S4 Fig). Future research should include more samples to capture the entire picture of the diet profiles in this study site.
Conclusion
In this study, we applied HTS-based DNA metabarcoding analyses targeting vertebrate 12S rRNA gene, invertebrate 16S rRNA gene, and fish COI gene sequences to characterize diet profiles of the Eurasian otters inhabiting a marshy estuary area of South Korea. In particular, we developed a Eurasian otter-specific blocking oligonucleotide (OBS1) for 12S rRNA gene sequencing for vertebrates. Overall, our study demonstrated that HTS-based DNA metabarcoding analyses were useful to provide in-depth information regarding their diet compositions. Our multi locus approach showed that the combination of 12S and 16S rRNA gene sequencing could cover a broad range of invertebrate and vertebrate prey animals of the Eurasian otters, including fishes detected by the COI gene sequencing. Our results suggest that the 12S and 16S rRNA gene sequencing analyses suffice for diet analyses of the Eurasian otters, with optional COI gene sequencing for fish-specific detection and identification. In this study, spraint samples were collected at one location in one season, and the diet profiles might differ spatiotemporally [3, 4]. By using HTS-based metabarcoding analyses, future research will explore detailed taxonomies of the otters’ diets across locations and seasons.
Supporting information
S1 Fig [tif]
Spraints collected in this study.
S2 Fig [tif]
DNA bands of the PCR amplicons by the Eurasian otter-specific primers LutcytF and LutcytR.
S3 Fig [a]
Rarefaction curves for observed unique sequences.
S4 Fig [tif]
Cumulative distributions of the numbers of detected genera.
S1 Table [xlsx]
DNA concentrations of extracts from spraint samples.
S2 Table [xlsx]
Number of high-quality vertebrate 12S rRNA gene sequence reads prepared with blocking oligonucleotides OBS1 and HomoB.
S3 Table [xlsx]
Number of high-quality vertebrate 12S rRNA gene sequence reads prepared without blocking oligonucleotides OBS1 and HomoB.
S4 Table [xlsx]
Number of high-quality invertebrate 16S rRNA gene sequence reads.
S5 Table [xlsx]
Number of high-quality COI gene sequence reads.
S6 Table [xlsx]
List of species in the reference database of 12S rRNA gene sequences.
S7 Table [xlsx]
List of species in the reference database of 16S rRNA gene sequences.
S8 Table [xlsx]
List of species in the reference database of COI sequences.
S9 Table [xlsx]
Comparison of fish genera detected by 12S rRNA gene and COI sequencing.
Zdroje
1. Roos A, Loy A, de Silva P, Hajkova P, Zemanová B. Lutra lutra. The IUCN Red List of Threatened Species 2015. 2015:e.T12419A21935287.
2. Won C, Smith KG. History and current status of mammals of the Korean Peninsula. Mammal Rev. 1999;29(1):3–36.
3. Krawczyk AJ, Bogdziewicz M, Majkowska K, Glazaczow A. Diet composition of the Eurasian otter Lutra lutra in different freshwater habitats of temperate Europe: a review and meta‐analysis. Mammal Rev. 2016;46(2):106–113.
4. Lanszki J, Lehoczky I, Kotze A, Somers MJ. Diet of otters (Lutra lutra) in various habitat types in the Pannonian biogeographical region compared to other regions of Europe. PeerJ. 2016;4:e2266. doi: 10.7717/peerj.2266 27602262
5. Mateo R, Millán J, Rodríguez-Estival J, Camarero PR, Palomares F, Ortiz-Santaliestra ME. Levels of organochlorine pesticides and polychlorinated biphenyls in the critically endangered Iberian lynx and other sympatric carnivores in Spain. Chemosphere. 2012;86(7):691–700. doi: 10.1016/j.chemosphere.2011.10.037 22099537
6. Guillaud E, Bearez P, Denys C, Raimond S. New data on fish diet and bone digestion of the Eurasian otter (Lutra lutra) (Mammalia: Mustelidae) in central France. Eur Zool J. 2017;84(1):226–237.
7. Park H-C, Han T-Y, Kim D-C, Min M-S, Han S-Y, Kim K-S, et al. Individual identification and sex determination of Eurasian otters (Lutra lutra) in Daegu city based on genetic analysis of otter spraint. Genes Genom. 2011;33(6):653–657.
8. Park C-S, Cho G-J. Individual identification of Eurasian otters (Lutra lutra) in South Korea (Sincheon River, Daegu) by microsatellite markers. J Vet Med Sci. 2017;79(6):1064–1067. doi: 10.1292/jvms.16-0563 28428484
9. Carss D. Foraging behaviour and feeding ecology of the otter Lutra lutra: a selective review. Hystrix It J Mamm. 1995;7(1–2):179–194.
10. Hong S, Gim J-S, Kim HG, Cowan P, Joo G-J. A molecular approach to identifying the relationship between resource use and availability in Eurasian otters (Lutra lutra). Can J Zool. 2019;97:797–804.
11. Shehzad W, Riaz T, Nawaz MA, Miquel C, Poillot C, Shah SA, et al. Carnivore diet analysis based on next-generation sequencing: Application to the leopard cat (Prionailurus bengalensis) in Pakistan. Mol Ecol. 2012;21(8):1951–1965. doi: 10.1111/j.1365-294X.2011.05424.x 22250784
12. Farrell LE, Roman J, Sunquist ME. Dietary separation of sympatric carnivores identified by molecular analysis of scats. Mol Ecol. 2000;9(10):1583–1590. doi: 10.1046/j.1365-294x.2000.01037.x 11050553
13. Palomares F, Godoy JA, Piriz A, O'Brien SJ. Faecal genetic analysis to determine the presence and distribution of elusive carnivores: design and feasibility for the Iberian lynx. Mol Ecol. 2002;11(10):2171–2182. doi: 10.1046/j.1365-294x.2002.01608.x 12296958
14. Vesterinen EJ, Ruokolainen L, Wahlberg N, Peña C, Roslin T, Laine VN, et al. What you need is what you eat? Prey selection by the bat Myotis daubentonii. Mol Ecol. 2016;25(7):1581–1594. doi: 10.1111/mec.13564 26841188
15. De Barba M, Miquel C, Boyer F, Mercier C, Rioux D, Coissac E, et al. DNA metabarcoding multiplexing and validation of data accuracy for diet assessment: application to omnivorous diet. Mol Ecol Resour. 2014;14(2):306–323. doi: 10.1111/1755-0998.12188 24128180
16. Vestheim H, Jarman SN. Blocking primers to enhance PCR amplification of rare sequences in mixed samples–a case study on prey DNA in Antarctic krill stomachs. Front Zool. 2008;5(1):12.
17. Riaz T, Shehzad W, Viari A, Pompanon F, Taberlet P, Coissac E. ecoPrimers: inference of new DNA barcode markers from whole genome sequence analysis. Nucleic Acids Res. 2011;39(21):e145–e145. doi: 10.1093/nar/gkr732 21930509
18. Ivanova NV, Zemlak TS, Hanner RH, Hebert PDN. Universal primer cocktails for fish DNA barcoding. Mol Ecol Notes. 2007;7(4):544–548.
19. An C, Okamoto Y, Xu S, Eo KY, Kimura J, Yamamoto N. Comparison of fecal microbiota of three captive carnivore species inhabiting Korea. J Vet Med Sci. 2017;79(3):542–546. doi: 10.1292/jvms.16-0472 28049922
20. Hospodsky D, Yamamoto N, Peccia J. Accuracy, precision, and method detection limits of quantitative PCR for airborne bacteria and fungi. Appl Environ Microbiol. 2010;76(21):7004–7012. doi: 10.1128/AEM.01240-10 20817798
21. Bolger AM, Lohse M, Usadel B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30(15):2114–2120. doi: 10.1093/bioinformatics/btu170 24695404
22. Boyer F, Mercier C, Bonin A, Le Bras Y, Taberlet P, Coissac E. obitools: a unix-inspired software package for DNA metabarcoding. Mol Ecol Resour. 2016;16(1):176–182. doi: 10.1111/1755-0998.12428 25959493
23. Ficetola GF, Coissac E, Zundel S, Riaz T, Shehzad W, Bessière J, et al. An in silico approach for the evaluation of DNA barcodes. BMC Genomics. 2010;11(1):434.
24. Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–7541. doi: 10.1128/AEM.01541-09 19801464
25. Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4:e2584. doi: 10.7717/peerj.2584 27781170
26. Wang Q, Garrity GM, Tiedje JM, Cole JR. Naïve Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73(16):5261–5267. doi: 10.1128/AEM.00062-07 17586664
27. Machida RJ, Leray M, Ho S-L, Knowlton N. Metazoan mitochondrial gene sequence reference datasets for taxonomic assignment of environmental samples. Sci Data. 2017;4:170027. doi: 10.1038/sdata.2017.27 28291235
28. Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: Community Ecology Package. R package version 2.5–6 [Internet]. 2019 [cited 2019 Sep 9]. Available from: https://CRAN.R-project.org/package=vegan.
29. Deagle BE, Chiaradia A, McInnes J, Jarman SN. Pyrosequencing faecal DNA to determine diet of little penguins: is what goes in what comes out? Conserv Genet. 2010;11(5):2039–2048.
30. Deagle BE, Kirkwood R, Jarman SN. Analysis of Australian fur seal diet by pyrosequencing prey DNA in faeces. Mol Ecol. 2009;18(9):2022–2038. doi: 10.1111/j.1365-294X.2009.04158.x 19317847
31. Miya M, Sato Y, Fukunaga T, Sado T, Poulsen JY, Sato K, et al. MiFish, a set of universal PCR primers for metabarcoding environmental DNA from fishes: detection of more than 230 subtropical marine species. Royal Soc Open Sci. 2(7):150088.
32. Stoeckle MY, Soboleva L, Charlop-Powers Z. Aquatic environmental DNA detects seasonal fish abundance and habitat preference in an urban estuary. PLOS ONE. 2017;12(4):e0175186. doi: 10.1371/journal.pone.0175186 28403183
33. Valentini A, Taberlet P, Miaud C, Civade R, Herder J, Thomsen PF, et al. Next-generation monitoring of aquatic biodiversity using environmental DNA metabarcoding. Mol Ecol. 2016;25(4):929–942. doi: 10.1111/mec.13428 26479867
34. Yoon J-D, Kim J-H, Byeon M-S, Yang H-J, Park J-Y, Shim J-H, et al. Distribution patterns of fish communities with respect to environmental gradients in Korean streams. Ann Limnol—Int J Lim. 2011;47:S63–S71.
35. Taastrøm H-M, Jacobsen L. The diet of otters (Lutra lutra L.) in Danish freshwater habitats: comparisons of prey fish populations. J Zool. 1999;248(1):1–13.
36. Copp GH, Roche K. Range and diet of Eurasian otters Lutra lutra (L.) in the catchment of the River Lee (south‐east England) since re‐introduction. Aquat Conserv: Mar Freshw Ecosyst. 2003;13(1):65–76.
37. Remonti L, Prigioni C, Balestrieri A, Sgrosso S, Priore G. Trophic flexibility of the otter (Lutra lutra) in southern Italy. Mamm Biol. 2008;73(4):293–302.
38. Jang MH, Kim JG, Park SB, Jeong KS, Cho GI, Joo GJ. The current status of the distribution of introduced fish in large river systems of South Korea. Int Rev Hydrobiol. 2002;87(2–3):319–328.
39. Jang M-H, Joo G-J, Lucas MC. Diet of introduced largemouth bass in Korean rivers and potential interactions with native fishes. Ecol Freshwat Fish. 2006;15(3):315–320.
40. Park SK, Davidson K, Pan M. Economic relationships between aquaculture and capture fisheries in the Republic of Korea. Aquacult Econ Manag. 2012;16(2):102–116.
41. Kawamura K, Yonekura R, Katano O, Taniguchi Y, Saitoh K. Origin and dispersal of bluegill sunfish, Lepomis macrochirus, in Japan and Korea. Mol Ecol. 2006;15(3):613–621. doi: 10.1111/j.1365-294X.2006.02823.x 16499689
42. Miranda R, Copp GH, Williams J, Beyer K, Gozlan RE. Do Eurasian otters Lutra lutra (L.) in the Somerset Levels prey preferentially on non-native fish species? Fundam Appl Limnol. 2008;172(4):339–347.
43. Ju J, Kim JH, Kim SH. Habitat fragmentation by a levee and its Impact on frog population in the civilian control zone. J Wetlands Res. 2016;18(2):113–120.
44. Park J, Song SJ, Ryu J, Kwon B-O, Hong S, Bae H, et al. Macrozoobenthos of Korean tidal flats: A review on species assemblages and distribution. Ocean Coast Manag. 2014;102:483–492.
45. Shih H-T, Suzuki H. Taxonomy, phylogeny, and biogeography of the endemic mudflat crab Helicel Chasmagnathus complex (Crustacea: Brachyura: Varunidae) from East Asia. Zool Stud. 2008;47(1):114–125.
46. Oh C-W, Suh H-L, Park K-Y, Ma C-W, Lim H-S. Growth and reproductive biology of the freshwater shrimp Exopalaemon modestus (Decapoda: Palaemonidae) in a lake of Korea. J Crustacean Biol. 2002;22(2):357–366.
47. Georgiev D, Stoycheva S. Freshwater crabs preyed on by the Eurasian Otter Lutra lutra in a river habitat of southern Bulgaria. Hystrix It J Mamm. 2007;17(2):129–135.
48. Freitas D, Gomes J, Luis TS, Madruga L, Marques C, Baptista G, et al. Otters and fish farms in the Sado estuary: ecological and socio-economic basis of a conflict. Hydrobiologia. 2007;587(1):51–62.
49. Gustafson GT, Miller KB. Systematics and evolution of the whirligig beetle tribe Dineutini (Coleoptera: Gyrinidae: Gyrininae). Zool J Linnean Soc. 2017;181(1):118–150.
50. Bowles E, Schulte PM, Tollit DJ, Deagle BE, Trites AW. Proportion of prey consumed can be determined from faecal DNA using real-time PCR. Mol Ecol Resour. 2011;11(3):530–540. doi: 10.1111/j.1755-0998.2010.02974.x 21481211
51. Elbrecht V, Leese F. Can DNA-based ecosystem assessments quantify species abundance? Testing primer bias and biomass—sequence relationships with an innovative metabarcoding protocol. PLOS ONE. 2015;10(7):e0130324. doi: 10.1371/journal.pone.0130324 26154168
52. Thomas AC, Deagle BE, Eveson JP, Harsch CH, Trites AW. Quantitative DNA metabarcoding: improved estimates of species proportional biomass using correction factors derived from control material. Mol Ecol Resour. 2016;16(3):714–726. doi: 10.1111/1755-0998.12490 26602877
53. Corse E, Meglécz E, Archambaud G, Ardisson M, Martin J-F, Tougard C, et al. A from-benchtop-to-desktop workflow for validating HTS data and for taxonomic identification in diet metabarcoding studies. Mol Ecol Resour. 2017;17(6):e146–e159. doi: 10.1111/1755-0998.12703 28776936
54. Simpson VR, Coxon KE. Intraspecific aggression, cannibalism and suspected infanticide in Otters. Br Wildl. 2000;11(6):423–426.
55. Corse E, Tougard C, Archambaud-Suard G, Agnèse J-F, Messu Mandeng FD, Bilong Bilong CF, et al. One-locus-several-primers: A strategy to improve the taxonomic and haplotypic coverage in diet metabarcoding studies. Ecol Evol. 2019;9(8):4603–4620. doi: 10.1002/ece3.5063 31031930
Článek vyšel v časopise
PLOS One
2019 Číslo 12
- Řekněte mi, co jíte, a já vám řeknu, jaká je úroveň vaší kognice
- Původní český test ALBA pro odhalení mírných kognitivních deficitů získal certifikaci MZ ČR
- Duševní onemocnění se u teenagerů přenáší v rámci školní třídy
- Od kostiček k lidské kůži. Jak si postavit 3D biotiskárnu?
- Prof. Libor Grubhoffer: Nová varianta mpox je příležitostí pokročit v poznání původce varioly
Nejčtenější v tomto čísle
- Methylsulfonylmethane increases osteogenesis and regulates the mineralization of the matrix by transglutaminase 2 in SHED cells
- Oregano powder reduces Streptococcus and increases SCFA concentration in a mixed bacterial culture assay
- The characteristic of patulous eustachian tube patients diagnosed by the JOS diagnostic criteria
- Parametric CAD modeling for open source scientific hardware: Comparing OpenSCAD and FreeCAD Python scripts