Genetic variation and population structure in native Americans

Dense SNP-based analyses complement forensic anthropology biogeographical ancestry assessments

Identification of unidentified human remains (UHRs) is crucial yet challenging, especially with traditional forensic techniques. Forensic anthropological examinations can yield ancestry estimations; however, the utility of these estimates is limited by the data points that can be collected from partial remains, complexities of admixture, and variation of phenotypic expression due to environmental effects. While it is generally known that anthropological estimates can be imprecise, the performance of these methods has not been studied at scale. Genome-wide SNP testing is an orthogonal approach for estimating ancestry and offers a unique opportunity to measure the magnitude of anthropological ancestry misattribution. Genomic ancestry inference leverages principal component analysis (PCA) and model-based clustering approaches. This study compares anthropologically determined ancestry with those estimated using genome-wide SNP markers. A dataset of 611 UHR samples with publicly available ancestry assessments from National Missing and Unidentified Persons System (NamUs) was analyzed. The genetic ancestry approach, validated against reference population samples, offers robust ancestry calculations for major population groups. Inconsistency between anthropological and genomic ancestry assignments were observed, particularly for admixed populations. Although forensic anthropological examinations remain valuable, their limitations emphasize the need for refinement and enhancement through the augmentation of SNP-based analyses. Further validation studies are crucial to define the uncertainty associated with both anthropological and genome-based ancestry estimates to resolve cases and aid law enforcement investigations. Additionally, current policy and practices for reporting ancestry for UHRs should be revisited to reduce potential misinformation.

Population expansion, larger, and more homogeneous native American ancestry among Mexican mestizo populations based on 10 X-chromosome STR loci (X-STR decaplex system)

OBJECTIVE

To evaluate the genetic diversity, admixture, genetic relationships, and sex-biased demographic processes in Mexican Mestizo (admixed) populations based on 10 X-chromosome STRs (X-STRs).

METHODS

We analyzed the X-STRs Decaplex system in 104 Mexican Native Americans to obtain the ancestral reference needed to complete the demographic analyses above mentioned. We included reported Iberian and Latin American (admixed) populations from Central and South America, as well as datasets from Mexican Mestizos based on Y-linked STRs (Y-STRs), autosomal STRs (A-STRs), and mtDNA.

RESULTS

Higher X-linked Native American ancestry was observed among Latin American populations regarding that reported from A-STRs and Y-STRs. The interpopulation differentiation based on ancestry among Mexican Mestizos diminished according to the inheritance pattern: Y-STRs (highest), A-STRs, X-STRs, and mtDNA (lowest). This finding is related to the peculiar admixture process that occurred during and after the Spanish Conquest of Mexico (and most of Latin America), involving a large number of Spanish men (Y-chromosomes) with a lesser proportion of X-chromosomes than autosomes; besides to the limited number of Spanish women (XX) arrived in the Americas in subsequent and shorter periods. Population expansion was detected in Mexican Mestizos from all the country, except those from the southeast region characterized by elevated indigenous ancestry, marginalization, and poorness.

CONCLUSIONS

Population growth was detected in most Mexican Mestizos, besides more homogeneous and larger Native American ancestry based on X-linked inheritance than that based on autosomal STRs and Y-STRs.

MULT: An allometric body mass index (ABMI) reference to assess nutritional status of multiethnic children and adolescents

OBJECTIVES

To develop an allometric body mass index (ABMI) reference that adjusts the weight in relation to height, taking into account the changes during development (MULT ABMI reference), and to compare it with international BMI references.

METHODS

The MULT ABMI reference was constructed through the LMS method, calculated with 65 644 ABMI observations of 17 447 subjects aged 5-22 years, from the United Kingdom, Ethiopia, India, Peru, Vietnam, Portugal, and Brazil. The M, S, and L curves of the MULT ABMI reference were compared with the curves of the MULT, World Health Organization (WHO), Centers for Disease Control and Prevention (CDC), International Obesity Task Force (IOTF), and Dutch Growth Study (DUTCH).

RESULTS

The greater differences in the M curve between MULT ABMI and WHO, CDC, IOTF, DUTCH, and MULT BMI references were around puberty (138 to 150 months for boys; 114 to 132 for girls). MULT ABMI presented S values similar to IOTF and DUTCH BMI references for boys 60 to 114 months and then became higher, approaching the MULT BMI S values from 198 to 240 months. For girls the MULT ABMI S values were close to the IOTF, CDC, and DUTCH from 60 to 110 months, and then became higher, approaching the MULT BMI S values until 240 months.

CONCLUSION

MULT ABMI presented an advantage in comparison to the existing BMI references because it takes into account the growth changes during puberty and is a new option to assess the nutritional status of multiethnic populations.

Genetic associations with disease in populations with Indigenous American ancestries

The genetic architecture of complex diseases affecting populations with Indigenous American ancestries is poorly understood due to their underrepresentation in genomics studies. While most of the genetic diversity associated with disease trait variation is shared among worldwide populations, a fraction of this component is expected to be unique to each continental group, including Indigenous Americans. Here, I describe the current state of knowledge from genome-wide association studies on Indigenous populations, as well as non-Indigenous populations with partial Indigenous ancestries from the American continent, focusing on disease susceptibility and anthropometric traits. While some studies identified risk alleles unique to Indigenous populations, their effects on trait variation are mostly small. I suggest that the associations rendered by many inter-population studies are probably inflated due to the absence of socio-cultural-economic covariates in the association models. I encourage the inclusion of admixed individuals in future GWAS studies to control for inter-ancestry differences in environmental factors. I suggest that some complex diseases might have arisen as trade-off costs of adaptations to past evolutionary selective pressures. Finally, I discuss how expanding panels with Indigenous ancestries in GWAS studies is key to accurately assess genetic risk in populations from the American continent, thus decreasing global health disparities.

Stature estimation equations for modern American Indians in the American Southwest

Stature estimation is a core component to the biological profile in forensic anthropology casework. Here we provide mathematical equations for estimating stature for contemporary American Indians (AI), which currently are lacking in forensic anthropology. Drawing on postmortem computed tomography data from the New Mexico Decedent Image Database we regressed cadaveric length on four long bone length measures of the tibia, femur, and humerus to produce 11 combinations of models. Separate regression models were calculated for the entire pooled sample, by sex, broad AI language groups, and age + sex subsamples and compared. Sex-specific models were statistically better than general models, which were more accurate than language group and age + sex models. Equations were created for general and sex-specific models. Application to an independent test sample demonstrates the equations are accurate for stature estimation with overestimates of less than 1 cm. The equations provide similar levels of precision to stature estimation programs like the FORDISC 3.0 module and other stature equations in the literature. We provide recommendations for equation use in casework based on our results. These equations are the first for estimating stature in contemporary AI. This paper demonstrates the appropriateness of these newly created stature equations for use in New Mexico and the surrounding region.

IKZF1PLUS alterations contribute to outcome disparities in Hispanic/Latino children with B-lymphoblastic leukemia

BACKGROUND

Compared to other ethnicities, Hispanics/Latinos (H/L) have a high incidence of acute lymphoblastic leukemia (ALL), enrichment of unfavorable ALL genetic subtypes, and worse outcomes, even after correcting for socioeconomic factors. We previously demonstrated increased incidence of the high-risk genetic drivers IKZF1 deletion and IGH::CRLF2 rearrangement in H/L compared to non-H/L children with B-ALL. Here in an expanded pediatric cohort, we sought to identify novel genetic drivers and secondary genetic alterations in B-ALL associated with H/L ethnicity.

