Estimating divergence times from DNA sequences

Persistent Gene Flow Suggests an Absence of Reproductive Isolation in an African Antelope Speciation Model

African antelope diversity is a globally unique vestige of a much richer world-wide Pleistocene megafauna. Despite this, the evolutionary processes leading to the prolific radiation of African antelopes are not well understood. Here, we sequenced 145 whole genomes from both subspecies of the waterbuck (Kobus ellipsiprymnus), an African antelope believed to be in the process of speciation. We investigated genetic structure and population divergence and found evidence of a mid-Pleistocene separation on either side of the eastern Great Rift Valley, consistent with vicariance caused by a rain shadow along the so-called "Kingdon's Line." However, we also found pervasive evidence of both recent and widespread historical gene flow across the Rift Valley barrier. By inferring the genome-wide landscape of variation among subspecies, we found 14 genomic regions of elevated differentiation, including a locus that may be related to each subspecies' distinctive coat pigmentation pattern. We investigated these regions as candidate speciation islands. However, we observed no significant reduction in gene flow in these regions, nor any indications of selection against hybrids. Altogether, these results suggest a pattern whereby climatically driven vicariance is the most important process driving the African antelope radiation and suggest that reproductive isolation may not set in until very late in the divergence process. This has a significant impact on taxonomic inference, as many taxa will be in a gray area of ambiguous systematic status, possibly explaining why it has been hard to achieve consensus regarding the species status of many African antelopes. Our analyses demonstrate how population genetics based on low-depth whole genome sequencing can provide new insights that can help resolve how far lineages have gone along the path to speciation.

Genetic and evolutionary divergence of harbour seals (Phoca vitulina) in Iliamna Lake, Alaska

Freshwater populations of typically marine species present unique opportunities to investigate biodiversity, evolutionary divergence, and the adaptive potential and niche width of species. A few pinniped species have populations that reside solely in freshwater. The harbour seals inhabiting Iliamna Lake, Alaska constitute one such population. Their remoteness, however, has long hindered scientific inquiry. We used DNA from seal scat and tissue samples provided by Indigenous hunters to screen for mitochondrial DNA and microsatellite variation within Iliamna Lake and eight regions across the Pacific Ocean. The Iliamna seals (i) were substantially and significantly discrete from all other populations ( [Formula: see text]F st-mtDNA = 0.544, [Formula: see text]Φ st - mtDNA = 0.541, [Formula: see text]F st-microsatellites = 0.308), (ii) formed a discrete genetic cluster separate from all marine populations (modal ∆k = 2, PC1 = 14.8%), had (iii) less genetic diversity (Hd, π, H exp), and (iv) higher inbreeding (F) than marine populations. These findings are both striking and unexpected revealing that Iliamna seals have likely been on a separate evolutionary trajectory for some time and may represent a unique evolutionary legacy for the species. Attention must now be given to the selective processes driving evolutionary divergence from harbour seals in marine habitats and to ensuring the future of the Iliamna seal.

Giraffe lineages are shaped by major ancient admixture events

Strong genetic structure has prompted discussion regarding giraffe taxonomy,1,2,3 including a suggestion to split the giraffe into four species: Northern (Giraffa c. camelopardalis), Reticulated (G. c. reticulata), Masai (G. c. tippelskirchi), and Southern giraffes (G. c. giraffa).4,5,6 However, their evolutionary history is not yet fully resolved, as previous studies used a simple bifurcating model and did not explore the presence or extent of gene flow between lineages. We therefore inferred a model that incorporates various evolutionary processes to assess the drivers of contemporary giraffe diversity. We analyzed whole-genome sequencing data from 90 wild giraffes from 29 localities across their current distribution. The most basal divergence was dated to 280 kya. Genetic differentiation, FST, among major lineages ranged between 0.28 and 0.62, and we found significant levels of ancient gene flow between them. In particular, several analyses suggested that the Reticulated lineage evolved through admixture, with almost equal contribution from the Northern lineage and an ancestral lineage related to Masai and Southern giraffes. These new results highlight a scenario of strong differentiation despite gene flow, providing further context for the interpretation of giraffe diversity and the process of speciation in general. They also illustrate that conservation measures need to target various lineages and sublineages and that separate management strategies are needed to conserve giraffe diversity effectively. Given local extinctions and recent dramatic declines in many giraffe populations, this improved understanding of giraffe evolutionary history is relevant for conservation interventions, including reintroductions and reinforcements of existing populations.

