Adélie penguin population diet monitoring by analysis of food DNA in scats

Global distribution patterns of siphonophores across horizontal and vertical oceanic gradients

Background

Siphonophores are diverse, globally distributed hydrozoans that play a central role in marine trophic webs worldwide. However, they still constitute an understudied fraction of the open ocean gelatinous taxa, mainly due to challenges related to siphonophore sampling and identification, which have led to a general knowledge gap about their diversity, distribution and abundance.

Methods

Here, we provide a global overview of the oceanic vertical distribution of siphonophores using DNA metabarcoding data from 77 bulk mesozooplankton samples collected at four different depth ranges (0-200, 200-500, 500-1000, 1000-3000 m depth) along the Atlantic, Pacific, and Indian Oceans during the MALASPINA-2010 circumnavigation expedition.

Results

We detected a total of 44 siphonophore species (which represents about one quarter of the described siphonophore species) from which 26 corresponded to Calycophores, 14 to Physonectae and 2 to Cystonectae. Our results suggest wider horizontal and vertical distributions of siphonophore species than previously described, including novel records of some species in certain oceanic basins. Also, we provide insights into the intraspecific variation of widely distributed species. Finally, we show a vertical structuring of siphonophores along the water column; Calycophores (siphonophores without pneumatophores) dominated the epipelagic (from the surface to 200 m depth) and upper mesopelagic layers (from 200 to 500 m depth), while the proportion Physonectids (siphonophores with pneumatophore) notably increased below 500 meters and were dominant at bathypelagic depths (>1000 m depth).

Conclusions

Our results support that the siphonophore community composition is vertically structured. Also, we provide insights into the potential existence of genetic variations within certain species that dominate some ocean basins or depth ranges. To our knowledge, this is the first time that DNA metabarcoding data is retrieved to study siphonophore distribution patterns, and the study provides evidence of the potential of molecular techniques to study the distribution of gelatinous organisms often destroyed in net sampling.

Molecular analysis of scats revealed diet and prey choice of grey wolves and Eurasian lynx in the contact zone between the Dinaric Mountains and the Alps

A comprehensive understanding of the dietary habits of carnivores is essential to get ecological insights into their role in the ecosystem, potential competition with other carnivorous species, and their effect on prey populations. Genetic analysis of non-invasive samples, such as scats, can supplement behavioural or microscopic diet investigations. The objective of this study was to employ DNA metabarcoding to accurately determine the prey species in grey wolf (Canis lupus) and Eurasian lynx (Lynx lynx) scat samples collected in the Julian Alps and the Dinaric Mountains, Slovenia. The primary prey of wolves were red deer (Cervus elaphus) (detected in 96% scat samples), European roe deer (Capreolus capreolus) (68%), and wild boar (Sus scrofa) (45%). A smaller portion of their diet consisted of mesocarnivores, small mammals, and domestic animals. In contrast, the lynx diet mostly consisted of European roe deer (82%) and red deer (64%). However, small mammals and domestic animals were also present in lynx diet, albeit to a lesser extent. Our findings indicate that the dietary habits of wolves and lynx are influenced by geographical location. Snapshot dietary analyses using metabarcoding are valuable for comprehending the behaviour and ecology of predators, and for devising conservation measures aimed at sustainable management of both their natural habitats and prey populations. However, to gain a more detailed understanding of wolf and lynx dietary habits and ecological impact, it would be essential to conduct long-term genetic monitoring of their diet.

VESPA: an optimized protocol for accurate metabarcoding-based characterization of vertebrate eukaryotic endosymbiont and parasite assemblages

Protocols for characterizing taxonomic assemblages by deep sequencing of short DNA barcode regions (metabarcoding) have revolutionized our understanding of microbial communities and are standardized for bacteria, archaea, and fungi. Unfortunately, comparable methods for host-associated eukaryotes have lagged due to technical challenges. Despite 54 published studies, issues remain with primer complementarity, off-target amplification, and lack of external validation. Here, we present VESPA (Vertebrate Eukaryotic endoSymbiont and Parasite Analysis) primers and optimized metabarcoding protocol for host-associated eukaryotic community analysis. Using in silico prediction, panel PCR, engineered mock community standards, and clinical samples, we demonstrate VESPA to be more effective at resolving host-associated eukaryotic assemblages than previously published methods and to minimize off-target amplification. When applied to human and non-human primate samples, VESPA enables reconstruction of host-associated eukaryotic endosymbiont communities more accurately and at finer taxonomic resolution than microscopy. VESPA has the potential to advance basic and translational science on vertebrate eukaryotic endosymbiont communities, similar to achievements made for bacterial, archaeal, and fungal microbiomes.

Effects of Sea-Ice Persistence on the Diet of Adélie Penguin (Pygoscelis adeliae) Chicks and the Trophic Differences between Chicks and Adults in the Ross Sea, Antarctica

In Antarctica, prey availability for the mesopredator Adélie penguin, Pygoscelis adeliae, depends on sea-ice dynamics. By affecting cycles of sea-ice formation and melt, climate change could thus affect penguin diet and recruitment. In the light of climate change, this raises concerns about the fate of this dominant endemic species, which plays a key role in the Antarctic food web. However, few quantitative studies measuring the effects of sea-ice persistence on the diet of penguin chicks have yet been conducted. The purpose of this study was to fill this gap by comparing penguin diets across four penguin colonies in the Ross Sea and evaluating latitudinal and interannual variation linked to different sea-ice persistence. Diet was evaluated by analysing the δ¹³C and δ¹⁵N values of penguin guano, and sea-ice persistence by means of satellite images. Isotopic values indicate that penguins consumed more krill in colonies with longer sea-ice persistence. In these colonies, the δ¹³C values of chicks were lower and closer to the pelagic chain than those of adults, suggesting that the latter apparently catch prey inshore for self-feeding and offshore for their chicks. The results indicate that sea-ice persistence is among the principal factors that influence the spatiotemporal variability of the penguins' diet.

New and legacy persistent organic pollutants (POPs) in breeding seabirds from the East Antarctic