PROCEDURE

Comprehensive clinicopathologic data from patients with B-ALL treated from 2016 to 2020 were analyzed. Subtype was determined from karyotype, fluorescence in situ hybridization (FISH), chromosome microarray (CMA), and our next-generation sequencing (NGS) panel (OncoKids). Non-driver genetic variants were also examined. p-Values less than .05 (Fisher's exact test) were considered significant.

RESULTS

Among patients with B-ALL at diagnosis (n = 273), H/L patients (189, 69.2%) were older (p = .018), more likely to present with CNS2 or CNS3 disease (p = .004), and NCI high-risk ALL (p = .014) compared to non-H/L patients. Higher incidence of IGH::CRLF2 rearrangement (B-ALL, BCR::ABL1-like, unfavorable; p = .016) and lower incidence of ETV6::RUNX1 rearrangement (favorable, p = .02) were also associated with H/L ethnicity. Among secondary (non-subtype-defining) genetic variants, B-ALL in H/L was associated with IKFZ1 deletion alone (p = .001) or with IGH::CRLF2 rearrangement (p = .003). The IKZF1PLUS profile (IKZF1 deletion plus CDKN2A/2Bdel, PAX5del, or P2RY8::CRLF2 rearrangement without DUX4 rearrangement) was identified as a novel high-risk feature enriched in H/L patients (p = .001).

CONCLUSIONS

Our study shows enrichment of high-risk genetic variants in H/L B-ALL and raises consideration for novel therapeutic targets.

The most frequent HLA alleles around the world: A fundamental synopsis

A comprehensive knowledge of human leukocyte antigen (HLA) molecular variation worldwide is essential in human population genetics research and disease association studies and is also indispensable for clinical applications such as allogeneic hematopoietic cell transplantation, where ensuring HLA compatibility between donors and recipients is paramount. Enormous progress has been made in this field thanks to several decades of HLA population studies allowing the development of helpful databases and bioinformatics tools. However, it is still difficult to appraise the global HLA population diversity in a synthetic way. We thus introduce here a novel approach, based on approximately 2000 data sets, to assess this complexity by providing a fundamental synopsis of the most frequent HLA alleles observed in different regions of the world. This new knowledge will be useful not only as a fundamental reference for basic research, but also as an efficient guide for clinicians working in the field of transplantation.

Comprehensive landscape of non-CODIS STRs in global populations provides new insights into challenging DNA profiles

The worldwide implementation of short tandem repeats (STR) profiles in forensic genetics necessitated establishing and expanding the CODIS core loci set to facilitated efficient data management and exchange. Currently, the mainstay CODIS STRs are adopted in most general-purpose forensic kits. However, relying solely on these loci failed to yield satisfactory results for challenging tasks, such as bio-geographical ancestry inference, complex DNA mixture profile interpretation, and distant kinship analysis. In this context, non-CODIS STRs are potent supplements to enhance the systematic discriminating power, particularly when combined with the high-throughput next-generation sequencing (NGS) technique. Nevertheless, comprehensive evaluation on non-CODIS STRs in diverse populations was scarce, hindering their further application in routine caseworks. To address this gap, we investigated genetic variations of 178 historically available non-CODIS STRs from ethnolinguistically different worldwide populations and studied their characteristics and forensic potentials via high-coverage whole genome sequencing (WGS) data. Initially, we delineated the genomic properties of these non-CODIS markers through sequence searching, repeat structure scanning, and manual inspection. Subsequent population genetics analysis suggested that these non-CODIS STRs had comparable polymorphism levels and forensic utility to CODIS STRs. Furthermore, we constructed a theoretical next-generation sequencing (NGS) panel comprising 108 STRs (20 CODIS STRs and 88 non-CODIS STRs), and evaluated its performance in inferring bio-geographical ancestry origins, deconvoluting complex DNA mixtures, and differentiating distant kinships using real and simulated datasets. Our findings demonstrated that incorporating supplementary non-CODIS STRs enabled the extrapolation of multidimensional information from a single STR profile, thereby facilitating the analysis of challenging forensic tasks. In conclusion, this study presents an extensive genomic landscape of forensic non-CODIS STRs among global populations, and emphasized the imperative inclusion of additional polymorphic non-CODIS STRs in future NGS-based forensic systems.

Clumppling: cluster matching and permutation program with integer linear programming

MOTIVATION

In the mixed-membership unsupervised clustering analyses commonly used in population genetics, multiple replicate data analyses can differ in their clustering solutions. Combinatorial algorithms assist in aligning clustering outputs from multiple replicates so that clustering solutions can be interpreted and combined across replicates. Although several algorithms have been introduced, challenges exist in achieving optimal alignments and performing alignments in reasonable computation time.

RESULTS

We present Clumppling, a method for aligning replicate solutions in mixed-membership unsupervised clustering. The method uses integer linear programming for finding optimal alignments, embedding the cluster alignment problem in standard combinatorial optimization frameworks. In example analyses, we find that it achieves solutions with preferred values of a desired objective function relative to those achieved by Pong and that it proceeds with less computation time than Clumpak. It is also the first method to permit alignments across replicates with multiple arbitrary values of the number of clusters K.

AVAILABILITY AND IMPLEMENTATION

Clumppling is available at https://github.com/PopGenClustering/Clumppling.

Germline mutations in pediatric cancer cohort with mixed-ancestry Mexicans

BACKGROUND

Childhood cancer is one of the primary causes of disease-related death in 5- to 14-year-old children and currently no prevention strategies exist to reduce the incidence of this disease. Childhood cancer has a larger hereditary component compared with cancer in adults. Few genetic studies have been conducted on children with cancer. Additionally, Latin American populations are underrepresented in genomic studies compared with other populations. Therefore, the aim of this study is to analyze germline mutations in a group of mixed-ancestry Mexican pediatric patients with solid and hematological cancers.

METHODS

We analyzed genetic variants from 40 Mexican childhood cancer patients and their relatives. DNA from saliva or blood samples was used for whole-exome sequencing. All variants were identified following GATK best practices.

RESULTS

We found that six patients (15%) were carriers of germline mutations in CDKN2A, CHEK2, DICER1, FANCA, MSH6, MUTYH, NF1, and SBDS cancer predisposition genes, and additional new variants predicted to be deleterious by in silico algorithms. A population genetics analysis detected five components consistent with the demographic models assumed for modern mixed-ancestry Mexicans.

CONCLUSIONS

This report identifies potential genetic risk factors and provides a better understanding of the underlying mechanisms of childhood cancer in this population.

Variants of IFNL4 Gene in Amazonian and Northern Brazilian Populations

Since the discovery of the polymorphic nature of the IFNL4 gene, its variants have been investigated and associated with several viral diseases, with an emphasis on hepatitis C. However, the impacts of these variants on mixed-race and native populations in the northern region of Brazil are scarce. We investigated three variants of the IFNL4 gene in populations from this location, which were among the 14 most frequent variants in worldwide populations, and compared the frequencies obtained to populational data from the 1000 Genomes Project, gnomAD and ABraOM databases. Our results demonstrate that mixed-race and native populations from the northern region of Brazil present frequencies like those of European and Asian groups for the rs74597329 and rs11322783 variants, and like all populations presented for the rs4803221 variant. These data reinforce the role of world populations in shaping the genetic profile of Brazilian populations, indicate patterns of illness according to the expressed genotype, and infer an individual predisposition to certain diseases.

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.

Craniovascular variation in four late Holocene human samples from southern South America

Craniovascular traits in the endocranium (traces of middle meningeal vessels and dural venous sinuses, emissary foramina) provide evidence of vascular anatomy in osteological samples. We investigate the craniovascular variation in four South American samples and the effect of artificial cranial modifications (ACM). CT scans of human adult crania from four archeological samples from southern South America (including skulls with ACM) are used for the analyses. The craniovascular features in the four samples are described, skulls with and without ACM are compared, and additionally, South Americans are compared to a previously analyzed sample of Europeans. Of the four South American samples, the Southern Patagonian differs the most, showing the most distinct cranial dimensions, no ACM, and larger diameters of the emissary foramina. Unlike previous studies, we did not find any major differences in craniovascular features between modified and non-modified skulls, except that the skulls with ACM present somewhat smaller foramina. South Americans significantly differed from Europeans, especially in the anteroposterior dominance of the middle meningeal artery, in the pattern of sinus confluence, in the occurrence of enlarged occipito-marginal sinuses, and in foramina frequencies and diameters. Craniovascular morphology is not affected by the cranial size, even in skulls with ACM, indicating a minor or null influence of structural topological factors. Concerning the samples from distinct geographic and climatic environments, it must be evaluated whether the craniovascular morphogenesis might be partially influenced by specific functions possibly associated with thermoregulation, intracranial pressure, and the maintenance of intracranial homeostasis.

An HLA map of the world: A comparison of HLA frequencies in 200 worldwide populations reveals diverse patterns for class I and class II

HLA frequencies show widespread variation across human populations. Demographic factors as well as selection are thought to have shaped HLA variation across continents. In this study, a worldwide comparison of HLA class I and class II diversity was carried out. Multidimensional scaling techniques were applied to 50 HLA-A and HLA-B (class I) as well as 13 HLA-DRB1 (class II) first-field frequencies in 200 populations from all continents. Our results confirm a strong effect of geography on the distribution of HLA class I allele groups, with principal coordinates analysis closely resembling geographical location of populations, especially those of Africa-Eurasia. Conversely, class II frequencies stratify populations along a continuum of differentiation less clearly correlated to actual geographic location. Double clustering analysis revealed finer intra-continental sub-clusters (e.g., Northern and Western Europe vs. South East Europe, North Africa and Southwest Asia; South and East Africa vs. West Africa), and HLA allele group patterns characteristic of these clusters. Ancient (Austronesian expansion) and more recent (Romani people in Europe) migrations, as well as extreme differentiation (Taiwan indigenous peoples, Native Americans), and interregional gene flow (Sámi, Egyptians) are also reflected by the results. Barrier analysis comparing D_ST and geographic location identified genetic discontinuities caused by natural barriers or human behavior explaining inter and intra-continental HLA borders for class I and class II. Overall, a progressive reduction in HLA diversity from African to Oceanian and Native American populations is noted. This analysis of HLA frequencies in a unique set of worldwide populations confirms previous findings on the remarkable similarity of class I frequencies to geography, but also shows a more complex development for class II, with implications for both human evolutionary studies and biomedical research.

Thyroidal and Extrathyroidal Requirements for Iodine and Selenium: A Combined Evolutionary and (Patho)Physiological Approach

Iodide is an antioxidant, oxidant and thyroid hormone constituent. Selenoproteins are needed for triiodothyronine synthesis, its deactivation and iodine release. They also protect thyroidal and extrathyroidal tissues from hydrogen peroxide used in the 'peroxidase partner system'. This system produces thyroid hormone and reactive iodine in exocrine glands to kill microbes. Exocrine glands recycle iodine and with high urinary clearance require constant dietary supply, unlike the thyroid. Disbalanced iodine-selenium explains relations between thyroid autoimmune disease (TAD) and cancer of thyroid and exocrine organs, notably stomach, breast, and prostate. Seafood is iodine unconstrained, but selenium constrained. Terrestrial food contains little iodine while selenium ranges from highly deficient to highly toxic. Iodine vs. TAD is U-shaped, but only low selenium relates to TAD. Oxidative stress from low selenium, and infection from disbalanced iodine-selenium, may generate cancer of thyroid and exocrine glands. Traditional Japanese diet resembles our ancient seashore-based diet and relates to aforementioned diseases. Adequate iodine might be in the milligram range but is toxic at low selenium. Optimal selenoprotein-P at 105 µg selenium/day agrees with Japanese intakes. Selenium upper limit may remain at 300-400 µg/day. Seafood combines iodine, selenium and other critical nutrients. It brings us back to the seashore diet that made us what we currently still are.

A genomic perspective on South American human history

It has generally been accepted that the current indigenous peoples of the Americas are derived from ancestors from northeastern Asia. The latter were believed to have spread into the American continent by the end of the Last Glacial Maximum. In this sense, a joint and in-depth study of the earliest settlement of East Asia and the Americas is required to elucidate these events accurately. The first Americans underwent an adaptation process to the Americas' vast environmental diversity, mediated by biological and cultural evolution and niche construction, resulting in enormous cultural diversity, a wealth of domesticated species, and extensive landscape modifications. Afterward, in the Late Holocene, the advent of intensive agricultural food production systems, sedentism, and climate change significantly reshaped genetic and cultural diversity across the continent, particularly in the Andes and Amazonia. Furthermore, starting around the end of the 15th century, European colonization resulted in massive extermination of indigenous peoples and extensive admixture. Thus, the present review aims to create a comprehensive picture of the main events involved in the formation of contemporary South American indigenous populations and the dynamics responsible for shaping their genetic diversity by integrating current genetic data with evidence from archeology, linguistics and other disciplines.

Disentangling Signatures of Selection Before and After European Colonization in Latin Americans

Throughout human evolutionary history, large-scale migrations have led to intermixing (i.e., admixture) between previously separated human groups. Although classical and recent work have shown that studying admixture can yield novel historical insights, the extent to which this process contributed to adaptation remains underexplored. Here, we introduce a novel statistical model, specific to admixed populations, that identifies loci under selection while determining whether the selection likely occurred post-admixture or prior to admixture in one of the ancestral source populations. Through extensive simulations, we show that this method is able to detect selection, even in recently formed admixed populations, and to accurately differentiate between selection occurring in the ancestral or admixed population. We apply this method to genome-wide SNP data of ∼4,000 individuals in five admixed Latin American cohorts from Brazil, Chile, Colombia, Mexico, and Peru. Our approach replicates previous reports of selection in the human leukocyte antigen region that are consistent with selection post-admixture. We also report novel signals of selection in genomic regions spanning 47 genes, reinforcing many of these signals with an alternative, commonly used local-ancestry-inference approach. These signals include several genes involved in immunity, which may reflect responses to endemic pathogens of the Americas and to the challenge of infectious disease brought by European contact. In addition, some of the strongest signals inferred to be under selection in the Native American ancestral groups of modern Latin Americans overlap with genes implicated in energy metabolism phenotypes, plausibly reflecting adaptations to novel dietary sources available in the Americas.

The Genomic Profile Associated with Risk of Severe Forms of COVID-19 in Amazonian Native American Populations

Genetic factors associated with COVID-19 disease outcomes are poorly understood. This study aimed to associate genetic variants in the SLC6A20, LZTFL1, CCR9, FYCO1, CXCR6, XCR1, and ABO genes with the risk of severe forms of COVID-19 in Amazonian Native Americans, and to compare the frequencies with continental populations. The study population was composed of 64 Amerindians from the Amazon region of northern Brazil. The difference in frequencies between the populations was analyzed using Fisher’s exact test, and the results were significant when p ≤ 0.05. We investigated 64 polymorphisms in 7 genes; we studied 47 genetic variants that were new or had impact predictions of high, moderate, or modifier. We identified 15 polymorphisms with moderate impact prediction in 4 genes (ABO, CXCR6, FYCO1, and SLC6A20). Among the variants analyzed, 18 showed significant differences in allele frequency in the NAM population when compared to others. We reported two new genetic variants with modifier impact in the Amazonian population that could be studied to validate the possible associations with COVID-19 outcomes. The genomic profile of Amazonian Native Americans may be associated with protection from severe forms of COVID-19. This work provides genomic data that may help forthcoming studies to improve COVID-19 outcomes.

Characterization of PCLO Gene in Amazonian Native American Populations

Genetic variations in PCLO have been associated with different pathologies in global literature, but there are no data regarding this gene in Native American populations. The Amazonian Native American populations have lower genetic diversity and are more different from other continental groups. We investigated 18 genetic variants in the PCLO gene in Amazonian indigenous and compared our results with the ones found in global populations, which were publicly available in the 1000 Genomes Project, gnmAD and ABraOM databases. The results demonstrated that the variants of the PCLO, especially rs17156844, rs550369696, rs61741659 and rs2877, have a significantly higher frequency in Amerindian populations in comparison with other continental populations. These data outline the singular genetic profile of the Native American population from the Brazilian Amazon region.

Association of Genetic Ancestry With the Molecular Subtypes and Prognosis of Childhood Acute Lymphoblastic Leukemia

IMPORTANCE

Racial and ethnic disparities persist in the incidence and treatment outcomes of childhood acute lymphoblastic leukemia (ALL). However, there is a paucity of data describing the genetic basis of these disparities, especially in association with modern ALL molecular taxonomy and in the context of contemporary treatment regimens.

OBJECTIVE

To evaluate the association of genetic ancestry with childhood ALL molecular subtypes and outcomes of modern ALL therapy.

DESIGN, SETTING, AND PARTICIPANTS

This multinational, multicenter genetic association study was conducted from March 1, 2000, to November 20, 2020, among 2428 children and adolescents with ALL enrolled in frontline trials from the United States, South East Asia (Singapore and Malaysia), and Latin America (Guatemala), representing diverse populations of European, African, Native American, East Asian, and South Asian descent. Statistical analysis was conducted from February 3, 2020, to April 19, 2021.

MAIN OUTCOMES AND MEASURES

Molecular subtypes of ALL and genetic ancestry were comprehensively characterized by performing RNA sequencing. Associations of genetic ancestries with ALL molecular subtypes and treatment outcomes were then evaluated.

RESULTS

Among the participants in the study, 1340 of 2318 (57.8%) were male, and the mean (SD) age was 7.8 (5.3) years. Of 21 ALL subtypes identified, 8 were associated with ancestry. East Asian ancestry was positively associated with the frequency of somatic DUX4 (odds ratio [OR], 1.30 [95% CI, 1.16-1.45]; P < .001) and ZNF384 (OR, 1.40 [95% CI, 1.18-1.66]; P < .001) gene rearrangements and negatively associated with BCR-ABL1-like ALL (OR, 0.79 [95% CI, 0.66-0.92]; P = .002) and T-cell ALL (OR, 0.80 [95% CI, 0.71-0.90]; P < .001). By contrast, occurrence of CRLF2 rearrangements was associated with Native American ancestry (OR, 1.48 [95% CI, 1.29-1.69]; P < .001). When the percentage of Native American ancestry increased, ETV6-RUNX1 fusion became less frequent (OR, 0.80 [95% CI, 0.70-0.91]; P < .001), with the opposite trend observed for ETV6-RUNX1-like ALL. There was a marked preponderance of T-cell ALL in children of African descent compared with those with a high percentage of Native American ancestry (African: OR, 1.22 [95% CI, 1.07-1.37]; P = .003; Native American: OR, 0.53 [95% CI, 0.40-0.67]; P < .001). African ancestry was also positively associated with the prevalence of TCF3-PBX1 (OR, 1.49 [95% CI, 1.25-1.76]; P < .001) and negatively associated with DUX4 rearrangements (OR, 0.70 [95% CI, 0.48-0.93]; P = .01) and hyperdiploidy (OR, 0.77 [95% CI, 0.68-0.86]; P < .001). African and Native American ancestries as continuous variables were both associated with poorer event-free survival (for every 25% increase in ancestry: hazard ratio [HR], 1.2; 95% CI, 1.1-1.4; P = .001 for African ancestry; HR, 1.3; 95% CI, 1.0-1.6; P = .04 for Native American ancestry) and overall survival (for every 25% increase in ancestry: HR, 1.2; 95% CI, 1.1-1.5; P = .01 for African ancestry; HR, 1.4; 95% CI, 1.0-1.8; P = .03 for Native American ancestry). Even after adjusting for biological subtypes and clinical features, Native American and African ancestries remained associated with poor prognosis.

CONCLUSIONS AND RELEVANCE

This study suggests that ALL molecular subtypes and prognosis are associated with genetic ancestry, potentially pointing to a genetic basis for some of the racial and ethnic disparities in ALL. Therefore, molecular subtype-driven treatment individualization is needed to help address racial and ethnic gaps in outcomes.

Remarkably Low KIR and HLA Diversity in Amerindians Reveals Signatures of Strong Purifying Selection Shaping the Centromeric KIR Region

The killer-cell immunoglobulin-like receptors (KIR) recognize human leukocyte antigen (HLA) molecules to regulate the cytotoxic and inflammatory responses of natural killer cells. KIR genes are encoded by a rapidly evolving gene family on chromosome 19 and present an unusual variation of presence and absence of genes and high allelic diversity. Although many studies have associated KIR polymorphism with susceptibility to several diseases over the last decades, the high-resolution allele-level haplotypes have only recently started to be described in populations. Here, we use a highly innovative custom next-generation sequencing method that provides a state-of-art characterization of KIR and HLA diversity in 706 individuals from eight unique South American populations: five Amerindian populations from Brazil (three Guarani and two Kaingang); one Amerindian population from Paraguay (Aché); and two urban populations from Southern Brazil (European and Japanese descendants from Curitiba). For the first time, we describe complete high-resolution KIR haplotypes in South American populations, exploring copy number, linkage disequilibrium, and KIR-HLA interactions. We show that all Amerindians analyzed to date exhibit the lowest numbers of KIR-HLA interactions among all described worldwide populations, and that 83-97% of their KIR-HLA interactions rely on a few HLA-C molecules. Using multiple approaches, we found signatures of strong purifying selection on the KIR centromeric region, which codes for the strongest NK cell educator receptors, possibly driven by the limited HLA diversity in these populations. Our study expands the current knowledge of KIR genetic diversity in populations to understand KIR-HLA coevolution and its impact on human health and survival.

The Heart of Silk Road "Xinjiang," Its Genetic Portray, and Forensic Parameters Inferred From Autosomal STRs

The Xinjiang Uyghur Autonomous Region of China (XUARC) harbors almost 50 ethnic groups including the Uyghur (UGR: 45.84%), Han (HAN: 40.48%), Kazakh (KZK: 6.50%), Hui (HUI: 4.51%), Kyrgyz (KGZ: 0.86%), Mongol (MGL: 0.81%), Manchu (MCH: 0.11%), and Uzbek (UZK: 0.066%), which make it one of the most colorful regions with abundant cultural and genetic diversities. In our previous study, we established allelic frequency databases for 14 autosomal short tandem repeats (STRs) for four minority populations from XUARC (MCH, KGZ, MGL, and UZK) using the AmpFlSTR® Identifiler PCR Amplification Kit. In this study, we genotyped 2,121 samples using the GoldenEye™ 20A Kit (Beijing PeopleSpot Inc., Beijing, China) amplifying 19 autosomal STR loci for four major ethnic groups (UGR, HAN, KZK, and HUI). These groups make up 97.33% of the total XUARC population. The total number of alleles for all the 19 STRs in these populations ranged from 232 (HAN) to 224 (KZK). We did not observe any departures from the Hardy-Weinberg equilibrium (HWE) in these populations after sequential Bonferroni correction. We did find minimal departure from linkage equilibrium (LE) for a small number of pairwise combinations of loci. The match probabilities for the different populations ranged from 1 in 1.66 × 1023 (HAN) to 6.05 × 1024 (HUI), the combined power of exclusion ranged from 0.999 999 988 (HUI) to 0.999 999 993 (UGR), and the combined power of discrimination ranged from 0.999 999 999 999 999 999 999 983 (HAN) to 0.999 999 999 999 999 999 999 997 (UGR). Genetic distances, principal component analysis (PCA), STRUCTURE analysis, and the phylogenetic tree showed that genetic affinity among studied populations is consistent with linguistic, ethnic, and geographical classifications.

HLA-G 3'UTR haplotype frequencies in highland and lowland South Native American populations

Human Leukocyte Antigen (HLA)-G participates in several biological processes, including reproduction, vascular remodeling, immune tolerance, and hypoxia response. HLA-G is a potential candidate gene for high altitude adaptation since its expression is modulated in both micro and macro environment under hypoxia and constant cellular stress. Besides the promoter region, the HLA-G 3'untranslated region (UTR) influences HLA-G expression patterns through several post-transcriptional mechanisms. Currently, the 3'UTR genetic diversity in terms of altitude adaptation of Native American populations is still unexplored, particularly at high altitude ecoregions. Here, we evaluated 288 Native Americans from 9 communities located in the Andes [highland (HL); ≥2,500 m (range = 2,838-4,433 m)] and 8 populations located in lowland (LL) regions [<2,500 m (range = 80-431 m); Amazonian tropical forest, Brazilian central plateau, and Chaco] of South America. In total, nine polymorphic sites and ten haplotypes were observed. The most frequent haplotypes (UTR-1, UTR-2, and UTR-3) accounted for ∼ 77% of haplotypes found in LL, while in the HL, the same haplotypes reach ∼ 93%. Also, a remarkable high frequency of putative ancestral UTR-5 haplotype was observed in LL (21.5%), while in HL UTR-2 reaches up to 47%. Further, UTR-2 frequency positively correlates with altitude-related variables, while a negative correlation for UTR-5 was observed. From an evolutionary perspective, we observed a tendency towards balancing selection in HL and LL populations thus suggesting that haplotypes of ancient and more derived alleles may have been co-opted for relatively recent adaptations such as those experienced by modern humans in the highland and lowland of South America. We also discuss how long-term balancing selection can be a reservoir of genetic variants that can be positively selected. Finally, our study provides some pieces of evidence that HLA-G 3'UTR haplotypes may have contributed to high altitude adaptation in the Andes.

Genetic differentiation between and within Northern Native American language groups: an argument for the expansion of the Native American CODIS database

The National Research Council recommends that genetic differentiation among subgroups of ethnic samples be lower than 3% of the total genetic differentiation within the ethnic sample to be used for estimating reliable random match probabilities for forensic use. Native American samples in the United States' Combined DNA Index System (CODIS) database represent four language families: Algonquian, Na-Dene, Eskimo-Aleut, and Salishan. However, a minimum of 27 Native American language families exists in the US, not including language isolates. Our goal was to ascertain whether genetic differences are correlated with language groupings and, if so, whether additional language families would provide a more accurate representation of current genetic diversity among tribal populations. The 21 short tandem repeat (STR) loci included in the Globalfiler® PCR Amplification Kit were used to characterize six indigenous language families, including three of the four represented in the CODIS database (i.e. Algonquian, Na-Dene, and Eskimo-Aleut), and two language isolates (Miwok and Seri) using major population genetic diversity metrics such as F statistics and Bayesian clustering analysis of genotype frequencies. Most of the genetic variation (97%) was found to be within language families instead of among them (3%). In contrast, when only the three of the four language families represented in both the CODIS database and the present study were considered, 4% of the genetic variation occurred among the language groups. Bayesian clustering resulted in a maximum posterior probability indicating three genetically distinct groups among the eight language families and isolates: (1) Eskimo, (2) Seri, and (3) all other language groups and isolates, thus confirming genetic subdivision among subgroups of the CODIS Native American database. This genetic structure indicates the need for an increased number of Native American populations based on language affiliation in the CODIS database as well as more robust sample sets for those language families. Supplemental data for this article is available online at https://doi.org/10.1080/20961790.2021.1963088 .

Performance of the CKD-EPI and MDRD equations for estimating glomerular filtration rate: a systematic review of Latin American studies

BACKGROUND

The most-used equations for estimating the glomerular filtration rate (GFR) are the CKD Epidemiology Collaboration (CKD-EPI) and Modification of Diet in Renal Disease (MDRD) equations. However, it is unclear which of these shows better performance in Latin America.

OBJECTIVE

To assess the performance of two equations for estimated GFR (eGFR) in Latin American countries.

DESIGN AND SETTING

Systematic review and meta-analysis in Latin American countries.

METHODS

We searched in three databases to identify studies that reported eGFR using both equations and compared them with measured GFR (mGFR) using exogenous filtration markers, among adults in Latin American countries. We performed meta-analyses on P30, bias (using mean difference [MD] and 95% confidence intervals [95% CI]), sensitivity and specificity; and evaluated the certainty of evidence using the GRADE methodology.

RESULTS

We included 12 papers, and meta-analyzed six (five from Brazil and one from Mexico). Meta-analyses that compared CKD-EPI using creatinine measured with calibration traceable to isotope dilution mass spectrometry (CKD-EPI-Cr IDMS) and using MDRD-4 IDMS did not show differences in bias (MD: 0.55 ml/min/1.73m2; 95% CI: -3.34 to 4.43), P30 (MD: 4%; 95% CI: -2% to 11%), sensitivity (76% and 75%) and specificity (91% and 89%), with very low certainty of evidence for bias and P30, and low certainty of evidence for sensitivity and specificity.

CONCLUSION

We found that the performances of CKD-EPI-Cr IDMS and MDRD-4 IDMS did not differ significantly. However, since most of the meta-analyzed studies were from Brazil, the results cannot be extrapolated to other Latin American countries.

REGISTRATION

PROSPERO (CRD42019123434) - https://www.crd.york.ac.uk/prospero/display_record.php?ID=CRD42019123434.

How natural selection shapes genetic differentiation in the MHC region: A case study with Native Americans

The Human Leukocyte Antigen (HLA) loci are extremely well documented targets of balancing selection, yet few studies have explored how selection affects population differentiation at these loci. In the present study we investigate genetic differentiation at HLA genes by comparing differentiation at microsatellites distributed genomewide to those in the MHC region. Our study uses a sample of 494 individuals from 30 human populations, 28 of which are Native Americans, all of whom were typed for genomewide and MHC region microsatellites. We find greater differentiation in the MHC than in the remainder of the genome (FST-MHC = 0.130 and FST-Genomic = 0.087), and use a permutation approach to show that this difference is statistically significant, and not accounted for by confounding factors. This finding lies in the opposite direction to the expectation that balancing selection reduces population differentiation. We interpret our findings as evidence that selection favors different sets of alleles in distinct localities, leading to increased differentiation. Thus, balancing selection at HLA genes simultaneously increases intra-population polymorphism and inter-population differentiation in Native Americans.

Exchanging fluids The sociocultural implications of microbial, cultural, and ethnic admixture in Latin America

Knowledge of evolutionary influences on patterns of human mating, social interactions, and differential health is increasing, yet these insights have rarely been applied to historical analyses of human population dynamics. The genetic and evolutionary forces behind biases in interethnic mating and in the health of individuals of different ethnic groups in Latin America and the Caribbean since the European colonization of America are still largely ignored. We discuss how historical and contemporary sociocultural interactions and practices are strongly influenced by population-level evolutionary forces. Specifically, we discuss the historical implications of functional (de facto) polygyny, sex-biased admixture, and assortative mating in Latin America. We propose that these three evolutionary mechanisms influenced mating patterns, shaping the genetic and cultural landscape across Latin America and the Caribbean. Further, we discuss how genetic differences between the original populations that migrated at different times into Latin America contributed to their accommodation to and survival in the different local ecologies and interethnic interactions. Relevant medical and social implications follow from the genetic and cultural changes reviewed.

Population structure and relatedness estimates in a Mexican sample

Population stratification (PS) is a confounding factor in genome-wide association studies (GWASs) and also an interesting process itself. Latin American populations have mixed genetic ancestry, which may account for PS. We have analyzed the relatedness, by means of the identity-by-descent (IBD) estimations, in a sample of 1805 individuals and 1.006.703 autosomal mutations from a case-control study of colorectal cancer in Mexico. When using the recommended protocol for quality control assessment, 402 should have been removed due to relatedness. Our purpose was to analyze this value in the context of an admixed population. For that aim, we reanalyzed the sample using two software designed for admixed populations, obtaining estimates of 110 and 70 related individuals to remove. The results showed that the first estimation of relatedness was an effect of the higher Native American contribution in part of the data samples, being a confounding factor for IBD estimations. We conclude in the importance of considering PS and genetic ancestry in order to avoid spurious results, not only in GWAS but also in relatedness analysis.

The evolution of altruism through war is highly sensitive to population structure and to civilian and fighter mortality

The importance of warfare in the evolution of human social behavior remains highly debated. One hypothesis is that intense warfare between groups favored altruism within groups, a hypothesis given some support by computational modeling and, in particular, the work of Choi and Bowles [J.-K. Choi, S. Bowles, Science 318, 636-640 (2007)]. The results of computational models are, however, sensitive to chosen parameter values and a deeper assessment of the plausibility of the parochial altruism hypothesis requires exploring this model in more detail. Here, I use a recently developed method to reexamine Choi and Bowles' model under a much broader range of conditions to those used in the original paper. Although the evolution of altruism is robust to perturbations in most of the default parameters, it is highly sensitive to group size and migration and to the lethality of war. The results show that the degree of genetic differentiation between groups (FST ) produced by Choi and Bowles' original model is much greater than empirical estimates of FST between hunter-gatherer groups. When FST in the model is close to empirically observed values, altruism does not evolve. These results cast doubt on the importance of war in the evolution of human sociality.

Use of craniometric data to facilitate migrant identifications at the United States/Mexico border

OBJECTIVES

Thousands of migrants have died along the United States/Mexico border and many remain unidentified. The purpose of this research is to test whether estimations of population affinity, derived from craniometric data, can facilitate identification of migrant remains and provide a geographic region of origin rather than the broad label Hispanic. The appropriateness of current forensic reference data will also be assessed.

METHODS

A case study combined with craniometric data from positively identified and unidentified migrants from the Pima County Office of the Medical Examiner (n = 489) in Arizona and operation identification (n = 201) in Texas were compared to skeletal data representing Guatemalan (n = 87) and Mexican (n = 65) Mayans. Biological distance and discriminant function analyses were used to assess overall population relationships and classificatory models for forensic anthropological application.

RESULTS

The majority of evidence indicates that estimations of population affinity can assist in the facilitation identification of migrant remains, even when a broad classification is used. Biological distances among the groups suggest that positively identified Guatemalan and Mexican migrants are similar to one another but differ from Guatemalan and Mexican Mayans.

CONCLUSIONS

Population affinity estimations can aid migrant identification, and current reference data used in forensic anthropological practice should be replaced with data from positive identifications. Estimates of geographic origin may be more useful than the broad generic term Hispanic for narrowing down the search for a missing person, but more data and research is needed to achieve this goal. Although, the utility of geographic origin estimates relies on transnational data centralization and sharing, which is not always the case.

Genetic diversity, structure, and admixture in Mayans from Guatemala and Mexico based on 15 short tandem repeats

OBJECTIVE

To analyze the genetic origin, relationships, structure, and admixture in Mayan Native American groups from Guatemala and Mexico based on 15 autosomal short tandem repeats (STRs) loci commonly used in human identification (HID).

METHODS

We genotyped 513 unrelated Mayan samples from Guatemala based on 15 STR loci (AmpFlSTR® Identifiler kit). Moreover, we included 4408 genotypes previously reported, as following: Mayas from Guatemala and Mexico (n = 1666) and from Latin American, European, and African (n = 2742) populations. Forensic parameters, genetic distances, admixture, and population structure were assessed.

RESULTS

Forensic parameters of the 15 STRs in different Mayan groups from Guatemala were reported. Low (Fst = 0.78%; p = 0.000) and non-significant differentiation (Fst = 1.8%; p = 0.108) were observed in Mayas from Guatemala and Mexico, respectively. The relative homogeneity observed among Mayan groups supported theories of extensive pre-Columbian gene flow and trade throughout the Mayan Empire. The distribution of the three Native American ancestries among these Mayan groups did not support the presumable Guatemalan origin of Tojolabal and Lacandon people (South, Mexico). The nonsignificant differentiation between Ladinos and Mayas suggests a relative panmixia in Guatemala. Mestizos from southeastern Mexico and Guatemala constitute a core of Native American ancestry in Latin America related to the Mayan Empire in Central America.

CONCLUSIONS

The higher European admixture and homogeneity in Mexican Mayas of the Yucatan Peninsula suggest more intensive post-Columbian gene flow in this region than in Guatemalan Mayas.

Identification of Variants (rs11571707, rs144848, and rs11571769) in the BRCA2 Gene Associated with Hereditary Breast Cancer in Indigenous Populations of the Brazilian Amazon

Estimates show that 5–10% of breast cancer cases are hereditary, caused by genetic variants in autosomal dominant genes; of these, 16% are due to germline mutations in the BRCA1 and BRCA2 genes. The comprehension of the mutation profile of these genes in the Brazilian population, particularly in Amazonian Amerindian groups, is scarce. We investigated fifteen polymorphisms in the BRCA1 and BRCA2 genes in Amazonian Amerindians and compared the results with the findings of global populations publicly available in the 1000 Genomes Project database. Our study shows that three variants (rs11571769, rs144848, and rs11571707) of the BRCA2 gene, commonly associated with hereditary breast cancer, had a significantly higher allele frequency in the Amazonian Amerindian individuals in comparison with the African, American, European, and Asian groups analyzed. These data outline the singular genetic profiles of the indigenous population from the Brazilian Amazon region. The knowledge about BRCA1 and BRCA2 variants is critical to establish public policies for hereditary breast cancer screening in Amerindian groups and populations admixed with them, such as the Brazilian population.

The genetic structure and adaptation of Andean highlanders and Amazonians are influenced by the interplay between geography and culture

Western South America was one of the worldwide cradles of civilization. The well-known Inca Empire was the tip of the iceberg of an evolutionary process that started 11,000 to 14,000 years ago. Genetic data from 18 Peruvian populations reveal the following: 1) The between-population homogenization of the central southern Andes and its differentiation with respect to Amazonian populations of similar latitudes do not extend northward. Instead, longitudinal gene flow between the northern coast of Peru, Andes, and Amazonia accompanied cultural and socioeconomic interactions revealed by archeology. This pattern recapitulates the environmental and cultural differentiation between the fertile north, where altitudes are lower, and the arid south, where the Andes are higher, acting as a genetic barrier between the sharply different environments of the Andes and Amazonia. 2) The genetic homogenization between the populations of the arid Andes is not only due to migrations during the Inca Empire or the subsequent colonial period. It started at least during the earlier expansion of the Wari Empire (600 to 1,000 years before present). 3) This demographic history allowed for cases of positive natural selection in the high and arid Andes vs. the low Amazon tropical forest: in the Andes, a putative enhancer in HAND2-AS1 (heart and neural crest derivatives expressed 2 antisense RNA1, a noncoding gene related to cardiovascular function) and rs269868-C/Ser1067 in DUOX2 (dual oxidase 2, related to thyroid function and innate immunity) genes and, in the Amazon, the gene encoding for the CD45 protein, essential for antigen recognition by T and B lymphocytes in viral-host interaction.

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.

