Development and assessment of an optimized next-generation DNA sequencing approach for the mtgenome using the Illumina MiSeq

Assessing the feasibility of free DNA for disaster victim identification and forensic applications

In tropical disaster victim identification (DVI) scenarios, challenging environmental conditions lead to accelerated DNA degradation in remains. To further enhance the utilization of leached DNA from tissue in the preservative solution (termed "free DNA") as an alternative source, we incorporated new results by assessing its integrity in postmortem and decomposing cadavers preserved in DNA/RNA Shield™ and modified TENT, with silica-based purification (QIAquick®) for faster processing. The psoas muscle tissues of one decomposed and ten cadavers were preserved in each solution at 25 °C and 35 °C for 3 months. Free DNA efficiency was compared with individual reference samples for reliable results in quantity, quality, and STR profiles. The findings revealed that DNA/RNA Shield™ effectively preserves free DNA integrity for extended storage, while modified TENT is more suitable for short-term storage due to higher degradation levels. Moreover, the use of free DNA samples with massive parallel sequencing displays potential for forensic DNA analysis. Successful amplification of the mtDNA control region enables variant calling and heteroplasmy analysis while also serving as quality control using ACTB and enabling differentiation within the 16S rRNA region for microbiome analysis. The simplicity of handling free DNA for PCR-based forensic analysis adds to its potential for various applications, including DVI and field-based analysis of biological evidence.

Discrimination of monozygotic twins using mtDNA heteroplasmy through probe capture enrichment and massively parallel sequencing

Differentiating between monozygotic (MZ) twins remains difficult because they have the same genetic makeup. Applying the traditional STR genotyping approach cannot differentiate one from the other. Heteroplasmy refers to the presence of two or more different mtDNA copies within a single cell and this phenomenon is common in humans. The levels of heteroplasmy cannot change dramatically during transmission in the female germ line but increase or decrease during germ-line transmission and in somatic tissues during life. As massively parallel sequencing (MPS) technology has advanced, it has shown the extraordinary quantity of mtDNA heteroplasmy in humans. In this study, a probe hybridization technique was used to obtain mtDNA and then MPS was performed with an average sequencing depth of above 4000. The results showed us that all ten pairs of MZ twins were clearly differentiated with the minor heteroplasmy threshold at 1.0%, 0.5%, and 0.1%, respectively. Finally, we used a probe that targeted mtDNA to boost sequencing depth without interfering with nuclear DNA and this technique can be used in forensic genetics to differentiate the MZ twins.

Mitochondrial DNA control region typing from highly degraded skeletal remains by single-multiplex next-generation sequencing

Poor nuclear DNA preservation from highly degraded skeletal remains is the most limiting factor for the genetic identification of individuals. Mitochondrial DNA (mtDNA) typing, and especially of the control region (CR), using next-generation sequencing (NGS), enables retrieval of valuable genetic information in forensic contexts where highly degraded human skeletal remains are the only source of genetic material. Currently, NGS commercial kits can type all mtDNA-CR in fewer steps than the conventional Sanger technique. The PowerSeq CRM Nested System kit (Promega Corporation) employs a nested multiplex-polymerase chain reaction (PCR) strategy to amplify and index all mtDNA-CR in a single reaction. Our study analyzes the success of mtDNA-CR typing of highly degraded human skeletons using the PowerSeq CRM Nested System kit. We used samples from 41 individuals from different time periods to test three protocols (M1, M2, and M3) based on modifications of PCR conditions. To analyze the detected variants, two bioinformatic procedures were compared: an in-house pipeline and the GeneMarker HTS software. The results showed that many samples were not analyzed when the standard protocol (M1) was used. In contrast, the M3 protocol, which includes 35 PCR cycles and longer denaturation and extension steps, successfully recovered the mtDNA-CR from highly degraded skeletal samples. Mixed base profiles and the percentage of damaged reads were both indicators of possible contamination and can provide better results if used together. Furthermore, our freely available in-house pipeline can provide variants concordant with the forensic software.

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.

Exploring rare differences in mitochondrial genome between MZ twins using Ion Torrent semiconductor sequencing

Monozygotic (MZ) twins are considered to be genetically identical in that they have the same genomic DNA sequences in theory, and thus cannot be differentiated using forensic standard STR-based DNA profiling. However, a recent study employed deep sequencing to explore extremely rare mutations in the nuclear genome and reported that the mutation analysis could be applied to differentiate between MZ twins. Compared with the nuclear genome, the mitochondrial DNA (mtDNA) exhibits higher mutation rates due to fewer DNA repair mechanisms in the mitochondrial genome (mtGenome) and the lack of proofreading capability of the mtDNA polymerase. In a previous study, we used Illumina ultra-deep sequencing to describe point heteroplasmy (PHP) and nucleotide variant of the mtGenomes in venous blood samples of MZ twins. In the present study, we characterized minor differences of the mtGenomes in three tissue samples from seven sets of MZ twins using Ion Torrent semiconductor sequencing (Thermo Fisher Ion S5 XL system) and commercialized mtGenome sequencing kit (Precision ID mtDNA Whole Genome Panel). PHP was observed in blood samples from one set of MZ twins and in saliva samples from two sets of twins, but it presented in hair shaft samples from all seven sets of MZ twins. Overall, the coding region of the mtGenome exhibits more PHPs than the control region. The results of this study have further attested the competence of mtGenome sequencing in differentiating between MZ twins, and that among the three kinds of samples tested, hair shaft is more likely to accumulate minor differences in the mtGenomes of MZ twins.

Massively parallel sequencing of Cannabis sativa chloroplast hotspots for forensic typing

BACKGROUND