Persistent organic pollutants (POPs) are pervasive and a significant threat to the environment worldwide. Yet, reports of POP levels in Antarctic seabirds based on blood are scarce, resulting in significant geographical gaps. Blood concentrations offer a snapshot of contamination within live populations, and have been used widely for Arctic and Northern Hemisphere seabird species but less so in Antarctica. This paper presents levels of legacy POPs (polychlorinated biphenyls (PCBs), organochlorine pesticides (OCPs) and polybrominated diphenyl ethers (PBDEs)) and novel brominated flame retardants (NBFRs) in the blood of five Antarctic seabird species breeding within Prydz Bay, East Antarctica. Legacy PCBs and OCPs were detected in all species sampled, with Adélie penguins showing comparatively high ∑PCB levels (61.1 ± 87.6 ng/g wet weight (ww)) compared to the four species of flying seabirds except the snow petrel (22.5 ± 15.5 ng/g ww), highlighting that legacy POPs are still present within Antarctic wildlife despite decades-long bans. Both PBDEs and NBFRs were detected in trace levels for all species and hexabromobenzene (HBB) was quantified in cape petrels (0.3 ± 0.2 ng/g ww) and snow petrels (0.2 ± 0.1 ng/g ww), comparable to concentrations found in Arctic seabirds. These results fill a significant data gap within the Antarctic region for POPs studies, representing a crucial step forward assessing the fate and impact of legacy POPs contamination in the Antarctic environment.

A review of applications of environmental DNA for reptile conservation and management

Reptile populations are in decline globally, with total reptile abundance halving in the past half century, and approximately a fifth of species currently threatened with extinction. Research on reptile distributions, population trends, and trophic interactions can greatly improve the accuracy of conservation listings and planning for species recovery, but data deficiency is an impediment for many species. Environmental DNA (eDNA) can detect species and measure community diversity at diverse spatio-temporal scales, and is especially useful for detection of elusive, cryptic, or rare species, making it potentially very valuable in herpetology. We aim to summarize the utility of eDNA as a tool for informing reptile conservation and management and discuss the benefits and limitations of this approach. A literature review was conducted to collect all studies that used eDNA and focus on reptile ecology, conservation, or management. Results of the literature search are summarized into key discussion points, and the review also draws on eDNA studies from other taxa to highlight methodological challenges and to identify future research directions. eDNA has had limited application to reptiles, relative to other vertebrate groups, and little use in regions with high species richness. eDNA techniques have been more successfully applied to aquatic reptiles than to terrestrial reptiles, and most (64%) of studies focused on aquatic habitats. Two of the four reptilian orders dominate the existing eDNA studies (56% Testudines, 49% Squamata, 5% Crocodilia, 0% Rhynchocephalia). Our review provides direction for the application of eDNA as an emerging tool in reptile ecology and conservation, especially when it can be paired with traditional monitoring approaches. Technologies associated with eDNA are rapidly advancing, and as techniques become more sensitive and accessible, we expect eDNA will be increasingly valuable for addressing key knowledge gaps for reptiles.

WOS, Steinberg DK, Swadling KM, Tarling GA, Thorpe SE, Veytia D, Ward P, Weldrick CK, Yang G. Status, Change, and Futures of Zooplankton in the Southern Ocean

In the Southern Ocean, several zooplankton taxonomic groups, euphausiids, copepods, salps and pteropods, are notable because of their biomass and abundance and their roles in maintaining food webs and ecosystem structure and function, including the provision of globally important ecosystem services. These groups are consumers of microbes, primary and secondary producers, and are prey for fishes, cephalopods, seabirds, and marine mammals. In providing the link between microbes, primary production, and higher trophic levels these taxa influence energy flows, biological production and biomass, biogeochemical cycles, carbon flux and food web interactions thereby modulating the structure and functioning of ecosystems. Additionally, Antarctic krill (Euphausia superba) and various fish species are harvested by international fisheries. Global and local drivers of change are expected to affect the dynamics of key zooplankton species, which may have potentially profound and wide-ranging implications for Southern Ocean ecosystems and the services they provide. Here we assess the current understanding of the dominant metazoan zooplankton within the Southern Ocean, including Antarctic krill and other key euphausiid, copepod, salp and pteropod species. We provide a systematic overview of observed and potential future responses of these taxa to a changing Southern Ocean and the functional relationships by which drivers may impact them. To support future ecosystem assessments and conservation and management strategies, we also identify priorities for Southern Ocean zooplankton research.

Molecular Diet Analysis of Adélie Penguins (Pygoscelis adeliae) in the Ross Sea Using Fecal DNA

Simple Summary

The diet of Adélie penguins, Pygoscelis adeliae, in the Ross Sea was studied applying quantitative polymerase chain reaction (qPCR) protocols to their feces. Two krill species (Euphausia superba and Euphausia crystallorophias) and notothenioid fish (mainly Pleuragramma antarctica and Pagothenia borchgrevinki) were among the most abundant components of the diet. The composition of the two krill species and notothenioid fish was found to be strongly related to the geographic characteristics of the Ross Sea.

Abstract

The diet of Adélie penguins, Pygoscelis adeliae, is a useful indicator in understanding the ecological conditions of their habitats. The diets of Adélie penguins were studied using metabarcoding and quantitative PCR (qPCR) analyses of fecal DNA from seven habitats along the Ross Sea region. Using metabarcoding analysis with dual universal primers (18Sv9 and miniFish), the overall diet composition and detailed information about piscine prey were clearly elucidated. It was found that two krill species (Euphausia superba and Euphausia crystallorophias) and notothenioid fish were the most abundant in the diets of Adélie penguins. Among the notothenioid prey, Pleuragramma antarctica (56.50%) and Pagothenia borchgrevinki (18.21%) were the two most abundant species. qPCR analysis showed a significant geographic difference in the composition of main prey. Penguins inhabiting outbound parts of the Ross Sea (Capes Adare (CA) and Duke of York Island (DY)) mainly preyed on E. superba, without any significant changes in prey composition. By contrast, those inhabiting the inbound parts of the Ross Sea (Edmonson Point (EP) and Inexpressible Island (II)) preyed on E. crystallorophias and notothenioid fish rather than E. superba. Compared with the outbound habitats, prey compositions for penguins inhabiting the inbound regions were significantly different year to year, which was presumably due to the food availability based on the annual environmental and meteorological conditions of the coastal region along with the inbound parts of the Ross Sea.

DNA-metabarcoding reveals the importance of gelatinous zooplankton in the diet of Pandalus borealis, a keystone species in the Arctic

Information about the dietary composition of a species is crucial to understanding their position and role in the food web. Increasingly molecular approaches such as DNA-metabarcoding are used in studying trophic relations, not least because they may alleviate problems such as low taxonomic resolution or underestimation of digestible taxa in the diet. Here, we used DNA-metabarcoding with universal primers for cytochrome c oxidase I (COI), to study the diet composition of the Northern shrimp (Pandalus borealis), an Arctic keystone species with large socio-economic importance. Across locations, jellyfish and chaetognaths were the most important components in the diet of P. borealis, jointly accounting for 40-60% of the total read abundance. This dietary importance of gelatinous zooplankton contrasts sharply with published results based on SCA. At the same time, diet composition differed between fjord and shelf locations, pointing to different food webs supporting P. borealis in these two systems. Our study underscores the potential of molecular approaches to provide new insights into the diet of marine invertebrates that are difficult to obtain with traditional methods, and calls for a revision of the role of gelatinous zooplankton in the diet of the key Arctic species P. borealis, and in extension, Arctic food webs.