Demographic history and selection at HLA loci in Native Americans

The American continent was the last to be occupied by modern humans, and native populations bear the marks of recent expansions, bottlenecks, natural selection, and population substructure. Here we investigate how this demographic history has shaped genetic variation at the strongly selected HLA loci. In order to disentangle the relative contributions of selection and demography process, we assembled a dataset with genome-wide microsatellites and HLA-A, -B, -C, and -DRB1 typing data for a set of 424 Native American individuals. We find that demographic history explains a sizeable fraction of HLA variation, both within and among populations. A striking feature of HLA variation in the Americas is the existence of alleles which are present in the continent but either absent or very rare elsewhere in the world. We show that this feature is consistent with demographic history (i.e., the combination of changes in population size associated with bottlenecks and subsequent population expansions). However, signatures of selection at HLA loci are still visible, with significant evidence selection at deeper timescales for most loci and populations, as well as population differentiation at HLA loci exceeding that seen at neutral markers.

Contrasting the ancestry patterns of three distinct population groups from the northernmost region of South America

Colombia, located in the north of the South American subcontinent is a country of great interest for population genetic studies given its high ethnic and cultural diversity represented by the admixed population, 102 indigenous peoples and African descent populations. In this study, an analysis of the genetic structure and ancestry was performed based on 46 ancestry informative INDEL markers (AIM-INDELs) and considering the genealogical and demographic variables of 451 unrelated individuals belonging to nine Native American, two African American, and four multiple ancestry populations. Measures of genetic diversity, ancestry components, and genetic substructure were analyzed to build a population model typical of the northernmost part of the South American continent. The model suggests three types of populations: Native American, African American, and multiple ancestry. The results support hypotheses posed by other authors about issues like the peopling of South America and the existence of two types of Native American ancestry. This last finding could be crucial for future research on the peopling of Colombia and South America in that a single origin of all indigenous communities should not be assumed. It then would be necessary to consider other events that could explain their genetic variability and complexity throughout the continent.

