Genetic variation and population structure in native Americans

We examined genetic diversity and population structure in the American landmass using 678 autosomal microsatellite markers genotyped in 422 individuals representing 24 Native American populations sampled from North, Central, and South America. These data were analyzed jointly with similar data available in 54 other indigenous populations worldwide, including an additional five Native American groups. The Native American populations have lower genetic diversity and greater differentiation than populations from other continental regions. We observe gradients both of decreasing genetic diversity as a function of geographic distance from the Bering Strait and of decreasing genetic similarity to Siberians--signals of the southward dispersal of human populations from the northwestern tip of the Americas. We also observe evidence of: (1) a higher level of diversity and lower level of population structure in western South America compared to eastern South America, (2) a relative lack of differentiation between Mesoamerican and Andean populations, (3) a scenario in which coastal routes were easier for migrating peoples to traverse in comparison with inland routes, and (4) a partial agreement on a local scale between genetic similarity and the linguistic classification of populations. These findings offer new insights into the process of population dispersal and differentiation during the peopling of the Americas.     

Continent-Wide Decoupling of Y-Chromosomal Genetic Variation from Language and Geography in Native South Americans

Numerous studies of human populations in Europe and Asia have revealed a concordance between their extant genetic structure and the prevailing regional pattern of geography and language. For native South Americans, however, such evidence has been lacking so far. Therefore, we examined the relationship between Y-chromosomal genotype on the one hand, and male geographic origin and linguistic affiliation on the other, in the largest study of South American natives to date in terms of sampled individuals and populations. A total of 1,011 individuals, representing 50 tribal populations from 81 settlements, were genotyped for up to 17 short tandem repeat (STR) markers and 16 single nucleotide polymorphisms (Y-SNPs), the latter resolving phylogenetic lineages Q and C. Virtually no structure became apparent for the extant Y-chromosomal genetic variation of South American males that could sensibly be related to their inter-tribal geographic and linguistic relationships. This continent-wide decoupling is consistent with a rapid peopling of the continent followed by long periods of isolation in small groups. Furthermore, for the first time, we identified a distinct geographical cluster of Y-SNP lineages C-M217 (C3*) in South America. Such haplotypes are virtually absent from North and Central America, but occur at high frequency in Asia. Together with the locally confined Y-STR autocorrelation observed in our study as a whole, the available data therefore suggest a late introduction of C3* into South America no more than 6,000 years ago, perhaps via coastal or trans-Pacific routes. Extensive simulations revealed that the observed lack of haplogroup C3* among extant North and Central American natives is only compatible with low levels of migration between the ancestor populations of C3* carriers and non-carriers. In summary, our data highlight the fact that a pronounced correlation between genetic and geographic/cultural structure can only be expected under very specific conditions, most of which are likely not to have been met by the ancestors of native South Americans.     

The Uto‐Aztecan premolar among North and South Amerindians: Geographic variation and genetics

The Uto‐Aztecan premolar (UAP) is a dental polymorphism characterized by an exaggerated distobuccal rotation of the paracone in combination with the presence of a fossa at the intersection of the distal occlusal ridge and distal marginal ridge of upper first premolars. This trait is important because, unlike other dental variants, it has been found exclusively in Native American populations. However, the trait's temporal and geographic variation has never been fully documented. The discovery of a Uto‐Aztecan premolar in a prehistoric skeletal series from northern South America calls into question the presumed linguistic and geographic limits of this trait. We examined published and unpublished data for this rare but highly distinctive trait in samples representing over 5,000 Native Americans from North and South America. Our findings in living Southwest Amerindian populations corroborate the notion that the variable goes beyond the bounds of the Uto‐Aztecan language family. It is found in prehistoric Native Americans from South America, eastern North America, Northern and Central Mexico, and in living and prehistoric populations in the American Southwest that are not members of the Uto‐Aztecan language stock. The chronology of samples, its geographic distribution, and trait frequencies suggests a North American origin (Southwest) for UAP perhaps between 15,000 BP and 4,000 BP and a rapid and widespread dispersal into South America during the late Holocene. Family data indicate that it may represent an autosomal recessive mutation that occurred after the peopling of the Americas as its geographic range appears to be limited to North and South Amerindian populations. Am J Phys Anthropol, 2010. © 2010 Wiley‐Liss, Inc.     

Genetic variability and level of differentiation among Brazilian pantanal wood stork populations

