Discovering hidden biodiversity: the use of complementary monitoring of fish diet based on DNA barcoding in freshwater ecosystems

Application of high-throughput sequencing to analyze the diet of black croaker (Atrobucca nibe) in southern coastal waters of Zhejiang

This study investigated the feeding habits of black croaker (Atrobucca nibe) to fully understand its role in the food web in southern coastal waters of Zhejiang, China. In total, 225 black croakers were collected in a bottom trawl survey, of which 195 were used for the stomach contents analysis and 30 for high-throughput sequencing (HTS). Both approaches showed that the main prey was shrimp, with Alpheus japonicus being the most predominant. In addition, black croaker also fed on other benthic species, suggesting that it was a demersal predator. Feeding intensity also varied with black croaker body length, as did prey preference, from small-sized species to larger species with increases in fish size; in addition, prey species were more diverse in the large-bodied prey group. Canonical correspondence analysis indicated that the feeding habit of black croaker was most significantly correlated with temperature and dissolved oxygen. HTS was shown to be a feasible high precision semi-quantitative analytical approach, although it must be combined with morphological methods to obtain more complete feeding data; thus, it could provide an effective scientific basis for feeding ecology studies of marine organisms.

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.

Aquatic Organisms Research with DNA Barcodes

Since their inception, DNA barcodes have become a powerful tool for understanding the biodiversity and biology of aquatic species, with multiple applications in diverse fields such as food security, fisheries, environmental DNA, conservation, and exotic species detection. Nevertheless, most aquatic ecosystems, from marine to freshwater, are understudied, with many species disappearing due to environmental stress, mostly caused by human activities. Here we highlight the progress that has been made in studying aquatic organisms with DNA barcodes, and encourage its further development in assisting sustainable use of aquatic resources and conservation.

eDNA metabarcoding for biodiversity assessment, generalist predators as sampling assistants

With an accelerating negative impact of anthropogenic actions on natural ecosystems, non-invasive biodiversity assessments are becoming increasingly crucial. As a consequence, the interest in the application of environmental DNA (eDNA) survey techniques has increased. The use of eDNA extracted from faeces from generalist predators, have recently been described as "biodiversity capsules" and suggested as a complementary tool for improving current biodiversity assessments. In this study, using faecal samples from two generalist omnivore species, the Eurasian badger and the red fox, we evaluated the applicability of eDNA metabarcoding in determining dietary composition, compared to macroscopic diet identification techniques. Subsequently, we used the dietary information obtained to assess its contribution to biodiversity assessments. Compared to classic macroscopic techniques, we found that eDNA metabarcoding detected more taxa, at higher taxonomic resolution, and proved to be an important technique to verify the species identification of the predator from field collected faeces. Furthermore, we showed how dietary analyses complemented field observations in describing biodiversity by identifying consumed flora and fauna that went unnoticed during field observations. While diet analysis approaches could not substitute field observations entirely, we suggest that their integration with other methods might overcome intrinsic limitations of single techniques in future biodiversity surveys.

Plant-herbivorous insect networks: who is eating what revealed by long barcodes using high-throughput sequencing and Trinity assembly

Interactions between plants and insects are among the most important life functions for all organism at a particular natural community. Usually a large number of samples are required to identify insect diets in food web studies. Previously, Sanger sequencing and next generation sequencing (NGS) with short DNA barcodes were used, resulting in low species-level identification; meanwhile the costs of Sanger sequencing are expensive for metabarcoding together with more samples. Here, we present a fast and effective sequencing strategy to identify larvae of Lepidoptera and their diets at the same time without increasing the cost on Illumina platform in a single HiSeq run, with long-multiplex-metabarcoding (COI for insects, rbcL, matK, ITS and trnL for plants) obtained by Trinity assembly (SHMMT). Meanwhile, Sanger sequencing (for single individuals) and NGS (for polyphagous) were used to verify the reliability of the SHMMT approach. Furthermore, we show that SHMMT approach is fast and reliable, with most high-quality sequences of five DNA barcodes of 63 larvae individuals (54 species) recovered (full length of 100% of the COI gene and 98.3% of plant DNA barcodes) using Trinity assembly (up-sized to 1015 bp). For larvae diets identification, 95% are reliable; the other 5% failed because their guts were empty. The diets identified by SHMMT approach are 100% consistent with the host plants that the larvae were feeding on during our collection. Our study demonstrates that SHMMT approach is reliable and cost-effective for insect-plants network studies. This will facilitate insect-host plant studies that generally contain a huge number of samples.

