Premium
Phylogenomics from low‐coverage whole‐genome sequencing
Author(s) -
Zhang Feng,
Ding Yinhuan,
Zhu ChaoDong,
Zhou Xin,
Orr Michael C.,
Scheu Stefan,
Luan YunXia
Publication year - 2019
Publication title -
methods in ecology and evolution
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 3.425
H-Index - 105
ISSN - 2041-210X
DOI - 10.1111/2041-210x.13145
Subject(s) - phylogenomics , genome , biology , whole genome sequencing , phylogenetic tree , computational biology , dna sequencing , genomics , sequence assembly , metagenomics , deep sequencing , hybrid genome assembly , pipeline (software) , transcriptome , genetics , gene , computer science , clade , gene expression , programming language
Abstract Phylogenetic studies are increasingly reliant on next‐generation sequencing. Transcriptomic and hybrid enrichment sequencing techniques remain the most prevalent methods for phylogenomic data collection due to their relatively low demands for computing powers and sequencing prices, compared to whole‐genome sequencing (WGS). However, the transcriptome‐based method is constrained by the availability of fresh materials and hybrid enrichment is limited by genomic resources necessary in probe designs, especially for non‐model organisms. We present a novel WGS‐based pipeline for extracting essential phylogenomic markers through rapid de novo genome assembling from low‐coverage genome data, employing a series of computationally efficient bioinformatic tools. We tested the pipeline on a Hexapoda dataset and a more focused Phthiraptera dataset (genome sizes 0.1–2 Gbp), and further investigated the effects of sequencing depth on target assembly success rate based on the raw data of six insect genomes (0.1–1 Gbp). Each genome assembly was completed in 2–24 hr on desktop PCs. We extracted 872–1,615 near‐universal single‐copy orthologs (Benchmarking Universal Single‐Copy Orthologs [BUSCOs]) per species. This method also enables the development of ultraconserved element (UCE) probe sets; we generated probes for Phthiraptera based on our WGS assemblies, containing 55,030 baits targeting 2,832 loci, from which we extracted 2,125–2,272 UCEs. Resulting phylogenetic trees all agreed with the currently accepted topologies, indicating that markers produced in our methods were valid for phylogenomic studies. We also showed that 10–20× sequencing coverage was sufficient to produce hundreds to thousands of targeted loci from BUSCO sets, and an even lower coverage (5×) was required for UCEs. Our study demonstrates the feasibility of conducting phylogenomics from low‐coverage WGS for a wide range of organisms without reference genomes. This new approach has major advantages in data collection, particularly in reducing sequencing cost and computing consumption, while expanding loci choices.