Exome Sequencing of Native Populations From the Amazon Reveals Patterns on the Peopling of South America

Studies on the peopling of South America have been limited by the paucity of sequence data from Native Americans, especially from the east part of the Amazon region. Here, we investigate the whole exome variation from 58 Native American individuals (eight different populations) from the Amazon region and draw insights into the peopling of South America. By using the sequence data generated here together with data from the public domain, we confirmed a strong genetic distinction between Andean and Amazonian populations. By testing distinct demographic models, our analysis supports a scenario of South America occupation that involves migrations along the Pacific and Atlantic coasts. Occupation of the southeast part of South America would involve migrations from the north, rather than from the west of the continent.

Comparative assessment of range-wide patterns of genetic diversity and structure with SNPs and microsatellites: A case study with Iberian amphibians

Reduced representation genome sequencing has popularized the application of single nucleotide polymorphisms (SNPs) to address evolutionary and conservation questions in nonmodel organisms. Patterns of genetic structure and diversity based on SNPs often diverge from those obtained with microsatellites to different degrees, but few studies have explicitly compared their performance under similar sampling regimes in a shared analytical framework. We compared range‐wide patterns of genetic structure and diversity in two amphibians endemic to the Iberian Peninsula: Hyla molleri and Pelobates cultripes, based on microsatellite (18 and 14 loci) and SNP (15,412 and 33,140 loci) datasets of comparable sample size and spatial extent. Model‐based clustering analyses with STRUCTURE revealed minor differences in genetic structure between marker types, but inconsistent values of the optimal number of populations (K) inferred. SNPs yielded more repeatable and less admixed ancestries with increasing K compared to microsatellites. Genetic diversity was weakly correlated between marker types, with SNPs providing a better representation of southern refugia and of gradients of genetic diversity congruent with the demographic history of both species. Our results suggest that the larger number of loci in a SNP dataset can provide more reliable inferences of patterns of genetic structure and diversity than a typical microsatellite dataset, at least at the spatial and temporal scales investigated.

