Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2021 Feb;590(7844):103-110.
doi: 10.1038/s41586-020-03053-2. Epub 2020 Dec 23.

A genetic history of the pre-contact Caribbean

Affiliations

A genetic history of the pre-contact Caribbean

Daniel M Fernandes et al. Nature. 2021 Feb.

Abstract

Humans settled the Caribbean about 6,000 years ago, and ceramic use and intensified agriculture mark a shift from the Archaic to the Ceramic Age at around 2,500 years ago1-3. Here we report genome-wide data from 174 ancient individuals from The Bahamas, Haiti and the Dominican Republic (collectively, Hispaniola), Puerto Rico, Curaçao and Venezuela, which we co-analysed with 89 previously published ancient individuals. Stone-tool-using Caribbean people, who first entered the Caribbean during the Archaic Age, derive from a deeply divergent population that is closest to Central and northern South American individuals; contrary to previous work4, we find no support for ancestry contributed by a population related to North American individuals. Archaic-related lineages were >98% replaced by a genetically homogeneous ceramic-using population related to speakers of languages in the Arawak family from northeast South America; these people moved through the Lesser Antilles and into the Greater Antilles at least 1,700 years ago, introducing ancestry that is still present. Ancient Caribbean people avoided close kin unions despite limited mate pools that reflect small effective population sizes, which we estimate to be a minimum of 500-1,500 and a maximum of 1,530-8,150 individuals on the combined islands of Puerto Rico and Hispaniola in the dozens of generations before the individuals who we analysed lived. Census sizes are unlikely to be more than tenfold larger than effective population sizes, so previous pan-Caribbean estimates of hundreds of thousands of people are too large5,6. Confirming a small and interconnected Ceramic Age population7, we detect 19 pairs of cross-island cousins, close relatives buried around 75 km apart in Hispaniola and low genetic differentiation across islands. Genetic continuity across transitions in pottery styles reveals that cultural changes during the Ceramic Age were not driven by migration of genetically differentiated groups from the mainland, but instead reflected interactions within an interconnected Caribbean world1,8.

PubMed Disclaimer

Figures

Extended Data Fig. 1:
Extended Data Fig. 1:. Temporal distribution of newly-reported individuals and overview of population structure.
(a) Numbers represent individuals from each site; thick lines denote direct 14C dates (95.4% calibrated confidence intervals); thin lines denote archaeological context dating; grey area identifies the first arrivals of ceramic-users in the Caribbean. Colors and labels are consistent with Fig. 1. (b) PCA plot with ancient individuals shown as solid squares or circles (Archaic- or Ceramic-associated individuals, respectively). Newly-reported individuals are outlined in black, genetic outliers are outlined in red, and individuals with <30,000 SNPs are outlined in blue. Individuals are separated by sub-clades, and three individuals from the site of Cueva Roja (Dominican Republic) who were excluded from clading analysis are labeled “Dominican Cueva Roja Archaic” and colored magenta. Individual PDI009, assessed elsewhere as an outlier, is denoted with an asterisk. Three previously-published ancient Caribbean individuals, are shown as inverted triangles outlined in gray and colored for the sub-clade that encompasses the geographic region with which they are associated. This plot focuses on ancient individuals and does not show some present-day populations; a full plot is provided as Fig. S17. (c) ADMIXTURE analysis best supports K=6 ancestral elements. Newly-reported and co-analyzed individuals are clustered by sub-clade; all newly-reported individuals are identified by a black bar to the side of the plot. The same three previously-published individuals, shown in Extended Data Fig. 1b are included, and three modern-day populations are shown for reference (Suruí, Cabécar, Piapoco).
Extended Data Fig. 2|
Extended Data Fig. 2|. FST distances.
Average pairwise FST distances and standard errors (x100) between (a) clades and (b) sites with more than two unrelated individuals, demonstrating both overall high levels of genetic similarity between the Caribbean_Ceramic sub-clades and the sites composing them, as well as the magnitude of genetic differentiation between those and the groups with Archaic- and Venezuela-related ancestries.
Extended Data Fig. 3:
Extended Data Fig. 3:. Maximum likelihood population tree from allele frequencies using Treemix.
The Caribbean_Ceramic sub-clades are shown on the same branch as modern Arawak-speaking groups (Palikur, Jamamadi). Orange arrows represent admixture events, although observations from other analyses (e.g., qpAdm admixture modeling) suggest that the indicated direction of admixture may be inaccurate (e.g., we believe it is more likely that there is GreaterAntilles_Archaic admixture into Haiti_Ceramic than the reverse scenario; Supplementary Information section 9).
Extended Data Fig. 4:
Extended Data Fig. 4:. Estimated effective population sizes.
(a) Estimates per site are based on ROH blocks 4–20 cM long using a likelihood model (Supplementary Information section 7). Colors as per sub-clades, numbers denote the count of analyzed individuals. Highly consanguineous individuals with a sum of ROH>20 above 50 cM were excluded. (b) Same as (a) but for IBD segments 8–20cM long shared on the X chromosome between all pairs of males. Closely related pairs of individuals with a sum of IBD X>20 above 25 cM were excluded. Numbers denote counts of all remaining pairs. In (a) and (b) points represent maximum likelihood estimate and vertical bars represent 95% CI.
Extended Data Fig. 5:
Extended Data Fig. 5:. Conditional heterozygosity by clade.
Conditional heterozygosity in the ancient Caribbean was similar to that of contemporaneous groups from Peru, except for the Archaic-associated groups and Venezuela_Ceramic. First- and second-degree relatives were excluded from the analysis, including the pair of related individuals representing Haiti_Ceramic. Colored circles represent point estimates (color scheme matching Fig. 1); bars represent three standard errors.
Extended Data Fig. 6:
Extended Data Fig. 6:. Pairwise kinship estimates for all individuals from sites where close relatives were identified using autosomal data.
Dotted lines identify family clusters and inter-site relationships; bottom rows correspond to relationships per individual.
Fig. 1:
Fig. 1:. Geography and significant genetic structure.
(a) Newly-reported data shown as large bordered shapes; co-analyzed data shown as small non-bordered shapes. Asterisk (*) denotes Archaic-associated site of Cueva Roja (excluded due to low-coverage); hash (#) denotes sites with admixed individuals. Andrés is represented as SECoastDR_Ceramic and Dominican_Archaic. Numbers of individuals and temporal distribution in Extended Data Fig. 1a. Map generated with the R package “maps” (R Core Team (2013). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL http://www.R-project.org/). (b) Relationships reconstructed from allele sharing (Supplementary Information section 8). Solid lines connect sub-groupings comprising a larger group; dashed lines represent admixture. Colored boxes represent final sub-clades with the color scheme matching Fig. 1a.
Fig. 2:
Fig. 2:. Genetic affinities of ancient Caribbean people.
(a) Outgroup f3-statistics measuring the relatedness of the clades GreaterAntilles_Archaic, Caribbean_Ceramic, and Venezuela_Ceramic to present-day populations (squares). Map generated with the R package “maps” (R Core Team (2013). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL http://www.R-project.org/). (b) We computed f4(Mbuti, Test; LanguageGroup1Pop, LanguageGroup2Pop) evaluating if each Test sub-clade is more closely related to populations belonging to one language family or another. Points represent the average Z-scores among all populations from each pair of language groups tested; horizontal lines show the range across such comparisons. Vertical lines represent a significance threshold corresponding to a 99.5% CI. (c) Admixture graph modelling of representative ancient Caribbean groupings and select non-Caribbean populations. We fit 12 groups, including the clades LesserAntilles_Ceramic and GreaterAntilles_Archaic, without mixture; the other three Caribbean_Ceramic sub-clades and the clade Venezuela_Ceramic fit as mixtures. The worst Z-score comparing observed to expected f-statistics is |3.6|, which is not significant after correcting for multiple hypothesis testing.
Fig. 3:
Fig. 3:. Estimates of effective population size from shared haplotypes.
Details in Supplementary Information section 7. (a) Number of generations since two chromosomes with a shared segment of a specific size shared a common ancestor, assuming a constant population size N=1000. (b) Average rate of ROH segments in different length bins after excluding highly consanguineous individuals (defined as having a sum of ROH>20 >50cM). (c) Rates of IBD segments shared on the X chromosome between pairs of males within length bins after excluding closely related individuals (defined as sum of IBD X>20 >25cM). For the Ne estimates quoted in the paper we use the pool of 12–20cM segments; for comparisons between the two major clades SECoastDR_Ceramic and EasternGreaterAntilles_Ceramic this gives Ne=3082 (95% CI 1530–8150). In (b) and (c) confidence intervals correspond to one standard deviation (68% coverage) assuming a Poisson distribution in each bin (vertical bars). Point estimates (circles) placed at the center of each 2cM bin, with jitter added for visual separation. Gray lines depict expectations for panmictic populations of various sizes.

Similar articles

  • Origins and genetic legacies of the Caribbean Taino.
    Schroeder H, Sikora M, Gopalakrishnan S, Cassidy LM, Maisano Delser P, Sandoval Velasco M, Schraiber JG, Rasmussen S, Homburger JR, Ávila-Arcos MC, Allentoft ME, Moreno-Mayar JV, Renaud G, Gómez-Carballa A, Laffoon JE, Hopkins RJA, Higham TFG, Carr RS, Schaffer WC, Day JS, Hoogland M, Salas A, Bustamante CD, Nielsen R, Bradley DG, Hofman CL, Willerslev E. Schroeder H, et al. Proc Natl Acad Sci U S A. 2018 Mar 6;115(10):2341-2346. doi: 10.1073/pnas.1716839115. Epub 2018 Feb 20. Proc Natl Acad Sci U S A. 2018. PMID: 29463742 Free PMC article.
  • Reconstructing the population genetic history of the Caribbean.
    Moreno-Estrada A, Gravel S, Zakharia F, McCauley JL, Byrnes JK, Gignoux CR, Ortiz-Tello PA, Martínez RJ, Hedges DJ, Morris RW, Eng C, Sandoval K, Acevedo-Acevedo S, Norman PJ, Layrisse Z, Parham P, Martínez-Cruzado JC, Burchard EG, Cuccaro ML, Martin ER, Bustamante CD. Moreno-Estrada A, et al. PLoS Genet. 2013 Nov;9(11):e1003925. doi: 10.1371/journal.pgen.1003925. Epub 2013 Nov 14. PLoS Genet. 2013. PMID: 24244192 Free PMC article.
  • Genetic continuity and change among the Indigenous peoples of California.
    Nakatsuka N, Holguin B, Sedig J, Langenwalter PE 2nd, Carpenter J, Culleton BJ, García-Moreno C, Harper TK, Martin D, Martínez-Ramírez J, Porcayo-Michelini A, Tiesler V, Villapando-Canchola ME, Valdes Herrera A, Callan K, Curtis E, Kearns A, Iliev L, Lawson AM, Mah M, Mallick S, Micco A, Michel M, Workman JN, Oppenheimer J, Qiu L, Zalzala F, Rohland N, Punzo Diaz JL, Johnson JR, Reich D. Nakatsuka N, et al. Nature. 2023 Dec;624(7990):122-129. doi: 10.1038/s41586-023-06771-5. Epub 2023 Nov 22. Nature. 2023. PMID: 37993721 Free PMC article.
  • Peopling of the Americas as inferred from ancient genomics.
    Willerslev E, Meltzer DJ. Willerslev E, et al. Nature. 2021 Jun;594(7863):356-364. doi: 10.1038/s41586-021-03499-y. Epub 2021 Jun 16. Nature. 2021. PMID: 34135521 Review.
  • Anthropological genetic insights on Caribbean population history.
    Nieves-Colón MA. Nieves-Colón MA. Evol Anthropol. 2022 May;31(3):118-137. doi: 10.1002/evan.21935. Epub 2022 Jan 21. Evol Anthropol. 2022. PMID: 35060661 Review.

Cited by

References

    1. Rouse I. The Tainos: Rise & Decline of the People who Greeted Columbus. (Yale University Press, 1992).
    1. Maggiolo MV La isla de Santo Domingo antes de Colón (Banco Central de la Republica Dominicana, 1993).
    1. Keegan WF & Hofman CL The Caribbean before Columbus. (Oxford University Press, 2017).
    1. Nägele K. et al. Genomic insights into the early peopling of the Caribbean. Science 369, 456–460 (2020). - PubMed
    1. Cook SF & Borah W. The Aboriginal Population of Hispaniola. vol. 1 376–410 (University of California Press, 1971).

Publication types