Alu insertion polymorphisms in Native Americans and related Asian populations

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.

Alu master copies serve as the drivers of differential SINE transposition in recent primate genomes

Alu elements, averaging ~300bp in length, are a family of primate-specific short intersperse nuclear elements (SINEs) with more than one million copies and contributing to ~11% of primate genomes. Despite mostly being shared among primates, our recent study revealed highly differential recent Alu transposition among the genomes of primates from Hominidae and Cercopithecidae families. To understand the underlying mechanism, we analyzed six primate genomes and revealed species- and lineage-specific Alu profile exclusively defined by AluY composition. Among all Alus from the 6 genomes, we identified 5401 Alu master copies with 99% being from the AluY subfamily. The numbers of Alu master copies are positively correlated to the number of AluY elements in the genomes with the baboon genome having the largest number of most recent Alu master copies at high activities, while the crab-eating macaque genome having a low number of Alu master copies with low activity. Furthermore, the expression level of Alu master copies is positively correlated with their transposition activity. Our results support the concept that Alu transposition in primate genomes is driven by a small number of master copies, the number and relative activity of which contribute to the differential Alu transposition in recent primate genomes.

Comparative Genomics Analysis Reveals High Levels of Differential Retrotransposition among Primates from the Hominidae and the Cercopithecidae Families

Mobile elements (MEs), making ∼50% of primate genomes, are known to be responsible for generating inter- and intra-species genomic variations and play important roles in genome evolution and gene function. Using a bioinformatics comparative genomics approach, we performed analyses of species-specific MEs (SS-MEs) in eight primate genomes from the families of Hominidae and Cercopithecidae, focusing on retrotransposons. We identified a total of 230,855 SS-MEs, with which we performed normalization based on evolutionary distances, and we also analyzed the most recent SS-MEs in these genomes. Comparative analysis of SS-MEs reveals striking differences in ME transposition among these primate genomes. Interesting highlights of our results include: 1) the baboon genome has the highest number of SS-MEs with a strong bias for SINEs, while the crab-eating macaque genome has a sustained extremely low transposition for all ME classes, suggesting the existence of a genome-wide mechanism suppressing ME transposition; 2) while SS-SINEs represent the dominant class in general, the orangutan genome stands out by having SS-LINEs as the dominant class; 3) the human genome stands out among the eight genomes by having the largest number of recent highly active ME subfamilies, suggesting a greater impact of ME transposition on its recent evolution; and 4) at least 33% of the SS-MEs locate to genic regions, including protein coding regions, presenting significant potentials for impacting gene function. Our study, as the first of its kind, demonstrates that mobile elements evolve quite differently among these primates, suggesting differential ME transposition as an important mechanism in primate evolution.

Genetic dissection of five ethnic groups from Punjab, North-West India-A study based on Autosomal Markers

The present study assessed the applicability of Alu insertion elements and Single Nucleotide Polymorphisms (SNPs) in forensic identification and estimated the extent of genetic variation in five major ethnic groups of Punjab, North-West India. A total of 1012 unrelated samples belonging to Banias, Brahmins, Jat Sikhs, Khatris and Scheduled Castes were genotyped for four Alu elements (ACE, APO, PLAT, D1) and six Single Nucleotide Polymorphisms [ESR (PvuII), LPL (PvuII), HTR2A (MspI), DRD2 Taq1A, Taq1B, Taq1D]. Allele frequencies observed heterozygosity and forensic efficacy parameters were determined. The data on the genetic affinity of the studied populations among themselves and with other populations of India was also analysed using a Neighbor-Joining tree and multidimensional scaling plot respectively. All the 10 loci were polymorphic and their average observed heterozygosity ranged from 0.3872 (Banias) to 0.4311 (Scheduled Castes). Allele frequency variation at the 9 out of 10 loci led to statistically significant pairwise differences among the five study population groups. The result from AMOVA, Structure analysis, and Phylogenetic tree suggests that these populations are homogenous. In the multidimensional scaling plot, the present study populations formed a compact cluster clearly separated from other populations, suggesting a unique genetic identity of the Punjab populations as a whole. All these observations suggest that either a recent common origin of these populations or extensive gene flow across the populations that dissolve the original genetic differences. The data generated in this study will be useful for forensic genetics, molecular anthropological and demographic studies.

