New Insights Into Mitochondrial DNA Reconstruction and Variant Detection in Ancient Samples

An easy-to-use pipeline to analyze amplicon-based Next Generation Sequencing results of human mitochondrial DNA from degraded samples

Genome and transcriptome examinations have become more common due to Next-Generation Sequencing (NGS), which significantly increases throughput and depth coverage while reducing costs and time. Mitochondrial DNA (mtDNA) is often the marker of choice in degraded samples from archaeological and forensic contexts, as its higher number of copies can improve the success of the experiment. Among other sequencing strategies, amplicon-based NGS techniques are currently being used to obtain enough data to be analyzed. There are some pipelines designed for the analysis of ancient mtDNA samples and others for the analysis of amplicon data. However, these pipelines pose a challenge for non-expert users and cannot often address both ancient and forensic DNA particularities and amplicon-based sequencing simultaneously. To overcome these challenges, a user-friendly bioinformatic tool was developed to analyze the non-coding region of human mtDNA from degraded samples recovered in archaeological and forensic contexts. The tool can be easily modified to fit the specifications of other amplicon-based NGS experiments. A comparative analysis between two tools, MarkDuplicates from Picard and dedup parameter from fastp, both designed for duplicate removal was conducted. Additionally, various thresholds of PMDtools, a specialized tool designed for extracting reads affected by post-mortem damage, were used. Finally, the depth coverage of each amplicon was correlated with its level of damage. The results obtained indicated that, for removing duplicates, dedup is a better tool since retains more non-repeated reads, that are removed by MarkDuplicates. On the other hand, a PMDS = 1 in PMDtools was the threshold that allowed better differentiation between present-day and ancient samples, in terms of damage, without losing too many reads in the process. These two bioinformatic tools were added to a pipeline designed to obtain both haplotype and haplogroup of mtDNA. Furthermore, the pipeline presented in the present study generates information about the quality and possible contamination of the sample. This pipeline is designed to automatize mtDNA analysis, however, particularly for ancient samples, some manual analyses may be required to fully validate results since the amplicons that used to be more easily recovered were the ones that had fewer reads with damage, indicating that special care must be taken for poor recovered samples.

Biomolecular analysis of the Epigravettian human remains from Riparo Tagliente in northern Italy

The Epigravettian human remains from Riparo Tagliente in northern Italy represent some of the earliest evidence of human occupation in the southern Alpine slopes after the Last Glacial Maximum. Genomic analyses of the 17,000-year-old Tagliente 2 mandible revealed the oldest presence of a genetic profile with affinities to the Near East in the Italian peninsula, which later became the most widespread hunter-gatherer ancestry across Europe. However, a comparable biomolecular characterization of the Tagliente 1 burial remains unavailable, preventing us from defining its biological relationships with Tagliente 2. Here, we apply paleogenomic, isotopic, and radiocarbon dating analyses on a femur fragment of Tagliente 1 and compare the reconstructed data with previously reported results from Tagliente 2. Despite their different isotopic signatures and non-overlapping radiocarbon dates, we reveal that the two human remains belong to the same male individual. We determine that the distinct isotopic values can be explained by different dietary practices during lifetime, whereas the non-overlapping radiocarbon dates can be caused by minimal radiocarbon contamination, possibly deriving from chemical treatments for conservation purposes. These findings highlight the importance of interdisciplinary biomolecular studies in offering new perspectives on the Palaeolithic fossil record and addressing long-standing bioarchaeological questions.

Revisiting mitogenome evolution in Medusozoa with eight new mitochondrial genomes

Mitogenomics has improved our understanding of medusozoan phylogeny. However, sequenced medusozoan mitogenomes remain scarce, and Medusozoa phylogeny studies often analyze mitogenomic sequences without incorporating mitogenome rearrangements. To better understand medusozoan evolution, we analyzed Medusozoa mitogenome phylogeny by sequencing and assembling eight mitogenomes from three classes (Cubozoa, Hydrozoa, and Scyphozoa). We reconstructed the mitogenome phylogeny using these mitogenomes and 84 other existing cnidarian mitogenomes to study mitochondrial gene rearrangements. All reconstructed mitogenomes had 13 mitochondrial protein-coding genes and two ribosomal genes typical for Medusozoa. Non-cubozoan mitogenomes were all linear and had typical gene orders, while arrangement of genes in the fragmented Cubozoa (Morbakka sp.) mitogenome differed from other Cubozoa mitogenomes. Gene order comparisons and ancestral state reconstruction suggest minimal rearrangements within medusozoan classes except for Hydrozoa. Our findings support a staurozoan ancestral medusozoan gene order, expand the pool of available medusozoan mitogenomes, and enhance our understanding of medusozoan phylogenetic relationships.

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

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