Pre-colonial Amerindian legacies in forest composition of southern Brazil

Past human societies have left persistent marks on forests worldwide. However, the degree to which pre-colonial Amerindian societies have affected forest structure is still not fully understood, especially in southern Brazil. This study investigated the influence of two distinct Amerindian groups (Southern-Jê and Guarani) over tree composition of forest fragments in the State of Santa Catarina. Vegetation data was obtained from the Santa Catarina Forest and Floristic Inventory (SCFFI): a statewide systematic vegetation sampling project. Archaeological data was collated from literature reviews as well as existing databases for archaeological sites occupied by Guarani and Southern-Jê groups. Using these sites of known Amerindian occupation, and corresponding environmental variables, ecological niche models were developed for each Amerindian group, predicting potential archaeological sites occupied by these groups across southern Brazil. Maps of these potential occupation sites of pre-colonial Amerindian groups were compared with 417 corresponding floristic inventory plots. Redundancy analysis (RDA) was used to identify floristic composition patterns linked to areas with a high probability of Southern-Jê or Guarani presence. Southern-Jê and Guarani pre-colonial occupations overlapped near main rivers; however, Southern-Jê groups generally occupied elevated areas whereas Guarani occupied mostly coastal areas. We observed differences in forest composition associated with the predicted occurrence of these pre-colonial Amerindian groups. Based on these results, we argue there is a relationship between tree species distribution and pre-colonial human occupation by these two Amerindian groups.

