A call for mtDNA data quality control in forensic science

Behind every good research there are data. What are they and their importance to forensic science

Data underpinning science have become one of the most precious assets in research, and while the principles of FAIR (Findable, Accessible, Interoperable and Reusable) have been put forward as a guide to how to approach data handling, data sharing and long-term storage still remain a challenge for many research areas including forensic science. The reporting and the sharing of data can be made easier by giving them structure, the use of suitable labels and the inclusion of descriptors collated into metadata prior to their deposition in repositories with persistent identifiers. Such a systematic approach would strengthen the quality and the integrity of research while providing greater transparency to published materials.

Cranio-facial bones evaluation based on clinical CT data for sex determination in Northwest Indian population

The skull and pelvis have been the first choice of bones for determination of unknown human remains. The goal of the present study was to derive discriminant function equations by using clinical CT scan data of cranio-facial bones for sex determination in Northwest Indian population. This study was conducted at Department of Radiology, by collecting the retrospective data of CT scan of 217 samples. In the data, 106 were males and 111 were females in the age group between 20 and 80 years. The total number of parameters under investigation were 10. All the selected variables were sexually dimorphic and showed significant values. 91.7% of original grouped cases were correctly classified to their sex category. The TEM, rTEM and R were under the acceptable limits. The univariate, multivariate and stepwise discriminant function analysis recorded an accuracy of 88.9%, 91.7% and 93.6% respectively. Multivariate direct discriminant function analysis stepwise method yielded the highest level of accuracy in differentiating males and females. All the variables reflected statistically significant difference between males and females (p less than 0.001). The best single parameter with highest level of sexual dimorphic trait was length of cranial base. This study aims to provide sex assessment using clinical data of CT scan in Northwest Indian population by incorporating the BIOFB cranio-facial parameter. The morphometric measurements taken on CT scan images can be utilized by forensic experts in identification process.

The Revolution of Animal Genomics in Forensic Sciences

Nowadays, the coexistence between humans and domestic animals (especially dogs and cats) has become a common scenario of daily life. Consequently, during a forensic investigation in civil or criminal cases, the biological material from a domestic animal could be considered "evidence" by law enforcement agencies. Animal genomics offers an important contribution in attacks and episodes of property destruction or in a crime scene where the non-human biological material is linked to the victim or perpetrator. However, only a few animal genetics laboratories in the world are able to carry out a valid forensic analysis, adhering to standards and guidelines that ensure the admissibility of data before a court of law. Today, forensic sciences focus on animal genetics considering all domestic species through the analysis of STRs (short tandem repeats) and autosomal and mitochondrial DNA SNPs (single nucleotide polymorphisms). However, the application of these molecular markers to wildlife seems to have gradually gained a strong relevance, aiming to tackle illegal traffic, avoid the loss of biodiversity, and protect endangered species. The development of third-generation sequencing technologies has glimmered new possibilities by bringing "the laboratory into the field", with a reduction of both the enormous cost management of samples and the degradation of the biological material.

Sequencing errors in Native American mitogenomes: impact on clade definitions, haplogroup assignation, and beyond

Available evidence allows the interpretation that some cases of absence of otherwise expected variation, based on phylogenetic expectations in mitogenomes of Native American origin, are due to artificial recombination rather than to homoplasy, while other more complex scenarios involving combination of original Cambridge Reference Sequence mistakes plus incomplete or incorrect scoring of variation are also showed. Several instances of mismatched control and coding regions as well as partially duplicated HV2 are observed in Peruvians, while intra-haplogroup chimaeras of different D1 subhaplogroups are referred to in Mexican Native Americans. A revised definition for haplogroup B2h is proposed, and preventive quality control measures are suggested.

Mitochondrial DNA insert into CD40 ligand gene-associated X-linked hyper-IgM syndrome

BACKGROUND

X-linked hyper-IgM (X-HIGM), which results from mutations in the CD40LG gene located on chromosome Xq26.3, is the most common form of HIGM. To date, more than 130 variants of the CD40L gene have been reported. We described a patient with novel de novo nuclear mitochondrial DNA sequences (NUMTs) in the CD40LG gene that have resulted in X-HIGM.

METHODS

Whole-exome sequencing (WES) analysis was used to screen for causal variants in the genome, and the candidate breakpoint was confirmed by Sanger sequencing.

RESULTS

A new mutation of CD40LG, which deletes A at position 17 followed by a 147-nucleotide from mitochondrial DNA copies insertion in exon 1, was detected in a 20-month-old boy harbouring an X-HIGM combined with immunodeficiency syndrome.

CONCLUSION

This is one of the few cases of a human genetic disease caused by nuclear mitochondrial DNA sequences (NUMTs). The presented data serve to demonstrate that de novo NUMT transfer of nucleic acid is a novel mechanism of X-HIGM.

The Prevalence and Impact of Model Violations in Phylogenetic Analysis

