Environmental DNA metabarcoding as a useful tool for evaluating terrestrial mammal diversity in tropical forests.

Innovative techniques, such as environmental DNA (eDNA) metabarcoding, are now promoting broader biodiversity monitoring at unprecedented scales, because of the reduction in time, presumably lower cost, and methodological efficiency. Our goal was to assess the efficiency of established inventory techniques (live-trapping grids, pitfall traps, camera trapping, mist netting) as well as eDNA for detecting Amazonian mammals. For terrestrial small mammals, we used 32 live-trapping grids based on Sherman and Tomahawk traps (total effort of 10,368 trap-nights); in addition to 16 pitfall traps (1,408 trap-nights). For bats, we used mist nets at 8 sites (4,800 net hours). For medium and large mammals, we used 72 camera trap stations (5,208 camera-days). We identified vertebrate and mammal taxa based on eDNA analysis (12S region, with V05 and Mamm01 markers) from water samples, including a total of 11 3-km transects for stagnant water sampling and seven small streams for running water sampling. A total of 106 mammal species were recorded. Building on sample-based rarefaction and extrapolation curves, both trapping grids and pitfall were successful, recording 91.16% and 82.1% of the expected species for these techniques (~22 and ~9 species), and 16.98% and 6.60% of the total recorded mammal species, respectively. Mist nets recorded 83.2% of the expected bat species (~48), and 34.91% of the total recorded species. Camera trapping recorded 99.2% of the predicted large- and medium-sized species (~31), and 33.02% of the total recorded species. eDNA recorded 75.4% of the expected mammal species for this technique (~68), and 47.0% of the total recorded species. eDNA resulted in a useful tool that saves on effort and reduces sampling costs. This study is among the first to show the large potential of eDNA metabarcoding for assessing Amazonian mammal communities, providing, in combination with conventional techniques, a rapid overview of mammal diversity with broad applications to monitoring, management and conservation. By including appropriate genetic markers and updated reference databases, eDNA metabarcoding method can be extended to the whole vertebrate community.

Evolutionary history of enigmatic bears in the Tibetan Plateau-Himalaya region and the identity of the yeti.

Although anecdotally associated with local bears (Ursus arctos and U. thibetanus), the exact identity of 'hominid'-like creatures important to folklore and mythology in the Tibetan Plateau-Himalaya region is still surrounded by mystery. Recently, two purported yeti samples from the Himalayas showed genetic affinity with an ancient polar bear, suggesting they may be from previously unrecognized, possibly hybrid, bear species, but this preliminary finding has been under question. We conducted a comprehensive genetic survey of field-collected and museum specimens to explore their identity and ultimately infer the evolutionary history of bears in the region. Phylogenetic analyses of mitochondrial DNA sequences determined clade affinities of the purported yeti samples in this study, strongly supporting the biological basis of the yeti legend to be local, extant bears. Complete mitochondrial genomes were assembled for Himalayan brown bear (U. a. isabellinus) and black bear (U. t. laniger) for the first time. Our results demonstrate that the Himalayan brown bear is one of the first-branching clades within the brown bear lineage, while Tibetan brown bears diverged much later. The estimated times of divergence of the Tibetan Plateau and Himalayan bear lineages overlap with Middle to Late Pleistocene glaciation events, suggesting that extant bears in the region are likely descendants of populations that survived in local refugia during the Pleistocene glaciations.

Next-generation monitoring of aquatic biodiversity using environmental DNA metabarcoding.

Global biodiversity in freshwater and the oceans is declining at high rates. Reliable tools for assessing and monitoring aquatic biodiversity, especially for rare and secretive species, are important for efficient and timely management. Recent advances in DNA sequencing have provided a new tool for species detection from DNA present in the environment. In this study, we tested whether an environmental DNA (eDNA) metabarcoding approach, using water samples, can be used for addressing significant questions in ecology and conservation. Two key aquatic vertebrate groups were targeted: amphibians and bony fish. The reliability of this method was cautiously validated in silico, in vitro and in situ. When compared with traditional surveys or historical data, eDNA metabarcoding showed a much better detection probability overall. For amphibians, the detection probability with eDNA metabarcoding was 0.97 (CI = 0.90-0.99) vs. 0.58 (CI = 0.50-0.63) for traditional surveys. For fish, in 89% of the studied sites, the number of taxa detected using the eDNA metabarcoding approach was higher or identical to the number detected using traditional methods. We argue that the proposed DNA-based approach has the potential to become the next-generation tool for ecological studies and standardized biodiversity monitoring in a wide range of aquatic ecosystems.

