Efficient COI barcoding using high throughput single-end 400 bp sequencing

Development of mitochondrial DNA cytochrome c oxidase subunit I primer sets to construct DNA barcoding library using next-generation sequencing

Insects are one of the most diverse eukaryotic groups on the planet, with one million or more species present, including those yet undescribed. The DNA barcoding system has been developed, which has aided in the identification of cryptic species and undescribed species. The mitochondrial cytochrome c oxidase I region (mtDNA COI) has been utilised for the barcoding analysis of insect taxa. Thereafter, next-generation sequencing (NGS) technology has been developed, allowing for rapid acquisition of massive amounts of sequence data for genetic analyses. Although NGS-based PCR primers designed to amplify the mtDNA COI region have been developed, their target regions were only a part of COI region and/or there were taxonomic bias for PCR amplification. As the mtDNA COI region is a traditional DNA marker for the DNA barcoding system, modified primers for this region would greatly contribute to taxonomic studies. In this study, we redesigned previously developed PCR primer sets that targetted the mtDNA COI barcoding region to improve amplification efficiency and to enable us to conduct sequencing analysis on NGS. As a result, the redesigned primer sets achieved a high success rate (> 85%) for species examined in this study, covering four insect orders (Coleoptera, Lepidoptera, Orthoptera and Odonata). Thus, by combining the primers with developed primer sets for 12S or 16S rRNA regions, we can conduct more detailed taxonomic, phylogeographic and conservation genetic studies using NGS.

Application and Comparison of Machine Learning and Database-Based Methods in Taxonomic Classification of High-Throughput Sequencing Data

The advent of high-throughput sequencing technologies has not only revolutionized the field of bioinformatics but has also heightened the demand for efficient taxonomic classification. Despite technological advancements, efficiently processing and analyzing the deluge of sequencing data for precise taxonomic classification remains a formidable challenge. Existing classification approaches primarily fall into two categories, database-based methods and machine learning methods, each presenting its own set of challenges and advantages. On this basis, the aim of our study was to conduct a comparative analysis between these two methods while also investigating the merits of integrating multiple database-based methods. Through an in-depth comparative study, we evaluated the performance of both methodological categories in taxonomic classification by utilizing simulated data sets. Our analysis revealed that database-based methods excel in classification accuracy when backed by a rich and comprehensive reference database. Conversely, while machine learning methods show superior performance in scenarios where reference sequences are sparse or lacking, they generally show inferior performance compared with database methods under most conditions. Moreover, our study confirms that integrating multiple database-based methods does, in fact, enhance classification accuracy. These findings shed new light on the taxonomic classification of high-throughput sequencing data and bear substantial implications for the future development of computational biology. For those interested in further exploring our methods, the source code of this study is publicly available on https://github.com/LoadStar822/Genome-Classifier-Performance-Evaluator. Additionally, a dedicated webpage showcasing our collected database, data sets, and various classification software can be found at http://lab.malab.cn/~tqz/project/taxonomic/.

Castaways: the Leeward Antilles endemic spider genus Papiamenta (Araneae: Pholcidae)

Ninetinae is a group of small to tiny short-legged spiders largely restricted to arid habitats. Among daddy-long-legs spiders (Pholcidae) this is by far the least diverse subfamily but this may partly be a result of inadequate collecting, poor representation in collections or scientific neglect. We build on a large recent collection of the ninetine genus Papiamenta Huber, 2000 from the Leeward Antilles and use cytochrome oxidase 1 (COI ) sequences, extensive scanning electron microscopy data, transmission electron microscopy data and karyotyping to analyse this geographically isolated and poorly known island genus. COI sequences support the split between the two morphologically distinct species on Curaçao but genetic distances between these are surprisingly low (7.4-9.8%; mean 8.6%). The type species P. levii (Gertsch, 1982) may include more than one species but COI and morphology suggest conflicting clade limits. A third species, P. bonay Huber sp. nov. is newly described from Bonaire. Our data on sperm ultrastructure and karyology are puzzling as these suggest different phylogenetic affinities of Papiamenta to other genera. Males transfer sperm as individual sperm (cleistosperm), agreeing with the putative closest relatives as suggested by molecular data, the North American genera Pholcophora and Tolteca . The sex chromosome system (X 1 X 2 X 3 Y ) of P. levii , however, is as in the South American Ninetinae genera Gertschiola and Nerudia but different from the putative closest relatives. ZooBank: urn:lsid:zoobank.org:pub:7A6A2E84-3A61-4637-AF6F-0E31A9FA79A8.

