Assessment of biodiversity using metabarcoding data, such as from bulk or environmental DNA sampling, is becoming increasingly relevant in ecology, biodiversity sciences and monitoring. Thereby, the taxonomic identification of species from their DNA sequences relies strongly on reference databases that link genetic sequences to taxonomic names. These databases vary in completeness and availability, depending on the taxonomic group studied and the genetic region targeted. The incompleteness of reference databases is an important argument to explain the nondetection by metabarcoding of species supposedly present. However, there exist further and generally overlooked problems with reference databases that can lead to false or inaccurate inferences of taxonomic assignment. Here, we synthesize all possible problems inherent to reference databases. In particular, we identify a complete, mutually nonexclusive list of seven classes of challenges when it comes to selecting, developing and using a reference database for taxonomic assignment. These are: (i) mislabelling, (ii) sequencing errors, (iii) sequence conflict, (iv) taxonomic conflict, (v) low taxonomic resolution, (vi) missing taxa and (vii) missing intraspecific variants. For each problem identified, we provide a description of possible consequences on the taxonomic assignment process. We illustrate the respective problem with examples taken from the literature or obtained by quantitative analyses of public databases, such as GenBank or BOLD. Finally, we discuss possible solutions to the identified problems and how to navigate them. Only by raising users' awareness of the limitations of metabarcoding data and DNA reference databases will adequate interpretations of these data be achieved.

Genetic markers have emerged as one of the most promising tools for species identification and geographic traceability in biodiversity conservation and international trade of biological products. However, traditional molecular markers rarely have sufficient resolution at lower taxonomic levels, especially for discriminating closely related forest tree species and their populations. In this study, we developed a panel of RNA-Seq based single nucleotide polymorphism (SNP) markers for tracing the geographic origin of an endangered conifer, Cathaya argyrophylla, which is a paleoendemic restricted to four mountain regions in subtropical China. A total of 69 individuals from five populations (DLS, SHS, HP, BMS, and DYS) covering the entire range were used for transcriptome sequencing. Based on these transcriptomic data, we evaluated genetic variation and population structure of C. argyrophylla, and found extremely low nucleotide diversity but strong population differentiation. We also screened 113 population-specific SNP loci, including 96 for BMS, eight for DYS, six for SHS, two for HP, and one for one of the three subpopulations from DLS. According to these geographically diagnostic SNPs, we designed four population-specific molecular barcodes for PCR amplification. To test the utility and efficiency of the four markers in geographic discrimination, double-blind experiment was performed using 157 individuals labelled without any locality information. We found that almost all tested individuals could be successfully assigned to their geographic localities. Our study not only sheds some new light on the genetic profile of C. argyrophylla, but also provides a practical and cost-efficient solution for geographic traceability using transcriptome-derived SNPs.

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.

Polyploids are cells or organisms with a genome consisting of more than two sets of homologous chromosomes. Polyploid plants have important traits that facilitate speciation and are thus often model systems for evolutionary, molecular ecology and agricultural studies. However, due to their unusual mode of inheritance and double-reduction, diploid models of population genetic analysis cannot properly be applied to autopolyploids. To overcome this problem, we developed a software package entitled vcfpop to perform a variety of population genetic analyses for autopolyploids, such as parentage analysis, analysis of molecular variance, principal coordinates analysis, hierarchical clustering analysis and Bayesian clustering. We used three data sets to evaluate the capability of vcfpop to analyse large data sets on a desktop computer. This software is freely available at http://github.com/huangkang1987/vcfpop.

Biomineralization-controlled exo-/endoskeleton growth contributes to body growth and body size diversity. Molluscan shells undergo ectopic biomineralization to form the exoskeleton and biocalcified "pearl" involved in invading defence. Notably, exo-/endoskeletons have a common ancestral origin, but their regulation and body growth are largely unknown. This study employed the pearl oyster, Pinctada fucata marntensii, a widely used experimental model for biomineralization in invertebrates, to perform whole-genome resequencing of 878 individuals from wild and breeding populations. This study characterized the genetic architecture of biomineralization-controlled growth and ectopic biomineralization. The insulin-like growth factor (IGF) endocrine signal interacted with ancient single-copy transcription factors to form the regulatory network. Moreover, the "cross-phylum" regulation of key long noncoding RNA (lncRNA) in bivalves and mammals indicated the conserved genetic and epigenetic regulation in exo-/endoskeleton growth. Thyroid hormone signal and apoptosis regulation in pearl oysters affected ectopic biomineralization in pearl oyster. These findings provide insights into the mechanism underlying the evolution and regulation of biomineralization in exo-/endoskeleton animals and ectopic biomineralization.

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.

Acipenseriformes (sturgeons and paddlefishes) are of substantial conservation concern, and development of genomic resources for these species is difficult due to past whole genome duplication. Development of disomic markers for polyploid organisms can be challenging due to difficulty in resolving alleles at a single locus from those among duplicated loci. In this study, we detail the development of disomic markers for the endangered pallid sturgeon (Scaphirhynchus albus) found in North America. One of the strategies for pallid sturgeon conservation is to stock U.S. rivers with offspring of pure pallid sturgeon, but introgression with the sympatric shovelnose sturgeon (S. platorynchus) threatens pallid sturgeon genetic integrity. Currently, 19 microsatellite loci are used to differentiate between both species and their hybrids, but the markers are insufficient to robustly identify backcrosses. We performed double digest restriction site-associated DNA sequencing (ddRADseq) on shovelnose sturgeon haploid gynogens to produce a reduced-representation genomic reference. Contiguous sequences that were heterozygous within a haploid individual were flagged as potentially encompassing multiple loci. Approximately 60 individuals of each species from two management units were sequenced, and reads were mapped to the haploid reference to identify single nucleotide polymorphisms (SNPs) at individual loci. The final data set contained 11,082 microhaplotyped loci which offer at least an order of magnitude greater resolution for species discrimination than the current panel of 19 microsatellites. These markers will be used to examine a larger sample of Scaphirhynchus individuals throughout their ranges to determine the extent and trajectory of hybridization.