Maternal ancestry and population history from whole mitochondrial genomes

Human mitochondrial DNA in public metagenomes: Opportunity or privacy threat?

Human DNA is unavoidably present in metagenomic analyses of human microbiomes. While current protocols remove human DNA before submission to public repositories, mitochondrial DNA (mtDNA) has been overlooked and frequently persists. We discuss the privacy risks and research opportunities associated with mtDNA, urging consideration by the scientific, ethics, and legal communities.

Picuris Pueblo oral history and genomics reveal continuity in US Southwest

Indigenous groups often encounter significant challenges when asserting ancestral claims and cultural affiliations based on oral histories, particularly in the USA where such narratives have historically been undervalued. Although ancient DNA offers a tool to complement traditional knowledge and address gaps in oral history, longstanding disregard for Indigenous sovereignty and beliefs has understandably led many Indigenous communities to distrust DNA studies1-4. Earlier research often focused on repatriation claims5-7, whereas more recent work has increasingly moved towards enhancing Tribal histories8,9. Here we present a collaborative study initiated by a federally recognized Native American tribe, the sovereign nation of Picuris Pueblo in the Northern Rio Grande region of New Mexico, USA, to address gaps in traditional knowledge and further their understanding of their population history and ancestry. We generated genomes from 16 ancient Picuris individuals and 13 present-day members of Picuris Pueblo, providing genomic data spanning the last millennium. We show genetic continuity between ancient and present-day Picuris, and more broadly with Ancestral Puebloans from Pueblo Bonito in Chaco Canyon10, 275 km to the west. This suggests a firm spatiotemporal link among these Puebloan populations of the North American Southwest. Furthermore, we see no evidence of population decline before European arrival11-13, and no Athabascan ancestry in individuals predating 1500 CE, challenging earlier migration hypotheses14-16. This work prioritizes Indigenous control of genetic data and brings together oral tradition, archaeology, ethnography and genetics.

Haplogroup Structure and Genetic Variation Analyses of Mitochondrial Genome SNPs in the Iranian Population

BACKGROUND

Mitochondrial DNA (mtDNA) is a valuable marker for population studies and forensic investigations. Recent advancements in massively parallel sequencing technologies enable whole mitochondrial genome sequencing. This study collected blood samples from unrelated Iranian participants from four ethnic groups: Persian, Kurd, Lur, and Azeri. We mapped mtDNA haplogroups according to genetic ancestry and investigated the ethnic similarities within the Iranian population.

METHODS

Complete mtDNA sequences were generated with targeted mtDNA sequencing method and haplogroups were determined on the base of mitogenome polymorphisms. Additionally, we used data from the whole exome sequencing (WES) of the current samples to compare the variants identified by two different mitochondrial testing methods. Principal component analysis (PCA) calculations were performed using the R software to determine diversity between unrelated individuals of various ethnicities.

RESULTS

A total of 129 sub-haplogroups were identified in 15 main haplogroups. The findings revealed high frequencies of haplogroups U and H (22.4% and 20.3%, respectively) in the Iranian population. The PCA scatter plots revealed overlapping diversity, with no distinct trends separating the groups in these four groups within the Iranian population. In the present samples, the WES method identified only 57.8% of the variants detected by the targeted mtDNA sequencing method.

CONCLUSION

Variant studies do not show much difference, which indicate a small genetic difference between the central ethnic groups of Iran. Furthermore, comparing the targeted whole mitochondrial genome to mitochondrial data from WES in our study samples highlights the notion that targeted entire mitochondrial genome is a gold standard method for variant detection.

Forensic DNA phenotyping: a review on SNP panels, genotyping techniques, and prediction models

In the past few years, forensic DNA phenotyping has attracted a strong interest in the forensic research. Among the increasing publications, many have focused on testing the available panels to infer biogeographical ancestry on less represented populations and understanding the genetic mechanisms underlying externally visible characteristics. However, there are currently no publications that gather all the existing panels limited to forensic DNA phenotyping and discuss the main technical limitations of the technique. In this review, we performed a bibliographic search in Scopus database of phenotyping-related literature, which resulted in a total of 48, 43, and 15 panels for biogeographical ancestry, externally visible characteristics, and both traits inference, respectively. Here we provide a list of commercial and non-commercial panels and the limitations regarding the lack of harmonization in terms of terminology (i.e., categorization and measurement of traits) and reporting, the lack of genetic knowledge and environment influence to select markers and develop panels, and the debate surrounding the selection of genotyping technologies and prediction models and algorithms. In conclusion, this review aims to be an updated guide and to present an overview of the current related literature.

Genealogical obscurement: mitochondrial replacement techniques and genealogical research

Mitochondrial replacement techniques (MRTs) are a new group of biotechnologies that aim to aid women whose eggs have disease-causing deleteriously mutated mitochondria to have genetically related healthy children. These techniques have also been used to aid women with poor oocyte quality and poor embryonic development, to have genetically related children. Remarkably, MRTs create humans with DNA from three sources: nuclear DNA from the intending mother and father, and mitochondrial DNA from the egg donor. In a recent publication Françoise Baylis argued that MRTs are detrimental for genealogical research via mitochondrial DNA because they would obscure the lines of individual descent. In this paper, I argue that MRTs do not obscure genealogical research, but rather that MRT-conceived children can have two mitochondrial lineages. I argue for this position by showing that MRTs are reproductive in nature and, thus, they create genealogy.

The evolution of contemporary livestock species: Insights from mitochondrial genome

The domestication of animals marks a pivotal moment in human history, profoundly influencing our demographic and cultural progress. This process has led to significant genetic, behavioral, and physical changes in livestock species compared to their wild ancestors. Understanding the evolutionary history and genetic diversity of livestock species is crucial, and mitochondrial DNA (mtDNA) has emerged as a robust marker for investigating molecular diversity in animals. Its highly conserved gene content across animal species, minimal duplications, absence of introns, and short intergenic regions make mtDNA analysis ideal for such studies. Mitochondrial DNA analysis has uncovered distinct cattle domestication events dating back to 8000 years BC in Southwestern Asia. The sequencing of water buffalo mtDNA in 2004 provided important insights into their domestication history. Caprine mtDNA analysis identified three haplogroups, indicating varied maternal origins. Sheep, domesticated 12,000 years ago, exhibit diverse mtDNA lineages, suggesting multiple domestication events. Ovine mtDNA studies revealed clades A, B, C, and a fourth lineage, group D. The origins of domestic pigs were traced to separate European and Asian events followed by interbreeding. In camels, mtDNA elucidated the phylogeographic structure and genetic differentiation between wild and domesticated species. Horses, domesticated around 3500 BC, show significant mtDNA variability, highlighting their diverse origins. Yaks exhibit unique adaptations for high-altitude environments, with mtDNA analysis providing insights into their adaptation. Chicken mtDNA studies supported a monophyletic origin from Southeast Asia's red jungle fowl, with evidence of multiple origins. This review explores livestock evolution and diversity through mtDNA studies, focusing on cattle, water buffalo, goat, sheep, pig, camel, horse, yak and chicken. It highlights mtDNA's significance in unraveling maternal lineages, genetic diversity, and domestication histories, concluding with insights into its potential application in improving livestock production and reproduction dynamics.

Repeat-Rich Regions Cause False-Positive Detection of NUMTs: A Case Study in Amphibians Using an Improved Cane Toad Reference Genome

