Skin and gut microbiomes of tadpoles vary differently with host and water environment: a short-term experiment using 16S metabarcoding.

The host-microbiome community is influenced by several host and environmental factors. In order to disentangle the individual effects of host and environment, we performed a laboratory experiment to assess the effects of the exposure to different water sources on the skin and gut microbiome of two amphibian species (Pelophylax perezi and Bufo spinosus). We observed that the bacterial communities greatly varied with water environment and host identity. Tadpoles of B. spinosus collected from a waterbody with poorer bacterial diversity exhibited a more diverse skin and gut microbiome after exposed to a richer water source. Tadpoles of P. perezi, originally collected from a richer water environment, exhibited less marked alterations in diversity patterns independently of the water source but showed alterations in gut composition. These results highlight that environment alterations, such as the water source, combined with the host effect, impact the microbiome of amphibian species in different ways; the population history (e.g., previous water environment and habitat) of the host species may also influence future alterations on tadpole microbiome.

Automating microsatellite screening and primer design from multi-individual libraries using Micro-Primers.

Analysis of intra- and inter-population diversity has become important for defining the genetic status and distribution patterns of a species and a powerful tool for conservation programs, as high levels of inbreeding could lead into whole population extinction in few generations. Microsatellites (SSR) are commonly used in population studies but discovering highly variable regions across species' genomes requires demanding computation and laboratorial optimization. In this work, we combine next generation sequencing (NGS) with automatic computing to develop a genomic-oriented tool for characterizing SSRs at the population level. Herein, we describe a new Python pipeline, named Micro-Primers, designed to identify, and design PCR primers for amplification of SSR loci from a multi-individual microsatellite library. By combining commonly used programs for data cleaning and microsatellite mining, this pipeline easily generates, from a fastq file produced by high-throughput sequencing, standard information about the selected microsatellite loci, including the number of alleles in the population subset, and the melting temperature and respective PCR product of each primer set. Additionally, potential polymorphic loci can be identified based on the allele ranges observed in the population, to easily guide the selection of optimal markers for the species. Experimental results show that Micro-Primers significantly reduces processing time in comparison to manual analysis while keeping the same quality of the results. The elapsed times at each step can be longer depending on the number of sequences to analyze and, if not assisted, the selection of polymorphic loci from multiple individuals can represent a major bottleneck in population studies.

Assessing changes in stream macroinvertebrate communities across ecological gradients using morphological versus DNA metabarcoding approaches.

Freshwater macroinvertebrates provide valuable indicators for biomonitoring ecosystem change in relation to natural and anthropogenic drivers. DNA metabarcoding is an efficient approach for estimating such indicators, but its results may differ from morphotaxonomic approaches traditionally used in biomonitoring. Here we test the hypothesis that despite differences in the number and identity of taxa recorded, both approaches may retrieve comparable patterns of community change, and detect similar ecological gradients influencing such changes. We compared results obtained with morphological identification at family level of macroinvertebrates collected at 80 streams under a Water Framework Directive biomonitoring program in Portugal, with results obtained with metabarcoding from the ethanol preserving the bulk samples, using either single (COI-M19BR2, 16S-Inse01, 18S-Euka02) or multiple markers. Metabarcoding recorded less families and different communities compared to morphotaxonomy, but community sensitivities to disturbance estimated with the IASPT index were more similar across approaches. Spatial variation in local community metrics and the factors influencing such variation were significantly correlated between morphotaxonomy and metabarcoding. After reducing random noise in the dissimilarity matrices, the spatial variation in community composition was also significantly correlated across methods. A dominant gradient of community change was consistently retrieved, and all methods identified a largely similar set of anthropogenic stressors strongly influencing such gradient. Overall, results confirm our initial hypothesis, suggesting that morphotaxonomy and metabarcoding can estimate consistent spatial patterns of community variation and their main drivers. These results are encouraging for macroinvertebrate biomonitoring using metabarcoding approaches, suggesting that they can be intercalibrated with morphotaxonomic approaches to recover equivalent spatial and temporal gradients of ecological change.

Modelling technical and biological biases in macroinvertebrate community assessment from bulk preservative using multiple metabarcoding markers.

