Next Article in Journal
Association between NTRK2 Polymorphisms, Hippocampal Volumes and Treatment Resistance in Major Depressive Disorder
Next Article in Special Issue
Genomic Signatures of Positive Selection in Human Populations of the OXT, OXTR, AVP, AVPR1A and AVR1B Gene Variants Related to the Regulation of Psychoemotional Response
Previous Article in Journal
A High-Quality Reference Genome Assembly of Prinsepia uniflora (Rosaceae)
Previous Article in Special Issue
High Mitochondrial Haplotype Diversity Found in Three Pre-Hispanic Groups from Colombia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unraveling the Genetic Threads of History: mtDNA HVS-I Analysis Reveals the Ancient Past of the Aburra Valley

by
Daniel Uricoechea Patiño
1,
Andrew Collins
2,
Oscar Julián Romero García
3,
Gustavo Santos Vecino
4,
Pablo Aristizábal Espinosa
5,
Jaime Eduardo Bernal Villegas
6,
Escilda Benavides Benitez
6,
Saray Vergara Muñoz
6 and
Ignacio Briceño Balcázar
1,*
1
Doctoral Program in Biosciences, Human Genetics Group, Faculty of Medicine, University of La Sabana, Chía 250001, Colombia
2
Human Genetics & Genomic Medicine, Faculty of Medicine, University of Southampton, Southampton SO16 6YD, UK
3
Genetics Group, National Institute of Legal Medicine and Forensic Sciences, Bogotá 110311, Colombia
4
Department of Anthropology, Faculty of Social and Human Science, Universidad de Antioquia, Medellín 050010, Colombia
5
School for Advanced Studies in the Social Sciences, Americanists Society of Paris, 75006 Paris, France
6
Faculty of Medicine, University of Sinú, Cartagena de Indias 130011, Colombia
*
Author to whom correspondence should be addressed.
Genes 2023, 14(11), 2036; https://doi.org/10.3390/genes14112036
Submission received: 9 August 2023 / Revised: 18 September 2023 / Accepted: 30 October 2023 / Published: 2 November 2023
(This article belongs to the Special Issue Population Structure and Human Genetic Diversity)

Abstract

:
This article presents a comprehensive genetic study focused on pre-Hispanic individuals who inhabited the Aburrá Valley in Antioquia, Colombia, between the tenth and seventeenth centuries AD. Employing a genetic approach, the study analyzed maternal lineages using DNA samples obtained from skeletal remains. The results illuminate a remarkable degree of biological diversity within these populations and provide insights into their genetic connections with other ancient and indigenous groups across the American continent. The findings strongly support the widely accepted hypothesis that the migration of the first American settlers occurred through Beringia, a land bridge connecting Siberia to North America during the last Ice Age. Subsequently, these early settlers journeyed southward, crossing the North American ice cap. Of particular note, the study unveils the presence of ancestral lineages from Asian populations, which played a pivotal role in populating the Americas. The implications of these results extend beyond delineating migratory routes and settlement patterns of ancient populations. They also enrich our understanding of the genetic diversity inherent in indigenous populations of the region. By revealing the genetic heritage of pre-Hispanic individuals from the Aburrá Valley, this study offers valuable insights into the history of human migration and settlement in the Americas. Furthermore, it enhances our comprehension of the intricate genetic tapestry that characterizes indigenous communities in the area.

1. Introduction

The Antioquia region has played a pivotal role in the pre-Hispanic history of South America, serving as a gateway for the initial human migrations to the continent and as a hub for the development of a sophisticated culture deeply rooted in the utilization of natural resources and advanced technologies [1]. Colombia boasts a multitude of archaeological sites that have enabled the tracking of the earliest human arrivals, dating back to approximately sixteen thousand four hundred (16,400) ± four hundred twenty (420) A.P. [2], as determined through radiocarbon dating of artifacts linked to ancient human activities. Among these sites, Pubenza, nestled in the lowlands of the middle-valley of the Magdalena River [3], and El Jordán in the municipality of Roncesvalles (Tolima), stand out as some of the oldest, alongside El Totumo, El Abra, Tibitó, Tequendama, and Nemocón. All are strategically situated in the Sabana de Bogotá within the Eastern Cordillera [4] (Table 1).
The earliest radiocarbon-dated evidence within the Antioquia Department is prominently concentrated in the Medellín/Porce River basin and the Medellín/Porce River high basin (Table 1) [5]. Early agricultural societies in the region exhibit distinct burial practices, including simple pit burials and cremations conducted in ceramic vessels. These burials, in some instances, feature accompanying vessels used as lids or funerary ornaments [6]. Two distinct funerary patterns have emerged within Colombia.
The first pattern corresponds to the Marrón Inciso ceramic style, spanning from the first century BC to the eighth century AD. The second pattern corresponds to the Ferrería ceramic style, ranging from the fourth century BC to the fourth century AD. The Marrón Inciso style extends throughout the Cauca mountain basin, while the Ferrería style is prominent in the central mountain range of Antioquia. Cerro El Volador, situated within the city of Medellín, has yielded Marrón Inciso-style burials, including an urn containing cremated human remains and a necklace bead resembling a praying mantis, linking ceramics to gold craftsmanship. The second funerary pattern, originating in Antioquia in the fourteenth century AD, is associated with the late ceramic style, spanning from the tenth century AD to the time of the Conquest. This pattern is characterized by burials in intricate and spacious funerary structures known as “pit tombs with a lateral chamber”, commonly found atop hills in cemeteries [7].
Cerro El Volador, for example, revealed twelve tombs featuring rectangular access shafts and expansive chambers of conical and elliptical configurations. These tombs, dated between the fifteenth and seventeenth centuries AD, and contained cremated human remains, indigenous pottery fragments, skeletal remains of cattle and horses, along with European artifacts. The architectural design of these burial chambers mirrors indigenous dwellings or huts. At La Colinita’s summit, three “side chamber well tombs” unveiled cremated and non-cremated human remains, spindle ruffles, gold nose-rings, tumbaga, vessels, and ceramic fragments. Another burial site atop Las Flores contained the skeletal remains of an adult individual, devoid of grave goods but adorned with engravings representing a hut. The “pit tombs with a lateral chamber” in the Aburrá Valley are intricate representations of indigenous dwellings, holding profound significance both in construction and cosmology [7] These are hills with privileged panoramic views, where accumulations of tombs or necropolises were discovered. It should be clarified that the tombs investigated in the three hills of the Aburrá Valley, where the city of Medellín is located, belong to this latter period.
Molecular biology, leveraging DNA analysis, has played a pivotal role in shedding light on hypotheses concerning the origins of human populations. DNA, the repository of hereditary information and essential protein-coding data, was first structurally elucidated in 1953 by Watson, Crick, and Franklin. Subsequently, Sanger’s sequencing method paved the way for reading and decoding genetic information. Recent advancements in data storage, analytical technologies, and informatics have catalyzed breakthroughs across multiple disciplines, including archaeology. Despite the inevitable degradation of DNA post-mortem, it can adhere to hydroxyapatite in bones, maintaining sufficient quality for reconstructing ancient genomes.
Mitochondrial DNA (mtDNA) analysis, particularly in tracing maternal lineages, has emerged as a powerful tool for unraveling population history and genetic ancestry. The exclusive maternal transmission of mtDNA enables the tracking of maternal lineages across generations, facilitating the reconstruction of ancient population movements and migrations.
The high copy-number of mtDNA within cells empowers the recovery and examination of genetic information from ancient remains, thereby aiding our comprehension of genetic diversity and population dynamics over time. Maternal lineages, a focal point of mtDNA analysis, are identified through the detection of characteristic mutations that serve as markers for distinct population groups. These findings have significantly contributed to the development of theories regarding the settlement of pre-Hispanic communities in Colombia.
Ancient DNA studies have sought to unravel the genetic structure of pre-Colombian populations and establish genetic links between historical and contemporary population groups. These investigations are instrumental in supporting models related to population replacements, migrations, and genetic admixtures. As of now, genetic ancestry data in the Aburrá Valley are scarce.
To gain a more comprehensive understanding of population dynamics, migration patterns, and genetic relationships within the region, future studies should consider incorporating complete nuclear and mitochondrial genome analysis. This advanced approach will enable a more detailed and precise assessment of the genetic history of the pre-Hispanic inhabitants in the Aburrá Valley; thereby, enriching our comprehension of the complex tapestry of human settlement and evolution in this pivotal region.

2. Materials and Methods

2.1. Sample Collection

This study was centered on the analysis of mitochondrial DNA HVS-I extracted from ten individuals originating from the Aburrá Valley in Antioquia. These individuals were associated with pit tombs characterized by lateral chambers. Gender determination was based on the anthropological examination of the skeletal remains.
We conducted research at three distinct archaeological sites: Cerro del Volador, which yielded one individual (bone); Alto de las Flores in Envigado, where we analyzed one individual (bone) in Medellín; and La Colinita, where we investigated eight individuals (molars) in Medellín (Figure 1). The samples selected for analysis date back to the late pre-Hispanic period and the contact period, specifically spanning from the eleventh to the eighteenth centuries AD [7].

2.2. Ancient DNA Extraction

The genetic analysis of individuals from bone remains was conducted at the Human Genetics Laboratory within the specialized area dedicated to archaic DNA procedures at the University of Sabana. Stringent authenticity criteria, designed to prevent contamination of ancient DNA as proposed by Cooper, Poinar, and Gilbert [8], were strictly adhered to during the entire process. Mitochondrial DNA sequences from the HVS-I region were detected through the polymerase chain reaction (PCR) amplification method. Subsequently, these sequences were aligned with the revised Cambridge Reference Sequence (rCRS) for detailed analysis [9].
To ensure the accuracy and reliability of the findings, the procedures were independently replicated by two researchers on separate occasions. Furthermore, a subset of the samples was sent to the Department of Genetic Epidemiology and Bioinformatics at the University of Southampton in England.

2.3. PCR Amplification and Sequencing