The Wood Stork (Mycteria americana) is a colonial wading bird of the tropical and lower subtropical zones. We assessed genetic structure within and among five stork colonies from the Brazilian Pantanal and compared our data with those from North American populations. Samples of 234 individuals were studied using protein electrophoresis to evaluate genetic variability and interpopulation differentiation. Of 22 loci examined, 7 were polymorphic (mean heterozygosity = 0.068). The low Fst value (0.005) indicated little intraspecific variation among breeding colonies. Estimated number of migrants per generation based on private alleles (Nm = 11.3) and on Fst (48.8) suggests high gene flow. Nei's genetic distance values among Pantanal colonies ranged from 0.0001 to 0.0034, demonstrating low genetic divergence among populations. Our data can be explained by supposing high gene flow levels among Pantanal colonies, and between North and South American populations, intermediated by a probable interbreeding population in Central America.     

Hierarchical modeling of genome‐wide Short Tandem Repeat (STR) markers infers native American prehistory

This study examines a genome‐wide dataset of 678 Short Tandem Repeat loci characterized in 444 individuals representing 29 Native American populations as well as the Tundra Netsi and Yakut populations from Siberia. Using these data, the study tests four current hypotheses regarding the hierarchical distribution of neutral genetic variation in native South American populations: (1) the western region of South America harbors more variation than the eastern region of South America, (2) Central American and western South American populations cluster exclusively, (3) populations speaking the Chibchan‐Paezan and Equatorial‐Tucanoan language stock emerge as a group within an otherwise South American clade, (4) Chibchan‐Paezan populations in Central America emerge together at the tips of the Chibchan‐Paezan cluster. This study finds that hierarchical models with the best fit place Central American populations, and populations speaking the Chibchan‐Paezan language stock, at a basal position or separated from the South American group, which is more consistent with a serial founder effect into South America than that previously described. Western (Andean) South America is found to harbor similar levels of variation as eastern (Equatorial‐Tucanoan and Ge‐Pano‐Carib) South America, which is inconsistent with an initial west coast migration into South America. Moreover, in all relevant models, the estimates of genetic diversity within geographic regions suggest a major bottleneck or founder effect occurring within the North American subcontinent, before the peopling of Central and South America. Am J Phys Anthropol 2010. © 2009 Wiley‐Liss, Inc.     

Demographic History of Indigenous Populations in Mesoamerica Based on mtDNA Sequence Data

The genetic characterization of Native American groups provides insights into their history and demographic events. We sequenced the mitochondrial D-loop region (control region) of 520 samples from eight Mexican indigenous groups. In addition to an analysis of the genetic diversity, structure and genetic relationship between 28 Native American populations, we applied Bayesian skyline methodology for a deeper insight into the history of Mesoamerica. AMOVA tests applying cultural, linguistic and geographic criteria were performed. MDS plots showed a central cluster of Oaxaca and Maya populations, whereas those from the North and West were located on the periphery. Demographic reconstruction indicates higher values of the effective number of breeding females (Nef) in Central Mesoamerica during the Preclassic period, whereas this pattern moves toward the Classic period for groups in the North and West. Conversely, Nef minimum values are distributed either in the Lithic period (i.e. founder effects) or in recent periods (i.e. population declines). The Mesomerican regions showed differences in population fluctuation as indicated by the maximum Inter-Generational Rate (IGRmax): i) Center-South from the lithic period until the Preclassic; ii) West from the beginning of the Preclassic period until early Classic; iii) North characterized by a wide range of temporal variation from the Lithic to the Preclassic. Our findings are consistent with the genetic variations observed between central, South and Southeast Mesoamerica and the North-West region that are related to differences in genetic drift, structure, and temporal survival strategies (agriculture versus hunter-gathering, respectively). Interestingly, although the European contact had a major negative demographic impact, we detect a previous decline in Mesoamerica that had begun a few hundred years before.     

Biological relationship between Central and South American Chibchan speaking populations: evidence from mtDNA

We examined mitochondrial DNA (mtDNA) haplogroup and haplotype diversity in 188 individuals from three Chibchan (Kogi, Arsario, and Ijka) populations and one Arawak (Wayuú) group from northeast Colombia to determine the biological relationship between lower Central American and northern South American Chibchan speakers. mtDNA haplogroups were obtained for all individuals and mtDNA HVS-I sequence data were obtained for 110 samples. Resulting sequence data were compared to 16 other Caribbean, South, and Central American populations using diversity measures, neutrality test statistics, sudden and spatial mismatch models, intermatch distributions, phylogenetic networks, and a multidimensional scaling plot. Our results demonstrate the existence of a shared maternal genetic structure between Central American Chibchan, Mayan populations and northern South American Chibchan-speakers. Additionally, these results suggest an expansion of Chibchan-speakers into South America associated with a shift in subsistence strategies because of changing ecological conditions that occurred in the region between 10,000-14,000 years before present.     

Molecular variation at functional genes and the history of human populations--data on candidate genes for cardiovascular risk in the Mediterranean

A screening of 22 DNA polymorphisms has been performed in western Mediterranean populations (Iberian Peninsula, Morocco, and Central Mediterranean Islands). The analyzed markers correspond to polymorphic sites in several candidate genes for cardiovascular disease including apolipopoteins and their receptors (APOA1, APOB, APOE, APOC1, APOC2, LPA, and LDLR), genes implied in the hemostasis regulation (Factor VII, alpha and beta-fibrinogen, alpha and beta platelet-integrin, tissue plasminogen activator, and plasminogen activator inhibitor-1), and the angiotensin converting enzyme gene. The results are presented of a partial analysis carried out in following population samples: 6 from the Iberian Peninsula, 2 from Morocco, and 3 from Central Islands. The degree of inter-population diversity was significant and consistent with data from other kind of genetic polymorphisms. The apportionment of the allele frequency variance supported a geographic structure into three main regions: Central Mediterranean Islands, the Iberia Peninsula and North Africa. The genetic distance pattern is compatible with a south-to-north North African influence in the Iberian Peninsula and a remarkable gene flow from sub-Saharan Africa into Morocco. Epidemiologically, North Africa is characterized by high frequencies of LPA PNR alleles with high number of repeats (protective for cardiovascular risk) and high frequencies of the APOE*E4 allele (risk factor) as compared with European populations.     

Population structure of Hispanics in the United States: the multi-ethnic study of atherosclerosis

Using ∼60,000 SNPs selected for minimal linkage disequilibrium, we perform population structure analysis of 1,374 unrelated Hispanic individuals from the Multi-Ethnic Study of Atherosclerosis (MESA), with self-identification corresponding to Central America (n = 93), Cuba (n = 50), the Dominican Republic (n = 203), Mexico (n = 708), Puerto Rico (n = 192), and South America (n = 111). By projection of principal components (PCs) of ancestry to samples from the HapMap phase III and the Human Genome Diversity Panel (HGDP), we show the first two PCs quantify the Caucasian, African, and Native American origins, while the third and fourth PCs bring out an axis that aligns with known South-to-North geographic location of HGDP Native American samples and further separates MESA Mexican versus Central/South American samples along the same axis. Using k-means clustering computed from the first four PCs, we define four subgroups of the MESA Hispanic cohort that show close agreement with self-identification, labeling the clusters as primarily Dominican/Cuban, Mexican, Central/South American, and Puerto Rican. To demonstrate our recommendations for genetic analysis in the MESA Hispanic cohort, we present pooled and stratified association analysis of triglycerides for selected SNPs in the LPL and TRIB1 gene regions, previously reported in GWAS of triglycerides in Caucasians but as yet unconfirmed in Hispanic populations. We report statistically significant evidence for genetic association in both genes, and we further demonstrate the importance of considering population substructure and genetic heterogeneity in genetic association studies performed in the United States Hispanic population.     

Disruption of historical patterns of isolation by distance in coastal steelhead

When natural populations exchange migrants at a rate proportional to their geographic distance, population genetics theory leads to the expectation of a pattern of isolation-by-distance (IBD), whereby geographic and genetic distance are correlated. However, the presence or absence of such patterns in modern populations may not fully reflect the historical relationships among those populations. Thus, historical samples, collected prior to modern human impacts, can often provide a critical baseline for comparison with modern populations. Steelhead, the anadromous form of rainbow trout, Oncorhynchus mykiss, are native to western North America and are endangered or threatened throughout most of California, near the southern extent of their native range. Population samples of steelhead collected in 1897 and 1909 in Central California rivers provided the opportunity to evaluate the historical genetic composition and population structure of these threatened fish. Here we show that these steelhead populations had a historically strong correlation between genetic and geographic distance that has been virtually erased in modern populations, suggesting that current relationships among modern steelhead populations are no longer reflective of natural migratory pathways. This demonstrates the critical role of migration in maintaining population relationships of threatened species and highlights the importance of natural history museums in providing historical baseline information.     

Molecular population genetics of the malaria vector Anopheles darlingi in Central and South America

To analyze the genetic relatedness and phylogeographic structure of Anopheles darlingi from 19 localities throughout Central and South America, we used a minimum spanning network, diversity measures, differentiation, neutrality tests, and mismatch distribution with mitochondrial cytochrome oxidase subunit I (COI) sequences. All the Central American haplotypes were separated by seven mutational steps from the South American haplotypes and the FST distance-based neighbor-joining tree showed a primary division between Central and South America, evidence for a putative gene pool division. More ancestral and diverse haplotypes were found in Amazonian and southern Brazil populations, suggesting that Central American populations may have originated in South America. The patterns of the mtDNA haplotype diversity and five of six tests for equilibrium implicate demographic expansion in the South American populations as the historical structure, but mismatch distribution depicts populations at mutation drift equilibrium (MDE). In South America, the departure from equilibrium was consistent with an expansion that occurred during the Pleistocene.     

Living in two worlds: Evolutionary mechanisms act differently in the native and introduced ranges of an invasive plant

Identifying the factors that influence spatial genetic structure among populations can provide insights into the evolution of invasive plants. In this study, we used the common reed (Phragmites australis), a grass native in Europe and invading North America, to examine the relative importance of geographic, environmental (represented by climate here), and human effects on population genetic structure and its changes during invasion. We collected samples of P. australis from both the invaded North American and native European ranges and used molecular markers to investigate the population genetic structure within and between ranges. We used path analysis to identify the contributions of each of the three factors—geographic, environmental, and human‐related—to the formation of spatial genetic patterns. Genetic differentiation was observed between the introduced and native populations, and their genetic structure in the native and introduced ranges was different. There were strong effects of geography and environment on the genetic structure of populations in the native range, but the human‐related factors manifested through colonization of anthropogenic habitats in the introduced range counteracted the effects of environment. The between‐range genetic differences among populations were mainly explained by the heterogeneous environment between the ranges, with the coefficient 2.6 times higher for the environment than that explained by the geographic distance. Human activities were the primary contributor to the genetic structure of the introduced populations. The significant environmental divergence between ranges and the strong contribution of human activities to the genetic structure in the introduced range suggest that invasive populations of P. australis have evolved to adapt to a different climate and to human‐made habitats in North America.     

Linguistic and maternal genetic diversity are not correlated in Native Mexicans

Mesoamerica, defined as the broad linguistic and cultural area from middle southern Mexico to Costa Rica, might have played a pivotal role during the colonization of the American continent. The Mesoamerican isthmus has constituted an important geographic barrier that has severely restricted gene flow between North and South America in pre-historical times. Although the Native American component has been already described in admixed Mexican populations, few studies have been carried out in native Mexican populations. In this study, we present mitochondrial DNA (mtDNA) sequence data for the first hypervariable region (HVR-I) in 477 unrelated individuals belonging to 11 different native populations from Mexico. Almost all of the Native Mexican mtDNAs could be classified into the four pan-Amerindian haplogroups (A2, B2, C1, and D1); only two of them could be allocated to the rare Native American lineage D4h3. Their haplogroup phylogenies are clearly star-like, as expected from relatively young populations that have experienced diverse episodes of genetic drift (e.g., extensive isolation, genetic drift, and founder effects) and posterior population expansions. In agreement with this observation, Native Mexican populations show a high degree of heterogeneity in their patterns of haplogroup frequencies. Haplogroup X2a was absent in our samples, supporting previous observations where this clade was only detected in the American northernmost areas. The search for identical sequences in the American continent shows that, although Native Mexican populations seem to show a closer relationship to North American populations, they cannot be related to a single geographical region within the continent. Finally, we did not find significant population structure in the maternal lineages when considering the four main and distinct linguistic groups represented in our Mexican samples (Oto-Manguean, Uto-Aztecan, Tarascan, and Mayan), suggesting that genetic divergence predates linguistic diversification in Mexico.     

Unexpected low genetic differentiation among Allopatric species of section Algarobia of Prosopis (Leguminosae)

The Chaqueña Biogeographic Province, in South America, is the main diversity centre of Prosopis. A group of sympatric species of Section Algarobia in this region constitutes a syngameon, characterised by frequent hybridization and introgression. These processes have been postulated as responsible for the low genetic differentiation observed among species within this group. In this study genetic variability and differentiation among geographically isolated species of the same section was analysed through isoenzyme electrophoresis. Variability parameters and fixation indices were estimated to determine the genetic structure of populations. Two Argentinean, P. ruscifolia (‘vinal’) and P. flexuosa (‘algarrobo amarillo’), and one North American species, P. glandulosa (honey mesquite) were studied. All of them showed, similarly to other species of the same section, high genetic variability and exhibit homozygote excess, probably due to population substructure and low rates of selfing. In contrast to our hypothesis, genetic similarity among species is not related to geographic distance. Genetic distances between P. glandulosa and South American species are similar to those observed among species of this subcontinent. The results obtained suggest that the high genetic similarity among the species of the section Algarobia studied is not due to hybridization.     

Relationships among cave swallow populations (Petrochelidon fulva) determined by comparisons of microsatellite and cytochrome b data

The cave swallow (Petrochelidon fulva) is a polytypic species with isolated populations in northwestern South America, southwestern North America, Yucatan, Greater Antilles, and Florida. We compared microsatellite genotypes of at least five individuals each from five populations and cytochrome b sequence data of two individuals each from seven populations plus two outgroups. Microsatellite allelic diversity was substantial, and the North American populations were about equally variable. In contrast, the Ecuadorian population had far less genetic variation. Gene flow was apparent among populations, especially between Texas and Florida. Genetic structure was greater than in widespread North American species but less than that of sedentary Neotropical birds. Microsatellite genetic distances indicated a close affinity between Ecuadorian and northern populations, especially Texas and Florida, but this finding was inconsistent with cytochrome b data, which indicated that the Ecuadorian population is the clear outgroup to the northern populations. Its outgroup status is consistent with recent classifications that designate South American populations as their own species (P. rufocollaris). The cytochrome b data further suggested that the northern populations are divided into two clades: Texas/Yucatan and Florida/Greater Antilles. The microsatellite data incorrectly measured the diversity and affinities of Ecuadorian birds apparently because of an ascertainment bias that results from the use of heterologous PCR primers. Despite these problems in measuring phylogenetic relationships, the microsatellite data appeared to work well as a population genetic marker in indicating population structure and gene flow.     

A striking example of the founder effect in the mollusc Littorina saxatilis

Two South African populations of Littorina saxatilis were examined by starch-gel electrophoresis at 16 enzyme loci and compared with 13 populations of North Atlantic saxatilis from both American and European coasts, and with six British populations of the closely related species Littorina arcana. The South African animals showed a severely reduced heterozygosity (– 0.052) compared with Atlantic populations of saxatilis ( = 0.181), and the mean genetic distance between the two areas was high ( = 0.203) compared with distances within the North Atlantic saxatilis populations (D = 0.034). In fact, the saxatilis from South Africa were genetically more distant from the North Atlantic samples of L. saxatilis than were the arcana from British shores. The reduced genetic heterozygosity and genetic divergence of the South African populations is attributed to founder effects following a postulated recent introduction by man.     

Phylogenetic relationships of Loxosceles and Sicarius spiders are consistent with Western Gondwanan vicariance

The modern geographic distribution of the spider family Sicariidae is consistent with an evolutionary origin on Western Gondwana. Both sicariid genera, Loxosceles and Sicarius are diverse in Africa and South/Central America. Loxosceles are also diverse in North America and the West Indies, and have species described from Mediterranean Europe and China. We tested vicariance hypotheses using molecular phylogenetics and molecular dating analyses of 28S, COI, 16S, and NADHI sequences. We recover reciprocal monophyly of African and South American Sicarius, paraphyletic Southern African Loxosceles and monophyletic New World Loxosceles within which an Old World species group that includes L. rufescens is derived. These patterns are consistent with a sicariid common ancestor on Western Gondwana. North American Loxosceles are monophyletic, sister to Caribbean taxa, and resolved in a larger clade with South American Loxosceles. With fossil data this pattern is consistent with colonization of North America via a land bridge predating the modern Isthmus of Panama.