Population genetic structure of traditional populations in the Peruvian Central Andes and implications for South American population history

Molecular-based characterizations of Andean peoples are traditionally conducted in the service of elucidating continent-level evolutionary processes in South America. Consequently, genetic variation among "western" Andean populations is often represented in relation to variation among "eastern" Amazon and Orinoco River Basin populations. This west-east contrast in patterns of population genetic variation is typically attributed to large-scale phenomena, such as dual founder colonization events or differing long-term microevolutionary histories. However, alternative explanations that consider the nature and causes of population genetic diversity within the Andean region remain underexplored. Here we examine population genetic diversity in the Peruvian Central Andes using data from the mtDNA first hypervariable region and Y-chromosome short tandem repeats among 17 newly sampled populations and 15 published samples. Using this geographically comprehensive data set, we first reassessed the currently accepted pattern of western versus eastern population genetic structure, which our results ultimately reject: mtDNA population diversities were lower, rather than higher, within Andean versus eastern populations, and only highland Y-chromosomes exhibited significantly higher within-population diversities compared with eastern groups. Multiple populations, including several highland samples, exhibited low genetic diversities for both genetic systems. Second, we explored whether the implementation of Inca state and Spanish colonial policies starting at about ad 1400 could have substantially restructured population genetic variation and consequently constitute a primary explanation for the extant pattern of population diversity in the Peruvian Central Andes. Our results suggest that Peruvian Central Andean population structure cannot be parsimoniously explained as the sole outcome of combined Inca and Spanish policies on the region's population demography: highland populations differed from coastal and lowland populations in mtDNA genetic structure only; highland groups also showed strong evidence of female-biased gene flow and/or effective sizes relative to other Peruvian ecozones. Taken together, these findings indicate that population genetic structure in the Peruvian Central Andes is considerably more complex than previously reported and that characterizations of and explanations for genetic variation may be best pursued within more localized regions and defined time periods.

Questioning the "melting pot": analysis of Alu inserts in three population samples from Uruguay

The way that immigrants integrate into recipient societies has been discussed for decades, mainly from the perspective of the social sciences. Uruguay, as other American countries, received diffferent waves of European immigrants, although the details of the process of assimilation, when it did occur, are unclear. In this study we used genetic markers to understand the process experienced by the Basques, one of the major migration waves that populated Uruguay, and their relation to other immigrants, as well as to Native American and African descendants. For this purpose, we analyzed the allele frequencies of 10 ALU loci (A25, ACE, APOA1, B65, D1, F13B, PV92, TPA25, HS2.43, and HS4.65) in three samples from Uruguay (two of Basque descendants, one of non-Basque descendants) from two locations: Montevideo and Trinidad. No departure from Hardy-Weinberg expectations was observed, with the exceptions of the APOA1 and D1 loci in the non-Basque descendants' samples. Our data show that the major genetic contribution in the three samples comes from Europe (78-88%), with minor African (10-15%) and Native American (0-10%) contributions. Genetic distances reveal that Basque descendants from Trinidad cluster with Europeans, whereas both Montevideo samples cluster together and are separate from other populations, showing two diffferent types of integration, related to the general characteristics of each regional population.

Alu polymorphisms in the Waorani tribe from the Ecuadorian Amazon reflect the effects of isolation and genetic drift

OBJECTIVES

The Amazon basin is inhabited by some of the most isolated human groups worldwide. Among them, the Waorani tribe is one of the most interesting Native American populations from the anthropological perspective. This study reports a genetic characterization of the Waorani based on autosomal genetic loci.

METHODS

We analyzed 12 polymorphic Alu insertions in 36 Waorani individuals from different communal longhouses settled in the Yasuní National Park.

RESULTS

The most notable finding was the strikingly reduced genetic diversity detected in the Waorani, corroborated by the existence of four monomorphic loci (ACE, APO, FXIIIB, and HS4.65), and of other four Alu markers that were very close to the fixation for the presence (PV92 and D1) or the absence (A25 and HS4.32) of the insertion. Furthermore, results of the centroid analysis supported the notion of the Waorani being one of the Amerindian groups less impacted by gene flow processes.

CONCLUSIONS

