A performance evaluation of Nextera XT and KAPA HyperPlus for rapid Illumina library preparation of long-range mitogenome amplicons

The effect of library preparation protocol on the efficiency of heteroplasmy detection in mitochondrial DNA using two massively parallel sequencing Illumina systems

Massively parallel sequencing (MPS) technology has become the gold standard in mitochondrial DNA research due to its high sensitivity in detecting mtDNA heteroplasmy, a prognostic marker in various medical applications. Various MPS technologies and platforms used for mtDNA analysis exist. Obtaining reliable and sensitive results requires deep and uniform coverage of the entire mtDNA sequence, which is heavily influenced by the choice of library preparation method and sequencing platform. Here, we present a comparison of the sequencing coverage and the ability to heteroplasmy detection using two library preparation protocols (Nextera XT DNA Library Preparation Kit and Nextera DNA Flex Library Preparation Kit) and two different (MiSeq FGx and ISeq 100) Illumina MPS platforms. Our study indicates that the Nextera DNA Flex Library protocol provides a more balanced coverage along the mitogenome and a reliable heteroplasmy detection with both MiSeq and iSeq Illumina MPS systems.

Capability of the iSeq 100 sequencing system from Illumina to detect low-level substitutions in the human mitochondrial genome

The significance of mtDNA heteroplasmy in forensic and medical genetics has increased recently because massively parallel sequencing (MPS) technologies enable more accurate and precise detection of minority nucleotide variants. Recent reports have shown that detection of low-level substitutions may depend on library preparation or sequencing protocol, and can vary for different MPS platforms. The MiSeq (Illumina) and Ion S5 (Thermo Fisher Scientific) are mainly used for heteroplasmy detection, but no data are available regarding the iSeq 100, an Illumina platform of the smallest throughput. Notably, unlike the other systems, the machine utilizes sequencing by synthesis one-channel chemistry to determine DNA sequences. Thus, it is important to verify the capability of the iSeq 100 system to determine mitochondrial haplotypes and detect heteroplasmic substitutions. In this study, previously determined entire mitochondrial genomes were sequenced with the iSeq 100 system. Each mitogenome was sequenced twice, giving approximately 2000x and 10,000x coverage. All homoplasmic mutations and minority variants above the 19 % level detected with the iSeq 100 system were also observed after dideoxy sequencing. Moreover, all heteroplasmic substitutions above the 2 % level were consistently detected with SBS one-channel chemistry. However, detection of low-level mtDNA variants may require additional, confirmatory experiments. In summary, the iSeq 100 system enables reproducible and accurate sequencing of human mitochondrial genomes. Detection of mtDNA minority variants depends on the laboratory protocol and sequencing platform used, but homoplasmic mutations and heteroplasmy above the 2 % level can be correctly detected with the iSeq 100 system.

Whole Mitochondrial Genome Detection and Analysis of Two- to Four-Generation Maternal Pedigrees Using a New Massively Parallel Sequencing Panel

Mitochondrial DNA (mtDNA) is an effective genetic marker in forensic practice, especially for aged bones and hair shafts. Detection of the whole mitochondrial genome (mtGenome) using traditional Sanger-type sequencing is laborious and time-consuming. Additionally, its ability to distinguish point heteroplasmy (PHP) and length heteroplasmy (LHP) is limited. The application of massively parallel sequencing in mtDNA detection helps researchers to study the mtGenome in-depth. The ForenSeq mtDNA Whole Genome Kit, which contains a total of 245 short amplicons, is one of the multiplex library preparation kits for the mtGenome. We used this system to detect the mtGenome in the blood samples and hair shafts of thirty-three individuals from eight two-generation pedigrees, one three-generation pedigree, and one four-generation pedigree. High-quality sequencing results were obtained. Ten unique mtGenome haplotypes were observed in the mothers from the ten pedigrees. A total of 26 PHPs were observed using the interpretation threshold of 6%. Eleven types of LHPs in six regions were evaluated in detail. When considering homoplasmic variants only, consistent mtGenome haplotypes were observed between the twice-sequenced libraries and between the blood and hair shafts from the same individual and among maternal relatives in the pedigrees. Four inherited PHPs were observed, and the remainder were de novo/disappearing PHPs in the pedigrees. Our results demonstrate the effective capability of the ForenSeq mtDNA Whole Genome Kit to generate the complete mtGenome in blood and hair shafts, as well as the complexity of mtDNA haplotype comparisons between different types of maternal relatives when heteroplasmy is considered.