First Gut Content Analysis of 4th Instar Midge Larvae (Diptera: Chronomidae) In Large-Scale Weirs Using a DNA Meta-Barcoding Approach

Chironomidae larvae play an important role in the food chain of river ecosystems in Korea, where it is dominant. However, detailed information on the diet of Chironomidae larvae are still lacking. The purpose of this study was to identify the gut contents of 4th instar larvae of a Chironomidae inhabiting four large-scale weirs (Sejong Weir, Juksan Weir, Gangjeong-Goryeong Weir, and Dalseong Weir) using a DNA meta-barcoding approach. We found that dominant Operational Taxonomic Unit (OUT) was assigned to Paractinolaimus sp. (Nematoda), and the sub-dominant OTU was assigned to Dicrotendipes fumidus (Chironomidae). The most common OTUs among the individuals included phytoplankton, such as Tetrahymena sp., D. armatus, Pseudopediastrum sp., Tetradesmus dimorphus, Biddulphia tridens, and Desmodesmus spp. We calculated the selectivity index (E’) and provided scientific evidence that Chironomidae larvae have a significant preference (E’ > 0.5) for Desmodesmus armatus, E. minima, and T. dimorphus, while it does not show preference for other species found in its gut. Differences in physico-chemical factors, such as water quality, nutrients, Chl-a, and carbon concentrations, resulting from anthropogenic impacts (i.e., construction of large-scale weirs) as well as the particle size of prey organisms (small-sized single cell) and effects of chemicals (chemokinesis) could affect the feeding behavior of Chironomidae larvae.

A molecular approach to identifying the relationship between resource use and availability in Eurasian otters (Lutra lutra)

In South Korea, the Eurasian otter (Lutra lutra (Linnaeus, 1758)), a semi-aquatic carnivore, is found mainly in lower order streams that tend to have a low abundance of preferred prey fish species. To investigate the relationship between resource use and availability, we used DNA barcoding to identify otter diet items in 24 otter spraints (faeces) from 16 sites along the Nakdong River basin from 4 to 6 June 2014. At these sites fish availability was assessed using scoop nets and casting nets. Fish formed the bulk of otter diet, which included also frogs, mammals, and reptiles. By DNA barcoding (success rate: 72.38%), we identified 79 prey items from 105 bone remains. The diet comprised mostly fish, but frogs, mammals, and reptiles were also identified. The fish fauna and otter diet composition differed significantly. Across the study sites, members of the Cyprinidae dominated in netted samples, but occurred less frequently in otter diet. Because most Cyprinidae are fast swimmers, otters also fed on benthic fishes and frogs, suggesting limited foraging flexibility in otters and specialization on more slowly moving prey.

Discrimination of Spatial Distribution of Aquatic Organisms in a Coastal Ecosystem Using eDNA

The nonlinearity and complexity of coastal ecosystems often cause difficulties when analyzing spatial and temporal patterns of ecological traits. Environmental DNA (eDNA) monitoring has provided an alternative to overcoming the aforementioned issues associated with classical monitoring. We determined aquatic community taxonomic composition using eDNA based on a meta-barcoding approach that characterizes the general ecological features in the Gwangyang Bay coastal ecosystem. We selected the V9 region of the 18S rDNA gene (18S V9), primarily because of its broad range among eukaryotes. Our results produced more detailed spatial patterns in the study area previously categorized (inner bay, main channel of the bay and outer bay) by Kim et al. (2019). Specifically, the outer bay zone was clearly identified by CCA using genus-level identification of aquatic organisms based on meta-barcoding data. We also found significant relationships between environmental factors. Therefore, eDNA monitoring based on meta-barcoding approach holds great potential as a complemental monitoring tool to identify spatial taxonomic distribution patterns in coastal areas.

Responses of fish assemblage structure to large-scale weir construction in riverine ecosystems