Advances and Limitations of Next Generation Sequencing in Animal Diet Analysis

Diet analysis is a critical content of animal ecology and the diet analysis methods have been constantly improving and updating. Contrary to traditional methods of high labor intensity and low resolution, the next generation sequencing (NGS) approach has been suggested as a promising tool for dietary studies, which greatly improves the efficiency and broadens the application range. Here we present a framework of adopting NGS and DNA metabarcoding into diet analysis, and discuss the application in aspects of prey taxa composition and structure, intra-specific and inter-specific trophic links, and the effects of animal feeding on environmental changes. Yet, the generation of NGS-based diet data and subsequent analyses and interpretations are still challenging with several factors, making it possible still not as widely used as might be expected. We suggest that NGS-based diet methods must be furthered, analytical pipelines should be developed. More application perspectives, including nutrient geometry, metagenomics and nutrigenomics, need to be incorporated to encourage more ecologists to infer novel insights on they work.

The diet of early birds based on modern and fossil evidence and a new framework for its reconstruction

Birds are some of the most diverse organisms on Earth, with species inhabiting a wide variety of niches across every major biome. As such, birds are vital to our understanding of modern ecosystems. Unfortunately, our understanding of the evolutionary history of modern ecosystems is hampered by knowledge gaps in the origin of modern bird diversity and ecosystem ecology. A crucial part of addressing these shortcomings is improving our understanding of the earliest birds, the non-avian avialans (i.e. non-crown birds), particularly of their diet. The diet of non-avian avialans has been a matter of debate, in large part because of the ambiguous qualitative approaches that have been used to reconstruct it. Here we review methods for determining diet in modern and fossil avians (i.e. crown birds) as well as non-avian theropods, and comment on their usefulness when applied to non-avian avialans. We use this to propose a set of comparable, quantitative approaches to ascertain fossil bird diet and on this basis provide a consensus of what we currently know about fossil bird diet. While no single approach can precisely predict diet in birds, each can exclude some diets and narrow the dietary possibilities. We recommend combining (i) dental microwear, (ii) landmark-based muscular reconstruction, (iii) stable isotope geochemistry, (iv) body mass estimations, (v) traditional and/or geometric morphometric analysis, (vi) lever modelling, and (vii) finite element analysis to reconstruct fossil bird diet accurately. Our review provides specific methodologies to implement each approach and discusses complications future researchers should keep in mind. We note that current forms of assessment of dental mesowear, skull traditional morphometrics, geometric morphometrics, and certain stable isotope systems have yet to be proven effective at discerning fossil bird diet. On this basis we report the current state of knowledge of non-avian avialan diet which remains very incomplete. The ancestral dietary condition in non-avian avialans remains unclear due to scarce data and contradictory evidence in Archaeopteryx. Among early non-avian pygostylians, Confuciusornis has finite element analysis and mechanical advantage evidence pointing to herbivory, whilst Sapeornis only has mechanical advantage evidence indicating granivory, agreeing with fossilised ingested material known for this taxon. The enantiornithine ornithothoracine Shenqiornis has mechanical advantage and pedal morphometric evidence pointing to carnivory. In the hongshanornithid ornithuromorph Hongshanornis only mechanical advantage evidence indicates granivory, but this agrees with evidence of gastrolith ingestion in this taxon. Mechanical advantage and ingested fish support carnivory in the songlingornithid ornithuromorph Yanornis. Due to the sparsity of robust dietary assignments, no clear trends in non-avian avialan dietary evolution have yet emerged. Dietary diversity seems to increase through time, but this is a preservational bias associated with a predominance of data from the Early Cretaceous Jehol Lagerstätte. With this new framework and our synthesis of the current knowledge of non-avian avialan diet, we expect dietary knowledge and evolutionary trends to become much clearer in the coming years, especially as fossils from other locations and climates are found. This will allow for a deeper and more robust understanding of the role birds played in Mesozoic ecosystems and how this developed into their pivotal role in modern ecosystems.

Multiple species-specific molecular markers using nanofluidic array as a tool to detect prey DNA from carnivore scats

Large carnivore feeding ecology plays a crucial role for management and conservation for predators and their prey. One of the keys to this kind of research is to identify the species composition in the predator diet, for example, prey determination from scat content. DNA-based methods applied to detect prey in predators' scats are viable alternatives to traditional macroscopic approaches, showing an increased reliability and higher prey detection rate. Here, we developed a molecular method for prey species identification in wolf (Canis lupus) scats using multiple species-specific marker loci on the cytochrome b gene for 18 target species. The final panel consisted of 80 assays, with a minimum of four markers per target species, and that amplified specifically when using a high-throughput Nanofluidic array technology (Fluidigm Inc.). As a practical example, we applied the method to identify target prey species DNA in 80 wolf scats collected in Sweden. Depending on the number of amplifying markers required to obtain a positive species call in a scat, the success in determining at least one prey species from the scats ranged from 44% to 92%. Although we highlight the need to evaluate the optimal number of markers for sensitive target species detection, the developed method is a fast and cost-efficient tool for prey identification in wolf scats and it also has the potential to be further developed and applied to other areas and large carnivores as well.

Metabarcoding confirms the opportunistic foraging behaviour of Atlantic bluefin tuna and reveals the importance of gelatinous prey

Studies of the diet, feeding habits and trophic activity of top marine predators are essential for understanding their trophodynamics. The main direct method used for such studies thus far has been morphological inventories of stomach contents. This approach presents limitations such as missing gelatinous prey, which are usually digested too quickly to be detectable. Here, we analysed the stomachs of 48 Atlantic bluefin tuna (Thunnus thynnus, approximately 15 to 60 kg, including juveniles and adult fishes) collected from the Mediterranean Sea through the metabarcoding of two gene regions (cytochrome c oxidase subunit I (COI) and the ribosomal 18S-V1V2 region). The identified prey taxa and their relative read abundances (RRAs) estimated using COI results were in line with the findings of morphologically based inventories simultaneously performed on the same set of tuna samples. In both cases (and with the same rankings), the prey taxa included anchovy (Engraulis encrasicolus, here detected in more than 80% of samples, RRA = 43%), sardine (Sardina pilchardus, also approximately 80%, RRA = 30%), sprat (Sprattus sprattus, approximately 66%, RRA = 8%), mackerel (Scomber colias, approximately 44%, RRA = 7%) and cephalopods (approximately 15%, RRA = 1.4%). Another striking result was the detection, based on 18S (with which vertebrates were detected as the most abundant group, RRA = 61.6%), of a high prevalence and diversity of gelatinous organisms (RRA = 27.1%), including cnidarians (6.7%), salps (11.7%), and ctenophores (8.7%), the latter increasing with the size of the predator. These results thus support the hypothesis of the role of gelatinous prey in the diet of Atlantic bluefin tuna, suggesting that this species is even more generalist and opportunistic than previously thought. This study further confirms that DNA metabarcoding can be a powerful tool for assessing the diet and trophodynamics of top marine predators.

