Complex Evolutionary History With Extensive Ancestral Gene Flow in an African Primate Radiation.

Understanding the drivers of speciation is fundamental in evolutionary biology, and recent studies highlight hybridization as an important evolutionary force. Using whole-genome sequencing data from 22 species of guenons (tribe Cercopithecini), one of the world's largest primate radiations, we show that rampant gene flow characterizes their evolutionary history and identify ancient hybridization across deeply divergent lineages that differ in ecology, morphology, and karyotypes. Some hybridization events resulted in mitochondrial introgression between distant lineages, likely facilitated by cointrogression of coadapted nuclear variants. Although the genomic landscapes of introgression were largely lineage specific, we found that genes with immune functions were overrepresented in introgressing regions, in line with adaptive introgression, whereas genes involved in pigmentation and morphology may contribute to reproductive isolation. In line with reports from other systems that hybridization might facilitate diversification, we find that some of the most species-rich guenon clades are of admixed origin. This study provides important insights into the prevalence, role, and outcomes of ancestral hybridization in a large mammalian radiation.

The evolution and changing ecology of the African hominid oral microbiome.

The oral microbiome plays key roles in human biology, health, and disease, but little is known about the global diversity, variation, or evolution of this microbial community. To better understand the evolution and changing ecology of the human oral microbiome, we analyzed 124 dental biofilm metagenomes from humans, including Neanderthals and Late Pleistocene to present-day modern humans, chimpanzees, and gorillas, as well as New World howler monkeys for comparison. We find that a core microbiome of primarily biofilm structural taxa has been maintained throughout African hominid evolution, and these microbial groups are also shared with howler monkeys, suggesting that they have been important oral members since before the catarrhine-platyrrhine split ca. 40 Mya. However, community structure and individual microbial phylogenies do not closely reflect host relationships, and the dental biofilms of Homo and chimpanzees are distinguished by major taxonomic and functional differences. Reconstructing oral metagenomes from up to 100 thousand years ago, we show that the microbial profiles of both Neanderthals and modern humans are highly similar, sharing functional adaptations in nutrient metabolism. These include an apparent Homo-specific acquisition of salivary amylase-binding capability by oral streptococci, suggesting microbial coadaptation with host diet. We additionally find evidence of shared genetic diversity in the oral bacteria of Neanderthal and Upper Paleolithic modern humans that is not observed in later modern human populations. Differences in the oral microbiomes of African hominids provide insights into human evolution, the ancestral state of the human microbiome, and a temporal framework for understanding microbial health and disease.

Ecology, Not Host Phylogeny, Shapes the Oral Microbiome in Closely Related Species.

Host-associated microbiomes are essential for a multitude of biological processes. Placed at the contact zone between external and internal environments, the little-studied oral microbiome has important roles in host physiology and health. Here, we investigate the roles of host evolutionary relationships and ecology in shaping the oral microbiome in three closely related gorilla subspecies (mountain, Grauer's, and western lowland gorillas) using shotgun metagenomics of 46 museum-preserved dental calculus samples. We find that the oral microbiomes of mountain gorillas are functionally and taxonomically distinct from the other two subspecies, despite close evolutionary relationships and geographic proximity with Grauer's gorillas. Grauer's gorillas show intermediate bacterial taxonomic and functional, and dietary profiles. Altitudinal differences in gorilla subspecies ranges appear to explain these patterns, suggesting a close connection between dental calculus microbiomes and the environment, likely mediated through diet. This is further supported by the presence of gorilla subspecies-specific phyllosphere/rhizosphere taxa in the oral microbiome. Mountain gorillas show a high abundance of nitrate-reducing oral taxa, which may promote adaptation to a high-altitude lifestyle by modulating blood pressure. Our results suggest that ecology, rather than evolutionary relationships and geographic distribution, shape the oral microbiome in these closely related species.

Isolated Grauer's gorilla populations differ in diet and gut microbiome.

The animal gut microbiome has been implicated in a number of key biological processes, ranging from digestion to behaviour, and has also been suggested to facilitate local adaptation. Yet studies in wild animals rarely compare multiple populations that differ ecologically, which is the level at which local adaptation may occur. Further, few studies simultaneously characterize diet and gut microbiome from the same sample, despite their probable interdependence. Here, we investigate the interplay between diet and gut microbiome in three geographically isolated populations of the critically endangered Grauer's gorilla (Gorilla beringei graueri), which we show to be genetically differentiated. We find population- and social group-specific dietary and gut microbial profiles and covariation between diet and gut microbiome, despite the presence of core microbial taxa. There was no detectable effect of age, and only marginal effects of sex and genetic relatedness on the microbiome. Diet differed considerably across populations, with the high-altitude population consuming a lower diversity of plants compared to low-altitude populations, consistent with plant availability constraining dietary choices. The observed pattern of covariation between diet and gut microbiome is probably a result of long-term social and environmental factors. Our study suggests that the gut microbiome is sufficiently plastic to support flexible food selection and hence contribute to local adaptation.

Genomic analyses show extremely perilous conservation status of African and Asiatic cheetahs (Acinonyx jubatus).

We live in a world characterized by biodiversity loss and global environmental change. The extinction of large carnivores can have ramifying effects on ecosystems like an uncontrolled increase in wild herbivores, which in turn can have knock-on impacts on vegetation regeneration and communities. Cheetahs (Acinonyx jubatus) serve important ecosystem functions as apex predators; yet, they are quickly heading towards an uncertain future. Threatened by habitat loss, human-wildlife conflict and illegal trafficking, there are only approximately 7100 individuals remaining in nature. We present the most comprehensive genome-wide analysis of cheetah phylogeography and conservation genomics to date, assembling samples from nearly the entire current and past species' range. We show that their phylogeography is more complex than previously thought, and that East African cheetahs (A. j. raineyi) are genetically distinct from Southern African individuals (A. j. jubatus), warranting their recognition as a distinct subspecies. We found strong genetic differentiation between all classically recognized subspecies, thus refuting earlier findings that cheetahs show only little differentiation. The strongest differentiation was observed between the Asiatic and all the African subspecies. We detected high inbreeding in the Critically Endangered Iranian (A. j. venaticus) and North-western (A. j. hecki) subspecies, and show that overall cheetahs, along with snow leopards, have the lowest genome-wide heterozygosity of all the big cats. This further emphasizes the cheetah's perilous conservation status. Our results provide novel and important information on cheetah phylogeography that can support evidence-based conservation policy decisions to help protect this species. This is especially relevant in light of ongoing and proposed translocations across subspecies boundaries, and the increasing threats of illegal trafficking.

The complex Y-chromosomal history of gorillas.

Studies of the evolutionary relationships among gorilla populations using autosomal and mitochondrial sequences suggest that male-mediated gene flow may have been important in the past, but data on the Y-chromosomal relationships among the gorilla subspecies are limited. Here, we genotyped blood and noninvasively collected fecal samples from 12 captives and 257 wild male gorillas of known origin representing all four subspecies (Gorilla gorilla gorilla, G. g. diehli, G. beringei beringei, and G. b. graueri) at 10 Y-linked microsatellite loci resulting in 102 unique Y-haplotypes for 224 individuals. We found that western lowland gorilla (G. g. gorilla) haplotypes were consistently more diverse than any other subspecies for all measures of diversity and comprised several genetically distinct groups. However, these did not correspond to geographical proximity and some closely related haplotypes were found several hundred kilometers apart. Similarly, our broad sampling of eastern gorillas revealed that mountain (G. b. beringei) and Grauer's (G. b. graueri) gorilla Y-chromosomal haplotypes did not form distinct clusters. These observations suggest structure in the ancestral population with subsequent mixing of differentiated haplotypes by male dispersal for western lowland gorillas, and postisolation migration or incomplete lineage sorting due to short divergence times for eastern gorillas.

Dental Calculus as a Tool to Study the Evolution of the Mammalian Oral Microbiome.

Dental calculus, the calcified form of the mammalian oral microbial plaque biofilm, is a rich source of oral microbiome, host, and dietary biomolecules and is well preserved in museum and archaeological specimens. Despite its wide presence in mammals, to date, dental calculus has primarily been used to study primate microbiome evolution. We establish dental calculus as a valuable tool for the study of nonhuman host microbiome evolution, by using shotgun metagenomics to characterize the taxonomic and functional composition of the oral microbiome in species as diverse as gorillas, bears, and reindeer. We detect oral pathogens in individuals with evidence of oral disease, assemble near-complete bacterial genomes from historical specimens, characterize antibiotic resistance genes, reconstruct components of the host diet, and recover host genetic profiles. Our work demonstrates that metagenomic analyses of dental calculus can be performed on a diverse range of mammalian species, which will allow the study of oral microbiome and pathogen evolution from a comparative perspective. As dental calculus is readily preserved through time, it can also facilitate the quantification of the impact of anthropogenic changes on wildlife and the environment.

The Genome of the Endangered Dryas Monkey Provides New Insights into the Evolutionary History of the Vervets.

Genomic data can be a powerful tool for inferring ecology, behavior, and conservation needs of highly elusive species, particularly, when other sources of information are hard to come by. Here, we focus on the Dryas monkey (Cercopithecus dryas), an endangered primate endemic to the Congo Basin with cryptic behavior and possibly <250 remaining adult individuals. Using whole-genome sequencing data, we show that the Dryas monkey represents a sister lineage to the vervets (Chlorocebus sp.) and has diverged from them ∼1.4 Ma with additional bidirectional gene flow ∼750,000-∼500,000 years ago that has likely involved the crossing of the Congo River. Together with evidence of gene flow across the Congo River in bonobos and okapis, our results suggest that the fluvial topology of the Congo River might have been more dynamic than previously recognized. Despite the presence of several homozygous loss-of-function mutations in genes associated with sperm mobility and immunity, we find high genetic diversity and low levels of inbreeding and genetic load in the studied Dryas monkey individual. This suggests that the current population carries sufficient genetic variability for long-term survival and might be larger than currently recognized. We thus provide an example of how genomic data can directly improve our understanding of highly elusive species.

The evolution of duplicate gene expression in mammalian organs.

Gene duplications generate genomic raw material that allows the emergence of novel functions, likely facilitating adaptive evolutionary innovations. However, global assessments of the functional and evolutionary relevance of duplicate genes in mammals were until recently limited by the lack of appropriate comparative data. Here, we report a large-scale study of the expression evolution of DNA-based functional gene duplicates in three major mammalian lineages (placental mammals, marsupials, egg-laying monotremes) and birds, on the basis of RNA sequencing (RNA-seq) data from nine species and eight organs. We observe dynamic changes in tissue expression preference of paralogs with different duplication ages, suggesting differential contribution of paralogs to specific organ functions during vertebrate evolution. Specifically, we show that paralogs that emerged in the common ancestor of bony vertebrates are enriched for genes with brain-specific expression and provide evidence for differential forces underlying the preferential emergence of young testis- and liver-specific expressed genes. Further analyses uncovered that the overall spatial expression profiles of gene families tend to be conserved, with several exceptions of pronounced tissue specificity shifts among lineage-specific gene family expansions. Finally, we trace new lineage-specific genes that may have contributed to the specific biology of mammalian organs, including the little-studied placenta. Overall, our study provides novel and taxonomically broad evidence for the differential contribution of duplicate genes to tissue-specific transcriptomes and for their importance for the phenotypic evolution of vertebrates.

Host-derived population genomics data provides insights into bacterial and diatom composition of the killer whale skin.

Recent exploration into the interactions and relationship between hosts and their microbiota has revealed a connection between many aspects of the host's biology, health and associated micro-organisms. Whereas amplicon sequencing has traditionally been used to characterize the microbiome, the increasing number of published population genomics data sets offers an underexploited opportunity to study microbial profiles from the host shotgun sequencing data. Here, we use sequence data originally generated from killer whale Orcinus orca skin biopsies for population genomics, to characterize the skin microbiome and investigate how host social and geographical factors influence the microbial community composition. Having identified 845 microbial taxa from 2.4 million reads that did not map to the killer whale reference genome, we found that both ecotypic and geographical factors influence community composition of killer whale skin microbiomes. Furthermore, we uncovered key taxa that drive the microbiome community composition and showed that they are embedded in unique networks, one of which is tentatively linked to diatom presence and poor skin condition. Community composition differed between Antarctic killer whales with and without diatom coverage, suggesting that the previously reported episodic migrations of Antarctic killer whales to warmer waters associated with skin turnover may control the effects of potentially pathogenic bacteria such as Tenacibaculum dicentrarchi. Our work demonstrates the feasibility of microbiome studies from host shotgun sequencing data and highlights the importance of metagenomics in understanding the relationship between host and microbial ecology.

Population-level assessment of genetic diversity and habitat fragmentation in critically endangered Grauer's gorillas.

OBJECTIVES: The critically endangered Grauer's gorilla (Gorilla beringei graueri) has experienced an estimated 77% population decline within a single generation. Although crucial for informed conservation decisions, there is no clear understanding about population structure and distribution of genetic diversity across the species' highly fragmented range. We fill this gap by studying several core and peripheral Grauer's gorilla populations throughout their distribution range. MATERIALS AND METHODS: We generated genetic profiles for a sampling of an unstudied population of Grauer's gorillas from within the species' core range at 13 autosomal microsatellite loci and combined them with previously published and newly generated data from four other Grauer's gorilla populations, two mountain gorilla populations, and one western lowland gorilla population. RESULTS: In agreement with previous studies, the genetic diversity of Grauer's gorillas is intermediate, falling between western lowland and mountain gorillas. Among Grauer's gorilla populations, we observe lower genetic diversity and high differentiation in peripheral compared with central populations, indicating a strong effect of genetic drift and limited gene flow among small, isolated forest fragments. DISCUSSION: Although genetically less diverse, peripheral populations are frequently essential for the long-term persistence of a species and migration between peripheral and core populations may significantly enrich the overall species genetic diversity. Thus, in addition to central Grauer's gorilla populations from the core of the distribution range that clearly deserve conservation attention, we argue that conservation strategies aiming to ensure long-term species viability should include preserving peripheral populations and enhancing habitat connectivity.

Phylogeny of the jacchus group of Callithrix marmosets based on complete mitochondrial genomes.

OBJECTIVES: Two subgroups make up the marmoset genus Callithrix. The "aurita" group is composed of two species, whereas evolutionary relationships among the four species of the "jacchus" group remain unclear. To uncover these relationships, we first sequenced mitochondrial genomes for C. kuhlii and C. penicillata to complement data available for congeners. We then constructed a phylogenetic tree based on mtDNA heavy chain protein coding genes from several primates to untangle species relationships and estimate divergence times of the jacchus group. MATERIALS AND METHODS: MtDNA genomes of C. kuhlii and C. penicillata were Sanger sequenced. These Callithrix mitogenomes were combined with other publically available primate mtDNA genomes. Phylogenies were produced using maximum likelihood and Bayesian inference. Finally, divergence times within the jacchus group of marmosets were estimated with Bayesian inference. RESULTS: In our phylogenetic tree, C. geoffroyi was the sister to all other jacchus group species, followed by C. kuhlii, while C. jacchus and C. penicillata diverged most recently. Bayesian inference showed that C. jacchus and C. penicillata diverged approximately 0.70 MYA and that the jacchus group radiated approximately 1.30 MYA. DISCUSSION: Callithrix nuclear and mtDNA phylogenies frequently result in polytomies and paraphyly. Here, we present a well-supported phylogenetic tree based on mitochondrial genome sequences, which facilitates the understanding of the divergence of the jacchus marmosets. Our results demonstrate how mitochondrial genomes can enrich Callithrix phylogenetic studies by alleviating some of the difficulties faced by previous mtDNA studies and allow formulation of hypotheses to test further under larger genomic-scale analyses.

Application of PE-RADSeq to the study of genomic diversity and divergence of two Brazilian marmoset species (Callithrix jacchus and C. penicillata).

Callithrix jacchus and C. penicillata are among the smallest anthropoid primates, are highly specialized tree gougers, and largely occupy Brazil's most extreme, semi-arid biomes. However, the underlying genomic factors that underpin the evolution of these species and their unique traits are under-investigated. Additionally, exotic populations of these two species are widely established throughout Brazil and hybridize with threatened native congers. Thus, both genomic and conservation factors call for a better understanding of C. jacchus and C. penicillata evolution. Here, we applied PE-RADseq to characterize genomic variation in these two species, using six C. jacchus and seven C. penicillata individuals. We identified an average of 7,463 and 5,180 SNPs/individual in C. penicillata and C. jacchus, respectively, and also found 1,395 variable sites that were represented in both species. C. penicillata showed overall higher levels of genetic diversity than C. jacchus at the variable sites present in both species. Additionally, among these variable sites, 106 showed relative interspecific divergence levels that were significantly higher than the genome-wide average. We further compared relative and absolute divergence for C. penicillata and C. jacchus between RAD loci associated with the 106 significantly diverged variable sites and all other RAD loci present in both species. The former RAD loci set showed significantly elevated relative and absolute divergence measures in comparison to the latter set. This convergence suggests that C. jacchus and C. penicillata may have diverged under a scenario of gene flow under secondary contact. Here, we demonstrate that RADseq is an efficient method to simultaneously discover and genotype a large number of markers and validate the utility of RADseq for examining Callithrix evolution.

Birth and expression evolution of mammalian microRNA genes.

MicroRNAs (miRNAs) are major post-transcriptional regulators of gene expression, yet their origins and functional evolution in mammals remain little understood due to the lack of appropriate comparative data. Using RNA sequencing, we have generated extensive and comparable miRNA data for five organs in six species that represent all main mammalian lineages and birds (the evolutionary outgroup) with the aim to unravel the evolution of mammalian miRNAs. Our analyses reveal an overall expansion of miRNA repertoires in mammals, with threefold accelerated birth rates of miRNA families in placentals and marsupials, facilitated by the de novo emergence of miRNAs in host gene introns. Generally, our analyses suggest a high rate of miRNA family turnover in mammals with many newly emerged miRNA families being lost soon after their formation. Selectively preserved mammalian miRNA families gradually evolved higher expression levels, as well as altered mature sequences and target gene repertoires, and were apparently mainly recruited to exert regulatory functions in nervous tissues. However, miRNAs that originated on the X chromosome evolved high expression levels and potentially diverse functions during spermatogenesis, including meiosis, through selectively driven duplication-divergence processes. Overall, our study thus provides detailed insights into the birth and evolution of mammalian miRNA genes and the associated selective forces.

Next-generation museomics disentangles one of the largest primate radiations.

Guenons (tribe Cercopithecini) are one of the most diverse groups of primates. They occupy all of sub-Saharan Africa and show great variation in ecology, behavior, and morphology. This variation led to the description of over 60 species and subspecies. Here, using next-generation DNA sequencing (NGS) in combination with targeted DNA capture, we sequenced 92 mitochondrial genomes from museum-preserved specimens as old as 117 years. We infer evolutionary relationships and estimate divergence times of almost all guenon taxa based on mitochondrial genome sequences. Using this phylogenetic framework, we infer divergence dates and reconstruct ancestral geographic ranges. We conclude that the extraordinary radiation of guenons has been a complex process driven by, among other factors, localized fluctuations of African forest cover. We find incongruences between phylogenetic trees reconstructed from mitochondrial and nuclear DNA sequences, which can be explained by either incomplete lineage sorting or hybridization. Furthermore, having produced the largest mitochondrial DNA data set from museum specimens, we document how NGS technologies can "unlock" museum collections, thereby helping to unravel the tree-of-life.

Recent divergences and size decreases of eastern gorilla populations.

Compared with other African apes, eastern gorillas (Gorilla beringei) have been little studied genetically. We used analysis of autosomal DNA genotypes obtained from non-invasively collected faecal samples to estimate the evolutionary histories of the two extant mountain gorilla populations and the closely related eastern lowland gorillas. Our results suggest that eastern lowland gorillas and mountain gorillas split beginning some 10 000 years ago, followed 5000 years ago by the split of the two mountain gorilla populations of Bwindi Impenetrable National Park and the Virungas Massif. All three populations have decreased in effective population size, with particularly substantial 10-fold decreases for the mountain gorillas. These dynamics probably reflect responses to habitat changes resulting from climate fluctuations over the past 20 000 years as well as increasing human influence in this densely populated region in the last several thousand years.

Comparative genomic analyses provide new insights into evolutionary history and conservation genomics of gorillas.

Genome sequencing is a powerful tool to understand species evolutionary history, uncover genes under selection, which could be informative of local adaptation, and infer measures of genetic diversity, inbreeding and mutational load that could be used to inform conservation efforts. Gorillas, critically endangered primates, have received considerable attention and with the recently sequenced Bwindi mountain gorilla population, genomic data is now available from all gorilla subspecies and both mountain gorilla populations. Here, we reanalysed this rich dataset with a focus on evolutionary history, local adaptation and genomic parameters relevant for conservation. We estimate a recent split between western and eastern gorillas of 150,000-180,000 years ago, with gene flow around 20,000 years ago, primarily between the Cross River and Grauer's gorilla subspecies. This gene flow event likely obscures evolutionary relationships within eastern gorillas: after excluding putatively introgressed genomic regions, we uncover a sister relationship between Virunga mountain gorillas and Grauer's gorillas to the exclusion of Bwindi mountain gorillas. This makes mountain gorillas paraphyletic. Eastern gorillas are less genetically diverse and more inbred than western gorillas, yet we detected lower genetic load in the eastern species. Analyses of indels fit remarkably well with differences in genetic diversity across gorilla taxa as recovered with nucleotide diversity measures. We also identified genes under selection and unique gene variants specific for each gorilla subspecies, encoding, among others, traits involved in immunity, diet, muscular development, hair morphology and behavior. The presence of this functional variation suggests that the subspecies may be locally adapted. In conclusion, using extensive genomic resources we provide a comprehensive overview of gorilla genomic diversity, including a so-far understudied Bwindi mountain gorilla population, identify putative genes involved in local adaptation, and detect population-specific gene flow across gorilla species.

Ancient reindeer mitogenomes reveal island-hopping colonisation of the Arctic archipelagos.

Climate warming at the end of the last glacial period had profound effects on the distribution of cold-adapted species. As their range shifted towards northern latitudes, they were able to colonise previously glaciated areas, including remote Arctic islands. However, there is still uncertainty about the routes and timing of colonisation. At the end of the last ice age, reindeer/caribou (Rangifer tarandus) expanded to the Holarctic region and colonised the archipelagos of Svalbard and Franz Josef Land. Earlier studies have proposed two possible colonisation routes, either from the Eurasian mainland or from Canada via Greenland. Here, we used 174 ancient, historical and modern mitogenomes to reconstruct the phylogeny of reindeer across its whole range and to infer the colonisation route of the Arctic islands. Our data shows a close affinity among Svalbard, Franz Josef Land and Novaya Zemlya reindeer. We also found tentative evidence for positive selection in the mitochondrial gene ND4, which is possibly associated with increased heat production. Our results thus support a colonisation of the Eurasian Arctic archipelagos from the Eurasian mainland and provide some insights into the evolutionary history and adaptation of the species to its High Arctic habitat.

Ghost admixture in eastern gorillas.

Archaic admixture has had a substantial impact on human evolution with multiple events across different clades, including from extinct hominins such as Neanderthals and Denisovans into modern humans. In great apes, archaic admixture has been identified in chimpanzees and bonobos but the possibility of such events has not been explored in other species. Here, we address this question using high-coverage whole-genome sequences from all four extant gorilla subspecies, including six newly sequenced eastern gorillas from previously unsampled geographic regions. Using approximate Bayesian computation with neural networks to model the demographic history of gorillas, we find a signature of admixture from an archaic 'ghost' lineage into the common ancestor of eastern gorillas but not western gorillas. We infer that up to 3% of the genome of these individuals is introgressed from an archaic lineage that diverged more than 3 million years ago from the common ancestor of all extant gorillas. This introgression event took place before the split of mountain and eastern lowland gorillas, probably more than 40 thousand years ago and may have influenced perception of bitter taste in eastern gorillas. When comparing the introgression landscapes of gorillas, humans and bonobos, we find a consistent depletion of introgressed fragments on the X chromosome across these species. However, depletion in protein-coding content is not detectable in eastern gorillas, possibly as a consequence of stronger genetic drift in this species.

Historical sampling reveals dramatic demographic changes in western gorilla populations.

BACKGROUND: Today many large mammals live in small, fragmented populations, but it is often unclear whether this subdivision is the result of long-term or recent events. Demographic modeling using genetic data can estimate changes in long-term population sizes while temporal sampling provides a way to compare genetic variation present today with that sampled in the past. In order to better understand the dynamics associated with the divergences of great ape populations, these analytical approaches were applied to western gorillas (Gorilla gorilla) and in particular to the isolated and Critically Endangered Cross River gorilla subspecies (G. g. diehli). RESULTS: We used microsatellite genotypes from museum specimens and contemporary samples of Cross River gorillas to infer both the long-term and recent population history. We find that Cross River gorillas diverged from the ancestral western gorilla population ~17,800 years ago (95% HDI: 760, 63,245 years). However, gene flow ceased only ~420 years ago (95% HDI: 200, 16,256 years), followed by a bottleneck beginning ~320 years ago (95% HDI: 200, 2,825 years) that caused a 60-fold decrease in the effective population size of Cross River gorillas. Direct comparison of heterozygosity estimates from museum and contemporary samples suggests a loss of genetic variation over the last 100 years. CONCLUSIONS: The composite history of western gorillas could plausibly be explained by climatic oscillations inducing environmental changes in western equatorial Africa that would have allowed gorilla populations to expand over time but ultimately isolate the Cross River gorillas, which thereafter exhibited a dramatic population size reduction. The recent decrease in the Cross River population is accordingly most likely attributable to increasing anthropogenic pressure over the last several hundred years. Isolation of diverging populations with prolonged concomitant gene flow, but not secondary admixture, appears to be a typical characteristic of the population histories of African great apes, including gorillas, chimpanzees and bonobos.