Sanger-type Sequencing Research Articles

FRI303 PRL-R Variants Are Not Only Associated With Prolactinomas But Also With Dopamine Agonist Resistance

Abstract Disclosure: A.R. Moreira: None. E.B. Trarbach: None. A. Glezer: None. G.A. Maciel: None. C.B. Bueno: None. A.S. Monteiro: None. M. Padula: None. I.P. Grande: None. M.D. Bronstein: None. Introduction: Knockout prolactin receptor gene (PRL-R) mice are animal models for prolactinomas and PRL acts via autocrine/paracrine inhibiting lactotroph proliferation. Recent studies identified variants in PRL-R gene in patients with prolactinoma and two of them Glu400(376)Gln/p.Asn516(492)Ile, were more common in patients compared to controls. Functional analysis of the p.Asn516(492)Ile variant showed an increase in Akt expression and cell proliferation induced by PRL. However, the association between PRL-R variants and response to cabergoline (CAB) treatment has not been investigated. Objective: This study aims to analyze the PRL-R gene in an extensive cohort of patients with sporadic prolactinomas and evaluate its association with clinical, laboratorial and tumor imaging characteristics, as well as CAB treatment resistance. Methods: DNA was extracted from peripheral blood samples of 178 patients with prolactinoma and 146 were classified according to cabergoline response. The PRL-R gene was analyzed by PCR techniques and automated Sanger-type sequencing. Four in silico tools were used for analysis. Results: We found six PRL-R variants: p.Ile100(76)Val (n=18), p.Ile170(146)Leu (n=6), p.Glu400(376)Gln/p.Asn516(492)Ile (n=3), p.Glu470Asp (n =1) and p.Ala591Pro (n=1), the last two described for the first time in prolactinomas´ patients. The variants p.Glu400(376)Gln/p.Asn516(492)Ile and p.Ala591Pro were more frequent amongst patients compared to genomic databases and the p.Asn516(492)Ile showed pathogenic potential using in silico analysis. PRL-R variants were associated with male sex (P=0.015), higher serum PRL levels (P=0.007), larger tumors (P=0.001), higher frequency of macroadenomas (P=0.031) and cabergoline resistance (P<0.001). Considering together only the PRL-R variants classified as pathogenic/probably pathogenic (p.Ile170(146)Leu, p.Glu400(376)Gln/Asn516(492) and p.Ala591Pro) versus the benign variants (p.Ile100(76)Val, p.Glu470Asp and wild type) we also observed a higher levels of serum PRL (P=0.008), largest tumor diameter (P=0.028) and higher frequency of CAB resistance (P=0.048).Conclusion: Our study showed a higher frequency of the PRL-R variants p.Glu400(376)Gln/p.Asn516(492)Ile, and p.Ala591Pro in patients with prolactinomas when compared to Brazilian and international databases, suggesting that these variants could play a role in the prolactinoma tumorigenesis. Patients with PRL-R variants presented male predominance, higher serum PRL levels, larger tumor size at diagnosis and higher frequency of CAB resistance. More studies evaluating germinative and somatic PRL-R variants are needed to clarify the PRL/PRL-R system role in the tumorigenesis and progression of prolactinomas and their interesting association with CAB resistance. Presentation: Friday, June 16, 2023

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.

Development of a multiplex assay for detection of autosomal and Y-chromosomal STRs, assessment of the degradation state of mitochondrial DNA and presence of mitochondrial length heteroplasmies