Metabarcoding options to study eukaryotic endoparasites of birds

There is growing interest in the study of avian endoparasite communities, and metabarcoding is a promising approach to complement more conventional or targeted methods. In the case of eukaryotic endoparasites, phylogenetic diversity is extreme, with parasites from 4 kingdoms and 11 phyla documented in birds. We addressed this challenge by comparing different primer sets across 16 samples from 5 bird species. Samples consisted of blood, feces, and controlled mixes with known proportions of bird and nematode DNA. Illumina sequencing revealed that a 28S primer set used in combination with a custom blocking primer allowed detection of various plasmodiid parasites and filarioid nematodes in the blood, coccidia in the feces, as well as two potentially pathogenic fungal groups. When tested on the controlled DNA mixes, these primers also increased the proportion of nematode DNA by over an order of magnitude. An 18S primer set, originally designed to exclude metazoan sequences, was the most effective at reducing the relative number of avian DNA sequences and was the only one to detect Trypanosoma in the blood. Expectedly, however, it did not allow nematode detection and also failed to detect avian malaria parasites. This study shows that a 28S set including a blocking primer allows detection of several major and very diverse bird parasite clades, while reliable amplification of all major parasite groups may require a combination of markers. It helps clarify options for bird parasite metabarcoding, according to priorities in terms of the endoparasite clades and the ecological questions researchers wish to focus on.

Fish as predators and prey: DNA-based assessment of their role in food webs

Fish are both consumers and prey, and as such part of a dynamic trophic network. Measuring how they are trophically linked, both directly and indirectly, to other species is vital to comprehend the mechanisms driving alterations in fish communities in space and time. Moreover, this knowledge also helps to understand how fish communities respond to environmental change and delivers important information for implementing management of fish stocks. DNA-based methods have significantly widened our ability to assess trophic interactions in both marine and freshwater systems and they possess a range of advantages over other approaches in diet analysis. In this review we provide an overview of different DNA-based methods that have been used to assess trophic interactions of fish as consumers and prey. We consider the practicalities and limitations, and emphasize critical aspects when analysing molecular derived trophic data. We exemplify how molecular techniques have been employed to unravel food web interactions involving fish as consumers and prey. In addition to the exciting opportunities DNA-based approaches offer, we identify current challenges and future prospects for assessing fish food webs where DNA-based approaches will play an important role.

High-resolution multi-marker DNA metabarcoding reveals sexual dietary differentiation in a bird with minor dimorphism

Although sexual dietary differentiation is well known in birds, it is usually linked with significant morphological dimorphism between males and females, with lower differentiation reported in sexually monomorphic or only slightly dimorphic species. However, this may be an artifact of poor taxonomic resolution achieved in most conventional dietary studies, which may be unable to detect subtle intraspecific differentiation in prey consumption. Here, we show the power of multi‐marker metabarcoding to address these issues, focusing on a slightly dimorphic generalist passerine, the black wheatear Oenanthe leucura. Using markers from four genomic regions (18S, 16S, COI, and trnL), we analyzed fecal droppings collected from 93 adult black wheatears during the breeding season. We found that sexes were rather similar in bill and body features, though males had a slightly thicker bill and longer wings and tail than females. Diet was dominated in both sexes by a very wide range of arthropod species and a few fleshy fruits, but the overall diet diversity was higher for males than females, and there was a much higher frequency of occurrence of ants in female (58%) than male (29%) diets. We hypothesize that the observed sexual differentiation was likely related to females foraging closer to their offspring on abundant prey, while males consumed a wider variety of prey while foraging more widely. Overall, our results suggest that dietary sexual differentiation in birds may be more widespread than recognized at present and that multi‐marker DNA metabarcoding is a particularly powerful tool to unveiling such differences.

Latitudinal, sex and inter-specific differences in mercury and other trace metal concentrations in Adélie and Emperor penguins in the Ross Sea, Antarctica

We sought to determine mercury (Hg) and other trace metal concentrations in Adélie (Pygoscelis adeliae) and emperor penguin (Aptenodytes forsteri) breast feathers from the Ross Sea, Antarctica, and relate those concentrations to the trophic position and the habitats in which each of these species forage. Adélie penguin feathers from the southern Ross Sea colonies were higher in Hg than those sampled further north in the Ross Sea, potentially due to greater exposure to local sources, such as volcanism. Female Adélie penguins had lower feather total Hg concentrations than males. This may reflect female penguin's capacity to eliminate Hg through the egg development and laying process, or the larger and/or older prey items that male birds can consume, reflected by their higher trophic position. Emperor penguins have higher Hg concentrations than Adélie penguins which is also partially explained by Adélie penguins feeding at lower trophic levels than emperor penguins.

Partitioning of diet between species and life history stages of sympatric and cryptic snappers (Lutjanidae) based on DNA metabarcoding

Lutjanus erythropterus and L. malabaricus are sympatric, sister taxa that are important to fisheries throughout the Indo-Pacific. Their juveniles are morphologically indistinguishable (i.e. cryptic). A DNA metabarcoding dietary study was undertaken to assess the diet composition and partitioning between the juvenile and adult life history stages of these two lutjanids. Major prey taxa were comprised of teleosts and crustaceans for all groups except adult L. erythropterus, which instead consumed soft bodied invertebrates (e.g. tunicates, comb jellies and medusae) as well as teleosts, with crustaceans being notably absent. Diet composition was significantly different among life history stages and species, which may be associated with niche habitat partitioning or differences in mouth morphology within adult life stages. This study provides the first evidence of diet partitioning between cryptic juveniles of overlapping lutjanid species, thus providing new insights into the ecological interactions, habitat associations, and the specialised adaptations required for the coexistence of closely related species. This study has improved our understanding of the differential contributions of the juvenile and adult diets of these sympatric species within food webs. The diet partitioning reported in this study was only revealed by the taxonomic resolution provided by the DNA metabarcoding approach and highlights the potential utility of this method to refine the dietary components of reef fishes more generally.

First genetic record of the non-native muzzled blenny Omobranchus punctatus (Teleostei: Blenniidae) in the Atlantic Coast of Central and South America