Mitochondrial DNA (mtDNA) has been widely used in genetics research for decades. Contamination from nuclear DNA of mitochondrial origin (NUMTs) can confound studies of phylogenetic relationships and mtDNA heteroplasmy. Homology searches with mtDNA are widely used to detect NUMTs in the nuclear genome. Nevertheless, false-positive detection of NUMTs is common when handling repeat-rich sequences, while fragmented genomes might result in missing true NUMTs. In this study, we investigated different NUMT detection methods and how the quality of the genome assembly affects them. We presented an improved nuclear genome assembly (aRhiMar1.3) of the invasive cane toad (Rhinella marina) with additional long-read Nanopore and 10× linked-read sequencing. The final assembly was 3.47 Gb in length with 91.3% of tetrapod universal single-copy orthologs (n = 5,310), indicating the gene-containing regions were well assembled. We used 3 complementary methods (NUMTFinder, dinumt, and PALMER) to study the NUMT landscape of the cane toad genome. All 3 methods yielded consistent results, showing very few NUMTs in the cane toad genome. Furthermore, we expanded NUMT detection analyses to other amphibians and confirmed a weak relationship between genome size and the number of NUMTs present in the nuclear genome. Amphibians are repeat-rich, and we show that the number of NUMTs found in highly repetitive genomes is prone to inflation when using homology-based detection without filters. Together, this study provides an exemplar of how to robustly identify NUMTs in complex genomes when confounding effects on mtDNA analyses are a concern.

Mitochondrial heterogeneity and crosstalk in aging: Time for a paradigm shift?

The hallmarks of aging have been influential in guiding the biology of aging research, with more recent and growing recognition of the interdependence of these hallmarks on age-related health outcomes. However, a current challenge is personalizing aging trajectories to promote healthy aging, given the diversity of genotypes and lived experience. We suggest that incorporating heterogeneity-including intrinsic (e.g., genetic and structural) and extrinsic (e.g., environmental and exposome) factors and their interdependence of hallmarks-may move the dial. This editorial perspective will focus on one hallmark, namely mitochondrial dysfunction, to exemplify how consideration of heterogeneity and interdependence or crosstalk may reveal new perspectives and opportunities for personalizing aging research. To this end, we highlight heterogeneity within mitochondria as a model.

Whole-genome analysis reveals phylogenetic and demographic history of Eurasian perch

The contemporary diversity and distribution of species are shaped by their evolutionary and ecological history. This can be deciphered with the help of phylogenetic and demographic analysis methods, ideally combining and supplementing information from mitochondrial and nuclear genomes. In this study, we investigated the demographic history of Eurasian perch (Perca fluviatilis), a highly adaptable teleost with a distribution range across Eurasia. We combined whole-genome resequencing data with available genomic resources to analyse the phylogeny, phylogeography, and demographic history of P. fluviatilis populations from Europe and Siberia. We identified five highly diverged evolutionary mtDNA lineages, three of which show a strong signal of admixture in the Baltic Sea region. The estimated mean divergence time between these lineages ranged from 0.24 to 1.42 million years. Based on nuclear genomes, two distinct demographic trajectories were observed in European and Siberian samples reflecting contrasting demographic histories ca. 30,000-100,000 years before the present. A comparison of mtDNA and nuclear DNA evolutionary trees and AMOVA revealed concordances, as well as incongruences, between the two types of data, most likely reflecting recent postglacial colonization and hybridization events. Overall, our findings demonstrate the power and usefulness of genome-wide information for delineating historical processes that have shaped the genome of P. fluviatilis. We also highlight the added value of data-mining existing transcriptomic resources to complement novel sequence data, helping to shed light on putative glacial refugia and postglacial recolonization routes.

Mitochondrial DNA of the Arabian Camel Camelus dromedarius

The Camelidae family, ranging from southwest Asia to north Africa, South America, and Australia, includes key domesticated species adapted to diverse environments. Among these, the Arabian camel (Camelus dromedarius) is vital to the cultural and economic landscape of the Arabian Peninsula. This review explores the mitochondrial DNA of the dromedary camel, focusing on the D-loop region to understand its genetic diversity, maternal inheritance, and evolutionary history. We aim to investigate the unique characteristics of Arabian camel mtDNA, analyze the D-loop for genetic diversity and maternal lineage patterns, and explore the implications of mitochondrial genomic studies for camel domestication and adaptation. Key findings on mtDNA structure and variation highlight significant genetic differences and adaptive traits. The D-loop, essential for mtDNA replication and transcription, reveals extensive polymorphisms and haplotypes, providing insights into dromedary camel domestication and breeding history. Comparative analyses with other camelid species reveal unique genetic signatures in the Arabian camel, reflecting its evolutionary and adaptive pathways. Finally, this review integrates recent advancements in mitochondrial genomics, demonstrating camel genetic diversity and potential applications in conservation and breeding programs. Through comprehensive mitochondrial genome analysis, we aim to enhance the understanding of Camelidae genetics and contribute to the preservation and improvement of these vital animals.

mtDNA "nomenclutter" and its consequences on the interpretation of genetic data

Population-based studies of human mitochondrial genetic diversity often require the classification of mitochondrial DNA (mtDNA) haplotypes into more than 5400 described haplogroups, and further grouping those into hierarchically higher haplogroups. Such secondary haplogroup groupings (e.g., "macro-haplogroups") vary across studies, as they depend on the sample quality, technical factors of haplogroup calling, the aims of the study, and the researchers' understanding of the mtDNA haplogroup nomenclature. Retention of historical nomenclature coupled with a growing number of newly described mtDNA lineages results in increasingly complex and inconsistent nomenclature that does not reflect phylogeny well. This "clutter" leaves room for grouping errors and inconsistencies across scientific publications, especially when the haplogroup names are used as a proxy for secondary groupings, and represents a source for scientific misinterpretation. Here we explore the effects of phylogenetically insensitive secondary mtDNA haplogroup groupings, and the lack of standardized secondary haplogroup groupings on downstream analyses and interpretation of genetic data. We demonstrate that frequency-based analyses produce inconsistent results when different secondary mtDNA groupings are applied, and thus allow for vastly different interpretations of the same genetic data. The lack of guidelines and recommendations on how to choose appropriate secondary haplogroup groupings presents an issue for the interpretation of results, as well as their comparison and reproducibility across studies. To reduce biases originating from arbitrarily defined secondary nomenclature-based groupings, we suggest that future updates of mtDNA phylogenies aimed for the use in mtDNA haplogroup nomenclature should also provide well-defined and standardized sets of phylogenetically meaningful algorithm-based secondary haplogroup groupings such as "macro-haplogroups", "meso-haplogroups", and "micro-haplogroups". Ideally, each of the secondary haplogroup grouping levels should be informative about different human population history events. Those phylogenetically informative levels of haplogroup groupings can be easily defined using TreeCluster, and then implemented into haplogroup callers such as HaploGrep3. This would foster reproducibility across studies, provide a grouping standard for population-based studies, and reduce errors associated with haplogroup nomenclatures in future studies.

Forensic Characterization, Genomic Variability and Ancestry Analysis of Six Populations from Odisha Using mtDNA SNPs and Autosomal STRs

Located on India's eastern coast, Odisha is known for its diverse tribes and castes. In the early days of genome sequencing technology, researchers primarily studied the Austroasiatic communities inhabiting this region to reconstruct the ancient origins and dispersal of this broad linguistic group. However, current research has shifted towards identifying population and individual-specific genome variation for forensic applications. This study aims to analyze the forensic efficiency and ancestry of six populations from Odisha. We assessed the SF mtDNA-SNP60™ PCR Amplification Kit by comparing it with PowerPlex® Fusion 6C System, a widely used autosomal STR (aSTR) kit, in an Indian cohort. Although the mtDNA SNP kit showed low discriminating power for individuals of a diverse population, it could identify deep lineage divergence. Also, we utilized mitochondrial and autosomal variation information to analyze the ancestry of six endogamous ethnic groups in Odisha. We observe two extremities-populations with higher West Asian affinity and those with East Asian affinity. This observation is in congruence with the existing information of their tribal and non-tribal affiliation. When compared with neighbouring populations from Central and Eastern India, multivariate analysis showed that the Brahmins clustered separately or with the Gopala, Kaibarta appeared as an intermediate, Pana and Kandha clustered with the Gonds, and Savara with the Munda tribes. Our findings indicate significant deep lineage stratification in the ethnic populations of Odisha and a gene flow from West and East Asia. The artefacts of unique deep lineage in such a diverse population will help in improving forensic identification. In addition, we conclude that the SF mtDNA-SNP60 PCR Amplification Kit may be used only as a supplementary tool for forensic analysis.