The current focus in most routine forensic casework is detection of autosomal or gonosomal Short Tandem Repeats (STRs). With increasing degradation, STR analysis tends to be less successful up to complete failure. For challenging samples such as telogen hair roots and shafts, touch DNA samples or skeletal remains, mitochondrial DNA (mtDNA) analysis provides a powerful tool. Determination of DNA quantity is an important part in the casework workflow. Several ready-to-use kits are commercially available for nuclear DNA targets. However, quantification of mtDNA targets requires the establishment of an in-house method. Some assays even contain assessment of degradation, which alleviates the choice of target enrichment for sequencing through medium or small amplicons. As Sanger-type Sequencing (STS) still remains the golden standard in many laboratories, identification of heteroplasmies in C-tract regions prior to the sequencing reaction is advantageous. Firstly, primer selection can be expanded with primers binding near the C-tract and secondly, determination of the dominant variant is straightforward. All those quantity (nuclear and mtDNA) and quality (degradation and length heteroplasmies) evaluations usually require at least two separate reactions. Therefore, the aim of this project was the combination of all these targets in one multiplex assay using capillary electrophoresis to spare valuable sample extract. Amplification of representative autosomal and Y-chromosomal STRs allows estimate of success of (Y-)STR analysis. Simultaneously, five length heteroplasmies in the mitochondrial control region are targeted as well as three conservative regions of differing fragment lengths for assessment of the mitochondrial degradation state. Based on the outcome of this assay, forensic examiners can decide if STR analysis may be suitable. In case of absent STR peaks, appropriate proceeding of mtDNA sequencing can be determined.

Clinical, biological, radiological, and genetic study of LPAC syndrome in Tunisian patients

Background and study aimsLow phospholipid-associated cholelithiasis (LPAC) syndrome is a form of cholelithiasis associated with the ABCB4 gene mutation. The defects of the protein ABCB4 encoded by this gene promote the formation of biliary cholesterol microcalculations. ABCB4 screening is negative in a significant proportion of patients. Patients and methodsAn analytical study of the epidemiological, clinical, biological, and radiological characteristics of 19 patients was conducted, followed by Sanger-type sequencing of the 27 exons encoding the ABCB4 gene. ResultsOur results showed a female predominance, symptomatic vesicular lithiasis predominance, and a high frequency of biliary complications in patients carrying an ABCB4 mutation. Normal ​​ liver enzyme values were found in 84.2% of the cases. Intrahepatic hyperechoic foci were present in 68.4%.Molecular analysis detected a pathogenic mutation of the ABCB4 gene in 31.57% of patients. The mutations found were a nonsense mutation and three missense mutations, including two new mutations. ConclusionOur epidemiological, clinical, and genetic results concord with previous studies of LPAC syndrome. Two of the mutations we found have never been detected in patients with LPAC. The low percentage of ABCB4 gene mutations can be explained by the absence of studies of other genes involved in bile acid homeostasis besides the ABCB4 gene and by the inclusion criteria used in this study.

Haplogroup Distribution of 309 Thais from Admixed Populations across the Country by HVI and HVII Sanger-Type Sequencing

The mitochondrial DNA (mtDNA) control region sequences for the hypervariable regions I (HVI) and II (HVII) of 309 Thai citizens were investigated using Sanger-type sequencing to generate an mtDNA reference dataset for forensic casework, and the haplogroup distribution within geographically proximal Asian populations was analyzed. The population sample set contained 264 distinct haplotypes and showed high haplotype diversity, low matching probability, and high powers of discrimination, at 0.9985, 0.4744%, and 0.9953, respectively, compared with previous reports. Subhaplogroup F1a showed the highest frequency in the Thai population, similar to Southeast Asian populations. The haplotype frequencies in the northern, northeastern, and southern populations of Thailand illustrate the relevance of social, religious, and historical factors in the biogeographical origin of the admixed Thai population as a whole. The HVI and HVII reference datasets will be useful for forensic casework applications, with improved genetic information content and discriminatory power compared to currently available techniques.

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.

The time is now for ubiquitous forensicmtMPSanalysis

AbstractMitochondrial (mt) DNA sequence analysis is useful for assessing ancestral origin and migration, identifying human remains, and examining evidentiary material in forensic casework. Conventional Sanger‐type sequencing (STS) has been used for more than three decades to address these interests. This paper reviews the methodologies and merits of using a massively parallel sequencing (MPS) approach for mtDNA testing in forensic laboratories, asThe Time is Now for Ubiquitous Forensic mtMPS Analysis.This article is categorized under:Forensic Biology > Haploid MarkersForensic Biology > Forensic DNA Technologies

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 FGxTM 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.

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.

