The use of coded PCR primers enables high-throughput sequencing of multiple homolog amplification products by 454 parallel sequencing

Effects of livestock on arthropod biodiversity in Iberian holm oak savannas revealed by metabarcoding

Increasing food production while avoiding negative impacts on biodiversity constitutes one of the main challenges of our time. Traditional silvopastoral systems like Iberian oak savannas ("dehesas") set an example, where free-range livestock has been reared for centuries while preserving a high natural value. Nevertheless, factors decreasing productivity need to be addressed, one being acorn losses provoked by pest insects. An increased and focalized grazing by livestock on infested acorns would kill the larvae inside and decrease pest numbers, but increased livestock densities could have undesired side effects on ground arthropod communities as a whole. We designed an experimental setup including areas under trees with livestock exclosures of different ages (short-term: 1-year exclusion, long-term: 10-year exclusion), along with controls (continuous grazing), using DNA metabarcoding (mitochondrial markers COI and 16S) to rapidly assess arthropod communities' composition. Livestock removal quickly increased grass cover and arthropod taxonomic richness and diversity, which was already higher in short-term (1-year exclosures) than beneath the canopies of control trees. Interestingly, arthropod diversity was not highest at long-term exclosures (≥10 years), although their community composition was the most distinct. Also, regardless of treatment, we found that functional diversity strongly correlated with the vegetation structure, being higher at trees beneath which there was higher grass cover and taller herbs. Overall, the taxonomic diversity peak at short term exclosures would support the intermediate disturbance hypothesis, which relates it with the higher microhabitat heterogeneity at moderately disturbed areas. Thus, we propose a rotatory livestock management in dehesas: plots with increased grazing should co-exist with temporal short-term exclosures. Ideally, a few long-term excluded areas should be also kept for the singularity of their arthropod communities. This strategy would make possible the combination of biological pest control and arthropod conservation in Iberian dehesas.

Diet components associated with specific bacterial taxa shape overall gut community compositions in omnivorous African viverrids

Gut bacterial communities provide flexibility to hosts during dietary changes. Despite the increasing number of studies exploring the associations between broader dietary guilds of mammalian hosts and their gut bacteria, it is generally unclear how diversity and variability in consumed diets link to gut bacterial taxa in wild non-primate mammals, particularly in omnivores. Here, we contribute to filling this gap by exploring consumed diets and gut bacterial community compositions with metabarcoding of faecal samples for two African mammals, Civettictis civetta and Genetta spp., from the family Viverridae. For each individual sample, we characterised bacterial communities and identified dietary taxa by sequencing vertebrate, invertebrate and plant markers. This led us to establish diet compositions that diverged from what has previously been found from visual identification methods. Specifically, while the two genera have been categorised into the same dietary guild, we detected more animal dietary items than plant items in C. civetta, while in Genetta spp., we observed the opposite. We further found that individuals with similar diets have similar gut bacterial communities within both genera. This association tended to be driven by specific links between dietary items and gut bacterial genera, rather than communities as a whole, implying diet-driven selection for specific gut microbes in individual wild hosts. Our findings underline the importance of molecular tools for improving characterisations of omnivorous mammalian diets and highlight the opportunities for simultaneously disentangling links between diets and gut symbionts. Such insights can inform robustness and flexibility in host-microbe symbioses to dietary change associated with seasonal and habitat changes.

Multiplexing PCR allows the identification of within-species genetic diversity in ancient eDNA

Sedimentary ancient DNA (sedaDNA) has rarely been used to obtain population-level data due to either a lack of taxonomic resolution for the molecular method used, limitations in the reference material or inefficient methods. Here, we present the potential of multiplexing different PCR primers to retrieve population-level genetic data from sedaDNA samples. Vaccinium uliginosum (Ericaceae) is a widespread species with a circumpolar distribution and three lineages in present-day populations. We searched 18 plastid genomes for intraspecific variable regions and developed 61 primer sets to target these. Initial multiplex PCR testing resulted in a final set of 38 primer sets. These primer sets were used to analyse 20 lake sedaDNA samples (11,200 cal. yr BP to present) from five different localities in northern Norway, the Alps and the Polar Urals. All known V. uliginosum lineages in these regions and all primer sets could be recovered from the sedaDNA data. For each sample on average 28.1 primer sets, representing 34.15 sequence variants, were recovered. All sediment samples were dominated by a single lineage, except three Alpine samples which had co-occurrence of two different lineages. Furthermore, lineage turnover was observed in the Alps and northern Norway, suggesting that present-day phylogeographical studies may overlook past genetic patterns. Multiplexing primer is a promising tool for generating population-level genetic information from sedaDNA. The relatively simple method, combined with high sensitivity, provides a scalable method which will allow researchers to track populations through time and space using environmental DNA.

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.

Host-gut microbiota interactions shape parasite infections in farmed Atlantic salmon

Animals and their associated microbiota share long evolutionary histories. However, it is not always clear how host genotype and microbiota interact to affect phenotype. We applied a hologenomic approach to explore how host-microbiota interactions shape lifetime growth and parasite infection in farmed Atlantic salmon (Salmo salar). Multi-omics data sets were generated from the guts of 460 salmon, 82% of which were naturally infected with an intestinal cestode. A single Mycoplasma bacterial strain, MAG01, dominated the gut metagenome of large, non-parasitized fish, consistent with previous studies showing high levels of Mycoplasma in the gut microbiota of healthy salmon. While small and/or parasitized salmon also had high abundance of MAG01, we observed increased alpha diversity in these individuals, driven by increased frequency of low-abundance Vibrionaceae and other Mycoplasma species that carried known virulence genes. Colonization by one of these cestode-associated Mycoplasma strains was associated with host individual genomic variation in long non-coding RNAs. Integrating the multi-omic data sets revealed coordinated changes in the salmon gut mRNA transcriptome and metabolome that correlated with shifts in the microbiota of smaller, parasitized fish. Our results suggest that the gut microbiota of small and/or parasitized fish is in a state of dysbiosis that partly depends on the host genotype, highlighting the value of using a hologenomic approach to incorporate the microbiota into the study of host-parasite dynamics.IMPORTANCEStudying host-microbiota interactions through the perspective of the hologenome is gaining interest across all life sciences. Intestinal parasite infections are a huge burden on human and animal health; however, there are few studies investigating the role of the hologenome during parasite infections. We address this gap in the largest multi-omics fish microbiota study to date using natural cestode infection of farmed Atlantic salmon. We find a clear association between cestode infection, salmon lifetime growth, and perturbation of the salmon gut microbiota. Furthermore, we provide the first evidence that the genetic background of the host may partly determine how the gut microbiota changes during parasite-associated dysbiosis. Our study therefore highlights the value of a hologenomic approach for gaining a more in-depth understanding of parasitism.

Interplay of competition and facilitation in grazing succession by migrant Serengeti herbivores

Competition, facilitation, and predation offer alternative explanations for successional patterns of migratory herbivores. However, these interactions are difficult to measure, leaving uncertainty about the mechanisms underlying body-size-dependent grazing-and even whether succession occurs at all. We used data from an 8-year camera-trap survey, GPS-collared herbivores, and fecal DNA metabarcoding to analyze the timing, arrival order, and interactions among migratory grazers in Serengeti National Park. Temporal grazing succession is characterized by a "push-pull" dynamic: Competitive grazing nudges zebra ahead of co-migrating wildebeest, whereas grass consumption by these large-bodied migrants attracts trailing, small-bodied gazelle that benefit from facilitation. "Natural experiments" involving intense wildfires and rainfall respectively disrupted and strengthened these effects. Our results highlight a balance between facilitative and competitive forces in co-regulating large-scale ungulate migrations.