Convergent Mutations and Single Nucleotide Variants in Mitochondrial Genomes of Modern Humans and Neanderthals

The genetic contributions of Neanderthals to the modern human genome have been evidenced by the comparison of present-day human genomes with paleogenomes. Neanderthal signatures in extant human genomes are attributed to intercrosses between Neanderthals and archaic anatomically modern humans (AMHs). Although Neanderthal signatures are well documented in the nuclear genome, it has been proposed that there is no contribution of Neanderthal mitochondrial DNA to contemporary human genomes. Here we show that modern human mitochondrial genomes contain 66 potential Neanderthal signatures, or Neanderthal single nucleotide variants (N-SNVs), of which 36 lie in coding regions and 7 result in nonsynonymous changes. Seven N-SNVs are associated with traits such as cycling vomiting syndrome, Alzheimer's disease and Parkinson's disease, and two N-SNVs are associated with intelligence quotient. Based on recombination tests, principal component analysis (PCA) and the complete absence of these N-SNVs in 41 archaic AMH mitogenomes, we conclude that convergent evolution, and not recombination, explains the presence of N-SNVs in present-day human mitogenomes.

Mitochondrial DNA variability and Covid-19 in the Slovak population

Recent studies have shown that mitochondria are involved in the pathogenesis of Covid-19. Mitochondria play a role in production of reactive oxygen species and induction of an innate immune response, both important during infections. Common variability of mitochondrial DNA (mtDNA) can affect oxidative phosphorylation and the risk or lethality of cardiovascular, neurodegenerative diseases and sepsis. However, it is unclear whether susceptibility of severe Covid-19 might be affected by mtDNA variation. Thus, we have analyzed mtDNA in a sample of 446 Slovak patients hospitalized due to Covid-19 and a control population group consisting of 1874 individuals. MtDNA variants in the HVRI region have been analyzed and classified into haplogroups at various phylogenetic levels. Binary logistic regression was used to assess the risk of Covid-19. Haplogroups T1, H11, K and variants 16256C > T, 16265A > C, 16293A > G, 16311 T > C and 16399A > G were associated with an increased Covid-19 risk. On contrary, Haplogroup J1, haplogroup clusters H + U5b and T2b + U5b, and the mtDNA variant 16189 T > C were associated with decreased risk of Covid-19. Following the application of the Bonferroni correction, statistical significance was observed exclusively for the cluster of haplogroups H + U5b. Unsurprisingly, the most significant factor contributing to the mortality of patients with Covid-19 is the age of patients. Our findings suggest that mtDNA haplogroups can play a role in Covid-19 pathogenesis, thus potentially useful in identifying susceptibility to its severe form. To confirm these associations, further studies taking into account the nuclear genome or other non-biological influences are needed.

Unraveling the mitochondrial phylogenetic landscape of Thailand reveals complex admixture and demographic dynamics

The evolutionary dynamics of mitochondrial DNA within the Thai population were comprehensively explored with a specific focus on the influence of South Asian admixture. A total of 166 samples were collected through randomized sampling, ensuring a diverse representation. Our findings unveil substantial genetic and haplogroup diversity within the Thai population. We have identified 164 haplotypes categorized into 97 haplogroups, with a notable inclusion of 20 novel haplogroups. The distribution of haplogroups exhibited variations across different populations and countries. The central Thai population displayed a high diversity of haplogroups from both the M and N clades. Maternal lineage affinities were discerned between several Mainland Southeast Asia (MSEA) and South Asian populations, implying ancestral genetic connections and a substantial influence of South Asian women in establishing these relationships. f4-statistics indicates the presence of a Tibeto-Burman genetic component within the Mon population from Thailand. New findings demonstrate two phases of population expansion occurring 22,000-26,000 and 2500-3800 years ago, coinciding with the Last Glacial Maximum, and Neolithic demographic transition, respectively. This research significantly enhances our understanding of the maternal genetic history of Thailand and MSEA, emphasizing the influence of South Asian admixture. Moreover, it underscores the critical role of prior information, such as mutation rates, within the Bayesian framework for accurate estimation of coalescence times and inferring demographic history.

Exploring the legacy of African and Indigenous Caribbean admixture in Puerto Rico

OBJECTIVES

From an anthropological genetic perspective, little is known about the ethnogenesis of African descendants in Puerto Rico. Furthermore, historical interactions between Indigenous Caribbean and African descendant peoples that may be reflected in the ancestry of contemporary populations are understudied. Given this dearth of genetic research and the precedence for Afro-Indigenous interactions documented by historical, archeological, and other lines of evidence, we sought to assess the biogeographic origins of African descendant Puerto Ricans and to query the potential for Indigenous ancestry within this community.

MATERIALS AND METHODS

Saliva samples were collected from 58 self-identified African descendant Puerto Ricans residing in Puerto Rico. We sequenced whole mitochondrial genomes and genotyped Y chromosome haplogroups for each male individual (n = 25). Summary statistics, comparative analyses, and network analysis were used to assess diversity and variation in haplogroup distribution between the sample and comparative populations.

RESULTS

As indicated by mitochondrial haplogroups, 66% had African, 5% had European, and 29% had Indigenous American matrilines. Along the Y chromosome, 52% had African, 28% had Western European, 16% had Eurasian, and, notably, 4% had Indigenous American patrilines. Both mitochondrial and Y chromosome haplogroup frequencies were significantly different from several comparative populations.

DISCUSSION

Biogeographic origins are consistent with historical accounts of African, Indigenous American, and European ancestry. However, this first report of Indigenous American paternal ancestry in Puerto Rico suggests distinctive features within African descendant communities on the island. Future studies expanding sampling and incorporating higher resolution genetic markers are necessary to more fully understand African descendant history in Puerto Rico.

High Mitochondrial Haplotype Diversity Found in Three Pre-Hispanic Groups from Colombia

The analysis of mitochondrial DNA (mtDNA) hypervariable region (HVR) sequence data from ancient human remains provides valuable insights into the genetic structure and population dynamics of ancient populations. mtDNA is particularly useful in studying ancient populations, because it is maternally inherited and has a higher mutation rate compared to nuclear DNA. To determine the genetic structure of three Colombian pre-Hispanic populations and compare them with current populations, we determined the haplotypes from human bone remains by sequencing several mitochondrial DNA segments. A wide variety of mitochondrial polymorphisms were obtained from 33 samples. Our results support a high population heterogeneity among pre-Hispanic populations in Colombia.

Selection and Interpretation of Molecular Diagnostics in Heart Transplantation

The number of heart transplants performed annually in the United States and worldwide continues to increase, but there has been little change in graft longevity and patient survival over the past 2 decades. The reference standard for diagnosis of acute cellular and antibody-mediated rejection includes histologic and immunofluorescence evaluation of endomyocardial biopsy samples, despite invasiveness and high interrater variability for grading histologic rejection. Circulating biomarkers and molecular diagnostics have shown substantial predictive value in rejection monitoring, and emerging data support their use in diagnosing other posttransplant complications. The use of genomic (cell-free DNA), transcriptomic (mRNA and microRNA profiling), and proteomic (protein expression quantitation) methodologies in diagnosis of these posttransplant outcomes has been evaluated with varying levels of evidence. In parallel, growing knowledge about the genetically mediated immune response leading to rejection (immunogenetics) has enhanced understanding of antibody-mediated rejection, associated graft dysfunction, and death. Antibodies to donor human leukocyte antigens and the technology available to evaluate these antibodies continues to evolve. This review aims to provide an overview of biomarker and immunologic tests used to diagnose posttransplant complications. This includes a discussion of pediatric heart transplantation and the disparate rates of rejection and death experienced by Black patients receiving a heart transplant. This review describes diagnostic modalities that are available and used after transplant and the landscape of future investigations needed to enhance patient outcomes after heart transplantation.

Establishing a Sequencing Method for the Whole Mitochondrial DNA of Domestic Dogs

