Phylogenomics from low-coverage whole-genome sequencing

View Researcher's Other Codes

Disclaimer: The provided code links for this paper are external links. Science Nest has no responsibility for the accuracy, legality or content of these links. Also, by downloading this code(s), you agree to comply with the terms of use as set out by the author(s) of the code(s).

Please contact us in case of a broken link from here

Authors Feng Zhang, Yinhuan Ding, Chao-Dong Zhu, Xin Zhou, Michael C. Orr, Stefan Scheu, Yun-Xia Luan
Journal/Conference Name Methods in Ecology and Evolution
Paper Category , ,
Paper 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.
Date of publication 2019
Code Programming Language Shell

Copyright Researcher 2022