Decrypting the mitochondrial gene pool of modern Panamanians

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

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.

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.

Mitochondrial DNA Footprints from Western Eurasia in Modern Mongolia

Mongolia is located in a strategic position at the eastern edge of the Eurasian Steppe. Nomadic populations moved across this wide area for millennia before developing more sedentary communities, extended empires, and complex trading networks, which connected western Eurasia and eastern Asia until the late Medieval period. We provided a fine-grained portrait of the mitochondrial DNA (mtDNA) variation observed in present-day Mongolians and capable of revealing gene flows and other demographic processes that took place in Inner Asia, as well as in western Eurasia. The analyses of a novel dataset (N = 2,420) of mtDNAs highlighted a clear matrilineal differentiation within the country due to a mixture of haplotypes with eastern Asian (EAs) and western Eurasian (WEu) origins, which were differentially lost and preserved. In a wider genetic context, the prevalent EAs contribution, larger in eastern and central Mongolian regions, revealed continuous connections with neighboring Asian populations until recent times, as attested by the geographically restricted haplotype-sharing likely facilitated by the Genghis Khan’s so-called Pax Mongolica. The genetic history beyond the WEu haplogroups, notably detectable on both sides of Mongolia, was more difficult to explain. For this reason, we moved to the analysis of entire mitogenomes (N = 147). Although it was not completely possible to identify specific lineages that evolved in situ, two major changes in the effective (female) population size were reconstructed. The more recent one, which began during the late Pleistocene glacial period and became steeper in the early Holocene, was probably the outcome of demographic events connected to western Eurasia. The Neolithic growth could be easily explained by the diffusion of dairy pastoralism, as already proposed, while the late glacial increase indicates, for the first time, a genetic connection with western Eurasian refuges, as supported by the unusual high frequency and internal sub-structure in Mongolia of haplogroup H1, a well-known post-glacial marker in Europe. Bronze Age events, without a significant demographic impact, might explain the age of some mtDNA haplogroups. Finally, a diachronic comparison with available ancient mtDNAs made it possible to link six mitochondrial lineages of present-day Mongolians to the timeframe and geographic path of the Silk Route.

Weaving Mitochondrial DNA and Y-Chromosome Variation in the Panamanian Genetic Canvas

The Isthmus of Panama was a crossroads between North and South America during the continent’s first peopling (and subsequent movements) also playing a pivotal role during European colonization and the African slave trade. Previous analyses of uniparental systems revealed significant sex biases in the genetic history of Panamanians, as testified by the high proportions of Indigenous and sub-Saharan mitochondrial DNAs (mtDNAs) and by the prevalence of Western European/northern African Y chromosomes. Those studies were conducted on the general population without considering any self-reported ethnic affiliations. Here, we compared the mtDNA and Y-chromosome lineages of a new sample collection from 431 individuals (301 males and 130 females) belonging to either the general population, mixed groups, or one of five Indigenous groups currently living in Panama. We found different proportions of paternal and maternal lineages in the Indigenous groups testifying to pre-contact demographic events and genetic inputs (some dated to Pleistocene times) that created genetic structure. Then, while the local mitochondrial gene pool was marginally involved in post-contact admixtures, the Indigenous Y chromosomes were differentially replaced, mostly by lineages of western Eurasian origin. Finally, our new estimates of the sub-Saharan contribution, on a more accurately defined general population, reduce an apparent divergence between genetic and historical data.

Evaluating the Impact of Sex-Biased Genetic Admixture in the Americas through the Analysis of Haplotype Data

A general imbalance in the proportion of disembarked males and females in the Americas has been documented during the Trans-Atlantic Slave Trade and the Colonial Era and, although less prominent, more recently. This imbalance may have left a signature on the genomes of modern-day populations characterised by high levels of admixture. The analysis of the uniparental systems and the evaluation of continental proportion ratio of autosomal and X chromosomes revealed a general sex imbalance towards males for European and females for African and Indigenous American ancestries. However, the consistency and degree of this imbalance are variable, suggesting that other factors, such as cultural and social practices, may have played a role in shaping it. Moreover, very few investigations have evaluated the sex imbalance using haplotype data, containing more critical information than genotypes. Here, we analysed genome-wide data for more than 5000 admixed American individuals to assess the presence, direction and magnitude of sex-biased admixture in the Americas. For this purpose, we applied two haplotype-based approaches, ELAI and NNLS, and we compared them with a genotype-based method, ADMIXTURE. In doing so, besides a general agreement between methods, we unravelled that the post-colonial admixture dynamics show higher complexity than previously described.

The Mitochondrial DNA Landscape of Modern Mexico