Nematode mitochondrial metagenomics: A new tool for biodiversity analysis

DNA barcoding approaches have greatly increased our understanding of biodiversity on the planet, and metabarcoding is widely used for classifying members of the phylum Nematoda. However, loci typically utilized in metabarcoding studies are often unable to resolve closely related species or are unable to recover all taxa present in a sample due to inadequate PCR primer binding. Mitochondrial metagenomics (mtMG) is an alternative approach utilizing shotgun sequencing of total DNA to recover the mitochondrial genomes of all species present in samples. However, this approach requires a comprehensive reference database for identification and currently available mitochondrial sequences for nematodes are highly dominated by sequences from the order Rhabditida, and excludes many clades entirely. Here, we analysed the efficacy of mtMG for the recovery of nematode taxa and the generation of mitochondrial genomes. We first developed a curated reference database of nematode mitochondrial sequences and expanded it with 40 newly sequenced taxa. We then tested the mito-metagenomics approach using a series of nematode mock communities consisting of morphologically identified nematode species representing various feeding traits, life stages, and phylogenetic relationships. We were able to identify all but two species through the de novo assembly of COX1 genes. We were also able to recover additional mitochondrial protein coding genes (PCGs) for 23 of the 24 detected species including a full array of 12 PCGs from five of the species. We conclude that mtMG offers a potential for the effective recovery of nematode biodiversity but remains limited by the breadth of the reference database.

Low-level point heteroplasmy detection in human mitogenomes amplified with different polymerases and sequenced on MiSeq FGx platform

Introduction

Massively parallel sequencing of mitogenomes usually requires prior amplification. The PCR step may influence the quality of the data obtained, especially when low-level heteroplasmy detection is applied.

Aim

The aim of this study was to compare the reliability of two different DNA polymerases in detecting homoplasmic and heteroplasmic substitutions in human mitogenomes.

Material and methods

Mitogenomes of five samples were amplified with Long PCR Enzyme Mix from Fermentas or TaKaRa LA Taq DNA Polymerase from TaKaRa. Then, NexteraTM XT DNA libraries were sequenced on MiSeq FGx platform (Illumina). mtDNA substitutions were called for alternative variants above the 1% level.

Results

All homoplasmic substitutions detected in amplicons generated with polymerases studied here and sequenced on MiSeq FGx system were consistently identified as homoplasmies with alternative sequencing methods. TaKaRa LA Taq DNA Polymerase was found to be less accurate in low-level heteroplasmy detection than Long PCR Enzyme Mix enzyme as more false negative and false positive results were observed for minority variants called above the 1% level. Nevertheless, both PCR systems studied can be successfully used to detect authentic mtDNA substitutions, for which minority variants exceed the 3.61% level assuming at least 10,000x coverage and sequencing Nextera XT DNA libraries on MiSeq FGx machine.

Conclusions

The accuracy and sensitivity of point heteroplasmy detection with the MiSeq FGx instrument varies on polymerase used for mtDNA amplification. Therefore, it is recommended to validate the laboratory protocols used for mtDNA substitution detection prior to their implementation for the forensic or medical genetics purposes.

Argopistes sexvittatus and Argopistes capensis (Chrysomelidae: Alticini): Mitogenomics and Phylogeny of Two Flea Beetles Affecting Olive Trees

The genus Argopistes (Chrysomelidae: Alticini) is the only group of flea beetles specialized in plant hosts in the family Oleaceae. In southern Africa, Argopistes are often found feeding on African Wild Olive (Olea europaea subsp. cuspidata) and European cultivated olive (O. e. subsp. europaea), and heavy infestations can be devastating to mature trees and compromise the development of young trees. Despite their negative agricultural impact, African Argopistes are an understudied group for which no genetic data were available. We assessed the species diversity of olive flea beetles in the Western Cape province of South Africa, the largest olive-producing region in sub-Saharan Africa, by collecting adult specimens on wild and cultivated olive trees between 2015 and 2017. Argopistes sexvittatus Bryant, 1922 (n = 289) dominated at all sampling sites, and Argopistes capensis Bryant, 1944 (n = 2) was found only once. Argopistes oleae Bryant, 1922, a third species previously reported in the region, was not found. The complete mitogenomes of one A. capensis and two A. sexvittatus (striped and black morphotypes) individuals were sequenced for phylogenetic reconstruction in the context of other 64 species. The two olive flea beetle species form a monophyletic clade with other Argopistes, supporting the hypothesis that the exclusive feeding habit on Oleaceae is an evolutionary adaptation in this genus.