In human beings, whole mitochondrial DNA (mtDNA) sequencing has been widely used in many research fields, including medicine, forensics, and genetics. With respect to the domestic dog (Canis lupus familiaris), which is commonly recognized as being an additional member of the traditional human family structure, research studies on mtDNA should be developed to expand and improve our collective knowledge of dog medicine and welfare as it seems that there is still room for further development in these areas. Moreover, a simple and robust method for sequencing whole mtDNA that can be applied to various dog breeds has not yet been described in the literature. In the present study, we aim to establish such a method for the whole mtDNA sequencing of the domestic dog. In the experiments we conducted, oral mucosa DNA samples obtained from six Japanese domestic dogs were used as a template. We designed four primer pairs that could amplify approximately 5 kbp from each region of the mtDNA and validated several PCR conditions. Subsequently, the PCR amplicons were pooled and subjected to library preparation. The sequencing of the libraries was performed using next-generation sequencing (NGS), followed by bioinformatics analysis. Our results demonstrate that the proposed method can be used to perform highly accurate resequencing. We believe that this method may be useful for future research conducted to better understand dog medicine and welfare.

Population analysis of complete mitogenomes for 334 samples from El Salvador

The use of mitochondrial DNA (mtDNA) in the field of forensic genetics is widely spread mainly due to its advantages when identifying highly degraded samples. In this sense, massive parallel sequencing has made the analysis of the whole mitogenome more accessible, noticeably increasing the informativeness of mtDNA haplotypes. The civil war (1980-1992) in El Salvador caused many deaths and disappearances (including children) all across the country and the economic and social instability after the war forced many people to emigration. For this reason, different organizations have collected DNA samples from relatives with the aim of identifying missing people. Thus, we present a dataset containing 334 complete mitogenomes from the Salvadoran general population. To the best of our knowledge, this is the first publication of a nationwide forensic-quality complete mitogenome database of any Latin American country. We found 293 different haplotypes, with a random match probability of 0.0041 and 26.6 mean pairwise differences, which is similar to other Latin American populations, and which represent a marked improvement from the values obtained with just control region sequences. These haplotypes belong to 54 different haplogroups, being 91% of them of Native American origin. Over a third (35.9%) of the individuals carried at least a heteroplasmic site (excluding length heteroplasmies). Ultimately, the present database aims to represent mtDNA haplotype diversity in the general Salvadoran populations as a basis for the identification of people that disappeared during or after the civil war.

Genetic features of Sri Lankan elephant, Elephas maximus maximus Linnaeus revealed by high throughput sequencing of mitogenome and ddRAD-seq

Elephas maximus maximus Linnaeus, the Sri Lankan subspecies is the largest and the darkest among Asian elephants. Patches of depigmented areas with no skin color on the ears, face, trunk, and belly morphologically differentiate it from the others. The elephant population in Sri Lanka is now limited to smaller areas and protected under Sri Lankan law. Despite its ecological and evolutionary importance, the relationship between Sri Lankan elephants and their phylogenetic position among Asian elephants remains controversial. While identifying genetic diversity is the key to any conservation and management strategies, limited data is currently available. To address such issues, we analyzed 24 elephants with known parental lineages with high throughput ddRAD-seq. The mitogenome suggested the coalescence time of the Sri Lankan elephant at ~0.2 million years, and sister to Myanmar elephants supporting the hypothesis of the movement of elephants in Eurasia. The ddRAD-seq approach identified 50,490 genome-wide SNPs among Sri Lankan elephants. The genetic diversity within Sri Lankan elephants assessed with identified SNPs suggests a geographical differentiation resulting in three main clusters; north-eastern, mid-latitude, and southern regions. Interestingly, though it was believed that elephants from the Sinharaja rainforest are of an isolated population, the ddRAD-based genetic analysis clustered it with the north-eastern elephants. The effect of habitat fragmentation on genetic diversity could be further assessed with more samples with specific SNPs identified in the current study.

HaploCart: Human mtDNA haplogroup classification using a pangenomic reference graph human mtDNA haplogroup inference

Current mitochondrial DNA (mtDNA) haplogroup classification tools map reads to a single reference genome and perform inference based on the detected mutations to this reference. This approach biases haplogroup assignments towards the reference and prohibits accurate calculations of the uncertainty in assignment. We present HaploCart, a probabilistic mtDNA haplogroup classifier which uses a pangenomic reference graph framework together with principles of Bayesian inference. We demonstrate that our approach significantly outperforms available tools by being more robust to lower coverage or incomplete consensus sequences and producing phylogenetically-aware confidence scores that are unbiased towards any haplogroup. HaploCart is available both as a command-line tool and through a user-friendly web interface. The C++ program accepts as input consensus FASTA, FASTQ, or GAM files, and outputs a text file with the haplogroup assignments of the samples along with the level of confidence in the assignments. Our work considerably reduces the amount of data required to obtain a confident mitochondrial haplogroup assignment.

Mitochondrial Deoxyribonucleic Acid (mtDNA), Maternal Inheritance, and Their Role in the Development of Cancers: A Scoping Review

Mitochondrial DNA (mtDNA) is a small, circular, double-stranded DNA inherited from the mother during fertilization. Evolutionary evidence supported by the endosymbiotic theory identifies mitochondria as an organelle that could have descended from prokaryotes. This may be the reason for the independent function and inheritance pattern shown by mtDNA. The unstable nature of mtDNA due to the lack of protective histones, and effective repair systems make it more vulnerable to mutations. The mtDNA and its mutations could be maternally inherited thereby predisposing the offspring to various cancers like breast and ovarian cancers among others. Although mitochondria are considered heteroplasmic wherein variations among the multiple mtDNA genomes are noticed, mothers can have mitochondrial populations that are homoplasmic for a given mitochondrial mutation. Homoplasmic mitochondrial mutations may be transmitted to all maternal offspring. However, due to the complex interplay between the mitochondrial and nuclear genomes, it is often difficult to predict disease outcomes, even with homoplasmic mitochondrial populations. Heteroplasmic mtDNA mutations can be maternally inherited, but the proportion of mutated alleles differs markedly between offspring within one generation. This led to the genetic bottleneck hypothesis, explaining the rapid changes in allele frequency witnessed during the transmission of mtDNA from one generation to the next. Although a physical reduction in mtDNA has been demonstrated in several species, a comprehensive understanding of the molecular mechanisms is yet to be demonstrated. Despite initially thought to be limited to the germline, there is evidence that blockages exist in different cell types during development, perhaps explaining why different tissues in the same organism contain different levels of mutated mtDNA. In this review, we comprehensively discuss the potential mechanisms through which mtDNA undergoes mutations and the maternal mode of transmission that contributes to the development of tumors, especially breast and ovarian cancers.

African mitochondrial haplogroup L7: a 100,000-year-old maternal human lineage discovered through reassessment and new sequencing

Archaeological and genomic evidence suggest that modern Homo sapiens have roamed the planet for some 300–500 thousand years. In contrast, global human mitochondrial (mtDNA) diversity coalesces to one African female ancestor (“Mitochondrial Eve”) some 145 thousand years ago, owing to the ¼ gene pool size of our matrilineally inherited haploid genome. Therefore, most of human prehistory was spent in Africa where early ancestors of Southern African Khoisan and Central African rainforest hunter-gatherers (RFHGs) segregated into smaller groups. Their subdivisions followed climatic oscillations, new modes of subsistence, local adaptations, and cultural-linguistic differences, all prior to their exodus out of Africa. Seven African mtDNA haplogroups (L0–L6) traditionally captured this ancient structure—these L haplogroups have formed the backbone of the mtDNA tree for nearly two decades. Here we describe L7, an eighth haplogroup that we estimate to be ~ 100 thousand years old and which has been previously misclassified in the literature. In addition, L7 has a phylogenetic sublineage L7a*, the oldest singleton branch in the human mtDNA tree (~ 80 thousand years). We found that L7 and its sister group L5 are both low-frequency relics centered around East Africa, but in different populations (L7: Sandawe; L5: Mbuti). Although three small subclades of African foragers hint at the population origins of L5'7, the majority of subclades are divided into Afro-Asiatic and eastern Bantu groups, indicative of more recent admixture. A regular re-estimation of the entire mtDNA haplotype tree is needed to ensure correct cladistic placement of new samples in the future.