Mexico is a rich source for anthropological and population genetic studies with high diversity in ethnic and linguistic groups. The country witnessed the rise and fall of major civilizations, including the Maya and Aztec, but resulting from European colonization, the population landscape has dramatically changed. Today, the majority of Mexicans do not identify themselves as Indigenous but as admixed, and appear to have very little in common with their pre-Columbian predecessors. However, when the maternally inherited mitochondrial (mt)DNA is investigated in the modern Mexican population, this is not the case. Control region sequences of 2021 samples deriving from all over the country revealed an overwhelming Indigenous American legacy, with almost 90% of mtDNAs belonging to the four major pan-American haplogroups A2, B2, C1, and D1. This finding supports a very low European contribution to the Mexican gene pool by female colonizers and confirms the effectiveness of employing uniparental markers as a tool to reconstruct a country’s history. In addition, the distinct frequency and dispersal patterns of Indigenous American and West Eurasian clades highlight the benefit such large and country-wide databases provide for studying the impact of colonialism from a female perspective and population stratification. The importance of geographical database subsets not only for forensic application is clearly demonstrated.

An isthmus of isolation: The likely elevated prevalence of genetic disease in ancient Panama and implications for considering rare diseases in paleopathology

OBJECTIVE

This study considers the evidence for elevated frequencies of "rare" diseases in ancient Panama. Indications of population isolation by multidisciplinary sources allow for the possibility that rare inherited conditions may have been maintained at relatively high prevalences in the region due to gene flow restriction.

MATERIALS

A sample of 267 skeletal human remains with diverse demographical characteristics from Pre-Columbian archaeological sites throughout Panama.

METHODS

Remains were analyzed macroscopically and hard tissue developmental anomalies were documented.

RESULTS

Frequencies of developmental anomalies and hard tissue changes consistent with specific rare genetic diseases, such as osteogenesis imperfecta, on the comparatively few human remains recovered from pre-Columbian archaeological sites are elevated as compared with global averages.

CONCLUSIONS

The paleopathological evidence is concordant with a scenario of isolation in Pre-Columbian times and with an increased presence of genetic disorders in the population.

SIGNIFICANCE

This study advocates for the special consideration of rare diseases by paleopathologists in regions where populations may have experienced prolonged geographical or social isolation in the past.

LIMITATIONS

A dearth of local modern epidemiological data and low sample sizes of preserved human remains in certain regions of the country limited the possibilities of spatiotemporal comparisons of rare disease prevalence.

SUGGESTIONS FOR FURTHER RESEARCH

Further scrutiny of developmental anomalies of genetic origin on ancient Panamanian remains and biomolecular testing of remains for specific disorders should be pursued to confirm the findings of this study.

Archaeogenomic distinctiveness of the Isthmo-Colombian area

The recently enriched genomic history of Indigenous groups in the Americas is still meager concerning continental Central America. Here, we report ten pre-Hispanic (plus two early colonial) genomes and 84 genome-wide profiles from seven groups presently living in Panama. Our analyses reveal that pre-Hispanic demographic events contributed to the extensive genetic structure currently seen in the area, which is also characterized by a distinctive Isthmo-Colombian Indigenous component. This component drives these populations on a specific variability axis and derives from the local admixture of different ancestries of northern North American origin(s). Two of these ancestries were differentially associated to Pleistocene Indigenous groups that also moved into South America, leaving heterogenous genetic footprints. An additional Pleistocene ancestry was brought by a still unsampled population of the Isthmus (UPopI) that remained restricted to the Isthmian area, expanded locally during the early Holocene, and left genomic traces up to the present day.

HLA allele and haplotype frequencies in the Panamanian population

In this study, we report for the first time HLA allele and haplotype frequencies in the modern Panamanian population at a two-field (four digits) resolution level. Reported frequencies were calculated from genotype data for the HLA-A, -B, -C, -DPB1, -DQB1 and -DRB1 loci of 462 healthy unrelated Panamanian adults of Hispanic ethnicity. In addition to providing new insights on the allelic structure of the Panamanian population and its origin, these data are critical for better planning of healthcare strategies in the country and for future research exploring the association with certain chronic and infectious diseases.

Ancient DNA Reconstructs the Genetic Legacies of Precontact Puerto Rico Communities

Indigenous peoples have occupied the island of Puerto Rico since at least 3000 BC. Due to the demographic shifts that occurred after European contact, the origin(s) of these ancient populations, and their genetic relationship to present-day islanders, are unclear. We use ancient DNA to characterize the population history and genetic legacies of precontact Indigenous communities from Puerto Rico. Bone, tooth, and dental calculus samples were collected from 124 individuals from three precontact archaeological sites: Tibes, Punta Candelero, and Paso del Indio. Despite poor DNA preservation, we used target enrichment and high-throughput sequencing to obtain complete mitochondrial genomes (mtDNA) from 45 individuals and autosomal genotypes from two individuals. We found a high proportion of Native American mtDNA haplogroups A2 and C1 in the precontact Puerto Rico sample (40% and 44%, respectively). This distribution, as well as the haplotypes represented, supports a primarily Amazonian South American origin for these populations and mirrors the Native American mtDNA diversity patterns found in present-day islanders. Three mtDNA haplotypes from precontact Puerto Rico persist among Puerto Ricans and other Caribbean islanders, indicating that present-day populations are reservoirs of precontact mtDNA diversity. Lastly, we find similarity in autosomal ancestry patterns between precontact individuals from Puerto Rico and the Bahamas, suggesting a shared component of Indigenous Caribbean ancestry with close affinity to South American populations. Our findings contribute to a more complete reconstruction of precontact Caribbean population history and explore the role of Indigenous peoples in shaping the biocultural diversity of present-day Puerto Ricans and other Caribbean islanders.