AmpliSAS and AmpliHLA: Web Server and Local Tools for MHC Typing of Non-model Species and Human Using NGS Data

AmpliSAS and AmpliHLA are tools for automatic genotyping of MHC genes from high-throughput sequencing data. AmpliSAS is designed specifically to analyze amplicon sequencing data from non-model species and it is able to perform de novo genotyping without any previous knowledge of the reference alleles. AmpliHLA is a human specific version; it performs HLA typing by comparing sequenced variants against human reference alleles from the IMGT/HLA database. Both tools are available in AmpliSAT web-server as well as scripts for local/server installation. Here we describe the installation and deployment of AmpliSAS and AmpliHLA Perl scripts and dependencies on a local or a server computer. We will show how to run them in the command line using as examples four genotyping protocols: the first two use amplicon sequencing data to genotype the MHC genes of a passerine bird and human respectively; the third and fourth present the HLA typing of a human cell line starting from RNA and exome sequencing data respectively.

Iron age genomic data from Althiburos - Tunisia renew the debate on the origins of African taurine cattle

The Maghreb is a key region for understanding the dynamics of cattle dispersal and admixture with local aurochs following their earliest domestication in the Fertile Crescent more than 10,000 years ago. Here, we present data on autosomal genomes and mitogenomes obtained for four archaeological specimens of Iron Age (∼2,800 cal BP-2,000 cal BP) domestic cattle from the Eastern Maghreb, i.e. Althiburos (El Kef, Tunisia). D-loop sequences were obtained for an additional eight cattle specimens from this site. Maternal lineages were assigned to the elusive R and ubiquitous African-T1 haplogroups found in two and ten Althiburos specimens, respectively. Our results can be explained by post-domestication hybridization of Althiburos cattle with local aurochs. However, we cannot rule out an independent domestication in North Africa considering the shared ancestry of Althiburos cattle with the pre-domestic Moroccan aurochs and present-day African taurine cattle.

The topological nature of tag jumping in environmental DNA metabarcoding studies

Metabarcoding of environmental DNA constitutes a state-of-the-art tool for environmental studies. One fundamental principle implicit in most metabarcoding studies is that individual sample amplicons can still be identified after being pooled with others-based on their unique combinations of tags-during the so-called demultiplexing step that follows sequencing. Nevertheless, it has been recognized that tags can sometimes be changed (i.e., tag jumping), which ultimately leads to sample crosstalk. Here, using four DNA metabarcoding data sets derived from the analysis of soils and sediments, we show that tag jumping follows very specific and systematic patterns. Specifically, we find a strong correlation between the number of reads in blank samples and their topological position in the tag matrix (described by vertical and horizontal vectors). This observed spatial pattern of artefactual sequences could be explained by polymerase activity, which leads to the exchange of the 3' tag of single stranded tagged sequences through the formation of heteroduplexes with mixed barcodes. Importantly, tag jumping substantially distorted our data sets-despite our use of methods suggested to minimize this error. We developed a topological model to estimate the noise based on the counts in our blanks, which suggested that 40%-80% of the taxa in our soil and sedimentary samples were likely false positives introduced through tag jumping. We highlight that the amount of false positive detections caused by tag jumping strongly biased our community analyses.

Genomic microbiome analyses of surface sand samples from the Kyzyl-Kum Desert (Uzbekistan): characterization and comparative study

The Kyzyl-Kum Desert extends over an area of 300,000 Km², in the region bordering Kazakhstan, Uzbekistan and Turkmenistan and is mainly covered by sand dunes. The Kyzyl-Kum desert is also known for its large deposits of minerals of economic interests, the exploitation of which is affecting the local ecosystem and increasing the desertification. We examined the bacterial biodiversity of surface sand samples from several sites from the Kyzyl-Kum desert using pyrosequencing of PCR amplified bacterial 16S rRNA genes from total extracted soil DNA. We also examined several physicochemical parameters of the sand samples to investigate any possible correlations between bacterial community structure and environmental drivers. The predominant bacterial phyla present in the samples were found to belong to members of the Actinobacteria, Proteobacteria and Bacteroidetes. The most abundant genera in our samples were found to belong to the Arthrobacter, Adhaeribacter and Roseomonas genera. We found that the relative abundance of members belonging to the Actinobacteria phylum, commonly found in desertic areas, increase in abundance in sites with higher content of organic matter and sulfur, whereas members of the Proteobacteria and Bacteroidetes phyla seems to diminish in abundance in coarse silt and fine-grained soils and those rich in magnesium, suggesting that those parameters might influence the bacterial community composition in this desertic area. This study is the first to provide new insights into the prokaryotic community composition from this unusual desert site.

A simplified protocol for DNA extraction from FTA cards for faecal microbiome studies

As metagenomic studies continue to increase in size and complexity, they are often required to incorporate data from geographically isolated locations or longitudinal time samples. This represents a technical challenge, given that many of the commonly used methods used for sample collection, storage, and DNA extraction are sensitive to differences related to the time, storage and chemistry involved. FTA cards have been previously proposed as a simple, reliable and cost-efficient method for the preservation of animal faecal microbiomes. In this study, we report a simplified extraction methodology for recovering microbiome DNA from faeces stored on FTA cards and compare its performance to a common alternative means of characterising such microbiomes; namely, immediate freezing of the faeces followed by DNA extraction using the Qiagen PowerSoil DNA isolation kit. Our results show that overall the application of our simplified DNA extraction methodology yields microbial community results that have higher diversity and an expanded core microbiome than that found using the PowerSoil methodology. This suggests that the FTA card extraction method presented here is a viable alternative for metagenomic studies using faecal material when traditional freeze-based storage methods are not feasible.

Sample Preparation for Fungal Community Analysis by High-Throughput Sequencing of Barcode Amplicons

Fungal species participate in vast numbers of processes in the landscape around us. However, their cryptic mycelial growth, inside various substrates and in highly diverse species assemblages, has been a major obstacle to thorough analysis of fungal communities, hampering exhaustive description of the fungal kingdom. Technological developments allowing rapid, high-throughput sequencing of mixed communities from many samples at once are currently having a tremendous impact in fungal community ecology. Universal DNA extraction followed by amplification and sequencing of fungal species-level barcodes such as the nuclear internal transcribed spacer (ITS) region now enables identification and relative quantification of fungal community members across well-replicated experimental settings.Here, we present the sample preparation procedure presently used in our laboratory for fungal community analysis by high-throughput sequencing of amplified ITS2 markers. We focus on the procedure optimized for studies of total fungal communities in humus-rich soils, wood, and litter. However, this procedure can be applied to other sample types and markers. We focus on the laboratory-based part of sample preparation, i.e., the procedure from the point where samples enter the laboratory until amplicons are submitted for sequencing. Our procedure comprises four main parts: (1) universal DNA extraction, (2) optimization of PCR conditions, (3) production of tagged ITS amplicons, and (4) preparation of the multiplexed amplicon pool to be sequenced. The presented procedure is independent of the specific high-throughput sequencing technology used, which makes it highly versatile.