Highly overlapping winter diet in two sympatric lemming species revealed by DNA metabarcoding.

Sympatric species are expected to minimize competition by partitioning resources, especially when these are limited. Herbivores inhabiting the High Arctic in winter are a prime example of a situation where food availability is anticipated to be low, and thus reduced diet overlap is expected. We present here the first assessment of diet overlap of high arctic lemmings during winter based on DNA metabarcoding of feces. In contrast to previous analyses based on microhistology, we found that the diets of both collared (Dicrostonyx groenlandicus) and brown lemmings (Lemmus trimucronatus) on Bylot Island were dominated by Salix while mosses, which were significantly consumed only by the brown lemming, were a relatively minor food item. The most abundant plant taxon, Cassiope tetragona, which alone composes more than 50% of the available plant biomass, was not detected in feces and can thus be considered to be non-food. Most plant taxa that were identified as food items were consumed in proportion to their availability and none were clearly selected for. The resulting high diet overlap, together with a lack of habitat segregation, indicates a high potential for resource competition between the two lemming species. However, Salix is abundant in the winter habitats of lemmings on Bylot Island and the non-Salix portion of the diets differed between the two species. Also, lemming grazing impact on vegetation during winter in the study area is negligible. Hence, it seems likely that the high potential for resource competition predicted between these two species did not translate into actual competition. This illustrates that even in environments with low primary productivity food resources do not necessarily generate strong competition among herbivores.

No evidence for the effect of MHC on male mating success in the brown bear.

Mate choice is thought to contribute to the maintenance of the spectacularly high polymorphism of the Major Histocompatibility Complex (MHC) genes, along with balancing selection from parasites, but the relative contribution of the former mechanism is debated. Here, we investigated the association between male MHC genotype and mating success in the brown bear. We analysed fragments of sequences coding for the peptide-binding region of the highly polymorphic MHC class I and class II DRB genes, while controlling for genome-wide effects using a panel of 18 microsatellite markers. Male mating success did not depend on the number of alleles shared with the female or amino-acid distance between potential mates at either locus. Furthermore, we found no indication of female mating preferences for MHC similarity being contingent on the number of alleles the females carried. Finally, we found no significant association between the number of MHC alleles a male carried and his mating success. Thus, our results provided no support for the role of mate choice in shaping MHC polymorphism in the brown bear.

Fifty thousand years of Arctic vegetation and megafaunal diet.

Although it is generally agreed that the Arctic flora is among the youngest and least diverse on Earth, the processes that shaped it are poorly understood. Here we present 50 thousand years (kyr) of Arctic vegetation history, derived from the first large-scale ancient DNA metabarcoding study of circumpolar plant diversity. For this interval we also explore nematode diversity as a proxy for modelling vegetation cover and soil quality, and diets of herbivorous megafaunal mammals, many of which became extinct around 10 kyr bp (before present). For much of the period investigated, Arctic vegetation consisted of dry steppe-tundra dominated by forbs (non-graminoid herbaceous vascular plants). During the Last Glacial Maximum (25-15 kyr bp), diversity declined markedly, although forbs remained dominant. Much changed after 10 kyr bp, with the appearance of moist tundra dominated by woody plants and graminoids. Our analyses indicate that both graminoids and forbs would have featured in megafaunal diets. As such, our findings question the predominance of a Late Quaternary graminoid-dominated Arctic mammoth steppe.

Fungal palaeodiversity revealed using high-throughput metabarcoding of ancient DNA from arctic permafrost.