Optimization of enzymatic fragmentation is crucial to maximize genome coverage: a comparison of library preparation methods for Illumina sequencing

Background

Novel commercial kits for whole genome library preparation for next-generation sequencing on Illumina platforms promise shorter workflows, lower inputs and cost savings. Time savings are achieved by employing enzymatic DNA fragmentation and by combining end-repair and tailing reactions. Fewer cleanup steps also allow greater DNA input flexibility (1 ng-1 μg), PCR-free options from 100 ng DNA, and lower price as compared to the well-established sonication and tagmentation-based DNA library preparation kits.

Results

We compared the performance of four enzymatic fragmentation-based DNA library preparation kits (from New England Biolabs, Roche, Swift Biosciences and Quantabio) to a tagmentation-based kit (Illumina) using low input DNA amounts (10 ng) and PCR-free reactions with 100 ng DNA. With four technical replicates of each input amount and kit, we compared the kits’ fragmentation sequence-bias as well as performance parameters such as sequence coverage and the clinically relevant detection of single nucleotide and indel variants. While all kits produced high quality sequence data and demonstrated similar performance, several enzymatic fragmentation methods produced library insert sizes which deviated from those intended. Libraries with longer insert lengths performed better in terms of coverage, SNV and indel detection. Lower performance of shorter-insert libraries could be explained by loss of sequence coverage to overlapping paired-end reads, exacerbated by the preferential sequencing of shorter fragments on Illumina sequencers. We also observed that libraries prepared with minimal or no PCR performed best with regard to indel detection.

Conclusions

The enzymatic fragmentation-based DNA library preparation kits from NEB, Roche, Swift and Quantabio are good alternatives to the tagmentation based Nextera DNA flex kit from Illumina, offering reproducible results using flexible DNA inputs, quick workflows and lower prices. Libraries with insert DNA fragments longer than the cumulative sum of both read lengths avoid read overlap, thus produce more informative data that leads to strongly improved genome coverage and consequently also increased sensitivity and precision of SNP and indel detection. In order to best utilize such enzymatic fragmentation reagents, researchers should be prepared to invest time to optimize fragmentation conditions for their particular samples.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08316-y.

Selective Chemical Labeling and Sequencing of 5-Hydroxymethylcytosine in DNA at Single-Base Resolution

5-Hydroxymethylcytosine (5hmC), the oxidative product of 5-methylcytosine (5mC) catalyzed by ten-eleven translocation enzymes, plays an important role in many biological processes as an epigenetic mediator. Prior studies have shown that 5hmC can be selectively labeled with chemically modified glucose moieties and enriched using click chemistry with biotin affinity approaches. Besides, DNA deaminases of the AID/APOBEC family can discriminate modified 5hmC bases from cytosine (C) or 5mC. Herein, we developed a method based on embryonic stem cell (ESC) whole-genome analysis, which could enrich 5hmC-containing DNA by selective chemical labeling and locate 5hmC sites at single-base resolution with enzyme-based deamination. The combination experimental design is an extension of previous methods, and we hope that this cost-effective single-base resolution 5hmC sequencing method could be used to promote the mechanism and diagnosis research of 5hmC.

Species Diversity and Phylogenetic Relationships of Olive Lace Bugs (Hemiptera: Tingidae) Found in South Africa