The history behind the mosaic of the Americas

Focusing on literature published in 2018-2020, we review inferences about: (i) how ancient DNA is contributing to clarify the peopling of the Americas and the dispersal of its first inhabitants, (ii) how the interplay between environmental diversity and culture has influenced the genetic structure and adaptation of Andean and Amazon populations, (iii) how genetics has contributed to our understanding of the Pre-Columbian Tupi expansion in Eastern South America, (iv) the subcontinental origins and dynamics of Post-Columbian admixture in the Americas, and finally, (v) episodes of adaptive natural selection in the American continent, particularly in the high altitudes of the Andes.

Major inconsistencies of inferred population genetic structure estimated in a large set of domestic horse breeds using microsatellites

STRUCTURE remains the most applied software aimed at recovering the true, but unknown, population structure from microsatellite or other genetic markers. About 30% of structure-based studies could not be reproduced (Molecular Ecology, 21, 2012, 4925). Here we use a large set of data from 2,323 horses from 93 domestic breeds plus the Przewalski horse, typed at 15 microsatellites, to evaluate how program settings impact the estimation of the optimal number of population clusters K opt that best describe the observed data. Domestic horses are suited as a test case as there is extensive background knowledge on the history of many breeds and extensive phylogenetic analyses. Different methods based on different genetic assumptions and statistical procedures (dapc, flock, PCoA, and structure with different run scenarios) all revealed general, broad-scale breed relationships that largely reflect known breed histories but diverged how they characterized small-scale patterns. structure failed to consistently identify K opt using the most widespread approach, the ΔK method, despite very large numbers of MCMC iterations (3,000,000) and replicates (100). The interpretation of breed structure over increasing numbers of K, without assuming a K opt, was consistent with known breed histories. The over-reliance on K opt should be replaced by a qualitative description of clustering over increasing K, which is scientifically more honest and has the advantage of being much faster and less computer intensive as lower numbers of MCMC iterations and repetitions suffice for stable results. Very large data sets are highly challenging for cluster analyses, especially when populations with complex genetic histories are investigated.