In this study we sequenced two mitochondrial (COI and 16S rRNA) and one nuclear (18S rRNA) gene fragment of an introduced muzzled blenny (Omobranchus punctatus) specimen collected from the Orinoco Delta (Gulf of Paria estuary) in Venezuela. This is the first genetic data generated for this species' introduced range in Central and South America, suggesting an introduction from the Indian Ocean.

Questioning the proverb 'more haste, less speed': classic versus metabarcoding approaches for the diet study of a remote island endemic gecko

Dietary studies can reveal valuable information on how species exploit their habitats and are of particular importance for insular endemics conservation as these species present higher risk of extinction. Reptiles are often neglected in island systems, principally the ones inhabiting remote areas, therefore little is known on their ecological networks. The Selvagens gecko Tarentola (boettgeri) bischoffi, endemic to the remote and integral reserve of Selvagens Archipelago, is classified as Vulnerable by the Portuguese Red Data Book. Little is known about this gecko’s ecology and dietary habits, but it is assumed to be exclusively insectivorous. The diet of the continental Tarentola species was already studied using classical methods. Only two studies have used next-generation sequencing (NGS) techniques for this genus thus far, and very few NGS studies have been employed for reptiles in general. Considering the lack of information on its diet and the conservation interest of the Selvagens gecko, we used morphological and DNA metabarcoding approaches to characterize its diet. The traditional method of morphological identification of prey remains in faecal pellets collected over a longer period was compared with metabarcoding of samples collected during rapid surveys. Molecular results revealed that this species is a generalist, feeding on invertebrate, plant and vertebrate items, whereas the morphological approaches were unable to detect the latter two. These results opened up new questions on the ecological role of the Selvagens gecko that deserves to be further explored, such as the possible predation on seabirds, plant services or trophic competition with the sympatric Madeira lizard Teira dugesii. Metabarcoding identified a greater diversity of dietary items at higher taxonomic resolution, but morphological identification enabled calculation of relative abundances and biomasses of ingested arthropods, and detected a dietary shift on invertebrate preys between seasons. Results of this study highlight the global applicability of rapid metabarcoding surveys for understudied taxa on remote islands that are difficult to access. We recommend using the metabarcoding approach, even if ‘speedy’ sampling only is possible, but we must highlight that disregarding long-term ecological data may lead to ‘hasty’ conclusion.

Species in the faeces: DNA metabarcoding as a method to determine the diet of the endangered yellow-eyed penguin

Abstract Context. Diet variability is a significant driver of seabird decline; however, data on seabird diet composition and trends have been affected by changes in precision and resolution owing to the evolution of different sampling methods over time. We investigated the effectiveness of applying a passive molecular diet method using faeces obtained from the endangered yellow-eyed penguin. Aims. To assess the feasibility of applying DNA metabarcoding methods to yellow-eyed penguin faeces to evaluate diet, and to compare the reliability of diet results derived from adults and chicks, and from latrine versus fresh faecal samples. Methods. We collected 313 faecal samples from yellow-eyed penguins resident on the Otago coast of New Zealand from October 2016 to August 2017. We used polymerase chain reaction (PCR) with mitochondrial 16S cephalopod and chordate primers to amplify prey DNA present in the faecal samples, and tested the completeness of our assembled reference databases based on previous diet research. Amplified prey DNA sequences were then assigned to taxa from our reference databases by using QIIME2. Key results. Mitochondrial 16S chordate PCR primers were effective at identifying 29 fish taxa, with 98.3% of amplified sequences being identified to species or genus level in 193 samples (61.7% collected). There was no significant difference in the number, occurrence or proportion of ray-finned fish prey DNA sequences derived from fresh samples or latrines. Mitochondrial 16S cephalopod PCR primers classified 1.98% of amplified DNA sequences as targets, with 96.5% of these target sequences being identified to species or genus level in 48 samples (15.3% collected), and five taxa identified. Conclusions. We recommend the collection of latrine samples to enable long-term monitoring of the diet of yellow-eyed penguins, which will optimise the trade-off between wildlife disturbance and dietary resolution. Further refinement is needed to identify cephalopod dietary components for yellow-eyed penguins, because our cephalopod primers were not as specific as those used for ray-finned fishes, amplifying a large number (>98%) of non-cephalopod species. Implications. DNA metabarcoding offers a robust and comprehensive alternative to other, more intrusive, seabird diet-assessment methods, but still requires parallel studies to provide critical information on prey size, true diet composition and diet quality.

Advancing the integration of multi-marker metabarcoding data in dietary analysis of trophic generalists

The application of DNA metabarcoding to dietary analysis of trophic generalists requires using multiple markers in order to overcome problems of primer specificity and bias. However, limited attention has been given to the integration of information from multiple markers, particularly when they partly overlap in the taxa amplified, and vary in taxonomic resolution and biases. Here, we test the use of a mix of universal and specific markers, provide criteria to integrate multi‐marker metabarcoding data and a python script to implement such criteria and produce a single list of taxa ingested per sample. We then compare the results of dietary analysis based on morphological methods, single markers, and the proposed combination of multiple markers. The study was based on the analysis of 115 faeces from a small passerine, the Black Wheatears (Oenanthe leucura). Morphological analysis detected far fewer plant taxa (12) than either a universal 18S marker (57) or the plant trnL marker (124). This may partly reflect the detection of secondary ingestion by molecular methods. Morphological identification also detected far fewer taxa (23) than when using 18S (91) or the arthropod markers IN16STK (244) and ZBJ (231), though each method missed or underestimated some prey items. Integration of multi‐marker data provided far more detailed dietary information than any single marker and estimated higher frequencies of occurrence of all taxa. Overall, our results show the value of integrating data from multiple, taxonomically overlapping markers in an example dietary data set.

DNA metabarcoding to assess diet partitioning and feeding strategies in generalist vertebrate predators: a case study on three syntopic lacertid lizards from Morocco

DNA metabarcoding is a fast and simple alternative to traditional microscopy methods, which have been the main tool for identification of prey in dietary studies of lizards. In this study, we applied a metabarcoding approach based on COI and 16S rRNA amplicons to assess diet partitioning and feeding strategies in three syntopic lizards from Taza, Morocco: Scelarcis perspicillata chabanaudi, Scelarcis perspicillata pellegrini and Podarcis vaucheri. In order to avoid competition, these lizards are expected to consume different prey species because they occupy distinct trophic niches, use different foraging strategies and express different dorsal pigmentation patterns. Given the spotted pattern of S. p. chabanaudi, we hypothesize a sit-and-wait foraging strategy with a less diverse diet and a higher consumption of mobile prey relative to the striped S. p. pellegrini and P. vaucheri which, as potential active foragers, are expected to have a higher diet diversity. Previous diet assessments using microscopy on faecal remains seem to contradict these expectations. Our results show that, as expected, the diet of S. p. chabanaudi is less diverse than the diet of S. p. pellegrini. Regarding P. vaucheri, our dietary data are consistent with the hypothesis that this species behaves as an active forager, owing to its high niche overlap with S. p. pellegrini. Advantages and limitations of molecular barcoding compared with the microscopy approach to the analysis of lizard diets are discussed.