taxalogue: a toolkit to create comprehensive CO1 reference databases

Background

Taxonomic identification through DNA barcodes gained considerable traction through the invention of next-generation sequencing and DNA metabarcoding. Metabarcoding allows for the simultaneous identification of thousands of organisms from bulk samples with high taxonomic resolution. However, reliable identifications can only be achieved with comprehensive and curated reference databases. Therefore, custom reference databases are often created to meet the needs of specific research questions. Due to taxonomic inconsistencies, formatting issues, and technical difficulties, building a custom reference database requires tremendous effort. Here, we present taxalogue, an easy-to-use software for creating comprehensive and customized reference databases that provide clean and taxonomically harmonized records. In combination with extensive geographical filtering options, taxalogue opens up new possibilities for generating and testing evolutionary hypotheses.

Methods

taxalogue collects DNA sequences from several online sources and combines them into a reference database. Taxonomic incongruencies between the different data sources can be harmonized according to available taxonomies. Dereplication and various filtering options are available regarding sequence quality or metadata information. taxalogue is implemented in the open-source Ruby programming language, and the source code is available at https://github.com/nwnoll/taxalogue. We benchmark four reference databases by sequence identity against eight queries from different localities and trapping devices. Subsamples from each reference database were used to compare how well another one is covered.

Results

taxalogue produces reference databases with the best coverage at high identities for most tested queries, enabling more accurate, reliable predictions with higher certainty than the other benchmarked reference databases. Additionally, the performance of taxalogue is more consistent while providing good coverage for a variety of habitats, regions, and sampling methods. taxalogue simplifies the creation of reference databases and makes the process reproducible and transparent. Multiple available output formats for commonly used downstream applications facilitate the easy adoption of taxalogue in many different software pipelines. The resulting reference databases improve the taxonomic classification accuracy through high coverage of the query sequences at high identities.

Development of nucleotide signatures for common poisonous organisms provides a new strategy for food poisoning diagnosis

DNA barcoding is widely used in toxic species authentication, but due to serious DNA degradation of forensic materials, the application of full-length barcode sequences in food poisoning diagnosis is greatly limited. Nucleotide signature, a shorter specific molecular marker, derived from traditional DNA barcoding has been proposed as an emerging tool of toxic species detection in deeply processed materials. In this study, to resolve the frequent food poisoning accidents with unknown origin, we envisioned developing a nucleotide signature data set of common poisonous organisms and combining high-throughput sequencing (HTS) to reveal the poisoning cause. Ninety-three individuals and 1093 DNA barcode sequences of twelve common poisonous plants, fish, mushrooms and their related species were collected. Through sequence alignment and screening, the nucleotide signatures were respectively developed and validated as their specific molecular markers. The sequence length varied from 19 bp to 38 bp. These fragments were conserved within the same species or genera, and the specificity between related species has been also demonstrated. To further evaluate the application potential of nucleotide signature in forensic diagnosis, simulated forensic specimens (SFS) containing different poisonous ingredients were sequenced by HTS with PCR-free libraries. As a result, the nucleotide signature was successfully captured from original HTS data without assembly and annotation, accompanied by a high detection sensitivity of 0.1 ng/µl in mixture system. Therefore, this method was suitable for the assay of forensic materials with serious DNA degradation. The present study undoubtedly provides a new perspective and strong support for the detection of toxic ingredients and the diagnosis of food poisoning.

Molecular parallelism in the evolution of a master sex-determining role for the anti-Mullerian hormone receptor 2 gene (amhr2) in Midas cichlids