Morphological variation of the early human remains from Quintana Roo, Yucatán Peninsula, Mexico: Contributions to the discussions about the settlement of the Americas

The human settlement of the Americas has been a topic of intense debate for centuries, and there is still no consensus on the tempo and mode of early human dispersion across the continent. When trying to explain the biological diversity of early groups across North, Central and South America, studies have defended a wide range of dispersion models that tend to oversimplify the diversity observed across the continent. In this study, we aim to contribute to this debate by exploring the cranial morphological affinities of four late Pleistocene/early Holocene specimens recovered from the caves of Quintana Roo, Mexico. The four specimens are among the earliest human remains known in the continent and permit the contextualization of biological diversity present during the initial millennia of human presence in the Americas. The specimens were compared to worldwide reference series through geometric morphometric analyses of 3D anatomical landmarks. Morphological data were analyzed through exploratory visual multivariate analyses and multivariate classification based on Mahalanobis distances. The results show very different patterns of morphological association for each Quintana Roo specimen, suggesting that the early populations of the region already shared a high degree of morphological diversity. This contrasts with previous studies of South American remains and opens the possibility that the initial populations of North America already had a high level of morphological diversity, which was reduced as populations dispersed into the southern continent. As such, the study of these rare remains illustrates that we are probably still underestimating the biological diversity of early Americans.

