Search

[Societal Impact Statement]: Crop archaeogenomics has rapidly flourished in recent years, leading to a new way of understanding the past and bringing answers to important questions about human history in relation to plant management and food production. Furthermore, the knowledge derived from the analysis of ancient crops can contribute to the development of a more sustainable future. However, the extant legal framework presents a number of challenges when applied to this research field, particularly in the current scenario of disparities in scientific outcomes between countries. We expose the uncertainties of the legal framework and the factors that maintain or exacerbate these inequalities, as well as possible solutions., [Summary]: Crop archaeogenomics is a flourishing field that has greatly benefited from next-generation sequencing technologies. Ancient and historical plant remains are currently considered genetic resources and as such are subject to legal frameworks like those implemented by the Nagoya Protocol. In addition to the challenges in complying with genetic resource regulations that crop archaeogenomics share with other basic plant research disciplines, there are additional difficulties specific to this interdisciplinary field that includes science and humanities, namely, the need to comply with two different legislations before accessing the samples (one for genetic resources and one for cultural heritage), along with a high risk of not obtaining DNA. As a result, most studies to date have been done on samples for which the laws regulating genetic resources did not apply, sometimes avoiding the need of reaching Access and Benefit Sharing agreements with the country that originally provided the samples. This phenomenon is likely to worsen in the future, as the archaeological record is a limited resource and competition between laboratories will only widen the gap between developed and developing economies. Because crop archaeogenomics is a new and promising scientific field, it is desirable to begin a dialogue with other basic biological research fields to facilitate the implementation of these agreements so that basic sciences can easily utilize these biological samples while ensuring the rights of all parties involved., A.I. is a FPI fellow (PRE2018-083529). L. B is a Ramón y Cajal Fellow. (RYC2018-024770-I) both fellowships funded by the Ministerio de Ciencia e Innovación—Agencia Estatal de Investigación/Fondo Social Europeo. We acknowledge financial support from the Spanish Ministry of Science and Innovation-State Research Agency (AEI), through the “Severo Ochoa Programme for Centres of Excellence in R&D” SEV-2015-0533 and CEX2019-000902-S. This work was also supported by the CERCA programme by the Generalitat de Catalunya.

The tree of life is the fundamental biological roadmap for navigating the evolution and properties of life on Earth, and yet remains largely unknown. Even angiosperms (flowering plants) are fraught with data gaps, despite their critical role in sustaining terrestrial life. Today, high-throughput sequencing promises to significantly deepen our understanding of evolutionary relationships. Here, we describe a comprehensive phylogenomic platform for exploring the angiosperm tree of life, comprising a set of open tools and data based on the 353 nuclear genes targeted by the universal Angiosperms353 sequence capture probes. The primary goals of this article are to (i) document our methods, (ii) describe our first data release, and (iii) present a novel open data portal, the Kew Tree of Life Explorer (https://treeoflife.kew.org). We aim to generate novel target sequence capture data for all genera of flowering plants, exploiting natural history collections such as herbarium specimens, and augment it with mined public data. Our first data release, described here, is the most extensive nuclear phylogenomic data set for angiosperms to date, comprising 3099 samples validated by DNA barcode and phylogenetic tests, representing all 64 orders, 404 families (96%⁠) and 2333 genera (17%⁠). A “first pass” angiosperm tree of life was inferred from the data, which totaled 824,878 sequences, 489,086,049 base pairs, and 532,260 alignment columns, for interactive presentation in the Kew Tree of Life Explorer. This species tree was generated using methods that were rigorous, yet tractable at our scale of operation. Despite limitations pertaining to taxon and gene sampling, gene recovery, models of sequence evolution and paralogy, the tree strongly supports existing taxonomy, while challenging numerous hypothesized relationships among orders and placing many genera for the first time. The validated data set, species tree and all intermediates are openly accessible via the Kew Tree of Life Explorer and will be updated as further data become available. This major milestone toward a complete tree of life for all flowering plant species opens doors to a highly integrated future for angiosperm phylogenomics through the systematic sequencing of standardized nuclear markers. Our approach has the potential to serve as a much-needed bridge between the growing movement to sequence the genomes of all life on Earth and the vast phylogenomic potential of the world’s natural history collections. [Angiosperms; Angiosperms353; genomics; herbariomics; museomics; nuclear phylogenomics; open access; target sequence capture; tree of life.], This work was supported by grants from the Calleva Foundation and the Sackler Trust to the Plant and Fungal Trees of Life project at the Royal Botanic Gardens, Kew. Additional support was received from the Garfield Weston Foundation, as part of the Global Tree Seed Bank Programme.

Tetraploid emmer wheat (Triticum turgidum ssp. dicoccon) is a progenitor of the world’s most widely grown crop, hexaploid bread wheat (Triticum aestivum), as well as the direct ancestor of tetraploid durum wheat (T. turgidum subsp. turgidum). Emmer was one of the first cereals to be domesticated in the old world; it was cultivated from around 9700 BC in the Levant1,2 and subsequently in south-western Asia, northern Africa and Europe with the spread of Neolithic agriculture3,4. Here, we report a whole-genome sequence from a museum specimen of Egyptian emmer wheat chaff, 14C dated to the New Kingdom, 1130–1000 BC. Its genome shares haplotypes with modern domesticated emmer at loci that are associated with shattering, seed size and germination, as well as within other putative domestication loci, suggesting that these traits share a common origin before the introduction of emmer to Egypt. Its genome is otherwise unusual, carrying haplotypes that are absent from modern emmer. Genetic similarity with modern Arabian and Indian emmer landraces connects ancient Egyptian emmer with early south-eastern dispersals, whereas inferred gene flow with wild emmer from the Southern Levant signals a later connection. Our results show the importance of museum collections as sources of genetic data to uncover the history and diversity of ancient cereals., M.F.S. and R.M. are supported by RCUK BBSRC grant BB/M011585/1. R.M. is also supported by RCUK BBSRC grant BB/P024726/1; L.R.B. by the Spanish Ministry of Economy and Competitiveness Severo Ochoa Programme for Centres of Excellence in R&D 2016-2019 (SEV-2015-0533) and CERCA Programme, Generalitat de Catalunya; M.G.T. and S.B. by a Wellcome Trust Senior Research Fellowship, grant 100719/Z/12/Z. D.Q.F. and C.S. are supported by the ERC ComPag project, grant number 323842. V.E.M. and S.B. are partially supported by the RCUK NERC grant NE/P012574/1. UCL computing infrastructure was supported by BBSRC grant BB/R01356X/1.

Sequencing of target-enriched libraries is an efficient and cost-effective method for obtaining DNA sequence data from hundreds of nuclear loci for phylogeny reconstruction. Much of the cost of developing targeted sequencing approaches is associated with the generation of preliminary data needed for the identification of orthologous loci for probe design. In plants, identifying orthologous loci has proven difficult due to a large number of whole-genome duplication events, especially in the angiosperms (flowering plants). We used multiple sequence alignments from over 600 angiosperms for 353 putatively single-copy protein-coding genes identified by the One Thousand Plant Transcriptomes Initiative to design a set of targeted sequencing probes for phylogenetic studies of any angiosperm group. To maximize the phylogenetic potential of the probes, while minimizing the cost of production, we introduce a k-medoids clustering approach to identify the minimum number of sequences necessary to represent each coding sequence in the final probe set. Using this method, 5–15 representative sequences were selected per orthologous locus, representing the sequence diversity of angiosperms more efficiently than if probes were designed using available sequenced genomes alone. To test our approximately 80,000 probes, we hybridized libraries from 42 species spanning all higher-order groups of angiosperms, with a focus on taxa not present in the sequence alignments used to design the probes. Out of a possible 353 coding sequences, we recovered an average of 283 per species and at least 100 in all species. Differences among taxa in sequence recovery could not be explained by relatedness to the representative taxa selected for probe design, suggesting that there is no phylogenetic bias in the probe set. Our probe set, which targeted 260 kbp of coding sequence, achieved a median recovery of 137 kbp per taxon in coding regions, a maximum recovery of 250 kbp, and an additional median of 212 kbp per taxon in flanking non-coding regions across all species. These results suggest that the Angiosperms353 probe set described here is effective for any group of flowering plants and would be useful for phylogenetic studies from the species level to higher-order groups, including the entire angiosperm clade itself., This work was supported by the Texas Tech University College of Arts and Sciences [to M.G.J.], the National Science Foundation [DEB-1239992, DEB-1342873 to N.J.W.], and by the Calleva Foundation, the Garfield Weston Foundation, and the Sackler Trust to the Royal Botanic Gardens, Kew.

Providing science and society with an integrated, up-to-date, high quality, open, reproducible and sustainable plant tree of life would be a huge service that is now coming within reach. However, synthesizing the growing body of DNA sequence data in the public domain and disseminating the trees to a diverse audience are often not straightforward due to numerous informatics barriers. While big synthetic plant phylogenies are being built, they remain static and become quickly outdated as new data are published and tree-building methods improve. Moreover, the body of existing phylogenetic evidence is hard to navigate and access for non-experts. We propose that our community of botanists, tree builders, and informaticians should converge on a modular framework for data integration and phylogenetic analysis, allowing easy collaboration, updating, data sourcing and flexible analyses. With support from major institutions, this pipeline should be re-run at regular intervals, storing trees and their metadata long-term. Providing the trees to a diverse global audience through user-friendly front ends and application development interfaces should also be a priority. Interactive interfaces could be used to solicit user feedback and thus improve data quality and to coordinate the generation of new data. We conclude by outlining a number of steps that we suggest the scientific community should take to achieve global phylogenetic synthesis.

A Correction to this paper has been published: https://doi.org/10.1038/s41586-021-03286-9.

This study investigates the influence of genetic differentiation in determining worldwide heterogeneity in osteoporosis-related hip fracture rates. The results indicate that global variation in fracture incidence exceeds that expected on the basis of random genetic variance. Worldwide, the incidence of osteoporotic hip fractures varies considerably. This variability is believed to relate mainly to non-genetic factors. It is conceivable, however, that genetic susceptibility indeed differs across populations. Here, we present the first quantitative assessment of the effects of genetic differentiation on global variability in hip fracture rates. We investigate the observed variance in publically reported age-standardized rates of hip fracture among 28 populations from around the world relative to the expected variance given the phylogenetic relatedness of these populations. The extent to which these variances are similar constitutes a “phylogenetic signal,” which was measured using the K statistic. Population genetic divergence was calculated using a robust array of genome-wide single nucleotide polymorphisms. While phylogenetic signal is maximized when K > 1, a K value of only 0.103 was detected in the combined-sex fracture rate pattern across the 28 populations, indicating that fracture rates vary more than expected based on phylogenetic relationships. When fracture rates for the sexes were analyzed separately, the degree of phylogenetic signal was also found to be small (females: K = 0.102; males: K = 0.081). The lack of a strong phylogenetic signal underscores the importance of factors other than stochastic genetic diversity in shaping worldwide heterogeneity in hip fracture incidence.

Next-generation sequencing technology has facilitated the discovery of millions of genetic variants in human genomes. A sizeable fraction of these variants are predicted to be deleterious. Here, we review the pattern of deleterious alleles as ascertained in genome sequencing data sets and ask whether human populations differ in their predicted burden of deleterious alleles - a phenomenon known as mutation load. We discuss three demographic models that are predicted to affect mutation load and relate these models to the evidence (or the lack thereof) for variation in the efficacy of purifying selection in diverse human genomes. We also emphasize why accurate estimation of mutation load depends on assumptions regarding the distribution of dominance and selection coefficients - quantities that remain poorly characterized for current genomic data sets.

Europe has played a major role in dog evolution, harbouring the oldest uncontested Palaeolithic remains and having been the centre of modern dog breed creation. Here we sequence the genomes of an Early and End Neolithic dog from Germany, including a sample associated with an early European farming community. Both dogs demonstrate continuity with each other and predominantly share ancestry with modern European dogs, contradicting a previously suggested Late Neolithic population replacement. We find no genetic evidence to support the recent hypothesis proposing dual origins of dog domestication. By calibrating the mutation rate using our oldest dog, we narrow the timing of dog domestication to 20,000–40,000 years ago. Interestingly, we do not observe the extreme copy number expansion of the AMY2B gene characteristic of modern dogs that has previously been proposed as an adaptation to a starch-rich diet driven by the widespread adoption of agriculture in the Neolithic., The European continent is thought to have played a major role in the origins of modern dogs. Here, analysing two ancient dog genomes from Germany, the authors find significant genetic continuity throughout the Neolithic period and time dog domestication to ∼20,000–40,000 years ago.

Botigué, Laura R. et al., Human genetic diversity in southern Europe is higher than in other regions of the continent. This difference has been attributed to postglacial expansions, the demic diffusion of agriculture from the Near East, and gene flow from Africa. Using SNP data from 2,099 individuals in 43 populations, we show that estimates of recent shared ancestry between Europe and Africa are substantially increased when gene flow from North Africans, rather than Sub-Saharan Africans, is considered. The gradient of North African ancestry accounts for previous observations of low levels of sharing with Sub-Saharan Africa and is independent of recent gene flow from the Near East. The source of genetic diversity in southern Europe has important biomedical implications; we find that most disease risk alleles from genome-wide association studies follow expected patterns of divergence between Europe and North Africa, with the principal exception of multiple sclerosis., M.H. and C.D.B. were supported by National Institutes of Health Grant 3R01HG003229. L.R.B. and D.C. were supported by Ministerio de Ciencia e Innovación Grant GL2010-14944/BOS and Generalitat de Catalunya Grant 2009SGR1101. C.F. was supported by Instituto de Salud Carlos III Grant PI11/00623. The Spanish National Institute for Bioinformatics supported this project.

The Out-of-Africa (OOA) dispersal ~50,000 years ago is characterized by a series of founder events as modern humans expanded into multiple continents. Population genetics theory predicts an increase of mutational load in populations undergoing serial founder effects during range expansions. To test this hypothesis, we have sequenced full genomes and high-coverage exomes from 7 geographically divergent human populations from Namibia, Congo, Algeria, Pakistan, Cambodia, Siberia and Mexico. We find that individual genomes vary modestly in the overall number of predicted deleterious alleles. We show via spatially explicit simulations that the observed distribution of deleterious allele frequencies is consistent with the OOA dispersal, particularly under a model where deleterious mutations are recessive. We conclude that there is a strong signal of purifying selection at conserved genomic positions within Africa, but that many predicted deleterious mutations have evolved as if they were neutral during the expansion out of Africa. Under a model where selection is inversely related to dominance, we show that OOA populations are likely to have a higher mutation load due to increased allele frequencies of nearly neutral variants that are recessive or partially recessive.

A population sample from Sao Tome e Principe (West Africa) was screened for the G6PD-deficient variants A– (376G/202A), Betica (376G/968C), and Santa Maria (376G/542T). G6PD locus haplotype diversity was also investigated using six intragenic RFLPs (FokI, PvuII, BspHI, PstI, BclI, NlaIII) and a (CTT)n microsatellite 18.61 kb within the G6PD locus. The estimated frequencies of the G6PD*B normal allele, the G6PD*A variant (376G), and the G6PD*A– allele were 0.698, 0.194, and 0.108, respectively. G6PD variants Betica and Santa Maria were not found. Similar levels of microsatellite diversity were found on variants G6PD*B and G6PD*A (H= 0.61 and 0.68, respectively), indicating a similar age for both alleles. All G6PD*A– alleles share the RFLP-microsatellite haplotype ++–+–+/195, the same haplotype described in nearly all the *A– alleles from sub-Saharan, Mexican Mestizo, and Portuguese populations, consistent with a single and recent origin of the G202A mutation on this *A haplotype.

The Levant is a region in the Near East with an impressive record of continuous human existence and major cultural developments since the Paleolithic period. Genetic and archeological studies present solid evidence placing the Middle East and the Arabian Peninsula as the first stepping-stone outside Africa. There is, however, little understanding of demographic changes in the Middle East, particularly the Levant, after the first Out-of-Africa expansion and how the Levantine peoples relate genetically to each other and to their neighbors. In this study we analyze more than 500,000 genome-wide SNPs in 1,341 new samples from the Levant and compare them to samples from 48 populations worldwide. Our results show recent genetic stratifications in the Levant are driven by the religious affiliations of the populations within the region. Cultural changes within the last two millennia appear to have facilitated/maintained admixture between culturally similar populations from the Levant, Arabian Peninsula, and Africa. The same cultural changes seem to have resulted in genetic isolation of other groups by limiting admixture with culturally different neighboring populations. Consequently, Levant populations today fall into two main groups: one sharing more genetic characteristics with modern-day Europeans and Central Asians, and the other with closer genetic affinities to other Middle Easterners and Africans. Finally, we identify a putative Levantine ancestral component that diverged from other Middle Easterners ~23,700-15,500 years ago during the last glacial period, and diverged from Europeans ~15,900-9,100 years ago between the last glacial warming and the start of the Neolithic. © 2013 Haber et al., This study was partly supported by the Lebanese American University and the National Geographic Society (The Genographic Project). CT-S was supported by grant number 098051 from The Wellcome Trust., Haber, Marc et al.

Campbell, Christopher L. et al., North African Jews constitute the second largest Jewish Diaspora group. However, their relatedness to each other; to European, Middle Eastern, and other Jewish Diaspora groups; and to their former North African non-Jewish neighbors has not been well defined. Here, genome-wide analysis of five North African Jewish groups (Moroccan, Algerian, Tunisian, Djerban, and Libyan) and comparison with other Jewish and non-Jewish groups demonstrated distinctive North African Jewish population clusters with proximity to other Jewish populations and variable degrees of Middle Eastern, European, and North African admixture. Two major subgroups were identified by principal component, neighbor joining tree, and identity-by-descent analysis - Moroccan/ Algerian and Djerban/Libyan - that varied in their degree of European admixture. These populations showed a high degree of endogamy and were part of a larger Ashkenazi and Sephardic Jewish group. By principal component analysis, these North African groups were orthogonal to contemporary populations from North and South Morocco, Western Sahara, Tunisia, Libya, and Egypt. Thus, this study is compatible with the history of North African Jews - founding during Classical Antiquity with proselytism of local populations, followed by genetic isolation with the rise of Christianity and then Islam, and admixture following the emigration of Sephardic Jews during the Inquisition., This work was supported in part by the Lewis and Rachel Rudin Foundation, the Iranian-American Jewish Federation of New York, the US–Israel Binational Science Foundation, National Institutes of Health Grant 5 U54 CA121852, and Ruth and Sidney Lapidus. L.R.B. and D.C. were supported by Ministerio de Ciencia e Innovación Grant CGL2010-14944/BOS.

Henn, Brenna M. et al., [Abstract] North African populations are distinct from sub-Saharan Africans based on cultural, linguistic, and phenotypic attributes; however, the time and the extent of genetic divergence between populations north and south of the Sahara remain poorly understood. Here, we interrogate the multilayered history of North Africa by characterizing the effect of hypothesized migrations from the Near East, Europe, and sub-Saharan Africa on current genetic diversity. We present dense, genome-wide SNP genotyping array data (730,000 sites) from seven North African populations, spanning from Egypt to Morocco, and one Spanish population. We identify a gradient of likely autochthonous Maghrebi ancestry that increases from east to west across northern Africa; this ancestry is likely derived from >back-to-Africa> gene flow more than 12,000 years ago (ya), prior to the Holocene. The indigenous North African ancestry is more frequent in populations with historical Berber ethnicity. In most North African populations we also see substantial shared ancestry with the Near East, and to a lesser extent sub-Saharan Africa and Europe. To estimate the time of migration from sub-Saharan populations into North Africa, we implement a maximum likelihood dating method based on the distribution of migrant tracts. In order to first identify migrant tracts, we assign local ancestry to haplotypes using a novel, principal component-based analysis of three ancestral populations. We estimate that a migration of western African origin into Morocco began about 40 generations ago (approximately 1,200 ya); a migration of individuals with Nilotic ancestry into Egypt occurred about 25 generations ago (approximately 750 ya). Our genomic data reveal an extraordinarily complex history of migrations, involving at least five ancestral populations, into North Africa. © 2012 Henn et al., [Author Summary] Proposed migrations between North Africa and neighboring regions have included Paleolithic gene flow from the Near East, an Arabic migration across the whole of North Africa 1,400 years ago (ya), and trans-Saharan transport of slaves from sub-Saharan Africa. Historical records, archaeology, and mitochondrial and Y-chromosome DNA have been marshaled in support of one theory or another, but there is little consensus regarding the overall genetic background of North African populations or their origin and expansion. We characterize the patterns of genetic variation in North Africa using ~730,000 single nucleotide polymorphisms from across the genome for seven populations. We observe two distinct, opposite gradients of ancestry: an east-to-west increase in likely autochthonous North African ancestry and an east-to-west decrease in likely Near Eastern Arabic ancestry. The indigenous North African ancestry may have been more common in Berber populations and appears most closely related to populations outside of Africa, but divergence between Maghrebi peoples and Near Eastern/Europeans likely precedes the Holocene (>12,000 ya). We also find significant signatures of sub-Saharan African ancestry that vary substantially among populations. These sub-Saharan ancestries appear to be a recent introduction into North African populations, dating to about 1,200 years ago in southern Morocco and about 750 years ago into Egypt, possibly reflecting the patterns of the trans-Saharan slave trade that occurred during this period., BMH, SG, and CDB were supported by NIH grant 3R01HG003229. LRB and DC were supported by MCINN grant CGL2010-14944/BOS and Generalitat de Catalunya grant 2009SGR1101. AM-E and AB were supported by NIH grant R01GM083606.

One of the main findings derived from the analysis of the Neandertal genome was the evidence for admixture between Neandertals and non-African modern humans. An alternative scenario is that the ancestral population of non-Africans was closer to Neandertals than to Africans because of ancient population substructure. Thus, the study of North African populations is crucial for testing both hypotheses. We analyzed a total of 780,000 SNPs in 125 individuals representing seven different North African locations and searched for their ancestral/derived state in comparison to different human populations and Neandertals. We found that North African populations have a significant excess of derived alleles shared with Neandertals, when compared to sub-Saharan Africans. This excess is similar to that found in non-African humans, a fact that can be interpreted as a sign of Neandertal admixture. Furthermore, the Neandertal's genetic signal is higher in populations with a local, pre-Neolithic North African ancestry. Therefore, the detected ancient admixture is not due to recent Near Eastern or European migrations. Sub-Saharan populations are the only ones not affected by the admixture event with Neandertals. © 2012 Sánchez-Quinto et al.