The taxonomic and ecological diversity of ancient fungal communities was assessed by combining next generation sequencing and metabarcoding of DNA preserved in permafrost. Twenty-six sediment samples dated 16 000-32 000 radiocarbon years old from two localities in Siberia were analysed for fungal ITS. We detected 75 fungal OTUs from 21 orders representing three phyla, although rarefaction analyses suggested that the full diversity was not recovered despite generating an average of 6677 ± 3811 (mean ± SD) sequences per sample and that preservation bias likely has considerable effect on the recovered DNA. Most OTUs (75.4%) represented ascomycetes. Due to insufficient sequencing depth, DNA degradation and putative preservation biases in our samples, the recovered taxa probably do not represent the complete historic fungal community, and it is difficult to determine whether the fungal communities varied geographically or experienced a composition shift within the period of 16 000-32 000 bp. However, annotation of OTUs to functional ecological groups provided a wealth of information on the historic communities. About one-third of the OTUs are presumed plant-associates (pathogens, saprotrophs and endophytes) typical of graminoid- and forb-rich habitats. We also detected putative insect pathogens, coprophiles and keratinophiles likely associated with ancient insect and herbivore faunas. The detection of putative insect pathogens, mycoparasites, aquatic fungi and endophytes broadens our previous knowledge of the diversity of fungi present in Beringian palaeoecosystems. A large group of putatively psychrophilic/psychrotolerant fungi was also detected, most likely representing a modern, metabolically active fungal community.

Demographic history and genetic diversity in West Indian Coereba flaveola populations.

The bananaquit (Coereba flaveola) has been well studied throughout the Caribbean region from a phylogenetic perspective. However, data concerning the population genetics and long-term demography of this bird species are lacking. In this study, we focused on three populations within the Lesser Antilles and one on Puerto Rico and assessed genetic and demographic processes, using five nuclear and two mitochondrial markers. We found that genetic diversity of bananaquits on Puerto Rico exceeds that on the smaller islands (Dominica, Guadeloupe and Grenada); this might reflect either successive founder events from Puerto Rico to Grenada, or more rapid drift in smaller populations subsequent to colonization. Population growth rate estimates showed no evidence of rapid expansion and migration was indicated only between populations from the closest islands of Dominica and Guadeloupe. Overall, our results suggest that a "demographic fission" model, considering only mutation and drift, but without migration, can be applied to these bananaquit populations in the West Indies.

New environmental metabarcodes for analysing soil DNA: potential for studying past and present ecosystems.

Metabarcoding approaches use total and typically degraded DNA from environmental samples to analyse biotic assemblages and can potentially be carried out for any kinds of organisms in an ecosystem. These analyses rely on specific markers, here called metabarcodes, which should be optimized for taxonomic resolution, minimal bias in amplification of the target organism group and short sequence length. Using bioinformatic tools, we developed metabarcodes for several groups of organisms: fungi, bryophytes, enchytraeids, beetles and birds. The ability of these metabarcodes to amplify the target groups was systematically evaluated by (i) in silico PCRs using all standard sequences in the EMBL public database as templates, (ii) in vitro PCRs of DNA extracts from surface soil samples from a site in Varanger, northern Norway and (iii) in vitro PCRs of DNA extracts from permanently frozen sediment samples of late-Pleistocene age (~16,000-50,000 years bp) from two Siberian sites, Duvanny Yar and Main River. Comparison of the results from the in silico PCR with those obtained in vitro showed that the in silico approach offered a reliable estimate of the suitability of a marker. All target groups were detected in the environmental DNA, but we found large variation in the level of detection among the groups and between modern and ancient samples. Success rates for the Pleistocene samples were highest for fungal DNA, whereas bryophyte, beetle and bird sequences could also be retrieved, but to a much lesser degree. The metabarcoding approach has considerable potential for biodiversity screening of modern samples and also as a palaeoecological tool.

Blocking human contaminant DNA during PCR allows amplification of rare mammal species from sedimentary ancient DNA.

Analyses of degraded DNA are typically hampered by contamination, especially when employing universal primers such as commonly used in environmental DNA studies. In addition to false-positive results, the amplification of contaminant DNA may cause false-negative results because of competition, or bias, during the PCR. In this study, we test the utility of human-specific blocking primers in mammal diversity analyses of ancient permafrost samples from Siberia. Using quantitative PCR (qPCR) on human and mammoth DNA, we first optimized the design and concentration of blocking primer in the PCR. Subsequently, 454 pyrosequencing of ancient permafrost samples amplified with and without the addition of blocking primer revealed that DNA sequences from a diversity of mammalian representatives of the Beringian megafauna were retrieved only when the blocking primer was added to the PCR. Notably, we observe the first retrieval of woolly rhinoceros (Coelodonta antiquitatis) DNA from ancient permafrost cores. In contrast, reactions without blocking primer resulted in complete dominance by human DNA sequences. These results demonstrate that in ancient environmental analyses, the PCR can be biased towards the amplification of contaminant sequences to such an extent that retrieval of the endogenous DNA is severely restricted. The application of blocking primers is a promising tool to avoid this bias and can greatly enhance the quantity and the diversity of the endogenous DNA sequences that are amplified.