DNA metabarcoding from the ethanol used to store macroinvertebrate bulk samples is a convenient methodological option in molecular biodiversity assessment and biomonitoring of aquatic ecosystems, as it preserves specimens and reduces problems associated with sample sorting. However, this method may be affected by errors and biases, which need to be thoroughly quantified before it can be mainstreamed into biomonitoring programmes. Here, we used 80 unsorted macroinvertebrate samples collected in Portugal under a Water Framework Directive monitoring programme, to compare community diversity and taxonomic composition metrics estimated through morphotaxonomy versus metabarcoding from storage ethanol using three markers (COI-M19BR2, 16S-Inse01 and 18S-Euka02) and a multimarker approach. A preliminary in silico analysis showed that the three markers were adequate for the target taxa, with detection failures related primarily to the lack of adequate barcodes in public databases. Metabarcoding of ethanol samples retrieved far less taxa per site (alpha diversity) than morphotaxonomy, albeit with smaller differences for COI-M19BR2 and the multimarker approach, while estimates of taxa turnover (beta diversity) among sites were similar across methods. Using generalized linear mixed models, we found that after controlling for differences in read coverage across samples, the probability of detection of a taxon was positively related to its proportional abundance, and negatively so to the presence of heavily sclerotized exoskeleton (e.g., Coleoptera). Overall, using our experimental protocol with different template dilutions, the COI marker showed the best performance, but we recommend the use of a multimarker approach to detect a wider range of taxa in freshwater macroinvertebrate samples. Further methodological development and optimization efforts are needed to reduce biases associated with body armouring and rarity in some macroinvertebrate taxa.

Have the cake and eat it: Optimizing nondestructive DNA metabarcoding of macroinvertebrate samples for freshwater biomonitoring.

DNA metabarcoding can contribute to improving cost-effectiveness and accuracy of biological assessments of aquatic ecosystems, but significant optimization and standardization efforts are still required to mainstream its application into biomonitoring programmes. In assessments based on freshwater macroinvertebrates, a key challenge is that DNA is often extracted from cleaned, sorted and homogenized bulk samples, which is time-consuming and may be incompatible with sample preservation requirements of regulatory agencies. Here, we optimize and evaluate metabarcoding procedures based on DNA recovered from 96% ethanol used to preserve field samples and thus including potential PCR inhibitors and nontarget organisms. We sampled macroinvertebrates at five sites and subsampled the preservative ethanol at 1 to 14 days thereafter. DNA was extracted using column-based enzymatic (TISSUE) or mechanic (SOIL) protocols, or with a new magnetic-based enzymatic protocol (BEAD), and a 313-bp COI fragment was amplified. Metabarcoding detected at least 200 macroinvertebrate taxa, including most taxa detected through morphology and for which there was a reference barcode. Better results were obtained with BEAD than SOIL or TISSUE, and with subsamples taken 7-14 than 1-7 days after sampling, in terms of DNA concentration and integrity, taxa diversity and matching between metabarcoding and morphology. Most variation in community composition was explained by differences among sites, with small but significant contributions of subsampling day and extraction method, and negligible contributions of extraction and PCR replication. Our methods enhance reliability of preservative ethanol as a potential source of DNA for macroinvertebrate metabarcoding, with a strong potential application in freshwater biomonitoring.

Differential effects of dietary protein on early life-history and morphological traits in natterjack toad (Epidalea calamita) tadpoles reared in captivity.

The production of high quality amphibian larvae through optimal diets is a critical component of amphibian conservation breeding programs. Larval period, survival, body weight and total length are frequently used as metrics of adequate nutrition. However, the effects of nutrition on tadpole and metamorph morphology are rarely tested in detail. In the present study, we analyzed the most common metrics and six other larval and post-metamorphic morphological traits in natterjack toads (Epidalea calamita) fed with three different commercial fish diets, varying in protein content (32.0%, 38.3%, and 46.2%). Our results suggest that early life-history (tadpole growth, development, and survival) and morphological traits of E. calamita tadpoles are differentially affected by the percentage of dietary protein. As protein content increased, tadpoles exhibited larger bodies along with shorter tail fins; however, with no significant differences in total length. Larval period was similar across treatments but mortality was lower in high-protein diet. At high-protein diets the metamorphs revealed significantly longer bodies, and wider heads and hind legs, but there was no significant difference in the average weight across all dietary treatments. Based on our results, feed containing 46.2% protein promotes growth, development and survival of E. calamita tadpoles better than either of the other two feeds tested. The use of other body measures beyond weight, tadpole total length, and snout-vent length in studies of amphibian nutrition in captivity may assist the selection of appropriate diets to optimize tadpole survival and metamorph fitness.