Species composition of shoreline wolf spider communities vary with salinity, but their diets vary with wrack inflow

Wolf spiders are typically the most common group of arthropod predators on both lake and marine shorelines because of the high prey availability in these habitats. However, shores are also harsh environments due to flooding and, in proximity to marine waters, to toxic salinity levels. Here, we describe the spider community, prey availabilities, and spider diets between shoreline sites with different salinities, albeit with comparatively small differences (5‰ vs. 7‰). Despite the small environmental differences, spider communities between lower and higher saline sites showed an almost complete species turnover. At the same time, differences in prey availability or spider gut contents did not match changes in spider species composition but rather changed with habitat characteristics within a region, where spiders collected at sites with thick wrack beds had a different diet than sites with little wrack. These data suggest that shifts in spider communities are due to habitat characteristics other than prey availabilities, and the most likely candidate restricting species in high salinity would be saline sensitivity. At the same time, species absence from low-saline habitats remains unresolved.

Deconstructing a Syndrome: Genomic Insights Into PCOS Causal Mechanisms and Classification

Polycystic ovary syndrome (PCOS) is among the most common disorders in women of reproductive age, affecting up to 15% worldwide, depending on the diagnostic criteria. PCOS is characterized by a constellation of interrelated reproductive abnormalities, including disordered gonadotropin secretion, increased androgen production, chronic anovulation, and polycystic ovarian morphology. It is frequently associated with insulin resistance and obesity. These reproductive and metabolic derangements cause major morbidities across the lifespan, including anovulatory infertility and type 2 diabetes (T2D). Despite decades of investigative effort, the etiology of PCOS remains unknown. Familial clustering of PCOS cases has indicated a genetic contribution to PCOS. There are rare Mendelian forms of PCOS associated with extreme phenotypes, but PCOS typically follows a non-Mendelian pattern of inheritance consistent with a complex genetic architecture, analogous to T2D and obesity, that reflects the interaction of susceptibility genes and environmental factors. Genomic studies of PCOS have provided important insights into disease pathways and have indicated that current diagnostic criteria do not capture underlying differences in biology associated with different forms of PCOS. We provide a state-of-the-science review of genetic analyses of PCOS, including an overview of genomic methodologies aimed at a general audience of non-geneticists and clinicians. Applications in PCOS will be discussed, including strengths and limitations of each study. The contributions of environmental factors, including developmental origins, will be reviewed. Insights into the pathogenesis and genetic architecture of PCOS will be summarized. Future directions for PCOS genetic studies will be outlined.

Strategies for sample labelling and library preparation in DNA metabarcoding studies

Metabarcoding of DNA extracted from environmental or bulk specimen samples is increasingly used to profile biota in basic and applied biodiversity research because of its targeted nature that allows sequencing of genetic markers from many samples in parallel. To achieve this, PCR amplification is carried out with primers designed to target a taxonomically informative marker within a taxonomic group, and sample-specific nucleotide identifiers are added to the amplicons prior to sequencing. The latter enables assignment of the sequences back to the samples they originated from. Nucleotide identifiers can be added during the metabarcoding PCR and during "library preparation", that is, when amplicons are prepared for sequencing. Different strategies to achieve this labelling exist. All have advantages, challenges and limitations, some of which can lead to misleading results, and in the worst case compromise the fidelity of the metabarcoding data. Given the range of questions addressed using metabarcoding, ensuring that data generation is robust and fit for the chosen purpose is critically important for practitioners seeking to employ metabarcoding for biodiversity assessments. Here, we present an overview of the three main workflows for sample-specific labelling and library preparation in metabarcoding studies on Illumina sequencing platforms; one-step PCR, two-step PCR, and tagged PCR. Further, we distill the key considerations for researchers seeking to select an appropriate metabarcoding strategy for their specific study. Ultimately, by gaining insights into the consequences of different metabarcoding workflows, we hope to further consolidate the power of metabarcoding as a tool to assess biodiversity across a range of applications.

Gut microbiota differences between paired intestinal wall and digesta samples in three small species of fish

The microbial gut communities of fish are receiving increased attention for their relevance, among others, in a growing aquaculture industry. The members of these communities are often split into resident (long-term colonisers specialised to grow in and adhere to the mucus lining of the gut) and transient (short-term colonisers originated from food items and the surrounding water) microorganisms. Separating these two communities in small fish are impeded by the small size and fragility of the gastrointestinal tract. With the aim of testing whether it is possible to recover two distinct communities in small species of fish using a simple sampling technique, we used 16S amplicon sequencing of paired intestinal wall and digesta samples from three small Cyprinodontiformes fish. We examined the diversity and compositional variation of the two recovered communities, and we used joint species distribution modelling to identify microbes that are most likely to be a part of the resident community. For all three species we found that the diversity of intestinal wall samples was significantly lower compared to digesta samples and that the community composition between sample types was significantly different. Across the three species we found seven unique families of bacteria to be significantly enriched in samples from the intestinal wall, encompassing most of the 89 ASVs enriched in intestinal wall samples. We conclude that it is possible to characterise two different microbial communities and identify potentially resident microbes through separately analysing samples from the intestinal wall and digesta from small species of fish. We encourage researchers to be aware that different sampling procedures for gut microbiome characterization will capture different parts of the microbiome and that this should be taken into consideration when reporting results from such studies on small species of fish.

Airborne environmental DNA for terrestrial vertebrate community monitoring

Biodiversity monitoring at the community scale is a critical element of assessing and studying species distributions, ecology, diversity, and movements, and it is key to understanding and tracking environmental and anthropogenic effects on natural ecosystems.1, 2, 3, 4 Vertebrates in terrestrial ecosystems are experiencing extinctions and declines in both population numbers and sizes due to increasing threats from human activities and environmental change.5, 6, 7, 8 Terrestrial vertebrate monitoring using existing methods is generally costly and laborious, and although environmental DNA (eDNA) is becoming the tool of choice to assess biodiversity, few sample types effectively capture terrestrial vertebrate diversity. We hypothesized that eDNA captured from air could allow straightforward collection and characterization of terrestrial vertebrate communities. We filtered air at three localities in the Copenhagen Zoo: a stable, outside between the outdoor enclosures, and in the Rainforest House. Through metabarcoding of airborne eDNA, we detected 49 vertebrate species spanning 26 orders and 37 families: 30 mammal, 13 bird, 4 fish, 1 amphibian, and 1 reptile species. These spanned animals kept at the zoo, species occurring in the zoo surroundings, and species used as feed in the zoo. The detected species comprise a range of taxonomic orders and families, sizes, behaviors, and abundances. We found shorter distance to the air sampling device and higher animal biomass to increase the probability of detection. We hereby show that airborne eDNA can offer a fundamentally new way of studying and monitoring terrestrial communities.

•

49 vertebrate species detected through metabarcoding of airborne eDNA from the zoo

•

Detections included 30 mammal, 13 bird, 4 fish, 1 amphibian, and 1 reptile species

•

6 to 21 vertebrate species were detected per air filtering sample

•

Shorter geographical distance and higher biomass increased probability of detection

The potential of aquatic bloodfeeding and nonbloodfeeding leeches as a tool for iDNA characterisation

