Next Article in Journal
Microarray-Based Prediction of Polycythemia after Exposure to High Altitudes
Next Article in Special Issue
Bioinformatic Characterization and Molecular Evolution of the Lucina pectinata Hemoglobins
Previous Article in Journal
TTC30A and TTC30B Redundancy Protects IFT Complex B Integrity and Its Pivotal Role in Ciliogenesis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 13
Issue 7
10.3390/genes13071192
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Testing Efficacy of Assembly-Free and Alignment-Free Methods for Species Identification Using Genome Skims, with Patellogastropoda as a Test Case

by
Tao Xu
1,
Lingfeng Kong
1,2,* and
Qi Li
1,2
1
Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, 5 Yushan Road, Qingdao 266003, China
2
Laboratory for Marine Fisheries Science and Food Production Processes, Qingdao National Laboratory for Marine Science and Technology, 5 Yushan Road, Qingdao 266003, China
*
Author to whom correspondence should be addressed.
Genes 2022, 13(7), 1192; https://doi.org/10.3390/genes13071192
Submission received: 17 April 2022 / Revised: 26 June 2022 / Accepted: 28 June 2022 / Published: 2 July 2022
(This article belongs to the Special Issue Genetic Evolution of Marine Shellfish)
Browse Figures
Versions Notes

Abstract

:
Most recently, species identification has leaped from DNA barcoding into shotgun sequencing-based “genome skimming” alternatives. Genome skims have mainly been used to assemble organelle genomes, which discards much of the nuclear genome. Recently, an alternative approach was proposed for sample identification, using unassembled genome skims, which can effectively improve phylogenetic signal and identification resolution. Studies have shown that the software Skmer and APPLES work well at estimating genomic distance and performing phylogenetic placement in birds and insects using low-coverage genome skims. In this study, we use Skmer and APPLES based on genome skims of 11 patellogastropods to perform assembly-free and alignment-free species identification and phylogenetic placement. Whether or not data corresponding to query species are present in the reference database, Skmer selects the best matching or closest species with COI barcodes under different sizes of genome skims except lacking species belonging to the same family as a query. APPLES cannot place patellogastropods in the correct phylogenetic position when the reference database is sparse. Our study represents the first attempt at assembly-free and alignment-free species identification of marine mollusks using genome skims, demonstrating its feasibility for patellogastropod species identification and flanking the necessity of establishing a database to share genome skims.

1. Introduction