The evolution of sex chromosomes and their differentiation from autosomes is a major event during genome evolution that happened many times in several lineages. The repeated evolution and lability of sex-determination mechanisms in fishes makes this a well-suited system to test for general patterns in evolution. According to current theory, differentiation is triggered by the suppression of recombination following the evolution of a new master sex-determining gene. However, the molecular mechanisms that establish recombination suppression are known from few examples, owing to the intrinsic difficulties of assembling sex-determining regions (SDRs). The development of forward-genetics and long-read sequencing have generated a wealth of data questioning central aspects of the current theory. Here, we demonstrate that sex in Midas cichlids is determined by an XY system, and identify and assemble the SDR by combining forward-genetics, long-read sequencing and optical mapping. We show how long-reads aid in the detection of artefacts in genotype-phenotype mapping that arise from incomplete genome assemblies. The male-specific region is restricted to a 100-kb segment on chromosome 4 that harbours transposable elements and a Y-specific duplicate of the anti-Mullerian receptor 2 gene, which has evolved master sex-determining functions repeatedly. Our data suggest that amhr2Y originated by an interchromosomal translocation from chromosome 20 to 4 pre-dating the split of Midas and Flier cichlids. In the latter, it is pseudogenized and translocated to another chromosome. Duplication of anti-Mullerian genes is a common route to establishing new sex determiners, highlighting the role of molecular parallelism in the evolution of sex determination.

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.

Testing cophylogeny between coral reef invertebrates and their bacterial and archaeal symbionts

Marine invertebrates harbour a complex suite of bacterial and archaeal symbionts, a subset of which are probably linked to host health and homeostasis. Within a complex microbiome it can be difficult to tease apart beneficial or parasitic symbionts from nonessential commensal or transient microorganisms; however, one approach is to detect strong cophylogenetic patterns between microbial lineages and their respective hosts. We employed the Procrustean approach to cophylogeny (PACo) on 16S rRNA gene derived microbial community profiles paired with COI, 18S rRNA and ITS1 host phylogenies. Second, we undertook a network analysis to identify groups of microbes that were co-occurring within our host species. Across 12 coral, 10 octocoral and five sponge species, each host group and their core microbiota (50% prevalence within host species replicates) had a significant fit to the cophylogenetic model. Independent assessment of each microbial genus and family found that bacteria and archaea affiliated to Endozoicomonadaceae, Spirochaetaceae and Nitrosopumilaceae have the strongest cophylogenetic signals. Further, local Moran's I measure of spatial autocorrelation identified 14 ASVs, including Endozoicomonadaceae and Spirochaetaceae, whose distributions were significantly clustered by host phylogeny. Four co-occurring subnetworks were identified, each of which was dominant in a different host group. Endozoicomonadaceae and Spirochaetaceae ASVs were abundant among the subnetworks, particularly one subnetwork that was exclusively comprised of these two bacterial families and dominated the octocoral microbiota. Our results disentangle key microbial interactions that occur within complex microbiomes and reveal long-standing, essential microbial symbioses in coral reef invertebrates.

Testing the utility of alternative metrics of branch support to address the ancient evolutionary radiation of tunas, stromateoids, and allies (Teleostei: Pelagiaria)

The use of high-throughput sequencing technologies to produce genome-scale datasets was expected to settle some long-standing controversies across the Tree of Life, particularly in areas where short branches occur at deep timescales. Instead, these datasets have often yielded many well-supported but conflicting topologies, and highly variable gene-tree distributions. A variety of branch-support metrics beyond the nonparametric bootstrap are now available to assess how robust a phylogenetic hypothesis may be, as well as new methods to quantify gene-tree discordance. We applied multiple branch support metrics to an ancient group of marine fishes (Teleostei: Pelagiaria) whose interfamilial relationships have proven difficult to resolve due to a rapid accumulation of lineages very early in its history. We analyzed hundreds of loci including published UCE data and newly generated exonic data along with their flanking regions to represent all 16 extant families for more than 150 out of 284 valid species in the group. Branch support was lower for interfamilial relationships (except the SH-like aLRT and aBayes methods) regardless of the type of marker used. Several nodes that were highly supported with bootstrap had very low site and gene-tree concordance, revealing underlying conflict. Despite this conflict, we were able to identify four consistent interfamilial clades, each comprised of two or three families. Combining exons with their flanking regions also produced increased branch lengths in the deep branches of the pelagiarian tree. Our results demonstrate the limitations of employing current metrics of branch support and species-tree estimation when assessing the confidence of ancient evolutionary radiations and emphasize the necessity to embrace alternative measurements to explore phylogenetic uncertainty and discordance in phylogenomic datasets.