To determine the mitochondrial haplotypes of the individuals in this study, we focused on a 388 bp fragment of HVS-I, specifically nucleotide positions 16024–16408. For amplification, we employed four pairs of overlapping PCR primers that were specifically developed for pre-Hispanic American populations [10]: H15995: GCT AAG ATT CTA ATT TAA ACT ATT CT (26 bp); L16174: GGA TTG GGT TTT TAT GTA CTA C (22 bp); H16105: TGC CAG CCA CCA TGA ATA TTG TAC (24 bp); L16256: GCT TTG GAG TTG CAG TTG ATG TGT (24 bp); H16194: ATG CTT ACA AGC AAG TAC AGC AA (23 bp); L16360: GAG AAG GGA TTT GAC TGT AAT GTG (24 bp); H16261: CCT CAC CCA CTA GGA TAC CAA CAA (24 bp); L16429: CCT CAC CCA CTA GGA TAC CAA CAA (21 bp). For the amplification of the HVS-I in human bone and molar remains containing ancient DNA, a standardized process was employed. This process involved a final volume of 10 µL and PCR amplification utilizing a hot start Polymerase and dNTPs. The resulting PCR product was verified using a 2% agarose gel, and the most suitable bands were selected for Sanger sequencing. Stringent alignment and analysis was performed using the QIAGEN CLC Genomics Workbench 20.0.4 program. Variants were independently and repeatedly identified by three investigators; each was assigned a GenBank accession number within the range BankIt 2736055: OR478541–OR478622. Haplogroups and haplotypes were determined using the HaploGrep2 program. Variations with a quality score exceeding 60% present in at least three duplicate samples were documented to ensure the reliability of our findings.

2.4. Statistical Analyses

To conduct individual analyses, we employed specific haplotypes and constructed phylogenetic trees using unique polymorphisms for each haplogroup. We determined common evolutionary relationships through several approaches, including neighbor-joining algorithms, UPGMA, direct networks using the NETWORK 10 program, and the HaploGrep 2 program [11].
Our study also encompassed intra- and inter-population analyses, performed with the Arlequin 3.5 program. We calculated various indices to discern population trends, and mean pairwise FST genetic distances were computed between populations based on their mitochondrial haplotype frequencies.
In this analysis, we considered twenty present-day populations from Colombia and seventy-seven pre-Hispanic American populations. These were categorized into eight groups based on previous research. We obtained distance matrices and visualized them using the MDS ALSCAL algorithm within the SPSS program. This method provides a balanced approach, addressing concerns regarding false positives and statistical power, as opposed to relying solely on phylogenetic trees [12,13,14,15].

3. Results

In our study, we identified a total of ten individuals, collectively representing the four founding haplogroups: A (60%), B (10%), C (20%), and D (10%) (Table 2). These haplotypes were defined by seventeen-point polymorphic sites, revealing thirty-one T↔C and thirteen A↔G transitions. It is worth noting that each individual possessed a unique haplotype, indicating the absence of a direct matrilineal relationship among the individuals under analysis.
To classify the mitochondrial haplogroups, we followed the Phylotree hierarchy [16], which begins with the common ancestor R, followed by N, M, and A. Among the individuals with the C16290T polymorphism, six were categorized into haplogroup A, with four of them belonging to the A2 subgroup due to distinct variants. Individual Abu10 was assigned to haplogroup B2s but also exhibited a variant characteristic of B4. Individuals Abu7 and Abu9 were designated as haplogroup C and shared specific variants, while the individual with the C16223T polymorphism was placed in haplogroup D, with individual Abu6 specifically categorized as D4j1a based on another unique variant.
The phylogenetic tree (Figure 2) illustrates the presence of various haplogroups and their corresponding samples. In the tree, we observe samples classified under haplogroup A, including A2, A2af1b1b, A, A2i, A2+(64)+ @16111, and A2ah. Additionally, there are samples classified under haplogroup C, namely C and C1c4. The tree also shows a sample categorized under haplogroup D, specifically D4j1a. Lastly, there is a sample denoted as B2s under haplogroup B2. These findings provide valuable insights into the genetic diversity and relationships among the studied population.
Median-joining (MJ) analysis was employed to explore maternal relationships between individuals in family graves discovered in the Aburrá Valley (Figure 3). The analysis did not reveal direct matrilineal relationships; however, it did show an ancestral connection between an Abu1 node and individual Abu9.
Our study’s results revealed the presence of ten individuals from the Aburrá Valley, representing four founding haplogroups (A, B, C, and D). Each individual exhibited a unique haplotype, indicating the absence of direct matrilineal relationships among them. Haplogroup assignments were made based on specific variants and polymorphisms identified in their sequences.
The haplotype network offers valuable insights into the genealogical relationships among haplotypes. It allows us to assess the degree of divergence between different lineages and estimate the distance to their most recent common ancestor. It is important to note, though, that the haplotype network does not directly establish maternal relationships between specific individuals. Instead, it primarily serves as a tool for understanding the genealogical connections and ancestral origins of different haplotypes. Its relevance lies in elucidating the historical relationships and evolutionary patterns within a population rather than establishing direct relationships between specific individuals.
Furthermore, the results of population diversity in the Aburrá Valley are presented in Table 3, with ‘n’ representing ten individuals and ‘K’ denoting ten sequences analyzed. The nucleotide diversity (π) for the study population is calculated as 0.012. Analysis of the Tajima index (-D) during pre-Hispanic periods indicates population expansion, with a result of −0.72937 and a p-value greater than 0.05; which suggests no significant differences between the analyzed samples [14].
In line with the research conducted by Raff [15], a multiscale analysis comparing seventy-eight ancient populations across the Americas with twenty contemporary indigenous Colombian populations highlights the highest genetic similarity in the Aburráes population (Figure 4).
This similarity is particularly striking with pre-Hispanic populations in Mesoamerica, as well as in the Arctic and Sub-Arctic regions. The study reveals the presence of haplogroups A, C, and D in both current and ancient populations throughout the Americas. Notably, haplogroup B is absent in the Caribbean, while haplogroup X is reported to be in low proportions in the northern part of the continent. Haplogroup M, on the other hand, is exclusively identified in an ancient population in the northwestern United States and not in current populations (Table 4).
The dataset analyzed includes archaeological sites categorized into various geographic zones, namely the Arctic/Subarctic, northwestern United States, Great Basin/southwestern United States, northeastern United States, Mesoamerica, Caribbean, and South America. This regional categorization, as presented in Table 4, is primarily based on the work of Raff [15]. However, it has been refined and expanded upon through recent research, with a particular focus on the inclusion of Colombia. The geographic distribution of these analyzed zones across the Americas is visually represented in Figure 5. Which employs a color scheme for easy identification of geographical locations and highlights the diverse array of regions covered in our study.
For more comprehensive information, including the haplogroup frequencies used for the multidimensional scaling (MDS) analysis, please consult Table 4. The specific matrix of FST values can be found in Table S1.
The distribution of mitochondrial haplogroups among the individuals revealed several patterns. Five individuals with the C16290T polymorphism were assigned to haplogroup A, with four of them specifically identified as A2. Notably, individual Abu10 was assigned to haplogroup B2s but also exhibited a variant typical of B4. Meanwhile, Individuals Abu4 and Abu7 were assigned to haplogroup C and shared specific variants. On the other hand, the individual with the C16223T polymorphism was placed in haplogroup D, with individual Abu6 classified as D4j1a, based on another specific variant (Table 5). These findings provide insights into the diversity of mitochondrial haplogroups in the studied population.

4. Discussion