Worldwide, increasing amounts of dams and weirs have been established in rivers in recent decades, often with drastic effects on their ecosystems. Between late 2009 and 2011, 16 large-scale dams were built in the main channels of the four largest rivers in South Korea, eight of these along the main channel of Nakdong River (300 km, 520 km in total). We studied the effect of these constructions on the fish community in the riparian zone based on intensive fish field surveys conducted in the littoral zone during 2007-2017, analysis of fishermen's catch data and molecular analysis of the diet of the keystone species, largemouth bass (Micropterus salmoides). Our results, based on RDA and GLM analyses, showed that environmental characteristics and fish species composition changed significantly with dam construction. Total fish abundance and biomass decreased after the start of the weir construction and fish richness decreased with time. The abundance and biomass of exotic fish decreased during construction but recovered afterwards, leading to higher proportions of exotics. Dominance and Shannon indices did not change significantly from before to after construction, while the Evenness index showed a significant decreasing trend. A comparison of the diet composition of largemouth bass showed decreasing genetic variation after construction. The changes in Functional Feeding Group (FFG) of the benthic communities at the study sites did not coincide with FFG changes in the prey items of largemouth bass, indicating a transition in feeding mode from before to after construction. In conclusion, the endemic and native fish species were most sensitive and showed lower resilience to disturbance by the large-scale dam construction than the translocated and exotic species, and the lake-like ecosystems after construction markedly improved the competitive capacity of these exotic fish over the native and endemic species in the riparian zone of the river.

Herbivore-Diet Analysis Based on Illumina MiSeq Sequencing: The Potential Use of an ITS2-Barcoding Approach to Establish Qualitative and Quantitative Predictions of Diet Composition of Mongolian Sheep