Colonialism and the co-evolution of ethnic and genetic structure in New Mexico

OBJECTIVE

Socially constructed ethnic identities are frequently rooted in beliefs about common descent that form when people with disparate cultures, languages, and biology come into contact. This study explores connections between beliefs about common descent, as represented by ethnic nomenclatures, and histories of migration and isolation ascertained from genomic data in New Mexicans of Spanish-speaking descent (NMS).

MATERIALS AND METHODS

We interviewed 507 NMS who further identified using one of seven ethnic terms that they associated with beliefs about connections to past ancestors. For groups of individuals who identified using each term, we estimated biogeographic ancestry, fit admixture models to ancestry distributions, and partitioned genetic distance into admixture and drift components.

RESULTS

Regardless of which ethnic term they used, all NMS had appreciable Native American (avg. 27%) and European ancestry (avg.71%). However, individuals who identified using terms associated with beliefs connecting them to colonial-period Spanish ancestors had significantly higher European ancestry than individuals who identified using terms associated with ancestral connections to post-colonial-period migrants from Mexico. Model-fitting analyses show that this ancestry difference reflects post-colonial gene flow with non-NMS European Americans, not colonial-period gene flow with Spaniards. Drift, not admixture, accounted for most of the genetic distance between NMS who expressed connections to Mexican versus Spanish ancestors, reflecting relative isolation of New Mexico and Mexico through the 19th century.

DISCUSSION

Patterns of genomic diversity in NMS are consistent with beliefs about common descent in showing that New Mexico was isolated for generations following initial colonization. They are inconsistent with these beliefs in showing that all NMS have substantial European and Native American ancestry, and in showing that a proportion of European ancestry derives from post-colonial-period admixture with non-NMS European Americans. Our findings provide insights into the construction of ethnic identity in contexts of migration and isolation in New Mexico and, potentially, throughout human prehistory.

Molecular analysis of an ancient Thule population at Nuvuk, Point Barrow, Alaska

OBJECTIVES

The North American archaeological record supports a Holocene origin of Arctic Indigenous peoples. Although the Paleo-Inuit were present for millennia, archaeological and genetic studies suggest that modern peoples descend from a second, more recent tradition known as the Neo-Inuit. Origins of the Neo-Inuit and their relations to the earlier and later Indigenous peoples are an area of active study. Here, we genetically analyze the maternal lineages present at Nuvuk, once the northernmost community in Alaska and located in a region identified as a possible origin point of the Neo-Inuit Thule. The cemetery at Nuvuk contains human remains representing a nearly one thousand year uninterrupted occupation from early Thule to post-contact Iñupiat.

MATERIALS AND METHODS

We selected 44 individuals from Nuvuk with calibrated dates between 981 AD and 1885 AD for molecular analysis. We amplified and sequenced the hypervariable segment I of the mitogenome. We compared the Nuvuk data with previously published sequences from 68 modern and ancient communities from across Asia and North America. Phylogeographic analyses suggest possible scenarios of Holocene Arctic and sub-Arctic population movements.

RESULTS

We successfully retrieved sequence data from 39 individuals. Haplogroup frequencies in Nuvuk were typed as 66.7% A2b1, 25.6% A2a, and 7.7% D4b1a2a1a. These results suggest that the population at Nuvuk was closest to the ancient Thule and modern Inuit of Canada, and to the Siberian Naukan people. We confirm that haplogroups A2a, A2b1, D2a, and D4b1a2a1a appear at high frequency in Arctic and sub-Arctic populations of North America and Chukotka. Sister clades D2b and D4b1a2a1b are present in Asian and Eastern European populations.

DISCUSSION

The ancient mitochondrial sequences from Nuvuk confirm the link between the North Slope and the Thule who later spread east, and the maternal discontinuity between the Neo-Inuit and Paleo-Inuit. We suggest haplogroups A2a, A2b, and D4b1a2a1a are linked to the ancestors of the Thule in eastern Beringia, whereas the D2 and D4b1a2a1 clades appear to have Asian Holocene origins. Further Siberian and Alaskan genomes are necessary to clarify these population migrations beyond a simple two-wave scenario of Neo-Inuit and Paleo-Inuit.