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

Assessment of ForenSeq mtDNA Whole Genome Kit for forensic application

Mitochondrial DNA (mtDNA) is an indispensable genetic marker in forensic genetics. The emergence and development of massively parallel sequencing (MPS) makes it possible to obtain complete mitochondrial genome sequences more quickly and accurately. The study evaluated the advantages and limitations of the ForenSeq mtDNA Whole Genome Kit in the practical application of forensic genetics by detecting human genomic DNA standards and thirty-three case samples. We used control DNA with different amount to determine sensitivity of the assay. Even when the input DNA is as low as 2.5 pg, most of the mitochondrial genome sequences could still be covered. For the detection of buccal swabs and aged case samples (bloodstains, bones, teeth), most samples could achieve complete coverage of mitochondrial genome. However, when ancient samples and hair samples without hair follicles were sequenced by the kit, it failed to obtain sequence information. In general, the ForenSeq mtDNA Whole Genome Kit has certain applicability to forensic low template and degradation samples, and these results provide the data basis for subsequent forensic applications of the assay. The overall detection process and subsequent analysis are easy to standardize, and it has certain application potential in forensic cases.

Comparison of Two Variant Analysis Programs for Next-Generation Sequencing Data of Whole Mitochondrial Genome

BACKGROUND

With advance of next-generation sequencing (NGS) techniques, the need for mitochondrial DNA analysis is increasing not only in the forensic area, but also in medical fields.

METHODS

Two commercial programs, Converge Software (CS) and Torrent Variant Caller for variant calling of NGS data, were compared with a considerable amount of sequence data of 50 samples with a homogeneous ethnicity.

RESULTS

About 2,300 variants were identified and the two programs showed about 90% of consistency. CS, a dedicated analysis program for mitochondrial DNA, showed some advantages for forensic use. By additional visual inspection, several causes of discrepancy in variant calling results were identified. Application of different notation rules for mitochondrial sequence and the minor allele frequency close to detection threshold were the two most significant reasons.

CONCLUSION

With prospective improvement of each program, researchers and practitioners should be aware of characteristics of the analysis program they use and prepare their own strategies to determine variants.

NuMY-A qPCR Assay Simultaneously Targeting Human Autosomal, Y-Chromosomal, and Mitochondrial DNA

The accurate quantification of DNA in forensic samples is of utmost importance. These samples are often present in limited amounts; therefore, it is indicated to use the appropriate analysis route with the optimum DNA amount (when possible). Also, DNA quantification can inform about the degradation stage and therefore support the decision on which downstream genotyping method to use. Consequently, DNA quantification aids in getting the best possible results from a forensic sample, considering both its DNA quantity and quality limitations. Here, we introduce NuMY, a new quantitative real-time PCR (qPCR) method for the parallel quantification of human nuclear (n) and mitochondrial (mt) DNA, assessing the male portion in mixtures of both sexes and testing for possible PCR inhibition. NuMY is based on previous work and follows the MIQE guidelines whenever applicable. Although quantification of nuclear (n)DNA by simultaneously analyzing autosomal and male-specific targets is available in commercial qPCR kits, tools that include the quantification of mtDNA are sparse. The quantification of mtDNA has proven relevant for samples with low nDNA content when conventional DNA fingerprinting techniques cannot be followed. Furthermore, the development and use of new massively parallel sequencing assays that combine multiple marker types, i.e., autosomal, Y-chromosomal, and mtDNA, can be optimized when precisely knowing the amount of each DNA component present in the input sample. For high-quality DNA extracts, NuMY provided nDNA results comparable to those of another quantification technique and has also proven to be a reliable tool for challenging, forensically relevant samples such as mixtures, inhibited, and naturally degraded samples.

Mitochondrial DNA in Human Diversity and Health: From the Golden Age to the Omics Era

Mitochondrial DNA (mtDNA) is a small fraction of our hereditary material. However, this molecule has had an overwhelming presence in scientific research for decades until the arrival of high-throughput studies. Several appealing properties justify the application of mtDNA to understand how human populations are-from a genetic perspective-and how individuals exhibit phenotypes of biomedical importance. Here, I review the basics of mitochondrial studies with a focus on the dawn of the field, analysis methods and the connection between two sides of mitochondrial genetics: anthropological and biomedical. The particularities of mtDNA, with respect to inheritance pattern, evolutionary rate and dependence on the nuclear genome, explain the challenges of associating mtDNA composition and diseases. Finally, I consider the relevance of this single locus in the context of omics research. The present work may serve as a tribute to a tool that has provided important insights into the past and present of humankind.

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

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