DNA-barcoding approaches to estimate the diet compositions of grazing animals have received significant attention, and particularly when combined with next-generation sequencing, these techniques have substantially improved in recent years. In this study, the identity and species composition of plant material ingested by Mongolian sheep were estimated through the use of 350 bp ITS2 gene sequences of the vegetation found in fecal samples. Four diets were formulated using varying amounts of eight plant species that are common in the grasslands of northern China. Sixteen Mongolian sheep were taken from pastures and randomly assigned to four groups, and each group received one of four diets. Each sheep was randomly assigned to one of 16 confinement pens and fed its respective diet for 12 consecutive days. Fecal samples were removed from each pen from days 7-12, preserved, and composited for each pen. All herbage species included in the daily diets were detected in each fecal sample, with the exception of Phragmites australis. Moreover, 12 additional different plant species were retrieved from feces of the experimental sheep. The obtained data provided preliminary support for the use of the ITS2 barcode to determine which plants were consumed. Moreover, the proportions of the herbage DNA sequences recovered from sheep feces and those of the herbage masses in the daily diets did not completely match. These results indicate that the non-Gramineae DNA sequences amplified with ITS2 primers (including those of Chenopodium album, Artemisia scoparia, Artemisia tanacetifolia, and Medicago sativa) far exceeded those of the Gramineae species (including Leymus chinensis and Puccinellia distans), which constitute the largest share of the experimental diets. A significant positive correlation (Spearman's ρ = 0.376, P = 0.003) between the actual herbage mass proportions in the experimental diets and the herbage-DNA-sequence proportions provided sufficiently favorable support for the further investigation of DNA barcoding for the quantification of plants in feces. A significant regression coefficient was found between the relative DNA-sequence proportions of L. chinensis ( R² = 0.82, P < 0.0001), P. distans ( R² = 0.64, P = 0.0017), and C. album ( R² = 0.98, P < 0.0001) and their respective herbage mass proportions. The quantitative relationship can be expressed by the linear-regression equations y = 0.90 x - 0.22, y = 0.98 x - 0.03, and y = 5.00 x - 0.25, respectively. Thus, these results demonstrate that dietary-DNA-barcoding methods exhibited potential in providing valuable quantitative information regarding food-item components. However, it should be noted that this explorative data needs to be further improved by using additional genes and by creating a sophisticated reference database, thus enhancing both quality and accuracy of the obtained results.

Next-generation sequencing analysis of Pardosa pseudoannulata's diet composition in different habitats

Spiders are the most common and predominant predators in terrestrial ecosystems. The predatory behavior of spiders affects the energy flow across the food web within an ecosystem. Traditiaonal methods for analyzing spider diets such as field observation, anatomy and faeces analysis are not suitable for spider experiments due to spiders’ special dietary behavior. The molecular method based on the specific primers of prey DNA seems to be inefficient either in spite of its wide application in diet analysis. As the next-generation sequencing (NGS) technology becomes prevalent in many different areas, several cases of the NGS-based analysis of mammal diets have been published. This study analyzed the diet differences of Pardosa pseudoannulata (Araneae: Lycosidae) in four habitats (a wetland, a tea plantation, an alpine meadow and a paddy field) by using the NGS technology, combined with the DNA barcode method. The results suggested that the Pardosa pseudoannulata feed on a broad range of prey, and 7 orders and 24 families of insects were detected in the four investigated habitats. Moreover, it is found that the diet diversity of Pardosa pseudoannulata is greatly influenced by their living environments and seasons. In a nutshell, this study established an NGS-based methodology for spider diets analysis, and the results provided some basic materials to inform the protection and utilization of the Pardosa pseudoannulata as a potential eco-friendly predator against pests.

A from-benchtop-to-desktop workflow for validating HTS data and for taxonomic identification in diet metabarcoding studies

The main objective of this work was to develop and validate a robust and reliable "from-benchtop-to-desktop" metabarcoding workflow to investigate the diet of invertebrate-eaters. We applied our workflow to faecal DNA samples of an invertebrate-eating fish species. A fragment of the cytochrome c oxidase I (COI) gene was amplified by combining two minibarcoding primer sets to maximize the taxonomic coverage. Amplicons were sequenced by an Illumina MiSeq platform. We developed a filtering approach based on a series of nonarbitrary thresholds established from control samples and from molecular replicates to address the elimination of cross-contamination, PCR/sequencing errors and mistagging artefacts. This resulted in a conservative and informative metabarcoding data set. We developed a taxonomic assignment procedure that combines different approaches and that allowed the identification of ~75% of invertebrate COI variants to the species level. Moreover, based on the diversity of the variants, we introduced a semiquantitative statistic in our diet study, the minimum number of individuals, which is based on the number of distinct variants in each sample. The metabarcoding approach described in this article may guide future diet studies that aim to produce robust data sets associated with a fine and accurate identification of prey items.

Application of DNA metabarcoding on faeces to identify European catfish Silurus glanis diet

In this study, the results of conventional stomach-content analysis are compared with the recent DNA metabarcoding approach on faeces to identify fish species consumed by non-native European catfish Silurus glanis in the Garonne River (south-western France), with a special emphasis on anadromous prey. Fourteen prey species were identified in the stomach contents or faeces, including four anadromous fish species. Despite higher intestine than stomach emptiness, more species were identified through faecal analysis (11 of 14) than through stomach-content analysis (five of 14) suggesting that DNA metabarcoding on faeces is an efficient, non-intrusive technique to study the diet of predatory fishes.

On the methodology of feeding ecology in fish

Feeding ecology explains predator’s preference to some preys over others in their habitat and their competitions thereof. The subject, as a functional and applied biology, is highly neglected, and in case of fish, a uniform and consistent methodology is absent. The currently practiced methods are largely centred on mathematical indices and highly erroneous because of non-uniform outcomes. Therefore, it requires a relook into the subject to elucidate functional contributions and to make it more comparable and comprehensive science. In this article, approachable methodological strategies have been forwarded in three hierarchical steps, namely, food occurrence, feeding biology and interpretative ecology. All these steps involve wide ranges of techniques, within the scope of ecology but not limited to, and traverse from narrative to functional evolutionary ecology. The first step is an assumption-observation practice to assess food of fish, followed by feeding biology that links morphological, histological, cytological, bacteriological or enzymological correlations to preferred food in the environment. Interpretative ecology is the higher level of analysis in which the outcomes are tested and discussed against evolutionary theories. A description of possible pedagogics on the methods of feeding ecological studies has also been forwarded.