Marijuana (Cannabis sativa) is the most commonly used illicit drug in the USA, and the use of DNA barcodes could assist drug trafficking investigations by indicating the biogeographical origin and crop type of a sample and providing a means for linking cases. Additionally, the legality of marijuana in the USA remains complicated with some states fully legalizing marijuana for recreational use while federally marijuana remains completely illegal. Massively parallel sequencing (MPS) offers distinct advantages over capillary electrophoresis (CE), including more comprehensive coverage of target loci, analysis of hundreds of markers simultaneously, and high throughput capabilities.

METHODS

This study reports on the development of a MiSeq FGx® assay targeting seven "hotspot" regions in the Cannabis sativa chloroplast genome that are highly polymorphic and informative in attempts to determine biogeographical origin and distinguishing between marijuana and hemp. Sequencing results were compared to previous studies that used CE-based genotyping methods.

RESULTS

A total of 49 polymorphisms were observed, 16 of which have not been previously reported. Additionally, sequence data revealed isoalleles at one locus, which were able to differentiate two samples that had the same haplotype using CE-based methods. This study reports preliminary results from sequencing 14 hemp and marijuana samples from different countries using the developed MPS assay.

CONCLUSION

Future studies should genotype a more comprehensive sample set from around the world to build a haplotype database, which could be used to provide investigative leads for law enforcement agencies investigating marijuana trafficking.

Exploring statistical weight estimates for mitochondrial DNA matches involving heteroplasmy

Massively parallel sequencing (MPS) of mitochondrial (mt) DNA allows forensic laboratories to report heteroplasmy on a routine basis. Statistical approaches will be needed to determine the relative frequency of observing an mtDNA haplotype when including the presence of a heteroplasmic site. Here, we examined 1301 control region (CR) sequences, collected from individuals in four major population groups (European, African, Asian, and Latino), and covering 24 geographically distributed haplogroups, to assess the rates of point heteroplasmy (PHP) on an individual and nucleotide position (np) basis. With a minor allele frequency (MAF) threshold of 2%, the data was similar across population groups, with an overall PHP rate of 37.7%, and the majority of heteroplasmic individuals (77.3%) having only one site of heteroplasmy. The majority (75.2%) of identified PHPs had an MAF of 2-10%, and were observed at 12.6% of the nps across the CR. Both the broad and phylogenetic testing suggested that in many cases the low number of observations of heteroplasmy at any one np results in a lack of statistical association. The posterior frequency estimates, which skew conservative to a degree depending on the sample size in a given haplogroup, had a mean of 0.152 (SD 0.134) and ranged from 0.031 to 0.83. As expected, posterior frequency estimates decreased in accordance with 1/n as the sample size (n) increased. This provides a proposed conservative statistical framework for assessing haplotype/heteroplasmy matches when applying an MPS technique in forensic cases and will allow for continual refinement as more data is generated, both within the CR and across the mitochondrial genome.

The frequency of the known mitochondrial variants associated with drug-induced toxicity in a Korean population

BACKGROUND

Few studies have annotated the whole mitochondrial DNA (mtDNA) genome associated with drug responses in Asian populations. This study aimed to characterize mtDNA genetic profiles, especially the distribution and frequency of well-known genetic biomarkers associated with diseases and drug-induced toxicity in a Korean population.

METHOD

Whole mitochondrial genome was sequenced for 118 Korean subjects by using a next-generation sequencing approach. The bioinformatic pipeline was constructed for variant calling, haplogroup classification and annotation of mitochondrial mutation.

RESULTS

A total of 681 variants was identified among all subjects. The MT-TRNP gene and displacement loop showed the highest numbers of variants (113 and 74 variants, respectively). The m.16189T > C allele, which is known to reduce the mtDNA copy number in human cells was detected in 25.4% of subjects. The variants (m.2706A > G, m.3010A > G, and m.1095T > C), which are associated with drug-induced toxicity, were observed with the frequency of 99.15%, 30.51%, and 0.08%, respectively. The m.2150T > A, a genotype associated with highly disruptive effects on mitochondrial ribosomes, was identified in five subjects. The D and M groups were the most dominant groups with the frequency of 34.74% and 16.1%, respectively.

CONCLUSIONS

Our finding was consistent with Korean Genome Project and well reflected the unique profile of mitochondrial haplogroup distribution. It was the first study to annotate the whole mitochondrial genome with drug-induced toxicity to predict the ADRs event in clinical implementation for Korean subjects. This approach could be extended for further study for validation of the potential ethnic-specific mitochondrial genetic biomarkers in the Korean population.

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.

Mitochondrial genome sequencing with short overlapping amplicons on MiSeq FGx system

With the development and maturation of massively parallel sequencing (MPS) technology, the mitochondrial genome (mitogenome) sequencing is increasingly applied in the forensic field. In this study, we employed the strategy of short overlapping amplicons for the whole mitogenome, library preparation with tagmentation using the Nextera® XT DNA Library Preparation Kit, sequencing on the MiSeq FGxTM Forensic Genomics System and analyzing data using the mitochondrial(mtDNA) MSR Plug-in and the mtDNA Variant Analyzer. A total of 27 libraries and 56 libraries were sequenced in a run using MiSeq Reagent Kit v2 and v3, respectively. Results showed more than 1800 × of averaged depth of coverage (DoC) at each position. Concordant haplotypes of 9947 A and 2800 M were obtained at 32 variants. Cross-reactivity was observed with 1 ng primate DNA and 10 ng non-primate DNA but could be easily distinguished. Full and accurate variants were obtained from at least 50 pg input DNA and from minor contributors between 19:1 and 1:19 mixed ratios with known reference profiles. More than 86% variants were detected from ≥200-bp degraded samples but its haplotype was assigned to more ancestral haplogroup. Further, a total of 3 962 variants were observed at 613 nucleotide positions from 103 Xibe mitogenomes with 25:1 ratio of transitions to transversions. Two new transversions (C13735A and A14755C) and two tri-alleles at nps 9824 and 16092 were identified. There were 103 unique mitogenome haplotypes from 103 Chinese Xibe that were assigned to 79 haplogroups. Haplogroup D was the preponderant top-level haplogroup in Xibe followed by F, B, M, A, N, G, C, Z, Y, HV and J. Random match probability (RMP) and haplotype diversity (HD) of the whole mitogenome was calculated as 0.0097 and 1.0000, respectively. Compared with HVS-I only, RMP decreased 33.56%, while the number of haplotypes and HD increased 15.73% and 0.49%, respectively. Principal component analysis (PCA) showed that Xibe was clustered to East and Southeast Asian. As a whole, this MPS strategy is suitable for the whole mitogenome sequencing especially for degraded samples and can facilitate generating mitogenome data to support the routine application in forensic sciences. EMP00726 is the first whole mitogenome dataset from Xibe contributed to the EMPOP. Supplemental data for this article are available online at.