Evolutionary plant-pollinator responses to anthropogenic land-use change: impacts on ecosystem services

Agricultural intensification at field and landscape scales, including increased use of agrochemicals and loss of semi-natural habitats, is a major driver of insect declines and other community changes. Efforts to understand and mitigate these effects have traditionally focused on ecological responses. At the same time, adaptations to pesticide use and habitat fragmentation in both insects and flowering plants show the potential for rapid evolution. Yet we lack an understanding of how such evolutionary responses may propagate within and between trophic levels with ensuing consequences for conservation of species and ecological functions in agroecosystems. Here, we review the literature on the consequences of agricultural intensification on plant and animal evolutionary responses and interactions. We present a novel conceptualization of evolutionary change induced by agricultural intensification at field and landscape scales and emphasize direct and indirect effects of rapid evolution on ecosystem services. We exemplify by focusing on economically and ecologically important interactions between plants and pollinators. We showcase available eco-evolutionary theory and plant-pollinator modelling that can improve predictions of how agricultural intensification affects interaction networks, and highlight available genetic and trait-focused methodological approaches. Specifically, we focus on how spatial genetic structure affects the probability of propagated responses, and how the structure of interaction networks modulates effects of evolutionary change in individual species. Thereby, we highlight how combined trait-based eco-evolutionary modelling, functionally explicit quantitative genetics, and genomic analyses may shed light on conditions where evolutionary responses impact important ecosystem services.

African bushpigs exhibit porous species boundaries and appeared in Madagascar concurrently with human arrival

Several African mammals exhibit a phylogeographic pattern where closely related taxa are split between West/Central and East/Southern Africa, but their evolutionary relationships and histories remain controversial. Bushpigs (Potamochoerus larvatus) and red river hogs (P. porcus) are recognised as separate species due to morphological distinctions, a perceived lack of interbreeding at contact, and putatively old divergence times, but historically, they were considered conspecific. Moreover, the presence of Malagasy bushpigs as the sole large terrestrial mammal shared with the African mainland raises intriguing questions about its origin and arrival in Madagascar. Analyses of 67 whole genomes revealed a genetic continuum between the two species, with putative signatures of historical gene flow, variable FST values, and a recent divergence time (<500,000 years). Thus, our study challenges key arguments for splitting Potamochoerus into two species and suggests their speciation might be incomplete. Our findings also indicate that Malagasy bushpigs diverged from southern African populations and underwent a limited bottleneck 1000-5000 years ago, concurrent with human arrival in Madagascar. These results shed light on the evolutionary history of an iconic and widespread African mammal and provide insight into the longstanding biogeographic puzzle surrounding the bushpig's presence in Madagascar.

Major Genetic Risk Factors for Dupuytren's Disease Are Inherited From Neandertals

Dupuytren's disease is characterized by fingers becoming permanently bent in a flexed position. Whereas people of African ancestry are rarely afflicted by Dupuytren's disease, up to ∼30% of men over 60 years suffer from this condition in northern Europe. Here, we meta-analyze 3 biobanks comprising 7,871 cases and 645,880 controls and find 61 genome-wide significant variants associated with Dupuytren's disease. We show that 3 of the 61 loci harbor alleles of Neandertal origin, including the second and third most strongly associated ones (P = 6.4 × 10-132 and P = 9.2 × 10-69, respectively). For the most strongly associated Neandertal variant, we identify EPDR1 as the causal gene. Dupuytren's disease is an example of how admixture with Neandertals has shaped regional differences in disease prevalence.