The prolonged isolation of the Waorani community, in conjunction with a historically low effective population size and high inbreeding levels, have resulted in the drastic reduction of their genetic diversity, because of the effects of severe genetic drift. Recurrent population bottlenecks most likely determined by certain deep-rooted sociocultural practices of the Waorani (characterized by violence, internal quarrels, and revenge killings until recent times) are likely responsible for this pattern of diversity. The findings of this study illustrate how sociocultural factors can shape the gene pool of human populations.

Towards a reconciling model about the initial peopling of America

The last two decades have seen numerous debates in the field of the initial settlement of America and noteworthy was the disagreement between physical and molecular anthropologists. Recently, it has been pointed out that this discordance could partly originate from the description methods and classification labels used in craniometry, which did not account fairly for the within-sample and within-group variance. From there, a federative model for the initial peopling of America has been designed which could now explain the biological variability found at both the craniofacial and genetic level. This is a major step in the study of the initial settlement of America, which deserved to be highlighted. The present paper recalls the two conflicting models that prevailed for the last 20 years of anthropological studies in America before browsing the newly accepted hypothesis about the origin of the first Amerindians as seen by its authors. Lastly, the article evokes some areas of investigations, which could furnish significant fallouts about the dynamics of the peopling of Americas in the future.

Microevolutionary processes due to landscape features in the province of Jujuy (Argentina)

OBJECTIVES

We seek to evaluate the influence of a diverse and rugged physical environment on the genetic background of human populations.

METHODS

We analyzed eight polymorphic Alu insertions in 226 individuals from Jujuy province (Argentina), which is composed of several regions with well-defined geographical features and marked contrasts between them associated with differences in altitude (range: 700-3300 m). This regional division was used to assess the spatial variation of the Alu diversity.

RESULTS

Deviations from Hardy-Weinberg Equilibrium expectations resulting from heterozygous deficit were found for FXIIIB and PV92 in the highest subpopulations. Several Alu elements showed genetic heterogeneity between the highest region (La Puna) and the lowest regions (Valle and Selva). Similarly, a decreasing trend of the average heterozygosity according to altitude was found. Both the centroid method and the admixture analysis unveiled a gene flow above the average in lowland populations, indicating a higher proportion of foreign genes introduced by immigrants of European and African ancestry. Furthermore, several Alu frequency clines fitting the orientation of the altitude gradient were detected.

CONCLUSIONS

Our study reveals a spatial patterning of the Alu diversity in Jujuy, most likely determined by disparities in landscape and environmental features between the different subregions. Differences in the physical environment would have drastically reduced the homogenizing effects of the gene flow and would have promoted genetic drift episodes in the highest subpopulations. Microevolutionary processes detected in Jujuy have played an important role in the shaping of the gene pool of the populations from this sub-Andean zone from Argentina.

Blood group O alleles in Native Americans: implications in the peopling of the Americas

All major ABO blood alleles are found in most populations worldwide, whereas the majority of Native Americans are nearly exclusively in the O group. O allele molecular characterization could aid in elucidating the possible causes of group O predominance in Native American populations. In this work, we studied exon 6 and 7 sequence diversity in 180 O blood group individuals from four different Mesoamerican populations. Additionally, a comparative analysis of genetic diversity and population structure including South American populations was performed. Results revealed no significant differences among Mesoamerican and South American groups, but showed significant differences within population groups attributable to previously detected differences in genetic drift and founder effects throughout the American continent. Interestingly, in all American populations, the same set of haplotypes O(1), O(1v), and O(1v(G542A)) was present, suggesting the following: (1) that they constitute the main genetic pool of the founding population of the Americas and (2) that they derive from the same ancestral source, partially supporting the single founding population hypothesis. In addition, the consistent and restricted presence of the G542A mutation in Native Americans compared to worldwide populations allows it to be employed as an Ancestry informative marker (AIM). Present knowledge of the peopling of the Americas allows the prediction of the way in which the G542A mutation could have emerged in Beringia, probably during the differentiation process of Asian lineages that gave rise to the founding population of the continent.

Autosomal and X chromosome Alu insertions in Bolivian Aymaras and Quechuas: two languages and one genetic pool

Thirty-two polymorphic Alu insertions (18 autosomal and 14 from the X chromosome) were studied in 192 individuals from two Amerindian populations of the Bolivian Altiplano (Aymara and Quechua speakers: the two main Andean linguistic groups), to provide relevant information about their genetic relationships and demographic processes. The main objective was to determine from genetic data whether the expansion of the Quechua language into Bolivia could be associated with demographic (Inca migration of Quechua-speakers from Peru into Bolivia) or cultural (language imposition by the Inca Empire) processes. Allele frequencies were used to assess the genetic relationships between these two linguistic groups. Our results indicated that the two Bolivian samples showed a high genetic similarity for both sets of markers and were clearly differentiated from the two Peruvian Quechua samples available in the literature. Additionally, our data were compared with the available literature to determine the genetic and linguistic structure, and East-West differentiation in South America. The close genetic relationship between the two Bolivian samples and their differentiation from the Quechua-speakers from Peru suggests that the Quechua language expansion in Bolivia took place without any important demographic contribution. Moreover, no clear geographical or linguistic structure was found for the Alu variation among South Amerindians.

The peopling of America: craniofacial shape variation on a continental scale and its interpretation from an interdisciplinary view

Twenty-two years ago, Greenberg, Turner and Zegura (Curr. Anthropol. 27,477-495, 1986) suggested a multidisciplinary model for the human settlement of the New World. Since their synthesis, several studies based mainly on partial evidence such as skull morphology and molecular genetics have presented competing, apparently mutually exclusive, settlement hypotheses. These contradictory views are represented by the genetic-based Single Wave or Out of Beringia models and the cranial morphology-based Two Components/Stocks model. Here, we present a geometric morphometric analysis of 576 late Pleistocene/early Holocene and modern skulls suggesting that the classical Paleoamerican and Mongoloid craniofacial patterns should be viewed as extremes of a continuous morphological variation. Our results also suggest that recent contact among Asian and American circumarctic populations took place during the Holocene. These results along with data from other fields are synthesized in a model for the settlement of the New World that considers, in an integrative and parsimonious way, evidence coming from genetics and physical anthropology. This model takes into account a founder population occupying Beringia during the last glaciation characterized by high craniofacial diversity, founder mtDNA and Y-chromosome lineages and some private autosomal alleles. After a Beringian population expansion, which could have occurred concomitant with their entry into America, more recent circumarctic gene flow would have enabled the dispersion of northeast Asian-derived characters and some particular genetic lineages from East Asia to America and vice versa.

A reevaluation of the Native American mtDNA genome diversity and its bearing on the models of early colonization of Beringia

The Americas were the last continents to be populated by humans, and their colonization represents a very interesting chapter in our species' evolution in which important issues are still contentious or largely unknown. One difficult topic concerns the details of the early peopling of Beringia, such as for how long it was colonized before people moved into the Americas and the demography of this occupation. A recent work using mitochondrial genome (mtDNA) data presented evidence for a so called "three-stage model" consisting of a very early expansion into Beringia followed by approximately 20,000 years of population stability before the final entry into the Americas. However, these results are in disagreement with other recent studies using similar data and methods. Here, we reanalyze their data to check the robustness of this model and test the ability of Native American mtDNA to discriminate details of the early colonization of Beringia. We apply the Bayesian Skyline Plot approach to recover the past demographic dynamic underpinning these events using different mtDNA data sets. Our results refute the specific details of the "three-stage model", since the early stage of expansion into Beringia followed by a long period of stasis could not be reproduced in any mtDNA data set cleaned from non-Native American haplotypes. Nevertheless, they are consistent with a moderate population bottleneck in Beringia associated with the Last Glacial Maximum followed by a strong population growth around 18,000 years ago as suggested by other recent studies. We suggest that this bottleneck erased the signals of ancient demographic history from recent Native American mtDNA pool, and conclude that the proposed early expansion and occupation of Beringia is an artifact caused by the misincorporation of non-Native American haplotypes.

Mitochondrial population genomics supports a single pre-Clovis origin with a coastal route for the peopling of the Americas

It is well accepted that the Americas were the last continents reached by modern humans, most likely through Beringia. However, the precise time and mode of the colonization of the New World remain hotly disputed issues. Native American populations exhibit almost exclusively five mitochondrial DNA (mtDNA) haplogroups (A-D and X). Haplogroups A-D are also frequent in Asia, suggesting a northeastern Asian origin of these lineages. However, the differential pattern of distribution and frequency of haplogroup X led some to suggest that it may represent an independent migration to the Americas. Here we show, by using 86 complete mitochondrial genomes, that all Native American haplogroups, including haplogroup X, were part of a single founding population, thereby refuting multiple-migration models. A detailed demographic history of the mtDNA sequences estimated with a Bayesian coalescent method indicates a complex model for the peopling of the Americas, in which the initial differentiation from Asian populations ended with a moderate bottleneck in Beringia during the last glacial maximum (LGM), around approximately 23,000 to approximately 19,000 years ago. Toward the end of the LGM, a strong population expansion started approximately 18,000 and finished approximately 15,000 years ago. These results support a pre-Clovis occupation of the New World, suggesting a rapid settlement of the continent along a Pacific coastal route.

Statistical evaluation of alternative models of human evolution

An appropriate model of recent human evolution is not only important to understand our own history, but it is necessary to disentangle the effects of demography and selection on genome diversity. Although most genetic data support the view that our species originated recently in Africa, it is still unclear if it completely replaced former members of the Homo genus, or if some interbreeding occurred during its range expansion. Several scenarios of modern human evolution have been proposed on the basis of molecular and paleontological data, but their likelihood has never been statistically assessed. Using DNA data from 50 nuclear loci sequenced in African, Asian and Native American samples, we show here by extensive simulations that a simple African replacement model with exponential growth has a higher probability (78%) as compared with alternative multiregional evolution or assimilation scenarios. A Bayesian analysis of the data under this best supported model points to an origin of our species approximately 141 thousand years ago (Kya), an exit out-of-Africa approximately 51 Kya, and a recent colonization of the Americas approximately 10.5 Kya. We also find that the African replacement model explains not only the shallow ancestry of mtDNA or Y-chromosomes but also the occurrence of deep lineages at some autosomal loci, which has been formerly interpreted as a sign of interbreeding with Homo erectus.

Genetic, geographic, and environmental correlates of human temporal bone variation

Temporal bone shape has been shown to reflect molecular phylogenetic relationships among hominoids and offers significant morphological detail for distinguishing taxa. Although it is generally accepted that temporal bone shape, like other aspects of morphology, has an underlying genetic component, the relative influence of genetic and environmental factors is unclear. To determine the impact of genetic differentiation and environmental variation on temporal bone morphology, we used three-dimensional geometric morphometric techniques to evaluate temporal bone variation in 11 modern human populations. Population differences were investigated by discriminant function analysis, and the strength of the relationships between morphology, neutral molecular distance, geographic distribution, and environmental variables were assessed by matrix correlation comparisons. Significant differences were found in temporal bone shape among all populations, and classification rates using cross-validation were relatively high. Comparisons of morphological distances to molecular distances based on short tandem repeats (STRs) revealed a significant correlation between temporal bone shape and neutral molecular distance among Old World populations, but not when Native Americans were included. Further analyses suggested a similar pattern for morphological variation and geographic distribution. No significant correlations were found between temporal bone shape and environmental variables: temperature, annual rainfall, latitude, or altitude. Significant correlations were found between temporal bone size and both temperature and latitude, presumably reflecting Bergmann's rule. Thus, temporal bone morphology appears to partially follow an isolation by distance model of evolution among human populations, although levels of correlation show that a substantial component of variation is unexplained by factors considered here.

Identification of a unique Alu-based polymorphism and its use in human population studies

Alu elements represent a family of short interspersed DNA elements (SINEs) found in primate genomes. These are members of a group of transposable elements that integrate into the genome by the process of retrotransposition. Recent integrations of Alu elements within the human genome have generated presence/absence variants useful as DNA markers in human population studies as well as in forensic and paternity analyses. Besides the ease of use, this type of marker is unique because the absence of the Alu represents the ancestral form. We have identified an Alu-based polymorphism that consists of four alleles in which we can predict the evolutionary order. Additionally, we have developed a simple PCR plus restriction endonuclease assay to readily distinguish the four alleles. We have thus far analyzed DNA from a small set of samples comprising ten different ethnic groups. The three populations of African descent exhibited a relatively low frequency of the absence allele in contrast to the other populations, as well as being the only populations in which all four alleles were identified. One presence allele was not found in both European Caucasian and South American populations that were sampled, whereas a different presence allele was not observed among the sampled Asian populations. Additionally, the four-allele system identified variations among populations not observed by simply scoring as presence/absence variants. Therefore, extending beyond the two-allele dimorphic Alu system further elucidates population variations. These features afford this marker as a unique tool in the study of both global and regional analyses of human populations.