Mitochondrial DNA control region variation in an Iraqi population sample.

Although Iraq has one of the world's oldest cultural histories and an important geographic location, forensic reference data on mitochondrial DNA (mtDNA) control region in Iraqi populations are scarce, particularly for populations residing in the southern part of Iraq. Mitochondrial DNA typing is an excellent tool for forensic investigations and in missing-person cases because of its unique qualities, such as mtDNA non-coding control region with specific genetic markers, high copy numbers in cells, maternal inheritance, and lack of recombination. Forensic analysis was performed on the entire mtDNA control region in 203 unrelated Iraqi individuals residing in Samawah City of Iraq. Polymorphisms in the mtDNA were detected using polymerase chain reaction and Sanger-type sequencing, and the sequences were aligned to compare with revised Cambridge Reference Sequence (rCRS). The sequencing results revealed 111 haplotypes characterized by 143 polymorphic positions. Of these haplotypes, 63 were unique and 48 were shared by more than one person. The haplotype data generated in this study will be available on EMPOP via accession number EMP00814.

Echinococcus granulosus genotypes verified in human hydatid disease around the world. Systematic review

The evidence regarding genotypic characteristics of Echinococcus granulosus infection in humans worldwide is scarce. To develop a synthesis of the available evidence regarding genotypes of E. granulosus verified in humans worldwide. Systematic review. Articles related with genotypes of E. granulosus, in humans, without language neither genotyped method restriction, published between 1990-2019 were included. A systematic in WoS, EMBASE, MEDLINE, SCOPUS, Trip Database, BIREME, SciELO, LILACS, IBECS, and PAHO-WHO was carried out. In study variables were year of publication, country, number of samples, host and parasite organs, genotype identified, molecular marker and genes. Descriptive statistics were applied. 701 related articles were identified; 62 fulfilled selection criteria, representing 1,511 samples. The existing evidence was published between 1994 and 2019; and mainly comes from Iran (45.2%). The most commonly used sequencing method was PCR amplification and Sanger type sequencing with partial or total genotyping of the cox1 gene. Genotyped method most frequently used was cox1 (79,0%). Genotypes most frequently identified were G1 and G1/G3 complex (49.1% and 32.2%). Publications related to genotypes of Eg in humans are scarce, heterogeneous, and presenting differing results. Eg G1/G3 accounts for most of the global burden worldwide.

Methods for the analysis of mitochondrial DNA

AbstractThe use of mitochondrial DNA (mtDNA) profiling as a tool for forming a link between an individual and a crime sample is a well‐established forensic biology technique. Though nuclear DNA typing is considered the gold standard method, certain evidence types have limited nuclear DNA available, such as telogen hairs and ancient or degraded remains, for which mtDNA analysis may be required. Traditionally mtDNA typing has been carried out using Sanger type sequencing, which is labor, time, and cost intensive. This has restricted the analysis of the mitochondrial genome to the control region, limiting the discrimination power of the technique. The introduction of massively parallel sequencing (MPS) into forensic laboratories has brought about the ability to analyze the whole mitochondrial genome of forensically relevant samples. The use of MPS for mtDNA analysis has resulted in an increase in the range of methods available for analysis, including: midi‐ and mini‐sized amplification strategies for whole genome analysis of casework‐type samples; whole genome single nucleotide polymorphism (SNP) multiplexes; and DNA capture methods allowing the analysis of highly degraded remains. The resulting increase in discriminatory power and the continued progression of the technology and associated methods increase its appeal as a DNA analysis method and may facilitate its implementation as a routinely used forensic tool. Here, we give an overview of current laboratory practices for the analysis of mtDNA using MPS.This article is categorized under: Forensic Biology > Haploid Markers Forensic Biology > Forensic DNA Technologies