Leeches play important roles in food webs due to their abundance, diversity and feeding habits. Studies using invertebrate-derived DNA (iDNA) extracted from leech gut contents to target vertebrate DNA have focused on the Indo-Pacific region and mainly leveraged the leech family Haemadipsidae, composed of bloodfeeding terrestrial leeches, while predatory, fluid/tissue-feeding and aquatic bloodfeeding species have been largely disregarded. While there is some general knowledge regarding the taxonomic groups that leeches prefer to feed on, detailed taxonomic resolution is missing and, therefore, their potential use for monitoring animals is unknown. In this study, 116 leeches from 12 species (six families) and spanning the three feeding habits were collected in Mexico and Canada. We used DNA metabarcoding to investigate their diet and assess their potential use for biodiversity monitoring. We detected vertebrates from five orders including fish, turtles and birds in the diet of aquatic bloodfeeding leeches; eight invertebrate orders of annelids, arthropods and molluscs in leeches that feed on body fluids and tissues; and 10 orders of invertebrates belonging to Arthropoda and Annelida, as well as one vertebrate and one parasitic nematode, in predatory leeches. These results show the potential use of iDNA from aquatic bloodfeeding leeches for retrieving vertebrate taxa, and from predatory and fluid-feeding leeches for invertebrates. Our study provides information about the dietary range of freshwater leeches and one terrestrial leech and contributes proof-of-concept for the use of these leeches for animal monitoring, expanding our knowledge of the use of iDNA from leech gut contents to North America.

Biological composition analysis of a natural medicine, Faeces Vespertilionis, with complex sources using DNA metabarcoding

Faeces Vespertilionis is a commonly used fecal traditional Chinese medicine. Traditionally, it is identified relying only on morphological characters. This poses a serious challenge to the composition analysis accuracy of this complex biological mixture. Thus, for quality control purposes, an accurate and effective method should be provided for taxonomic identification of Faeces Vespertilionis. In this study, 26 samples of Faeces Vespertilionis from ten provinces in China were tested using DNA metabarcoding. Seven operational taxonomic units (OTUs) were detected as belonging to bats. Among them, Hipposideros armiger (Hodgson, 1835) and Rhinolophus ferrumequinum (Schober and Grimmberger, 1997) were the main host sources of Faeces Vespertilionis samples, with average relative abundances of 59.3% and 24.1%, respectively. Biodiversity analysis showed that Diptera and Lepidoptera were the most frequently consumed insects. At the species level, 19 taxa were clearly identified. Overall, our study used DNA metabarcoding to analyze the biological composition of Faeces Vespertilionis, which provides a new idea for the quality control of this special traditional Chinese medicine.

Biodiversity and vector-borne diseases: host dilution and vector amplification occur simultaneously for Amazonian leishmaniases

Changes in biodiversity may impact infectious disease transmission through multiple mechanisms. We explored the impact of biodiversity changes on the transmission of Amazonian leishmaniases, a group of wild zoonoses transmitted by phlebotomine sand flies (Psychodidae), which represent an important health burden in a region where biodiversity is both rich and threatened. Using molecular analyses of sand fly pools and blood-fed dipterans, we characterized the disease system in forest sites in French Guiana undergoing different levels of human-induced disturbance. We show that the prevalence of Leishmania parasites in sand flies correlates positively with the relative abundance of mammal species known as Leishmania reservoirs. In addition, Leishmania reservoirs tend to dominate in less diverse mammal communities, in accordance with the dilution effect hypothesis. This results in a negative relationship between Leishmania prevalence and mammal diversity. On the other hand, higher mammal diversity is associated with higher sand fly density, possibly because more diverse mammal communities harbor higher biomass and more abundant feeding resources for sand flies, although more research is needed to identify the factors that shape sand fly communities. As a consequence of these antagonistic effects, decreased mammal diversity comes with an increase of parasite prevalence in sand flies, but has no detectable impact on the density of infected sand flies. These results represent additional evidence that biodiversity changes may simultaneously dilute and amplify vector-borne disease transmission through different mechanisms that need to be better understood before drawing generalities on the biodiversity-disease relationship.

DNA metabarcoding reveals introduced species predominate in the diet of a threatened endemic omnivore, Telfair's skink (Leiolopisma telfairii)

Introduced species can exert disproportionately negative effects on island ecosystems, but their potential role as food for native consumers is poorly studied. Telfair's skinks are endemic omnivores living on Round Island, Mauritius, a globally significant site of biodiversity conservation. We aimed to determine the dietary diversity and key trophic interactions of Telfair's skinks, whether introduced species are frequently consumed, and if diet composition changes seasonally between male and female skinks. We used DNA metabarcoding of skink fecal samples to identify animals (COI) and plants (ITS2) consumed by skinks. There were 389 dietary presence counts belonging to 77 dietary taxa found across the 73 Telfair's skink fecal samples. Introduced taxa were cumulatively consumed more frequently than other categories, accounting for 49.4% of all detections, compared to cryptogenic (20.6%), native (20.6%), and endemic taxa (9.5%). The most frequently consumed introduced species was the ant, Pheidole megacephala, present in 40% of samples. Blue latan palm, Latania loddigesii, was the most frequently consumed endemic species, present in 33% of samples but was only detected in the dry season, when fruits are produced. We found a strong seasonal difference in diet composition explained by the presence of certain plant species solely or primarily in one season and a marked increase in the consumption of animal prey in the dry season. Male and female skinks consumed several taxa at different frequencies. These results present a valuable perspective on the role of introduced species in the trophic network of their invaded ecosystem. Both native and introduced species provide nutritional resources for skinks, and this may have management implications in the context of species conservation and island restoration.

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.

Plant Sedimentary Ancient DNA From Far East Russia Covering the Last 28,000 Years Reveals Different Assembly Rules in Cold and Warm Climates

Woody plants are expanding into the Arctic in response to the warming climate. The impact on arctic plant communities is not well understood due to the limited knowledge about plant assembly rules. Records of past plant diversity over long time series are rare. Here, we applied sedimentary ancient DNA metabarcoding targeting the P6 loop of the chloroplast trnL gene to a sediment record from Lake Ilirney (central Chukotka, Far Eastern Russia) covering the last 28 thousand years. Our results show that forb-rich steppe-tundra and dwarf-shrub tundra dominated during the cold climate before 14 ka, while deciduous erect-shrub tundra was abundant during the warm period since 14 ka. Larix invasion during the late Holocene substantially lagged behind the likely warmest period between 10 and 6 ka, where the vegetation biomass could be highest. We reveal highest richness during 28–23 ka and a second richness peak during 13–9 ka, with both periods being accompanied by low relative abundance of shrubs. During the cold period before 14 ka, rich plant assemblages were phylogenetically clustered, suggesting low genetic divergence in the assemblages despite the great number of species. This probably originates from environmental filtering along with niche differentiation due to limited resources under harsh environmental conditions. In contrast, during the warmer period after 14 ka, rich plant assemblages were phylogenetically overdispersed. This results from a high number of species which were found to harbor high genetic divergence, likely originating from an erratic recruitment process in the course of warming. Some of our evidence may be of relevance for inferring future arctic plant assembly rules and diversity changes. By analogy to the past, we expect a lagged response of tree invasion. Plant richness might overshoot in the short term; in the long-term, however, the ongoing expansion of deciduous shrubs will eventually result in a phylogenetically more diverse community.

The role of the gut microbiota in the dietary niche expansion of fishing bats

Background

Due to its central role in animal nutrition, the gut microbiota is likely a relevant factor shaping dietary niche shifts. We analysed both the impact and contribution of the gut microbiota to the dietary niche expansion of the only four bat species that have incorporated fish into their primarily arthropodophage diet.