Epidemiologic Characteristics of Suicide in Panama, 2007–2016

Background and objectives: We aim to describe the demographic characteristics associated with suicide in Panama, to estimate the suicide mortality rate and years of potential life lost (YPLL) to suicide, and to explore the correlation of suicide rates with the Multidimensional Poverty Index (MPI). We present a descriptive retrospective epidemiological report of suicide-related mortality (Panama, 2007–2016). Materials and Methods: Data were matched-merged to calculate unadjusted suicide mortality rates (overall, and by sex, age groups, and administrative region), YPLL, and coefficients (r) for the correlation of MPI and suicide rates. Results: There were 1475 deaths by suicide (86% among men, 47% between 20 and 39 years). The average mortality rate was estimated at 3.91 per 100,000 population with an average YPLL rate of 3.79 per 1000 population. There was a statistically significant trend to reduce YPLL over time (r = −0.93; p< 0.001). Exploratory analyses did not show a significant correlation between the MPI and suicide rates. Our study showed a 6:1 male-to-female ratio of suicide, mostly affecting the age groups of 20–29 and over 80 years. Conclusions: Exploratory analyses on the correlation of the MPI and the suicide rates did not achieve statistical significance, and alternative explanations, such as access to pesticides and alcohol, were further explored to inform potential interventions.

Social determinants and cardiovascular disease mortality in Panama, 2012-2016

Background

The aims of this study were to determine the geographic and time variation of social determinants of health (SDH) and cardiovascular disease (CVD) mortality in Panama from 2012 to 2016, and to identify which of the SDH has the strongest correlation with a socioeconomic index (SEI).

Methods

We conducted an ecological study obtaining mortality from the National Mortality Register and socioeconomic variables derived from the National Household Survey (NHS). The International Classification of Diseases 10th revision codes I20–I25 and I60–I69 were used for ischemic heart disease (IHD) and stroke, respectively. Standardized age-adjusted mortality rates were calculated by direct method. Mortality rates and socioeconomic variables were evaluated together in a panel data model. A SEI was developed from factorial analysis by principal components with a polychoric correlation matrix. Provinces and regions were categorized in tertiles according to median value of the SEI score.

Results

The NHS evaluated an average of 15,919 households per year. The mean of age throughout the study period was 41 years. The average monthly income increased, from US$ (SD) 331.94 (5.38) in 2012, to 406.24 (5.81) in 2016, whereas the social security health coverage remained in a range of 57–58%. The mean number of school years was twelve. Significant geographical and temporal variations in social determinants and mortality rates were observed throughout the country. Colon, categorized in the middle tertile according to the SEI, presented higher IHD mortality rates. Darién (in the lowest SEI tertile) Colón and Herrera had higher stroke mortality rates. The SEI categorized indigenous territories in the lowest tertile. Total years of education was the strongest correlated variable with the SEI, when we excluded the population living in indigenous territories. However, when this population was included, social security coverage had the strongest correlation with the SEI.

Conclusion

We observed geographical and temporal disparities in SDH and CVD mortality rates. Further epidemiological studies are warranted in the provinces of Colón, Darien, Herrera and Los Santos to explore in-depth the higher CVD mortality rates observed in these provinces.

Electronic supplementary material

The online version of this article (10.1186/s12889-019-6508-8) contains supplementary material, which is available to authorized users.

Análisis de ADN mitocondrial en una muestra de restos óseos precolombinos de Norte de Santander, Colombia (Área Cultural Chitarera)

Los análisis de ADN antiguo (ADNa) han incrementado en los últimos años permitiendo conocer la diversidad genética de las poblaciones precolombinas. En Colombia, existen pocos registros arqueológicos de la población prehispánica del Norte de Santander habitada en el siglo XVI por el grupo Chitarero. Por este motivo, nos propusimos analizar la diversidad genética a partir de secuencias de la región HVRI del ADNmt y determinar sus posibles relaciones con otras comunidades tanto antiguas como contemporáneas. Se analizaron siete individuos precolombinos asociados a este grupo prehispánico, recuperados en los municipios de Cácota y Silos en el departamento de Norte de Santander de los Andes Orientales colombianos, siguiendo criterios estrictos de autenticidad para el ADNa. En todos los individuos se logró identificar el haplogrupo B caracterizado por el polimorfismo en la posición 16217C, siendo éste uno de los más frecuentes en comunidades precolombinas y contemporáneas de los Andes Suramericanos. Este hallazgo indica que este grupo poblacional se encuentra estrechamente emparentado por línea materna, con posibles índices de endogamia, con una probable densidad demográfica baja y una baja diversidad genética, similares a lo observado en comunidades pertenecientes a periodos anteriores como el Formativo. Este grupo precolombino exhibe una de las diversidades genéticas más bajas reportadas en las poblaciones pertenecientes a la familia lingüistica Chibcha. Estos resultados genéticos coinciden con los planteamientos sobre el grupo Chitarero de pertenecer a comunidades pequeñas independientes, con asentamientos dispersos, apartados unos de otros. 