In phylogenetic inference, we commonly use models of substitution which assume that sequence evolution is stationary, reversible, and homogeneous (SRH). Although the use of such models is often criticized, the extent of SRH violations and their effects on phylogenetic inference of tree topologies and edge lengths are not well understood. Here, we introduce and apply the maximal matched-pairs tests of homogeneity to assess the scale and impact of SRH model violations on 3,572 partitions from 35 published phylogenetic data sets. We show that roughly one-quarter of all the partitions we analyzed (23.5%) reject the SRH assumptions, and that for 25% of data sets, tree topologies inferred from all partitions differ significantly from topologies inferred using the subset of partitions that do not reject the SRH assumptions. This proportion increases when comparing trees inferred using the subset of partitions that rejects the SRH assumptions, to those inferred from partitions that do not reject the SRH assumptions. These results suggest that the extent and effects of model violation in phylogenetics may be substantial. They highlight the importance of testing for model violations and possibly excluding partitions that violate models prior to tree reconstruction. Our results also suggest that further effort in developing models that do not require SRH assumptions could lead to large improvements in the accuracy of phylogenomic inference. The scripts necessary to perform the analysis are available in https://github.com/roblanf/SRHtests, and the new tests we describe are available as a new option in IQ-TREE (http://www.iqtree.org).

A large-scale species level dated angiosperm phylogeny for evolutionary and ecological analyses

Phylogenies are a central and indispensable tool for evolutionary and ecological research. Even though most angiosperm families are well investigated from a phylogenetic point of view, there are far less possibilities to carry out large-scale meta-analyses at order level or higher. Here, we reconstructed a large-scale dated phylogeny including nearly 1/8th of all angiosperm species, based on two plastid barcoding genes, matK (incl. trnK) and rbcL. Novel sequences were generated for several species, while the rest of the data were mined from GenBank. The resulting tree was dated using 56 angiosperm fossils as calibration points. The resulting megaphylogeny is one of the largest dated phylogenetic tree of angiosperms yet, consisting of 36,101 sampled species, representing 8,399 genera, 426 families and all orders. This novel framework will be useful for investigating different broad scale research questions in ecological and evolutionary biology.

Molecular markers used in forensic genetics

Forensic genetics is a field that has become subject to increasing interest in recent years. Both the technology and the markers used for forensic purposes have changed since the 1980s. The minisatellite sequences used in the famous Pitchfork case introduced genetics to the forensic sciences. Minisatellite sequences have now been replaced by more sensitive microsatellite markers, which have become the basis for the creation of genetic profile databases. Modern molecular methods also exploit single nucleotide polymorphisms, which are often the only way to identify degraded DNA samples. The same type of variation is taken into consideration in attempting to establish the ethnicity of a perpetrator and to determine phenotypic traits such as the eye or hair colour of the individual who is the source of the genetic material. This paper contains a review of the techniques and molecular markers used in human and animal forensic genetics, and also presents the potential trends in forensic genetics such as phenotyping.

Next generation database search algorithm for forensic mitogenome analyses