Characterization and Phylogenetic Analyses of the Complete Mitochondrial Genome of Sugarcane (Saccharum spp. Hybrids) Line A1

Modern sugarcane cultivars are highly polyploid with complex nuclear genomic genetic background, while their mitochondrion (mt) genomes are much simpler, smaller and more manageable and could provide useful phylogenetic information. In this study, the mt genome of a modern commercial cultivar A1 was sequenced via Illumina Hiseq XTen and PacBio Sequel platform. The assembled and annotated mitochondrial genomes of A1 were composed of two circular DNA molecules, one large and one small, which were named Chromosome 1 and Chromosome 2. The two distinct circular chromosomes of mitogenome construct is consisted with other sugarcane cultivars i.e., Saccharum officinarum Khon Kaen 3 and Saccharum spp. hybrids ROC22 and FN15. The Chromosome 1 of A1 mitogenome is 300,822 bp in length with the GC content of 43.94%, and 7.14% of Chromosome 1 sequences (21,468 nucleotides) are protein coding genes (PCGs) while 92.86% (279,354 nucleotides) are intergenic region. The length of Chromosome 2 is 144,744 bp with the GC content of 43.57%, and 8.20% of Chromosome 2 sequences (11,865 nucleotides) are PCGs while 91.80% (132,879 nucleotides) are intergenic region. A total of 43 genes are located on Chromosome 1, which contains 22 PCGs (six nad genes, four rps genes, four atp genes, three ccm genes, three cox genes, one mat gene and one mtt gene) and 21 non-coding genes including 15 tRNAs and 6 rRNAs. Chromosome 2 includes 18 genes in total, which contains 13 PCGs (four nad genes, three rps genes, two atp genes, one ccm gene, one cob gene, one cox gene and one rpl gene) and five non-coding genes (tRNA genes). Analysis of codon usage of 35 PCGs showed that codon ending in A/U was preferred. Investigation of gene composition indicated that the types and copy numbers of CDS genes, tRNAs and rRNAs of A1 and FN15 were identical. The cox1 gene has two copies and the trnP gene has one copy in A1, FN15 and ROC22 three lines, while there is only one copy of cox1 and two copies of trnP in S. officinarum Khon Kaen 3. In addition, S. officinarum Khon Kaen 3 have no nad1 gene and rps7 gene. 100 sequence repeats, 38 SSRs and 444 RNA editing sites in A1 mt genome were detected. Moreover, the maximum likelihood phylogenetic analysis found that A1 were more closely related to S. spp. hybrid (ROC22 and FN15) and S. officinarum (Khon Kaen 3). Herein, the complete mt genome of A1 will provide essential DNA molecular information for further phylogenetic and evolutionary analysis for Saccharum and Poaceae.

Mitochondrial DNA Consensus Calling and Quality Filtering for Constructing Ancient Human Mitogenomes: Comparison of Two Widely Applied Methods

Retrieving high-quality endogenous ancient DNA (aDNA) poses several challenges, including low molecular copy number, high rates of fragmentation, damage at read termini, and potential presence of exogenous contaminant DNA. All these factors complicate a reliable reconstruction of consensus aDNA sequences in reads from high-throughput sequencing platforms. Here, we report findings from a thorough evaluation of two alternative tools (ANGSD and schmutzi) aimed at overcoming these issues and constructing high-quality ancient mitogenomes. Raw genomic data (BAM/FASTQ) from a total of 17 previously published whole ancient human genomes ranging from the 14th to the 7th millennium BCE were retrieved and mitochondrial consensus sequences were reconstructed using different quality filters, with their accuracy measured and compared. Moreover, the influence of different sequence parameters (number of reads, sequenced bases, mean coverage, and rate of deamination and contamination) as predictors of derived sequence quality was evaluated. Complete mitogenomes were successfully reconstructed for all ancient samples, and for the majority of them, filtering substantially improved mtDNA consensus calling and haplogroup prediction. Overall, the schmutzi pipeline, which estimates and takes into consideration exogenous contamination, appeared to have the edge over the much faster and user-friendly alternative method (ANGSD) in moderate to high coverage samples (>1,000,000 reads). ANGSD, however, through its read termini trimming filter, showed better capabilities in calling the consensus sequence from low-quality samples. Among all the predictors of overall sample quality examined, the strongest correlation was found for the available number of sequence reads and bases. In the process, we report a previously unassigned haplogroup (U3b) for an Early Chalcolithic individual from Southern Anatolia/Northern Levant.

A new linear combination method of haplogroup distribution central vectors to model population admixtures

We introduce a novel population genetic approach suitable to model the origin and relationships of populations, using new computation methods analyzing Hg frequency distributions. Hgs were selected into groups which show correlated frequencies in subsets of populations, based on the assumption that correlations were established in ancient separation, migration and admixture processes. Populations are defined with this universal Hg database, then using unsupervised artificial intelligence, central vectors (CVs) are determined from local condensations of the Hg-distribution vectors in the multidimensional point system. Populations are clustered according to their proximity to CVs. We show that CVs can be regarded as approximations of ancient populations and real populations can be modeled as weighted linear combinations of the CVs using a new linear combination algorithm based on a gradient search for the weights. The efficacy of the method is demonstrated by comparing Copper Age populations of the Carpathian Basin to Middle Age ones and modern Hungarians. Our analysis reveals significant population continuity since the Middle Ages, and the presence of a substrate component since the Copper Age.

Circum-Saharan Prehistory through the Lens of mtDNA Diversity

African history has been significantly influenced by the Sahara, which has represented a barrier for migrations of all living beings, including humans. Major exceptions were the gene flow events that took place between North African and sub-Saharan populations during the so-called African Humid Periods, especially in the Early Holocene (11.5 to 5.5 thousand years ago), and more recently in connection with trans-Saharan commercial routes. In this study, we describe mitochondrial DNA (mtDNA) diversity of human populations from both sides of the Sahara Desert, i.e., both from North Africa and the Sahel/Savannah belt. The final dataset of 7213 mtDNA sequences from 134 African populations encompasses 470 newly collected and 6743 previously published samples, which were analyzed using descriptive methods and Bayesian statistics. We completely sequenced 26 mtDNAs from sub-Saharan samples belonging to the Eurasian haplogroup N1. Analyses of these N1 mitogenomes revealed their possible routes to the Sahel, mostly via Bab el-Mandab. Our results indicate that maternal gene flow must have been important in this circum-Saharan space, not only within North Africa and the Sahel/Savannah belt but also between these two regions.

Interactions between nuclear and mitochondrial SNPs and Parkinson's disease risk

Interactions between the products of the nuclear and mitochondrial genomes are critical for the function of most eukaryotic cells. Recently the introduction of mitochondrial replacement therapy has raised the question of incompatibilities between mitochondrial and nuclear variants, and their potential influence on the genetic makeup of human populations. Such interactions could also contribute to the variability of the penetrance of pathogenic DNA variants. This led us to investigate the frequencies of combinations of nuclear and mitochondrial SNP alleles (mitonuclear combinations) in healthy individuals (n = 5375) and in a cohort of patients with Parkinson's disease (PD, n = 2210). In the unaffected population, we were not able to find associations between nuclear and mitochondrial variants with a false discovery rate below 0.05 after accounting for multiple testing (i.e., the number of combinations examined). However, in the PD cohort, five combinations surpassed this threshold. Next, after combining both cohorts, we investigated whether these associations were modulated by disease status. All five combinations were significant (p < 10-3 for all tests). These combinations also showed significant evidence for an effect of the interaction between the mitochondrial and nuclear variants on disease risk. Their nuclear components mapped to TBCA, NIBAN3, and GLT25D1 and an uncharacterised intergenic region. In summary, starting from a single cohort design we identified combinations of nuclear and mitochondrial variants affecting PD disease risk.

Structural variation of mitochondrial genomes sheds light on evolutionary history of soybeans