DNA and protein analyses of hair in forensic genetics

Hair is one of the most common pieces of biological evidence found at a crime scene and plays an essential role in forensic investigation. Hairs, especially non-follicular hairs, are usually found at various crime scenes, either by natural shedding or by forcible shedding. However, the genetic material in hairs is usually highly degraded, which makes forensic analysis difficult. As a result, the value of hair has not been fully exploited in forensic investigations and trials. In recent years, with advances in molecular biology, forensic analysis of hair has achieved remarkable strides and provided crucial clues in numerous cases. This article reviews recent developments in DNA and protein analysis of hair and attempts to provide a comprehensive solution to improve forensic hair analysis.

Cumulative UV Exposure or a Modified SCINEXA™-Skin Aging Score Do Not Play a Substantial Role in Predicting the Risk of Developing Keratinocyte Cancers after Solid Organ Transplantation-A Case Control Study

The risk of keratinocyte cancer is determined by intrinsic and extrinsic factors, which also influence skin aging. Few studies have linked skin aging and UV exposure with the incidence of non-melanoma skin cancer (NMSC). We evaluated signs of actinic skin damage and aging, individual UV burden, and melanocortin-1 receptor (MC1R) variants. A total of 194 organ transplant recipients (OTR) who suffered from NMSC were compared to 194 tumor-free controls matched for gender, age, type of transplanted organ, post-transplantation (TX) period, and immunosuppressive therapy. Compared with the cases, the controls scored higher in all skin aging scores and there were no differences in UV burden except for intentional whole-body UV exposure for specific UV scenarios and periods of life in favor of cases. The number of NMSCs correlated with all types of skin aging scores, the extent of intentional sun exposure, older age, longer post-TX period, shorter interval from TX to first NMSC, and specific MC1R risk groups. Multivariable models revealed a 7.5-fold risk of developing NMSC in individuals with actinic keratosis; 4.1- or 3.6-fold in those with green or blue eyes, respectively; and a 1.9-fold increased risk in the MC1R medium- + high-risk group. In the absence of skin aging contributing to NMSC development, certain MC1R risk types may identify OTR at risk for high tumor burden.

Recent advances in forensic biology and forensic DNA typing: INTERPOL review 2019-2022

This review paper covers the forensic-relevant literature in biological sciences from 2019 to 2022 as a part of the 20th INTERPOL International Forensic Science Managers Symposium. Topics reviewed include rapid DNA testing, using law enforcement DNA databases plus investigative genetic genealogy DNA databases along with privacy/ethical issues, forensic biology and body fluid identification, DNA extraction and typing methods, mixture interpretation involving probabilistic genotyping software (PGS), DNA transfer and activity-level evaluations, next-generation sequencing (NGS), DNA phenotyping, lineage markers (Y-chromosome, mitochondrial DNA, X-chromosome), new markers and approaches (microhaplotypes, proteomics, and microbial DNA), kinship analysis and human identification with disaster victim identification (DVI), and non-human DNA testing including wildlife forensics. Available books and review articles are summarized as well as 70 guidance documents to assist in quality control that were published in the past three years by various groups within the United States and around the world.

Whole mitochondrial genome analysis of Tai-Kadai-speaking populations in Southwest China

As a single matrilineal gene, human mitochondrial DNA plays a very important role in the study of population genetics. The whole mitogenome sequences of 287 individuals of the Tai-Kadai-speaking population in Guizhou were obtained. It was discovered that there were 82, 104, and 94 haplotypes in 83 Bouyei individuals, 107 Dong individuals, and 97 Sui individuals, respectively; and the haplotype diversity in Bouyei, Dong, and Sui groups was 1.000 ± 0.02, 0.9993 ± 0.0015, and 0.999 ± 0.002, respectively. The result of neutrality tests of the Tai-Kadai-speaking population in Guizhou showed significant negative values, and the analysis of mismatch distribution showed an obvious unimodal distribution. The results implied that Guizhou Tai-Kadai-speaking populations had high genetic diversities and may have experienced recent population expansion. In addition, the primary haplogroups of studied populations were M*, F, B, D, and R*, implying that they may origin from Southern China. The matrilineal genetic structure of the Tai-Kadai-speaking populations in Guizhou was analyzed by merging the mitogenome data of 79 worldwide populations as reference data. The results showed that there were close relationships between studied populations and other Tai-Kadai as well as some Austronesian populations in East and Southeast Asia. Overall, the mitogenome data generated in this study will provide important data for the study of genetic structure of Tai-Kadai speaking populations.