First-tier detection of intragenomic 16S rRNA gene variation in culturable endophytic bacteria from cacao seeds.

BackgroundIntragenomic variability in 16S rDNA is a limiting factor for taxonomic and diversity characterization of Bacteria, and studies on its occurrence in natural/environmental populations are scarce. In this work, direct DNA amplicon sequencing coupled with frequent-cutter restriction analysis allowed detection of intragenomic 16S rDNA variation in culturable endophytic bacteria from cacao seeds in a fast and attractive manner.MethodsTotal genomic DNA from 65 bacterial strains was extracted and the 16S rDNA hyper variable V5–V9 regions were amplified for enzyme digestion and direct Sanger-type sequencing. The resulting electropherograms were visually inspected and compared to the corresponding AluI-restriction profiles, as well as to complete genome sequences in databases. Restriction analysis were employed to substitute the need of amplicon cloning and re-sequencing. A specifically improved polyacrylamide-gradient electrophoresis allowed to resolve 5-bp differences in restriction fragment sizes. Chi-square analysis on 2 × 2 contingency table tested for the independence between the ‘number of AluI bands’ and ‘type of eletropherogram’.ResultsTwo types of electropherograms were obtained: unique template, with single peaks per base (clean chromatograms), and heterogeneous template, with various levels of multiple peaks per base (mixed chromatograms). Statistics revealed significant interaction between number of restriction fragments and type of electropherogram for the same amplicons: clean or mixed ones associated to ≤5 or ≥6 bands, respectively. The mixed-template pattern combined with the AluI-restriction profiles indicated a high proportion of 49% of the culturable endophytes from a tropical environment showing evidence of intragenomic 16S rDNA heterogeneity.ConclusionThe approach presented here was useful for a rapid, first-tier detection of intragenomic variation in culturable isolates, which can be applied in studies of other natural populations; a preliminary view of intragenomic heterogeneity levels can complement culture-dependent and -independent methods. Consequences of these findings in taxonomic and diversity studies in complex bacterial communities are discussed.

Impact of the sequencing method on the detection and interpretation of mitochondrial DNA length heteroplasmy

Advancements in sequencing technologies allow for rapid and efficient analysis of mitochondrial DNA (mtDNA) in forensic laboratories, which is particularly beneficial for specimens with limited nuclear DNA. Next generation sequencing (NGS) offers higher throughput and sensitivity over traditional Sanger-type sequencing (STS) as well as the ability to quantitatively analyze the data. Changes in sample preparation, sequencing method and analysis required for NGS may alter the mtDNA haplotypes compared to previously generated STS data. Thus, the present study aimed to characterize the impact of different sequencing workflows on the detection and interpretation of length heteroplasmy (LHP), a particularly complicated aspect of mtDNA analysis. Whole mtDNA genome (mitogenome) data were generated for 16 high-quality samples using well-established Illumina and Ion methods, and the NGS data were compared to previously-generated STS mtDNA control region data. Although the mitogenome haplotypes were concordant with the exception of length and low-level variants (<30 % variant frequency), LHP in the hypervariable segment (HVS) polycytosine regions (C-tracts) differed across sequencing methods. Consistent with previous studies, LHP in HVS1 was observed in samples with nine or more consecutive cytosines (Cs) and eight Cs in the HVS2 region in the STS data. The Illumina data produced a similar pattern of LHP as the STS data, whereas the Ion data were noticeably different. More complex LHP (i.e. more length molecules) was observed in the Ion data, as length variation occurred in multiple homopolymer stretches within the targeted HVS regions. Further, the STS dominant or major molecule (MM) differed from the Ion MM in 11 (37 %) of the 30 regions evaluated and six instances (20 %) in Illumina data. This is of particular interest, as the MM is used by many forensic laboratories to report the HVS C-tract in the mtDNA haplotype. In general, the STS MMs were longer than the Illumina MMs, while the Ion MMs were the shortest. The higher rate of homopolymer indels in Ion data likely contributed to these differences. Supplemental analysis with alternative approaches demonstrated that the LHP pattern may also be altered by the bioinformatic tool and workflow used for data interpretation. The broader application of NGS in forensic laboratories will undoubtedly result in the use of varying sample preparation and sequencing methods. Based on these findings, minor LHP differences are expected across sequencing workflows, and it will be important that C-tract indels continue to be ignored for forensic queries and comparisons.