DNA metabarcoding assays reveal a diverse prey assemblage for Mobula rays in the Bohol Sea, Philippines

Diet studies provide base understanding of trophic structure and are a valuable initial step for many fields of marine ecology, including conservation and fisheries biology. Considerable complexity in marine trophic structure can exist due to the presence of highly mobile species with long life spans. Mobula rays are highly mobile, large, planktivorous elasmobranchs that are frequently caught either directly or as bycatch in fisheries, which, combined with their conservative life history strategy, makes their populations susceptible to decline in intensely fished regions. Effective management of these iconic and vulnerable species requires an understanding of the diets that sustain them, which can be difficult to determine using conventional sampling methods. We use three DNA metabarcode assays to identify 44 distinct taxa from the stomachs (n = 101) of four sympatric Mobula ray species (Mobula birostris, Mobula tarapacana, Mobula japanica, and Mobula thurstoni) caught over 3 years (2013-2015) in a direct fishery off Bohol in the Philippines. The diversity and incidence of bony fishes observed in ray diets were unprecedented. Nevertheless, rays showed dietary overlap, with krill (Euphausia) dominating their diet. Our results provide a more detailed assessment of sympatric ray diets than was previously described and reveal the complexity that can exist in food webs at critical foraging habitats.

Multi-locus DNA metabarcoding of zooplankton communities and scat reveal trophic interactions of a generalist predator

To understand the ecosystem dynamics that underpin the year-round presence of a large generalist consumer, the Bryde's whale (Balaenoptera edeni brydei), we use a DNA metabarcoding approach and systematic zooplankton surveys to investigate seasonal and regional changes in zooplankton communities and if whale diet reflects such changes. Twenty-four zooplankton community samples were collected from three regions throughout the Hauraki Gulf, New Zealand, over two temperature regimes (warm and cool seasons), as well as 20 samples of opportunistically collected Bryde's whale scat. Multi-locus DNA barcode libraries were constructed from 18S and COI gene fragments, representing a trade-off between identification and resolution of metazoan taxa. Zooplankton community OTU occurrence and relative read abundance showed regional and seasonal differences based on permutational analyses of variance in both DNA barcodes, with significant changes in biodiversity indices linked to season in COI only. In contrast, we did not find evidence that Bryde's whale diet shows seasonal or regional trends, but instead indicated clear prey preferences for krill-like crustaceans, copepods, salps and ray-finned fishes independent of prey availability. The year-round presence of Bryde's whales in the Hauraki Gulf is likely associated with the patterns of distribution and abundance of these key prey items.

Biostatistical Options for Quantitative Diet Analysis

Sufficient statistics knowledge is crucial for the correct design of a research plan. The elucidations of results are interpretive only if appropriate statistical methods are applied. Statistical strategies are a particular approach to demonstrate complicated information in broad and explicable conclusions. The emergence of biostatistical approaches for diet evaluation has improved the accuracy of diet estimation, and different methodologies of data integration promise to magnify our understanding of ecological communities. The present study aimed to compile multiple statistical methods used for diet analysis. More specifically, the significant analysis used in diet assessment, central expectations, and preferences related to each measure was conceptualized. In addition, the ability of each test to evaluate diversity, richness, differentiation, fluctuation, similarity, and quantification of multiple diet items was summarized. Moreover, different options were proposed for researchers to select the appropriate statistical tests. This study covers a framework, aim, and understanding of the statistical test methods of diet analysis.

Evaluation of plant contamination in metabarcoding diet analysis of a herbivore

Fecal DNA metabarcoding is currently used in various fields of ecology to determine animal diets. Contamination of non-food DNA from complex field environments is a considerable challenge to the reliability of this method but has rarely been quantified. We evaluated plant DNA contamination by sequencing the chloroplast trnL P6 loop region from food-controlled geese feces. The average percentage of contaminant sequences per sample was 1.86%. According to the results of generalized linear models, the probability of contamination was highest in samples placed in wet soil. The proportion of contaminant sequences was lowest at the earliest sampling point and was slightly higher in samples placed in open conditions. Exclusion of rare OTUs (operational taxonomic units) was effective for obtaining reliable dietary data from the obtained sequences, and a 1% cutoff reduced the percentage of contaminated samples to less than 30%. However, appropriate interpretation of the barcoding results considering inevitable contamination is an important issue to address. We suggest the following procedures for fecal sampling and sequence data treatment to increase the reliability of DNA metabarcoding diet analyses: (i) Collect samples as soon as possible after deposition, (ii) avoid samples from deposits on wet soil, and (iii) exclude rare OTUs from diet composition estimations.

Sex-specific prey partitioning in breeding piscivorous birds examined via a novel, noninvasive approach

Piscivorous birds frequently display sex-specific differences in their hunting and feeding behavior, which lead to diverging impacts on prey populations. Cormorants (Phalacrocoracidae), for example, were previously studied to examine dietary differences between the sexes and males were found to consume larger fish in coastal areas during autumn and winter. However, information on prey partitioning during breeding and generally on sex-specific foraging in inland waters is missing. Here, we assess sex-specific prey choice of Great Cormorants (Phalacrocorax carbo) during two subsequent breeding seasons in the Central European Alpine foreland, an area characterized by numerous stagnant and flowing waters in close proximity to each other. We developed a unique, noninvasive approach and applied it to regurgitated pellets: molecular cormorant sexing combined with molecular fish identification and fish-length regression analysis performed on prey hard parts. Altogether, 364 pellets delivered information on both, bird sex, and consumed prey. The sexes differed significantly in their overall prey composition, even though Perca fluviatilis, Rutilus rutilus, and Coregonus spp. represented the main food source for both. Albeit prey composition did not indicate the use of different water bodies by the sexes, male diet was characterized by higher prey diversity within a pellet and the consumption of larger fish. The current findings show that female and male cormorants to some extent target the available prey spectrum at different levels. Finally, the comprehensive and noninvasive approach has great potential for application in studies of other piscivorous bird species.

Genetic and genomic monitoring with minimally invasive sampling methods