Helena's Many Daughters: More Mitogenome Diversity behind the Most Common West Eurasian mtDNA Control Region Haplotype in an Extended Italian Population Sample

The high number of matching haplotypes of the most common mitochondrial (mt)DNA lineages are considered to be the greatest limitation for forensic applications. This study investigates the potential to solve this constraint by massively parallel sequencing a large number of mitogenomes that share the most common West Eurasian mtDNA control region (CR) haplotype motif (263G 315.1C 16519C). We augmented a pilot study on 29 to a total of 216 Italian mitogenomes that represents the largest set of the most common CR haplotype compiled from a single country. The extended population sample confirmed and extended the huge coding region diversity behind the most common CR motif. Complete mitogenome sequencing allowed for the detection of 163 distinct haplotypes, raising the power of discrimination from 0 (CR) to 99.6% (mitogenome). The mtDNAs were clustered into 61 named clades of haplogroup H and did not reveal phylogeographic trends within Italy. Rapid individualization approaches for investigative purposes are limited to the most frequent H clades of the dataset, viz. H1, H3, and H7.

Genetic Background of Kirgiz Ethnic Group From Northwest China Revealed by Mitochondrial DNA Control Region Sequences on Massively Parallel Sequencing

The mitochondrial DNA (mtDNA) has been used to trace population evolution and apply to forensic identification due to the characteristics including lack of recombination, higher copy number and matrilineal inheritance comparing with nuclear genome DNA. In this study, mtDNA control region sequences of 91 Kirgiz individuals from the Northwest region of China were sequenced to identify genetic polymorphisms and gain insight into the genetic background of the Kirgiz ethnic group. MtDNA control region sequences of Kirgiz individuals presented relatively high genetic polymorphisms. The 1,122 bp sequences of mtDNA control region could differ among unrelated Kirgiz individuals, which suggested the mtDNA control region sequences have a good maternal pedigree tracing capability among different Kirgiz individuals. The neutrality test, mismatch distribution, Bayesian phylogenetic inference, Bayesian skyline analysis, and the median network analyses showed that the Kirgiz group might occurred population expansion, and the expansion could be observed at about ∼53.41 kilo years ago (kya) when ancestries of modern humans began to thrive in Eurasia. The pairwise population comparisons, principal component analyses, and median network analyses were performed based on haplogroup frequencies or mtDNA control region sequences of 5,886 individuals from the Kirgiz group and the 48 reference populations all over the world. And the most homologous haplotypes were found between Kirgiz individuals and the East Asian individuals, which indicated that the Kirgiz group might have gene exchanges with the East Asian populations.

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 Ancestry of Eastern Paraguay: A Typical South American Profile with a Unique Pattern of Admixture

Immigrants from diverse origins have arrived in Paraguay and produced important demographic changes in a territory initially inhabited by indigenous Guarani. Few studies have been performed to estimate the proportion of Native ancestry that is still preserved in Paraguay and the role of females and males in admixture processes. Therefore, 548 individuals from eastern Paraguay were genotyped for three marker sets: mtDNA, Y-SNPs and autosomal AIM-InDels. A genetic homogeneity was found between departments for each set of markers, supported by the demographic data collected, which showed that only 43% of the individuals have the same birthplace as their parents. The results show a sex-biased intermarriage, with higher maternal than paternal Native American ancestry. Within the native mtDNA lineages in Paraguay (87.2% of the total), most haplogroups have a broad distribution across the subcontinent, and only few are concentrated around the Paraná River basin. The frequency distribution of the European paternal lineages in Paraguay (92.2% of the total) showed a major contribution from the Iberian region. In addition to the remaining legacy of the colonial period, the joint analysis of the different types of markers included in this study revealed the impact of post-war migrations on the current genetic background of Paraguay.

From Forensics to Clinical Research: Expanding the Variant Calling Pipeline for the Precision ID mtDNA Whole Genome Panel