The findings of this study offer valuable genetic insights into the population history of pre-Hispanic inhabitants in the Aburrá Valley region. This research primarily focused on the analysis of mitochondrial DNA (mtDNA), which has provided significant information about migratory patterns, interactions, and genetic relationships among ancient American populations.
Studying ancient populations is a multidisciplinary endeavor that combines genetics and archaeology to unravel the complex tapestry of human history. Haplogroups, which are genetic markers passed down through generations, serve as invaluable tools for understanding ancient population dynamics. In this research, we present a comprehensive table featuring archaeological sites across distinct geographic zones, each holding a wealth of genetic information in the form of haplogroups A, B, C, D, X, M, and H.
Variations in haplogroup distribution are evident across these geographic zones, indicating that different regions were inhabited by populations with unique genetic profiles. By examining these archaeological periods, we can observe how haplogroup distributions evolved over millennia; shedding light on genetic adaptations and the ebb and flow of ancient populations. Comparative analysis of haplogroups allows us to reconstruct population movements, migrations, and genetic relationships, providing insights into the interactions between ancient human groups. However, it is crucial to consider sample sizes, as they impact the statistical robustness of genetic findings. Larger sample sizes lead to more reliable insights.
Regarding genetic diversity, the Aburráes population in Colombia exhibited a genetic diversity level of 0.644. When compared to various populations from different geographic regions and archaeological periods, this falls within the range of genetic diversity observed across these populations. Some populations displayed higher genetic-diversity values, while others had lower values. Thus, reflecting the intricate tapestry of human genetic variation across different time periods and geographical locations. It is important to note that, in general, a larger number of samples provide a more accurate reflection of the true diversity of the population or community being studied. Therefore, to obtain reliable and representative results, it is recommended to use a significant number of samples in the diversity analysis.
It is important to emphasize that the haplotype network does not directly establish maternal relationships between specific individuals. Instead, its primary function is to elucidate genealogical connections and ancestral origins among various haplotypes. Its significance lies in unveiling historical relationships and evolutionary patterns within a population rather than confirming direct relationships between specific individuals.
Regarding the Colinita tomb, it’s worth noting that the eight adult individuals found within it are not from the same family. This information holds significant value for archaeologists and anthropologists as it provides insights into burial rituals and helps determine whether the tombs are related to a single family or not.
The meticulous examination of various haplogroups, including A2, B2, C1, and D4j1a, has provided valuable insights into the intricate demographic processes that shaped the region. Notably, Haplogroup A2, renowned for its prevalence in the high latitudes of the Americas and its historical ties to Siberian and Eskimo tribes, has emerged as a pivotal contributor to the initial waves of migration into the continent [16,74]. Moreover, the identification of three distinct sub-haplogroups within A2 among pre-Hispanic Antioquia populations signifies local genetic diversification, underscoring the region’s role as a genetic melting-pot.
This research delves deeply into the dynamics of these populations. The discovery of A2af1b1b within contemporary populations in Panama and Costa Rica carries profound implications, pointing to significant demographic events approximately ten thousand years ago. While this specific haplogroup is presently observed in only two individuals—one from Chiriqui, Panama, and the other from Costa Rica—it originated 4700 years before the present, with a time to most recent common ancestor (TMRCA) of 1000 years before the present [78].
Most intriguingly, the study identifies the presence of the A2af1b1b sub-haplogroup in an individual from the Valle de Aburra. This finding suggests that A2af1b1b has deep historical roots in the region, indicating its presence within ancient populations. Furthermore, the fact that this haplogroup is still present in contemporary populations emphasizes its enduring legacy [78]. This underscores the notion that A2af1b1b is not just a relic of the past but has persisted through generations and is still part of the genetic makeup of current populations in the region.
The enduring presence of Haplogroup A2ah for over five millennia across various regions offers a compelling narrative of historical population movements. This relatively rare lineage among Native Americans, discerned through comprehensive analyses of complete mitochondrial genomes, is characterized by the uncommon combination of tandem transitions, T16097C and A16098G. Which is rarely encountered in mitochondrial DNA databases [94].
Within the larger Haplogroup A2, A2ah emerges as a distinct branch within the Native American mitochondrial family. This research marks a significant milestone in the field of genetics. Notably, for the first time, researchers have unveiled a mitochondrial lineage in the Andean region that was previously believed to be exclusive to the South American lowlands. Specifically, they have identified the A2ah lineage in the Andean area, even within a pre-Hispanic individual. The implications of this discovery are profound, raising intriguing questions about the potential gene flow between the Andean and lowland regions before the Spanish arrival [95].
Archaeological evidence indicates interactions and trade among diverse populations across various eco-regions during pre-Hispanic times. The presence of trade networks, such as those connected to the Quebrada de La Cueva ravine, which housed the Pukara de La Cueva, suggests extensive pre-Hispanic connectivity. Yet, whether these networks also facilitated gene-flow remains uncertain. This underscores the significance of ancient DNA studies in shedding light on this aspect of pre-Hispanic population dynamics [96].
If indeed gene-flow was a regular occurrence between highland and lowland pre-Hispanic populations, it becomes plausible that certain genetic variants from one region could be found in the other. The identification of an infrequent mitochondrial haplotype from the South American lowlands in the Andean area in this study lends support to this hypothesis. To further investigate this proposition, the researchers recommend expanding genetic analyses to encompass other individuals retrieved from Pukara de La Cueva at the population level. Additionally, isotopic analysis could offer insights into the region where the analyzed individual resided during their early years [97].
The research suggests an approximate origin dating back to 5200 years ago, with some within-clade variations. Geographically, this lineage appears to be primarily concentrated in Bolivia and its adjacent regions with a limited distribution [98,99,100].
The discovery of the A2ah mitochondrial lineage in the Andean region, as previously discussed, presents an intriguing parallel to the genetic findings in the Valle de Aburra. In both instances, researchers have uncovered the existence of uncommon mitochondrial haplogroups in regions that defied earlier expectations.
This shared discovery of sub-haplogroup A2ah among both ancient and contemporary Native American populations emphasizes its profound historical importance and its enduring legacy within the tapestry of mitochondrial DNA. The presence of this lineage in diverse temporal contexts underscores its resilience and persistence through time.
By delving deeper into the study of this unique lineage, researchers can unlock invaluable insights into the historical narratives, migration patterns, and genetic tapestry of Native American communities. This connection between seemingly distant regions challenges preconceived notions and beckons for further exploration into the intricate web of human genetic heritage across the Americas.
The genetic narrative of the A2i haplogroup, explored in Achilli et al.’s 2008 study within the Valle de Aburra [75], finds resonance in the genetic tapestry of Alaska. Recent studies involving mitogenomes from Alaska have brought to light a range of haplogroups, including A2, A2a, and A2i. The intriguing presence of A2i in both locales speaks to the remarkable genetic continuity and adaptation across distinct regions and epochs [85,86,87].
Within the complex genetic framework of lineage C1, a thorough network analysis has unveiled a compelling historical narrative. This analysis, which considered sequences from diverse regions such as the Amur River, Japan, and the Americas, lends support to the notion that the foundational haplotype for Native American C1 finds its roots in the Amur River region, echoing previous hypotheses proposed by Starikovskaya [86].
The C1c subclade, an exclusive genetic branch within haplogroup C1 exclusive to Native American populations, stands as a testament to the deep ancestral ties between indigenous peoples and the Americas. C1c plays a pivotal role in genetic and phylogenetic investigations, offering a unique perspective into the evolutionary history and migration patterns of specific groups. This subclade provides invaluable insights into the genetic diversity of Native American populations and their ancestral journeys [87].
The research also delved into broader colonization patterns. The presence of haplogroup B2 across diverse Native American populations adds weight to the idea of an early expansion from North to South America. Additionally, the widespread distribution of haplogroup C1 and its subclades suggests a unified migratory wave that played a crucial role in populating the Americas. The identification of variant 16298C in Abu4 within haplogroup C underscores the distinctive genetic legacy of pre-Hispanic Native American settlers.
Moreover, the presence of Haplogroup C1c, with its older coalescent age and evidence of intra-haplogroup variation, in both Mexican Americans and ancient Valle de Aburra populations hints at a captivating connection. This shared genetic heritage suggests that the genetic diversity within Haplogroup C1c, as well as its sub-branches, was firmly established among ancestral populations before their migration and expansion into the Americas. The exploration of this haplogroup’s distribution and variation across populations promises profound insights into historical movements and genetic relationships among these groups, ultimately unraveling the intricate history of Native American origins and migrations.
Within the realm of Native American haplogroups, B2 emerges as a distinct lineage exclusive to indigenous populations of the Americas. Setting it apart from the other major founding haplogroups (A2, C1, D1, and X2a), the basal B2 lineage lacks defining variants in the control region (CR). This unique feature presents a challenge in distinguishing B2 from closely related haplogroups, a critical task for distinguishing Native American from Asian haplotypes [91].
Compelling evidence suggests that the initial wave of people entering what is now Mexican territory carried with them the A2 and B2 haplogroups around 15–18-thousand years ago. The distribution of individuals with higher frequencies of B2 or C1 compared to A2 indicates a specific migration pattern. A2, B2, and C1 appear to have been primarily dispersed along a Pacific coast route, spanning approximately 800 km along the Neovolcanic axis region. Subsequently, they made their way southward along the Pacific coast route, eventually reaching the Yucatan Peninsula via the Gulf coast route.
It is noteworthy that the Neovolcanic axis region, characterized by the presence of significant rivers like Lerma and Balsas, acts as a natural barrier against adverse Pacific weather conditions. This geographical feature likely played a pivotal role in facilitating the spread of these haplogroups [86,91].
Sub-haplotype D4j1a stands out from the norm when it comes to substitutions observed in mitochondrial DNA clusters associated with the southern distribution of aboriginal Siberian peoples. These substitutions are typically linked to major mitochondrial DNA haplogroups prevalent in East and South Asia. The intriguing presence of sub-haplotype D4j1a among both aboriginal Siberian populations and the ancient inhabitants of the Valle de Aburra unveils a fascinating genetic connection spanning vast geographical distances. Importantly, this sub-haplotype diverges from the usual mitochondrial DNA clusters found in aboriginal Siberian populations, aligning instead with major mitochondrial DNA haplogroups prevalent in East and South Asia [87,88,89,90,91,92,93].
The genetic connection between Siberia and the Valle de Aburra holds the potential to illuminate ancient migration patterns and the interactions between diverse populations. To gain a deeper understanding of the historical context and the implications of this shared genetic marker, further genetic and archaeological investigations are imperative (Figure 5).
Understanding the genetic connection between Siberia and Valle de Aburra offers a glimpse into the intricate journey of ancient human migrations. It suggests that populations from distant regions might have come into contact, possibly through migratory movements, trade routes, or cultural exchanges.
To provide a meaningful context for this genetic link, it is crucial to highlight the fundamental role of archaeological investigations, which, in combination with genetics, offer a more accurate representation of ancient life. The genetic relationship between the Valle de Aburra and Siberia, along with its proximity to Mesoamerica, suggests the diffusion of genetic markers across populations. Mesoamerica’s genetic diversity, a result of its complex history, can provide insights into ancient cultural exchanges and trade routes. A multidisciplinary approach, encompassing genetics and archaeological excavations, is essential for a comprehensive understanding of the history and implications of this shared genetic marker.
The revelation of a genetic connection weaving through Siberia, Mesoamerica, and Valle de Aburra is nothing short of enthralling. This revelation acts as a beckoning call, compelling us to embark on an expedition of profound historical significance. Expanding our investigation into Mesoamerica opens doors to unraveling ancient ties, migratory routes, and the intricate web of interactions that wove together diverse populations. In the realm of collaborative genetic and archaeological exploration, we stand on the precipice of unearthing the narratives of our forebears, shedding brilliant light upon the enigmatic pathways of our shared past.
Comparisons with other populations across the Americas revealed that the Aburráes population showed the highest genetic similarity with pre-Hispanic populations in America. The results presented in our study contribute to the understanding of the genetic diversity and ancestral connections of the Aburráes population. These findings have implications for the population history and migration patterns in the Americas.
Furthermore, the study sheds light on population relationships across regions. The genetic affinity between the Aburráes and Mesoamerican and Arctic–sub-Arctic populations showcase intricate gene flow and migratory dynamics. These findings align with the hypothesis of a coastal migration route and potential interactions between the indigenous Chibchas population and the Mayan population.
The primary limitation of this study arises from its exclusive focus on the HVS-I of ancient mtDNA. This choice was necessitated by the inherent challenges in working with degraded and ancient DNA samples. While the results obtained are indeed intriguing and informative, it is imperative to recognize and acknowledge this limitation.
Likewise, it is advisable to approach the study’s conclusions with a degree of caution. This caution is particularly warranted due to the relatively small sample size, consisting of only 10 individuals, and the restriction of the analysis to a single mitochondrial region (HVS-I).
Despite these challenges, it is noteworthy that the study was successful in elucidating important haplogroups. These findings make a significant contribution to our understanding of the genetic history of the ancient population under investigation.

5. Conclusions

In conclusion, the present study has provided valuable insights into the genetic diversity and population history of pre-Hispanic inhabitants in the Aburrá Valley region. The identification of various mitochondrial haplogroups, particularly A2 and its sub-lineages (A2af1b1b, A2ah, and A2i), Haplogroup B2, Haplogroup C1, and Haplogroup D4j1a, has shed light on the complex migratory patterns and interactions of ancient populations in the Americas.
One of the most striking findings of this research is the presence of distinct sub-haplogroups within the larger haplogroups, such as A2af1b1b, A2ah, and A2i. These sub-haplogroups not only reveal local genetic diversification but also highlight the enduring legacy of certain genetic lineages in the region.
While the genetic data presented in this study contribute valuable insights, it represents only a starting point in the investigation of the region’s population history. The pursuit of more extensive and refined research will undoubtedly bring us closer to a deeper appreciation of the rich and complex human history that unfolded in the Aburrá Valley throughout the ages.
Ongoing research efforts are vital to fully unravel the region’s population history. By continuously exploring and integrating diverse data sources, we can achieve a more nuanced understanding of the region’s past, its connections to broader human migrations in the Americas, and the factors shaping its distinctive genetic landscape.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes14112036/s1, Table S1: Genetic Distance Matrix Using Frequency Data from Seven Prehispanic Populations.

Author Contributions