Sequencing of full mitochondrial genomes for NIST population samples

The U.S. National Institute of Standards and Technology in collaboration with the Armed Forces Medical Examiner System’s Armed Forces DNA Identification Laboratory has undertaken a population sequencing study of 607 mitochondrial genomes (mtGenomes). The goal of the study is to increase the availability of high-quality whole mtGenome sequences for haplotype frequency estimations, allowing the increased discrimination power of the entire extranuclear genome to be used. Whole mtGenome sequencing is enabled by next generation sequencing workflows, which are vastly simplified compared to Sanger-type sequencing. Here we evaluated the relative performance of mtGenome versus sequencing the control region of the mtGenome alone.

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.

Optimized mtDNA Control Region Primer Extension Capture Analysis for Forensically Relevant Samples and Highly Compromised mtDNA of Different Age and Origin.

The analysis of mitochondrial DNA (mtDNA) has proven useful in forensic genetics and ancient DNA (aDNA) studies, where specimens are often highly compromised and DNA quality and quantity are low. In forensic genetics, the mtDNA control region (CR) is commonly sequenced using established Sanger-type Sequencing (STS) protocols involving fragment sizes down to approximately 150 base pairs (bp). Recent developments include Massively Parallel Sequencing (MPS) of (multiplex) PCR-generated libraries using the same amplicon sizes. Molecular genetic studies on archaeological remains that harbor more degraded aDNA have pioneered alternative approaches to target mtDNA, such as capture hybridization and primer extension capture (PEC) methods followed by MPS. These assays target smaller mtDNA fragment sizes (down to 50 bp or less), and have proven to be substantially more successful in obtaining useful mtDNA sequences from these samples compared to electrophoretic methods. Here, we present the modification and optimization of a PEC method, earlier developed for sequencing the Neanderthal mitochondrial genome, with forensic applications in mind. Our approach was designed for a more sensitive enrichment of the mtDNA CR in a single tube assay and short laboratory turnaround times, thus complying with forensic practices. We characterized the method using sheared, high quantity mtDNA (six samples), and tested challenging forensic samples (n = 2) as well as compromised solid tissue samples (n = 15) up to 8 kyrs of age. The PEC MPS method produced reliable and plausible mtDNA haplotypes that were useful in the forensic context. It yielded plausible data in samples that did not provide results with STS and other MPS techniques. We addressed the issue of contamination by including four generations of negative controls, and discuss the results in the forensic context. We finally offer perspectives for future research to enable the validation and accreditation of the PEC MPS method for final implementation in forensic genetic laboratories.

Hairy matters: MtDNA quantity and sequence variation along and among human head hairs

Hairs from the same donor have been found to differ in mtDNA sequence within and among themselves and from other tissues, which impacts interpretation of results obtained in a forensic setting. However, little is known on the magnitude of this phenomenon and published data on systematic studies are scarce. We addressed this issue by generating mtDNA control region (CR) profiles of >450 hair fragments from 21 donors by Sanger-type sequencing (STS). To mirror forensic scenarios, we compared hair haplotypes from the same donors to each other, to the corresponding buccal swab reference haplotypes and analyzed several fragments of individual hairs. We also investigated the effects of hair color, donor sex and age, mtDNA haplogroup and chemical treatment on mtDNA quantity, amplification success and variation. We observed a wide range of individual CR sequence variation. The reference haplotype was the only or most common (≥75%) hair haplotype for most donors. However, in two individuals, the reference haplotype was only found in about a third of the investigated hairs, mainly due to differences at highly variable positions. Similarly, most hairs revealed the reference haplotype along their entire length, however, about a fifth of the hairs contained up to 71% of segments with deviant haplotypes, independent of the longitudinal position. Variation affected numerous positions, typically restricted to the individual hair and in most cases heteroplasmic, but also fixed (i.e. homoplasmic) substitutions were observed. While existing forensic mtDNA interpretation guidelines were found still sufficient for all comparisons to reference haplotypes, some comparisons between hairs from the same donor could yield false exclusions when those guidelines are strictly followed. This study pinpoints the special care required when interpreting mtDNA results from hair in forensic casework.