Genetic Evidence of Authenticity of a Hair Shaft Relic from the Portrait of Tsesarevich Alexei, Son of the Last Russian Emperor

To determine the value of an art piece, authenticity of the artwork must be verified. We demonstrate here a genetic approach to determine origin of a historic relic in the museum piece. We tested two hair shafts of unknown origin framed into a watercolor portrait of Tsesarevich Alexei Romanov, son of the last Russian Tzar Nicholas II, which is a unique item kept in the State Historical Museum. Genetic identification of the hair shafts was performed by analysis of mitochondrial DNA (mtDNA) markers using both massive parallel genomic sequencing and multiplex targeted PCR, followed by Sanger sequencing. In previous works, we reconstructed the complete mtDNA sequence inherited to Alexei Romanov through the Queen Victoria lineage [Rogaev et al. (2009) Proc. Natl. Acad. Sci. USA, 106, 5258-5263]. DNA extracts were obtained from the two thin hair shafts and used for comparative genetic analysis. Despite the very low quantity and quality of the DNA templates retrieved from the historical single hair shaft specimen, informative mtDNA sequences were determined. The mtDNA haplotype in the hair shafts corresponds to the mtDNA haplotype of Tsarevich Alexei, his sisters, and his mother, Empress Alexandra Feodorovna. This haplotype remains unique in the currently available mtDNA databases. Our results reveal that the hair relic from the portrait is associated with the family of the last Russian Emperor Nicholas II. The study is an example of first application of the genetic methodology for verification of the value of museum artwork items.

iPSC culture expansion selects against putatively actionable mutations in the mitochondrial genome

Therapeutic application of induced pluripotent stem cell (iPSC) derivatives requires comprehensive assessment of the integrity of their nuclear and mitochondrial DNA (mtDNA) to exclude oncogenic potential and functional deficits. It is unknown, to which extent mtDNA variants originate from their parental cells or from de novo mutagenesis, and whether dynamics in heteroplasmy levels are caused by inter- and intracellular selection or genetic drift. Sequencing of mtDNA of 26 iPSC clones did not reveal evidence for de novo mutagenesis, or for any selection processes during reprogramming or differentiation. Culture expansion, however, selected against putatively actionable mtDNA mutations. Altogether, our findings point toward a scenario in which intracellular selection of mtDNA variants during culture expansion shapes the mutational landscape of the mitochondrial genome. Our results suggest that intercellular selection and genetic drift exert minor impact and that the bottleneck effect in context of the mtDNA genetic pool might have been overestimated.

•

Expansion culture selects against putatively actionable mtDNA mutations in iPSCs

•

Intracellular selection on mtDNA molecules shapes the mutational landscape

•

Random genetic drift and intercellular selection exert minor impact

•

Selection acts during culture expansion but not during reprogramming or differentiation

Benchmarking laboratory processes to characterise low-biomass respiratory microbiota

The low biomass of respiratory samples makes it difficult to accurately characterise the microbial community composition. PCR conditions and contaminating microbial DNA can alter the biological profile. The objective of this study was to benchmark the currently available laboratory protocols to accurately analyse the microbial community of low biomass samples. To study the effect of PCR conditions on the microbial community profile, we amplified the 16S rRNA gene of respiratory samples using various bacterial loads and different number of PCR cycles. Libraries were purified by gel electrophoresis or AMPure XP and sequenced by V2 or V3 MiSeq reagent kits by Illumina sequencing. The positive control was diluted in different solvents. PCR conditions had no significant influence on the microbial community profile of low biomass samples. Purification methods and MiSeq reagent kits provided nearly similar microbiota profiles (paired Bray–Curtis dissimilarity median: 0.03 and 0.05, respectively). While profiles of positive controls were significantly influenced by the type of dilution solvent, the theoretical profile of the Zymo mock was most accurately analysed when the Zymo mock was diluted in elution buffer (difference compared to the theoretical Zymo mock: 21.6% for elution buffer, 29.2% for Milli-Q, and 79.6% for DNA/RNA shield). Microbiota profiles of DNA blanks formed a distinct cluster compared to low biomass samples, demonstrating that low biomass samples can accurately be distinguished from DNA blanks. In summary, to accurately characterise the microbial community composition we recommend 1. amplification of the obtained microbial DNA with 30 PCR cycles, 2. purifying amplicon pools by two consecutive AMPure XP steps and 3. sequence the pooled amplicons by V3 MiSeq reagent kit. The benchmarked standardized laboratory workflow presented here ensures comparability of results within and between low biomass microbiome studies.

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.

Interpreting NUMTs in forensic genetics: Seeing the forest for the trees

Nuclear mitochondrial DNA (mtDNA) segments (NUMTs) were discovered shortly after sequencing the first human mitochondrial genome. They have earlier been considered to represent archaic elements of ancient insertion events, but modern sequencing technologies and growing databases of mtDNA and NUMT sequences confirm that they are abundant and some of them phylogenetically young. Here, we build upon mtDNA/NUMT review articles published in the mid 2010 s and focus on the distinction of NUMTs and other artefacts that can be observed in aligned sequence reads, such as mixtures (contamination), point heteroplasmy, sequencing error and cytosine deamination. We show practical examples of the effect of the mtDNA enrichment method on the representation of NUMTs in the mapped sequence data and discuss methods to bioinformatically filter NUMTs from mtDNA reads.

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.