The Paleo-Indian Entry into South America According to Mitogenomes

Recent and compelling archaeological evidence attests to human presence ∼14.5 ka at multiple sites in South America and a very early exploitation of extreme high-altitude Andean environments. Considering that, according to genetic evidence, human entry into North America from Beringia most likely occurred ∼16 ka, these archeological findings would imply an extremely rapid spread along the double continent. To shed light on this issue from a genetic perspective, we first completely sequenced 217 novel modern mitogenomes of Native American ancestry from the northwestern area of South America (Ecuador and Peru); we then evaluated them phylogenetically together with other available mitogenomes (430 samples, both modern and ancient) from the same geographic area and, finally, with all closely related mitogenomes from the entire double continent. We detected a large number (N = 48) of novel subhaplogroups, often branching into further subclades, belonging to two classes: those that arose in South America early after its peopling and those that instead originated in North or Central America and reached South America with the first settlers. Coalescence age estimates for these subhaplogroups provide time boundaries indicating that early Paleo-Indians probably moved from North America to the area corresponding to modern Ecuador and Peru over the short time frame of ∼1.5 ka comprised between 16.0 and 14.6 ka.

Successful reconstruction of whole mitochondrial genomes from ancient Central America and Mexico

The northern and southern peripheries of ancient Mesoamerica are poorly understood. There has been speculation over whether borderland cultures such as Greater Nicoya and Casas Grandes represent Mesoamerican outposts in the Isthmo-Colombian area and the Greater Southwest, respectively. Poor ancient DNA preservation in these regions challenged previous attempts to resolve these questions using conventional genetic techniques. We apply advanced in-solution mitogenome capture and high-throughput sequencing to fourteen dental samples obtained from the Greater Nicoya sites of Jícaro and La Cascabel in northwest Costa Rica (n = 9; A.D. 800–1250) and the Casas Grandes sites of Paquimé and Convento in northwest Mexico (n = 5; A.D. 1200–1450). Full mitogenome reconstruction was successful for three individuals from Jícaro and five individuals from Paquimé and Convento. The three Jícaro individuals belong to haplogroup B2d, a haplogroup found today only among Central American Chibchan-speakers. The five Paquimé and Convento individuals belong to haplogroups C1c1a, C1c5, B2f and B2a which, are found in contemporary populations in North America and Mesoamerica. We report the first successfully reconstructed ancient mitogenomes from Central America, and the first genetic evidence of ancestry affinity of the ancient inhabitants of Greater Nicoya and Casas Grandes with contemporary Isthmo-Columbian and Greater Southwest populations, respectively.

Enclaves of genetic diversity resisted Inca impacts on population history

The Inca Empire is claimed to have driven massive population movements in western South America, and to have spread Quechua, the most widely-spoken language family of the indigenous Americas. A test-case is the Chachapoyas region of northern Peru, reported as a focal point of Inca population displacements. Chachapoyas also spans the environmental, cultural and demographic divides between Amazonia and the Andes, and stands along the lowest-altitude corridor from the rainforest to the Pacific coast. Following a sampling strategy informed by linguistic data, we collected 119 samples, analysed for full mtDNA genomes and Y-chromosome STRs. We report a high indigenous component, which stands apart from the network of intense genetic exchange in the core central zone of Andean civilization, and is also distinct from neighbouring populations. This unique genetic profile challenges the routine assumption of large-scale population relocations by the Incas. Furthermore, speakers of Chachapoyas Quechua are found to share no particular genetic similarity or gene-flow with Quechua speakers elsewhere, suggesting that here the language spread primarily by cultural diffusion, not migration. Our results demonstrate how population genetics, when fully guided by the archaeological, historical and linguistic records, can inform multiple disciplines within anthropology.

[Mitochondrial DNA diversity in prehispanic bone remains on the eastern Colombian Andes]

INTRODUCTION

DNA extracted from ancient human bones allows to analyze the genetic makeup of pre-Columbian populations and to determine the dynamics that gave rise to the diversity of contemporary populations.

OBJECTIVE

To determine the genetic diversity of skeletal remains associated with the Templo del Sol (Sun Temple) and their relationship with other contemporary and ancient communities of America.

MATERIALS AND METHODS

We analyzed 13 individuals belonging to the pre-Columbian Muisca Period (IX-XVI centuries AD) from the vicinities of the Templo del Sol (Sun Temple) (Sogamoso, Boyacá) in the eastern Colombian Andes. Mitochondrial DNA was amplified and RFLPs were performed in order to type the four traditional Amerindian haplogroups (A, B, C and D). In addition, autosomal markers including amelogenin and Y-chromosome STRs were amplified.