Conceptualization, D.U.P., I.B.B. and J.E.B.V.; data curation, D.U.P. and A.C.; formal analysis, D.U.P. and O.J.R.G.; funding acquisition, I.B.B., J.E.B.V., E.B.B. and S.V.M.; investigation, D.U.P. and O.J.R.G.; methodology, D.U.P. and A.C.; project administration, I.B.B., A.C. and J.E.B.V.; resources, I.B.B., G.S.V. and P.A.E.; software, D.U.P. and A.C.; supervision, I.B.B. and A.C.; validation, D.U.P. and O.J.R.G.; visualization, D.U.P.; writing—original draft, D.U.P., I.B.B., E.B.B., S.V.M. and J.E.B.V.; writing—review and editing, D.U.P., I.B.B., J.E.B.V. and A.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Universidad de La Sabana through the Convocatoria interna de 2013 (Internal Call of 2013) as per acta 809 de la comisión de asuntos generales, proyecto MED-168-2013. The APC (Article Processing Charge) was funded by the Universidad de La Sabana.

Institutional Review Board Statement

Ethic Committee Name: Comité de Ética en Investigación, Dirección de Investigación, Universidad de La Sabana, Approval Code: Acta No. 29 del 25 de mayo de 2012, Approval Date: 25 May 2012.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data supporting the reported results can be found in GenBank flat file: LOCUS OR478541 483 bp DNA linear PRI 17-OCT-2023 DEFINITION Homo sapiens isolate M1-P1_H15995_P1 D-loop, partial sequence; mitochondrial.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Aceituno Bocanegra, F.J.; Uriarte González, A. Conectando un territorio: Simulación de rutas de movilidad entre cazadores-recolectores y primeros cultivadores. El caso del Cauca medio (Macizo Volcánico 2019, Colombia). Trab. Prehist. 2019, 76, 219–235. [Google Scholar] [CrossRef]
  2. Correal, G.; Van der Hammen, T. Supervivencia de Mastodontes, Megaterios y presencia del Hombre en el Valle del Magdalena (Colombia) entre 6000 y 5000 AP. Rev. Acad. Colomb. De Cienc. 2003, 27, 159–164. [Google Scholar] [CrossRef]
  3. Díaz-Matallana, M.; Gómez, A.; Briceño, I.; Rodríguez, J.V. Genetic analysis of Paleo-Colombians from Nemocón, Cundinamarca provides insights on the early peopling of Northwestern South America. Rev. La Acad. Colomb. Cienc. Exactas Físicas Nat. 2016, 40, 461–483. [Google Scholar] [CrossRef]
  4. Delgado, M.; Rodríguez, F.; Kassadjikova, K.; Fehren-Schmitz, L. A paleogenetic perspective of the Sabana de Bogotá (Northern South America) population history over the Holocene (9000–550 cal BP). bioRxiv 2020. bioRxiv:2020.01.24.918425. [Google Scholar] [CrossRef]
  5. Gustavo, V.S.; Marín, C.A.M.; Salas, L.V.C. Alteration of tropical forest vegetation from the Pleistocene–Holocene transition and plant cultivation from the end of early Holocene through Middle Holocene in Northwest Colombia. Quat. Int. 2015, 363, 28–42. [Google Scholar]
  6. Aristizábal Espinosa, P. Los Aburráes. In Tras Los Rastros de Nuestros Ancestros, 1st ed.; Secretaría de Cultura Ciudadana, Alcaldía de Medellín: Medellín, Colombia, 2015. [Google Scholar]
  7. Santos Vecino, G. Las prácticas funerarias prehispánicas de la región central de Antioquia. In Catálogo Región Andina; Universidad de Antioquia: Medellín, Colombia, 2017; pp. 41–71. [Google Scholar]
  8. Cooper, A.; Poinar, H. Ancient DNA: Do it right or not at All. Science 2000, 289, 1139. [Google Scholar] [CrossRef]
  9. Anderson, S.; Bankier, A.; Barrell, B.; de Bruijn, M.; Coulson, A.; Drouin, J.; Eperon, I.; Nierlich, D.; Roe, B.; Sanger, F.; et al. Sequence and organization of the human mitochondrial genome. Nature 1981, 290, 457–465. [Google Scholar] [CrossRef] [PubMed]
  10. Fehren-Schmitz, L.; Warnberg, O.; Reindel, M.; Seidenberg, V.; Tomasto-Cagigao, E.; Isla-Cuadrado, J.; Hummel, S.; Herrmann, B. Diachronic investigations of mitochondrial and Y-chromosomal genetic markers in Pre-Columbian Andean highlanders from South Peru. Ann. Hum. Genet. 2011, 75, 266–283. [Google Scholar] [CrossRef]
  11. Excoffier, L.; Lischer, H. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 2010, 10, 564–567. [Google Scholar] [CrossRef]
  12. Phillip, M.; Briceño, I.; Gómez, A.; Devor, E.; Bernal, J.E.; Crawford, M. Biological relationship between central and South American Chibchan speaking populations: Evidence from mtDNA. Am. J. Phys. Anthropol. 2007, 133, 753–770. [Google Scholar] [CrossRef]
  13. Claudia, N.S.M.; Anderson, C.E.; Uricoechea, D.; Durán, C.; Briceño-Balcázar, I.; Villegas, J.B. Mitochondrial DNA analysis suggests a Chibchan migration into Colombia. Univ. Sci. 2015, 20, 261–278. [Google Scholar] [CrossRef]
  14. Casas-Vargas, A.; Romero, L.; Usaquén, W.; Zea, S.; Silva, M.; Briceño, I.; Gómez, A.; Rodríguez, J.V. Diversidad del ADN mitocondrial en restos óseos prehispánicos. Biomédica 2017, 37, 548–560. [Google Scholar] [CrossRef] [PubMed]
  15. Raff., A.J.; Bolnick, D.; Tackney, J.; O’Rourke, D. Ancient DNA perspectives on American colonization and population history. Am. J. Phys. Anthropol. 2011, 146, 503–514. [Google Scholar] [CrossRef] [PubMed]
  16. Starikovskaya, Y.B.; Sukernik, R.I.; Schurr, T.G.; Kogelnik, A.M.; Wallace, D.C. mtDNA diversity in Chukchi and Siberian Eskimos: Implications for the genetic history of Ancient Beringia and the peopling of the New World. Am. J. Hum. Genet. 1998, 63, 1473–1491. [Google Scholar] [CrossRef] [PubMed]
  17. Smith, S.E.; Hayes, M.G.; Cabana, G.S.; Huff, C.; Coltrain, J.B.; O’rourke, D.H. Inferring population continuity versus replacement with aDNA: A cautionary tale from the Aleutian Islands. Hum. Biol. 2009, 81, 407–426. [Google Scholar] [CrossRef] [PubMed]
  18. Raff, J.; Tackney, J.; O’Rourke, D.H. South from Alaska: A pilot aDNA study of genetic history on the Alaska Peninsula and the eastern Aleutians. Hum. Biol. 2010, 82, 677–693. [Google Scholar] [CrossRef]
  19. Malhi, R.S.; Kemp, B.M.; Eshleman, J.A.; Cybulski, J.; Smith, D.G.; Cousins, S.; Harry, H. Mitochondrial haplogroup M discovered in prehistoric North Americans. J. Archaeol. Sci. 2007, 34, 642–648. [Google Scholar] [CrossRef]
  20. Monsalve, M.V.; Stone, A.C.; Lewis, C.M.; Rempel, A.; Richards, M.; Straathof, D.; Devine, D.V. Brief communication: Molecular analysis of the Kwäday Dän Ts’ finchi ancient remains found in a glacier in Canada. Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol. 2002, 119, 288–291. [Google Scholar] [CrossRef] [PubMed]
  21. Hayes, M.G. Paleogenetic Assessments of Human Migration and Population Replacement in North American Arctic Prehistory; The University of Utah: Salt Lake City, UT, USA, 2002. [Google Scholar]
  22. Gilbert, M.T.P.; Kivisild, T.; Grønnow, B.; Andersen, P.K.; Metspalu, E.; Reidla, M.; Tamm, E.; Axelsson, E.; Campos, P.F.; Willerslev, E.; et al. Paleo-Eskimo mtDNA genome reveals matrilineal discontinuity in Greenland. Science 2008, 320, 1787–1789. [Google Scholar] [CrossRef]
  23. Malhi, R.S.; Mortensen, H.M.; Eshleman, J.A.; Kemp, B.M.; Lorenz, J.G.; Kaestle, F.A.; Johnson, J.R.; Gorodezky, C.; Smith, D.G. Native American mtDNA prehistory in the American southwest. Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol. 2003, 120, 108–124. [Google Scholar] [CrossRef]
  24. Kaestle, F.A.; Horsburgh, K.A. Ancient DNA in anthropology: Methods, applications, and ethics. Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol. 2002, 119, 92–130. [Google Scholar] [CrossRef] [PubMed]
  25. Stone, A.C.; Stoneking, M. mtDNA analysis of a prehistoric Oneota population: Implications for the peopling of the New World. Am. J. Hum. Genet. 1998, 62, 1153–1170. [Google Scholar] [CrossRef]
  26. Snow, M.H.; Durand, K.R.; Smith, D.G. Ancestral Puebloan mtDNA in context of the greater southwest. J. Archaeol. Sci. 2010, 37, 1635–1645. [Google Scholar] [CrossRef] [PubMed]
  27. Carlyle, S.W.; Parr, R.L.; Hayes, M.G.; O’Rourke, D.H. Context of maternal lineages in the Greater Southwest. Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol. 2000, 113, 85–101. [Google Scholar] [CrossRef]
  28. LeBlanc, S.A.; Cobb Kreisman, L.S.; Kemp, B.M.; Smiley, F.E.; Carlyle, S.W.; Dhody, A.N.; Benjamin, T. Quids and aprons: Ancient DNA from artifacts from the American Southwest. J. Field Archaeol. 2007, 32, 161–175. [Google Scholar] [CrossRef]
  29. Schultz Shook, B.A.; Smith, D.G. Using ancient mtDNA to reconstruct the population history of northeastern North America. Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol. 2008, 137, 14–29. [Google Scholar] [CrossRef] [PubMed]
  30. Bolnick, D.A.; Smith, D.G. Migration and social structure among the Hopewell: Evidence from ancient DNA. Am. Antiq. 2007, 72, 627–644. [Google Scholar] [CrossRef]
  31. Raff, J.A. An Ancient DNA Perspective on the Prehistory of the Lower Illinois Valley. Doctoral Dissertation, Indiana University, Bloomington, Indiana, 2008. [Google Scholar]
  32. Mills, L.A. Mitochondrial DNA Analysis of the Ohio Hopewell of the Hopewell Mound Group; The Ohio State University: Columbus, OH, USA, 2003. [Google Scholar]
  33. Kemp, B.M.; Reséndez, A.; Román Berrelleza, J.A.; Malhi, R.S.; Smith, D.G. An analysis of ancient Aztec mtDNA from Tlatelolco: Pre-Columbian relations and the spread of Uto-Aztecan. Biomol. Archaeol. Genet. Approaches Past 2005, 32, 22–42. [Google Scholar]
  34. González-Oliver, A.; Márquez-Morfín, L.; Jiménez, J.C.; Torre-Blanco, A. Founding amerindian mitochondrial DNA lineages in ancient Maya from Xcaret, Quintana Roo. Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol. 2001, 116, 230–235. [Google Scholar] [CrossRef]
  35. Merriwether, D.A.; Reed, D.M.; Ferrell, R.E. Ancient and Contemporary Mitochondrial DNA Variation in the Maya; Smithsonian Institute Press: Washington, DC, USA, 1997; pp. 208–217. [Google Scholar]
  36. Lalueza-Fox, C.; Calderon, F.L.; Calafell, F.; Morera, B.; Bertranpetit, J. MtDNA from extinct Tainos and the peopling of the Caribbean. Ann. Hum. Genet. 2001, 65, 137–151. [Google Scholar] [CrossRef]
  37. Lalueza-Fox, C.; Gilbert, M.T.P.; Martínez-Fuentes, A.J.; Calafell, F.; Bertranpetit, J. Mitochondrial DNA from pre-Columbian Ciboneys from Cuba and the prehistoric colonization of the Caribbean. Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol. 2003, 121, 97–108. [Google Scholar] [CrossRef] [PubMed]
  38. Torroni, A.; Sukernik, R.I.; Schurr, T.G.; Starikorskaya, Y.B.; Cabell, M.F.; Crawford, M.H.; Comuzzie, A.G.; Wallace, D.C. mtDNA variation of aboriginal Siberians reveals distinct genetic affinities with Native Americans. Am. J. Hum. Genet. 1993, 53, 591. [Google Scholar] [PubMed]
  39. Santos, M.; Ward, R.H.; Barrantes, R. mtDNA variation in the Chibcha Amerindian Huetar from Costa Rica. Hum. Biol. 1994, 66, 963–977. [Google Scholar] [PubMed]
  40. Kolman, C.J.; Bermingham, E.; Cooke, R.; Ward, R.H.; Arias, T.D.; Guionneau-Sinclair, F. Reduced mtDNA diversity in the Ngöbé Amerinds of Panamá. Genetics 1995, 140, 275–283. [Google Scholar] [CrossRef] [PubMed]
  41. Batista, O.; Kolman, C.J.; Bermingham, E. Mitochondrial DNA diversity in the Kuna Amerinds of Panama. Hum. Mol. Genet. 1995, 4, 921–929. [Google Scholar] [CrossRef] [PubMed]
  42. Casas-Vargas, A.; Gómez, A.; Briceño, I.; Díaz-Matallana, M.; Bernal, J.E.; Rodríguez, J.V. High genetic diversity on a sample of pre-Columbian bone remains from Guane territories in northwestern Colombia. Am. J. Phys. Anthropol. 2011, 146, 637–649. [Google Scholar] [CrossRef] [PubMed]
  43. Uricoechea Patiño, D.; Collins, A.; García, O.J.R.; Santos Vecino, G.; Cuenca, J.V.R.; Bernal, J.E.; Benavides Benítez, E.; Vergara Muñoz, S.; Briceño Balcázar, I. High Mitochondrial Haplotype Diversity Found in Three Pre-Hispanic Groups from Colombia. Genes 2023, 14, 1853. [Google Scholar] [CrossRef]
  44. Fernandez, C. La Arqueología Molecular Aplicada a la Solución de Problemas Prehistóricos: Análisis de ADN Mitocondrial en Momias y Restos Óseos Prehispánicos. Undergraduate Thesis, Universidad Nacional de Colombia, Bogotá, Colombia, 1999. [Google Scholar]
  45. Jara, N.P.; Díaz, M.; Villegas, V.; López de Mesa, C.; Torres, D.; Bernal, J.; Gómez, A.; Briceño, I. Aplication of authenticity criteria in mitochondrial studies on archaic bone remains from a prehispanic Muisca population. Colomb. Médica 2010, 41, 306–314. [Google Scholar]
  46. Sánchez, C. Secuenciación de ADN Mitocondrial a Partir de Fragmentos Óseos Prehispánicos Hallados en el Sector de Candelaria La Nueva en Bogotá; Pontificia Universidad Javeriana: Bogotá, Colombia, 2007. [Google Scholar]
  47. Noguera-Santamaría, M.C.; Rivera-Sandoval, J.; Martín, J.G.; Briceño-Balcázar, I.; Gómez-Gutiérrez, A. Análisis genético de restos humanos precolombinos del Bajo Magdalena sugiere una ruta migratoria y continuidad genética matrilineal en el norte de Suramérica. Rev. La Acad. Colomb. Cienc. Exactas Físicas Nat. 2020, 44, 704–715. [Google Scholar] [CrossRef]
  48. Romero Murillo, L.M.; Silva Montaña, V.M.; Zea Montoya, S.; Usaquén Martínez, W.; Briceño Balcázar, I.; Gómez Gutiérrez, A.; Casas Vargas, L.A. Diversidad del ADN Mitocondrial en Restos óseos Prehispánicos Asociados al Templo del Sol en los Andes Orientales Colombianos; Pontificia Universidad Javeriana: Bogotá, Colombia, 2017. [Google Scholar]
  49. Silva, A.; Briceño, I.; Burgos, J.; Torres, D.; Villegas, V.; Gómez, A.; Rodríguez Cuenca, J.V. Mitochondrial DNA analysis on pre-Columbian bone remains of the Herrera period. Biomédica 2008, 28, 569–577. [Google Scholar] [CrossRef]
  50. Miller, M.J.; Agarwal, S.C.; Aristizabal, L.; Langebaek, C. The daily grind: Sex-and age-related activity patterns inferred from cross-sectional geometry of long bones in a pre-Columbian muisca population from Tibanica, Colombia. Am. J. Phys. Anthropol. 2018, 167, 311–326. [Google Scholar] [CrossRef]
  51. Monsalve, M.V.; Cardenas, F.; Guhl, F.; Delaney, A.D.; Devine, D.V. Phylogenetic analysis of mtDNA lineages in South American mummies. Ann. Hum. Genet. 1996, 60, 293–303. [Google Scholar] [CrossRef] [PubMed]
  52. Usme-Romero, S.; Alonso, M.; Hernandez-Cuervo, H.; Yunis, E.J.; Yunis, J.J. Genetic differences between Chibcha and Non-Chibcha speaking tribes based on mitochondrial DNA (mtDNA) haplogroups from 21 Amerindian tribes from Colombia. Genet. Mol. Biol. 2013, 36, 149–157. [Google Scholar] [CrossRef] [PubMed]
  53. Hoopes, J.W.; Fonseca, O. Goldwork and Chibchan identity: Endogenous change and diffuse unity in the Isthmo-Colombian area. In Gold and Power in Ancient Costa Rica, Panama, and Colombia; Dumbarton Oaks Research Library and Collection: Washington, DC, USA, 2003; pp. 49–89. [Google Scholar]
  54. Ward, R.H.; Redd, A.; Valencia, D.; Frazier, B.; Pääbo, S. Genetic and linguistic differentiation in the Americas. Proc. Natl. Acad. Sci. USA 1993, 90, 10663–10667. [Google Scholar] [CrossRef] [PubMed]
  55. Keyeux, G.; Rodas, C.; Gelvez, N.; Carter, D. Possible migration routes into South America deduced from mitochondrial DNA studies in Colombian Amerindian populations. Hum. Biol. 2002, 74, 211–233. [Google Scholar] [CrossRef] [PubMed]
  56. Schurr, T.G.; Ballinger, S.W.; Gan, Y.Y.; Hodge, J.A.; Merriwether, D.A.; Lawrence, D.N.; Knowler, W.C.; Weiss, K.M.; Wallace, D.C. Amerindian mitochondrial DNAs have rare Asian mutations at high frequencies, suggesting they derived from four primary maternal lineages. Am. J. Hum. Genet. 1990, 46, 613. [Google Scholar] [PubMed]
  57. Xavier, C.; Builes, J.J.; Gomes, V.; Ospino, J.M.; Aquino, J.; Parson, W.; Amorim, A.; Gusmão, L.; Goios, A. Admixture and genetic diversity distribution patterns of non-recombining lineages of Native American ancestry in Colombian populations. PLoS ONE 2015, 10, e0120155. [Google Scholar] [CrossRef] [PubMed]
  58. Shimada, I.; Shinoda, K.; Farnum, J.; Corruccini, R.; Watanabe, H. An integrated analysis of pre-hispanic mortuary practices: A middle Sicn case study. Curr. Anthropol. 2004, 45, 369–402. [Google Scholar] [CrossRef]
  59. Shinoda, K.I.; Adachi, N.; Guillen, S.; Shimada, I. Mitochondrial DNA analysis of ancient Peruvian highlanders. Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol. 2006, 131, 98–107. [Google Scholar] [CrossRef]
  60. Kemp, B.M.; Tung, T.A.; Summar, M.L. Genetic continuity after the collapse of the Wari empire: Mitochondrial DNA profiles from Wari and post-Wari populations in the ancient Andes. Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol. 2009, 140, 80–91. [Google Scholar] [CrossRef]
  61. Luciani, S.; Fornaciari, G.; Rickards, O.; Labarga, C.M.; Rollo, F. Molecular characterization of a pre-Columbian mummy and in situ coprolite. Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol. 2006, 129, 620–629. [Google Scholar] [CrossRef] [PubMed]
  62. Valverde, G.; Barreto Romero, M.I.; Flores Espinoza, I.; Cooper, A.; Fehren-Schmitz, L.; Llamas, B.; Haak, W. Ancient DNA analysis suggests negligible impact of the Wari empire expansion in Peru’s central coast during the Middle Horizon. PLoS ONE 2016, 11, e0155508. [Google Scholar] [CrossRef] [PubMed]
  63. Fehren-Schmitz, L.; Reindel, M.; Cagigao, E.T.; Hummel, S.; Herrmann, B. Pre-Columbian population dynamics in coastal southern Peru: A diachronic investigation of mtDNA patterns in the Palpa region by ancient DNA analysis. Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol. 2010, 141, 208–221. [Google Scholar] [CrossRef] [PubMed]
  64. Lewis, C.M.; Buikstra, J.E.; Stone, A.C. Ancient DNA and genetic continuity in the South Central Andes. Lat. Am. Antiq. 2007, 18, 145–160. [Google Scholar] [CrossRef]
  65. Moraga, M.; Santoro, C.M.; Standen, V.G.; Carvallo, P.; Rothhammer, F. Microevolution in prehistoric Andean populations: Chronologic mtDNA variation in the desert valleys of northern Chile. Am. J. Phys. Anthropol. Off. Publ. Am. Assoc. Phys. Anthropol. 2005, 127, 170–181. [Google Scholar] [CrossRef] [PubMed]
  66. Carnese, F.R.; Mendisco, F.; Keyser, C.; Dejean, C.B.; Dugooujon, J.-M.; Bravi, C.M.; Ludes, B.; Crubézy, E. Paleogenetical study of pre-Columbian samples from Pampa Grande (Salta Argentina). Am. J. Phys. Anthropol. 2010, 141, 452–462. [Google Scholar] [CrossRef] [PubMed]
  67. Manríquez Soto, G.R.; Moraga Vergara, M.; Santoro, C.; Aspillaga Fontaine, E.; Arriaza, B.T.; Rothhammer Engel, F. Morphometric and mtdna analyses of archaic skeletal remains from southwestern South America. Chungara Rev. De Antropol. Chil. 2011, 43, 283–292. [Google Scholar] [CrossRef]
  68. Ginther, C.; Corach, D.; Penacino, G.A.; Rey, J.A.; Carnese, F.R.; Hutz, M.H.; Anderson, A.; Just, J.; Salzano, F.M.; King, M.C. Genetic variation among the Mapuche Indians from the Patagonian region of Argentina: Mitochondrial DNA sequence variation and allele frequencies of several nuclear genes. In DNA Fingerprinting: State of the Science; Springer: Berlin/Heidelberg, Germany, 1993; pp. 211–219. [Google Scholar]
  69. Lewis, C.M.; Lizárraga, B.; Tito, R.Y.; López, P.W.; Iannacone, G.C.; Medina, A.; Martínez, R.; Polo, S.I.; Augusto, F.; Cáceres, A.M.; et al. Mitochondrial DNA and the peopling of South America. Hum. Biol. 2007, 79, 159–178. [Google Scholar] [CrossRef]
  70. Wang, S.; Lewis, C.M., Jr.; Jakobsson, M.; Ramachandran, S.; Ray, N.; Bedoya, G.; Rojas, W.; Parra, M.V.; Molina, J.A.; Ruiz-Linares, A. Genetic variation and population structure in Native Americans. PLoS Genet. 2007, 3, e185. [Google Scholar] [CrossRef]
  71. Torroni, A.; Schurr, T.G.; Cabell, M.F.; Brown, M.D.; Neel, J.V.; Larsen, M.; Smith, D.G.; Vullo, C.M.; Wallace, D.C. Asian affinities and continental radiation of the four founding Native American mtDNAs. Am. J. Hum. Genet. 1993, 53, 563. [Google Scholar]
  72. Ward, R.H.; Salzano, F.M.; Bonatto, S.L.; Hutz, M.H.; Coimbra, C.E.A., Jr.; Santos, R.V. Mitochondrial DNA polymorphism in three Brazilian Indian tribes. Am. J. Hum. Biol. Off. J. Hum. Biol. Assoc. 1996, 8, 317–323. [Google Scholar] [CrossRef]
  73. Huber, N.; Parson, W.; Dür, A. Next generation database search algorithm for forensic mitogenome analyses. Forensic Sci. International. Genet. 2018, 37, 204–214. [Google Scholar] [CrossRef] [PubMed]
  74. Capodiferro, M.R.; Aram, B.; Raveane, A.; Migliore, N.R.; Colombo, G.; Ongaro, L.; Rivera, J.; Mendizábal, T.; Hernández-Mora, I.; Achilli, A.; et al. Archaeogenomic distinctiveness of the Isthmo-Colombian area. Cell 2021, 184, 1706–1723. [Google Scholar] [CrossRef] [PubMed]
  75. Achilli, A.; Perego, U.A.; Bravi, C.M.; Coble, M.D.; Kong, Q.P.; Woodward, S.R.; Salas, A.; Torroni, A.; Bandelt, H.J. The phylogeny of the four pan-American MtDNA haplogroups: Implications for evolutionary and disease studies. PLoS ONE 2008, 3, e1764. [Google Scholar] [CrossRef] [PubMed]
  76. Duggan, A.T.; Harris, A.J.T.; Marciniak, S.; Marshall, I.; Kuch, M.; Kitchen, A.; Renaud, G.; Southon, J.; Fuller, B.; Young, J.; et al. Genetic Discontinuity between the Maritime Archaic and Beothuk Populations in Newfoundland, Canada. Curr. Biol. CB 2017, 27, 3149–3156.e11. [Google Scholar] [CrossRef] [PubMed]
  77. Lindo, J.; Achilli, A.; Perego, U.A.; Archer, D.; Valdiosera, C.; Petzelt, B.; Mitchell, J.; Worl, R.; Dixon, E.J.; Malhi, R.S.; et al. Ancient individuals from the North American Northwest Coast reveal 10,000 years of regional genetic continuity. Proc. Natl. Acad. Sci. USA 2017, 114, 4093–4098. [Google Scholar] [CrossRef] [PubMed]
  78. Perego, U.A.; Lancioni, H.; Tribaldos, M.; Angerhofer, N.; Ekins, J.E.; Olivieri, A.; Woodward, S.R.; Pascale, J.M.; Cooke, R.; Motta, J.; et al. Decrypting the mitochondrial gene pool of modern Panamanians. PLoS ONE 2012, 7, e38337. [Google Scholar] [CrossRef]
  79. Gayden, T.; Perez, A.; Persad, P.J.; Bukhari, A.; Chennakrishnaiah, S.; Simms, T.; Maloney, T.; Rodriguez, K.; Herrera, R.J. The Himalayas: Barrier and conduit for gene flow. Am. J. Phys. Anthropol. 2013, 151, 169–182. [Google Scholar] [CrossRef]
  80. Kumar, S.; Ravuri, R.R.; Koneru, P.; Urade, B.P.; Sarkar, B.N.; Chandrasekar, A.; Rao, V.R. Reconstructing Indian-Australian phylogenetic link. BMC Evol. Biol. 2009, 9, 173. [Google Scholar] [CrossRef]
  81. Chandrasekar, A.; Kumar, S.; Sreenath, J.; Sarkar, B.N.; Urade, B.P.; Mallick, S.; Bandopadhyay, S.S.; Barua, P.; Barik, S.S.; Rao, V.R.; et al. Updating phylogeny of mitochondrial DNA macrohaplogroup m in India: Dispersal of modern human in South Asian corridor. PLoS ONE 2009, 4, e7447. [Google Scholar] [CrossRef]
  82. Bhandari, S.; Zhang, X.; Cui, C.; Bianba Liao, S.; Peng, Y.; Zhang, H.; Xiang, K.; Shi, H.; Qi, X.; Su, B.; et al. Genetic evidence of a recent Tibetan ancestry to Sherpas in the Himalayan region. Sci. Rep. 2015, 5, 16249. [Google Scholar] [CrossRef] [PubMed]
  83. Lee, H.Y.; Song, I.; Ha, E.; Cho, S.B.; Yang, W.I.; Shin, K.J. mtDNAmanager: A Web-based tool for the management and quality analysis of mitochondrial DNA control-region sequences. BMC Bioinform. 2008, 9, 483. [Google Scholar] [CrossRef] [PubMed]
  84. Malyarchuk, B.A.; Derenko, M.V. Evaluating the role of selection in the evolution of mitochondrial genomes of aboriginal peoples of Siberia. Vavilov J. Genet. Breed. 2023, 27, 218. [Google Scholar] [CrossRef] [PubMed]
  85. Khan, K.; Siddiqi, M.H.; Ali, S.; Ali, S.; Sabar, M.F. Mitochondrial DNA control region variants analysis in Balti population of Gilgit-Baltistan, Pakistan. Meta Gene 2020, 23, 100630. [Google Scholar] [CrossRef]
  86. Starikovskaya, E.B.; Sukernik, R.I.; Derbeneva, O.A.; Volodko, N.V.; Ruiz-Pesini, E.; Torroni, A.; Brown, M.D.; Lott, M.T.; Hosseini, S.H.; Wallace, D.C.; et al. Mitochondrial DNA diversity in indigenous populations of the southern extent of Siberia, and the origins of Native American haplogroups. Ann. Hum. Genet. 2005, 69, 67–89. [Google Scholar] [CrossRef] [PubMed]
  87. Der Sarkissian, C.; Brotherton, P.; Balanovsky, O.; Templeton, J.E.; Llamas, B.; Soubrier, J.; Moiseyev, V.; Khartanovich, V.; Cooper, A.; Genographic Consortium; et al. Mitochondrial genome sequencing in Mesolithic North East Europe Unearths a new sub-clade within the broadly distributed human haplogroup C1. PLoS ONE 2014, 9, e87612. [Google Scholar]
  88. Russo, M.G.; Dejean, C.B.; Avena, S.A.; Seldes, V.; Ramundo, P. Mitochondrial lineage A2ah found in a pre-Hispanic individual from the Andean region. Am. J. Hum. Biol. 2018, 30, e23134. [Google Scholar] [CrossRef]
  89. Russo, M.G.; Arencibia, V.; Emery, M.; Bettera Marcat, G.; Seldes, V.; Mercolli, P.; Soria, S.; Maldonado, L.; Kamenetzky, L.; Stone, A.C.; et al. Ancient mitochondrial genome diversity in South America: Contributions from Quebrada del Toro, Northwestern Argentina. Am. J. Biol. Anthropol. 2023, 181, 597–610. [Google Scholar] [CrossRef]
  90. Bandelt, H.J.; Herrnstadt, C.; Yao, Y.G.; Kong, Q.P.; Kivisild, T.; Rengo, C.; Scozzari, R.; Richards, M.; Villems, R.; Macaulay, V.; et al. Identification of Native American founder mtDNAs through the analysis of complete mtDNA sequences: Some caveats. Ann. Hum. Genet. 2003, 67, 512–524. [Google Scholar] [CrossRef]
  91. Sala, A.; Caputo, M.; Ginart, S.; Theiler, G.; Parolin, M.L.; Carnese, R.F.; Fainboim, L.; Corach, D. Historical records under the genetic evidence: “Chiriguano” tribe genesis as a test case. Mol. Biol. Rep. 2018, 45, 987–1000. [Google Scholar] [CrossRef]
  92. Cui, Y.; Lindo, J.; Hughes, C.E.; Johnson, J.W.; Hernandez, A.G.; Kemp, B.M.; Ma, J.; Cunningham, R.; Petzelt, B.; Mitchell, J.; et al. Ancient DNA Analysis of Mid-Holocene Individuals from the Northwest Coast of North America Reveals Different Evolutionary Paths for Mitogenomes. PLoS ONE 2013, 8, e66948. [Google Scholar] [CrossRef] [PubMed]
  93. Wood, M.R.; Sturk-Andreaggi, K.; Ring, J.D.; Huber, N.; Bodner, M.; Crawford, M.H.; Parson, W.; Marshall, C. Resolving mitochondrial haplogroups B2 and B4 with next-generation mitogenome sequencing to distinguish Native American from Asian haplotypes. Forensic Sci. Int. Genet. 2019, 43, 102143. [Google Scholar] [CrossRef] [PubMed]
  94. Flores-Huacuja, M.; Snow, M.; Ramos-Madrigal, J.; Contreras-Cubas, C.; Barajas-Olmos, F.; Gonzalez-Oliver, A.; Mendoza-Caamal, E.; Ciceron-Arellano, I.; Centeno-Cruz, F.; Orozco, L.; et al. Whole mitogenome analysis highlights demographic history and shared connections among distal Indigenous groups of Mexico Complete mitogenome sequencing from 60 Mexican Native American groups. bioRxiv 2023. bioRxiv:2023–09. [Google Scholar]
  95. Aqil, A.; Gill, S.; Gokcumen, O.; Malhi, R.S.; Reese, E.A.; Smith, J.L.; Heaton, T.T.; Lindqvist, C. A paleogenome from a Holocene individual supports genetic continuity in Southeast Alaska. Iscience 2023, 26. [Google Scholar] [CrossRef] [PubMed]
  96. Malhi, R.S.; Breece, K.E.; Shook, B.A.S.; Kaestle, F.A.; Chatters, J.C.; Hackenberger, S.; Smith, D.G. Patterns of mtDNA diversity in northwestern North America. Hum. Biol. 2004, 76, 33–54. [Google Scholar] [CrossRef] [PubMed]
  97. Mannis, O.; Kayser, M. Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum. Mutat. 2009, 30, E386–E394. [Google Scholar] [CrossRef]
  98. Taboada-Echalar, P.; Alvarez-Iglesias, V.; Heinz, T.; Vidal-Bralo, L.; Gómez-Carballa, A.; Catelli, L.; Pardo-Seco, J.; Pastoriza, A.; Carracedo, A.; Torres-Balanza, A.; et al. The genetic legacy of the pre-colonial period in contemporary Bolivians. PLoS ONE 2013, 8, e58980. [Google Scholar] [CrossRef]
  99. Behar, D.M.; Rosset, S.; Blue-Smith, J.; Balanovsky, O.; Tzur, S.; Comas, D.; Mitchell, R.J.; Quintana-Murci, L.; Tyler-Smith, C.; Genographic Consortium; et al. The Genographic Project public participation mitochondrial DNA database. PLoS Genet. 2007, 3, e104. [Google Scholar] [CrossRef]
  100. Kong, Q.P.; Yao, Y.G.; Sun, C.; Bandelt, H.J.; Zhu, C.L.; Zhang, Y.P. Phylogeny of East Asian mitochondrial DNA lineages inferred from complete sequences. Am. J. Hum. Genet. 2003, 73, 671–676. [Google Scholar] [CrossRef]