Analyzing Low-Level mtDNA Heteroplasmy-Pitfalls and Challenges from Bench to Benchmarking

Massive parallel sequencing technologies are promising a highly sensitive detection of low-level mutations, especially in mitochondrial DNA (mtDNA) studies. However, processes from DNA extraction and library construction to bioinformatic analysis include several varying tasks. Further, there is no validated recommendation for the comprehensive procedure. In this study, we examined potential pitfalls on the sequencing results based on two-person mtDNA mixtures. Therefore, we compared three DNA polymerases, six different variant callers in five mixtures between 50% and 0.5% variant allele frequencies generated with two different amplification protocols. In total, 48 samples were sequenced on Illumina MiSeq. Low-level variant calling at the 1% variant level and below was performed by comparing trimming and PCR duplicate removal as well as six different variant callers. The results indicate that sensitivity, specificity, and precision highly depend on the investigated polymerase but also vary based on the analysis tools. Our data highlight the advantage of prior standardization and validation of the individual laboratory setup with a DNA mixture model. Finally, we provide an artificial heteroplasmy benchmark dataset that can help improve somatic variant callers or pipelines, which may be of great interest for research related to cancer and aging.

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.

Searching the undetected mtDNA variants in forensic MPS data

The efficiency of MPS in forensic mtDNA analysis has been thoroughly proven, although a reliable and well established data evaluation still remains a critical point. Numerous bioinformatics tools have been developed, but most of them require specific operating systems and high costs, while free open-source programs with user-friendly interfaces are few. In this study, 43 full mtGenomes were sequenced using the Ion Personal Genome Machine™ (PGM™) System and analyzed utilizing the plug-in Variant Caller (TVC) of the Ion Torrent Software Suite and the mtDNA-Server (mDS), a free web-based mitochondrial analysis tool for MPS data. The outcomes of these two different analysis tools were compared to variants noted after manual inspection of the aligned reads performed using Integrative Genomics Viewer (IGV). The comparison highlighted the presence of thirty-nine discordant variant calls, which were resolved by Sanger sequencing that confirmed the presence of all variants, except for 7 deletions. The combined adoption of IGV and Sanger type sequencing confirmatory steps, in addition of TVC and mDS analysis, resulted in a more accurate variants assignment with the detection of 32 additional true polymorphisms, which were noted in the final dataset. Regarding the heteroplasmy issue, out of a total of thirty heteroplasmic variants, twenty-eight were detected by the TVC, while the mDS detected twenty-two. Overall, none of the used bioinformatics tools were the perfect choice and a secondary analysis with an expert's opinion in complete mtGenome MPS data evaluation is still required in forensic genetic analysis.

Damage patterns observed in mtDNA control region MPS data for a range of template concentrations and when using different amplification approaches

Massively parallel sequencing (MPS) of mitochondrial (mt) DNA allows practitioners the ability to fully resolve heteroplasmic sites. In forensic DNA analysis, identifying heteroplasmy (a naturally occurring mixture of two mtDNA profiles) can provide additional mtDNA profile information which can lead to an increase in the discrimination potential of an mtDNA match between an evidentiary sample and reference source. Forensic samples such as hair and skeletal remains, especially older, more compromised samples, can often exhibit DNA damage. Because both damage and heteroplasmy can manifest as a mixture of two nucleotides, it is important to differentiate between the two conditions when interpreting mtDNA MPS data. In this study, DNA damage was applied under controlled conditions to samples containing a range of template concentrations, including some with identified heteroplasmy. Damage was applied via storage in water at room temperature on samples diluted before or after storage to mimic low template scenarios. Damage was assessed with respect to the following areas: mtDNA quantification and degradation ratios, MPS read depth, MPS profile results, overall damage rates, and the interpretation of heteroplasmy. Datasets were generated to assess and compare two different amplification and library preparation strategies: the Promega PowerSeq™ CRM Nested System kit and a 1.16 kb target amplicon of the entire mtDNA control region followed by a Nextera® XT library preparation. The results of this study provide an evaluation of the Promega 10-plex MPS procedure as an improved process to mitigate the impact of mtDNA damage on low template samples. Some of the negative effects of damage observed in this study were a decrease in mtDNA yield by 20-30% and lower quality MPS sequencing results. These effects were observed more frequently when samples were diluted prior to inducing damage, illustrating that low template samples are more susceptible to damage. The findings of this study will assist forensic laboratories in differentiating between damage and heteroplasmy, which is essential when developing robust mtDNA MPS interpretation guidelines such as setting appropriate reporting thresholds.

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.

Massive parallel sequencing in forensics: advantages, issues, technicalities, and prospects

In the last decade, next-generation sequencing (NGS) technology, alternatively massive parallel sequencing (MPS), was applied to all fields of biological research. Its introduction to the field of forensics was slower, mainly due to lack of accredited sequencers, kits, and relatively higher sequencing error rates as compared with standardized Sanger sequencing. Currently, a majority of the problematic issues have been solved, which is proven by the body of reports in the literature. Here, we discuss the utility of NGS sequencing in forensics, emphasizing the advantages, issues, the technical aspects of the experiments, commercial solutions, and the potentially interesting applications of MPS.

Characterization of background noise in MiSeq MPS data when sequencing human mitochondrial DNA from various sample sources and library preparation methods