RESULTS

Among the observed mitochondrial lineages, haplogroup A was the most frequent, followed by haplogroups B and C; no evidence of haplogroup D was found. The genetic variation analysis indicated a similar diversity of pre-ColumbianMuiscas to that of contemporary populations belonging to the Chibcha linguistic family from Colombia and Central America. Molecular sexing was accomplished and it was compared to osteological data. With only one exception, anthropological and molecular data were consistent.

CONCLUSIONS

Our results contribute new genetic elements supporting the hypothesis of Central American origin of the Chibcha groups of the Cundiboyacense plateau, and allowed sex typing and kinship evaluations.

Exploring the Y Chromosomal Ancestry of Modern Panamanians

Geologically, Panama belongs to the Central American land-bridge between North and South America crossed by Homo sapiens >14 ka ago. Archaeologically, it belongs to a wider Isthmo-Colombian Area. Today, seven indigenous ethnic groups account for 12.3% of Panama’s population. Five speak Chibchan languages and are characterized by low genetic diversity and a high level of differentiation. In addition, no evidence of differential structuring between maternally and paternally inherited genes has been reported in isthmian Chibchan cultural groups. Recent data have shown that 83% of the Panamanian general population harbour mitochondrial DNAs (mtDNAs) of Native American ancestry. Considering differential male/female mortality at European contact and multiple degrees of geographical and genetic isolation over the subsequent five centuries, the Y-chromosome Native American component is expected to vary across different geographic regions and communities in Panama. To address this issue, we investigated Y-chromosome variation in 408 modern males from the nine provinces of Panama and one indigenous territory (the comarca of Kuna Yala). In contrast to mtDNA data, the Y-chromosome Native American component (haplogroup Q) exceeds 50% only in three populations facing the Caribbean Sea: the comarca of Kuna Yala and Bocas del Toro province where Chibchan languages are spoken by the majority, and the province of Colón where many Kuna and people of mixed indigenous-African-and-European descent live. Elsewhere the Old World component is dominant and mostly represented by western Eurasian haplogroups, which signal the strong male genetic impact of invaders. Sub-Saharan African input accounts for 5.9% of male haplotypes. This reflects the consequences of the colonial Atlantic slave trade and more recent influxes of West Indians of African heritage. Overall, our findings reveal a local evolution of the male Native American ancestral gene pool, and a strong but geographically differentiated unidirectional sex bias in the formation of local modern Panamanian populations.

Genetic Diversity in the Lesser Antilles and Its Implications for the Settlement of the Caribbean Basin

Historical discourses about the Caribbean often chronicle West African and European influence to the general neglect of indigenous people's contributions to the contemporary region. Consequently, demographic histories of Caribbean people prior to and after European contact are not well understood. Although archeological evidence suggests that the Lesser Antilles were populated in a series of northward and eastern migratory waves, many questions remain regarding the relationship of the Caribbean migrants to other indigenous people of South and Central America and changes to the demography of indigenous communities post-European contact. To explore these issues, we analyzed mitochondrial DNA and Y-chromosome diversity in 12 unrelated individuals from the First Peoples Community in Arima, Trinidad, and 43 unrelated Garifuna individuals residing in St. Vincent. In this community-sanctioned research, we detected maternal indigenous ancestry in 42% of the participants, with the remainder having haplotypes indicative of African and South Asian maternal ancestry. Analysis of Y-chromosome variation revealed paternal indigenous American ancestry indicated by the presence of haplogroup Q-M3 in 28% of the male participants from both communities, with the remainder possessing either African or European haplogroups. This finding is the first report of indigenous American paternal ancestry among indigenous populations in this region of the Caribbean. Overall, this study illustrates the role of the region's first peoples in shaping the genetic diversity seen in contemporary Caribbean populations.

Trend Analysis of Cancer Mortality and Incidence in Panama, Using Joinpoint Regression Analysis

Cancer is one of the leading causes of death worldwide and its incidence is expected to increase in the future. In Panama, cancer is also one of the leading causes of death. In 1964, a nationwide cancer registry was started and it was restructured and improved in 2012. The aim of this study is to utilize Joinpoint regression analysis to study the trends of the incidence and mortality of cancer in Panama in the last decade. Cancer mortality was estimated from the Panamanian National Institute of Census and Statistics Registry for the period 2001 to 2011. Cancer incidence was estimated from the Panamanian National Cancer Registry for the period 2000 to 2009. The Joinpoint Regression Analysis program, version 4.0.4, was used to calculate trends by age-adjusted incidence and mortality rates for selected cancers. Overall, the trend of age-adjusted cancer mortality in Panama has declined over the last 10 years (-1.12% per year). The cancers for which there was a significant increase in the trend of mortality were female breast cancer and ovarian cancer; while the highest increases in incidence were shown for breast cancer, liver cancer, and prostate cancer. Significant decrease in the trend of mortality was evidenced for the following: prostate cancer, lung and bronchus cancer, and cervical cancer; with respect to incidence, only oral and pharynx cancer in both sexes had a significant decrease. Some cancers showed no significant trends in incidence or mortality. This study reveals contrasting trends in cancer incidence and mortality in Panama in the last decade. Although Panama is considered an upper middle income nation, this study demonstrates that some cancer mortality trends, like the ones seen in cervical and lung cancer, behave similarly to the ones seen in high income countries. In contrast, other types, like breast cancer, follow a pattern seen in countries undergoing a transition to a developed economy with its associated lifestyle, nutrition, and body weight changes.