Figure 1. Geographic overview of excavation sites for pre-Hispanic human remains in the department of Antioquia, municipalities of Medellín, and Envigado.
Figure 1. Geographic overview of excavation sites for pre-Hispanic human remains in the department of Antioquia, municipalities of Medellín, and Envigado.
Genes 14 02036 g001
Figure 2. Phylogenetic tree (samples classified in haplogroups and haplotypes). All sorted samples are combined into a resulting tree that includes all rCRS-related polymorphisms. The polymorphisms marked in blue are private polymorphisms for this group, already known by Phylotree. Individuals are marked in green. Based on Phylotree [16], the phylogenetic tree was constructed using the neighbor-joining algorithm. The samples names have been abbreviated for clarity.
Figure 2. Phylogenetic tree (samples classified in haplogroups and haplotypes). All sorted samples are combined into a resulting tree that includes all rCRS-related polymorphisms. The polymorphisms marked in blue are private polymorphisms for this group, already known by Phylotree. Individuals are marked in green. Based on Phylotree [16], the phylogenetic tree was constructed using the neighbor-joining algorithm. The samples names have been abbreviated for clarity.
Genes 14 02036 g002
Figure 3. Network analysis of ten HVS-I sequences was performed using the reduced median algorithm. The star contraction option was applied to condense extensive data by identifying phylogenetic clusters into a single ancestral type.
Figure 3. Network analysis of ten HVS-I sequences was performed using the reduced median algorithm. The star contraction option was applied to condense extensive data by identifying phylogenetic clusters into a single ancestral type.
Genes 14 02036 g003
Figure 4. MDS analysis comparing haplotype frequencies of pre-Hispanic populations. AP: ancient population; CP: current population (the stress value for the matrix is 0.05841, with an RSQ of 0.98776).
Figure 4. MDS analysis comparing haplotype frequencies of pre-Hispanic populations. AP: ancient population; CP: current population (the stress value for the matrix is 0.05841, with an RSQ of 0.98776).
Genes 14 02036 g004
Figure 5. Geographic locations of analyzed populations in the Americas. These zones are color-coded to simplify identification: blue: Arctic and Subarctic; green: northwestern United States; orange: Great Basin; yellow: southwestern United States; purple: Mesoamerican regions; gray: Caribbean; dark blue: Colombia and light green: South America.
Figure 5. Geographic locations of analyzed populations in the Americas. These zones are color-coded to simplify identification: blue: Arctic and Subarctic; green: northwestern United States; orange: Great Basin; yellow: southwestern United States; purple: Mesoamerican regions; gray: Caribbean; dark blue: Colombia and light green: South America.
Genes 14 02036 g005
Table 1. Notable Colombian archaeological sites dated between the late Pleistocene and early Holocene.
Table 1. Notable Colombian archaeological sites dated between the late Pleistocene and early Holocene.
PlaceRegionDates BP Not Calibrated
Tocaima-CundinamarcaPubenza16,400 ± 420
Zipaquirá-CundinamarcaEl Abra12,400 ± 160
Tocancipá-CundinamarcaTibitó 111,740 ± 110
Tequerdama-CundinamarcaTequendama 1 y 29740 ± 135
Nemocón-CundinamarcaChecua8500–3000
Roncesvalles -TolimaJordan12.910 ± 60–9760 ± 160
Porce 045-AntioquiaPorce 7.080 ± 130
La Morena Envigado-AntioquiaPorce River High Basin10,060 ± 60–9680 ± 60
Porce, Primavera II-AntioquiaPorce River Low Basin7730 ± 170
Medellín Guayabal-Antioquia La Colinita880 ± 30
Medellín, Robledo-AntioquiaCerro El Volador T10 950 ± 70
Medellín, Robledo-AntioquiaCerro El Volador T5 530 ± 80
Medellín, Robledo-AntioquiaCerro El Volador T13 480 ± 60
Medellín, Robledo-AntioquiaCerro El Volador T9 420 ± 50
Medellín, Robledo-AntioquiaCerro El Volador T8 330 ± 60
Envigado, Loma del Barro-AntioquiaAlto de las Flores140 ± 30
Table 2. Results of haplogroups and haplotypes identification based on the HVS-I sequence for each of the analyzed samples with the gender assignment, the origin of the excavation site and type of tissue used for DNA extraction.
Table 2. Results of haplogroups and haplotypes identification based on the HVS-I sequence for each of the analyzed samples with the gender assignment, the origin of the excavation site and type of tissue used for DNA extraction.
Sample GenderSite DescriptionRangeHaplogroup QualityPolymorphisms
Abu 1Indet.La Colinita Molar 16,024–16,569A, A286%16111T 16223T 16290T
Abu 2Indet.La Colinita Molar 16,024–16,569; 1–576A, A2af1b1b60%16168T 16223T 16290T 16362C
Abu 3Indet.La Colinita Molar 16,024–16,569A94%16129A 16223T 16290T 16319A 16362C
Abu 4Indet.La Colinita Molar 16,024–16,569C97%16223T 16288C 16298C 16311C 16327T
Abu 5Indet.La Colinita Molar 16,024–16,569A, A2i100%16111T 16223T 16290T 16319A 16325C 16362C
Abu 6Indet.La Colinita Molar 16,024–16,569; 1–576D, D4j1a65%16086C 16183C 16223T
Abu 7Indet.La Colinita Molar 16,024–16,569;C, C1c494%16223T 16274A 16327Y
Abu 8Indet.La Colinita Molar 16,024–16,569A, A2+(64)+ @1611182%16290T 16319A 16362C
Abu 9Female Cerro el VoladorBone remains16,024–16,569; 1–576A, A2ah81%16098G 16111T 16223T 16290T
Abu 10Indet.Alto de las Flores Bone remains16,024–16,569; 1–576B2, B2s62%16217C 16325C
Table 3. Measures of diversity in the Aburráes population from the mitochondrial sequences of the HVS-I. Where k number of different sequences; H, haplotype/genetic diversity; π, nucleotide diversity and Tajima’s index (-D) and p value.
Table 3. Measures of diversity in the Aburráes population from the mitochondrial sequences of the HVS-I. Where k number of different sequences; H, haplotype/genetic diversity; π, nucleotide diversity and Tajima’s index (-D) and p value.
PopulationKHπTajima’sp-Value
Aburráes100.6440.012−0.729370.254
Table 4. Table of haplogroups, showing genetic diversity across ancient and contemporary populations throughout the American continent. The genetic diversity (h) for various archaeological sites in different geographic zones, including the Aburráes region in Colombia. Sites with “N.C.” values indicate that there is no genetic diversity data available for those sites. AP = 0 for modern indigenous populations.
Table 4. Table of haplogroups, showing genetic diversity across ancient and contemporary populations throughout the American continent. The genetic diversity (h) for various archaeological sites in different geographic zones, including the Aburráes region in Colombia. Sites with “N.C.” values indicate that there is no genetic diversity data available for those sites. AP = 0 for modern indigenous populations.
Position on the MapGeographic Zone Archaeological SiteAPNABCDXMHReference
1Arctic, SubarcticEarly Aleuts3500–1200118003000.436[17]
2Arctic, SubarcticLate Pre-Contact Aleuts1000–40052120040000.362[17]
3Arctic, SubarcticBrooks River938–131885201000.607[18]
4Arctic, SubarcticHot Springs20703100200N.C.[18]
5Arctic, SubarcticOn-Your-Knees Cave10,3001000100N.C.[19]
6Arctic, SubarcticMink Island1029–121561005000.333[18]
7Arctic, SubarcticTatshenshini-Alsek Glacier5501100000N.C.[20]
8Arctic, SubarcticThule1130–6281515000000[21]
9Arctic, SubarcticSadlermiut977–6821810008000.523[21]
10Arctic, SubarcticDorset2260–12162000200N.C.[21]
11Northwest United StatesChina Lake and Big Bar Lake49753100002N.C.[19]
12Northwest United StatesPlateau Salish200113611000.673[19]
13Northwest United StatesPlateau Sahaptian20080422000.714[19]
14Northwest United StatesWishram2003371709000.635[19]
15Northwest United StatesVantage500–150072013100.81[19]
16Northwest United StatesPaisley 5 Mile Point Caves14,0002110000N.C.[22]
17Great Basin, Southwest United StatesCecil 3600–2860 160196000.575[23]
18Great Basin, Southwest United StatesCook 200023121010000.64[23]
19Great Basin, Southwest United StatesApplegate 1765–2055 60240000.533[23]
20Great Basin, Southwest United StatesPyramid Lake 860–5905 1826010000.601[24]
21Great Basin, Southwest United StatesPyramid Lake (Wizards Beach) 92001001000N.C.[24]
22Great Basin, Southwest United StatesStillwater Marsh 290–3290 2118012000.552[24]
23Great Basin, Southwest United StatesHourglass Cave 80001010000N.C.[25]
24Great Basin, Southwest United StatesTommy Site 850–1150 3612555000.492[26]
25Great Basin, Southwest United StatesMine Canyon 650–850 127410000.591[26]
26Great Basin, Southwest United StatesFremont 500–1500 3002442000.349[27]
27Great Basin, Southwest United StatesAnasazi 1010–2010 3842770000.462[27]
28Great Basin, Southwest United StatesWestern Basketmaker II 2500–1300 2331811000.383[27,28]
29Northeast United StatesGreat Western Park 80062040000.533[29]
30Northeast United StatesGlacial Kame 29001831131000.601[29]
31Northeast United StatesMorse 270091350000.639[29]
32Northeast United StatesOrendorf 800115033000.709[29]
33Northeast United StatesNorris Farms 7001083413469600.702[25]
34Northeast United StatesPete Klunk Mound Group 18253995195100.694[30]
35Northeast United StatesSchild Mississippian 90047186114800.762[31]
36Northeast United StatesSchild Late Woodland 1200195141800.743[31]
37Northeast United StatesOhio Hopewell Mound Group 170034143107000.715[32]
38MesoamericaTlatelolco Post-Classic Aztec 500–6752315314000.574[33]
39MesoamericaMaya (Xcaret) 480–1400 2420120000.236[34]
40MesoamericaMaya (Copan) 750–1300 90081000.222[35]
41CaribbeanLa Caleta (Tainos) 1330–320 2400186000.391[36]
42CaribbeanCuba (Ciboneys) 4700–1620 151095000.562[37]
43CaribbeanBoruca01431001000.473[38]
44CaribbeanHuetar02719107000.453[39]
45CaribbeanNgöbe046311500000.449[40]
46CaribbeanTainos01900154000.351[41]
47Colombia La Purnia (Guane)1090176704000.875[42]
48Colombia Jerico (Lache)70052300000.600[43]
49Colombia Andes Orientales Colombianos (Muisca, Lache, Guane, Checua, Aguazuque) 1200–700131741000.654[44]
50Colombia Madrid (Muisca)200060600000[45]
51Colombia Candelaria (muisca)9001010000000[46]
52Colombia Antioquia (Aburraes)800106121000.644[43]
53Colombia Malambo 300060600000[47]
54Colombia Sogamoso (Muisca)800–1500 139310000.500[48]
55Colombia Nemocon (Checua)800053110000.700[3]
56Colombia Herrera (Muisca)1400–10077000000[49]
57Colombia Tibanica (Muisca)8001812204000.523[50]
58Colombia Agroalfarero (early Muisca)9002413641000.641[51]
59Colombia Páez 2 0361031012100.747[52]
60Colombia Valle del Cauca (Calima)800–1200175822000.706[43]
61Colombia Ovando (Quimbaya)1500–700 70601000.286[53]
62Colombia Desano0203395000.726[52]
63Colombia Curripaco0221983100.710[52]
64Colombia Wayuu03088130100.692[12,54]
65Colombia Guane0176704000.691[42]
66Colombia Waunana040219117100.682[55]
67Colombia Ingano0481917111000.679[52]
68Colombia Tule-Cuna03015860100.660[55,56]
69Colombia Tucano0141318100.659[52]
70Colombia Embera02121162000.657[52]
71Colombia Antioquia0381810010000.654[57]
72Colombia Maya02714642000.652[56]
73Colombia Embera02171012000.647[55]
74Colombia Tatuyo0104051000.644[52]
75Colombia Zenu034511171000.642[55]
76Colombia Chibcha 0156242100.781[52]
77Colombia Puinave061531204100.633[52]
78Colombia Cauca060115385100.560[57]
79Colombia Páez 1031182110000.551[55]
80Colombia Guane-Butar03342108000.538[55]
81Colombia Cubeo 0248276100.764[52]
82Colombia Guambiano080565910000.436[57]
83Colombia Arsario02820080000.423[12]
84Colombia Kuna063451800000.415[12]
85Colombia Kogui02117040000.324[12,54]
86Colombia Arhuaco02119110000.186[52]
87Colombia Ijka03128120000.185[12]
88Colombia Chimila05650112100.174[52]
89South AmericaNorth Peruvian Coast 1000216717000.729[58]
90South AmericaPeruvian Highlanders 550–4503532381000.523[59]
91South AmericaConchapata 1400–1200 144721000.692[60]
92South AmericaCuzco 1020–830 1000100N.C.[61]
93South AmericaHuari900–600 1834101000.647[62]
94South AmericaParacas (Peninsula) 2800–2200 100037000.467[63]
95South AmericaChen Chen 1215–1000 239941000.692[64]
96South AmericaParacas (Palpa) 2800–2200 2820422000.370[63]
97South AmericaNasca-Rural (Palpa) 2200–1400 3714824000.535[63]
98South AmericaNasca-Urban (Palpa) 2200–1400 28051211000.653[63]
99South AmericaMiddle Horizon (Palpa) 1400–1000 110344000.727[63]
100South AmericaPacapaccari (Highlands)8201601150000.458[63]
101South AmericaYacotogia 118725114100000.547[63]
102South AmericaOcoro 1400–1000 51310000.700[63]
103South AmericaBotigiriayocc 1000–600 120840000.485[63]
104South AmericaHuayuncalla 97851220000.800[63]
105South AmericaLayuni 1400–1000 90531000.639[63]
106South AmericaNorth Chile Late Archaic 6000–3900 147511000.659[65]
107South AmericaNorth Chile Middle Horizon 1650–1000 196850000.690[65]
108South AmericaNorth Chile Late Intermediate 1000–500 153831000.676[65]
109South AmericaPampa Grande 1600–1350 192908000.620[66]
110South AmericaBaño nuevo, Tagua Tagua, Camarones9000–75003091092000.729[67]
111South AmericaMapuche039615810000.740[68]
112South AmericaYungay02011054000.679[69]
113South AmericaQuechua02321371000.605[70]
114South AmericaYanomami053053115300.593[71]
115South AmericaZoro02960418100.571[72]
116South AmericaXavante02442000000.290[72]
Table 5. Distribution of mitochondrial haplogroups in the sampled population. Offering an insightful breakdown of the genetic haplogroups assigned to each individual in the study. This comprehensive presentation aids in the categorization and enhances our understanding of the genetic diversity within the sampled population.
Table 5. Distribution of mitochondrial haplogroups in the sampled population. Offering an insightful breakdown of the genetic haplogroups assigned to each individual in the study. This comprehensive presentation aids in the categorization and enhances our understanding of the genetic diversity within the sampled population.
Sample IDLocalityHaplogroupHaplotypeHaplotype Global Distribution
Abu 1La Colinita A216111T 16223T 16290TFound across the entire continent of the Americas [73]
Abu 2La Colinita A2af1b1b16168T 16223T 16290T 16362CSpecifically present in Panama and Costa Rica [74,75]
Abu 3La Colinita A16129A 16223T 16290T 16319A 16362CFound across the entire continent of the Americas [73]
Abu 4La Colinita C16223T 16288C 16298C 16311C 16327TPrevalent in most parts of the American continent [73]
Abu 5La Colinita A2i16111T 16223T 16290T 16319A 16325C 16362CConcentrated in central North America, including Nebraska, Minnesota, and South Dakota [73], found in ancient populations in Canada [76,77,78]
Abu 6La Colinita D4j1a16086C 16183C 16223TFound in a small number of individuals in Thailand and Tibet [73,79,80,81,82,83,84,85]
Abu 7La Colinita C1c416223T 16274A 16327YLocated in the southern region of Mexico [86,87]
Abu 8La Colinita A216290T 16319A 16362CFound across the entire continent of the Americas [73]
Abu 9Cerro el VoladorA2ah16098G 16111T 16223T 16290TPresent in North America but more concentrated in South America [73,88,89,90,91]
Abu 10Alto de las Flores B2s16217C 16325CSpecifically found in North America [73,92,93]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Uricoechea Patiño, D.; Collins, A.; Romero García, O.J.; Santos Vecino, G.; Aristizábal Espinosa, P.; Bernal Villegas, J.E.; Benavides Benitez, E.; Vergara Muñoz, S.; Briceño Balcázar, I. Unraveling the Genetic Threads of History: mtDNA HVS-I Analysis Reveals the Ancient Past of the Aburra Valley. Genes 2023, 14, 2036. https://doi.org/10.3390/genes14112036