The architecture and genetic diversity of mitogenome (mtDNA) are largely unknown in cultivated soybean (Glycine max), which is domesticated from the wild progenitor, Glycine soja, 5000 years ago. Here, we de novo assembled the mitogenome of the cultivar 'Williams 82' (Wm82_mtDNA) with Illumina PE300 deep sequencing data, and verified it with polymerase chain reaction (PCR) and Southern blot analyses. Wm82_mtDNA maps as two autonomous circular chromosomes (370 871-bp Chr-m1 and 62 661-bp Chr-m2). Its structure is extensively divergent from that of the mono-chromosomal mitogenome reported in the landrace 'Aiganhuang' (AGH_mtDNA). Synteny analysis showed that the structural variations (SVs) between two genomes are mainly attributed to ectopic and illegitimate recombination. Moreover, Wm82_mtDNA and AGH_mtDNA each possess six and four specific regions, which are absent in their counterparts and likely result from differential sequence-loss events. Mitogenome SV was further studied in 39 wild and 182 cultivated soybean accessions distributed world-widely with PCR/Southern analyses or a comparable in silico analysis. The results classified both wild and cultivated soybeans into five cytoplasmic groups, named as GSa-GSe and G1-G5; 'Williams 82' and 'Aiganhuang' belong to G1 and G5, respectively. Notably, except for members in GSe and G5, all accessions carry a bi-chromosomal mitogenome with a common Chr-m2. Phylogenetic analyses based on mtDNA structures and chloroplast gene sequences both inferred that G1-G3, representing >90% of cultigens, likely inherited cytoplasm from the ancestor of domestic soybean, while G4 and G5 likely inherited cytoplasm from wild soybeans carrying GSa- and GSe-like cytoplasm through interspecific hybridization, offering new insights into soybean cultivation history.

Mitochondrial DNA haplogroup study: residents of Sulaymaniyah city in the Iraqi Kurdistan Region may be genetically closer to European lineage

Background

Being the native inhabitants of the Neolithic Fertile Crescent, Kurds were included in several maternal lineage studies concerning the Eurasian population. However, no study was performed on the Kurdish population of Sulaymaniyah city (latitude 33.314690 and longitude 44.376759). This study was carried out on a sample of Sorani Kurds living in Sulaymaniyah for the identification of population-related single nucleotide polymorphisms (SNPs) and modes of maternal lineage.

Results

In this study, 36 randomly selected healthy unrelated Kurdish subjects were enrolled. Whole mitochondrial DNA sequencing was performed. HaploGrep 2.0 and neutrality test (Tajima’s D) were employed for haplogroup identification and historical demography determination. When the outcomes were compared with previous studies in Kurds and the neighbouring nations, the identified haplogroups in the sample of study were members of the Western Eurasian haplogroups with a predominance of haplogroup H.

Conclusions

The whole mitochondrial DNA sequence is superior to the traditional analysis of the non-coding (control) region. Our study indicates a stronger relation of the studied group to the European lineage than to their neighbouring nations.

Genome-wide autosomal, mtDNA, and Y chromosome analysis of King Bela III of the Hungarian Arpad dynasty

The ancient Hungarians, "Madzsars", established their control of the Carpathian Basin in the late ninth century and founded the Hungarian Kingdom around 1000AD. The origin of the Magyars as a tribal federation has been much debated in the past. From the time of the conquest to the early fourteenth century they were ruled by descendants of the Arpad family. In order to learn more about the genetic origin of this family, we here analyzed the genome of Bela III one of the most prominent members of the early Hungarian dynasty that ruled the Hungarian Kingdom from 1172 to 1196. The Y-Chromosome of Bela III belongs to haplogroup R1a-Z2123 that is today found in highest frequency in Central Asia, supporting a Central Asian origin for the ruling lineage of the Hungarian kingdom. The autosomal DNA profile of Bela III, however, falls within the genetic variation of present-day east European populations. This is further supported through his mtDNA genome that belongs to haplogroup H, the most common European maternal lineage, but also found in Central Asia. However, we didn't find an exact haplotype match for Bela III. The typical autosomal and maternal Central Eastern European ancestry among Bela III autosomes might be best explained by consecutive intermarriage with local European ruling families.

Pilot Screening of Cell-Free mtDNA in NIPT: Quality Control, Variant Calling, and Haplogroup Determination

Clinical tests based on whole-genome sequencing are generally focused on a single task approach, testing one or several parameters, although whole-genome sequencing (WGS) provides us with large data sets that can be used for many supportive analyses. In spite of low genome coverage, data of WGS-based non-invasive prenatal testing (NIPT) contain fully sequenced mitochondrial DNA (mtDNA). This mtDNA can be used for variant calling, ancestry analysis, population studies and other approaches that extend NIPT functionality. In this study, we analyse mtDNA pool from 645 cell-free DNA (cfDNA) samples of pregnant women from different regions of Russia, explore the effects of transportation and storing conditions on mtDNA content, analyse effects, frequency and location of mitochondrial variants called from samples and perform haplogroup analysis, revealing the most common mitochondrial superclades. We have shown that, despite the relatively low sequencing depth of unamplified mtDNA from cfDNA samples, the mtDNA analysis in these samples is still an informative instrument suitable for research and screening purposes.

Association of whole mtDNA, an NADPH G11914A variant, and haplogroups with high physical performance in an elite military troop

Physical performance is a multifactorial and complex trait influenced by environmental and hereditary factors. Environmental factors alone have been insufficient to characterize all outstanding phenotypes. Recent advances in genomic technologies have enabled the investigation of whole nuclear and mitochondrial genome sequences, increasing our ability to understand interindividual variability in physical performance. Our objective was to evaluate the association of mitochondrial polymorphic loci with physical performance in Brazilian elite military personnel. Eighty-eight male military personnel who participated in the Command Actions Course of the Army were selected. Total DNA was obtained from blood samples and a complete mitochondrial genome (mtDNA) was sequenced using Illumina MiSeq platform. Twenty-nine subjects completed the training program (FINISHED, 'F'), and fifty-nine failed to complete (NOT_FINISHED, 'NF'). The mtDNA from NF was slightly more similar to genomes from African countries frequently related to endurance level. Twenty-two distinct mtDNA haplogroups were identified corroborating the intense genetic admixture of the Brazilian population, but their distribution was similar between the two groups (FST=0.0009). Of 745 polymorphisms detected in the mtDNA, the position G11914A within the NADPH gene component of the electron transport chain, was statistically different between F and NF groups (P=0.011; OR: 4.286; 95%CI: 1.198-16.719), with a higher frequency of the G allele in group F individuals). The high performance of military personnel may be mediated by performance-related genomic traits. Thus, mitochondrial genetic markers such as the ND4 gene may play an important role on physical performance variability.

Ancient and modern mitogenomes from Central Argentina: new insights into population continuity, temporal depth and migration in South America

The inverted triangle shape of South America places Argentina territory as a geographical crossroads between the two principal peopling streams that followed either the Pacific or the Atlantic coasts, which could have then merged in Central Argentina (CA). Although the genetic diversity from this region is therefore crucial to decipher past population movements in South America, its characterization has been overlooked so far. We report 92 modern and 22 ancient mitogenomes spanning a temporal range of 5000 years, which were compared with a large set of previously reported data. Leveraging this dataset representative of the mitochondrial diversity of the subcontinent, we investigate the maternal history of CA populations within a wider geographical context. We describe a large number of novel clades within the mitochondrial DNA tree, thus providing new phylogenetic interpretations for South America. We also identify several local clades of great temporal depth with continuity until the present time, which stem directly from the founder haplotypes, suggesting that they originated in the region and expanded from there. Moreover, the presence of lineages characteristic of other South American regions reveals the existence of gene flow to CA. Finally, we report some lineages with discontinuous distribution across the Americas, which suggest the persistence of relic lineages likely linked to the first population arrivals. The present study represents to date the most exhaustive attempt to elaborate a Native American genetic map from modern and ancient complete mitochondrial genomes in Argentina and provides relevant information about the general process of settlement in South America.

Mitochondrial genome architecture and phylogenetic relationships of Odontesthes argentinensis within Atherinomorpha