Results

We first compared the taxonomic and functional features of the gut microbiota of the four piscivorous bats to that of 11 strictly arthropodophagous species using 16S rRNA targeted amplicon sequencing. Second, we increased the resolution of our analyses for one of the piscivorous bat species, namely Myotis capaccinii, and analysed multiple populations combining targeted approaches with shotgun sequencing. To better understand the origin of gut microorganisms, we also analysed the gut microbiota of their fish prey (Gambusia holbrooki). Our analyses showed that piscivorous bats carry a characteristic gut microbiota that differs from that of their strict arthropodophagous counterparts, in which the most relevant bacteria have been directly acquired from their fish prey. This characteristic microbiota exhibits enrichment of genes involved in vitamin biosynthesis, as well as complex carbohydrate and lipid metabolism, likely providing their hosts with an enhanced capacity to metabolise the glycosphingolipids and long-chain fatty acids that are particularly abundant in fish.

Conclusions

Our results depict the gut microbiota as a relevant element in facilitating the dietary transition from arthropodophagy to piscivory.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00137-w.

Dung beetles as samplers of mammals in Malaysian Borneo-a test of high throughput metabarcoding of iDNA

Invertebrate-derived DNA (iDNA) sampling in biodiversity surveys is becoming increasingly widespread, with most terrestrial studies relying on DNA derived from the gut contents of blood-feeding invertebrates, such as leeches and mosquitoes. Dung beetles (superfamily Scarabaeoidea) primarily feed on the faecal matter of terrestrial vertebrates and offer several potential benefits over blood-feeding invertebrates as samplers of vertebrate DNA. Importantly, these beetles can be easily captured in large numbers using simple, inexpensive baited traps, are globally distributed, and occur in a wide range of habitats. To build on the few existing studies demonstrating the potential of dung beetles as sources of mammalian DNA, we subjected the large-bodied, Bornean dung beetle (Catharsius renaudpauliani) to a controlled feeding experiment. We analysed DNA from gut contents at different times after feeding using qPCR techniques. Here, we first describe the window of DNA persistence within a dung beetle digestive tract. We found that the ability to successfully amplify cattle DNA decayed over relatively short time periods, with DNA copy number decreasing by two orders of magnitude in just 6 h. In addition, we sampled communities of dung beetles from a lowland tropical rainforest in Sabah, Malaysia, in order to test whether it is possible to identify vertebrate sequences from dung beetle iDNA. We sequenced both the gut contents from large dung beetle species, as well as whole communities of smaller beetles. We successfully identified six mammalian species from our samples, including the bearded pig (Sus barbatus) and the sambar deer (Rusa unicolor)—both vulnerable species on the IUCN red list. Our results represent the first use of dung beetle iDNA to sample Southeast Asian vertebrate fauna, and highlight the potential for dung beetle iDNA to be used in future biodiversity monitoring surveys.

Towards reproducible metabarcoding data: Lessons from an international cross-laboratory experiment

Advances in high-throughput sequencing (HTS) are revolutionizing monitoring in marine environments by enabling rapid, accurate and holistic detection of species within complex biological samples. Research institutions worldwide increasingly employ HTS methods for biodiversity assessments. However, variance in laboratory procedures, analytical workflows and bioinformatic pipelines impede the transferability and comparability of results across research groups. An international experiment was conducted to assess the consistency of metabarcoding results derived from identical samples and primer sets using varying laboratory procedures. Homogenized biofouling samples collected from four coastal locations (Australia, Canada, New Zealand and the USA) were distributed to 12 independent laboratories. Participants were asked to follow one of two HTS library preparation workflows. While DNA extraction, primers and bioinformatic analyses were purposefully standardized to allow comparison, many other technical variables were allowed to vary among laboratories (amplification protocols, type of instrument used, etc.). Despite substantial variation observed in raw results, the primary signal in the data was consistent, with the samples grouping strongly by geographical origin for all data sets. Simple post hoc data clean-up by removing low-quality samples gave the best improvement in sample classification for nuclear 18S rRNA gene data, with an overall 92.81% correct group attribution. For mitochondrial COI gene data, the best classification result (95.58%) was achieved after correction for contamination errors. The identified critical methodological factors that introduced the greatest variability (preservation buffer, sample defrosting, template concentration, DNA polymerase, PCR enhancer) should be of great assistance in standardizing future biodiversity studies using metabarcoding.

Cascading effects of moth outbreaks on subarctic soil food webs

The increasing severity and frequency of natural disturbances requires a better understanding of their effects on all compartments of biodiversity. In Northern Fennoscandia, recent large-scale moth outbreaks have led to an abrupt change in plant communities from birch forests dominated by dwarf shrubs to grass-dominated systems. However, the indirect effects on the belowground compartment remained unclear. Here, we combined eDNA surveys of multiple trophic groups with network analyses to demonstrate that moth defoliation has far-reaching consequences on soil food webs. Following this disturbance, diversity and relative abundance of certain trophic groups declined (e.g., ectomycorrhizal fungi), while many others expanded (e.g., bacterivores and omnivores) making soil food webs more diverse and structurally different. Overall, the direct and indirect consequences of moth outbreaks increased belowground diversity at different trophic levels. Our results highlight that a holistic view of ecosystems improves our understanding of cascading effects of major disturbances on soil food webs.

Structure and diversity of bacterial community in semiarid soils cultivated with prickly-pear cactus (Opuntia ficus-indica (L.) Mill.)

Studies evaluating the structure and diversity of bacterial communities in arid environments including the rhizosphere of local and adapted plant species are important. Therefore, we used a sequencing of the 16S ribosomal RNA gene for describing the structure and diversity of soil bacterial community in three zones: Agreste, Transition and Sertão. The bacterial community was clustered in 9,838 OTUs in Agreste, 8,388 OTUs in the transition, and 14,849 OTUs for Sertão. Among the most abundant phyla, Proteobacteria and Acidobacteria were abundant in Agreste and Sertão, respectively, while Actinobacteria were abundant in Transition and Sertão. Specific taxa of Proteobacteria, in Agreste, and Actinobacteria, in Sertão, exhibited differences according to biotic and abiotic conditions. Thus, the structure and diversity of bacterial community were different in these areas and were influenced by environmental and soil conditions.

Microbial community diversity and composition in river sediments contaminated with tetrabromobisphenol A and copper

The microbial community composition in aquatic ecosystems have received increased attention. However, the knowledge gap relative to the responses of bacterial, archaeal and fungal communities in co-contaminated river sediments remain poorly studied. Here, we investigated the changes of tetrabromobisphenol A (TBBPA) and copper (Cu) concentrations and the responses of microbial communities in co-contaminated sediments during long-term incubation. TBBPA concentrations significantly decreased over time, whereas Cu concentrations remained relatively stable over the 60-day incubation. Abundances of the bacterial 16S rRNA, archaeal 16S rRNA and fungal ITS genes ranged from 6.53 × 10⁶ to 1.26 × 10⁹ copies g^-1, 1.12 × 10⁶ to 5.47 × 10⁶ copies g^-1 and 5.33 × 10³ to 7.51 × 10⁴ copies g^-1 in the samples, respectively. A total of 11, 6 and 5 bacterial, archaeal and fungal phyla were identified across all samples. Bacterial, archaeal and fungal communities mainly consisted of members from the phyla Proteobacteria and Acidobacteria, Methanomicrobia and Woesearchaeia as well as Agaricales and Helotiales, respectively. Fungal communities showed a stronger response to pollutant addition after a long incubation compared with bacterial and archaeal communities. The variance analysis results revealed that the bacterial, archaeal and fungal microbial communities of all treatments were distinctly distributed into two separated clusters according to incubation time. However, the three microbial communities did not significantly change in response to pollutant types, which was consistent with variation in relative abundances of the three microbial communities. These findings improve our understanding of the ecotoxicological effects of co-exposure on sediment microbial communities.