Despite a multitude of methods for the sample preparation, sequencing, and data analysis of mitochondrial DNA (mtDNA), the demand for innovation remains, particularly in comparison with nuclear DNA (nDNA) research. The Applied Biosystems™ Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific, USA) is an innovative library preparation kit suitable for degraded samples and low DNA input. However, its bioinformatic processing occurs in the enterprise Ion Torrent Suite™ Software (TSS), yielding BAM files aligned to an unorthodox version of the revised Cambridge Reference Sequence (rCRS), with a heteroplasmy threshold level of 10%. Here, we present an alternative customizable pipeline, the PrecisionCallerPipeline (PCP), for processing samples with the correct rCRS output after Ion Torrent sequencing with the Precision ID library kit. Using 18 samples (3 original samples and 15 mixtures) derived from the 1000 Genomes Project, we achieved overall improved performance metrics in comparison with the proprietary TSS, with optimal performance at a 2.5% heteroplasmy threshold. We further validated our findings with 50 samples from an ongoing independent cohort of stroke patients, with PCP finding 98.31% of TSS's variants (TSS found 57.92% of PCP's variants), with a significant correlation between the variant levels of variants found with both pipelines.

Internal validation and improvement of mitochondrial genome sequencing using the Precision ID mtDNA Whole Genome Panel

With the recent advances in next-generation sequencing (NGS), mitochondrial whole-genome sequencing has begun to be applied to the field of the forensic biology as an alternative to the traditional Sanger-type sequencing (STS). However, experimental workflows, commercial solutions, and output data analysis must be strictly validated before being implemented into the forensic laboratory. In this study, we performed an internal validation for an NGS-based typing of the entire mitochondrial genome using the Precision ID mtDNA Whole Genome Panel (Thermo Fisher Scientific) on the Ion S5 sequencer (Thermo Fisher Scientific). Concordance, repeatability, reproducibility, sensitivity, and heteroplasmy detection analyses were assessed using the 2800 M and 9947A standard control DNA as well as typical casework specimens, and results were compared with conventional Sanger sequencing and another NGS sequencer in a different laboratory. We discuss the strengths and limitations of this approach, highlighting some issues regarding noise thresholds and heteroplasmy detection, and suggesting solutions to mitigate these effects and improve overall data interpretation. Results confirmed that the Precision ID Whole mtDNA Genome Panel is highly reproducible and sensitive, yielding useful full mitochondrial DNA sequences also from challenging DNA specimens, thus providing further support for its use in forensic practice.

Massively parallel sequencing of human skeletal remains in Vietnam using the precision ID mtDNA control region panel on the Ion S5™ system

Mitochondrial DNA (mtDNA) analysis using Sanger sequencing has been a routine practice for the identification of human skeletal remains. However, this process is usually challenging since DNA from the remains is highly degraded and at low concentration. Recently, the advent and implementation of massively parallel sequencing (MPS) have been offered the ability to improve mtDNA sequence data for forensic analysis. To assess the utility of the Ion S5™ system - an MPS platform for mtDNA analysis in challenging samples, we sequenced the mitochondrial control region of 52 age-old skeletal remains. Using the Precision ID mtDNA Control Region Panel, 50 full and two partial control region haplotypes at relatively high mean coverage of 2494 × were achieved for variant calling. Further variant analysis at 10% threshold for point heteroplasmy showed high degradation degree in terms of DNA damage in our bone samples. A higher point heteroplasmy threshold of 20% was required to diminish most of background noise caused by the damage. The results from this study indicated the potential application of the Ion S5™ system in sequencing degraded samples in Vietnam and provided valuable data sources for forensic analyses in the future.

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.

Evidence for multi-copy Mega-NUMTs in the human genome

The maternal mode of mitochondrial DNA (mtDNA) inheritance is central to human genetics. Recently, evidence for bi-parental inheritance of mtDNA was claimed for individuals of three pedigrees that suffered mitochondrial disorders. We sequenced mtDNA using both direct Sanger and Massively Parallel Sequencing in several tissues of eleven maternally related and other affiliated healthy individuals of a family pedigree and observed mixed mitotypes in eight individuals. Cells without nuclear DNA, i.e. thrombocytes and hair shafts, only showed the mitotype of haplogroup (hg) V. Skin biopsies were prepared to generate ρ° cells void of mtDNA, sequencing of which resulted in a hg U4c1 mitotype. The position of the Mega-NUMT sequence was determined by fluorescence in situ hybridization and two different quantitative PCR assays were used to determine the number of contributing mtDNA copies. Thus, evidence for the presence of repetitive, full mitogenome Mega-NUMTs matching haplogroup U4c1 in various tissues of eight maternally related individuals was provided. Multi-copy Mega-NUMTs mimic mixtures of mtDNA that cannot be experimentally avoided and thus may appear in diverse fields of mtDNA research and diagnostics. We demonstrate that hair shaft mtDNA sequencing provides a simple but reliable approach to exclude NUMTs as source of misleading results.