Olive lace bugs (Hemiptera: Tingidae) are small sap-sucking insects that feed on wild and cultivated Olea europaea. The diversity of olive lace bug species in South Africa, the most important olive producer on the continent, has been incompletely surveyed. Adult specimens were collected in the Western Cape province for morphological and DNA-based species identification, and sequencing of complete mitogenomes. Cysteochila lineata, Plerochila australis, Neoplerochila paliatseasi and Neoplerochila sp. were found at 12 sites. Intra- and interspecific genetic divergences and phylogenetic clustering in 30 species in 18 genera of Tingidae using new and publicly available DNA barcodes showed high levels of congruity between taxonomic and genetic data. The phylogenetic position of the four species found in South Africa was inferred using new and available mitogenomes of Tingidae. Notably, olive lace bugs formed a cluster of closely related species. However, Cysteochila was non-monophyletic as C. lineata was recovered as a sister species to P. australis whereas Cysteochila chiniana, the other representative of the genus, was grouped with Trachypeplus jacobsoni and Tingis cardui in a different cluster. This result suggests that feeding on O. europaea may have a common origin in Tingidae and warrants future research on potential evolutionary adaptations of olive lace bugs to this plant host.

Benchmarking and optimization of a high-throughput sequencing based method for transgene sequence variant analysis in biotherapeutic cell line development

In recent years, High-Throughput Sequencing (HTS) based methods to detect mutations in biotherapeutic transgene products have become a key quality step deployed during the development of manufacturing cell line clones. Previously we reported on a higher throughput, rapid mutation detection method based on amplicon sequencing (targeting transgene RNA) and detailed its implementation to facilitate cell line clone selection. By gaining experience with our assay in a diverse set of cell line development programs, we improved the computational analysis as well as experimental protocols. Here we report on these improvements as well as on a comprehensive benchmarking of our assay. We evaluated assay performance by mixing amplicon samples of a verified mutated antibody clone with a non-mutated antibody clone to generate spike-in mutations from ∼60% down to ∼0.3% frequencies. We subsequently tested the effect of 16 different sample and HTS library preparation protocols on the assay's ability to quantify mutations and on the occurrence of false-positive background error mutations (artifacts). Our evaluation confirmed assay robustness, established a high confidence limit of detection of ∼0.6%, and identified protocols that reduce error levels thereby significantly reducing a source of false positives that bottlenecked the identification of low-level true mutations.

Human Mitochondrial Control Region and mtGenome: Design and Forensic Validation of NGS Multiplexes, Sequencing and Analytical Software

Forensic mitochondrial DNA (mtDNA) analysis conducted using next-generation sequencing (NGS), also known as massively parallel sequencing (MPS), as compared to Sanger-type sequencing brings modern advantages, such as deep coverage per base (herein referred to as read depth per base pair (bp)), simultaneous sequencing of multiple samples (libraries) and increased operational efficiencies. This report describes the design and developmental validation, according to forensic quality assurance standards, of end-to-end workflows for two multiplexes, comprised of ForenSeq mtDNA control region and mtDNA whole-genome kits the MiSeq FGx^TM instrument and ForenSeq universal analysis software (UAS) 2.0/2.1. Polymerase chain reaction (PCR) enrichment and a tiled amplicon approach target small, overlapping amplicons (60–150 bp and 60–209 bp for the control region and mtGenome, respectively). The system provides convenient access to data files that can be used outside of the UAS if desired. Studies assessed a range of environmental and situational variables, including but not limited to buccal samples, rootless hairs, dental and skeletal remains, concordance of control region typing between the two multiplexes and as compared to orthogonal data, assorted sensitivity studies, two-person DNA mixtures and PCR-based performance testing. Limitations of the system and implementation considerations are discussed. Data indicated that the two mtDNA multiplexes, MiSeq FGx and ForenSeq software, meet or exceed forensic DNA quality assurance (QA) guidelines with robust, reproducible performance on samples of various quantities and qualities.

Capture enrichment and massively parallel sequencing for human identification

In the past decade, hybridization capture has gained attention within the forensic field for its possible use in human identification. One of the primary benefits to capture enrichment is its applicability to degraded DNA fragments that, due to their reduced size, are not amenable to traditional PCR enrichment techniques. Hybridization capture is typically introduced after genomic library preparation of extracted DNA templates for the subsequent enrichment of mitochondrial DNA or single nucleotide polymorphisms within the nuclear genome. The enriched molecules are then subjected to massively parallel sequencing (MPS) for sensitive and high-throughput DNA sequence generation. Bioinformatic analysis of capture product removes PCR duplicates that were introduced during the laboratory workflow in order to characterize the original DNA template molecules. In the case of aged and degraded skeletal remains, the fraction of endogenous human DNA may be very low; therefore low-coverage sequence analysis may be required. This review contains an overview of current capture methodologies and the primary literature on hybridization capture as evaluated for forensic applications.

Assessment of Illumina® Human mtDNA Genome assay: workflow evaluation with development of analysis and interpretation guidelines

Mitochondrial DNA (mtDNA) is a small but significant part of the human genome, whose applicability potential has gradually increased with the advent of massively parallel sequencing (MPS) technology. Knowledge of the particular workflow, equipment, and reagents used, along with extensive usage of negative controls to monitor all preparation steps constitute the prerequisites for confident reporting of results. In this study, we performed an assessment of Illumina® Human mtDNA Genome assay on MiSeq FGx™ instrument. Through analysis of several types of negative controls, as well as mtDNA positive controls, we established thresholds for data analysis and interpretation, consisting of several components: minimum read depth (220 reads), minimum quality score (41), percentage of minor allele sufficient for analysis (3.0%), percentage of minor allele sufficient for interpretation (6.0%), and percentage of major allele sufficient for homoplasmic variant call (97.0%). Based on these criteria, we defined internal guidelines for analysis and interpretation of mtDNA results obtained by MPS. Our study shows that the whole mtDNA assay on MiSeq FGx™ produces repeatable and reproducible results, independent of the analyst, which are also concordant with Sanger-type sequencing results for mtDNA control region, as well as with MPS results produced by NextSeq®. Overall, established thresholds and interpretation guidelines were successfully applied for the sequencing of complete mitochondrial genomes from high-quality samples. The underlying principles and proposed methodology on the definition of internal laboratory guidelines for analysis and interpretation of MPS results may be applicable to similar MPS workflows, e.g. targeting good-quality samples in forensic genetics and molecular diagnostics.

Platinum-Quality Mitogenome Haplotypes from United States Populations

A total of 1327 platinum-quality mitochondrial DNA haplotypes from United States (U.S.) populations were generated using a robust, semi-automated next-generation sequencing (NGS) workflow with rigorous quality control (QC). The laboratory workflow involved long-range PCR to minimize the co-amplification of nuclear mitochondrial DNA segments (NUMTs), PCR-free library preparation to reduce amplification bias, and high-coverage Illumina MiSeq sequencing to produce an average per-sample read depth of 1000 × for low-frequency (5%) variant detection. Point heteroplasmies below 10% frequency were confirmed through replicate amplification, and length heteroplasmy was quantitatively assessed using a custom read count analysis tool. Data analysis involved a redundant, dual-analyst review to minimize errors in haplotype reporting with additional QC checks performed by EMPOP. Applying these methods, eight sample sets were processed from five U.S. metapopulations (African American, Caucasian, Hispanic, Asian American, and Native American) corresponding to self-reported identity at the time of sample collection. Population analyses (e.g., haplotype frequencies, random match probabilities, and genetic distance estimates) were performed to evaluate the eight datasets, with over 95% of haplotypes unique per dataset. The platinum-quality mitogenome haplotypes presented in this study will enable forensic statistical calculations and thereby support the usage of mitogenome sequencing in forensic laboratories.

Characterization and mitigation of fragmentation enzyme-induced dual stranded artifacts

High-throughput short-read sequencing relies on fragmented DNA for optimal sampling of input nucleic acid. Several vendors now offer proprietary enzyme cocktails as a cheaper and more streamlined method of fragmentation when compared to acoustic shearing. We have discovered that these enzymes induce the formation of library molecules containing regions of nearby DNA from opposite strands. Sequencing reads derived from these molecules can lead to artifact-derived variant calls appearing at variant allele frequencies <5%. We present Fragmentation Artifact Detection and Elimination (FADE), software to remove these artifacts from mapped reads and mitigate artifact-related effects on downstream analysis. We find that the artifacts principally affect downstream analyses that are sensitive to a 1-3% artifact bias in the sequencing reads, such as targeted resequencing and rare variant discovery.

Next generation sequencing of STR artifacts produced from historical bone samples

STR artifacts are commonly observed in electrophoretic data and can complicate interpretation of the profiles produced. Even when a consensus approach is applied, reproducible artifacts have the potential to convolute the analysis. DNA obtained from historical bone samples is often heavily degraded and damaged, requiring the use of more sensitive procedures to increase allele recovery. Additionally, skeletal remains exposed to environmental conditions may be afflicted with microbial DNA contamination that cross-reacts with the primers during short tandem repeat (STR) multiplex amplification. STR artifacts manifested as a result of these circumstances can be sourced and characterized using new sequencing technologies to potentially ease the analysis burden. For this study, PCR product from 17 low-quality bone samples exhibiting reproducible autosomal and Y-chromosomal STR (Y-STR) artifacts in capillary electrophoresis (CE) data were sequenced with next generation sequencing (NGS). Sequenced reads were bioinformatically sorted using STRait Razor to determine the authenticity of alleles and confirm the profile generated by CE. Sequence data from the PCR products and a subset of the associated extracts were further analyzed with Kaiju to classify the microbial species present and identify potential sources of artifact peaks. A suspected Y-STR artifact was similar in sequence to Pseudomonas sp. BAY1663, a species ubiquitously found in soil. Regions of homology were observed between the Pseudomonas genome and the presumed primer binding locations for Y-STRs included in the AmpFlSTR Y-Filer STR kit. Characterization of such supposed artifact peaks may aid in interpretation of CE data and ultimately lead to increased confidence in the reported results.

Mitogenome phylogenetics in the genus Palaemon (Crustacea: Decapoda) sheds light on species crypticism in the rockpool shrimp P. elegans

The genus Palaemon comprises worldwide marine and freshwater shrimps and prawns, and some of them are ecologically or commercially important species. Palaemon is not currently a monophyletic group, so phylogenetics and systematics are constantly changing. Species crypticism has been pointed out in several Palaemon species, being the clearest evidence in the European rockpool shrimp P. elegans. Here we sequenced and described seven European Palaemon mitochondrial genomes. The mitochondrial protein-coding genes were used, along with those of three other Palaemon species, to perform mitogenome phylogenetic analyses to clarify the evolutionary relationships within the genus, and particularly to shed light on the cryptic species found within P. elegans. The Messinian Salinity Crisis (5.3-5.9 Ma, late Miocene) was proposed to be the origin of this cryptic species and it was used as aged constraint for calibration analysis. We provide the largest and the first time-calibrated mitogenome phylogeny of the genus Palaemon and mitogenome substitution rate was estimated (1.59% per million years) in Decapoda for the first time. Our results highlighted the need for future systematics changes in Palaemon and crypticism in P. elegans was confirmed. Mitochondrial genome and cox1 (1.41%) substitution rate estimates matched those published elsewhere, arguing that the Messinian Salinity Crisis was a plausible event driving the split between P. elegans and its cryptic species. Molecular dating suggested that Pleistocene glaciations were likely involved in the differentiation between the Atlantic and Mediterranean populations of P. elegans. On the contrary, the divergence between the Atlantic and Mediterranean populations of the common littoral shrimp P. serratus was greater and dated to be much older (4.5-12.3 Ma, Plio-Miocene), so we considered that they could represent two separated species. Therefore, species crypticism in the genus Palaemon seems to be a common phenomenon.

A universal and independent synthetic DNA ladder for the quantitative measurement of genomic features

Standard units of measurement are required for the quantitative description of nature; however, few standard units have been established for genomics to date. Here, we have developed a synthetic DNA ladder that defines a quantitative standard unit that can measure DNA sequence abundance within a next-generation sequencing library. The ladder can be spiked into a DNA sample, and act as an internal scale that measures quantitative genetics features. Unlike previous spike-ins, the ladder is encoded within a single molecule, and can be equivalently and independently synthesized by different laboratories. We show how the ladder can measure diverse quantitative features, including human genetic variation and microbial abundance, and also estimate uncertainty due to technical variation and improve normalization between libraries. This ladder provides an independent quantitative unit that can be used with any organism, application or technology, thereby providing a common metric by which genomes can be measured.

Comparison of Illumina MiSeq and the Ion Torrent PGM and S5 platforms for whole-genome sequencing of picornaviruses and caliciviruses

Next-generation sequencing is a powerful tool for virological surveillance. While Illumina® and Ion Torrent® sequencing platforms are used extensively for generating viral RNA genome sequences, there is limited data comparing different platforms. The Illumina MiSeq, Ion Torrent PGM and Ion Torrent S5 platforms were evaluated using a panel of sixteen specimens containing picornaviruses and human caliciviruses (noroviruses and sapoviruses). The specimens were processed, using combinations of three library preparation and five sequencing kits, to assess the quality and completeness of assembled viral genomes, and an estimation of cost per sample to generate the data was calculated. The choice of library preparation kit and sequencing platform was found to impact the breadth of genome coverage and accuracy of consensus viral genomes. The Ion Torrent S5 510 chip runs produced more reads at a lower cost per sample than the highest output Ion Torrent PGM 318 chip run, and generated the highest proportion of reads for enterovirus D68 samples. However, indels at homopolymer regions impacted the accuracy of consensus genome sequences. For lower throughput sequencing runs (i.e., Ion Torrent 510 and Illumina MiSeq Nano V2), the cost per sample was lower on the MiSeq platform, whereas with higher throughput runs (Ion Torrent 530 and Illumina MiSeq V2) there is less of a difference in the cost per sample between the two sequencing platforms ($5.47-$10.25 more per sample for an Ion Torrent 530 chip run when multiplexing 24 samples). These findings suggest that the Ion Torrent S5 and Illumina MiSeq platforms are both viable options for genomic sequencing of RNA viruses, each with specific advantages and tradeoffs.

Genome sequence data of Mangrovimonas sp. strain CR14 isolated from mangrove forest at Tanjung Piai National Park, Malaysia

Mangrovimonas sp. strain CR14 is a halophilic bacterium affiliated with family Flavobacteriaceae which was successfully isolated from mangrove soil samples obtained from Tanjung Piai National Park, Johor. The whole genome of strain CR14 was sequenced on an Illumina HiSeq 2500 platform (2 × 150 bp paired end). Herein, we report the genome sequence of Mangrovimonas sp. strain CR14 in which its assembled genome consisted 20 contigs with a total size of 3,590,195 bp, 3209 coding sequences, and an average 36.08% G + C content. Genome annotation and gene mining revealed that this bacterium demonstrated proteolytic activity which could be potentially applied in detergent industry. This whole-genome shotgun data of Mangrovimonas sp. strain CR14 has been deposited at DDBJ/ENA/GenBank under the accession JAAFZY000000000. The version described in this paper is version JAAFZY010000000.

Mechanisms of nuclear content loading to exosomes

Exosome cargoes are highly varied and include proteins, small RNAs, and genomic DNA (gDNA). The presence of gDNA suggests that different intracellular compartments contribute to exosome loading, resulting in distinct exosome subpopulations. However, the loading of gDNA and other nuclear contents into exosomes (nExo) remains poorly understood. Here, we identify the relationship between cancer cell micronuclei (MN), which are markers of genomic instability, and nExo formation. Imaging flow cytometry analyses reveal that 10% of exosomes derived from cancer cells and <1% of exosomes derived from blood and ascites from patients with ovarian cancer carry nuclear contents. Treatment with genotoxic drugs resulted in increased MN and nExos both in vitro and in vivo. We observed that multivesicular body precursors and exosomal markers, such as the tetraspanins, directly interact with MN. Collectively, this work provides new insights related to nExos, which have implications for cancer biomarker development.

Determination of materials present in skeletonized human remains and the associated DNA: Development of a GC/MS protocol

DNA testing of skeletonized human remains is considered to be challenging, especially when the remains have been exposed to inhibitory materials during decomposition. Inhibitors affect the processing of DNA, either by preventing efficient extraction or interfering with down-stream PCR-based processes. Limited studies have been performed on real-world samples that have been exposed to such inhibitors. This paper presents the development of a gas chromatography/mass spectrometry (GC/MS) protocol for the evaluation of materials present in skeletonized human remains recovered from the field, as well as the DNA extracted from the same materials. Twenty-one bone samples and seventeen DNA extracts were evaluated across three solvents and multiple GC/MS parameters to determine the optimal conditions for the recovery of trace materials present. The aim of this work is to provide a technique that can determine the presence of inhibitors prior to DNA extraction, allowing analysts to optimize removal of inhibitory materials.