ITS as an environmental DNA barcode for fungi: an in silico approach reveals potential PCR biases.

BACKGROUND: During the last 15 years the internal transcribed spacer (ITS) of nuclear DNA has been used as a target for analyzing fungal diversity in environmental samples, and has recently been selected as the standard marker for fungal DNA barcoding. In this study we explored the potential amplification biases that various commonly utilized ITS primers might introduce during amplification of different parts of the ITS region in samples containing mixed templates ('environmental barcoding'). We performed in silico PCR analyses with commonly used primer combinations using various ITS datasets obtained from public databases as templates. RESULTS: Some of the ITS primers, such as ITS1-F, were hampered with a high proportion of mismatches relative to the target sequences, and most of them appeared to introduce taxonomic biases during PCR. Some primers, e.g. ITS1-F, ITS1 and ITS5, were biased towards amplification of basidiomycetes, whereas others, e.g. ITS2, ITS3 and ITS4, were biased towards ascomycetes. The assumed basidiomycete-specific primer ITS4-B only amplified a minor proportion of basidiomycete ITS sequences, even under relaxed PCR conditions. Due to systematic length differences in the ITS2 region as well as the entire ITS, we found that ascomycetes will more easily amplify than basidiomycetes using these regions as targets. This bias can be avoided by using primers amplifying ITS1 only, but this would imply preferential amplification of 'non-dikarya' fungi. CONCLUSIONS: We conclude that ITS primers have to be selected carefully, especially when used for high-throughput sequencing of environmental samples. We suggest that different primer combinations or different parts of the ITS region should be analyzed in parallel, or that alternative ITS primers should be searched for.

New perspectives in diet analysis based on DNA barcoding and parallel pyrosequencing: the trnL approach.

The development of DNA barcoding (species identification using a standardized DNA sequence), and the availability of recent DNA sequencing techniques offer new possibilities in diet analysis. DNA fragments shorter than 100-150 bp remain in a much higher proportion in degraded DNA samples and can be recovered from faeces. As a consequence, by using universal primers that amplify a very short but informative DNA fragment, it is possible to reliably identify the plant taxon that has been eaten. According to our experience and using this identification system, about 50% of the taxa can be identified to species using the trnL approach, that is, using the P6 loop of the chloroplast trnL (UAA) intron. We demonstrated that this new method is fast, simple to implement, and very robust. It can be applied for diet analyses of a wide range of phytophagous species at large scales. We also demonstrated that our approach is efficient for mammals, birds, insects and molluscs. This method opens new perspectives in ecology, not only by allowing large-scale studies on diet, but also by enhancing studies on resource partitioning among competing species, and describing food webs in ecosystems.

The dynamic evolutionary history of the bananaquit (Coereba flaveola) in the Caribbean revealed by a multigene analysis.

BACKGROUND: The bananaquit (Coereba flaveola) is a small nectivorous and frugivorous emberizine bird (order Passeriformes) that is an abundant resident throughout the Caribbean region. We used multi-gene analyses to investigate the evolutionary history of this species throughout its distribution in the West Indies and in South and Middle America. We sequenced six mitochondrial genes (3744 base pairs) and three nuclear genes (2049 base pairs) for forty-four bananaquits and three outgroup species. We infer the ancestral area of the present-day bananaquit populations, report on the species' phylogenetic, biogeographic and evolutionary history, and propose scenarios for its diversification and range expansion. RESULTS: Phylogenetic concordance between mitochondrial and nuclear genes at the base of the bananaquit phylogeny supported a West Indian origin for continental populations. Multi-gene analysis showing genetic remnants of successive colonization events in the Lesser Antilles reinforced earlier research demonstrating that bananaquits alternate periods of invasiveness and colonization with biogeographic quiescence. Although nuclear genes provided insufficient information at the tips of the tree to further evaluate relationships of closely allied but strongly supported mitochondrial DNA clades, the discrepancy between mitochondrial and nuclear data in the population of Dominican Republic suggested that the mitochondrial genome was recently acquired by introgression from Jamaica. CONCLUSION: This study represents one of the most complete phylogeographic analyses of its kind and reveals three patterns that are not commonly appreciated in birds: (1) island to mainland colonization, (2) multiple expansion phases, and (3) mitochondrial genome replacement. The detail revealed by this analysis will guide evolutionary analyses of populations in archipelagos such as the West Indies, which include islands varying in size, age, and geological history. Our results suggest that multi-gene phylogenies will permit improved comparative analysis of the evolutionary histories of different lineages in the same geographical setting, which provide replicated "natural experiments" for testing evolutionary hypotheses.