Mitochondrial DNA Copy Number and Heteroplasmy in Monozygotic Twins Discordant for Schizophrenia

Schizophrenia (SZ) is a severe, complex, and common mental disorder with high heritability (80%), an adult age of onset, and high discordance (∼50%) in monozygotic twins (MZ). Extensive studies on familial and non-familial cases have implicated a number of segregating mutations and de novo changes in SZ that may include changes to the mitochondrial genome. Yet, no single universally causal variant has been identified, highlighting its extensive genetic heterogeneity. This report specifically focuses on the assessment of changes in the mitochondrial genome in a unique set of monozygotic twins discordant (MZD) for SZ using blood. Genomic DNA from six pairs of MZD twins and two sets of parents (N = 16) was hybridized to the Affymetrix Human SNP Array 6.0 to assess mitochondrial DNA copy number (mtDNA-CN). Whole genome sequencing (WGS) and quantitative polymerase chain reaction (qPCR) was performed for a subset of MZD pairs and their parents and was also used to derive mtDNA-CN estimates. The WGS data were further analyzed to generate heteroplasmy (HP) estimates. Our results show that mtDNA-CN estimates for within-pair and mother-child differences were smaller than comparisons involving unrelated individuals, as expected. MZD twins showed discordance in mtDNA-CN estimates and displayed concordance in directionality of differences for mtDNA-CN across all technologies. Further, qPCR performed better than Affymetrix in estimating mtDNA-CN based on relatedness. No reliable differences in HP were detected between MZD twins. The within-MZD differences in mtDNA-CN observed represent postzygotic somatic changes that may contribute to discordance of MZ twins for diseases, including SZ.

Using loop-primer mediated PCR to enhance the detection of poorly preserved DNA

Ancient DNA is vitally important in evolutionary research, and obtaining authentic ancient DNA sequences is critical for a proper analysis. However, it is difficult to acquire amplicons accurately and efficiently from ancient DNA templates using current techniques. Here, we established a loop-primer-mediated amplification method (L-PCR) to obtain target ancient DNA sequences with high accuracy and efficiency. The method was tested using 66 ancient samples (including 27 pig bones or teeth and 39 chicken bones) and serially diluted modern animal DNA templates. Compared to nested PCR, L-PCR was proven to be more efficient and accurate and could obtain more amplicons from both ancient pig samples and chicken bones and detect as low as 10−3 ng/μl modern pig template DNA. The efficiency was at least 100-fold that of the nested PCR. The results suggest that L-PCR is advantageous for obtaining authentic DNA sequences from poorly preserved or recalcitrant ancient specimens.

Genetic structure and differentiation from early bronze age in the mediterranean island of sicily: Insights from ancient mitochondrial genomes

Sicily is one of the main islands of the Mediterranean Sea, and it is characterized by a variety of archaeological records, material culture and traditions, reflecting the history of migrations and populations' interaction since its first colonization, during the Paleolithic. These deep and complex demographic and cultural dynamics should have affected the genomic landscape of Sicily at different levels; however, the relative impact of these migrations on the genomic structure and differentiation within the island remains largely unknown. The available Sicilian modern genetic data gave a picture of the current genetic structure, but the paucity of ancient data did not allow so far to make predictions about the level of historical variation. In this work, we sequenced and analyzed the complete mitochondrial genomes of 36 individuals from five different locations in Sicily, spanning from Early Bronze Age to Iron Age, and with different cultural backgrounds. The comparison with coeval groups from the Mediterranean Basin highlighted structured genetic variation in Sicily since Early Bronze Age, thus supporting a demic impact of the cultural transitions within the Island. Explicit model testing through Approximate Bayesian Computation allowed us to make predictions about the origin of Sicanians, one of the three indigenous peoples of Sicily, whose foreign origin from Spain, historically attributed, was not confirmed by our analysis of genetic data. Sicilian modern mitochondrial data show a different, more homogeneous, genetic composition, calling for a recent genetic replacement in the Island of pre-Iron Age populations, that should be further investigated.