The decreasing cost and increasing scope and power of emerging genomic technologies are reshaping the field of molecular ecology. However, many modern genomic approaches (e.g., RAD-seq) require large amounts of high-quality template DNA. This poses a problem for an active branch of conservation biology: genetic monitoring using minimally invasive sampling (MIS) methods. Without handling or even observing an animal, MIS methods (e.g., collection of hair, skin, faeces) can provide genetic information on individuals or populations. Such samples typically yield low-quality and/or quantities of DNA, restricting the type of molecular methods that can be used. Despite this limitation, genetic monitoring using MIS is an effective tool for estimating population demographic parameters and monitoring genetic diversity in natural populations. Genetic monitoring is likely to become more important in the future as many natural populations are undergoing anthropogenically driven declines, which are unlikely to abate without intensive adaptive management efforts that often include MIS approaches. Here, we profile the expanding suite of genomic methods and platforms compatible with producing genotypes from MIS, considering factors such as development costs and error rates. We evaluate how powerful new approaches will enhance our ability to investigate questions typically answered using genetic monitoring, such as estimating abundance, genetic structure and relatedness. As the field is in a period of unusually rapid transition, we also highlight the importance of legacy data sets and recommend how to address the challenges of moving between traditional and next-generation genetic monitoring platforms. Finally, we consider how genetic monitoring could move beyond genotypes in the future. For example, assessing microbiomes or epigenetic markers could provide a greater understanding of the relationship between individuals and their environment.

Combined effects of ocean acidification and temperature on planula larvae of the moon jellyfish Aurelia coerulea

Rapidly rising levels of atmospheric CO₂ have caused two environmental stressors, ocean acidification and seawater temperature increases, which represent major abiotic threats to marine organisms. Here, we investigated for the first time the combined effects of ocean acidification and seawater temperature increases on the behavior, survival, and settlement of the planula larvae of Aurelia coerulea, which is considered a nuisance species around the world. Three pH levels (8.1, 7.7 and 7.3) and two temperature levels (24 °C and 27 °C) were used in the present study. There were no interactive effects of temperature and pH on the behavior, survival, and settlement of planula larvae of A. coerulea. We found that the swimming speed and mortality of the planula larvae of A. coerulea were significantly affected by temperature, and low pH significantly affected settlement. Planula larvae of A. coerulea from the elevated temperature treatment moved faster and showed higher mortality than those at the control temperature. The settlement rate of A. coerulea planulae was significantly higher at the pH level of 7.3 than at other pH levels. These results suggest that seawater temperature increase, rather than reduced pH, was the main stress factor affecting the survival of A. coerulea planulae. Overall, the planula larvae of the common jellyfish A. coerulea appeared to be resistant to ocean acidification, but may be negatively affected by future seawater temperature increases.

Review: the energetic value of zooplankton and nekton species of the Southern Ocean

Understanding the energy flux through food webs is important for estimating the capacity of marine ecosystems to support stocks of living resources. The energy density of species involved in trophic energy transfer has been measured in a large number of small studies, scattered over a 40-year publication record. Here, we reviewed energy density records of Southern Ocean zooplankton, nekton and several benthic taxa, including previously unpublished data. Comparing measured taxa, energy densities were highest in myctophid fishes (ranging from 17.1 to 39.3 kJ g^-1 DW), intermediate in crustaceans (7.1 to 25.3 kJ g^-1 DW), squid (16.2 to 24.0 kJ g^-1 DW) and other fish families (14.8 to 29.9 kJ g^-1 DW), and lowest in jelly fish (10.8 to 18.0 kJ g^-1 DW), polychaetes (9.2 to 14.2 kJ g^-1 DW) and chaetognaths (5.0-11.7 kJ g^-1 DW). Data reveals differences in energy density within and between species related to size, age and other life cycle parameters. Important taxa in Antarctic food webs, such as copepods, squid and small euphausiids, remain under-sampled. The variability in energy density of Electrona antarctica was likely regional rather than seasonal, although for many species with limited data it remains difficult to disentangle regional and seasonal variability. Models are provided to estimate energy density more quickly using a species' physical parameters. It will become increasingly important to close knowledge gaps to improve the ability of bioenergetic and food web models to predict changes in the capacity of Antarctic ecosystems to support marine life.

Counting with DNA in metabarcoding studies: How should we convert sequence reads to dietary data?

Advances in DNA sequencing technology have revolutionized the field of molecular analysis of trophic interactions, and it is now possible to recover counts of food DNA sequences from a wide range of dietary samples. But what do these counts mean? To obtain an accurate estimate of a consumer's diet should we work strictly with data sets summarizing frequency of occurrence of different food taxa, or is it possible to use relative number of sequences? Both approaches are applied to obtain semi-quantitative diet summaries, but occurrence data are often promoted as a more conservative and reliable option due to taxa-specific biases in recovery of sequences. We explore representative dietary metabarcoding data sets and point out that diet summaries based on occurrence data often overestimate the importance of food consumed in small quantities (potentially including low-level contaminants) and are sensitive to the count threshold used to define an occurrence. Our simulations indicate that using relative read abundance (RRA) information often provides a more accurate view of population-level diet even with moderate recovery biases incorporated; however, RRA summaries are sensitive to recovery biases impacting common diet taxa. Both approaches are more accurate when the mean number of food taxa in samples is small. The ideas presented here highlight the need to consider all sources of bias and to justify the methods used to interpret count data in dietary metabarcoding studies. We encourage researchers to continue addressing methodological challenges and acknowledge unanswered questions to help spur future investigations in this rapidly developing area of research.

DNA barcoding reveals seasonal shifts in diet and consumption of deep-sea fishes in wedge-tailed shearwaters

The foraging ecology of pelagic seabirds is difficult to characterize because of their large foraging areas. In the face of this difficulty, DNA metabarcoding may be a useful approach to analyze diet compositions and foraging behaviors. Using this approach, we investigated the diet composition and its seasonal variation of a common seabird species on the Ogasawara Islands, Japan: the wedge-tailed shearwater Ardenna pacifica. We collected fecal samples during the prebreeding (N = 73) and rearing (N = 96) periods. The diet composition of wedge-tailed shearwater was analyzed by Ion Torrent sequencing using two universal polymerase chain reaction primers for the 12S and 16S mitochondrial DNA regions that targeted vertebrates and mollusks, respectively. The results of a BLAST search of obtained sequences detected 31 and 1 vertebrate and mollusk taxa, respectively. The results of the diet composition analysis showed that wedge-tailed shearwaters frequently consumed deep-sea fishes throughout the sampling season, indicating the importance of these fishes as a stable food resource. However, there was a marked seasonal shift in diet, which may reflect seasonal changes in food resource availability and wedge-tailed shearwater foraging behavior. The collected data regarding the shearwater diet may be useful for in situ conservation efforts. Future research that combines DNA metabarcoding with other tools, such as data logging, may provide further insight into the foraging ecology of pelagic seabirds.