The mitochondrial DNA makeup of Romanians: A forensic mtDNA control region database and phylogenetic characterization

To evaluate the pattern of Romanian population from a mitochondrial perspective and to establish an appropriate mtDNA forensic database, we generated a high-quality mtDNA control region dataset from 407 Romanian subjects belonging to four major historical regions: Moldavia, Transylvania, Wallachia and Dobruja. The entire control region (CR) was analyzed by Sanger-type sequencing assays and the resulting 306 different haplotypes were classified into haplogroups according to the most updated mtDNA phylogeny. The Romanian gene pool is mainly composed of West Eurasian lineages H (31.7%), U (12.8%), J (10.8%), R (10.1%), T (9.1%), N (8.1%), HV (5.4%),K (3.7%), HV0 (4.2%), with exceptions of East Asian haplogroup M (3.4%) and African haplogroup L (0.7%).The pattern of mtDNA variation observed in this study indicates that the mitochondrial DNA pool is geographically homogeneous across Romania and that the haplogroup composition reveals signals of admixture of populations of different origin. The PCA scatterplot supported this scenario, with Romania located in southeastern Europe area, close to Bulgaria and Hungary, and as a borderland with respect to east Mediterranean and other eastern European countries.High haplotype diversity (0.993) and nucleotide diversity indices (0.00838±0.00426), together with low random match probability (0.0087) suggest the usefulness of this control region dataset as a forensic database in routine forensic mtDNA analysis and in the investigation of maternal genetic lineages in the Romanian population.

Concordance and reproducibility of a next generation mtGenome sequencing method for high-quality samples using the Illumina MiSeq

Sanger-type sequencing (STS) of mitochondrial DNA (mtDNA), specifically the control region (CR), is routinely employed in forensics in human identification and missing persons scenarios. Yet next-generation sequencing (NGS) has the potential to overcome some of the major limitations of STS processing, permitting reasonable paths forward for full mitochondrial genome (mtGenome) sequencing, while also offering higher-throughput and higher sensitivity capabilities. To establish the accuracy and reproducibility of NGS for the development of mtDNA data, 90 DNA extracts that were previously used to generate forensic quality full mtGenomes using STS were sequenced using Nextera XT library preparation and the Illumina MiSeq. Using the same amplicon product, replicate library sets were generated and sequenced at different laboratories, and analysis was performed in replicate using the CLC Genomics Workbench. Both sequencing sets resulted in 99.998% of positions with greater than 10X coverage when 96 samples (including controls) were multiplexed. Overall, 99.9996% concordance was observed between the NGS data and the STS data for the full mtGenome. The only “discordant” calls involved low level point heteroplasmies, with the differences resulting from stochastic variation and/or the increased sensitivity of NGS. Higher sensitivity also allowed for the detection of a mixed sample previously not detected with STS. Additionally, variant calls were reproducible between sequencing sets and between software analysis versions with the variant frequency only differing by 0.23% and 0.01%, respectively. Further validation studies and specialized software functionality tailored to forensic practice should facilitate the incorporation of NGS processing into standard casework applications. The data herein comprise the largest, and likely most thoroughly examined, complete mtGenome STS-NGS concordance dataset available.