Silversides are a widely distributed group across South America, with several species occupying marine, freshwater and estuarine environments. Several authors suggest main transitions among these environments took place during Pleistocene, and were accompanied with rapid speciation events. This scenario produced very limited genetic and morphological differentiation among the species. However, most of these surveys have an incomplete coverage of the intraspecific genetic diversity of the taxa studied. In this work, we reconstructed six mitochondrial genomes of O. argentinensis using transcriptomic data, and used them-in combination with several nuclear markers retrieved from the same transcriptomes-to explore the effect of additional coverage of intraspecific diversity of this species in phylogenetic reconstructions. Unlike previous works, phylogenetic analyses failed to identify O. argentinensis as a monophyletic group in relation with closely related taxa. Our results suggest that several species of the genus, especially those related to O. argentinensis, need further taxonomic revision.

Ancestry and admixture of a southernmost Chilean population: The reflection of a migratory history

OBJECTIVES

Punta Arenas is a Chilean city situated on ancestral Aönikenk territory. The city was founded by 19th- and 20th-century colonists from Chile (Chiloé) and Europe (Croatia). This work uses uniparental and ancestry-informative markers (AIMs) to explore the effects of historic migratory and admixture patterns on the current genetic composition of Punta Arenas.

METHODS

We analyzed mitochondrial DNA (mtDNA), Y-chromosome single-nucleotide polymorphisms (SNPs), and 141 AIMs obtained from 129 DNA samples from male residents with regional ancestry. After characterizing uniparental lineages and ancestry proportions, multivariate analysis was used to explore relationships among the various types of data.

RESULTS

Punta Arenas has an admixed population with three main genetic components: European (56.5%), northern Native (11.3%), and south-central Native (28.6%). The Native component is preponderant in the mtDNA (83.76%), while the foreign component predominates in the Y-chromosome (92.25%). Non-Native mtDNA lineages are associated with European genetic ancestry, and Native mtDNA lineages originated mainly in the southern and southernmost regions of Chile. Most non-Native Y-chromosome SNPs originated in Spain, and secondly, in Croatia.

CONCLUSIONS

The population of Punta Arenas is mainly of Chilote origin with south-central Native and Spanish ancestral components, as well as some Croatian components. The persistence of local Native lineages is notable, suggesting continuity with the ancestral populations of the region such as the Kawésqar, Aönikenk, Yámana, or Selknam peoples. This study contributes to our knowledge of local history and its links to national and global developments in genetic ancestry.

Maternal Lineages from 10-11th Century Commoner Cemeteries of the Carpathian Basin

Nomadic groups of conquering Hungarians played a predominant role in Hungarian prehistory, but genetic data are available only from the immigrant elite strata. Most of the 10-11th century remains in the Carpathian Basin belong to common people, whose origin and relation to the immigrant elite have been widely debated. Mitogenome sequences were obtained from 202 individuals with next generation sequencing combined with hybridization capture. Median joining networks were used for phylogenetic analysis. The commoner population was compared to 87 ancient Eurasian populations with sequence-based (Fst) and haplogroup-based population genetic methods. The haplogroup composition of the commoner population markedly differs from that of the elite, and, in contrast to the elite, commoners cluster with European populations. Alongside this, detectable sub-haplogroup sharing indicates admixture between the elite and the commoners. The majority of the 10-11th century commoners most likely represent local populations of the Carpathian Basin, which admixed with the eastern immigrant groups (which included conquering Hungarians).

Population inference based on mitochondrial DNA control region data by the nearest neighbors algorithm

Population and geographic assignment are frequently undertaken using DNA sequences on the mitochondrial genome. Assignment to broad continental populations is common, although finer resolution to subpopulations can be less accurate due to shared genetic ancestry at a local level and members of different ancestral subpopulations cohabiting the same geographic area. This study reports on the accuracy of population and subpopulation assignment by using the sequence data obtained from the 3070 mitochondrial genomes and applying the K-nearest neighbors (KNN) algorithm. These data also included training samples used for continental and population assignment comprised of 1105 Europeans (including Austria, France, Germany, Spain, and England and Caucasian countries), 374 Africans (including North and East Africa and non-specific area (Pan-Africa)), and 1591 Asians (including Japan, Philippines, and Taiwan). Subpopulations included in this study were 1153 mitochondrial DNA (mtDNA) control region sequences from 12 subpopulations in Taiwan (including Han, Hakka, Ami, Atayal, Bunun, Paiwan, Puyuma, Rukai, Saisiyat, Tsou, Tao, and Pingpu). Additionally, control region sequence data from a further 50 samples, obtained from the Sigma Company, were included after they were amplified and sequenced. These additional 50 samples acted as the "testing samples" to verify the accuracy of the population. In this study, based on genetic distances as genetic metric, we used the KNN algorithm and the K-weighted-nearest neighbors (KWNN) algorithm weighted by genetic distance to classify individuals into continental populations, and subpopulations within the same continent. Accuracy results of ethnic inferences at the level of continental populations and of subpopulations among KNN and KWNN algorithms were obtained. The training sample set achieved an overall accuracy of 99 to 82% for assignment to their continental populations with K values from 1 to 101. Population assignment for subpopulations with K assignments from 1 to 5 reached an accuracy of 77 to 54%. Four out of 12 Taiwanese populations returned an accuracy of assignment of over 60%, Ami (66%), Atayal (67%), Saisiyat (66%), and Tao (80%). For the testing sample set, results of ethnic prediction for continental populations with recommended K values as 5, 10, and 35, based on results of the training sample set, achieved overall an accuracy of 100 to 94%. This study provided an accurate method in population assignment for not only continental populations but also subpopulations, which can be useful in forensic and anthropological studies.

Delineation of Mitochondrial DNA Variants From Exome Sequencing Data and Association of Haplogroups With Obesity in Kuwait

BACKGROUND/OBJECTIVES

Whole-exome sequencing is a valuable tool to determine genetic variations that are associated with rare and common health conditions. A limited number of studies demonstrated that mitochondrial DNA can be captured using whole-exome sequencing. Previous studies have suggested that mitochondrial DNA variants and haplogroup lineages are associated with obesity. Therefore, we investigated the role of mitochondrial variants and haplogroups contributing to the risk of obesity in Arabs in Kuwait using exome sequencing data.

SUBJECTS/METHODS

Indirect mitochondrial genomes were extracted from exome sequencing data from 288 unrelated native Arab individuals from Kuwait. The cohort was divided into obese [body mass index (BMI) ≥ 30 kg/m2] and non-obese (BMI < 30 kg/m2) groups. Mitochondrial variants were identified, and haplogroups were classified and compared with other sequencing technologies. Statistical analysis was performed to determine associations and identify mitochondrial variants and haplogroups affecting obesity.

RESULTS

Haplogroup R showed a protective effect on obesity [odds ratio (OR) = 0.311; P = 0.006], whereas haplogroup L individuals were at high risk of obesity (OR = 2.285; P = 0.046). Significant differences in mitochondrial variants between the obese and non-obese groups were mainly haplogroup-defining mutations and were involved in processes in energy generation. The majority of mitochondrial variants and haplogroups extracted from exome were in agreement with technical replica from Sanger and whole-genome sequencing.

CONCLUSIONS

This is the first to utilize whole-exome data to extract entire mitochondrial haplogroups to study its association with obesity in an Arab population.

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

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

Platinum-Quality Mitogenome Haplotypes from United States Populations

A total of 1327 platinum-quality mitochondrial DNA haplotypes from United States (U.S.) populations were generated using a robust, semi-automated next-generation sequencing (NGS) workflow with rigorous quality control (QC). The laboratory workflow involved long-range PCR to minimize the co-amplification of nuclear mitochondrial DNA segments (NUMTs), PCR-free library preparation to reduce amplification bias, and high-coverage Illumina MiSeq sequencing to produce an average per-sample read depth of 1000 × for low-frequency (5%) variant detection. Point heteroplasmies below 10% frequency were confirmed through replicate amplification, and length heteroplasmy was quantitatively assessed using a custom read count analysis tool. Data analysis involved a redundant, dual-analyst review to minimize errors in haplotype reporting with additional QC checks performed by EMPOP. Applying these methods, eight sample sets were processed from five U.S. metapopulations (African American, Caucasian, Hispanic, Asian American, and Native American) corresponding to self-reported identity at the time of sample collection. Population analyses (e.g., haplotype frequencies, random match probabilities, and genetic distance estimates) were performed to evaluate the eight datasets, with over 95% of haplotypes unique per dataset. The platinum-quality mitogenome haplotypes presented in this study will enable forensic statistical calculations and thereby support the usage of mitogenome sequencing in forensic laboratories.