Metagenomics: A viable tool for reconstructing herbivore diet

Metagenomics can generate data on the diet of herbivores, without the need for primer selection and PCR enrichment steps as is necessary in metabarcoding. Metagenomic approaches to diet analysis have remained relatively unexplored, requiring validation of bioinformatic steps. Currently, no metagenomic herbivore diet studies have utilized both chloroplast and nuclear markers as reference sequences for plant identification, which would increase the number of reads that could be taxonomically informative. Here, we explore how in silico simulation of metagenomic data sets resembling sequences obtained from faecal samples can be used to validate taxonomic assignment. Using a known list of sequences to create simulated data sets, we derived reliable identification parameters for taxonomic assignments of sequences. We applied these parameters to characterize the diet of western capercaillies (Tetrao urogallus) located in Norway, and compared the results with metabarcoding trnL P6 loop data generated from the same samples. Both methods performed similarly in the number of plant taxa identified (metagenomics 42 taxa, metabarcoding 43 taxa), with no significant difference in species resolution (metagenomics 24%, metabarcoding 23%). We further observed that while metagenomics was strongly affected by the age of faecal samples, with fresh samples outperforming old samples, metabarcoding was not affected by sample age. On the other hand, metagenomics allowed us to simultaneously obtain the mitochondrial genome of the western capercaillies, thereby providing additional ecological information. Our study demonstrates the potential of utilizing metagenomics for diet reconstruction but also highlights key considerations as compared to metabarcoding for future utilization of this technique.

Evaluation of V2 18S rDNA barcode marker and assessment of sample collection and DNA extraction methods for metabarcoding of autotrophic euglenids

Even though the interest in metabarcoding in environmental studies is growing, euglenids are still underrepresented in both sea and freshwater bodies researches. The reason for this situation could be the unsuitability of universal eukaryotic DNA barcodes and primers as well as the lack of a verified protocol, suitable to assess euglenid diversity. In this study, using specific primers for the V2 hypervariable region of 18S rDNA for metabarcoding resulted in obtaining a high fraction (85%) of euglenid reads and species‐level identification of almost 90% of them. Fifty species were detected by the metabarcoding method, including almost all species observed using a light microscope. We investigated three biomass harvesting methods (filtering, centrifugation and scraping the side of a collection vessel) and determined that centrifugation and filtration outperformed scrapes, but the choice between them is not crucial for the reliability of the analysis. In addition, eight DNA extraction methods were evaluated. We compared five commercially available DNA isolation kits, two CTAB‐based protocols and a chelating resin. For this purpose, the efficiency of extraction, quality of obtained DNA, preparation time and generated costs were taken into consideration. After examination of the aforementioned criteria, we chose the GeneMATRIX Soil DNA Purification Kit as the most suitable for DNA isolation.

Vegetation Changes in Southeastern Siberia During the Late Pleistocene and the Holocene

Relationships between climate, species composition, and species richness are of particular importance for understanding how boreal ecosystems will respond to ongoing climate change. This study aims to reconstruct changes in terrestrial vegetation composition and taxa richness during the glacial Late Pleistocene and the interglacial Holocene in the sparsely studied southeastern Yakutia (Siberia) by using pollen and sedimentary ancient DNA (sedaDNA) records. Pollen and sedaDNA metabarcoding data using the trnL g and h markers were obtained from a sediment core from Lake Bolshoe Toko. Both proxies were used to reconstruct the vegetation composition, while metabarcoding data were also used to investigate changes in plant taxa richness. The combination of pollen and sedaDNA approaches allows a robust estimation of regional and local past terrestrial vegetation composition around Bolshoe Toko during the last ∼35,000 years. Both proxies suggest that during the Late Pleistocene, southeastern Siberia was covered by open steppe-tundra dominated by graminoids and forbs with patches of shrubs, confirming that steppe-tundra extended far south in Siberia. Both proxies show disturbance at the transition between the Late Pleistocene and the Holocene suggesting a period with scarce vegetation, changes in the hydrochemical conditions in the lake, and in sedimentation rates. Both proxies document drastic changes in vegetation composition in the early Holocene with an increased number of trees and shrubs and the appearance of new tree taxa in the lake’s vicinity. The sedaDNA method suggests that the Late Pleistocene steppe-tundra vegetation supported a higher number of terrestrial plant taxa than the forested Holocene. This could be explained, for example, by the “keystone herbivore” hypothesis, which suggests that Late Pleistocene megaherbivores were able to maintain a high plant diversity. This is discussed in the light of the data with the broadly accepted species-area hypothesis as steppe-tundra covered such an extensive area during the Late Pleistocene.

Microbial community compositions in breast implant biofilms associated with contracted capsules

Subclinical bacterial infections (biofilms) are strongly implicated in breast augmentation failure due to capsular contracture, and while these infections are generally ascribed to common skin commensals, this remains largely unsubstantiated through robust cultivation independent analyses. To determine capsule biofilm microbial community compositions, we employed amplicon sequencing of the 16S rRNA gene using DNA extracted from breast implant capsule samples. These cultivation independent analyses revealed that capsule associated biofilms are more diverse than canonical single-species infections, but have relatively low diversity (~ <100 species) compared to many host-associated microbial communities. In addition to taxa commonly associated with capsular contracture, the biofilms analyzed comprised a number of taxa that escaped detection in cultivation-dependent work. We have also isolated several key taxa identified through the culture-independent analyses. Together our analyses reveal that capsule biofilms are more diverse than cultivation studies suggest and can be heterogeneous within an individual capsule, between breasts of the same patient, across similar implant types, and over a range in severity of contracture. The complex nature of these communities requires further study across a broader suite of patients in addition to higher resolution analyses including metagenomics to better assess the fundamental role of microorganisms in capsular contracture.

A library preparation optimized for metagenomics of RNA viruses

Our understanding of the viral communities associated to animals has not yet reached the level attained on the bacteriome. This situation is due to, among others, technical challenges in adapting metagenomics using high-throughput sequencing to the study of RNA viromes in animals. Although important developments have been achieved in most steps of viral metagenomics, there is yet a key step that has received little attention: the library preparation. This situation differs from bacteriome studies in which developments in library preparation have largely contributed to the democratisation of metagenomics. Here, we present a library preparation optimized for metagenomics of RNA viruses from insect vectors of viral diseases. The library design allows a simple PCR-based preparation, such as those routinely used in bacterial metabarcoding, that is adapted to shotgun sequencing as required in viral metagenomics. We first optimized our library preparation using mock viral communities and then validated a full metagenomic approach incorporating our preparation in two pilot studies with field-caught insect vectors; one including a comparison with a published metagenomic protocol. Our approach provided a fold increase in virus-like sequences compared to other studies, and nearly-full genomes from new virus species. Moreover, our results suggested conserved trends in virome composition within a population of a mosquito species. Finally, the sensitivity of our approach was compared to a commercial diagnostic PCR for the detection of an arbovirus in field-caught insect vectors. Our approach could facilitate studies on viral communities from animals and the democratization of metagenomics in community ecology of viruses.