A Continuous Statistical Phasing Framework for the Analysis of Forensic Mitochondrial DNA Mixtures

Despite the benefits of quantitative data generated by massively parallel sequencing, resolving mitotypes from mixtures occurring in certain ratios remains challenging. In this study, a bioinformatic mixture deconvolution method centered on population-based phasing was developed and validated. The method was first tested on 270 in silico two-person mixtures varying in mixture proportions. An assortment of external reference panels containing information on haplotypic variation (from similar and different haplogroups) was leveraged to assess the effect of panel composition on phasing accuracy. Building on these simulations, mitochondrial genomes from the Human Mitochondrial DataBase were sourced to populate the panels and key parameter values were identified by deconvolving an additional 7290 in silico two-person mixtures. Finally, employing an optimized reference panel and phasing parameters, the approach was validated with in vitro two-person mixtures with differing proportions. Deconvolution was most accurate when the haplotypes in the mixture were similar to haplotypes present in the reference panel and when the mixture ratios were neither highly imbalanced nor subequal (e.g., 4:1). Overall, errors in haplotype estimation were largely bounded by the accuracy of the mixture's genotype results. The proposed framework is the first available approach that automates the reconstruction of complete individual mitotypes from mixtures, even in ratios that have traditionally been considered problematic.

Developmental Validation of a MPS Workflow with a PCR-Based Short Amplicon Whole Mitochondrial Genome Panel

For the adoption of massively parallel sequencing (MPS) systems by forensic laboratories, validation studies on specific workflows are needed to support the feasibility of implementation and the reliability of the data they produce. As such, the whole mitochondrial genome sequencing methodology—Precision ID mtDNA Whole Genome Panel, Ion Chef, Ion S5, and Converge—has been subjected to a variety of developmental validation studies. These validation studies were completed in accordance with the Scientific Working Group on DNA Analysis Methods (SWGDAM) validation guidelines and assessed reproducibility, repeatability, accuracy, sensitivity, specificity to human DNA, and ability to analyze challenging (e.g., mixed, degraded, or low quantity) samples. Intra- and inter-run replicates produced an average maximum pairwise difference in variant frequency of 1.2%. Concordance with data generated with traditional Sanger sequencing and an orthogonal MPS platform methodology was used to assess accuracy, and generation of complete and concordant haplotypes at DNA input levels as low as 37.5 pg of nuclear DNA or 187.5 mitochondrial genome copies illustrated the sensitivity of the system. Overall, data presented herein demonstrate that highly accurate and reproducible results were generated for a variety of sample qualities and quantities, supporting the reliability of this specific whole genome mitochondrial DNA MPS system for analysis of forensic biological evidence.

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.

Distinguishing mitochondrial DNA and NUMT sequences amplified with the precision ID mtDNA whole genome panel

Nuclear mitochondrial DNA segments (NUMTs) are generated via transfer of portions of the mitochondrial genome into the nuclear genome. Given their common origin, there is the possibility that both the mitochondrial and NUMT segments may co-amplify using the same set of primers. Thus, analysis of the variation of the mitochondrial genome must take into account this co-amplification of mitochondrial and NUMT sequences. The study herein builds on data from the study by Strobl et al. (Strobl et al., 2019), in which multiple point heteroplasmies were called with an "N" to prevent labeling NUMT sequences mimicking mitochondrial heteroplasmy and being interpreted as true mitochondrial in origin sequence variants. Each of these point heteroplasmies was studied in greater detail, both molecularly and bioinformatically, to determine whether NUMT or true mitochondrial DNA variation was present. The bioinformatic and molecular tools available to help distinguish between NUMT and mitochondrial DNA and the effect of NUMT sequences on interpretation were discussed.

Numt identification and removal with RtN!

Motivation

Assays in mitochondrial genomics rely on accurate read mapping and variant calling. However, there are known and unknown nuclear paralogs that have fundamentally different genetic properties than that of the mitochondrial genome. Such paralogs complicate the interpretation of mitochondrial genome data and confound variant calling.