Improved resolution of massively parallel sequencing (MPS) allows for the characterization of mitochondrial (mt) DNA heteroplasmy to levels previously unattainable with traditional sequencing approaches. An essential criterion for the reporting of heteroplasmy is the ability of the MPS method to distinguish minor sequence variants (MSVs) from system noise, or error. Therefore, an assessment of the background noise in the MPS method is desirable to identify the point at which reliable data can be reported. Substitution and sequence specific error (SSE) was evaluated for a variety of sample types and two library preparations. Substitution error rates ranged from 0.18 to 0.49 per 100 nucleotides with C positions generally having the highest rate of misincorporation. Comparison of error rates across sample types indicated a significant increase for samples with damaged DNA. The positions of error were varied across datasets (pairwise concordance 0-68%), but had greater consistency within the damaged samples (80-96%). The most commonly observed motif preceding error in forward reads was CCG, while GGT was most common in reverse reads, both consistent with previous findings. The findings illustrate that for datasets containing samples with damaged DNA, reporting thresholds for heteroplasmy may have to be modified and individual sites with error levels exceeding thresholds should be scrutinized. Collectively, the shifting error profiles observed across the various sample types and library preparation methods demonstrates the need for an assessment of error under these varying circumstances. Characterization of the applicable background noise will help to ensure that thresholds are reliably set for detection of true MSVs.

Characterization of sequence variation at 30 autosomal STRs in Chinese Han and Tibetan populations

Massively parallel sequencing (MPS) technologies have the ability to reveal sequence variations within STR alleles as well as their nominal allele lengths, which have traditionally been detected by CE instruments. Recently, Thermo Fisher Scientific has updated the MPS-STR panel, named the Precision ID GlobalFiler next-generation sequencing (NGS) STR Panel version 2, with primers redesigned to add two pentanucleotide tandem repeat loci and profile interpretation supported by the Converge software. Using the Ion Chef System, the Ion S5XL System, and the Converge software, genetic variations were characterized within STR repeat and flanking regions of 30 autosomal STR markers in 115 unrelated individuals from two Chinese population groups (58 Tibetans and 57 Hans). Nineteen STRs demonstrated a relative increase in diversity with the variant sequence alleles compared with those of traditional nominal length alleles. In total, 390 alleles were identified by their sequences compared with 258 alleles that were identified by length. Of these 92 sequence variants found within the STR repeat regions, 40 variants were located in STR flanking regions. Additionally, the agreement of the results with CE data was evaluated, as was the ability of this new MPS panel to analyze case-type (11 samples) and artificially degraded samples (seven samples in triplicate). The results generated from this study illustrate that extensive sequence variation exists in commonly used STR markers in the selected population samples and indicate that this NGS STR panel has the potential to be used as an effective tool for human forensics.

Comparison of whole mitochondrial genome variants between hair shafts and reference samples using massively parallel sequencing

Hair shafts are one of the most common types of evidence at crime scenes, and mitochondrial DNA (mtDNA) has been analyzed as a valuable genetic marker for hair shafts in forensic casework. However, the mtDNA analysis strategy may vary according to the quantity and quality of DNA extracted from a forensic sample and the available massively parallel sequencing (MPS) platform in laboratories. Forensic practitioners often have to interpret mtDNA sequences exhibiting point heteroplasmy (PHP) that are analyzed using different analytical methods. In the present study, the whole mitochondrial genome (mtGenome) variants of hair shaft samples obtained from 20 donors, which were sampled in duplicate and stored at room temperature for > 1 year, were analyzed using the Precision ID mtDNA Whole Genome Panel and Ion S5 system. The whole mtGenome variants of 20 blood and 20 buccal swab samples (reference samples) from the hair shaft donors were analyzed using the Nextera XT DNA Library Prep Kit and MiSeq System. A total of 20 unique mtGenome haplotypes were observed, and 56 PHP variants were identified across the 4 sets of tissue. When the major nucleotide of PHP was considered, 16 of 20 haplotypes of the hair shaft samples matched those of the corresponding blood and buccal swab samples. In four donors, the major nucleotide of PHP was inverted at one nucleotide position between the hair shaft and reference samples. However, the data obtained on MPS, showing high PHP resolution, provided substantial information to avoid false exclusion when comparing two haplotypes containing PHP with inverted major nucleotides. In conclusion, the present study demonstrates the utility of MPS in forensic casework in the comparative analysis of mtGenome variants containing PHP.

Estimation of the number of contributors to mixed samples of DNA by mitochondrial DNA analyses using massively parallel sequencing

We evaluated whether the number of contributors to mixed DNA samples can be estimated by analyzing the D-loop of mitochondrial DNA using massively parallel sequencing. The A- (positions 16,209-16,400) and B- (positions 30-284) amplicons in hypervariable regions 1 and 2, respectively, were sequenced using MiSeq with 2 × 251 cycles. Sequence extraction and trimming were performed using CLC Genomics Workbench 11 and the number of observed haplotypes was counted for each amplicon type using Microsoft Excel. The haplotype ratios were calculated by dividing the number of counted reads of the corresponding haplotype by the total number of sequence reads. Haplotypes that were over the threshold (5%) were defined as positive haplotypes. The number of larger positive haplotypes in either of the two amplicon types was defined as the number of contributors. Samples were collected from seven individuals. Seventeen mixed samples were prepared by mixing DNA from two to five contributors at various ratios. The number of contributors was correctly estimated from almost all of the mixed samples containing equal amounts of DNA from two to five people. In mixed samples of two or three people, the minor components were detected down to a ratio of 20:1 or 8:2:1. However, heteroplasmy, base deletions, and sharing of the same haplotypes caused incorrect estimations of the number of contributors. Although this method still has room for improvement, it may be useful for estimating the number of contributors in a mixed sample, as it does not rely on forensic mathematics.

Unraveling heteroplasmy patterns with NOVOPlasty

