Comparative performance of the BGISEQ-500 vs Illumina HiSeq2500 sequencing platforms for palaeogenomic sequencing

Abstract

Ancient DNA research has been revolutionized following development of next-generation sequencing platforms. Although a number of such platforms have been applied to ancient DNA samples, the Illumina series are the dominant choice today, mainly because of high production capacities and short read production. Recently a potentially attractive alternative platform for palaeogenomic data generation has been developed, the BGISEQ-500, whose sequence output are comparable with the Illumina series. In this study, we modified the standard BGISEQ-500 library preparation specifically for use on degraded DNA, then directly compared the sequencing performance and data quality of the BGISEQ-500 to the Illumina HiSeq2500 platform on DNA extracted from 8 historic and ancient dog and wolf samples. The data generated were largely comparable between sequencing platforms, with no statistically significant difference observed for parameters including level (P = 0.371) and average sequence length (P = 0718) of endogenous nuclear DNA, sequence GC content (P = 0.311), double-stranded DNA damage rate (v. 0.309), and sequence clonality (P = 0.093). Small significant differences were found in single-strand DNA damage rate (δS; slightly lower for the BGISEQ-500, P = 0.011) and the background rate of difference from the reference genome (θ; slightly higher for BGISEQ-500, P = 0.012). This may result from the differences in amplification cycles used to polymerase chain reaction–amplify the libraries. A significant difference was also observed in the mitochondrial DNA percentages recovered (P = 0.018), although we believe this is likely a stochastic effect relating to the extremely low levels of mitochondria that were sequenced from 3 of the samples with overall very low levels of endogenous DNA. Although we acknowledge that our analyses were limited to animal material, our observations suggest that the BGISEQ-500 holds the potential to represent a valid and potentially valuable alternative platform for palaeogenomic data generation that is worthy of future exploration by those interested in the sequencing and analysis of degraded DNA.

The first chromosome‐level genome for a marine mammal as a resource to study ecology and evolution

“…Using sequencing data from WGS libraries (200 bp, 400 bp, 2 kb, 5 kb and 10 kb), we first assembled the genome using soapdenovo (version 2.04) (Luo et al, ) with parameter settings “‐K 49 –d 4 –D 4 –R,” and gaps in scaffolds were then filled in using kgf (version 1.19) and gapcloser (version 1.10) (Luo et al, ) with default parameters. For assembly of the 10 × Genomics Chromium library data, the clean fastq files were converted so as to be recognized by 10 × genomics supernova (version 1.2.0) (Weisenfeld, Kumar, Shah, Church, & Jaffe, ) using an in‐house script.…”

MobiSeq: De novo SNP discovery in model and non‐model species through sequencing the flanking region of transposable elements

“…Size‐selected DNA was then quantified using Qubit HS reagents, as previously described, and loaded in a TapeStation High Sensitivity (TS‐HS) 2200 (Agilent) to obtain a length profile of the fragmented material. Illumina library building was based on a recently developed blunt‐end single‐tube protocol (Carøe et al, ) with some modifications as in Mak et al () (full protocol available as Supporting Information Appendix S1). In particular, the ligated adapters differed from the ones in the original protocol by excluding the use of a P7 adapter (only using the mP5 adapter) and excluding the adapter fill‐in reaction.…”

Tagsteady: a metabarcoding library preparation protocol to avoid false assignment of sequences to samples

“…End-repair without T4 DNA Polymerase: for treatments -/+ and -/-, an end-repair mastermix was made by combining 4 μl T4 DNA ligase reaction buffer (New England Biolabs, NEB, Ipswich, Massachusetts, US), 0.5 μl dATP (10mM) (Thermo-Fisher), 2 μl reaction booster mix (consisting of 25 % PEG-4000 (Sigma Aldrich, 50%), 2 mg/ml BSA (Thermo-Fisher) and 400 mM NaCl) (Mak et al 2017) , 2 μl T4 PNK (NEB, cat#M0201S, 10 U/μl) and 1.5 μl Klenow Fragment (3'->5' exo-) (NEB, cat#M0212S, 5 U/μl) per amplicon pool reaction. Ten μl of this mastermix was then added to each 30 μl amplicon pool, mixed well by pipetting and incubated for 30 minutes at 37°C followed by 30 minutes at 65°C and finally cooled to 4°C.…”