The deterioration of the global ecosystems has accelerated biodiversity loss in recent years, including undiscovered germplasm resources. Rapid and inexpensive taxonomic identification to discover and protect biodiversity has become a hot topic among taxonomists. Currently, the standardized and broadly used method of species identification is DNA barcoding (e.g., COI, 12/16S, matK, ITS) [1,2,3,4], which is more efficient and accurate than previously widely used morphological identification methods. There is no doubt that DNA barcoding is the fundamental pillar of many current and future studies [5]. However, the high-quality DNA required for PCR amplification limits the available specimens, and the limited phylogenetic signals prevent barcoding from distinguishing all species (e.g., wasp) [6]. Moreover, although the Barcode of Life Data System (BOLD) (https://www.boldsystems.org/, accessed on 12 January 2022) provides a reliable database to assign an identity to query samples using a reference database of taxonomically preidentified vouchers, the barcode sequences contained in BOLD are only part of the biodiversity [7].
Recently, the decreasing costs of shotgun sequencing have led researchers to propose an alternative method, that is, using low-pass sequencing to generate genome skims [8,9], typically around 1–5 Gb [8,10], providing 0.1–10× coverage, and usually not enough to recover the nuclear genome [11]. Reconstructing organelle genomes from low-pass sequencing data to perform de novo reconstruction (e.g., maximum likelihood (ML) tree) is the most common application of genome skimming. There is no doubt that organelle genomes will provide a greater phylogenetic signal and identification resolution than barcoding in species identification [12,13,14,15]. However, the relatively time-consuming manual curation steps (i.e., assembly and alignment) and limited scalability prevent this method from being applied to ultra-large trees with higher accuracy and resolution. Furthermore, reconstructed phylogenetic relationships based on the mitochondrial genome are easily affected by long-branch attraction (LBA), which leads to the clustering of rapidly evolving taxa and increases the difficulty of reconstructing deep (Cambrian) molluscan evolution, especially with limited taxon sampling [16]. Most importantly, the above approach discards a vast proportion of the nuclear genome (as much as 99% of the sequence data) [5,17], making low-pass shotgun data underutilized.
The nuclear genome represents the evolutionary history of any nonclonal organism [5]. As the desirable approach to identifying samples, acquiring fully assembled nuclear genomes requires higher sequencing depth (at least 50×) and computational power [5,18]. In addition to the high sequencing cost and computing requirements, repeat elements of nuclear genomes will prevent unambiguous assembly when longer reads are sequenced. Although there are existing solutions to this problem (e.g., construction of mate-pair/large insert libraries for short-read technologies, high molecular weight DNA extraction, and long-read sequencing using single-molecule sequencing), the methods limit the available samples and require complex equipment and skills [5]. Therefore, a fully assembled nuclear genome is not presently suitable for extensive use in species identification.
Nevertheless, given that nuclear genome sequences represent the ultimate source of information for taxonomic assignment, some recent studies have proposed using all unassembled reads from low-pass sequencing data to perform taxonomic assignments [5,8,17]. It is argued that for a genome of size n and ignoring repeats, the probability of finding a k-mer with sufficient size (log4 n) in another genome relates directly to the evolutionary distance to the other genome [17]. Therefore, assembly-free and alignment-free species identification using genome skims might be a viable alternative to DNA barcoding. There is already some assembly-free or alignment-free software available [19,20,21], but this either requires high enough coverage (e.g., Co-phylog, Mash [21] and Simka) or the accuracy of the results is not high enough (e.g., AAF [20]). Recently, Sarmashghi et al. proposed a new assembly-free and alignment-free method for species identification and developed Skmer. This software can accurately compute the genomic distance from low-coverage genome skims [17]. Subsequently, Balaban et al. developed the software APPLES (Accurate Phylogenetic Placement using Least Squares) for distance-based phylogenetic placement, which can improve identification accuracy [22]. APPLES can find the optimal position for a new query species on an existing backbone (or reference) tree, which relies on the assembly-free estimates of genomic distance estimated from low coverage genome skims by Skmer or other alternatives. Compared to Maximum likelihood (ML) methods, APPLES is more memory efficient, more scalable, and runs faster.
Although both software performed well in the tests of birds and insects [17,22], broader testing is currently lacking. Patellogastropoda has been recognized as the most ‘primitive’ group of living gastropods. Due to historical vicariance and dispersal of the Pangean supercontinent, patellogastropod species with low dispersal ability have undergone geographic isolation and diverged on isolated islands, resulting in abundant cryptic species [23,24].
In this study, we sequenced the low-pass whole genome data of 11 patellogastropod species using genome skimming. Our aims were (1) to test the assembly-free and alignment-free species identification effectiveness of Skmer whether data correspond to query species in the reference database or not; (2) to test the phylogenetic placement effectiveness of APPLES in the above two cases; (3) to determine the minimum size of genome skims needed to obtain reliable species identification and phylogenetic placement results.

2. Material and Methods

2.1. DNA Extraction and Sequencing

The collecting site of each specimen is shown in Table 1. All specimens were immediately preserved in 95% ethanol following collection. The total genomic DNA was extracted using the TIANamp Marine Animals DNA Kit (TIANGEN Biotech Beijing Co., Ltd., Beijing, China), following the manufacturer’s protocols. Genomic DNA was sequenced by Novogene Technology Co., Ltd. (Beijing, China) on the Illumina NovaSeq 6000 platform using a PE150 protocol. To compare the species identification effectiveness of Skmer and APPLES between whether there is data corresponding to query species in the reference database or not, we selected Scutellastra flexuosa, Cellana toreuma (HN), and Patelloida conulus as the query species and other patellogastropod species constituted the reference library. S. flexuosa is the only species that belongs to Patellidae, while other families have at least three species in our study. C. toreuma (HN) and C. toreuma (GD) are specimens collected from different sampling sites. P. conulus, Patelloida ryukyuensis, and Patelloida saccharina lanx belong to the genus Patelloida. Query species from different scenarios allow us to better evaluate the application of assembly-free and alignment-free species identification using Skmer and APPLES in marine shellfish.

2.2. Genomic Reads Subsampling and Preprocessing

First, genome skims with 0.1 Gb, 0.5 Gb, 1 Gb, 2 Gb, and 4 Gb of data were generated from the paired-end reads using BBTools [25] by randomly subsampling the reads. As L. goshimai has the smallest size (7.3 Gb) of clean data in the study (Table S1), we used BBTools to subsample the paired-end reads of the remaining species to 7.3 Gb and called this the largest data in this article. Then we used fastp [26] to filter low-quality reads and remove adapters. Kraken II; [27] was used to filter possible extraneous reads such as bacteria, archaea, viral and human contamination. After filtering, BBtools was used to clean up deduplicate reads and merge paired-end reads. The statistics of the processing results for each step are shown in Table S3. Through this pipeline, we obtained the six datasets for the 11 patellogastropod species, and Trochus maculatus was selected as the outgroup.

2.3. Distance Calculation and Phylogenetic Placement

The process of Skmer calculation distance is as follows: First, Jelly Fish [28] was used to compute the frequency profiles of the subset of genome skims (i.e., 0.1 Gb, 0.5 Gb, 1 Gb, 2 Gb, 4 Gb, or largest data), and then to estimate the coverage, genome length, error rate, and read length, which can help analyze the test results further. Second, we used the hashing technique of Mash to retain a subset of the above subset of genome skims, and then it was used to compute the Jaccard index. Finally, we used these estimates to compute the genomic distance between the query and reference. To compare with the DNA barcoding method, COI barcodes for each species were downloaded from the BOLD database. MEGA v. 5.1 [29] was used to align all barcodes and calculate the pairwise p-distance.
After obtaining the distance matrix of reference samples, we transformed the genomic distances to Jukes–Cantor (JC) distances. Then, we used FastME [30] to infer the backbone tree, which is the necessary input for APPLES. Finally, APPLES used Treeswift [31] to place the query on the optimal position of the backbone tree based on the distance matrices mentioned above and used the JC69 model to independently compute phylogenetic distances [32] without the Γ model of rate variation for all pairs. Phylogenetic trees were visualized in FigTree v. 1.4.4 [33]. To further evaluate the effectiveness of the method, the reference mitogenome trees using exactly the same species and the outgroup were reconstructed based on 13 mitochondrial protein-coding genes (PCGs, refer to Xu et al. [34] for the specific method).

3. Results

Our subsampled genome skims ranged in coverage from 0.50× to 1.68×, 0.69× to 1.13×, 0.63× to 1.15×, 0.69× to 1.92×, and 0.93× to 2.78× for 0.5 Gb, 1 Gb, 2 Gb, 4 Gb, and largest data, respectively (Figure 1; Table S2). However, Skmer was unable to estimate the coverage of 0.1 Gb data. We then compared estimated distances computed from COI and unassembled shotgun sequence data.
First, we select C. toreuma (HN) as the query species. Skmer correctly identified the best match under different sizes of genome skims with very small differences in the genomic distance except for 0.1 Gb (Table 2). Then, when S. flexuosa was used as the query species, the calculation results from different sizes of genome skims agreed that P. saccharina lanx was the closest species to the query skims in our reference samples. However, the pairwise distance of all COI barcodes calculated by MEGA v. 5.1 showed that C. grata and C. toreuma have a smaller genetic distance with S. flexuosa than P. saccharina lanx (Table 3). For P. conulus, both COI and genome skims agreed that the closest species was P. ryukyuensis (Table 4).
Although Nacellidae and Lottiidae were a monophyletic group and clustered into a clade in most backbone trees (Figures S1 and S2), subtle differences existed in the phylogenetic relationships of backbone trees inferred by FastME based on specific sizes of genome skims. Specifically, the backbone trees of 0.5 Gb and 1 Gb data placed Nacellidae as a paraphyletic grade at the base of patellogastropods in Figure S1. Nacellidae nested in the family Lottiidae in the backbone tree inferred based on 2 Gb genome skims. In Figure S2, Lottiidae nested within Nacellidae in the reference tree of 0.5 Gb. However, the backbone trees with the same phylogenetic relationship had different branch lengths due to different pairwise distances to query samples.
APPLES placed C. toreuma (HN) into a clade with C. toreuma (GD) except for the 0.5 Gb data (Figure 2). The phylogenetic position of S. flexuosa was variable in the results obtained based on different sizes of genome skims (Figure 3). APPLES placed S. flexuosa at or near the base of Patellogastropoda except for the phylogenetic results of 4 Gb and the largest data, which placed the query into a sister clade with L. cassis and N. radula. When P. conulus was selected as the query species, it formed a clade with P. ryukyuensis or P. ryukyuensis + L. cassis under 0.1 Gb, 0.5 Gb, and 1 Gb data, while it was sister to a clade formed by N. radula and L. cassis in the placement results based on data volumes of 2 Gb and above (Figure 4).

4. Discussion

We present here the first assembly-free and alignment-free species identification using genome skims on Patellogastropoda. Skmer consistently selected the best matching or closest species under different data amounts, whether there were data corresponding to the query in the reference database or not. However, when species belonging to the same family were lacking, the taxon selection of COI barcodes was different from that of genome skims. APPLES placed query species in the accurate phylogenetic position only when the corresponding species data were available in the reference database.

4.1. Coverage

In our study, the coverage of genome skims estimated by Skmer was less than 3×, even in the largest data. The coverage was not uniform; it was randomly distributed and did not increase exponentially as the genome skims increased, which may be due to the overrepresentation of mitochondrial sequences (Figure 1, Table S2) [17].

4.2. Selection of the Best Matched Species

In our study, Skmer selected the same species that best matched the query sequence in the reference datasets even based on different sizes of genome skims. For C. toreuma (HN), Skmer selected the correct best matching species. The closest species of P. conulus selected based on different sizes of genome skims and COI was consistent. Interestingly, P. saccharina lanx, which belongs to Lottiidae, was selected by Skmer as the best matching species for S. flexuosa, but the result of COI barcodes showed that C. grata and C. toreuma had a much closer genomics distance with S. flexuosa than P. saccharina lanx.
According to the estimated distances, we found no correlation between the results of COI and genome skims even though both of them selected the same best matching species (Table 2, Table 3 and Table 4). There are several possible reasons for this discrepancy between COI barcoding and genome skims. First, the length of COI barcoding is around 658 bp, which can only provide limited sequence information and low phylogenetic resolution [35]. In contrast, genome skims not only include all of the different ‘standard’ animal barcoding regions (e.g., COI, 12S, 16S), but also provide sequence data from other loci [8,36], which hold valuable information that can further achieve the goals of species identification. In genome skimming approaches, as much as 99% of the sequence data is from the nuclear genome [5,17,37]. Considering this fact, we speculate that the second reason could be the great difference in sequence information provided between organelle and nuclear genome, which may result in gene tree/species tree discordances. For example, Patellogastropoda and Vetigastropoda were recovered as the sister clade of the remaining gastropods based on transcriptomes [38], while Patellogastropoda was the sister lineage of the remaining gastropods when reconstructed based on the mitochondrial genome [16]. The third possible reason is that Patellogastropoda might have complex mutations, such as large-scale repetition, especially S. flexuosa. When calculating distance based on genome skims, Skmer will simplify the evolutionary process, such as ignoring repeats and assuming that mutations are uniformly distributed [5,17], which might reduce the accuracy of the results of the calculation about the patellogastropods. While in the calculation with COI barcodes, the p-distance or Kimura 2-parameter model is usually chosen according to the affinity between the query species and the reference species. The potential inherent biases in the methods may also explain the observed differences between COI and genome skims.
In addition, a crucial step before the application of assembly-free and alignment-free species identification is to use Kraken-II to remove possible contamination reads after subsampling. Several facts about Kraken-II might affect the accuracy of calculated distance. Firstly, studies have shown that Kraken-II can effectively reduce the adverse effects of contamination only when the contaminants and the contaminant reference library have a match, which is within a 5–10% genomic distance [5,11]. Therefore, an incomplete contaminant reference library and unconfident matching to the contaminant database may affect our estimation results. Secondly, impure query or reference skims can lead to underestimating the accurate distance by Skmer, especially when the impurity of the query skims is similar to that of the reference skims [5,11]. Thirdly, Kraken-II has been shown to occasionally over-correct errors, which might lead to an overestimation of the true distance [11].

4.3. Phylogenetic Placement

According to the phylogenetic placement result, APPLES can accurately determine the location of the query species above 1 Gb genome skims when the corresponding species data are available in the reference database. Although stable phylogenetic positions were obtained based on 4 Gb and above data, the results were irrational when the reference database did not have corresponding query data. The phylogenetic tree reconstructed by mitogenomes (Figure S3; [34]), mitochondrial and nuclear genes [23] showed that P. conulus has a closer relationship with P. ryukyuensis. However, APPLES placed it in the sister clade of N. radula and L. cassis. S. flexuosa, which lacks species belonging to the Patellidae family in the backbone tree, was nested in Lottiidae and then formed a clade with N. radula and L. cassis. However, in the phylogenetic reconstruction based on the mitochondrial genome, S. flexuosa was clustered into a clade with Nacellidae (Figure S4; [34]) and was placed in the most basal position in the molecular phylogeny, including the most extensive sampling of specimens [23]. The results indicated that APPLES could not place P. conulus and S. flexuosa in the exact phylogenetic position.
In addition to the possible influencing factors mentioned above in the Skmer section, we considered that the misplacement of P. conulus and S. flexuosa might also be influenced by sparse taxon sampling, which might decrease the accuracy of APPLES. Similarly, in the phylogenetic inference of mitochondrial genomes, the evolutionary relationship of L. digitalis suffers from long-branch attraction (LBA), which results in inconsistent relationships among different studies [16,39]. Nevertheless, it has recently been proven that improved taxon sampling can effectively alleviate the LBA [34]. Future studies using a denser taxon sampling, especially in Patelloida and Patellidae, can further confirm our hypothesis.
The above evidence shows that APPLES cannot place patellogastropods in a reliable phylogenetic position on a sparse backbone tree.

5. Conclusions

In our study, the assembly-free and alignment-free methods for species identification using genome skims perform well in Patellogastropoda, meaning that Skmer has good potential for application in more taxa besides insects and birds. Problems that existed in the phylogenetic placement of APPLES might be affected by limited taxon sampling and need to be further discussed by increasing the number of species in Patelloida and Patellidae. More comparative studies covering denser sampling and different groups of mollusks should be implemented in the future. Furthermore, to better apply this new species identification method, we propose establishing a large reference database to store the processed shotgun sequencing data.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes13071192/s1, Figure S1: The backbone tree of C. toreuma (HN) and S. flexuosa inferred by FastME under different size of genome skims; Figure S2: The backbone tree of P. conulus inferred by FastME under different size of genome skims; Figure S3: Phylogenetic relationships of P. conulus (query species) inferred based on concatenated amino acids of 13 mitochondrial protein-coding genes. Numbers at nodes are statistical support values for BI (posterior probabilities)/ML (bootstrap proportions in percentage); Figure S4: Phylogenetic relationships of S. flexuosa (query species) inferred based on concatenated amino acids of 13 mitochondrial protein-coding genes. Numbers at nodes are statistical support values for BI (posterior probabilities)/ML (bootstrap proportions in percentage); Table S1: The original paired-end reads; Table S2: The parameters of different size of genome skims calculated by Skmer; Table S3: The statistics of the processing results for BBTools, fastp and Kraken II.

Author Contributions

Writing—original draft preparation, visualization, writing—review and editing, T.X.; supervisor: L.K. and Q.L.; funding acquisition: L.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China (2019YFD0900403); National Natural Science Foundation of China: 31772414; Fundamental Research Funds for the Central Universities: 201964001.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original shotgun sequencing data used in this study were deposited in the NCBI Sequence Read Archive (BioProject PRJNA766309, accessed on 7 June 2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hebert, P.; Cywinska, A.; Ball, S.; Dewaard, J. Biological identifications through DNA barcodes. Proc. R Soc. B Biol. Sci. 2003, 270, 313–321. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Savolainen, V.; Cowan, R.S.; Vogler, A.P.; Roderick, G.K.; Lane, R. Towards writing the encyclopaedia of life: An introduction to DNA barcoding. Phil. Trans. R. Soc. B 2005, 360, 1805–1811. [Google Scholar] [CrossRef] [PubMed]
  3. Schindel, D.E.; Miller, S.E. DNA barcoding a useful tool for taxonomists. Nature 2005, 435, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Taberlet, P.; Coissac, E.; Pompanon, F.; Brochmann, C.; Willerslev, E. Towards next-generation biodiversity assessment using DNA metabarcoding. Mol. Ecol. 2012, 21, 2045–2050. [Google Scholar] [CrossRef]
  5. Bohmann, K.; Mirarab, S.; Bafna, V.; Gilbert, M.T.P. Beyond DNA barcoding: The unrealized potential of genome skim data in sample identification. Mol. Ecol. 2020, 29, 2521–2534. [Google Scholar] [CrossRef]
  6. Quicke, D.L.J.; Smith, M.A.; Janzen, D.H.; Hallwachs, W.; Fernandez-Triana, J.; Laurenne, N.M.; Zaldívar-Riverón, A.; Shaw, M.R.; Broad, G.R.; Klopfstein, S.; et al. Utility of the DNA barcoding gene fragment for parasitic wasp phylogeny (Hymenoptera: Ichneumonoidea): Data release and new measure of taxonomic congruence. Mol. Ecol. Resour. 2012, 12, 676–685. [Google Scholar] [CrossRef]
  7. Ratnasingham, S.; Hebert, P.D.N. BOLD: The barcode of life data system (www.Barcodinglife.Org). Mol. Ecol. Notes 2007, 7, 355–364. [Google Scholar] [CrossRef] [Green Version]
  8. Coissac, E.; Hollingsworth, P.M.; Lavergne, S.; Taberlet, P. From barcodes to genomes: Extending the concept of DNA barcoding. Mol. Ecol. 2016, 25, 1423–1428. [Google Scholar] [CrossRef] [Green Version]
  9. Straub, S.C.K.; Parks, M.; Weitemier, K.; Fishbein, M.; Cronn, R.C.; Liston, A. Navigating the tip of the genomic iceberg: Next-generation sequencing for plant systematics. Am. J. Bot. 2012, 99, 349–364. [Google Scholar] [CrossRef] [Green Version]
  10. Dodsworth, S. Genome skimming for next-generation biodiversity analysis. Trends Plant Sci. 2015, 20, 525–527. [Google Scholar] [CrossRef]
  11. Rachtman, E.; Balaban, M.; Bafna, V.; Mirarab, S. The impact of contaminants on the accuracy of genome skimming and the effectiveness of exclusion read filters. Mol. Ecol. Resour. 2020, 20, 649–661. [Google Scholar] [CrossRef] [PubMed]
  12. Yang, Y.; Li, Q.; Kong, L.; Yu, H. Mitogenomic phylogeny of Nassarius (Gastropoda: Neogastropoda). Zool. Scr. 2019, 48, 302–312. [Google Scholar] [CrossRef]
  13. Uribe, J.E.; Williams, S.T.; Templado, J.; Abalde, S.; Zardoya, R. Denser mitogenomic sampling improves resolution of the phylogeny of the superfamily Trochoidea (Gastropoda: Vetigastropoda). J. Molluscan Stud. 2017, 83, 111–118. [Google Scholar] [CrossRef]
  14. Irisarri, I.; Uribe, J.E.; Eernisse, D.J.; Zardoya, R. A mitogenomic phylogeny of chitons (Mollusca: Polyplacophora). BMC Evol. Biol. 2020, 20, 22. [Google Scholar] [CrossRef]
  15. Zardoya, R.; Meyer, A. Phylogenetic performance of mitochondrial protein-coding genes in resolving relationships among vertebrates. Mol. Biol. Evol. 1996, 13, 933–942. [Google Scholar] [CrossRef] [Green Version]
  16. Uribe, J.E.; Irisarri, I.; Templado, J.; Zardoya, R. New patellogastropod mitogenomes help counteracting long-branch attraction in the deep phylogeny of gastropod mollusks. Mol. Phylogenetics Evol. 2019, 133, 12–23. [Google Scholar] [CrossRef]
  17. Sarmashghi, S.; Bohmann, K.; Gilbert, M.T.P.; Bafna, V.; Mirarab, S. Skmer: Assembly-free and alignment-free sample identification using genome skims. Genome Biol. 2019, 20, 34. [Google Scholar] [CrossRef] [Green Version]
  18. Sohn, J.I.; Nam, J.W. The present and future of de novo whole-genome assembly. Brief Bioinform. 2018, 19, 23–40. [Google Scholar] [CrossRef] [Green Version]
  19. Benoit, G.; Peterlongo, P.; Mariadassou, M.; Drezen, E.; Schbath, S.; Lavenier, D.; Lemaitre, C. Multiple comparative metagenomics using multiset k-mer counting. PeerJ Comput. Sci. 2016, 2, e94. [Google Scholar] [CrossRef] [Green Version]
  20. Fan, H.; Ives, A.R.; Surget-Groba, Y.; Cannon, C.H. An assembly and alignment-free method of phylogeny reconstruction from next-generation sequencing data. BMC Genom. 2015, 16, 522. [Google Scholar] [CrossRef] [Green Version]
  21. Ondov, B.D.; Treangen, T.J.; Melsted, P.; Mallonee, A.B.; Bergman, N.H.; Koren, S.; Phillippy, A.M. Mash: Fast genome and metagenome distance estimation using MinHash. Genome Biol. 2016, 17, 132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Balaban, M.; Sarmashghi, S.; Mirarab, S. APPLES: Scalable distance-based phylogenetic placement with or without alignments. Syst. Biol. 2020, 69, 566–578. [Google Scholar] [CrossRef] [PubMed]
  23. Nakano, T.; Ozawa, T. Worldwide phylogeography of limpets of the order Patellogastropoda: Molecular, morphological and palaeontological evidence. J. Molluscan Stud. 2007, 73, 79–99. [Google Scholar] [CrossRef]
  24. Nakano, T.; Sasaki, T. Recent advances in molecular phylogeny, systematics and evolution of patellogastropod limpets. J. Molluscan Stud. 2011, 77, 203–217. [Google Scholar] [CrossRef]
  25. Bushnell, B. BBTools Software Package. Available online: http://sourceforge.Net/projects/bbmap (accessed on 1 December 2021).
  26. Chen, S.F.; Zhou, Y.Q.; Chen, Y.R.; Gu, J. Fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  27. Wood, D.E.; Lu, J.; Langmead, B. Improved metagenomic analysis with kraken 2. Genome Biol. 2019, 20, 257. [Google Scholar] [CrossRef] [Green Version]
  28. Marçais, G.; Kingsford, C. A fast, lock-free approach for efficient parallel counting of occurrences of k-mers. Bioinformatics 2011, 27, 764–770. [Google Scholar] [CrossRef] [Green Version]
  29. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [Green Version]
  30. Lefort, V.; Desper, R.; Gascuel, O. FastME 2.0: A comprehensive, accurate, and fast distance-based phylogeny inference program. Mol. Biol. Evol. 2015, 32, 2798–2800. [Google Scholar] [CrossRef] [Green Version]
  31. Moshiri, N. TreeSwift: A massively scalable Python tree package. SoftwareX 2020, 11, 100436. [Google Scholar] [CrossRef]
  32. Jukes, T.H.; Cantor, C.R. Evolution of protein molecules. In Mammalian Protein Metabolism; Munro, H.N., Ed.; Academic Press: New York, NY, USA, 1969; Volume III, pp. 21–132. ISBN 978-1-4832-3211-9. [Google Scholar]
  33. Rambaut, A. FigTree v1.4.4. Available online: http://tree.bio.ed.ac.uk/software/figtree (accessed on 15 December 2021).
  34. Xu, T.; Qi, L.; Kong, L.F.; Li, Q. Mitogenomics reveals phylogenetic relationships of Patellogastropoda (Mollusca, Gastropoda) and dynamic gene rearrangements. Zool. Scr. 2022, 51, 147–160. [Google Scholar] [CrossRef]
  35. Hickerson, M.J.; Meyer, C.P.; Moritz, C.; Hedin, M. DNA Barcoding will often fail to discover new animal species over broad parameter space. Syst. Biol. 2006, 55, 729–739. [Google Scholar] [CrossRef] [PubMed]
  36. Besnard, G.; Christin, P.A.; Male, P.-J.G.; Lhuillier, E.; Lauzeral, C.; Coissac, E.; Vorontsova, M.S. From museums to genomics: Old herbarium specimens shed light on a C3 to C4 transition. J. Exp. Bot. 2014, 65, 6711–6721. [Google Scholar] [CrossRef] [PubMed]
  37. Liu, S.; Wang, X.; Xie, L.; Tan, M.; Li, Z.; Su, X.; Zhang, H.; Misof, B.; Kjer, K.M.; Tang, M.; et al. Mitochondrial capture enriches mito-DNA 100-fold, enabling PCR-free mitogenomics biodiversity analysis. Mol. Ecol. Resour. 2016, 16, 470–479. [Google Scholar] [CrossRef]
  38. Zapata, F.; Wilson, N.G.; Howison, M.; Andrade, S.C.S.; Jörger, K.M.; Schroedl, M.; Goetz, F.E.; Giribet, G.; Dunn, C.W. Phylogenomic analyses of deep gastropod relationships reject Orthogastropoda. Proc. R Soc. B 2014, 281, 20141739. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Gaitán-Espitia, J.D.; González-Wevar, C.A.; Poulin, E.; Cardenas, L. Antarctic and sub-Antarctic Nacella limpets reveal novel evolutionary characteristics of mitochondrial genomes in Patellogastropoda. Mol. Phylogenetics Evol. 2019, 131, 1–7. [Google Scholar] [CrossRef] [PubMed]
Genes 13 01192 g001 550
Figure 1. Coverage distribution of P. saccharina lanx (Psl), P. conulus (Pco), P. ryukyuensis (Pry), C. toreuma (GD) (Cto-GD), C. toreuma (HN) (Cto-HN), C. grata (Cgr), C. nigrolineata (Cni), L. goshimai (Lgo), L. cassis (Lca), N. radula (Nra), and S. flexuosa (Sfl) under different sizes of genome skims.
Figure 1. Coverage distribution of P. saccharina lanx (Psl), P. conulus (Pco), P. ryukyuensis (Pry), C. toreuma (GD) (Cto-GD), C. toreuma (HN) (Cto-HN), C. grata (Cgr), C. nigrolineata (Cni), L. goshimai (Lgo), L. cassis (Lca), N. radula (Nra), and S. flexuosa (Sfl) under different sizes of genome skims.
Genes 13 01192 g001
Genes 13 01192 g002 550
Figure 2. Phylogenetic placement of C. toreuma (query species) under different sizes of genome skims.
Figure 2. Phylogenetic placement of C. toreuma (query species) under different sizes of genome skims.
Genes 13 01192 g002
Genes 13 01192 g003 550
Figure 3. Phylogenetic placement of S. flexuosa (query species) under different sizes of genome skims.
Figure 3. Phylogenetic placement of S. flexuosa (query species) under different sizes of genome skims.
Genes 13 01192 g003
Genes 13 01192 g004 550
Figure 4. Phylogenetic placement of P. conulus (query species) under different sizes of genome skims.
Figure 4. Phylogenetic placement of P. conulus (query species) under different sizes of genome skims.
Genes 13 01192 g004
Table
Table 1. List of species used in this study.
Table 1. List of species used in this study.
SubclassFamilySpeciesLocality
Patellogastropoda
Nacellidae
Cellana toreumaYangjiang, Guangdong, China
Wenchang, Hainan, China
Cellana nigrolineataJeju Island, South Korea
Cellana grataNingde, Fujian, China
Patellidae
Scutellastra flexuosaSansha, Hainan, China
Lottiidae
Nipponacmea radulaWeihai, Shandong, China
Patelloida ryukyuensisWeihai, Shandong, China
Patelloida saccharina lanxWenchang, Hainan, China
Patelloida conulusWeihai, Shandong, China
Lottia cassisWeihai, Shandong, China
Lottia goshimaiQingdao, Shandong, China
Vetigastropoda
Trochidae
Trochus maculatusSanya, Hainan, China
Table
Table 2. The calculation distance from COI and different sizes of genome skims between C. toreuma (HN) and reference species in our study. Color shows the distance ranking between reference species and query species, that is, the darker the color, the farther the relationship.
Table 2. The calculation distance from COI and different sizes of genome skims between C. toreuma (HN) and reference species in our study. Color shows the distance ranking between reference species and query species, that is, the darker the color, the farther the relationship.
COI0.1 Gb0.5 Gb1 Gb2 Gb4 GbLargest Data
C. toreuma (GD)0.0030.09180.00950.01150.01310.01450.0159
C. grata0.1740.19810.12720.13340.13590.13820.1406
C. nigrolineata0.1870.19870.12360.12990.13370.13730.1397
P. ryukyuensis0.3680.27270.19860.22610.23180.25310.2554
P. conulus0.3840.26300.19850.21060.22540.23950.2410
L. goshimai0.3920.25410.18830.19550.20420.21740.2229
P. saccharina lanx0.3990.24800.18070.20410.21120.22300.228
L. cassis0.6760.27050.20780.23110.23420.24790.2609
N. radula0.6840.24390.20160.21120.21160.23200.2389
Color from light to dark:Genes 13 01192 i001.
Table
Table 3. The calculation distance from COI and different sizes of genome skims between S. flexuosa and reference species in our study. Color shows the distance ranking between reference species and query species, that is, the darker the color, the farther the relationship.
Table 3. The calculation distance from COI and different sizes of genome skims between S. flexuosa and reference species in our study. Color shows the distance ranking between reference species and query species, that is, the darker the color, the farther the relationship.
COI0.1 Gb0.5 Gb1 Gb2 Gb4 GbLargest Data
C. grata0.2160.2541 0.19530.21380.2174 0.2359 0.2408
C. toreuma0.2160.2439 0.1697 0.1806 0.1974 0.2139 0.2230
C. nigrolineata0.2240.2582 0.1926 0.2077 0.2173 0.2365 0.2470
L. goshimai0.3760.2503 0.1806 0.1850 0.2037 0.2179 0.2229
P. ryukyuensis0.3790.27050.1905 0.2120 0.2194 0.23840.2379
P. conulus0.3880.2480 0.1885 0.2011 0.2112 0.2218 0.2299
P. saccharina lanx0.4080.2330 0.1637 0.1788 0.1881 0.2003 0.2043
L. cassis0.6540.27270.19610.21410.24070.24030.2526
N. radula0.6730.2480 0.1863 0.2030 0.22170.2317 0.2360
Color from light to dark:Genes 13 01192 i002.
Table
Table 4. The calculation distance from COI and different sizes of genome skims between P. conulus and reference species in our study. Color shows the distance ranking between reference species and query species, that is, the darker the color, the farther the relationship.
Table 4. The calculation distance from COI and different sizes of genome skims between P. conulus and reference species in our study. Color shows the distance ranking between reference species and query species, that is, the darker the color, the farther the relationship.
COI0.1 Gb0.5 Gb1 Gb2 Gb4 GbLargest Data
P. ryukyuensis0.1520.11550.05010.05700.06560.07290.0778
P. saccharina lanx0.2340.19000.13080.13750.14810.15810.1627
L. goshimai0.2450.17440.10440.11220.12160.12630.1288
L. cassis0.3060.21760.1520.15740.16660.17870.1797
N. radula0.3250.19510.13100.14340.15280.15970.1614
C. nigrolineata0.3590.27750.23530.24340.23940.24580.2524
C. grata0.3620.28020.22250.24870.24450.24350.2501
C. toreuma0.3650.26300.19860.20410.23000.24910.2467
Color from light to dark:Genes 13 01192 i003.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Xu, T.; Kong, L.; Li, Q. Testing Efficacy of Assembly-Free and Alignment-Free Methods for Species Identification Using Genome Skims, with Patellogastropoda as a Test Case. Genes 2022, 13, 1192. https://doi.org/10.3390/genes13071192

AMA Style

Xu T, Kong L, Li Q. Testing Efficacy of Assembly-Free and Alignment-Free Methods for Species Identification Using Genome Skims, with Patellogastropoda as a Test Case. Genes. 2022; 13(7):1192. https://doi.org/10.3390/genes13071192

Chicago/Turabian Style

Xu, Tao, Lingfeng Kong, and Qi Li. 2022. "Testing Efficacy of Assembly-Free and Alignment-Free Methods for Species Identification Using Genome Skims, with Patellogastropoda as a Test Case" Genes 13, no. 7: 1192. https://doi.org/10.3390/genes13071192

APA Style

Xu, T., Kong, L., & Li, Q. (2022). Testing Efficacy of Assembly-Free and Alignment-Free Methods for Species Identification Using Genome Skims, with Patellogastropoda as a Test Case. Genes, 13(7), 1192. https://doi.org/10.3390/genes13071192

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop