A Comparison and Integration of MiSeq and MinION Platforms for Sequencing Single Source and Mixed Mitochondrial Genomes

Forensic use of human mitochondrial DNA: A review

In forensics, genetic human identification is generally achieved by nuclear STR DNA typing. However, forensic samples often yield DNA in exiguous quantity and low quality, impairing the generation of conclusive DNA profiles by STR typing. In such cases, mitochondrial DNA (mtDNA) can be used as an alternative solution in forensic human identification. The high copy number, small circular DNA, high mutation rate, maternal inheritance, and absence of recombination are mtDNA's key features in forensics. In this work, we review mtDNA characteristics, forensic applications, sequencing methodologies and present some relevant examples in the forensic science literature.

Research progress and application of the third-generation sequencing technologies in forensic medicine

Third-generation sequencing technologies, exemplified by single-molecule real-time sequencing and nanopore sequencing, provide a constellation of advantages, including long read lengths, high throughput, real-time sequencing capabilities, and remarkable portability. These cutting-edge methodologies have provided new tools for genomic analysis in forensic medicine. To gain a comprehensive understanding of the current applications and cutting-edge trends of third-generation sequencing technologies in forensic medicine, this study retrieved relevant literature from the China National Knowledge Infrastructure (CNKI) database and the Web of Science (WOS) database. Using bibliometric software CiteSpace 6.1.R6, the study visualized publication volume, countries, and keywords related to the application of third-generation sequencing technologies in forensic medicine from 2014 to 2023. The review then summarized the foundational principles, characteristics, and promising prospects of third-generation sequencing technologies in forensic medicine. Notably, it highlights their remarkable contributions in forensic individual identification, body fluid identification, forensic epigenetic analysis, microbial analysis and forensic species identification.

A proof-of-principle study: The potential application of MiniHap biomarkers in ancestry inference based on the QNome nanopore sequencing

Haplotyped SNPs convey forensic-related information, and microhaplotypes (MHs), as the most representative of this kind of marker, have proved the potential value for human forensics. In recent years, nanopore sequencing technology has developed rapidly, with its outstanding ability to sequence long continuous DNA fragments and obtain phase information, making the detection of longer haplotype marker possible. In this proof-of-principle study, we proposed a new type of forensic marker, MiniHap, based on five or more SNPs within a molecular distance less than 800 bp, and investigated the haplotype data of 56 selected MiniHaps in five Chinese populations using the QNome nanopore sequencing. The sequencing performance, allele (haplotype) frequencies, forensic parameters, effective number of alleles (Ae), and informativeness (In) were subsequently calculated. In addition, we performed principal component analysis (PCA), phylogenetic tree, and structure analysis to investigate the population genetic relationships and ancestry components among the five investigated populations and 26 worldwide populations. MiniHap-04 exhibited remarkable forensic efficacy, with 148 haplotypes reported and the Ae was 66.9268. In addition, the power of discrimination (PD) was 0.9934, the probability of exclusion (PE) was 0.9898, and the In value was 0.7893. Of the 56 loci, 85.71% had PD values above 0.85, 66.07% had PE values above 0.54, 67.86% had Ae values over 7.0%, and 55.36% were with In values above 0.2 across all samples, indicating that most of the MiniHaps are suitable for individual identification, paternity testing, mixture deconvolution, and ancestry inference. Moreover, the results of PCA, phylogenetic tree and structure analysis demonstrated that this MiniHap panel had the competency in continental population ancestry inference, but the differentiation within intracontinental/linguistically restricted subpopulations was not ideal. Such findings suggested that the QNome device for MiniHap detection was feasible and this novel marker has the potential in ancestry inference. Yet, the establishment of a more comprehensive database with sufficient reference population data remains necessary to screen more suitable MiniHaps.

Genomic diversity, pathogenicity and antimicrobial resistance of Escherichia coli isolated from poultry in the southern United States

Escherichia coli (E. coli) are typically present as commensal bacteria in the gastro-intestinal tract of most animals including poultry species, but some avian pathogenic E. coli (APEC) strains can cause localized and even systematic infections in domestic poultry. Emergence and re-emergence of antimicrobial resistant isolates (AMR) constrain antibiotics usage in poultry production, and development of an effective vaccination program remains one of the primary options in E. coli disease prevention and control for domestic poultry. Thus, understanding genetic and pathogenic diversity of the enzootic E. coli isolates, particularly APEC, in poultry farms is the key to designing an optimal vaccine candidate and to developing an effective vaccination program. This study explored the genomic and pathogenic diversity among E. coli isolates in southern United States poultry. A total of nine isolates were recovered from sick broilers from Mississippi, and one from Georgia, with epidemiological variations among clinical signs, type of housing, and bird age. The genomes of these isolates were sequenced by using both Illumina short-reads and Oxford Nanopore long-reads, and our comparative analyses suggested data from both platforms were highly consistent. The 16 s rRNA based phylogenetic analyses showed that the 10 bacteria strains are genetically closer to each other than those in the public database. However, whole genome analyses showed that these 10 isolates encoded a diverse set of reported virulence and AMR genes, belonging to at least nine O:H serotypes, and are genetically clustered with at least five different groups of E. coli isolates reported by other states in the United States. Despite the small sample size, this study suggested that there was a large extent of genomic and serological diversity among E. coli isolates in southern United States poultry. A large-scale comprehensive study is needed to understand the overall genomic diversity and the associated virulence, and such a study will be important to develop a broadly protective E. coli vaccine.

Benchmarking Low-Frequency Variant Calling With Long-Read Data on Mitochondrial DNA

Background: Sequencing quality has improved over the last decade for long-reads, allowing for more accurate detection of somatic low-frequency variants. In this study, we used mixtures of mitochondrial samples with different haplogroups (i.e., a specific set of mitochondrial variants) to investigate the applicability of nanopore sequencing for low-frequency single nucleotide variant detection.

Methods: We investigated the impact of base-calling, alignment/mapping, quality control steps, and variant calling by comparing the results to a previously derived short-read gold standard generated on the Illumina NextSeq. For nanopore sequencing, six mixtures of four different haplotypes were prepared, allowing us to reliably check for expected variants at the predefined 5%, 2%, and 1% mixture levels. We used two different versions of Guppy for base-calling, two aligners (i.e., Minimap2 and Ngmlr), and three variant callers (i.e., Mutserve2, Freebayes, and Nanopanel2) to compare low-frequency variants. We used F₁ score measurements to assess the performance of variant calling.

Results: We observed a mean read length of 11 kb and a mean overall read quality of 15. Ngmlr showed not only higher F₁ scores but also higher allele frequencies (AF) of false-positive calls across the mixtures (mean F₁ score = 0.83; false-positive allele frequencies < 0.17) compared to Minimap2 (mean F₁ score = 0.82; false-positive AF < 0.06). Mutserve2 had the highest F₁ scores (5% level: F₁ score >0.99, 2% level: F₁ score >0.54, and 1% level: F₁ score >0.70) across all callers and mixture levels.

Conclusion: We here present the benchmarking for low-frequency variant calling with nanopore sequencing by identifying current limitations.

Estimating individual mtDNA haplotypes in mixed DNA samples by combining MinION and MiSeq

We tried to estimate individual mtDNA haplotypes in mixed DNA samples by combining MinION and MiSeq. The BAM files produced by MiSeq were viewed using Integrative Genomics Viewer (IGV) to verify mixed bases. By sorting the reads according to base type for each mixed base, partial haplotypes were determined. Then, the BAM files produced by MinKNOW were viewed using IGV. To determine haplotypes with IGV, only mixed bases determined by MiSeq were used as target bases. By sorting the reads according to base type for each target base, each contributor's haplotype was estimated. In mixed samples from two contributors, even a haplotype with a minor contribution of 5% could be distinguished from the haplotype of the major contributor. In mixed samples of three contributors (mixture ratios of 1:1:1 and 4:2:1), each haplotype could also be distinguished. Sequences of C-stretches were determined very inaccurately in the MinION analysis. Although the analysis method was simple, each haplotype was correctly detected in all mixed samples with two or three contributors in various mixture ratios by combining MinION and MiSeq. This should be useful for identifying contributors to mixed samples.

Forensic nanopore sequencing of microhaplotype markers using QitanTech's QNome

In recent years, extraordinary progress has been made in genome sequencing technologies, which has led to a decrease in cost and an increase in the diversity of sequenced genomes. Nanopore sequencing is one of the latest genome sequencing technologies. It aims to sequence longer contiguous pieces of DNA, which are essential for resolving structurally complex regions, and provides a new approach for forensic genetics to detect longer markers in real time. To date, multiple studies have been conducted to sequence forensic markers using MinION from Oxford Nanopore Technologies (ONT), and the results indicate that nanopore sequencing holds promise for forensic applications. Qitan Technology (QitanTech) recently launched its first commercial nanopore genome sequencer, QNome. It could achieve a read length of more than 150 kbp, and could generate approximately 500 Mb of data in 8 h. In this pilot study, we explored and validated this alternative nanopore sequencing device for microhaplotype (MH) profiling using a custom set of 15 MH loci. Seventy single-contributor samples were divided into 7 batches, each of which included 10 samples and control DNA 9947A and was sequenced by QNome. MH genotypes generated from QNome were compared to those from Ion Torrent sequencing (Ion S5XL system) to evaluate the accuracy and stability. Twelve samples randomly selected from the last three batches and Control DNA 9947A were also subjected to ONT MinION sequencing (with R9.4 flow cell) for parallel comparison. Based on MHtyper, a bioinformatics workflow developed for automated MH designation, all MH loci can be genotyped and reliably phased using the QNome data, with an overall accuracy of 99.83% (4 errors among 2310 genotypes). Three occurred near or in the region of homopolymer sequences, and one existed within 50 bp of the start of the sequencing reaction. In the last 15 samples (12 individual samples and 3 replicates of control DNA 9947A), two SNPs located at 4-mer homopolymers failed to obtain reliable genotypes on the MinION data. This study shows the potential of state-of-the-art nanopore sequencing methods to analyze forensic MH markers. Given the rapid pace of change, sporadic and nonrepetitive errors presented in this study are expected to be resolved by further developments of nanopore technologies and analysis tools.

Comparing two seized drug workflows for the analysis of synthetic cannabinoids, cathinones, and opioids

As the challenges faced by drug chemists persist, due to the presence of emerging drugs, laboratories continue to look for new solutions, ranging from existing methods to implementation of entirely new technology. A common barrier for making workflow changes is a lack of pre-existing data demonstrating the potential impact of these changes. In this study, we compare, qualitatively and quantitatively, an existing workflow for seized drug analysis to an experimental workflow. Four chemists were asked to analyze a total of 50 mock case samples across the two workflows. The existing workflow employed color tests for screening alongside general purpose GC-FID and GC-MS analyses for confirmation. The experimental workflow combined DART-MS screening with class-specific (targeted) GC-MS analysis for confirmation. Comparison of the workflows showed that screening by DART-MS required the same amount of time as color tests but yielded more accurate and specific information. Confirmation using the existing workflow required more than twice the amount of instrument time and data interpretation time while also presenting other analytical challenges that prevented compound confirmation in select samples. Targeted GC-MS methods simplified data interpretation, reduced consumption of reference materials, and addressed almost all limitations of general-purpose methods. While the experimental workflow requires modifications and answering of additional research questions, this study shows how rethinking analytical workflows for seized drug analysis could reduce turnaround times, backlogs, and standards consumption. It also demonstrates the potential impact of being able to investigate workflow changes prior to implementation.

An inducible transposon mutagenesis approach for the intracellular human pathogen Chlamydia trachomatis

Background: Chlamydia trachomatis is a prolific human pathogen that can cause serious long-term conditions if left untreated. Recent developments in Chlamydia genetics have opened the door to conducting targeted and random mutagenesis experiments to identify gene function. In the present study, an inducible transposon mutagenesis approach was developed for C. trachomatis using a self-replicating vector to deliver the transposon-transposase cassette - a significant step towards our ultimate aim of achieving saturation mutagenesis of the Chlamydia genome. Methods: The low transformation efficiency of C. trachomatis necessitated the design of a self-replicating vector carrying the transposon mutagenesis cassette (i.e. the Himar-1 transposon containing the beta lactamase gene as well as a hyperactive transposase gene under inducible control of the tet promoter system with the addition of a riboswitch). Chlamydia transformed with this vector (pSW2-RiboA-C9Q) were induced at 24 hours post-infection. Through dual control of transcription and translation, basal expression of transposase was tightly regulated to stabilise the plasmid prior to transposition. Results: Here we present the preliminary sequencing results of transposon mutant pools of both C. trachomatis biovars, using two plasmid-free representatives: urogenital strain   C. trachomatis SWFP- and the lymphogranuloma venereum isolate L2(25667R). DNA sequencing libraries were generated and analysed using Oxford Nanopore Technologies' MinION technology. This enabled 'proof of concept' for the methods as an initial low-throughput screen of mutant libraries; the next step is to employ high throughput sequencing to assess saturation mutagenesis. Conclusions: This significant advance provides an efficient method for assaying C. trachomatis gene function and will enable the identification of the essential gene set of C. trachomatis. In the long-term, the methods described herein will add to the growing knowledge of chlamydial infection biology leading to the discovery of novel drug or vaccine targets.

Enhanced mixture interpretation with macrohaplotypes based on long-read DNA sequencing

Deconvoluting mixture samples is one of the most challenging problems confronting DNA forensic laboratories. Efforts have been made to provide solutions regarding mixture interpretation. The probabilistic interpretation of Short Tandem Repeat (STR) profiles has increased the number of complex mixtures that can be analyzed. A portion of complex mixture profiles, particularly for mixtures with a high number of contributors, are still being deemed uninterpretable. Novel forensic markers, such as Single Nucleotide Variants (SNV) and microhaplotypes, also have been proposed to allow for better mixture interpretation. However, these markers have both a lower discrimination power compared with STRs and are not compatible with CODIS or other national DNA databanks worldwide. The short-read sequencing (SRS) technologies can facilitate mixture interpretation by identifying intra-allelic variations within STRs. Unfortunately, the short size of the amplicons containing STR markers and sequence reads limit the alleles that can be attained per STR. The latest long-read sequencing (LRS) technologies can overcome this limitation in some samples in which larger DNA fragments (including both STRs and SNVs) with definitive phasing are available. Based on the LRS technologies, this study developed a novel CODIS compatible forensic marker, called a macrohaplotype, which combines a CODIS STR and flanking variants to offer extremely high number of haplotypes and hence very high discrimination power per marker. The macrohaplotype will substantially improve mixture interpretation capabilities. Based on publicly accessible data, a panel of 20 macrohaplotypes with sizes of ~ 8 k bp and the maximum high discrimination powers were designed. The statistical evaluation demonstrates that these macrohaplotypes substantially outperform CODIS STRs for mixture interpretation, particularly for mixtures with a high number of contributors, as well as other forensic applications. Based on these results, efforts should be undertaken to build a complete workflow, both wet-lab and bioinformatics, to precisely call the variants and generate the macrohaplotypes based on the LRS technologies.

A benchmarking of human mitochondrial DNA haplogroup classifiers from whole-genome and whole-exome sequence data

The mitochondrial genome (mtDNA) is of interest for a range of fields including evolutionary, forensic, and medical genetics. Human mitogenomes can be classified into evolutionary related haplogroups that provide ancestral information and pedigree relationships. Because of this and the advent of high-throughput sequencing (HTS) technology, there is a diversity of bioinformatic tools for haplogroup classification. We present a benchmarking of the 11 most salient tools for human mtDNA classification using empirical whole-genome (WGS) and whole-exome (WES) short-read sequencing data from 36 unrelated donors. We also assessed the best performing tool in third-generation long noisy read WGS data obtained with nanopore technology for a subset of the donors. We found that, for short-read WGS, most of the tools exhibit high accuracy for haplogroup classification irrespective of the input file used for the analysis. However, for short-read WES, Haplocheck and MixEmt were the most accurate tools. Based on the performance shown for WGS and WES, and the accompanying qualitative assessment, Haplocheck stands out as the most complete tool. For third-generation HTS data, we also showed that Haplocheck was able to accurately retrieve mtDNA haplogroups for all samples assessed, although only after following assembly-based approaches (either based on a referenced-based assembly or a hybrid de novo assembly). Taken together, our results provide guidance for researchers to select the most suitable tool to conduct the mtDNA analyses from HTS data.

Nanopore Single-Molecule Sequencing for Mitochondrial DNA Methylation Analysis: Investigating Parkin-Associated Parkinsonism as a Proof of Concept

Objective: To establish a workflow for mitochondrial DNA (mtDNA) CpG methylation using Nanopore whole-genome sequencing and perform first pilot experiments on affected Parkin biallelic mutation carriers (Parkin-PD) and healthy controls. Background: Mitochondria, including mtDNA, are established key players in Parkinson's disease (PD) pathogenesis. Mutations in Parkin, essential for degradation of damaged mitochondria, cause early-onset PD. However, mtDNA methylation and its implication in PD is understudied. Herein, we establish a workflow using Nanopore sequencing to directly detect mtDNA CpG methylation and compare mtDNA methylation between Parkin-related PD and healthy individuals. Methods: To obtain mtDNA, whole-genome Nanopore sequencing was performed on blood-derived from five Parkin-PD and three control subjects. In addition, induced pluripotent stem cell (iPSC)-derived midbrain neurons from four of these patients with PD and the three control subjects were investigated. The workflow was validated, using methylated and unmethylated 897 bp synthetic DNA samples at different dilution ratios (0, 50, 100% methylation) and mtDNA without methylation. MtDNA CpG methylation frequency (MF) was detected using Nanopolish and Megalodon. Results: Across all blood-derived samples, we obtained a mean coverage of 250.3X (SD ± 80.5X) and across all neuron-derived samples 830X (SD ± 465X) of the mitochondrial genome. We detected overall low-level CpG methylation from the blood-derived DNA (mean MF ± SD = 0.029 ± 0.041) and neuron-derived DNA (mean MF ± SD = 0.019 ± 0.035). Validation of the workflow, using synthetic DNA samples showed that highly methylated DNA molecules were prone to lower Guppy Phred quality scores and thereby more likely to fail Guppy base-calling. CpG methylation in blood- and neuron-derived DNA was significantly lower in Parkin-PD compared to controls (Mann-Whitney U-test p < 0.05). Conclusion: Nanopore sequencing is a useful method to investigate mtDNA methylation architecture, including Guppy-failed reads is of importance when investigating highly methylated sites. We present a mtDNA methylation workflow and suggest methylation variability across different tissues and between Parkin-PD patients and controls as an initial model to investigate.

A New Method for Sequencing the Mitochondrial Genome by Using Long Read Technology

We describe a protocol to prepare a multiplexed mtDNA library from a blood sample for performing a long read sequencing of the mitochondrial genome. All steps are carefully described to get a high enrichment of mtDNA relative to total DNA extracted from the blood sample. The obtained mutiplexed library allows the production of long sequence mtDNA reads up to 16.5 kbp with a quality enabling variant-calling by using a portable sequencer (MinION, Oxford Nanopore Technologies).

International Wildlife Trafficking: A perspective on the challenges and potential forensic genetics solutions

International wildlife trafficking (IWT) is a thriving and pervasive illegal enterprise that adversely affects modern societies. Yet, despite being globally recognized as a threat to biodiversity, national security, economy, and biosecurity, IWT remains largely unabated and is proliferating at an alarming rate. The increase in IWT is generally attributed to a lack of prioritization to curb wildlife crime through legal and scientific infrastructure. This review: (1) lays out the damaging scope and influence of IWT; (2) discusses the potential of DNA marker systems, barcodes, and emerging molecular technologies, such as long-read portable sequencing, to facilitate rapid, in situ identification of species and individuals; and (3) encourages initiatives that promote quality and innovation. Interdisciplinary collaboration promises to be one of the most effective ways forward to surmounting the complex scientific and legal challenges posed by IWT.

Progress in forensic bone DNA analysis: Lessons learned from ancient DNA

Research on ancient and forensic DNA is related in many ways, and the two fields must deal with similar obstacles. Therefore, communication between these two communities has the potential to improve results in both research fields. Here, we present the insights gained in the ancient DNA community with regard to analyzing DNA from aged skeletal material and the potential use of the developed protocols in forensic work. We discuss the various steps, from choosing samples for DNA extraction to deciding between classical PCR amplification and massively parallel sequencing approaches. Based on the progress made in ancient DNA analyses combined with the requirements of forensic work, we suggest that there is substantial potential for incorporating ancient DNA approaches into forensic protocols, a process that has already begun to a considerable extent. However, taking full advantage of the experiences gained from ancient DNA work will require comparative studies by the forensic DNA community to tailor the methods developed for ancient samples to the specific needs of forensic studies and case work. If successful, in our view, the benefits for both communities would be considerable.

A new method for long-read sequencing of animal mitochondrial genomes: application to the identification of equine mitochondrial DNA variants

BACKGROUND

Mitochondrial DNA is remarkably polymorphic. This is why animal geneticists survey mitochondrial genomes variations for fundamental and applied purposes. We present here an approach to sequence whole mitochondrial genomes using nanopore long-read sequencing. Our method relies on the selective elimination of nuclear DNA using an exonuclease treatment and on the amplification of circular mitochondrial DNA using a multiple displacement amplification step.

RESULTS

We optimized each preparative step to obtain a 100 million-fold enrichment of horse mitochondrial DNA relative to nuclear DNA. We sequenced these amplified mitochondrial DNA using nanopore sequencing technology and obtained mitochondrial DNA reads that represented up to half of the sequencing output. The sequence reads were 2.3 kb of mean length and provided an even coverage of the mitochondrial genome. Long-reads spanning half or more of the whole mtDNA provided a coverage that varied between 118X and 488X. We evaluated SNPs identified using these long-reads by Sanger sequencing as ground truth and found a precision of 100.0%; a recall of 93.1% and a F1-score of 0.964 using the Twilight horse mtDNA reference. The choice of the mtDNA reference impacted variant calling efficiency with F1-scores varying between 0.947 and 0.964.

CONCLUSIONS

Our method to amplify mtDNA and to sequence it using the nanopore technology is usable for mitochondrial DNA variant analysis. With minor modifications, this approach could easily be applied to other large circular DNA molecules.

Comparative sequencing data analysis of Ion Torrent and MinION sequencing platforms using a clinical diagnostic haematology panel

INTRODUCTION

Currently, two second-generation sequencing platforms-Ion Torrent and Illumina are being widely used for clinical testing and reporting of human samples. However, third-generation long read platforms like single molecule (PacBio) or direct sequencing (Nanopore) are gaining widespread interest in clinical genomics.

AIMS AND OBJECTIVES

The current study attempts to analyse and compare MinION sequencing data with Ion Torrent data, using a haematology ampliseq panel, to shed light on its current standing in reporting of clinical diagnostic samples.

METHODOLOGY

A custom targeted Next-generation sequencing ampliseq panel comprising of 33 genes related to detection of inherited bone marrow failure syndrome cases was used to sequence five samples on both Ion Torrent and Oxford MinION platforms. The resulting data were analysed for output, quality and variant metrics across both platforms independently.

RESULTS

Overall, MinION produced longer reads (range 108-7227 bp) than Ion Torrent (25-580 bp). Moreover, it generated more reads than Ion Torrent in high %GC content (P < .005) or repeat regions. But Ion Torrent had much lower error rate of 1.59% compared with MinION's 20.31%. Despite high error, MinION platform was able to identify and report the pathogenic variant in all samples.

DISCUSSION AND CONCLUSION

The extremely long read lengths of MinION sequencers and better coverage in difficult to sequence regions give it an edge in generating contig free whole-genome sequences. However, the pore technology and chemistry of MinION needs further tuning to reduce error rate before it can be incorporated for clinical testing and reporting of human samples.

Shotgun metagenome data of a defined mock community using Oxford Nanopore, PacBio and Illumina technologies

Metagenomic sequence data from defined mock communities is crucial for the assessment of sequencing platform performance and downstream analyses, including assembly, binning and taxonomic assignment. We report a comparison of shotgun metagenome sequencing and assembly metrics of a defined microbial mock community using the Oxford Nanopore Technologies (ONT) MinION, PacBio and Illumina sequencing platforms. Our synthetic microbial community BMock12 consists of 12 bacterial strains with genome sizes spanning 3.2–7.2 Mbp, 40–73% GC content, and 1.5–7.3% repeats. Size selection of both PacBio and ONT sequencing libraries prior to sequencing was essential to yield comparable relative abundances of organisms among all sequencing technologies. While the Illumina-based metagenome assembly yielded good coverage with few misassemblies, contiguity was greatly improved by both, Illumina + ONT and Illumina + PacBio hybrid assemblies but increased misassemblies, most notably in genomes with high sequence similarity to each other. Our resulting datasets allow evaluation and benchmarking of bioinformatics software on Illumina, PacBio and ONT platforms in parallel.

Measurement(s)	metagenomic data • sequence_assembly
Technology Type(s)	ONT MinION • Illumina sequencing • PacBio RS II
Factor Type(s)	sequencing platform
Sample Characteristic - Organism	Bacteria
Sample Characteristic - Environment	mock community

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.10260740

Next-generation sequencing for the clinical management of hepatitis C virus infections: does one test fits all purposes?

While the prospect of viral cure is higher than ever for individuals infected with the hepatitis C virus (HCV) due to ground-breaking progress in antiviral treatment, success rates are still negatively influenced by HCV's high genetic variability. This genetic diversity is represented in the circulation of various genotypes and subtypes, mixed infections, recombinant forms and the presence of numerous drug resistant variants among infected individuals. Common misclassifications by commercial genotyping assays in combination with the limitations of currently used targeted population sequencing approaches have encouraged researchers to exploit alternative methods for the clinical management of HCV infections. Next-generation sequencing (NGS), a revolutionary and powerful tool with a variety of applications in clinical virology, can characterize viral diversity and depict viral dynamics in an ultra-wide and ultra-deep manner. The level of detail it provides makes it the method of choice for the diagnosis and clinical assessment of HCV infections. The sequence library provided by NGS is of a higher magnitude and sensitivity than data generated by conventional methods. Therefore, these technologies are helpful to guide clinical practice and at the same time highly valuable for epidemiological studies. The decreasing costs of NGS to determine genotypes, mixed infections, recombinant strains and drug resistant variants will soon make it feasible to employ NGS in clinical laboratories, to assist in the daily care of patients with HCV.

Using the Plasmodium mitochondrial genome for classifying mixed-species infections and inferring the geographical origin of P. falciparum parasites imported to the U.S

The ability to identify mixed-species infections and track the origin of Plasmodium parasites can further enhance the development of treatment and prevention recommendations as well as outbreak investigations. Here, we explore the utility of using the full Plasmodium mitochondrial genome to classify Plasmodium species, detect mixed infections, and infer the geographical origin of imported P. falciparum parasites to the United States (U.S.). Using the recently developed standardized, high-throughput Malaria Resistance Surveillance (MaRS) protocol, the full Plasmodium mitochondrial genomes of 265 malaria cases imported to the U.S. from 2014-2017 were sequenced and analyzed. P. falciparum infections were found in 94.7% (251/265) of samples. Five percent (14/265) of samples were identified as mixed- Plasmodium species or non-P. falciparum, including P. vivax, P. malariae, P. ovale curtisi, and P. ovale wallikeri. P. falciparum mitochondrial haplotypes analysis revealed greater than eighteen percent of samples to have at least two P. falciparum mitochondrial genome haplotypes, indicating either heteroplasmy or multi-clonal infections. Maximum-likelihood phylogenies of 912 P. falciparum mitochondrial genomes with known country origin were used to infer the geographical origin of thirteen samples from persons with unknown travel histories as: Africa (country unspecified) (n = 10), Ghana (n = 1), Southeast Asia (n = 1), and the Philippines (n = 1). We demonstrate the utility and current limitations of using the Plasmodium mitochondrial genome to classify samples with mixed-infections and infer the geographical origin of imported P. falciparum malaria cases to the U.S. with unknown travel history.

Nanopore sequencing: An enrichment-free alternative to mitochondrial DNA sequencing

Mitochondrial DNA sequence data are often utilized in disease studies, conservation genetics and forensic identification. The current approaches for sequencing the full mtGenome typically require several rounds of PCR enrichment during Sanger or MPS protocols followed by fairly tedious assembly and analysis. Here we describe an efficient approach to sequencing directly from genomic DNA samples without prior enrichment or extensive library preparation steps. A comparison is made between libraries sequenced directly from native DNA and the same samples sequenced from libraries generated with nine overlapping mtDNA amplicons on the Oxford Nanopore MinION™ device. The native and amplicon library preparation methods and alternative base calling strategies were assessed to establish error rates and identify trends of discordance between the two library preparation approaches. For the complete mtGenome, 16 569 nucleotides, an overall error rate of approximately 1.00% was observed. As expected with mtDNA, the majority of error was detected in homopolymeric regions. The use of a modified basecaller that corrects for ambiguous signal in homopolymeric stretches reduced the error rate for both library preparation methods to approximately 0.30%. Our study indicates that direct mtDNA sequencing from native DNA on the MinION™ device provides comparable results to those obtained from common mtDNA sequencing methods and is a reliable alternative to approaches using PCR-enriched libraries.

Recent progress, methods and perspectives in forensic epigenetics

Forensic epigenetics, i.e., investigating epigenetics variation to resolve forensically relevant questions unanswerable with standard forensic DNA profiling has been gaining substantial ground over the last few years. Differential DNA methylation among tissues and individuals has been proposed as useful resource for three forensic applications i) determining the tissue type of a human biological trace, ii) estimating the age of an unknown trace donor, and iii) differentiating between monozygotic twins. Thus far, forensic epigenetic investigations have used a wide range of methods for CpG marker discovery, prediction modelling and targeted DNA methylation analysis, all coming with advantages and disadvantages when it comes to forensic trace analysis. In this review, we summarize the most recent literature on these three main topics of current forensic epigenetic investigations and discuss limitations and practical considerations in experimental design and data interpretation, such as technical and biological biases. Moreover, we provide future perspectives with regard to new research questions, new epigenetic markers and recent technological advances that - as we envision - will move the field towards forensic epigenomics in the near future.

Massively parallel sequencing-enabled mixture analysis of mitochondrial DNA samples

The mitochondrial genome has a number of characteristics that provide useful information to forensic investigations. Massively parallel sequencing (MPS) technologies offer improvements to the quantitative analysis of the mitochondrial genome, specifically the interpretation of mixed mitochondrial samples. Two-person mixtures with nuclear DNA ratios of 1:1, 5:1, 10:1, and 20:1 of individuals from different and similar phylogenetic backgrounds and three-person mixtures with nuclear DNA ratios of 1:1:1 and 5:1:1 were prepared using the Precision ID mtDNA Whole Genome Panel and Ion Chef, and sequenced on the Ion PGM or Ion S5 sequencer (Thermo Fisher Scientific, Waltham, MA, USA). These data were used to evaluate whether and to what degree MPS mixtures could be deconvolved. Analysis was effective in identifying the major contributor in each instance, while SNPs from the minor contributor's haplotype only were identified in the 1:1, 5:1, and 10:1 two-person mixtures. While the major contributor was identified from the 5:1:1 mixture, analysis of the three-person mixtures was more complex, and the mixed haplotypes could not be completely parsed. These results indicate that mixed mitochondrial DNA samples may be interpreted with the use of MPS technologies.

Analysing Microbial Community Composition through Amplicon Sequencing: From Sampling to Hypothesis Testing

Microbial ecology as a scientific field is fundamentally driven by technological advance. The past decade's revolution in DNA sequencing cost and throughput has made it possible for most research groups to map microbial community composition in environments of interest. However, the computational and statistical methodology required to analyse this kind of data is often not part of the biologist training. In this review, we give a historical perspective on the use of sequencing data in microbial ecology and restate the current need for this method; but also highlight the major caveats with standard practices for handling these data, from sample collection and library preparation to statistical analysis. Further, we outline the main new analytical tools that have been developed in the past few years to bypass these caveats, as well as highlight the major requirements of common statistical practices and the extent to which they are applicable to microbial data. Besides delving into the meaning of select alpha- and beta-diversity measures, we give special consideration to techniques for finding the main drivers of community dissimilarity and for interaction network construction. While every project design has specific needs, this review should serve as a starting point for considering what options are available.