Heteroplasmy, the existence of multiple mitochondrial haplotypes within an individual, has been studied across different scientific fields. Mitochondrial genome polymorphisms have been linked to multiple severe disorders and are of interest to evolutionary studies and forensic science. Before the development of massive parallel sequencing (MPS), most studies of mitochondrial genome variation were limited to short fragments and to heteroplasmic variants associated with a relatively high frequency (>10%). By utilizing ultra-deep sequencing, it has now become possible to uncover previously undiscovered patterns of intra-individual polymorphisms. Despite these technological advances, it is still challenging to determine the origin of the observed intra-individual polymorphisms. We therefore developed a new method that not only detects intra-individual polymorphisms within mitochondrial and chloroplast genomes more accurately, but also looks for linkage among polymorphic sites by assembling the sequence around each detected polymorphic site. Our benchmark study shows that this method is capable of detecting heteroplasmy more accurately than any method previously available and is the first tool that is able to completely or partially reconstruct the sequence for each mitochondrial haplotype (allele). The method is implemented in our open source software NOVOPlasty that can be downloaded at https://github.com/ndierckx/NOVOPlasty.

Rare human mitochondrial HV lineages spread from the Near East and Caucasus during post-LGM and Neolithic expansions

Of particular significance to human population history in Eurasia are the migratory events that connected the Near East to Europe after the Last Glacial Maximum (LGM). Utilizing 315 HV*(xH,V) mitogenomes, including 27 contemporary lineages first reported here, we found the genetic signatures for distinctive movements out of the Near East and South Caucasus both westward into Europe and eastward into South Asia. The parallel phylogeographies of rare, yet widely distributed HV*(xH,V) subclades reveal a connection between the Italian Peninsula and South Caucasus, resulting from at least two (post-LGM, Neolithic) waves of migration. Many of these subclades originated in a population ancestral to contemporary Armenians and Assyrians. One such subclade, HV1b-152, supports a postexilic, northern Mesopotamian origin for the Ashkenazi HV1b2 lineages. In agreement with ancient DNA findings, our phylogenetic analysis of HV12 and HV14, the two exclusively Asian subclades of HV*(xH,V), point to the migration of lineages originating in Iran to South Asia before and during the Neolithic period. With HV12 being one of the oldest HV subclades, our results support an origin of HV haplogroup in the region defined by Western Iran, Mesopotamia, and the South Caucasus, where the highest prevalence of HV has been found.

Validation of NGS for mitochondrial DNA casework at the FBI Laboratory

As a first step towards integrating next generation sequencing (NGS) technology into the FBI Laboratory's operational casework, the PowerSeq™ CRM Nested System, an NGS-based mitochondrial DNA (mtDNA) control region assay, was developmentally and internally validated. The validation studies were conducted in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) Validation Guidelines for Forensic DNA Analysis Methods, and the FBI's Quality Assurance Standards (QAS) for Forensic DNA Testing Laboratories. The assay was shown to be highly reproducible, with variant frequencies across intra and inter-run replicates of the same sample differing, on average, by just 0.3% for substitutions and point heteroplasmies and 1.5% for insertions and deletions. The assay was also shown to be extremely sensitive, yielding complete control region sequence data from as few as 2000 copies of mtDNA. This is a more than 20-fold increase in sensitivity when compared to the FBI Laboratory's current Sanger sequencing-based protocols and, based on mtDNA quantitation values of samples routinely encountered in mtDNA casework, suggests that the percentage of questioned samples from which full control region data can be recovered will increase from our current 20% to approximately 90% success with NGS technology. In addition, the assay requires on average only 30% of the extract volume typically required to develop control region profiles from degraded samples via Sanger sequencing. Overall, these studies establish the reliability of the PowerSeq™ CRM Nested System for accurate mtDNA control region typing and can serve as a model for laboratories seeking to validate NGS protocols for forensic mtDNA analysis.

Antibiotic-Induced, Increased Conjugative Transfer Is Common to Diverse Naturally Occurring ESBL Plasmids in Escherichia coli

Previously, we showed that cefotaxime (CTX) exposure increases conjugative transfer of a bla CTX-M- 1 encoding IncI1 plasmid (IncI1/pST49/CTX-M-1) in Escherichia coli in a SOS-independent manner. This study aimed at investigating whether the observation was unique for that plasmid/strain/antibiotic combination or whether antibiotic-induced plasmid transfer (PT) is a more general phenomenon among plasmids in E. coli. Whole genome sequences of 25 E. coli strains were analyzed to identify different extended spectrum beta-lactamases (ESBL) plasmids enabling selection of a diverse collection of plasmids. Experiments were performed following exposure of these strains to 1/2 minimal inhibitory concentration (MIC) of CTX, ampicillin (AMP), or ciprofloxacin (CIP) before conjugation experiments. The frequency of PT was measured and compared to that of donors not exposed to antibiotics. Reverse-transcribed-quantitative real time polymerase chain reaction (RT-qPCR) was used to measure mRNA levels of five PT genes and two SOS response genes in donors exposed to antibiotics. The PT of eight strains (30.8% of strains tested) with IncI1/pST7/CTX-M-1, IncI1/pST49/CTX-M-1, IncI1/pST3/CTX-M-1, IncI1/pST293/CTX-M-1, IncI1/pST295/CTX-M-1, IncI1/pST16/CTX-M-55, and IncFII/CTX-M-14 (n = 2) plasmids was significantly increased following antibiotic exposure. CTX increased PT in all of these eight strain/plasmid combinations, AMP and CIP increased the PT in six and three strains, respectively. RT-qPCR showed that PT genes were up-regulated in the presence of the three antibiotics, whereas SOS-response genes were up-regulated only following CIP exposure. Our findings reveal that antibiotics can increase PT in E. coli strains with various ESBL plasmids. Thus, antibiotic-induced conjugative transfer of ESBL plasmids appears to be a common phenomenon in E. coli, having important implications for assessing the risks of antibiotic use.