Development of a triplex mtDNA qPCR assay to assess quantification, degradation, inhibition, and amplification target copy numbers

A hybrid absolute/relative qPCR assay which provides information regarding the condition of mitochondrial DNA (mtDNA) in a DNA sample is described. MtDNA concentration (copy number/μL) is determined via absolute quantification using a standard curve of a synthetic duplex DNA previously described (Kavlick et al., 2011). The state of mtDNA degradation is determined via the relative quantification of a mtDNA target found within the 16 s rRNA gene which is 3× longer than that of the short target in the former duplex assay, using the delta, delta Ct (ΔΔCt) method. The presence or absence of PCR inhibitors in the sample is qualitatively determined using a custom internal positive control (IPC) system which targets a unique and non-naturally occurring duplex DNA sequence. This IPC effectively detected inhibition by humic acid, tannic acid, melanin, and EDTA. All three assay components utilize sensitive and specific hydrolysis probes. The utility of ΔΔCt method was demonstrated in a series of experiments involving laboratory-fragmented DNA. Also described is a method for estimating copy number of any mtDNA target longer than the two targets amplified. The described triplex assay works well for intact and for fragmented or degraded mtDNA and therefore may be useful in forensic and ancient DNA disciplines as well as in biomedical research or practice.

A Microfluidics Workflow for Sample Preparation for Next-Generation DNA Sequencing

Next-generation sequencing technology requires amplified, short DNA fragments with known end sequences. Samples must undergo processing steps, including extraction and purification of genomic DNA (gDNA), fragmentation, end repair, adapter ligation, and amplification, to prepare a sequencing library. The process of sample preparation requires careful control of temperature and buffer conditions, as well as the stringent removal of contaminants. As a result, library preparation methods are often plagued by sample loss, long protocol times, numerous manual steps, and high cost. We attempt to understand and optimize each step of sample preparation on a microfluidic platform using magnetic bead motion through channels containing immiscible phases. Our platform integrates all steps associated with library preparation with no buffer exchanges and utilizes just 30-60 µL of reagents. Our chip shows a sixfold improvement in yield compared with an affinity spin column when capturing gDNA from samples of ~50 ± 4 MCF-7 cells. Finally, we show whole-genome shotgun sequencing results from 660 pg of human gDNA, in which >93 ± 1% of reads map to a reference genome at or above 99.9% confidence, matching a commercially available sample preparation kit optimized for low-cell-count samples.

Developmental validation of a Nextera XT mitogenome Illumina MiSeq sequencing method for high-quality samples

Generating mitochondrial genome (mitogenome) data from reference samples in a rapid and efficient manner is critical to harnessing the greater power of discrimination of the entire mitochondrial DNA (mtDNA) marker. The method of long-range target enrichment, Nextera XT library preparation, and Illumina sequencing on the MiSeq is a well-established technique for generating mitogenome data from high-quality samples. To this end, a validation was conducted for this mitogenome method processing up to 24 samples simultaneously along with analysis in the CLC Genomics Workbench and utilizing the AQME (AFDIL-QIAGEN mtDNA Expert) tool to generate forensic profiles. This validation followed the Federal Bureau of Investigation's Quality Assurance Standards (QAS) for forensic DNA testing laboratories and the Scientific Working Group on DNA Analysis Methods (SWGDAM) validation guidelines. The evaluation of control DNA, non-probative samples, blank controls, mixtures, and nonhuman samples demonstrated the validity of this method. Specifically, the sensitivity was established at ≥25 pg of nuclear DNA input for accurate mitogenome profile generation. Unreproducible low-level variants were observed in samples with low amplicon yields. Further, variant quality was shown to be a useful metric for identifying sequencing error and crosstalk. Success of this method was demonstrated with a variety of reference sample substrates and extract types. These studies further demonstrate the advantages of using NGS techniques by highlighting the quantitative nature of heteroplasmy detection. The results presented herein from more than 175 samples processed in ten sequencing runs, show this mitogenome sequencing method and analysis strategy to be valid for the generation of reference data.