Results

Remove the Numts! (RtN!) was developed to categorize reads from massively parallel sequencing data not based on the expected properties and sequence identities of paralogous nuclear encoded mitochondrial sequences, but instead using sequence similarity to a large database of publicly available mitochondrial genomes. RtN! removes low-level sequencing noise and mitochondrial paralogs while not impacting variant calling, while competing methods were shown to remove true variants from mitochondrial mixtures.

Availability and implementation

https://github.com/Ahhgust/RtN

Supplementary information

Supplementary data are available at Bioinformatics online.

The lot-to-lot variability in the mitochondrial genome of controls

Current research in the biomedical field has illustrated how cell lines used as reference standards can change over time and, more importantly, can affect research and diagnostic results obtained from these cell lines. With the use of increasingly sensitive and highly resolving technologies (e.g., massively parallel sequencing), forensic scientists must be aware of and account for potential variability in the cell lines used as controls in their validation studies and day-to-day casework. In this study, multiple lot numbers from four commonly-used control cell line DNAs were sequenced with massively parallel sequencing on the Ion S5. The variability among these different lots was evaluated, and the effect on forensic laboratory work discussed.

Massively parallel sequencing of mitogenome sequences reveals the forensic features and maternal diversity of tai-kadai-speaking hlai islanders

As a single maternally inherited locus, human mitochondrial DNA (mtDNA) is geographically arranged and plays a key role in forensic applications. Hlai population has been evidenced as the most typical and unmixed representative of the Tai-Kadai-speaking populations via genome-wide analyses. However, forensic features and maternal diversity of the complete mitogenomes in this Tai-Kadai ancestrally related population are scarce. Thus, we sequenced the complete mitogenomes in 127 Hainan Hlais and found 109 distinct haplotypes belonging to 43 terminal haplogroups resulting in the haplotype diversity of 0.9970. Our results of comprehensive population comparisons showed that Hlai islanders had a close genetic affinity with Tai-Kadai-speaking populations from Southeast Asia, which is consistent with the back-migration of Chinese Neolithic farmers into this region via the inland route. Besides, maternally genetic evidence further revealed a close genetic relationship between Tai-Kadai-speaking and Austronesian-speaking populations when only East Asian dataset was considered, which is consistent with the common origin from Yangtze rice farmers and then spread southward along the inland and coastal routes, respectively. In the reconstructed phylogenetic tree and median-joining networks, the vast majority of Hlais were clustered in exclusive clades, which demonstrated that Hlai people probably had undergone founder effect or genetic bottleneck in their history, and remained genetically isolated for a long time. Collectively, Hainan Hlai did not exhibit detectable maternal gene flow from surrounding or incoming populations. Mitogenome information generated in this study is a contribution in mitigating the underrepresentation of Chinese data in forensic mitogenetics and will assist geography-, metapopulation-, as well as phylogeny-based queries.

Forensically relevant phylogeographic evaluation of mitogenome variation in the Basque Country

The Basque Country has been the focus of population (genetic) and evolutionary studies for decades, as it represents an interesting evolutionary feature: it is the only European country where a non-Indo-European language is still spoken today and, for which there are no known living or extinct relatives. Early studies that were based on anatomical and serological methods, along with subsequent molecular genetic investigations, contain controversial interpretations of their data. Additionally, the analysis of mitochondrial DNA, which is maternally inherited and thus suitable for the examination of the maternal phylogeny of the population, was the focus of some studies. Early mtDNA studies were however restricted to the information provided by the control region or its hypervariable segments only. These are known to harbour little phylogenetic information, particularly for haplogroup H that is dominant in Westeurasian populations including the Basques. Later studies analysed complete mitogenome sequences. Their information content is however limited, either because the number of samples was low, or because these studies only considered particular haplogroups. In this study we present the full mitogenome sequences of 178 autochthonous Basque individuals that were carefully selected based on their familial descent and discuss the observed phylogenetic signals in the light of earlier published findings. We confirm the presence of Basque-specific mtDNA lineages and extend the knowledge of these lineages by providing data on their distribution in comparison to other Basque and non-Basque populations. This dataset improves our understanding of the Basque mtDNA phylogeny and serves as a high-quality dataset that is provided via EMPOP for forensic genetic purposes.