Genomic insights on the ethno-history of the Maya and the 'Ladinos' from Guatemala

Background

Guatemala is a multiethnic and multilingual country located in Central America. The main population groups separate ‘Ladinos’ (mixed Native American-African-Spanish), and Native indigenous people of Maya descent. Among the present-day Guatemalan Maya, there are more than 20 different ethnic groups separated by different languages and cultures. Genetic variation of these communities still remains largely unexplored. The principal aim of this study is to explore the genetic variability of the Maya and ‘Ladinos’ from Guatemala by means of uniparental and ancestry informative markers (AIMs).

Results

Analyses of uniparental genetic markers indicate that Maya have a dominant Native American ancestry (mitochondrial DNA [mtDNA]: 100%; Y-chromosome: 94%). ‘Ladino’, however, show a clear gender-bias as indicated by the large European ancestry observed in the Y-chromosome (75%) compared to the mtDNA (0%). Autosomal polymorphisms (AIMs) also mirror this marked gender-bias: (i) Native American ancestry: 92% for the Maya vs. 55% for the ‘Ladino’, and (ii) European ancestry: 8% for the Maya vs. 41% for the ‘Ladino’. In addition, the impact of the Trans-Atlantic slave trade on the present-day Guatemalan population is very low (and only occurs in the ‘Ladino’; mtDNA: 9%; AIMs: 4%), in part mirroring the fact that Guatemala has a predominant orientation to the Pacific Ocean instead of a Caribbean one. Sequencing of entire Guatemalan mitogenomes has led to improved Native American phylogeny via the addition of new haplogroups that are mainly observed in Mesoamerica and/or the North of South America.

Conclusions

The data reveal the existence of a fluid gene flow in the Mesoamerican area and a predominant unidirectional flow towards South America, most likely occurring during the Pre-Classic (1800 BC-200 AD) and the Classic (200–1000 AD) Eras of the Mesoamerican chronology, coinciding with development of the most distinctive and advanced Mesoamerican civilization, the Maya. Phylogenetic features of mtDNA data also suggest a demographic scenario that is compatible with moderate local endogamy and isolation in the Maya combined with episodes of gene exchange between ethnic groups, suggesting an ethno-genesis in the Guatemalan Maya that is recent and supported on a cultural rather than a biological basis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1339-1) contains supplementary material, which is available to authorized users.

Indian signatures in the westernmost edge of the European Romani diaspora: new insight from mitogenomes

In agreement with historical documentation, several genetic studies have revealed ancestral links between the European Romani and India. The entire mitochondrial DNA (mtDNA) of 27 Spanish Romani was sequenced in order to shed further light on the origins of this population. The data were analyzed together with a large published dataset (mainly hypervariable region I [HVS-I] haplotypes) of Romani (N=1,353) and non-Romani worldwide populations (N>150,000). Analysis of mitogenomes allowed the characterization of various Romani-specific clades. M5a1b1a1 is the most distinctive European Romani haplogroup; it is present in all Romani groups at variable frequencies (with only sporadic findings in non-Romani) and represents 18% of their mtDNA pool. Its phylogeographic features indicate that M5a1b1a1 originated 1.5 thousand years ago (kya; 95% CI: 1.3-1.8) in a proto-Romani population living in Northwest India. U3 represents the most characteristic Romani haplogroup of European/Near Eastern origin (12.4%); it appears at dissimilar frequencies across the continent (Iberia: ≈ 31%; Eastern/Central Europe: ≈ 13%). All U3 mitogenomes of our Iberian Romani sample fall within a new sub-clade, U3b1c, which can be dated to 0.5 kya (95% CI: 0.3-0.7); therefore, signaling a lower bound for the founder event that followed admixture in Europe/Near East. Other minor European/Near Eastern haplogroups (e.g. H24, H88a) were also assimilated into the Romani by introgression with neighboring populations during their diaspora into Europe; yet some show a differentiation from the phylogenetically closest non-Romani counterpart. The phylogeny of Romani mitogenomes shows clear signatures of low effective population sizes and founder effects. Overall, these results are in good agreement with historical documentation, suggesting that cultural identity and relative isolation have allowed the Romani to preserve a distinctive mtDNA heritage, with some features linking them unequivocally to their ancestral Indian homeland.