Mitochondrial DNA (mtDNA) variation is being reported relative to the corrected version of the first sequenced human mitochondrial genome. A review of the existing literature across disciplines that employ mtDNA demonstrates that insertions and deletions are not reported in a standardized way. This may lead to false exclusions of identical sequences, unidentified matches in missing persons mtDNA databases, biased mtDNA database frequency estimates and overestimation of the genetic evidence. Seven years ago we introduced alignment-free database search software (SAM) and implemented it into the mtDNA database EMPOP (https://empop.online) to produce reliable and conservative frequency estimates that are required in the forensic context. However, ambiguity remained in how laboratories have been reporting mitotypes, as often more than one single alignment of a given mtDNA sequence was feasible. In order to overcome this limitation we here describe a concept and provide software for producing stable, harmonized phylogenetic alignment of mtDNA sequences for database searches. The new software SAM 2 will be made available via EMPOP and provide the user with the already established conservative frequency estimates. In addition, SAM 2 offers the rCRS-coded haplotype of a given mtDNA sequence following the established and widely accepted phylogenetic alignment. This provides the user with feedback on how mitotypes are stored in EMPOP and how they should be reported in order to harmonize nomenclature. Finally, this approach does not only permit reliable mtDNA nomenclature in forensics but invites related disciplines to take advantage of a standardized way of reporting mtDNA variation, thus closing the ranks between different genetic fields and supporting dialogue and collaboration between mtDNA scholars from various disciplines.

Numtogenesis as a mechanism for development of cancer

Transfer of genetic material from cytoplasmic organelles to the nucleus, an ongoing process, has implications in evolution, aging, and human pathologies such as cancer. The transferred mitochondrial DNA (mtDNA) fragments in the nuclear genome are called nuclear mtDNA or NUMTs. We have named the process numtogenesis, defining the term as the transfer of mtDNA into the nuclear genome, or, less specifically, the transfer of mitochondria or mitochondrial components into the nucleus. There is increasing evidence of the involvement of NUMTs in human biology and pathology. Although information pertaining to NUMTs and numtogenesis is sparse, the role of this aspect of mitochondrial biology to human cancers is apparent. In this review, we present available knowledge about the origin and mechanisms of numtogenesis, with special emphasis on the role of NUMTs in human malignancies. We describe studies undertaken in our laboratory and in others and discuss the influence of NUMTs in tumor initiation and progression and in survival of cancer patients. We describe suppressors of numtogenesis and evolutionary conserved mechanisms underlying numtogenesis in cancer. An understanding the emerging field of numtogenesis should allow comprehension of this process in various malignancies and other diseases and, more generally, in human health.

Assessment of sex in a modern Turkish population using cranial anthropometric parameters

The utilization of radiological imaging methods in anthropometric studies is being expanded by the application of modern imaging methods, leading to a decrease in costs, a decrease in the time required for analysis and the ability to create three-dimensional images. This retrospective study investigated 400 patients within the 18-45-years age group (mean age: 30.7±11.2years) using cranial computed tomography images. We measured 14 anthropometric parameters (basion-bregma height, basion-prosthion length, maximum cranial length and cranial base lengths, maximum cranial breadth, bizygomatic diameter, upper facial breadth, bimastoid diameter, orbital breadth, orbital length, biorbital breadth, interorbital breadth, foramen magnum breadth and foramen magnum length) of cranial measurements. The intra- and inter-observer repeatability and consistency were good. From the results of logistic regression analysis using morphometric measurements, the most conspicuous measurements in terms of dimorphism were maximum cranial length, bizygomatic diameter, basion-bregma height, and cranial base length. The most dimorphic structure was the bizygomatic diameter with an accuracy rate of 83% in females and 77% in males. In this study, 87.5% of females and 87.0% of males were classified accurately by this model including four parameters with a sensitivity of 91.5% and specificity of 85.0%. In conclusion, CT cranial morphometric analysis may be reliable for the assessment of sex in the Turkish population and is recommended for comparison of data of modern populations with those of former populations. Additionally, cranial morphometric data that we obtained from modern Turkish population may reveal population specific data, which may help current criminal investigations and identification of disaster victims.

Massively parallel sequencing of the entire control region and targeted coding region SNPs of degraded mtDNA using a simplified library preparation method

The application of next-generation sequencing (NGS) to forensic genetics is being explored by an increasing number of laboratories because of the potential of high-throughput sequencing for recovering genetic information from multiple markers and multiple individuals in a single run. A cumbersome and technically challenging library construction process is required for NGS. In this study, we propose a simplified library preparation method for mitochondrial DNA (mtDNA) analysis that involves two rounds of PCR amplification. In the first-round of multiplex PCR, six fragments covering the entire mtDNA control region and 22 fragments covering interspersed single nucleotide polymorphisms (SNPs) in the coding region that can be used to determine global haplogroups and East Asian haplogroups were amplified using template-specific primers with read sequences. In the following step, indices and platform-specific sequences for the MiSeq(®) system (Illumina) were added by PCR. The barcoded library produced using this simplified workflow was successfully sequenced on the MiSeq system using the MiSeq Reagent Nano Kit v2. A total of 0.4 GB of sequences, 80.6% with base quality of >Q30, were obtained from 12 degraded DNA samples and mapped to the revised Cambridge Reference Sequence (rCRS). A relatively even read count was obtained for all amplicons, with an average coverage of 5200 × and a less than three-fold read count difference between amplicons per sample. Control region sequences were successfully determined, and all samples were assigned to the relevant haplogroups. In addition, enhanced discrimination was observed by adding coding region SNPs to the control region in in silico analysis. Because the developed multiplex PCR system amplifies small-sized amplicons (<250 bp), NGS analysis using the library preparation method described here allows mtDNA analysis using highly degraded DNA samples.

The Use and Effectiveness of Triple Multiplex System for Coding Region Single Nucleotide Polymorphism in Mitochondrial DNA Typing of Archaeologically Obtained Human Skeletons from Premodern Joseon Tombs of Korea

Previous study showed that East Asian mtDNA haplogroups, especially those of Koreans, could be successfully assigned by the coupled use of analyses on coding region SNP markers and control region mutation motifs. In this study, we tried to see if the same triple multiplex analysis for coding regions SNPs could be also applicable to ancient samples from East Asia as the complementation for sequence analysis of mtDNA control region. By the study on Joseon skeleton samples, we know that mtDNA haplogroup determined by coding region SNP markers successfully falls within the same haplogroup that sequence analysis on control region can assign. Considering that ancient samples in previous studies make no small number of errors in control region mtDNA sequencing, coding region SNP analysis can be used as good complimentary to the conventional haplogroup determination, especially of archaeological human bone samples buried underground over long periods.

Complete mitochondrial genome database and standardized classification system for Canis lupus familiaris

To contribute to the complete mitogenome database of the species Canis lupus familiaris and shed more light on its origin, we have sequenced mitochondrial genomes of 120 modern dogs from worldwide populations. Together with all the previously published mitogenome sequences of acceptable quality, we have reconstructed a global phylogenetic tree of 555 C. l. familiaris mitogenomes and standardized haplogroup nomenclature. The phylogenetic tree presented here and available online at http://clf.mtdna.tree.cm.umk.pl/ could be further used by forensic and evolutionary geneticists as well cynologists, for data quality control and unambiguous haplogroup classification. Our in-depth phylogeographic analysis of all C. l. familiaris mitogenomes confirmed that domestic dogs may have originated in East Asia during the Mesolithic and Upper Paleolithic time periods and started to expand to other parts of the world during Neolithic times.

Mitochondrial DNA mutations in single human blood cells

Determination mitochondrial DNA (mtDNA) sequences from extremely small amounts of DNA extracted from tissue of limited amounts and/or degraded samples is frequently employed in medical, forensic, and anthropologic studies. Polymerase chain reaction (PCR) amplification followed by DNA cloning is a routine method, especially to examine heteroplasmy of mtDNA mutations. In this review, we compare the mtDNA mutation patterns detected by three different sequencing strategies. Cloning and sequencing methods that are based on PCR amplification of DNA extracted from either single cells or pooled cells yield a high frequency of mutations, partly due to the artifacts introduced by PCR and/or the DNA cloning process. Direct sequencing of PCR product which has been amplified from DNA in individual cells is able to detect the low levels of mtDNA mutations present within a cell. We further summarize the findings in our recent studies that utilized this single cell method to assay mtDNA mutation patterns in different human blood cells. Our data show that many somatic mutations observed in the end-stage differentiated cells are found in hematopoietic stem cells (HSCs) and progenitors within the CD34(+) cell compartment. Accumulation of mtDNA variations in the individual CD34+ cells is affected by both aging and family genetic background. Granulocytes harbor higher numbers of mutations compared with the other cells, such as CD34(+) cells and lymphocytes. Serial assessment of mtDNA mutations in a population of single CD34(+) cells obtained from the same donor over time suggests stability of some somatic mutations. CD34(+) cell clones from a donor marked by specific mtDNA somatic mutations can be found in the recipient after transplantation. The significance of these findings is discussed in terms of the lineage tracing of HSCs, aging effect on accumulation of mtDNA mutations and the usage of mtDNA sequence in forensic identification.

Improved visibility of character conflicts in quasi-median networks with the EMPOP NETWORK software

Aim

To provide a valuable tool for graphical representation of mitochondrial DNA (mtDNA) data that enables visual emphasis on complex substructures within the network to highlight possible ambiguities and errors.

Method

We applied the new NETWORK graphical user interface, available via EMPOP (European DNA Profiling Group Mitochondrial DNA Population Database; www.empop.org) by means of two mtDNA data sets that were submitted for quality control.

Results

The quasi-median network torsi of the two data sets resulted in complex reticulations, suggesting ambiguous data. To check the corresponding raw data, accountable nodes and connecting branches of the network could be identified by highlighting induced subgraphs with concurrent dimming of their complements. This is achieved by accentuating the relevant substructures in the network: mouse clicking on a node displays a list of all mtDNA haplotypes included in that node; the selection of a branch specifies the mutation(s) connecting two nodes. It is indicated to evaluate these mutations by means of the raw data.

Conclusion

Inspection of the raw data confirmed the presence of phantom mutations due to suboptimal electrophoresis conditions and data misinterpretation. The network software proved to be a powerful tool to highlight problematic data and guide quality control of mtDNA data tables.

Haplotype diversity in mitochondrial DNA hypervariable region in a population of southeastern Brazil

Brazilian population derives from Native Amerindians, Europeans, and Africans. Southeastern Brazil is the most populous region of the country. The present study intended to characterize the maternal genetic ancestry of 290 individuals from southeastern (Brazil) population. Thus, we made the sequencing of the three hypervariable regions (HV1, HV2, and HV3) of the mitochondrial DNA (mtDNA). The statistical analyses were made using Arlequin software, and the median-joining haplotype networks were generated using Network software. The analysis of three hypervariable regios showed 230 (79.3 %) unique haplotypes and the most common haplotype was "263G" carried by 12 (4.1 %) individuals. The strikingly high variability generated by intense gene flow is mirrored in a high sequence diversity (0.9966 ± 0.0010), and the probability of two random individuals showing identical mtDNA haplotypes were 0.0068. The analysis of haplogroup distribution revealed that 36.9 % (n = 107) presented Amerindian haplogroups, 35.2 % (n = 102) presented African haplogroups, 27.6 % (n = 80) presented European haplogroups, and one (0.3 %) individual presented East Asian haplogroup, evidencing that the southeastern population is extremely heterogeneous and the coexistence of matrilineal lineages with three different phylogeographic origins. The genetic diversity found in the mtDNA control region in the southeastern Brazilian population reinforces the importance of increased national database in order to be important and informative in forensic cases.

Reviewing population studies for forensic purposes: Dog mitochondrial DNA

The identification of dog hair through mtDNA analysis has become increasingly important in the last 15 years, as it can provide associative evidence connecting victims and suspects. The evidential value of an mtDNA match between dog hair and its potential donor is determined by the random match probability of the haplotype. This probability is based on the haplotype’s population frequency estimate. Consequently, implementing a population study representative of the population relevant to the forensic case is vital to the correct evaluation of the evidence. This paper reviews numerous published dog mtDNA studies and shows that many of these studies vary widely in sampling strategies and data quality. Therefore, several features influencing the representativeness of a population sample are discussed. Moreover, recommendations are provided on how to set up a dog mtDNA population study and how to decide whether or not to include published data. This review emphasizes the need for improved dog mtDNA population data for forensic purposes, including targeting the entire mitochondrial genome. In particular, the creation of a publicly available database of qualitative dog mtDNA population studies would improve the genetic analysis of dog traces in forensic casework.

Assessing DNA barcoding as a tool for species identification and data quality control

In recent years, the number of sequences of diverse species submitted to GenBank has grown explosively and not infrequently the data contain errors. This problem is extensively recognized but not for invalid or incorrectly identified species, sample mixed-up, and contamination. DNA barcoding is a powerful tool for identifying and confirming species and one very important application involves forensics. In this study, we use DNA barcoding to detect erroneous sequences in GenBank by evaluating deep intraspecific and shallow interspecific divergences to discover possible taxonomic problems and other sources of error. We use the mitochondrial DNA gene encoding cytochrome b (Cytb) from turtles to test the utility of barcoding for pinpointing potential errors. This gene is widely used in phylogenetic studies of the speciose group. Intraspecific variation is usually less than 2.0% and in most cases it is less than 1.0%. In comparison, most species differ by more than 10.0% in our dataset. Overlapping intra- and interspecific percentages of variation mainly involve problematic identifications of species and outdated taxonomies. Further, we detect identical problems in Cytb from Insectivora and Chiroptera. Upon applying this strategy to 47,524 mammalian CoxI sequences, we resolve a suite of potentially problematic sequences. Our study reveals that erroneous sequences are not rare in GenBank and that the DNA barcoding can serve to confirm sequencing accuracy and discover problems such as misidentified species, inaccurate taxonomies, contamination, and potential errors in sequencing.

A one step multiplex PCR assay for rapid screening of East Asian mtDNA haplogroups on forensic samples

The mitochondrial DNA (mtDNA) haplogroup typing has become an essential tool to study human evolutionary history and to infer the matrilineal bio-geographic ancestry. In forensic field, the screening of mtDNA haplogroups by genotyping of mtDNA single nucleotide polymorphisms (SNPs) can help guarantee the quality of mtDNA sequence data as well as can reduce the need to sequence samples that do not match. Here, a multiplex mutagenically separated (MS) polymerase chain reaction (PCR) system was developed for simultaneous rapid detection of 14 coding region SNPs and one deletion motif representing common mtDNA haplogroups of East Asia. The multiplex MS PCR system we developed has the advantage of being a one step procedure that requires only a single PCR amplification with allele-specific primers and allowing straightforward designation of haplogroups along the branches of the phylogenetic tree. Therefore, it would be a simple, rapid, and reliable detection method useful for large-scale screening of mtDNA variations to determine East Asian mtDNA haplogroups.

Capillary electrophoresis of Big-Dye terminator sequencing reactions for human mtDNA Control Region haplotyping in the identification of human remains

Cycle sequencing reaction with Big-Dye terminators provides the methodology to analyze mtDNA Control Region amplicons by means of capillary electrophoresis. DNA sequencing with ddNTPs or terminators was developed by (1). The progressive automation of the method by combining the use of fluorescent-dye terminators with cycle sequencing has made it possible to increase the sensibility and efficiency of the method and hence has allowed its introduction into the forensic field. PCR-generated mitochondrial DNA products are the templates for sequencing reactions. Different set of primers can be used to generate amplicons with different sizes according to the quality and quantity of the DNA extract providing sequence data for different ranges inside the Control Region.

Mitochondria in anthropology and forensic medicine

Mitochondria's role in crucial metabolic pathways is probably the first answer which comes to our minds for the question: what do these tiny organelles serve for? However, specific features of their DNA made them extremely useful also in the field of anthropology and forensics. MtDNA analyses became a milestone in the complex task of unraveling earliest human migrations. Evidence provided by these experiments left no doubts on modern humans origins pointing to Africa being our cradle. It also contributed to interpretation of putative ways of our dispersal around Asia and Americas thousands years ago. On the other hand, analysis of mtDNA is well established and valuable tool in forensic genetics. When other definitely more popular markers give no answer on identity, it is the time to employ information carried by mitochondria. This chapter summarizes not only current reports on the role of mitochondria in forensics and reconstruction of modern humans phylogeny, but also calls one's attention to a broad range of difficulties and constraints associated with mtDNA analyses.

Current genetic methodologies in the identification of disaster victims and in forensic analysis

This review presents the basic problems and currently available molecular techniques used for genetic profiling in disaster victim identification (DVI). The environmental conditions of a mass disaster often result in severe fragmentation, decomposition and intermixing of the remains of victims. In such cases, traditional identification based on the anthropological and physical characteristics of the victims is frequently inconclusive. This is the reason why DNA profiling became the gold standard for victim identification in mass-casualty incidents (MCIs) or any forensic cases where human remains are highly fragmented and/or degraded beyond recognition. The review provides general information about the sources of genetic material for DNA profiling, the genetic markers routinely used during genetic profiling (STR markers, mtDNA and single-nucleotide polymorphisms [SNP]) and the basic statistical approaches used in DNA-based disaster victim identification. Automated technological platforms that allow the simultaneous analysis of a multitude of genetic markers used in genetic identification (oligonucleotide microarray techniques and next-generation sequencing) are also presented. Forensic and population databases containing information on human variability, routinely used for statistical analyses, are discussed. The final part of this review is focused on recent developments, which offer particularly promising tools for forensic applications (mRNA analysis, transcriptome variation in individuals/populations and genetic profiling of specific cells separated from mixtures).

Determination of sex from juvenile crania by means of discriminant function analysis

This study provides evidence of craniofacial growth variation between the sexes in juveniles of European descent. Data were collected from lateral cephalometric radiographs belonging to the Michigan Craniofacial Growth Study. The collection consists of longitudinal lateral radiographs that represent individuals 5-16 years of age. Each radiograph was manually traced on hyprint vellum from which eight craniometric points were identified. From these points, 20 craniofacial measurements were recorded and then analyzed by means of a canonical discriminant function analysis. Sex classification equations were then created by applying a backward stepwise procedure to the discriminant functions. The analysis demonstrates the presence of sexually dimorphic differences in craniofacial growth. The neurocranium is the most sexually dimorphic region of the juvenile craniofacial skeleton, until the onset of puberty. Size is the main source of variation with males having taller and longer heads than females. Overall, sex classification in the sample ranges from 78 to 89% accuracy.

Uncovering the profile of somatic mtDNA mutations in Chinese colorectal cancer patients

In the past decade, a high incidence of somatic mitochondrial DNA (mtDNA) mutations has been observed, mostly based on a fraction of the molecule, in various cancerous tissues; nevertheless, some of them were queried due to problems in data quality. Obviously, without a comprehensive understanding of mtDNA mutational profile in the cancerous tissue of a specific patient, it is unlikely to disclose the genuine relationship between somatic mtDNA mutations and tumorigenesis. To achieve this objective, the most straightforward way is to directly compare the whole mtDNA genome variation among three tissues (namely, cancerous tissue, para-cancerous tissue, and distant normal tissue) from the same patient. Considering the fact that most of the previous studies on the role of mtDNA in colorectal tumor focused merely on the D-loop or partial segment of the molecule, in the current study we have collected three tissues (cancerous, para-cancerous and normal tissues) respectively recruited from 20 patients with colorectal tumor and completely sequenced the mitochondrial genome of each tissue. Our results reveal a relatively lower incidence of somatic mutations in these patients; intriguingly, all somatic mutations are in heteroplasmic status. Surprisingly, the observed somatic mutations are not restricted to cancer tissues, for the para-cancer tissues and distant normal tissues also harbor somatic mtDNA mutations with a lower frequency than cancerous tissues but higher than that observed in the general population. Our results suggest that somatic mtDNA mutations in cancerous tissues could not be simply explained as a consequence of tumorigenesis; meanwhile, the somatic mtDNA mutations in normal tissues might reflect an altered physiological environment in cancer patients.

The role of the mitochondrial genome in ageing and carcinogenesis

Mitochondrial DNA mutations and polymorphisms have been the focus of intensive investigations for well over a decade in an attempt to understand how they affect fundamental processes such as cancer and aging. Initial interest in mutations occurring in mitochondrial DNA of cancer cells diminished when most were found to be the same mutations which occurred during the evolution of human mitochondrial haplogroups. However, increasingly correlations are being found between various mitochondrial haplogroups and susceptibility to cancer or diseases in some cases and successful aging in others.

MitoTool: a web server for the analysis and retrieval of human mitochondrial DNA sequence variations

MitoTool, a web-based bioinformatics platform, is designed for deciphering human mitochondrial DNA (mtDNA) data in batch mode. The platform has advantages in (i) parsing diverse types of mtDNA data; (ii) automatically classifying haplogroup according to mtDNA sequences or variants; (iii) discovering possibly missing variants of the samples with claimed haplogroups status; (iv) estimating the evolutionary conservation index, protein coding effect and potential pathogenicity of certain substitutions; (v) performing statistical analysis for haplogroup distribution frequency between case and control groups. Furthermore, it offers an integrated database for retrieving five types of mitochondrion-related information. The MitoTool is freely accessed at http://www.mitotool.org.

A reduced number of mtSNPs saturates mitochondrial DNA haplotype diversity of worldwide population groups

Background

The high levels of variation characterising the mitochondrial DNA (mtDNA) molecule are due ultimately to its high average mutation rate; moreover, mtDNA variation is deeply structured in different populations and ethnic groups. There is growing interest in selecting a reduced number of mtDNA single nucleotide polymorphisms (mtSNPs) that account for the maximum level of discrimination power in a given population. Applications of the selected mtSNP panel range from anthropologic and medical studies to forensic genetic casework.

Methodology/Principal Findings

This study proposes a new simulation-based method that explores the ability of different mtSNP panels to yield the maximum levels of discrimination power. The method explores subsets of mtSNPs of different sizes randomly chosen from a preselected panel of mtSNPs based on frequency. More than 2,000 complete genomes representing three main continental human population groups (Africa, Europe, and Asia) and two admixed populations (“African-Americans” and “Hispanics”) were collected from GenBank and the literature, and were used as training sets. Haplotype diversity was measured for each combination of mtSNP and compared with existing mtSNP panels available in the literature. The data indicates that only a reduced number of mtSNPs ranging from six to 22 are needed to account for 95% of the maximum haplotype diversity of a given population sample. However, only a small proportion of the best mtSNPs are shared between populations, indicating that there is not a perfect set of “universal” mtSNPs suitable for all population contexts. The discrimination power provided by these mtSNPs is much higher than the power of the mtSNP panels proposed in the literature to date. Some mtSNP combinations also yield high diversity values in admixed populations.

Conclusions/Significance

The proposed computational approach for exploring combinations of mtSNPs that optimise the discrimination power of a given set of mtSNPs is more efficient than previous empirical approaches. In contrast to precedent findings, the results seem to indicate that only few mtSNPs are needed to reach high levels of discrimination power in a population, independently of its ancestral background.

eCOMPAGT integrates mtDNA: import, validation and export of mitochondrial DNA profiles for population genetics, tumour dynamics and genotype-phenotype association studies

Background

Mitochondrial DNA (mtDNA) is widely being used for population genetics, forensic DNA fingerprinting and clinical disease association studies. The recent past has uncovered severe problems with mtDNA genotyping, not only due to the genotyping method itself, but mainly to the post-lab transcription, storage and report of mtDNA genotypes.

Description

eCOMPAGT, a system to store, administer and connect phenotype data to all kinds of genotype data is now enhanced by the possibility of storing mtDNA profiles and allowing their validation, linking to phenotypes and export as numerous formats. mtDNA profiles can be imported from different sequence evaluation programs, compared between evaluations and their haplogroup affiliations stored. Furthermore, eCOMPAGT has been improved in its sophisticated transparency (support of MySQL and Oracle), security aspects (by using database technology) and the option to import, manage and store genotypes derived from various genotyping methods (SNPlex, TaqMan, and STRs). It is a software solution designed for project management, laboratory work and the evaluation process all-in-one.

Conclusions

The extended mtDNA version of eCOMPAGT was designed to enable error-free post-laboratory data handling of human mtDNA profiles. This software is suited for small to medium-sized human genetic, forensic and clinical genetic laboratories. The direct support of MySQL and the improved database security options render eCOMPAGT a powerful tool to build an automated workflow architecture for several genotyping methods. eCOMPAGT is freely available at http://dbis-informatik.uibk.ac.at/ecompagt.

Human mitochondrial DNA complete amplification and sequencing: a new validated primer set that prevents nuclear DNA sequences of mitochondrial origin co-amplification

To date, there are no published primers to amplify the entire mitochondrial DNA (mtDNA) that completely prevent the amplification of nuclear DNA (nDNA) sequences of mitochondrial origin. The main goal of this work was to design, validate and describe a set of primers, to specifically amplify and sequence the complete human mtDNA, allowing the correct interpretation of mtDNA heteroplasmy in healthy and pathological samples. Validation was performed using two different approaches: (i) Basic Local Alignment Search Tool and (ii) amplification using isolated nDNA obtained from sperm cells by differential lyses. During the validation process, two mtDNA regions, with high similarity with nDNA, represent the major problematic areas for primer design. One of these could represent a non-published nuclear DNA sequence of mitochondrial origin. For two of the initially designed fragments, the amplification results reveal PCR artifacts that can be attributed to the poor quality of the DNA. After the validation, nine overlapping primer pairs to perform mtDNA amplification and 22 additional internal primers for mtDNA sequencing were obtained. These primers could be a useful tool in future projects that deal with mtDNA complete sequencing and heteroplasmy detection, since they represent a set of primers that have been tested for the non-amplification of nDNA.

Contamination and sample mix-up can best explain some patterns of mtDNA instabilities in buccal cells and oral squamous cell carcinoma

The study of somatic DNA instabilities constitutes a debatable topic because different causes can lead to seeming DNA alteration patterns between different cells or tissues from the same individual. Carcinogenesis or the action of a particular toxic could generate such patterns, and this is in fact the leitmotif of a number of studies on mitochondrial DNA (mtDNA) instability. Patterns of seeming instabilities could also arise from technical errors at any stage of the analysis (DNA extraction, amplification, mutation screening/sequencing, and documentation). Specifically, inadvertent DNA contamination or sample mixing would yield mosaic variation that could be erroneously interpreted as real mutation differences (instabilities) between tissues from the same individual. From the very beginning, mtDNA studies comparing cancerous to non-cancerous tissues have suffered from such mosaic results. We demonstrate here that the phylogenetic linkage of whole arrays of mtDNA mutations provides strong evidence of artificial recombination in previous studies on buccal cells and oral squamous cell carcinoma.

Median network analysis of defectively sequenced entire mitochondrial genomes from early and contemporary disease studies

Sequence analysis of the mitochondrial genome has become a routine method in the study of mitochondrial diseases. Quite often, the sequencing efforts in the search of pathogenic or disease-associated mutations are affected by technical and interpretive problems, caused by sample mix-up, contamination, biochemical problems, incomplete sequencing, misdocumentation and insufficient reference to previously published data. To assess data quality in case studies of mitochondrial diseases, it is recommended to compare any mtDNA sequence under consideration to their phylogenetically closest lineages available in the Web. The median network method has proven useful for visualizing potential problems with the data. We contrast some early reports of complete mtDNA sequences to more recent total mtDNA sequencing efforts in studies of various mitochondrial diseases. We conclude that the quality of complete mtDNA sequences generated in the medical field in the past few years is somewhat unsatisfactory and may even fall behind that of pioneer manual sequencing in the early nineties. Our study provides a paradigm for an a posteriori evaluation of sequence quality and for detection of potential problems with inferring a pathogenic status of a particular mutation.

Analysis of mitochondrial DNA polymorphisms in Guangdong Han Chinese

Previous investigations on Chinese mitochondrial DNA (mtDNA) variation revealed that the matrilineal gene pool of southern Han Chinese is rather complex, with much higher genetic diversity and more basal/ancient lineages than the northern Hans. The extreme case is Guangdong Han populations, among which pronounced (matrilineal) differentiation has been observed, indicative of complex demography of the region. To get more insights into the maternal makeup of southern Han Chinese, mtDNA variation of a total of 106 individuals sampled from Dongguan, Guangdong Province, China, was analyzed in this study. With the aid of the information from control-region hypervariable segments I and II (HVS-I and -II) as well as some necessary coding-region segments, the phylogenetic status of all mtDNAs under examination were determined according to the reconstructed East Asian mtDNA tree. In this way, the mtDNAs have been classified into various haplogroups or sub-haplogroups. The southern-prevalent haplogroups, such as R9 (20.8%), B (17.9%), M7b (14.2%), show relatively high distribution frequencies in Dongguan Hans; whereas the frequencies of Northern-prevalent haplogroups (with the exception of D) are quite low: C (1.9%), G2 (1.9%) and Z (1.9%), indicating the southern-origin of Dongguan Hans.

Genetic signatures of parental contribution in black and white populations in Brazil

Two hundred and three individuals classified as white were tested for 11 single nucleotide polymorphisms plus two insertion/deletions in their Y-chromosomes. A subset of these individuals (n = 172) was also screened for sequences in the first hypervariable segment of their mitochondrial DNA (mtDNA). In addition, complementary studies were done for 11 of the 13 markers indicated above in 54 of 107 black subjects previously investigated in this southern Brazilian population. The prevalence of Y-chromosome haplogroups among whites was similar to that found in the Azores (Portugal) or Spain, but not to that of other European countries. About half of the European or African mtDNA haplogroups of these individuals were related to their places of origin, but not their Amerindian counterparts. Persons classified in these two categories of skin color and related morphological traits showed distinct genomic ancestries through the country. These findings emphasize the need to consider in Brazil, despite some general trends, a notable heterogeneity in the pattern of admixture dynamics within and between populations/groups.

A homogenous nature of native Chinese duck matrilineal pool

Background

China, with around 30 unique breeds, has a diverse duck genetic pool. Currently, there is no systematic report which investigates the genetic diversity, phylogenetic relationship, and matrilineal genetic structure of these domestic breeds and wild mallards (Anas platyrhynchos).

Results

In this study, we sequenced the mitochondrial DNA (mtDNA) control region segments in 278 domestic ducks (Anas platyrhynchos domestica) from 19 indigenous breeds/populations and 70 wild mallard samples and analyzed them together with the 101 control region sequences from published sources. Fifty-two samples were then sequenced for a cytochrome b (Cyt b) gene fragment to solidify the pattern emerged from the control region sequences. All domestic duck and wild mallard haplotypes were essentially indistinguishable and were clustered together in the phylogenetic tree. There was no geographic differentiation and breed/population-specific distribution of duck lineages.

Conclusion

Our results showed that unlike other domesticated farm animals in China such as chicken, cattle, goat, and yak with multiple matrilineal components, the matrilineal pool of Chinese ducks was homogenous.

Distilling artificial recombinants from large sets of complete mtDNA genomes

Background

Large-scale genome sequencing poses enormous problems to the logistics of laboratory work and data handling. When numerous fragments of different genomes are PCR amplified and sequenced in a laboratory, there is a high immanent risk of sample confusion. For genetic markers, such as mitochondrial DNA (mtDNA), which are free of natural recombination, single instances of sample mix-up involving different branches of the mtDNA phylogeny would give rise to reticulate patterns and should therefore be detectable.

Methodology/Principal Findings

We have developed a strategy for comparing new complete mtDNA genomes, one by one, to a current skeleton of the worldwide mtDNA phylogeny. The mutations distinguishing the reference sequence from a putative recombinant sequence can then be allocated to two or more different branches of this phylogenetic skeleton. Thus, one would search for two (or three) near-matches in the total mtDNA database that together best explain the variation seen in the recombinants. The evolutionary pathway from the mtDNA tree connecting this pair together with the recombinant then generate a grid-like median network, from which one can read off the exchanged segments.

Conclusions

We have applied this procedure to a large collection of complete human mtDNA sequences, where several recombinants could be distilled by our method. All these recombinant sequences were subsequently corrected by de novo experiments – fully concordant with the predictions from our data-analytical approach.

Mitochondrial DNA control region polymorphism in the population of Alagoas state, north-eastern Brazil

The sequences of the two hypervariable (HV) segments of the mitochondrial DNA (mtDNA) control region were determined in 167 randomly selected, unrelated individuals living in the state of Alagoas, north-eastern Brazil. One hundred and forty-five different haplotypes, associated with 139 variable positions, were determined. More than 95% of the mtDNA sequences could be allocated to specific mtDNA haplogroups according to the mutational motifs. Length heteroplasmy in the C-stretch HV1 and HV2 regions was observed in 22 and 11%, respectively, of the population sample. The genetic diversity was estimated to be 0.9975 and the probability of two random individuals presenting identical mtDNA haplotypes was 0.0084. The most frequent haplotype was shared by six individuals. All sequences showed high-quality values and phantom mutations were not detected. The diversity revealed in the mitochondrial control region indicates the importance of this locus for forensic casework and population studies within Alagoas, Brazil.

Application of DNA-based methods in forensic entomology

A forensic entomological investigation can benefit from a variety of widely practiced molecular genotyping methods. The most commonly used is DNA-based specimen identification. Other applications include the identification of insect gut contents and the characterization of the population genetic structure of a forensically important insect species. The proper application of these procedures demands that the analyst be technically expert. However, one must also be aware of the extensive list of standards and expectations that many legal systems have developed for forensic DNA analysis. We summarize the DNA techniques that are currently used in, or have been proposed for, forensic entomology and review established genetic analyses from other scientific fields that address questions similar to those in forensic entomology. We describe how accepted standards for forensic DNA practice and method validation are likely to apply to insect evidence used in a death or other forensic entomological investigation.

Somatic mutations of mitochondrial genome in early stage breast cancer

The complete mitochondrial genomes of the primary cancerous, matched paracancerous normal and distant normal tissues from 10 early-stage breast cancer patients were analyzed in this study, with special attempt (i) to investigate whether the reported high frequency of mitochondrial DNA (mtDNA) somatic mutations in breast cancer could be repeated under a stringent data quality control, and (ii) to characterize the spectrum of mtDNA somatic mutations in Chinese breast cancer patients and evaluate their potential significance in early cancer diagnosis. Two heteroplasmic somatic transitions (T2275C and A8601G) were identified in our samples. The transition A8601G was present in the primary cancerous and paracancerous normal tissues from patient no. 3. Transition T2275C was found in the primary cancerous tissue but not in other normal tissues from patient no. 6; this transition has been reported in the colonic crypts and is located at a highly conserved site in the 16S rRNA gene. Subsequent cloning sequencing confirmed the absence of both mutations in the distant normal tissues from the 2 patients. The overall rate of somatic mutations in our patients was much lower than those of previous studies of breast cancer. Our results gave support to the recent claim that the high frequency of mtDNA somatic mutations in cancer studies is overestimated. Based on the mtDNA mutation pattern in early stage breast cancer observed in this study, we cautioned the enthusiasm and efforts to look for somatic mutations that were of diagnostic value in cancer early detection.