Evaluation of mitogenome sequence concordance, heteroplasmy detection, and haplogrouping in a worldwide lineage study using the Precision ID mtDNA Whole Genome Panel

The emergence of Massively Parallel Sequencing technologies enabled the analysis of full mitochondrial (mt)DNA sequences from forensically relevant samples that have, so far, only been typed in the control region or its hypervariable segments. In this study, we evaluated the performance of a commercially available multiplex-PCR-based assay, the Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific), for the amplification and sequencing of the entire mitochondrial genome (mitogenome) from even degraded forensic specimens. For this purpose, more than 500 samples from 24 different populations were selected to cover the vast majority of established superhaplogroups. These are known to harbor different signature sequence motifs corresponding to their phylogenetic background that could have an effect on primer binding and, thus, could limit a broad application of this molecular genetic tool. The selected samples derived from various forensically relevant tissue sources and were DNA extracted using different methods. We evaluated sequence concordance and heteroplasmy detection and compared the findings to conventional Sanger sequencing as well as an orthogonal MPS platform. We discuss advantages and limitations of this approach with respect to forensic genetic workflow and analytical requirements.

Oral Administration of Lactobacillus curvatus HY7601 and Lactobacillus plantarum KY1032 with Cinnamomi Ramulus Extract Reduces Diet-Induced Obesity and Modulates Gut Microbiota

Obesity is a major health issue worldwide, and is associated with many diseases including type 2 diabetes. In this study, we evaluated the anti-obesity effects of combinations of two lactic acid bacteria (LAB), Lactobacillus curvatus HY7601 and Lactobacillus plantarum KY1032, and Cinnamomi Ramulus (CR) extract, and explored the mechanism through which they modulate gut microbiota using diet-induced obese mice. Male C57BL/6J mice were randomly divided into five groups that received a high-fat diet (HFD), HFD and LAB (HFD+LAB), HFD and CR extract (HFD+CR), HFD with LAB and CR extract (HFD+LAB+CR), or normal diet for 10 weeks. The mice in the HFD+LAB+CR group showed significant reductions in body weight gain, in particular epididymal fat and liver, blood leptin levels, and an increase in the levels of blood adiponectin. In addition, the LAB and CR extract altered the gut microbiota, mainly increasing the alpha diversity. These results demonstrate that a mixture of two LAB (Lactobacillus curvatus HY7601 and Lactobacillus plantarum KY1032) and CR extracts alleviate HFD-induced obesity, and has potential of being used as a strategy for the treatment of obesity.

Diversity of the Gut Microbiota in Dihydrotestosterone-Induced PCOS Rats and the Pharmacologic Effects of Diane-35, Probiotics, and Berberine

Polycystic ovary syndrome (PCOS) is a frequent endocrine and metabolic syndrome in reproductive-age women. Recently, emerging evidence has shown that gut microbiota is closely related to metabolic diseases such as type 2 diabetes, obesity and PCOS. In the present study, we established dihydrotestosterone (DHT)-induced PCOS rats and used Illumina MiSeq sequencing (PE300) to examine the composition, diversity, and abundance of the gut microbiota in PCOS. We compared the effects of three PCOS treatments: Diane-35 (estrogen and progesterone), probiotics and berberine. The DHT-induced rats showed constant estrous cycles, the loss of mature ovarian follicles, insulin resistance and obesity. The reproductive and metabolic functions in the PCOS rats were improved by treatment with Diane-35 and probiotics. Diane-35 and probiotics could restore the diversity of the gut microbiota, and the recovery of gut microbiota disorders improved the reproductive function in PCOS-like rats. However, berberine drastically reduced the species diversity and amount of gut microbiota and showed no improvement in PCOS. The findings of this study will help us to better understand the influence of the gut microbiota in the metabolic and reproductive alterations in PCOS as well as suggest opportunities for future personal dietary guidance for PCOS.

Bioinformatics Tools and Databases to Assess the Pathogenicity of Mitochondrial DNA Variants in the Field of Next Generation Sequencing

The development of next generation sequencing (NGS) has greatly enhanced the diagnosis of mitochondrial disorders, with a systematic analysis of the whole mitochondrial DNA (mtDNA) sequence and better detection sensitivity. However, the exponential growth of sequencing data renders complex the interpretation of the identified variants, thereby posing new challenges for the molecular diagnosis of mitochondrial diseases. Indeed, mtDNA sequencing by NGS requires specific bioinformatics tools and the adaptation of those developed for nuclear DNA, for the detection and quantification of mtDNA variants from sequence alignment to the calling steps, in order to manage the specific features of the mitochondrial genome including heteroplasmy, i.e., coexistence of mutant and wildtype mtDNA copies. The prioritization of mtDNA variants remains difficult, relying on a limited number of specific resources: population and clinical databases, and in silico tools providing a prediction of the variant pathogenicity. An evaluation of the most prominent bioinformatics tools showed that their ability to predict the pathogenicity was highly variable indicating that special efforts should be directed at developing new bioinformatics tools dedicated to the mitochondrial genome. In addition, massive parallel sequencing raised several issues related to the interpretation of very low mtDNA mutational loads, discovery of variants of unknown significance, and mutations unrelated to patient phenotype or the co-occurrence of mtDNA variants. This review provides an overview of the current strategies and bioinformatics tools for accurate annotation, prioritization and reporting of mtDNA variations from NGS data, in order to carry out accurate genetic counseling in individuals with primary mitochondrial diseases.

Mind the gap! The mitochondrial control region and its power as a phylogenetic marker in echinoids

BACKGROUND

In Metazoa, mitochondrial markers are the most commonly used targets for inferring species-level molecular phylogenies due to their extremely low rate of recombination, maternal inheritance, ease of use and fast substitution rate in comparison to nuclear DNA. The mitochondrial control region (CR) is the main non-coding area of the mitochondrial genome and contains the mitochondrial origin of replication and transcription. While sequences of the cytochrome oxidase subunit 1 (COI) and 16S rRNA genes are the prime mitochondrial markers in phylogenetic studies, the highly variable CR is typically ignored and not targeted in such analyses. However, the higher substitution rate of the CR can be harnessed to infer the phylogeny of closely related species, and the use of a non-coding region alleviates biases resulting from both directional and purifying selection. Additionally, complete mitochondrial genome assemblies utilizing next generation sequencing (NGS) data often show exceptionally low coverage at specific regions, including the CR. This can only be resolved by targeted sequencing of this region.

RESULTS

Here we provide novel sequence data for the echinoid mitochondrial control region in over 40 species across the echinoid phylogenetic tree. We demonstrate the advantages of directly targeting the CR and adjacent tRNAs to facilitate complementing low coverage NGS data from complete mitochondrial genome assemblies. Finally, we test the performance of this region as a phylogenetic marker both in the lab and in phylogenetic analyses, and demonstrate its superior performance over the other available mitochondrial markers in echinoids.

CONCLUSIONS

Our target region of the mitochondrial CR (1) facilitates the first thorough investigation of this region across a wide range of echinoid taxa, (2) provides a tool for complementing missing data in NGS experiments, and (3) identifies the CR as a powerful, novel marker for phylogenetic inference in echinoids due to its high variability, lack of selection, and high compatibility across the entire class, outperforming conventional mitochondrial markers.

Evaluation of the precision ID whole MtDNA genome panel for forensic analyses

Mitochondrial DNA (mtDNA) amplification and Massively Parallel Sequencing (MPS) using an early access version of the Precision ID Whole MtDNA Genome Panel (Thermo Fisher Scientific) and the Ion Personal Genome Machine (PGM) were evaluated using 15 forensically relevant samples. Samples were selected to represent typical forensic specimens for mtDNA analysis including hairs, hair shafts, swabs and ancient solid tissue samples (bones and teeth) that were stored in the freezer for up to several years after having been typed with conventional Sanger-type Sequencing and Capillary Electrophoresis. The MPS haplotypes confirmed the earlier results in all samples and provided additional sequence information that improved discrimination power and haplogroup estimation. The results raised the appetite for further experiments to validate and apply the new technology in forensic practice.

Deep-Coverage MPS Analysis of Heteroplasmic Variants within the mtGenome Allows for Frequent Differentiation of Maternal Relatives

Distinguishing between maternal relatives through mitochondrial (mt) DNA sequence analysis has been a longstanding desire of the forensic community. Using a deep-coverage, massively parallel sequencing (DCMPS) approach, we studied the pattern of mtDNA heteroplasmy across the mtgenomes of 39 mother-child pairs of European decent; haplogroups H, J, K, R, T, U, and X. Both shared and differentiating heteroplasmy were observed on a frequent basis in these closely related maternal relatives, with the minor variant often presented as 2–10% of the sequencing reads. A total of 17 pairs exhibited differentiating heteroplasmy (44%), with the majority of sites (76%, 16 of 21) occurring in the coding region, further illustrating the value of conducting sequence analysis on the entire mtgenome. A number of the sites of differentiating heteroplasmy resulted in non-synonymous changes in protein sequence (5 of 21), and to changes in transfer or ribosomal RNA sequences (5 of 21), highlighting the potentially deleterious nature of these heteroplasmic states. Shared heteroplasmy was observed in 12 of the 39 mother-child pairs (31%), with no duplicate sites of either differentiating or shared heteroplasmy observed; a single nucleotide position (16093) was duplicated between the data sets. Finally, rates of heteroplasmy in blood and buccal cells were compared, as it is known that rates can vary across tissue types, with similar observations in the current study. Our data support the view that differentiating heteroplasmy across the mtgenome can be used to frequently distinguish maternal relatives, and could be of interest to both the medical genetics and forensic communities.

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.

Assessment of mitochondrial DNA heteroplasmy detected on commercial panel using MPS system with artificial mixture samples

Mitochondrial DNA (mtDNA) heteroplasmy is a potential genetic marker for forensic mtDNA analysis as well as phylogenic studies. Frequency of mtDNA heteroplasmy has been investigated in different populations through massively parallel sequencing (MPS) analysis, revealing various levels of frequency based on different MPS systems. For accurate heteroplasmy identification, it is essential to explore reliable detection threshold on various MPS systems. In addition, software solutions and pipelines need to be evaluated to analyze data effectively. In this study, heteroplasmy analysis was conducted on a commercially available mtDNA analysis system developed for forensic caseworks with artificially mixed DNA samples known for ratios and variant positions for assessment. mtDNA heteroplasmy > 10% was detectable with Torrent Variant Caller (TVC) while lower levels were identified using GeneMarker® HTS specialized software for minor variant detection. This study implies that analytical parameters and tools need to be optimized and evaluated for low-level heteroplasmy identification. Automated system with simple and efficient workflow is needed for forensic caseworks.

A review of bioinformatic methods for forensic DNA analyses

Short tandem repeats, single nucleotide polymorphisms, and whole mitochondrial analyses are three classes of markers which will play an important role in the future of forensic DNA typing. The arrival of massively parallel sequencing platforms in forensic science reveals new information such as insights into the complexity and variability of the markers that were previously unseen, along with amounts of data too immense for analyses by manual means. Along with the sequencing chemistries employed, bioinformatic methods are required to process and interpret this new and extensive data. As more is learnt about the use of these new technologies for forensic applications, development and standardization of efficient, favourable tools for each stage of data processing is being carried out, and faster, more accurate methods that improve on the original approaches have been developed. As forensic laboratories search for the optimal pipeline of tools, sequencer manufacturers have incorporated pipelines into sequencer software to make analyses convenient. This review explores the current state of bioinformatic methods and tools used for the analyses of forensic markers sequenced on the massively parallel sequencing (MPS) platforms currently most widely used.