Phylogenetic relationships of three Italian merino-derived sheep breeds evaluated through a complete mitogenome analysis

In Italy, the crisis of the wool industry triggered the necessity to reconvert the two traditional Merino-derived breeds, Gentile di Puglia and Sopravissana, to meat production, by creating the Merinizzata Italiana. The aim of the present study was to assess the genetic diversity of these three Italian Merino-derived (IMd) breeds by examining the molecular information encoded in the maternally-inherited mitochondrial DNA (mtDNA). A parallel molecular investigation was performed on the putative paternal and maternal breeds, the Merino from Spain and the Appenninica from Italy, respectively, as well as on three unrelated dairy breeds (Sarda and Comisana from Italy, and Lacaune from France). Firstly, the mtDNA control region of 291 samples was analyzed. When comparing the overall genetic distances among the eight stocks, the three IMd breeds clustered together close to the Appenninica, thus confirming its parental role. Among the 90 IMd samples, 82 different haplotypes were observed, almost all belonging to haplogroup B, and only one to A. For 23 mtDNAs, including nine IMd, the analysis was then brought to the level of entire mitogenomes. Three distinct sub-haplogroups within B were found to encompass the IMd samples, with one clade (B1a2a1) apparently restricted to those sheep. Thus, despite experiencing a drastic reduction in number (mainly due to changes in breeding practices driven by the economy), the IMd breeds still represent a reservoir of distinctive mitochondrial variants, which could potentially contribute to the development of conservation and management programs of Italian sheep breeds.

The first peopling of South America: new evidence from Y-chromosome haplogroup Q

Recent progress in the phylogenetic resolution of the Y-chromosome phylogeny permits the male demographic dynamics and migratory events that occurred in Central and Southern America after the initial human spread into the Americas to be investigated at the regional level. To delve further into this issue, we examined more than 400 Native American Y chromosomes (collected in the region ranging from Mexico to South America) belonging to haplogroup Q – virtually the only branch of the Y phylogeny observed in modern-day Amerindians of Central and South America – together with 27 from Mongolia and Kamchatka. Two main founding lineages, Q1a3a1a-M3 and Q1a3a1-L54(xM3), were detected along with novel sub-clades of younger age and more restricted geographic distributions. The first was also observed in Far East Asia while no Q1a3a1-L54(xM3) Y chromosome was found in Asia except the southern Siberian-specific sub-clade Q1a3a1c-L330. Our data not only confirm a southern Siberian origin of ancestral populations that gave rise to Paleo-Indians and the differentiation of both Native American Q founding lineages in Beringia, but support their concomitant arrival in Mesoamerica, where Mexico acted as recipient for the first wave of migration, followed by a rapid southward migration, along the Pacific coast, into the Andean region. Although Q1a3a1a-M3 and Q1a3a1-L54(xM3) display overlapping general distributions, they show different patterns of evolution in the Mexican plateau and the Andean area, which can be explained by local differentiations due to demographic events triggered by the introduction of agriculture and associated with the flourishing of the Great Empires.

Reconciling migration models to the Americas with the variation of North American native mitogenomes

In this study we evaluated migration models to the Americas by using the information contained in native mitochondrial genomes (mitogenomes) from North America. Molecular and phylogeographic analyses of B2a mitogenomes, which are absent in Eskimo-Aleut and northern Na-Dene speakers, revealed that this haplogroup arose in North America ∼11-13 ka from one of the founder Paleo-Indian B2 mitogenomes. In contrast, haplogroup A2a, which is typical of Eskimo-Aleuts and Na-Dene, but also present in the easternmost Siberian groups, originated only 4-7 ka in Alaska, led to the first Paleo-Eskimo settlement of northern Canada and Greenland, and contributed to the formation of the Na-Dene gene pool. However, mitogenomes also show that Amerindians from northern North America, without any distinction between Na-Dene and non-Na-Dene, were heavily affected by an additional and distinctive Beringian genetic input. In conclusion, most mtDNA variation (along the double-continent) stems from the first wave from Beringia, which followed the Pacific coastal route. This was accompanied or followed by a second inland migratory event, marked by haplogroups X2a and C4c, which affected all Amerindian groups of Northern North America. Much later, the ancestral A2a carriers spread from Alaska, undertaking both a westward migration to Asia and an eastward expansion into the circumpolar regions of Canada. Thus, the first American founders left the greatest genetic mark but the original maternal makeup of North American Natives was subsequently reshaped by additional streams of gene flow and local population dynamics, making a three-wave view too simplistic.

Microevolution, migration, and the population structure of five Amerindian populations from Nicaragua and Costa Rica

OBJECTIVE

This research examines the coevolution of languages and uniparental genetic marker (mitochondrial DNA [mtDNA] and nonrecombining Y-chromosome [NRY]) variation within five Lower Central American (Rama, Chorotega, Maléku, Zapatón-Huetar, and Abrojo-Guaymí) Amerindian groups. This pattern occurred since European contact.

METHODS

