Phylogenomics from low‐coverage whole‐genome sequencing. Issue 4 (30th January 2019)

• Record Type: Journal Article
• Authors: Zhang, Feng; Ding, Yinhuan; Zhu, Chao‐Dong; Zhou, Xin; Orr, Michael C.; Scheu, Stefan; Luan, Yun‐Xia
• Editors: Matschiner, Michael
• Is Part Of: Methods in ecology and evolution. Volume 10:Issue 4(2019)
• Page Start: 507
• Page End: 517
• Publication Date: 2019-01-30
• Subjects: desktop PC -- genome assembly -- hybrid enrichment -- single‐copy orthologs -- ultraconserved elements; Ecology -- Periodicals; Evolution -- Periodicals; 577
• DOI: 10.1111/2041-210X.13145
• Languages: English
• ISSNs: 2041-210X

Abstract

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.