Targeted Sequencing of Mitochondrial Genes Reveals Signatures of Molecular Adaptation in a Nearly Panmictic Small Pelagic Fish Species

Ongoing climatic changes, with predictable impacts on marine environmental conditions, are expected to trigger organismal responses. Recent evidence shows that, in some marine species, variation in mitochondrial genes involved in the aerobic conversion of oxygen into ATP at the cellular level correlate with gradients of sea surface temperature and gradients of dissolved oxygen. Here, we investigated the adaptive potential of the European sardine Sardina pilchardus populations offshore the Iberian Peninsula. We performed a seascape genetics approach that consisted of the high throughput sequencing of mitochondria’s ATP6, COI, CYTB and ND5 and five microsatellite loci on 96 individuals coupled with environmental information on sea surface temperature and dissolved oxygen across five sampling locations. Results show that, despite sardines forming a nearly panmictic population around Iberian Peninsula, haplotype frequency distribution can be explained by gradients of minimum sea surface temperature and dissolved oxygen. We further identified that the frequencies of the most common CYTB and ATP6 haplotypes negatively correlate with minimum sea surface temperature across the sampled area, suggestive of a signature of selection. With signatures of selection superimposed on highly connected populations, sardines may be able to follow environmental optima and shift their distribution northwards as a response to the increasing sea surface temperatures.

DNA-Based Arthropod Diversity Assessment in Amazonian Iron Mine Lands Show Ecological Succession Towards Undisturbed Reference Sites

Human activities change natural landscapes, and in doing so endanger biodiversity and associated ecosystem services. To reduce the net impacts of these activities, such as mining, disturbed areas are rehabilitated and restored. During this process, monitoring is important to ensure that desired trajectories are maintained. In the Carajás region of the Brazilian Amazon, exploration for iron ores has transformed the original ecosystem; natural forest and a savanna formation with lateritic iron duricrust outcrops named canga. Here, native vegetation is logged and topsoil removed and deposited in waste piles along with mine waste. During rehabilitation, these waste piles are hydroseeded with non-native plant species to achieve rapid revegetation. Further, seeds of native canga and forest plant species are planted to point ecological succession towards natural ecosystems. In this study, we investigate diversity and composition of the arthropod community along a post-mining rehabilitation and restoration gradient, taking seasonality and primer bias into account. We use DNA metabarcoding of bulk arthropod samples collected in both the dry and rainy seasons from waste-pile benches at various stages of revegetation: non-revegetated exposed soils, initial stage with one-to-three-year-old stands, intermediate stage with four-to-five-year-old stands, and advanced stage with six-to-seven-year-old stands. We use samples from undisturbed cangas and forests as reference sites. In addition, we vegetation diversity and structure were measured to investigate relations between arthropod community and vegetation structure. Our results show that, over time, the arthropod community composition of the waste piles becomes more similar to the reference forests, but not to the reference cangas. Nevertheless, even the communities in the advanced-stage waste piles are different from the reference forests, and full restoration in these highly diverse ecosystems is not achieved, even after 6 to 7 years. Finally, our results show seasonal variation in arthropod communities and primer bias.

Leech blood-meal invertebrate-derived DNA reveals differences in Bornean mammal diversity across habitats

The application of metabarcoding to environmental and invertebrate‐derived DNA (eDNA and iDNA) is a new and increasingly applied method for monitoring biodiversity across a diverse range of habitats. This approach is particularly promising for sampling in the biodiverse humid tropics, where rapid land‐use change for agriculture means there is a growing need to understand the conservation value of the remaining mosaic and degraded landscapes. Here we use iDNA from blood‐feeding leeches (Haemadipsa picta) to assess differences in mammalian diversity across a gradient of forest degradation in Sabah, Malaysian Borneo. We screened 557 individual leeches for mammal DNA by targeting fragments of the 16S rRNA gene and detected 14 mammalian genera. We recorded lower mammal diversity in the most heavily degraded forest compared to higher quality twice logged forest. Although the accumulation curves of diversity estimates were comparable across these habitat types, diversity was higher in twice logged forest, with more taxa of conservation concern. In addition, our analysis revealed differences between the community recorded in the heavily logged forest and that of the twice logged forest. By revealing differences in mammal diversity across a human‐modified tropical landscape, our study demonstrates the value of iDNA as a noninvasive biomonitoring approach in conservation assessments.

Tagsteady: A metabarcoding library preparation protocol to avoid false assignment of sequences to samples

Metabarcoding of environmental DNA (eDNA) and DNA extracted from bulk specimen samples is a powerful tool in studies of biodiversity, diet and ecological interactions as its inherent labelling of amplicons allows sequencing of taxonomically informative genetic markers from many samples in parallel. However, the occurrence of so-called 'tag-jumps' can cause incorrect assignment of sequences to samples and artificially inflate diversity. Two steps during library preparation of pools of 5' nucleotide-tagged amplicons have been suggested to cause tag-jumps: (a) T4 DNA polymerase blunt-ending in the end-repair step and (b) postligation PCR amplification of amplicon libraries. The discovery of tag-jumps has led to recommendations to only carry out metabarcoding PCR amplifications with primers carrying twin-tags to ensure that tag-jumps cannot result in false assignments of sequences to samples. As this increases both cost and workload, a metabarcoding library preparation protocol which circumvents the two steps that causes tag-jumps is needed. Here, we demonstrate Tagsteady, a PCR-free metabarcoding Illumina library preparation protocol for pools of nucleotide-tagged amplicons that enables efficient and cost-effective generation of metabarcoding data with virtually no tag-jumps. We use pools of twin-tagged amplicons to investigate the effect of T4 DNA polymerase blunt-ending and postligation PCR on the occurrence of tag-jumps and demonstrate that both blunt-ending and postligation PCR, alone or together, can result in detrimental amounts of tag-jumps (here, up to ca. 49% of total sequences), while leaving both steps out (the Tagsteady protocol) results in amounts of sequences carrying new combinations of used tags (tag-jumps) comparable to background contamination.

An appetite for pests: Synanthropic insectivorous bats exploit cotton pest irruptions and consume various deleterious arthropods

Conservation biological control (CBC) seeks to minimize the deleterious effects of agricultural pests by enhancing the efficiency of natural enemies. Despite the documented potential of insectivorous bats to consume pests, many synanthropic bat species are still underappreciated as beneficial species. We investigated the diet of Kuhl's pipistrelle (Pipistrellus kuhlii), a common synanthropic insectivorous bat that forages in urban and agricultural areas, to determine whether it may function as a natural enemy in CBC. Faecal samples of P. kuhlii were collected throughout the cotton-growing season from five roost sites near cotton fields located in a Mediterranean agroecosystem, Israel, and analyzed using DNA metabarcoding. Additionally, data on estimated abundance of major cotton pests were collected. We found that the diet of P. kuhlii significantly varied according to sites and dates and comprised 27 species of agricultural pests that were found in 77.2% of the samples, including pests of key economic concern. The dominant prey was the widespread cotton pest, the pink bollworm, Pectinophora gossypiella, found in 31% of the samples and in all the roosts. Pink bollworm abundance was positively correlated with its occurrence in the bat diet. Furthermore, the bats' dietary breadth narrowed, while temporal dietary overlap increased, in relation to increasing frequencies of pink bollworms in the diet. This suggests that P. kuhlii exploits pink bollworm irruptions by opportunistic feeding. We suggest that synanthropic bats provide important pest suppression services, may function as CBC agents of cotton pests and potentially contribute to suppress additional deleterious arthropods found in their diet in high frequencies.