Comprehensive analysis of mitochondrial and nuclear DNA variations in patients affected by hemoglobinopathies: A pilot study

The hemoglobin disorders are the most common single gene disorders in the world. Previous studies have suggested that they are deeply geographically structured and a variety of genetic determinants influences different clinical phenotypes between patients inheriting identical β-globin gene mutations. In order to get new insights into the heterogeneity of hemoglobin disorders, we investigated the molecular variations on nuclear genes (i.e. HBB, HBG2, BCL11A, HBS1L and MYB) and mitochondrial DNA control region. This pilot study was carried out on 53 patients belonging to different continents and molecularly classified in 4 subgroup: β-thalassemia (β+/β+, β0/β0 and β+/β0)(15), sickle cell disease (HbS/HbS)(20), sickle cell/β-thalassemia (HbS/β+ or HBS/β0)(10), and non-thalassemic compound heterozygous (HbS/HbC, HbO-Arab/HbC)(8). This comprehensive phylogenetic analysis provided a clear separation between African and European patients either in nuclear or mitochondrial variations. Notably, informing on the phylogeographic structure of affected individuals, this accurate genetic stratification, could help to optimize the diagnostic algorithm for patients with uncertain or unknown origin.

Skim-Sequencing Based Genotyping Reveals Genetic Divergence of the Wild and Domesticated Population of Black Tiger Shrimp (Penaeus monodon) in the Indo-Pacific Region

The domestication of a wild-caught aquatic animal is an evolutionary process, which results in genetic discrimination at the genomic level in response to strong artificial selection. Although black tiger shrimp (Penaeus monodon) is one of the most commercially important aquaculture species, a systematic assessment of genetic divergence and structure of wild-caught and domesticated broodstock populations of the species is yet to be documented. Therefore, we used skim sequencing (SkimSeq) based genotyping approach to investigate the genetic structure of 50 broodstock individuals of P. monodon species, collected from five sampling sites (n = 10 in each site) across their distribution in Indo-Pacific regions. The wild-caught P. monodon broodstock population were collected from Malaysia (MS) and Japan (MJ), while domesticated broodstock populations were collected from Madagascar (MMD), Hawaii, HI, USA (MMO), and Thailand (MT). After various filtering process, a total of 194,259 single nucleotide polymorphism (SNP) loci were identified, in which 4983 SNP loci were identified as putatively adaptive by the pcadapt approach. In both datasets, pairwise F_ST estimates high genetic divergence between wild and domesticated broodstock populations. Consistently, different spatial clustering analyses in both datasets categorized divergent genetic structure into two clusters: (1) wild-caught populations (MS and MJ), and (2) domesticated populations (MMD, MMO and MT). Among 4983 putatively adaptive SNP loci, only 50 loci were observed to be in the coding region. The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses suggested that non-synonymous mutated genes might be associated with the energy production, metabolic functions, respiration regulation and developmental rates, which likely act to promote adaptation to the strong artificial selection during the domestication process. This study has demonstrated the applicability of SkimSeq in a highly duplicated genome of P. monodon specifically, across a range of genetic backgrounds and geographical distributions, and would be useful for future genetic improvement program of this species in aquaculture.

Genetic Consequences of the Transatlantic Slave Trade in the Americas

According to historical records of transatlantic slavery, traders forcibly deported an estimated 12.5 million people from ports along the Atlantic coastline of Africa between the 16th and 19th centuries, with global impacts reaching to the present day, more than a century and a half after slavery's abolition. Such records have fueled a broad understanding of the forced migration from Africa to the Americas yet remain underexplored in concert with genetic data. Here, we analyzed genotype array data from 50,281 research participants, which-combined with historical shipping documents-illustrate that the current genetic landscape of the Americas is largely concordant with expectations derived from documentation of slave voyages. For instance, genetic connections between people in slave trading regions of Africa and disembarkation regions of the Americas generally mirror the proportion of individuals forcibly moved between those regions. While some discordances can be explained by additional records of deportations within the Americas, other discordances yield insights into variable survival rates and timing of arrival of enslaved people from specific regions of Africa. Furthermore, the greater contribution of African women to the gene pool compared to African men varies across the Americas, consistent with literature documenting regional differences in slavery practices. This investigation of the transatlantic slave trade, which is broad in scope in terms of both datasets and analyses, establishes genetic links between individuals in the Americas and populations across Atlantic Africa, yielding a more comprehensive understanding of the African roots of peoples of the Americas.

Mitochondrial DNA haplogrouping to assist with the identification of unknown service members from the World War II Battle of Tarawa

The World War II Battle of Tarawa, 1943, was a devastating conflict that resulted in losses of more than 1100 American and 4690 Japanese troops. The United States government aims to identify and repatriate the remains of all missing American service members through the Defense Prisoner of War/Missing in Action (POW/MIA) Accounting Agency (DPAA) and its partners such as the Armed Forces Medical Examiner System's Armed Forces DNA Identification Laboratory (AFMES-AFDIL). Remains associated with the Battle of Tarawa have been recovered from field excavations conducted by History Flight, a DPAA strategic partner, as well as from the National Memorial Cemetery of the Pacific (NMCP) in Hawaii where unknowns have been disinterred for identification. DNA testing at the AFMES-AFDIL has produced mitochondrial DNA (mtDNA) sequences from 1027 case samples to date. Haplogroup assignments indicate that more than one third (36.2 %) of field-collected samples are likely of Asian maternal ancestry. Therefore the field collections from the Tarawa battlefield comprise the remains of American service members but also those of foreign nationals from Asia. The mtDNA of the NMCP unknowns is similar in ancestry proportion to the family reference sample distribution. The DPAA uses the ancestry information gleaned from mtDNA sequence data in conjunction with anthropological evidence to make foreign national determinations. In this way, mtDNA haplogrouping is used to sort the commingled and fragmentary remains recovered from Tarawa between Americans and foreign nationals, which are then repatriated to their country of origin.

Population genetics of two chromatic morphs of the Chagas disease vector Rhodnius pallescens Barber, 1932 in Panamá

Rhodnius pallescens is the principal vector of Chagas disease in Panama. Recently a dark chromatic morph has been discovered in the highlands of Veraguas Province. Limited genetic studies have been conducted with regards to the population structure and dispersal potential of Triatominae vectors, particularly in R. pallescens. Next generation sequencing methods such as RADseq and complete mitochondrial DNA (mtDNA) genome sequencing have great potential for examining vector biology across space and time. Here we utilize a RADseq method (3RAD), along with complete mtDNA sequencing, to examine the population structure of the two chromatic morpho types of R. pallescens in Panama. We sequenced 105 R. pallescens samples from five localities in Panama. We generated a 2216 SNP dataset and 6 complete mtDNA genomes. RADseq showed significant differentiation among the five localities (F_CT = 0.695; P = .004), but most of this was between localities with the dark vs. light chromatic morphs (Veraguas vs. Panama Oeste). The mtDNA genomes showed a 97-98% similarity between dark and light chromatic morphs across all genes and a 502 bp insert in light morphs. Thus, both the RADseq and mtDNA data showed highly differentiated clades with essentially no gene flow between the dark and light chromatic morphs from Veraguas and central Panama respectively. We discuss the growing evidence showing clear distinctions between these two morpho types with the possibility that these are separate species, an area of research that requires further investigation. Finally, we discuss the cost-effectiveness of 3RAD which is a third of the cost compared to other RADseq methods used recently in Chagas disease vector research.