Dropping Hints: Estimating the diets of livestock in rangelands using DNA metabarcoding of faeces

The introduction of domesticated animals into new environments can lead to considerable ecological disruption, and it can be difficult to predict their impact on the new ecosystem. In this study, we use faecal metabarcoding to characterize the diets of three ruminant taxa in the rangelands of south-western New South Wales, Australia. Our study organisms included goats (Capraaegagrushircus) and two breeds of sheep (Ovisaries): Merinos, which have been present in Australia for over two hundred years, and Dorpers, which were introduced in the 1990s. We used High-Throughput Sequencing methods to sequence therbcLand ITS2 genes of plants in the faecal samples, and identified the samples using the GenBank and BOLD online databases, as well as a reference collection of sequences from plants collected in the study area. We found that the diets of all three taxa were dominated by the family Malvaceae, and that the Dorper diet was more diverse than the Merino diet at both the family and the species level. We conclude that Dorpers, like Merinos, are potentially a threat to some vulnerable species in the rangelands of New South Wales.

Evaluation of an automated protocol for efficient and reliable DNA extraction of dietary samples

Molecular techniques have become an important tool to empirically assess feeding interactions. The increased usage of next‐generation sequencing approaches has stressed the need of fast DNA extraction that does not compromise DNA quality. Dietary samples here pose a particular challenge, as these demand high‐quality DNA extraction procedures for obtaining the minute quantities of short‐fragmented food DNA. Automatic high‐throughput procedures significantly decrease time and costs and allow for standardization of extracting total DNA. However, these approaches have not yet been evaluated for dietary samples. We tested the efficiency of an automatic DNA extraction platform and a traditional CTAB protocol, employing a variety of dietary samples including invertebrate whole‐body extracts as well as invertebrate and vertebrate gut content samples and feces. Extraction efficacy was quantified using the proportions of successful PCR amplifications of both total and prey DNA, and cost was estimated in terms of time and material expense. For extraction of total DNA, the automated platform performed better for both invertebrate and vertebrate samples. This was also true for prey detection in vertebrate samples. For the dietary analysis in invertebrates, there is still room for improvement when using the high‐throughput system for optimal DNA yields. Overall, the automated DNA extraction system turned out as a promising alternative to labor‐intensive, low‐throughput manual extraction methods such as CTAB. It is opening up the opportunity for an extensive use of this cost‐efficient and innovative methodology at low contamination risk also in trophic ecology.

The diet of a nocturnal pelagic predator, the Bulwer's petrel, across the lunar cycle

The lunar cycle is believed to strongly influence the vertical distribution of many oceanic taxa, with implications for the foraging behaviour of nocturnal marine predators. Most studies to date testing lunar effects on foraging have focused on predator activity at-sea, with some birds and marine mammals demonstrating contrasting behavioural patterns, depending on the lunar-phase. However, to date no study has focused on how the lunar cycle might actually affect predator-prey interactions in the upper layers of the ocean. Here, we tested whether the diet of the predominantly nocturnal pelagic predator, the Bulwer’s petrel (Bulweria bulwerii) would change throughout the lunar cycle, using molecular analysis to augment detection and taxonomic resolution of prey collected from stomach-contents. We found no evidence of dietary shifts in species composition or diversity, with Bulwer’s petrel always consuming a wide range of mesopelagic species. Other co-variables potentially affecting light availability at-sea, such as percentage of cloud cover, did not confound our results. Moreover, many of the species found are thought not to reach the sea-surface. Our findings reveal that nocturnal predators are probably more specialized than previously assumed, irrespective of ambient-light, but also reveal deficiencies in our current understanding of species vertical distribution and predation-dynamics at-sea.

Diet analysis in piscivorous birds: What can the addition of molecular tools offer?

In trophic studies on piscivorous birds, it is vital to know which kind of dietary sample provides the information of interest and how the prey can be identified reliably and efficiently. Often, noninvasively obtained dietary samples such as regurgitated pellets, feces, and regurgitated fish samples are the preferred source of information. Fish prey has usually been identified via morphological analysis of undigested hard parts, but molecular approaches are being increasingly used for this purpose. What remains unknown, however, is which dietary sample type is best suited for molecular diet analysis and how the molecular results compare to those obtained by morphological analysis. Pellets, feces, and regurgitated fish samples of Great Cormorants (Phalacrocorax carbo sinensis) were examined for prey using both morphological hard part analysis and molecular prey identification. The sample types and methods were compared regarding number of species detected (overall and per sample) as well as the prey species composition and its variability among individual samples. Via molecular analysis, significantly higher numbers of prey species were detected in pellets, feces, and fish samples. Of the three sample types, pellets contained the most comprehensive trophic information and could be obtained with the lowest sampling effort. Contrastingly, dietary information obtained from feces was least informative and most variable. For all sample types, the molecular approach outperformed morphological hard part identification regarding the detectable prey spectrum and prey species composition. We recommend the use of pellets in combination with molecular prey identification to study the diet of piscivorous birds.

Effect of marker choice and thermal cycling protocol on zooplankton DNA metabarcoding studies

DNA metabarcoding is a promising approach for rapidly surveying biodiversity and is likely to become an important tool for measuring ecosystem responses to environmental change. Metabarcoding markers need sufficient taxonomic coverage to detect groups of interest, sufficient sequence divergence to resolve species, and will ideally indicate relative abundance of taxa present. We characterized zooplankton assemblages with three different metabarcoding markers (nuclear 18S rDNA, mitochondrial COI, and mitochondrial 16S rDNA) to compare their performance in terms of taxonomic coverage, taxonomic resolution, and correspondence between morphology‐ and DNA‐based identification. COI amplicons sequenced on separate runs showed that operational taxonomic units representing >0.1% of reads per sample were highly reproducible, although slightly more taxa were detected using a lower annealing temperature. Mitochondrial COI and nuclear 18S showed similar taxonomic coverage across zooplankton phyla. However, mitochondrial COI resolved up to threefold more taxa to species compared to 18S. All markers revealed similar patterns of beta‐diversity, although different taxa were identified as the greatest contributors to these patterns for 18S. For calanoid copepod families, all markers displayed a positive relationship between biomass and sequence reads, although the relationship was typically strongest for 18S. The use of COI for metabarcoding has been questioned due to lack of conserved primer‐binding sites. However, our results show the taxonomic coverage and resolution provided by degenerate COI primers, combined with a comparatively well‐developed reference sequence database, make them valuable metabarcoding markers for biodiversity assessment.