We examined mtDNA sequence variation from the hypervariable region 1 (HVS-1) and NRY genetic variation using short tandem repeat (STR) loci (DYS19, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, and DYS439) and NRY haplogroups (Q1a3a, Q1a3*, C3b, R1b1b2, E1b1, G2a2, and I) identified through single-nucleotide polymorphisms. Phylogenetic analysis included multidimensional scaling (MDS), heterozygosity versus rii , and analysis of molecular variance (AMOVA).

RESULTS

Eighteen mtDNA haplotypes were characterized in 131 participants with 94.6% of these assigned to the Amerindian mtDNA subclades, A2 and B2. The Amerindian NRY haplogroup, Q1a3a, was present in all five groups and ranged from 85% (Zapatón-Huetar) to 35% (Chorotega). Four populations (Rama, Chorotega, Zapatón-Huetar, and Abrojo-Guaymí) were also characterized by the presence of NRY haplogroup R1b1b2 indicative of western European admixture. Seventy NRY STR haplotypes were identified of which 69 (97%) were population specific. MDS plots demonstrated genetic similarities between Mesoamericans and northern Chibchan Amerindian populations, absent in mtDNA analyses, which is further supported by heterozygosity versus rii results.

CONCLUSIONS

We conclude that although these linguistically related populations in geographic proximity demonstrate a high degree of paternal genetic differentiation, recent demographic events have dramatically altered the paternal genetic structure of the regions Amerindian populations.

The genetic legacy of the pre-colonial period in contemporary Bolivians

Only a few genetic studies have been carried out to date in Bolivia. However, some of the most important (pre)historical enclaves of South America were located in these territories. Thus, the (sub)-Andean region of Bolivia was part of the Inca Empire, the largest state in Pre-Columbian America. We have genotyped the first hypervariable region (HVS-I) of 720 samples representing the main regions in Bolivia, and these data have been analyzed in the context of other pan-American samples (>19,000 HVS-I mtDNAs). Entire mtDNA genome sequencing was also undertaken on selected Native American lineages. Additionally, a panel of 46 Ancestry Informative Markers (AIMs) was genotyped in a sub-set of samples. The vast majority of the Bolivian mtDNAs (98.4%) were found to belong to the main Native American haplogroups (A: 14.3%, B: 52.6%, C: 21.9%, D: 9.6%), with little indication of sub-Saharan and/or European lineages; however, marked patterns of haplogroup frequencies between main regions exist (e.g. haplogroup B: Andean [71%], Sub-Andean [61%], Llanos [32%]). Analysis of entire genomes unraveled the phylogenetic characteristics of three Native haplogroups: the pan-American haplogroup B2b (originated ∼21.4 thousand years ago [kya]), A2ah (∼5.2 kya), and B2o (∼2.6 kya). The data suggest that B2b could have arisen in North California (an origin even in the north most region of the American continent cannot be disregarded), moved southward following the Pacific coastline and crossed Meso-America. Then, it most likely spread into South America following two routes: the Pacific path towards Peru and Bolivia (arriving here at about ∼15.2 kya), and the Amazonian route of Venezuela and Brazil southwards. In contrast to the mtDNA, Ancestry Informative Markers (AIMs) reveal a higher (although geographically variable) European introgression in Bolivians (25%). Bolivia shows a decreasing autosomal molecular diversity pattern along the longitudinal axis, from the Altiplano to the lowlands. Both autosomes and mtDNA revealed a low impact (1-2%) of a sub-Saharan component in Bolivians.

Arrival of Paleo-Indians to the southern cone of South America: new clues from mitogenomes

With analyses of entire mitogenomes, studies of Native American mitochondrial DNA (mtDNA) variation have entered the final phase of phylogenetic refinement: the dissection of the founding haplogroups into clades that arose in America during and after human arrival and spread. Ages and geographic distributions of these clades could provide novel clues on the colonization processes of the different regions of the double continent. As for the Southern Cone of South America, this approach has recently allowed the identification of two local clades (D1g and D1j) whose age estimates agree with the dating of the earliest archaeological sites in South America, indicating that Paleo-Indians might have reached that region from Beringia in less than 2000 years. In this study, we sequenced 46 mitogenomes belonging to two additional clades, termed B2i2 (former B2l) and C1b13, which were recently identified on the basis of mtDNA control-region data and whose geographical distributions appear to be restricted to Chile and Argentina. We confirm that their mutational motifs most likely arose in the Southern Cone region. However, the age estimate for B2i2 and C1b13 (11–13,000 years) appears to be younger than those of other local clades. The difference could reflect the different evolutionary origins of the distinct South American-specific sub-haplogroups, with some being already present, at different times and locations, at the very front of the expansion wave in South America, and others originating later in situ, when the tribalization process had already begun. A delayed origin of a few thousand years in one of the locally derived populations, possibly in the central part of Chile, would have limited the geographical and ethnic diffusion of B2i2 and explain the present-day occurrence that appears to be mainly confined to the Tehuelche and Araucanian-speaking groups.