Distortion of Population Statistics due to the Use of Different Methodological Approaches to the Construction of Genomic DNA Libraries

Several different methods of DNA library preparation for paleogenetic studies are now available. However, the chemical reactions underlying each of them can affect the primary sequence of ancient DNA (aDNA) in the libraries and taint the results of a statistical analysis. In this paper, we compare the results of a sequencing of the aDNA libraries of a Bronze Age sample from burials of the Caucasian burial ground Klady, prepared using three different approaches: (1) shotgun sequencing, (2) strategies for selecting target genomic regions, and (3) strategies for selecting target genomic regions, including DNA pre-treatment with a mixture of uracil-DNA glycosylase (UDG) and endonuclease VIII. The impact of the studied approaches to genomic library preparation on the results of a secondary analysis of the statistical data, namely F4 statistics, ADMIXTURE, and principal component analysis (PCA), was analyzed. It was shown that preparation of genomic libraries without the use of UDG can result in distorted statistical data due to postmortem chemical modifications of the aDNA. This distortion can be alleviated by analyzing only the single nucleotide polymorphisms caused by transversions in the genome.

Warthog Genomes Resolve an Evolutionary Conundrum and Reveal Introgression of Disease Resistance Genes

African wild pigs have a contentious evolutionary and biogeographic history. Until recently, desert warthog (Phacochoerus aethiopicus) and common warthog (P. africanus) were considered a single species. Molecular evidence surprisingly suggested they diverged at least 4.4 million years ago, and possibly outside of Africa. We sequenced the first whole-genomes of four desert warthogs and 35 common warthogs from throughout their range. We show that these two species diverged much later than previously estimated, 400,000–1,700,000 years ago depending on assumptions of gene flow. This brings it into agreement with the paleontological record. We found that the common warthog originated in western Africa and subsequently colonized eastern and southern Africa. During this range expansion, the common warthog interbred with the desert warthog, presumably in eastern Africa, underlining this region’s importance in African biogeography. We found that immune system–related genes may have adaptively introgressed into common warthogs, indicating that resistance to novel diseases was one of the most potent drivers of evolution as common warthogs expanded their range. Hence, we solve some of the key controversies surrounding warthog evolution and reveal a complex evolutionary history involving range expansion, introgression, and adaptation to new diseases.

Philippine Ayta possess the highest level of Denisovan ancestry in the world

Multiple lines of evidence show that modern humans interbred with archaic Denisovans. Here, we report an account of shared demographic history between Australasians and Denisovans distinctively in Island Southeast Asia. Our analyses are based on ∼2.3 million genotypes from 118 ethnic groups of the Philippines, including 25 diverse self-identified Negrito populations, along with high-coverage genomes of Australopapuans and Ayta Magbukon Negritos. We show that Ayta Magbukon possess the highest level of Denisovan ancestry in the world-∼30%-40% greater than that of Australians and Papuans-consistent with an independent admixture event into Negritos from Denisovans. Together with the recently described Homo luzonensis, we suggest that there were multiple archaic species that inhabited the Philippines prior to the arrival of modern humans and that these archaic groups may have been genetically related. Altogether, our findings unveil a complex intertwined history of modern and archaic humans in the Asia-Pacific region, where distinct Islander Denisovan populations differentially admixed with incoming Australasians across multiple locations and at various points in time.

Reconstructing 50,000 years of human history from our DNA: lessons from modern genomics

The advent of high throughput sequencing approaches and ancient DNA techniques have enabled reconstructing the history of human populations at an unprecedented level of resolution. The symposium from the French Academy of Sciences "50,000 ans d'épopée humaine dans notre ADN" has reviewed some of the latest contributions from the fields of genomics, archaeology, and linguistics to our understanding of >300,000 years of human history. DNA has revealed the richness of the human journey, from the deep divergences between human populations in Africa, to the first encounters of Homo Sapiens with other hominins on their way to Eurasia and the peopling of Remote Oceania. The symposium has also emphasized how migrations, cultural practices, and environmental pathogens have contributed to shape the genetic diversity of modern humans, through admixture, genetic drift or genetic adaptation. Finally, special attention was also given to how human behaviours have shaped the genome of other species, through the spreading of microbes and pathogens, as in the case of Yersinia Pestis, or through domestication, as elegantly demonstrated for dogs, horses, and apples. Altogether, this conference illustrated how the complex history of human populations is tightly linked with their contemporary genetic diversity that, in turn, has direct effects on their identity and health.

Estimation of coalescence probabilities and population divergence times from SNP data

We present a method called the G(A|B) method for estimating coalescence probabilities within population lineages from genome sequences when one individual is sampled from each population. Population divergence times can be estimated from these coalescence probabilities if additional assumptions about the history of population sizes are made. Our method is based on a method presented by Rasmussen et al. (2014) to test whether an archaic genome is from a population directly ancestral to a present-day population. The G(A|B) method does not require distinguishing ancestral from derived alleles or assumptions about demographic history before population divergence. We discuss the relationship of our method to two similar methods, one introduced by Green et al. (2010) and called the F(A|B) method and the other introduced by Schlebusch et al. (2017) and called the TT method. When our method is applied to individuals from three or more populations, it provides a test of whether the population history is treelike because coalescence probabilities are additive on a tree. We illustrate the use of our method by applying it to three high-coverage archaic genomes, two Neanderthals (Vindija and Altai) and a Denisovan.

Multiple migrations to the Philippines during the last 50,000 years

A key link to understand human history in Island Southeast Asia is the Philippine archipelago and its poorly investigated genetic diversity. We analyzed the most comprehensive set of population-genomic data for the Philippines: 1,028 individuals covering 115 indigenous communities. We demonstrate that the Philippines were populated by at least five waves of human migration. The Cordillerans migrated into the Philippines prior to the arrival of rice agriculture, where some remain as the least admixed East Asians carrying an ancestry shared by all Austronesian-speaking populations, thereby challenging an exclusive out-of-Taiwan model of joint farming–language–people dispersal. Altogether, our findings portray the Philippines as a crucial gateway, with a multilayered history, that ultimately changed the genetic landscape of the Asia-Pacific region.

Island Southeast Asia has recently produced several surprises regarding human history, but the region’s complex demography remains poorly understood. Here, we report ∼2.3 million genotypes from 1,028 individuals representing 115 indigenous Philippine populations and genome-sequence data from two ∼8,000-y-old individuals from Liangdao in the Taiwan Strait. We show that the Philippine islands were populated by at least five waves of human migration: initially by Northern and Southern Negritos (distantly related to Australian and Papuan groups), followed by Manobo, Sama, Papuan, and Cordilleran-related populations. The ancestors of Cordillerans diverged from indigenous peoples of Taiwan at least ∼8,000 y ago, prior to the arrival of paddy field rice agriculture in the Philippines ∼2,500 y ago, where some of their descendants remain to be the least admixed East Asian groups carrying an ancestry shared by all Austronesian-speaking populations. These observations contradict an exclusive “out-of-Taiwan” model of farming–language–people dispersal within the last four millennia for the Philippines and Island Southeast Asia. Sama-related ethnic groups of southwestern Philippines additionally experienced some minimal South Asian gene flow starting ∼1,000 y ago. Lastly, only a few lowlanders, accounting for <1% of all individuals, presented a low level of West Eurasian admixture, indicating a limited genetic legacy of Spanish colonization in the Philippines. Altogether, our findings reveal a multilayered history of the Philippines, which served as a crucial gateway for the movement of people that ultimately changed the genetic landscape of the Asia-Pacific region.