Changes in Bacterial Communities in Seawater-Flooded Soil in the Four Years After the 2011 Tohoku Tsunami in Japan

The 2011 Tohoku tsunami had a serious impact, such as an increase in harmful substances and salinity over a large area. Herein, we evaluated transitions in bacterial communities in agricultural fields in the four years after the 2011 Tohoku tsunami. Bacterial communities were compared across four different types of soil—unflooded field (UF) soil, soil flooded for a short term (ST), soil flooded for the long term (LT), soil flooded long term and cultivated fields (LTC), and marine environmental materials (bay sediment, sea sand and sea water), using a polymerase chain reaction (PCR) and pyrosequencing of 16S ribosomal RNA genes. In the soil bacterial communities that were flooded by the 2011 Tohoku tsunami, these effects were not seen after 2013. Although the difference in bacterial communities between LT and UF became smaller during the four years, the bacterial communities in LT were different from those in UF in several ways, such as a higher tendency frequency of sulfur-oxidizing bacteria (SOB) and the presence of halotolerant SOB. Therefore, it is thought that the Tohoku tsunami affected the microbial communities in the soil for more than four years. Especially genus Halothiobacillus, which is Halotolerant SOB in flooded soils, was detected neither in unflooded soil nor in the marine environment. Therefore, it is thought that inundation by a tsunami produces a unique environment with bacterial communities to form in soil. Further, SOB structure, especially halotolerant, might serve as a good indicators of the impacts of inundation on bacterial communities in agricultural fields over the long term.

Trophic rewilding revives biotic resistance to shrub invasion

Trophic rewilding seeks to rehabilitate degraded ecosystems by repopulating them with large animals, thereby re-establishing strong top-down interactions. Yet there are very few tests of whether such initiatives can restore ecosystem structure and functions, and on what timescales. Here we show that war-induced collapse of large-mammal populations in Mozambique's Gorongosa National Park exacerbated woody encroachment by the invasive shrub Mimosa pigra-considered one of the world's 100 worst invasive species-and that one decade of concerted trophic rewilding restored this invasion to pre-war baseline levels. Mimosa occurrence increased between 1972 and 2015, a period encompassing the near extirpation of large herbivores during the Mozambican Civil War. From 2015 to 2019, mimosa abundance declined as ungulate biomass recovered. DNA metabarcoding revealed that ruminant herbivores fed heavily on mimosa, and experimental exclosures confirmed the causal role of mammalian herbivory in containing shrub encroachment. Our results provide mechanistic evidence that trophic rewilding has rapidly revived a key ecosystem function (biotic resistance to a notorious woody invader), underscoring the potential for restoring ecological health in degraded protected areas.

The Soil Microbiome of the Laurel Forest in Garajonay National Park (La Gomera, Canary Islands): Comparing Unburned and Burned Habitats after a Wildfire

The evergreen laurel forest is a relic of ancient subtropical/tropical forests, of which the best remnant in the Canary Islands is in Garajonay National Park, on La Gomera island. The soil microbiome associated with a mature undisturbed (unburned) laurel forest was characterized at two locations at different topographical positions on the mountain: The slope and the ridge crest. Given the unusual circumstance of an intense wildfire that severely affected part of this forest, the burned soils were also studied. The soil in undisturbed areas was relatively uniform. The bacterial community composition was dominated by bacteria from phyla Proteobacteria, Acidobacteria, and Actinobacteria. The wildfire changed the composition of the bacterial communities. The Acidobacteria, Actinobacteria, and Alphaproteobacteria (dominant class in unburned forests) significantly decreased in burned soils along with a parallel high increase in Betaproteobacteria, Bacteroidetes, and Firmicutes. We further showed the dramatic effect of a wildfire on the soil microbiome of the laurel forest, appearing as a loss of species richness and diversity, species dominance, and changes in the composition of the bacterial communities.

Adapterama I: universal stubs and primers for 384 unique dual-indexed or 147,456 combinatorially-indexed Illumina libraries (iTru & iNext)

Massively parallel DNA sequencing offers many benefits, but major inhibitory cost factors include: (1) start-up (i.e., purchasing initial reagents and equipment); (2) buy-in (i.e., getting the smallest possible amount of data from a run); and (3) sample preparation. Reducing sample preparation costs is commonly addressed, but start-up and buy-in costs are rarely addressed. We present dual-indexing systems to address all three of these issues. By breaking the library construction process into universal, re-usable, combinatorial components, we reduce all costs, while increasing the number of samples and the variety of library types that can be combined within runs. We accomplish this by extending the Illumina TruSeq dual-indexing approach to 768 (384 + 384) indexed primers that produce 384 unique dual-indexes or 147,456 (384 × 384) unique combinations. We maintain eight nucleotide indexes, with many that are compatible with Illumina index sequences. We synthesized these indexing primers, purifying them with only standard desalting and placing small aliquots in replicate plates. In qPCR validation tests, 206 of 208 primers tested passed (99% success). We then created hundreds of libraries in various scenarios. Our approach reduces start-up and per-sample costs by requiring only one universal adapter that works with indexed PCR primers to uniquely identify samples. Our approach reduces buy-in costs because: (1) relatively few oligonucleotides are needed to produce a large number of indexed libraries; and (2) the large number of possible primers allows researchers to use unique primer sets for different projects, which facilitates pooling of samples during sequencing. Our libraries make use of standard Illumina sequencing primers and index sequence length and are demultiplexed with standard Illumina software, thereby minimizing customization headaches. In subsequent Adapterama papers, we use these same primers with different adapter stubs to construct amplicon and restriction-site associated DNA libraries, but their use can be expanded to any type of library sequenced on Illumina platforms.

New insights on lake sediment DNA from the catchment: importance of taphonomic and analytical issues on the record quality

Over the last decade, an increasing number of studies have used lake sediment DNA to trace past landscape changes, agricultural activities or human presence. However, the processes responsible for lake sediment formation and sediment properties might affect DNA records via taphonomic and analytical processes. It is crucial to understand these processes to ensure reliable interpretations for “palaeo” studies. Here, we combined plant and mammal DNA metabarcoding analyses with sedimentological and geochemical analyses from three lake-catchment systems that are characterised by different erosion dynamics. The new insights derived from this approach elucidate and assess issues relating to DNA sources and transfer processes. The sources of eroded materials strongly affect the “catchment-DNA” concentration in the sediments. For instance, erosion of upper organic and organo-mineral soil horizons provides a higher amount of plant DNA in lake sediments than deep horizons, bare soils or glacial flours. Moreover, high erosion rates, along with a well-developed hydrographic network, are proposed as factors positively affecting the representation of the catchment flora. The development of open and agricultural landscapes, which favour the erosion, could thus bias the reconstructed landscape trajectory but help the record of these human activities. Regarding domestic animals, pastoral practices and animal behaviour might affect their DNA record because they control the type of source of DNA (“point” vs. “diffuse”).