AMA Style

Uricoechea Patiño D, Collins A, Romero García OJ, Santos Vecino G, Aristizábal Espinosa P, Bernal Villegas JE, Benavides Benitez E, Vergara Muñoz S, Briceño Balcázar I. Unraveling the Genetic Threads of History: mtDNA HVS-I Analysis Reveals the Ancient Past of the Aburra Valley. Genes. 2023; 14(11):2036. https://doi.org/10.3390/genes14112036

Chicago/Turabian Style

Uricoechea Patiño, Daniel, Andrew Collins, Oscar Julián Romero García, Gustavo Santos Vecino, Pablo Aristizábal Espinosa, Jaime Eduardo Bernal Villegas, Escilda Benavides Benitez, Saray Vergara Muñoz, and Ignacio Briceño Balcázar. 2023. "Unraveling the Genetic Threads of History: mtDNA HVS-I Analysis Reveals the Ancient Past of the Aburra Valley" Genes 14, no. 11: 2036. https://doi.org/10.3390/genes14112036

APA Style

Uricoechea Patiño, D., Collins, A., Romero García, O. J., Santos Vecino, G., Aristizábal Espinosa, P., Bernal Villegas, J. E., Benavides Benitez, E., Vergara Muñoz, S., & Briceño Balcázar, I. (2023). Unraveling the Genetic Threads of History: mtDNA HVS-I Analysis Reveals the Ancient Past of the Aburra Valley. Genes, 14(11), 2036. https://doi.org/10.3390/genes14112036

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop