Mitochondrial DNA studies of Native Americans: Conceptions and misconceptions of the population prehistory of the Americas

Molecular characterization of mitochondrial Amerindian haplogroups and the amelogenin gene in human ancient DNA from three archaeological sites in Lambayeque - Peru

Important pre-Inca civilizations, known by their great political and religious structures, inhabited the northern coast of Peru. Archeological and anthropological studies have shown that people from these villages have hierarchical strata, but the genetic structure has been poorly studied. Here, we aimed to perform a molecular characterization of the Amerindian maternal lineages and the amelogenin gene in skeletons collected from three archeological sites in Lambayeque. Ancient DNA (aDNA) samples were analyzed with conventional PCR to assess the nine-base pair (9 bp) deletion corresponding to mitochondrial haplogroup B and the identification of haplogroups A, C, and D were obtained with PCR-RFLP experiments. The sex was characterized via amplification of the AMEL(X/Y) locus. Haplogroup frequencies were compared with available data from other ancient and modern civilizations from the Peruvian coast and highlands using statistical methods. Our results showed that haplogroup C had the highest frequency, while haplogroup B showed variable diversity in the analyzed populations. The meta-analysis revealed a positive correlation among some coastal villages. We concluded that ancient populations analyzed in our study showed the presence of four Amerindian mitochondrial haplogroups, which is consistent with previous studies.

Annotated Draft Genome Assemblies for the Northern Bobwhite (Colinus virginianus) and the Scaled Quail (Callipepla squamata) Reveal Disparate Estimates of Modern Genome Diversity and Historic Effective Population Size

Northern bobwhite (Colinus virginianus; hereafter bobwhite) and scaled quail (Callipepla squamata) populations have suffered precipitous declines across most of their US ranges. Illumina-based first- (v1.0) and second- (v2.0) generation draft genome assemblies for the scaled quail and the bobwhite produced N50 scaffold sizes of 1.035 and 2.042 Mb, thereby producing a 45-fold improvement in contiguity over the existing bobwhite assembly, and ≥90% of the assembled genomes were captured within 1313 and 8990 scaffolds, respectively. The scaled quail assembly (v1.0 = 1.045 Gb) was ∼20% smaller than the bobwhite (v2.0 = 1.254 Gb), which was supported by kmer-based estimates of genome size. Nevertheless, estimates of GC content (41.72%; 42.66%), genome-wide repetitive content (10.40%; 10.43%), and MAKER-predicted protein coding genes (17,131; 17,165) were similar for the scaled quail (v1.0) and bobwhite (v2.0) assemblies, respectively. BUSCO analyses utilizing 3023 single-copy orthologs revealed a high level of assembly completeness for the scaled quail (v1.0; 84.8%) and the bobwhite (v2.0; 82.5%), as verified by comparison with well-established avian genomes. We also detected 273 putative segmental duplications in the scaled quail genome (v1.0), and 711 in the bobwhite genome (v2.0), including some that were shared among both species. Autosomal variant prediction revealed ∼2.48 and 4.17 heterozygous variants per kilobase within the scaled quail (v1.0) and bobwhite (v2.0) genomes, respectively, and estimates of historic effective population size were uniformly higher for the bobwhite across all time points in a coalescent model. However, large-scale declines were predicted for both species beginning ∼15-20 KYA.

A draft de novo genome assembly for the northern bobwhite (Colinus virginianus) reveals evidence for a rapid decline in effective population size beginning in the Late Pleistocene

Wild populations of northern bobwhites (Colinus virginianus; hereafter bobwhite) have declined across nearly all of their U.S. range, and despite their importance as an experimental wildlife model for ecotoxicology studies, no bobwhite draft genome assembly currently exists. Herein, we present a bobwhite draft de novo genome assembly with annotation, comparative analyses including genome-wide analyses of divergence with the chicken (Gallus gallus) and zebra finch (Taeniopygia guttata) genomes, and coalescent modeling to reconstruct the demographic history of the bobwhite for comparison to other birds currently in decline (i.e., scarlet macaw; Ara macao). More than 90% of the assembled bobwhite genome was captured within <40,000 final scaffolds (N50 = 45.4 Kb) despite evidence for approximately 3.22 heterozygous polymorphisms per Kb, and three annotation analyses produced evidence for >14,000 unique genes and proteins. Bobwhite analyses of divergence with the chicken and zebra finch genomes revealed many extremely conserved gene sequences, and evidence for lineage-specific divergence of noncoding regions. Coalescent models for reconstructing the demographic history of the bobwhite and the scarlet macaw provided evidence for population bottlenecks which were temporally coincident with human colonization of the New World, the late Pleistocene collapse of the megafauna, and the last glacial maximum. Demographic trends predicted for the bobwhite and the scarlet macaw also were concordant with how opposing natural selection strategies (i.e., skewness in the r-/K-selection continuum) would be expected to shape genome diversity and the effective population sizes in these species, which is directly relevant to future conservation efforts.

Genetic differences between Chibcha and Non-Chibcha speaking tribes based on mitochondrial DNA (mtDNA) haplogroups from 21 Amerindian tribes from Colombia

We analyzed the frequency of four mitochondrial DNA haplogroups in 424 individuals from 21 Colombian Amerindian tribes. Our results showed a high degree of mtDNA diversity and genetic heterogeneity. Frequencies of mtDNA haplogroups A and C were high in the majority of populations studied. The distribution of these four mtDNA haplogroups from Amerindian populations was different in the northern region of the country compared to those in the south. Haplogroup A was more frequently found among Amerindian tribes in northern Colombia, while haplogroup D was more frequent among tribes in the south. Haplogroups A, C and D have clinal tendencies in Colombia and South America in general. Populations belonging to the Chibcha linguistic family of Colombia and other countries nearby showed a strong genetic differentiation from the other populations tested, thus corroborating previous findings. Genetically, the Ingano, Paez and Guambiano populations are more closely related to other groups of south eastern Colombia, as also inferred from other genetic markers and from archeological data. Strong evidence for a correspondence between geographical and linguistic classification was found, and this is consistent with evidence that gene flow and the exchange of customs and knowledge and language elements between groups is facilitated by close proximity.

Calculating expected DNA remnants from ancient founding events in human population genetics

Background

Recent advancements in sequencing and computational technologies have led to rapid generation and analysis of high quality genetic data. Such genetic data have achieved wide acceptance in studies of historic human population origins and admixture. However, in studies relating to small, recent admixture events, genetic factors such as historic population sizes, genetic drift, and mutation can have pronounced effects on data reliability and utility. To address these issues we conducted genetic simulations targeting influential genetic parameters in admixed populations.

Results

We performed a series of simulations, adjusting variable values to assess the affect of these genetic parameters on current human population studies and what these studies infer about past population structure. Final mean allele frequencies varied from 0.0005 to over 0.50, depending on the parameters.

Conclusion

The results of the simulations illustrate that, while genetic data may be sensitive and powerful in large genetic studies, caution must be used when applying genetic information to small, recent admixture events. For some parameter sets, genetic data will not be adequate to detect historic admixture. In such cases, studies should consider anthropologic, archeological, and linguistic data where possible.

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.

Beringian standstill and spread of Native American founders

Native Americans derive from a small number of Asian founders who likely arrived to the Americas via Beringia. However, additional details about the intial colonization of the Americas remain unclear. To investigate the pioneering phase in the Americas we analyzed a total of 623 complete mtDNAs from the Americas and Asia, including 20 new complete mtDNAs from the Americas and seven from Asia. This sequence data was used to direct high-resolution genotyping from 20 American and 26 Asian populations. Here we describe more genetic diversity within the founder population than was previously reported. The newly resolved phylogenetic structure suggests that ancestors of Native Americans paused when they reached Beringia, during which time New World founder lineages differentiated from their Asian sister-clades. This pause in movement was followed by a swift migration southward that distributed the founder types all the way to South America. The data also suggest more recent bi-directional gene flow between Siberia and the North American Arctic.

Deciphering host migrations and origins by means of their microbes

Mitochondrial DNA and microsatellite sequences are powerful genetic markers for inferring the genealogy and the population genetic structure of animals but they have only limited resolution for organisms that display low genetic variability due to recent strong bottlenecks. An alternative source of data for deciphering migrations and origins in genetically uniform hosts can be provided by some of their microbes, if their evolutionary history correlates closely with that of the host. In this review, we first discuss how a variety of viruses, and the bacterium Helicobacter pylori, can be used as genetic tracers for one of the most intensively studied species, Homo sapiens. Then, we review statistical problems and limitations that affect the calculation of particular population genetic parameters for these microbes, such as mutation rates, with particular emphasis on the effects of recombination, selection and mode of transmission. Finally, we extend the discussion to other host-parasite systems and advocate the adoption of an integrative approach to both sampling and analysis.

Rapid deployment of the five founding Amerind mtDNA haplogroups via coastal and riverine colonization

Numerous studies of variation in mtDNA in Amerindian populations established that four haplogroups are present throughout both North and South America. These four haplogroups (A, B, C, and D) and perhaps a fifth (X) in North America are postulated to be present in the initial founding migration to the Americas. Furthermore, studies of ancient mtDNA in North America suggested long-term regional continuity of the frequencies of these founding haplogroups. Present-day tribal groups possess high frequencies of private mtDNA haplotypes (variants within the major haplogroups), consistent with early establishment of local isolation of regional populations. Clearly these patterns have implications for the mode of colonization of the hemisphere. Recently, the earlier consensus among archaeologists for an initial colonization by Clovis hunters arriving through an ice-free corridor and expanding in a "blitzkrieg " wave was shown to be inconsistent with extensive genetic variability in Native Americans; a coastal migration route avoids this problem. The present paper demonstrates through a computer simulation model how colonization along coasts and rivers could have rapidly spread the founding lineages widely through North America.

South Amerindian craniofacial morphology: Diversity and implications for Amerindian evolution

The most compelling models concerning the peopling of the Americas consider that modern Amerindians share a common biological pattern, showing affinities with populations of the Asian Northeast. The aim of the present study was to assess the degree of variation of craniofacial morphology of South American Amerindians in a worldwide context. Forty-three linear variables were analyzed on crania derived from American, Asian, Australo-Melanesian, European, South-Saharan African, and Polynesian regions. South America was represented by seven Amerindian samples. In order to understand morphologic diversity among Amerindians of South America, variation was estimated using regions and local populations as units of analysis. Variances and F(ST) values were calculated for each unit, respectively. Both analyses indicated that morphologic variation in Southern Amerindians is extremely high: an F(ST) of 0.01531 was obtained for Southern Amerindians, and values from 0.0371-0.1205 for other world regions. Some aspects linked to the time and mode of the peopling of the Americas and various microevolutionary processes undergone by Amerindians are discussed. Some of the alternatives proposed to explain this high variation include: a greater antiquity of the peopling than what is mostly accepted, a peopling by several highly differentiated waves, an important effect of genetic drift, and gene flow with Paleoamericans. A combination of some of these alternatives explains at least some of the variation.

Identification of Native American founder mtDNAs through the analysis of complete mtDNA sequences: some caveats

In this study, a detailed analysis of both previously published and new data was performed to determine whether complete, or almost complete, mtDNA sequences can resolve the long-debated issue of which Asian mtDNAs were founder sequences for the Native American mtDNA pool. Unfortunately, we now know that coding region data and their analysis are not without problems. To obtain and report reasonably correct sequences does not seem to be a trivial task, and to discriminate between Asian and Native American mtDNA ancestries may be more complex than previously believed. It is essential to take into account the effects of mutational hot spots in both the control and coding regions, so that the number of apparent Native American mtDNA founder sequences is not erroneously inflated. As we report here, a careful analysis of all available data indicates that there is very little evidence that more than five founder mtDNA sequences entered Beringia before the Last Glacial Maximum and left their traces in the current Native American mtDNA pool.

African Haplogroup L mtDNA sequences show violations of clock-like evolution

A set of 96 complete mtDNA sequences that belong to the three major African haplogroups (L1, L2, and L3) was analyzed to determine if mtDNA has evolved as a molecular clock. Likelihood ratio tests (LRTs) were carried out with each of the haplogroups and with combined haplogroup sequence sets. Evolution has not been clock-like, neither for the coding region nor for the control region, in combined sets of African haplogroup L mtDNA sequences. In tests of individual haplogroups, L2 mtDNAs showed violations of a molecular clock under all conditions and in both the control and coding regions. In contrast, haplogroup L1 and L3 sequences, both for the coding and control regions, show clock-like evolution. In clock tests of individual L2 subclades, the L2a sequences showed a marked violation of clock-like evolution within the coding region. In addition, the L2a and L2c branch lengths of both the coding and control regions were shorter relative to those of the L2b and L2d sequences, a result that indicates lower levels of sequence divergence. Reduced median network analyses of the L2a sequences indicated the occurrence of marked homoplasy at multiple sites in the control region. After exclusion of the L2a and L2c sequences, African mtDNA coding region evolution has not significantly departed from a molecular clock, despite the results of neutrality tests that indicate the mitochondrial coding region has evolved under nonneutral conditions. In contrast, control region evolution is clock-like only at the haplogroup level, and it thus appears to have evolved essentially independently from the coding region. The results of the clock tests, the network analyses, and the branch length comparisons all caution against the use of simple mtDNA clocks.