Are islands the end of the colonization road?

Ecologists have, up to now, widely regarded colonization of islands from continents as a one-way journey, mainly because of widely accepted assertions that less diverse island communities are easier to invade. However, continents present large targets and island species should be capable of making the reverse journey, considering they are the direct descendants of successful colonists and provided that they have not lost their dispersal abilities. Recent mapping of geography onto molecular phylogenies has revealed several cases of 'reverse colonization' (from islands to continents). We suggest this phenomenon warrants closer attention in ecology and biogeography. Assessing its significance will contribute to understanding the role of dispersal and establishment in biogeographic distributions and the assembly of natural biotas.

Genetic estimates of annual reproductive success in male brown bears: the effects of body size, age, internal relatedness and population density.

1. We studied male yearly reproductive success (YRS) and its determinants (phenotypic characteristics, age, population density) in two Scandinavian brown bear populations, using molecular techniques to determine paternity. 2. We found a significant difference in male YRS between the study areas, with lower YRS in the south than in the north. 3. In general, older and larger males had higher YRS. Older males may be more experienced in competition for reproduction (male dominance). Large body size is of direct benefit in male-male competition and of advantage in endurance competition for the access to females. 4. Age was relatively more important for YRS in the north and body size was more important in the south, due perhaps to differences in male age structure due to illegal killing. A single old male dominated the reproduction in the north during the study, which resulted most probably in the relatively higher importance of age in the north. In the south, with a more even male age structure, no single male was able to dominate, probably resulting in a more intense competition among males, with body size as the deciding factor. 5. Male YRS was correlated positively with population density. This may be related to the structure of the expanding bear population, with female densities declining towards the population edge. 6. Internal relatedness, a measure of genetic heterozygosity, was correlated negatively with YRS, suggesting that outbred individuals have a higher YRS. Individual heterozygosity at key or many loci may reflect male physical qualities and condition-sensitive traits, which may benefit males directly in contest or in sperm competition.

The dilemma of female mate selection in the brown bear, a species with sexually selected infanticide.

Because of differential investment in gametes between sexes, females tend to be the more selective sex. Based on this concept, we investigate mate selection in a large carnivore: the brown bear (Ursus arctos). We hypothesize that, in this species with sexually selected infanticide (SSI), females may be faced with a dilemma: either select a high-quality partner based on phenotypic criteria, as suggested by theories of mate choice, or rather mate with future potentially infanticidal males as a counter-strategy to SSI. We evaluated which male characteristics were important in paternity assignment. Among males available in the vicinity of the females, the largest, most heterozygous and less inbred and also the geographically closest males were more often the fathers of the female's next litter. We suggest that female brown bears may select the closest males as a counter-strategy to infanticide and exercise a post-copulatory cryptic choice, based on physical attributes, such as a large body size, reflecting male genetic quality. However, male-male competition either in the form of fighting before copulation or during the post-copulatory phase, in the form of sperm competition, cannot entirely be ruled out.

Genotyping errors: causes, consequences and solutions.

Although genotyping errors affect most data and can markedly influence the biological conclusions of a study, they are too often neglected. Errors have various causes, but their occurrence and effect can be limited by considering these causes in the production and analysis of the data. Procedures that have been developed for dealing with errors in linkage studies, forensic analyses and non-invasive genotyping should be applied more broadly to any genetic study. We propose a protocol for estimating error rates and recommend that these measures be systemically reported to attest the reliability of published genotyping studies.