Disjunct distributions of freshwater snails testify to a central role of the Congo system in shaping biogeographical patterns in Africa

Phylogenomic analyses reveal incongruences between divergence times and fossil records of freshwater snails in East Asia

Fossils provide important insight into our understanding of phylogenetic history by serving as calibration points for divergence time estimation. However, uncertainties in the fossil record due to parallel evolution and convergent evolution can critically affect estimates of node ages. Here, we compare and contrast estimates of phylogenetic divergence with geologic and fossil history for two freshwater snail genera of the family Viviparidae in East Asia (Cipangopaludina and Margarya). Cipangopaludina species are commonly widely distributed species in East Asia, but extant Margarya species are endemic to the ancient lakes in Yunnan, China. According to some previous studies, parallel evolution or convergent evolution of shell morphology has occurred in the family several times which may affect divergence time estimation using fossil records. In this study, we used SNP data derived from ddRAD-seq loci to investigate population demographic history of both genera. Our results show a common pattern of lake endemic lineages diversifying from widely distributed lineages in the Miocene, and multiple colonization to a single ancient lake occurred in the Pleistocene. Our results indicate substantial incongruence among estimated phylogenomic divergence times, some fossil records, and formation ages of ancient lakes. These findings suggest some fossil records may be misidentified in these groups and highlight the need to carefully evaluate geological evidence and fossil records when using these for divergence time estimation.

Population Genomic Evidence for the Diversification of Bellamya aeruginosa in Different River Systems in China

Clarifying the genetic structure can facilitate the understanding of a species evolution history. It is crucial for the management of germplasm resources and providing useful guidance for effective selective breeding. Bellamya is an economically and ecologically important freshwater snail for fish, birds and even humans. Population genetic structures of the Bellamya species, however, were unknown in previous studies. Population genomics approaches with tens to hundreds of thousands of single nucleotide polymorphisms (SNPs) make it possible to detect previously unidentified structures. The population genomic study of seven populations of B. aeruginosa across three river systems (Yellow River, Yangtze River and Pearl River) in China was conducted by SLAF-seq. SLAF-seq obtained a total of 4737 polymorphisms SLAF-tags and 25,999 high-consistency genome-wide SNPs. The population genetic structure showed a clear division among populations from the Yellow River basin (YH and WL) and the Pearl River basin (QSH and LB), as well as population YC from the Yangtze River basin using the SNPs data. However, there existed no distinct population structure using the mitochondrial DNA (mtDNA). The anthropogenic translocation from the Yangtze River basin to the Pearl River basin and the passive dispersion from the Yangtze River basin to the Yellow River basin by flooding have weakened the phylogeographic pattern of B. aeruginosa. The divergence of B. aeruginosa in the three river systems suggests that the anthropogenic dispersal for aquaculture and breeding requires serious consideration of the population structure for the preservation of genetic diversity and effective utilization of germplasm resources.

Reanalysis on Phylogeographic Pattern of Sharpbelly Hemiculter leucisculus (Cyprinidae: Cultrinae) in China: A Review and the Implications for Conservation

Hemiculter leucisclus, as a widely distributed freshwater fish in China, provides an interesting model to explore the impact of drainage evolution and geologic history in the Pleistocene on diversification patterns. We collected the mitochondrial cytochrome b (cytb) gene and the recombination activating gene 2 (RAG2) from 1,070 individuals from 59 sampling locations. Phylogenetic and population genetic approaches were used to describe the phylogeographic pattern and to test how the geological and climatic factors on diversification. The results suggested that there existed four sublineages of the H. leucisclus across six river systems, among which two sublineages, showing strongly indigenous characteristics, are constrained to particular geographical regions in China. The molecular data and ancestral states demonstrated that the H. leucisclus possibly originated from the Pearl River basin during the later Pliocene. The phylogeographic pattern in H. leucisclus appears to have been driven by palaeoenvironmental perturbations rather than anthropogenic translocations. The geographically constrained sublineages A in the middle and lower Pearl River basin and sublineage B in the upper Yangtze River basin deserves special protection.

Patterns of variation of mutation rates of mitochondrial and nuclear genes of gastropods

BACKGROUND

Although mitochondrial DNA (mtDNA) of many animals tends to mutate at higher rates than nuclear DNA (nuDNA), a recent survey of mutation rates of various animal groups found that the gastropod family Bradybaenidae (suborder Helicina) shows a nearly 40-fold difference in mutation rates of mtDNA ([Formula: see text]_m) and nuDNA ([Formula: see text]_n), while other gastropod taxa exhibit only two to five-fold differences. To determine if Bradybaenidae represents an outlier within Gastropoda, I compared estimated values of [Formula: see text]_m/[Formula: see text]_n of additional gastropod groups. In particular, I reconstructed mtDNA and nuDNA gene trees of 121 datasets that include members of various clades contained within the gastropod subclasses Caenogastropoda, Heterobranchia, Patellogastropoda, and Vetigastropoda and then used total branch length estimates of these gene trees to infer [Formula: see text]_m/[Formula: see text]_n.

RESULTS

Estimated values of [Formula: see text]_m/[Formula: see text]_n range from 1.4 to 91.9. Datasets that exhibit relatively large values of [Formula: see text]_m/[Formula: see text]_n (i.e., > 20), however, show relatively lower estimates of [Formula: see text]_n (and not elevated [Formula: see text]_m) in comparison to groups with lower values. These datasets also tend to contain sequences of recently diverged species. In addition, datasets with low levels of phylogenetic breadth (i.e., contain members of single genera or families) exhibit higher values of [Formula: see text]_m/[Formula: see text]_n than those with high levels (i.e., those that contain representatives of single superfamilies or higher taxonomic ranks).

CONCLUSIONS

Gastropods exhibit considerable variation in estimates of [Formula: see text]_m/[Formula: see text]_n. Large values of [Formula: see text]_m/[Formula: see text]_n that have been calculated for Bradybaenidae and other gastropod taxa may be overestimated due to possible sampling artifacts or processes that depress estimates of total molecular divergence of nuDNA in groups that recently diversified.

DNA barcoding of Iranian radicine freshwater snails begins to untangle the taxonomy and phylogeography of intermediate hosts of schistosomiasis and fasciolosis from the Middle East and across Central Asia

In the Middle East radicine snails are of considerable medical and veterinary importance acting as vectors of trematodes. In Iran, such snails are responsible for the transmission of the zoonotic trematodes Schistosoma turkestanicum and Fasciola gigantica. Historically, Radix gedrosiana has been incriminated as an important intermediate host for both trematodes, however, controversy remains over the snail's true taxonomic status. This species has been determined using morphological characters that has resulted in erroneous identification of species, affecting understanding of population biology, and ultimately affecting vector incrimination. In this current study DNA barcoding using cox1 and phylogenetic analyses revealed that snails identified as R. gedrosiana from Iran split into two separate species, Radix euphratica and Ampullaceana sp. The cox1 also provided useful insights into the evolutionary history of R. euphratica populations. Phylogeographic analyses indicated that R. euphratica had an Iraqi/Iranian origin approximately 3.3 MYA and exists as a large stable population across the Middle East and Central Asia, and a lack of genetic differentiation between geographical isolates. Such molecular barcoding techniques are crucial for the identification of radicine snails of Iran being invaluable for the monitoring of zoonotic flukes, understanding the distribution of infection and the accurate incrimination of snail vectors.

Global Diversification Dynamics Since the Jurassic: Low Dispersal and Habitat-Dependent Evolution Explain Hotspots of Diversity and Shell Disparity in River Snails (Viviparidae)

The Viviparidae, commonly known as River Snails, is a dominant group of freshwater snails with a nearly worldwide distribution that reaches its highest taxonomic and morphological diversity in Southeast Asia. The rich fossil record is indicative of a probable Middle Jurassic origin on the Laurasian supercontinent where the group started to diversify during the Cretaceous. However, it remains uncertain when and how the biodiversity hotspot in Southeast Asia was formed. Here, we used a comprehensive genetic data set containing both mitochondrial and nuclear markers and comprising species representing 24 out of 28 genera from throughout the range of the family. To reconstruct the spatiotemporal evolution of viviparids on a global scale, we reconstructed a fossil-calibrated phylogeny. We further assessed the roles of cladogenetic and anagenetic events in range evolution. Finally, we reconstructed the evolution of shell features by estimating ancestral character states to assess whether the appearance of sculptured shell morphologies was driven by major habitat shifts. The molecular phylogeny supports the monophyly of the three subfamilies, the Bellamyinae, Lioplacinae, and Viviparinae, but challenges the currently accepted genus-level classification in several cases. The almost global distribution of River Snails has been influenced both by comparatively ancient vicariance and more recent founder events. In Southeast Asia, Miocene dispersal was a main factor in shaping the modern species distributions. A recurrent theme across different viviparid taxa is that many species living in lentic waters exhibit sculptured shells, whereas only one strongly sculptured species is known from lotic environments. We show that such shell sculpture is habitat-dependent and indeed evolved several times independently in lentic River Snails. Considerably high transition rates between shell types in lentic habitats probably caused the co-occurrence of morphologically distinct shell types in several lakes. In contrast, directional evolution toward smooth shells in lotic habitats, as identified in the present analyses, explains why sculptured shells are rarely found in these habitats. However, the specific factors that promoted changes in shell morphology require further work. [biogeographical analyses; fossil-calibrated phylogeny; fossil-constrained analyses; Southeast Asia; stochastic character mapping.].

Population collapse in viviparid gastropods of the Lake Victoria ecoregion started before the Last Glacial Maximum

Ecosystems of Lake Victoria and riparian communities have been strongly disrupted by the introduction of the invasive Nile perch and its fishing industry. Beyond this invasion and other recent anthropogenic stressors, the Lake Victoria ecoregion also underwent phases of pronounced aridity over the Late Pleistocene, lastly during the Last Glacial Maximum (LGM). The consequences of recent and historic environmental change have been canvassed for the adaptive radiation of haplochromine cichlids occupying the ecoregion, but their effect on freshwater invertebrate diversity remains largely unknown. Here, we use 15 microsatellite loci and approximate Bayesian computation to test whether viviparid gastropods experienced a population bottleneck during the LGM, as did cichlids. Clustering analyses support three viviparid gene pools in the Lake Victoria ecoregion, gathering specimens from 1) Lake Albert and the White Nile, 2) the Victoria Nile and Lake Kyoga and 3) Lake Victoria and tributaries. The last group contains the highest genetic diversity, but all groups have a considerable number of private alleles and are inferred to predate the LGM. Examinations of demographic history reveal a 190- to 500-fold population decline that started ~ 125-150 ka ago, thus substantially before the LGM bottleneck documented in haplochromine cichlids. Population collapses in viviparids are an order of magnitude more severe than declines in cichlids and have not been halted by the re-establishment of freshwater ecosystems since the LGM. Recent anthropogenic ecosystem deterioration is causing homogenization of previously diversified microhabitats, which may contribute to (local) extinction and enhanced gene flow among species within gene pools.

Checklist of the fresh and brackish water snails (Mollusca, Gastropoda) of Bénin and adjacent West African ecoregions

Currently no comprehensive checklist of fresh and brackish water gastropods from Bénin exists, and those for adjacent West African areas are outdated. Yet, such checklists provide essential biodiversity information and a consistent taxonomic and nomenclatural framework for that biodiversity. Here a first checklist of the fresh and brackish water gastropods from Bénin and adjacent West African ecoregions is presented, based on an extensive literature review and field surveys between September 2014 and June 2019 in six major fresh and brackish water ecosystems in Bénin. This inventory includes information on synonymy, species distribution in West Africa, habitats, and conservation status. The fresh and brackish water gastropod fauna includes 60 species, belonging to 28 genera and 16 families. Pachychilidae, Ampullariidae, Neritidae, and Bulinidae were the most diverse families with 9, 8, 7, and 7 species, respectively. However, literature and field data indicated that 23 species observed in West African basins that extend to Bénin do not occur in the territory of Bénin. These species were not detected in our field surveys, most likely because they are rare at collecting sites. Of the 60 species included, five are classified as “Data Deficient”, 43 as “Least Concern”, two as “Nearly Threatened”, one as “Vulnerable”, and six as “Endangered” by the IUCN, whereas the remaining three species were not evaluated. Because the taxonomy of fresh and brackish water gastropods in West Africa is still largely based on morphology, comparative molecular and taxonomic studies may result in substantial revisions of this checklist over the coming years.

Where ichthyofaunal provinces meet: the fish fauna of the Lake Edward system, East Africa

Based on literature, museum collections and three recent expeditions, an annotated species list of the Lake Edward, East Africa, drainage system is presented, excluding the endemic haplochromines. A total of 34 non-Haplochromis species belonging to 10 families and 21 genera are recorded from the system. Three of these are endemic and two others have been introduced in the region. Six species are new records for the Lake Edward system. A species accumulation curve indicates that we probably covered most of the non-Haplochromis species in the area sampled during the recent expeditions, but undetected species might still be present in the Congolese part of the system, which is poorly sampled. A comparison of the species list with those of neighbouring basins confirmed the placement of the Lake Edward system within the east-coast ichthyofaunal province.

Diversification dynamics of freshwater bivalves (Unionidae: Parreysiinae: Coelaturini) indicate historic hydrographic connections throughout the East African Rift System

Invertebrates are exceptionally diverse, but many are in decline because of anthropogenic changes to their habitat. This situation is particularly problematic for taxa that are not well monitored or taxonomically poorly understood, because the lack of knowledge hampers conservation. Despite their important functional role in freshwater ecosystems, African bivalves of the family Unionidae remain poorly studied compared to their highly threatened relatives in Europe, the U.S.A. and Canada. To resolve relationships and to study diversification dynamics in space and time, we performed time-calibrated phylogenetic studies and biogeographical modeling on the unionids from the East African Rift System and surroundings, including representatives of all currently recognized Afrotropical genera except for Brazzaea (and Unio from southern Africa). Our analyses indicate that all sampled Afrotropical unionids belong to the tribe Coelaturini (subfamily Parreysiinae), as does the genus Moncetia from Lake Tanganyika, which is currently attributed to the family Iridinidae. Colonization of Africa from Eurasia by Parreysiinae occurred ~17 Ma ago, and the subsequent diversification of Coelaturini in Africa continued at a steady pace, although net diversification decreased over time as more niches and ecoregions became occupied. Clades in Coelaturini largely reflect drainage basins, with the oldest lineages and highest regional diversity occurring in Lake Tanganyika, followed by the Congo Basin watershed in general. The species assemblage of Lake Tanganyika reflects multiple independent events of colonization and intralacustrine diversification since the Late Miocene or Early Pliocene. The clades of other regions, including that containing the species from Lake Malawi, are comparatively young. Biogeographical analyses indicate that the colonization history was mainly driven by cladogenesis in sympatry, whereas few anagenetic events contributed to the modern distribution of Coelaturini. Ancestral range estimations demonstrate that Coelaturini originated in the Victoria and/or Tanganyika ecoregions, and that the Congo Basin played an essential role in the colonization of Africa by Coelaturini.

Enigmatic incongruence between mtDNA and nDNA revealed by multi-locus phylogenomic analyses in freshwater snails

Phylogenetic incongruence has frequently been encountered among different molecular markers. Recent progress in molecular phylogenomics has provided detailed and important information for evolutionary biology and taxonomy. Here we focused on the freshwater viviparid snails (Cipangopaludina chinensis chinensis and C. c. laeta) of East Asia. We conducted phylogenetic analyses and divergence time estimation using two mitochondrial markers. We also performed population genetic analyses using genome-wide SNPs. We investigated how and which phylogenetic patterns reflect shell morphology. The results showed these two species could be separated into four major mitochondrial clades, whereas the nuclear clusters supported two groups. The phylogenetic patterns of both mtDNA and nDNA largely reflected the geographical distribution. Shell morphology reflected the phylogenetic clusters based on nDNA. The findings also showed these two species diversified in the Pliocene to early Pleistocene era, and occurred introgressive hybridisation. The results also raise the taxonomic issue of the two species.

Phylogeography of Bellamya (Mollusca: Gastropoda: Viviparidae) snails on different continents: contrasting patterns of diversification in China and East Africa

BACKGROUND

Species diversity is determined by both local environmental conditions that control differentiation and extinction and the outcome of large-scale processes that affect migration. The latter primarily comprises climatic change and dynamic landscape alteration. In the past few million years, both Southeast Asia and Eastern Africa experienced drastic climatic and geological oscillations: in Southeast Asia, especially in China, the Tibetan Plateau significantly rose up, and the flow of the Yangtze River was reversed. In East Africa, lakes and rivers experienced frequent range expansions and regressions due to the African mega-droughts. To test how such climatic and geological histories of both regions relate to their respective regional species and genetic diversity, a large scale comparative phylogeographic study is essential. Bellamya, a species rich freshwater snail genus that is widely distributed across China and East Africa, represents a suitable model system to address this question. We sequenced mitochondrial and nuclear DNA for members of the genus from China and used published sequences from Africa and some other locations in Asia to investigate their phylogeny and distribution of genetic diversity.

RESULTS

Our phylogenetic analysis revealed two monophyletic groups, one in China and one in East Africa. Within the Chinese group, Bellamya species show little genetic differentiation. In contrast, we observe fairly deep divergence among the East African lakes with almost every lake possessing its unique clade. Our results show that strong divergence does not necessarily depend on intrinsic characteristics of a species, but rather is related to the landscape dynamics of a region.

CONCLUSION

Our phylogenetic results suggest that the Bellamya in China and East Africa are independent phylogenetic clades with different evolutionary trajectories. The different climate and geological histories likely contributed to the diverging evolutionary patterns. Repeated range expansions and regressions of lakes likely contributed to the great divergence of Bellamya in East Africa, while reversal of the river courses and intermingling of different lineages had an opposite effect on Bellamya diversification in China.

Phylogeography of freshwater planorbid snails reveals diversification patterns in Eurasian continental islands

Background

Islands have traditionally been the centre of evolutionary biological research, but the dynamics of immigration and differentiation at continental islands have not been well studied. Therefore, we focused on the Japanese archipelago, the continental islands located at the eastern end of the Eurasian continent. While the Japanese archipelago is characterised by high biodiversity and rich freshwater habitats, the origin and formation mechanisms of its freshwater organisms are not clear. In order to clarify the history of the planorbid gastropod fauna, we conducted phylogenetic analysis, divergence time estimation, ancestral state reconstruction, and lineage diversity estimations.

Results

Our analyses revealed the formation process of the planorbid fauna in the Japanese archipelago. Most lineages in the Japanese archipelago have closely related lineages on the continent, and the divergence within the Japanese lineages presumably occurred after the late Pliocene. In addition, each lineage is characterised by different phylogeographical patterns, suggesting that immigration routes from the continent to the Japanese archipelago differ among lineages. Furthermore, a regional lineage diversity plot showed that the present diversity in the Japanese archipelago potentially reflects the differentiation of lineages within the islands after the development of the Japanese archipelago.

Conclusions

Although additional taxon sampling and genetic analysis focused on each lineage are needed, our results suggest that immigration from multiple routes just prior to the development of the Japanese archipelago and subsequent diversification within the islands are major causes of the present-day diversity of the Japanese planorbid fauna.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1273-3) contains supplementary material, which is available to authorized users.

Analysis of genetic diversity and population structure of Bellamya quadrata from lakes of middle and lower Yangtze River

As an endemic species of freshwater gastropods in China, Bellamya quadrata plays an important role in ecosystem service provision and commercial importance. However, the species is overharvested and its natural habitats are under severe threat due to fragmentation and loss. To estimate the genetic diversity and population structure of B. quadrata, 285 individuals from eight lake populations across middle and lower Yangtze River were sampled. Seven microsatellite loci were genotyped. Our results showed that (i) the genetic diversity of B. quadrata was high in most of the studied populations, yet effective population sizes appear to be rather small in some populations; (ii) low levels of genetic differentiation exists among populations but gene flow was generally high; (iii) no clear geographic or genetic structure was observed in the studied region, implying mechanisms (zoochoric dispersal and anthropogenic translocations) that enhance dispersal and gene flow have promoted population connectivity. However, the comparatively high genetic diversity of B. quadrata could be attributed to a lag phase, suggesting that the genetic diversity of this species may be lost in the future and the priorities for conservation of B. quadrata are necessary.

Recurrent camouflaged invasions and dispersal of an Asian freshwater gastropod in tropical Africa

BACKGROUND

Non-indigenous taxa currently represent a large fraction of the species and biomass of freshwater ecosystems. The accumulation of invasive taxa in combination with other stressors in these ecosystems may alter the habitats to which native taxa are adapted, which could elicit evolutionary changes in native populations and their ecological interactions. Assessing ecological and evolutionary consequences of invasions simultaneously may therefore be the most effective approach to study taxa with complex invasion histories. Here we apply such an integrated approach to the cerithioid gastropod Melanoides tuberculata, a model system in invasion biology.

RESULTS

Molecular phylogenetics and ancestral range reconstructions allowed us to identify several independent Asian invasions in Lakes Malawi and Tanganyika, the Congo River, Nigeria and Cameroon. Some invasive M. tuberculata populations display much variation in shell morphology, and overlap in morphospace with M. tuberculata populations native to Africa. Experiments confirmed great ecophenotyic plasticity in some invasive populations, which, in combination with the overlap in disparity with native populations, masks invaders and their dispersal through Africa. Finally, the results of geographic modeling indicate that cryptic M. tuberculata invasions occurred primarily in densely populated areas.

CONCLUSIONS

We reveal the continental nature of invasions of Asian M. tuberculata to Africa. Several of the affected ecosystems have high endemicity in Cerithioidea: Lake Tanganyika has an unparalleled diversity in freshwater cerithioids (>10 endemic genera) and the Congo Basin and Lake Malawi are home to the two largest endemic species clusters of Melanoides in Africa (~12 and ~8 species, respectively). Cerithioids perform ecologically important functions in the benthic ecosystems of African freshwaters, but invaders and ecosystem change pose risks to their native diversity. We draw suggestions for more effective conservation strategies from our integrated approach.

Multilocus Phylogeny of the Afrotropical Freshwater Crab Fauna Reveals Historical Drainage Connectivity and Transoceanic Dispersal Since the Eocene

Phylogenetic reconstruction, divergence time estimations and ancestral range estimation were undertaken for 66% of the Afrotropical freshwater crab fauna (Potamonautidae) based on four partial DNA loci (12S rRNA, 16S rRNA, cytochrome oxidase one [COI], and histone 3). The present study represents the most comprehensive taxonomic sampling of any freshwater crab family globally, and explores the impact of paleodrainage interconnectivity on cladogenesis among freshwater crabs. Phylogenetic analyses of the total evidence data using maximum-likelihood (ML), maximum parsimony (MP), and Bayesian inference (BI) produced a robust statistically well-supported tree topology that reaffirmed the monophyly of the Afrotropical freshwater crab fauna. The estimated divergence times suggest that the Afrotropical Potamonautidae diverged during the Eocene. Cladogenesis within and among several genera occurred predominantly during the Miocene, which was associated with major tectonic and climatic ameliorations throughout the region. Paleodrainage connectivity was observed with specimens from the Nilo-Sudan and East African coast proving to be sister to specimens from the Upper Guinea Forests in West Africa. In addition, we observed strong sister taxon affinity between specimens from East Africa and the Congo basin, including specimens from Lake Tanganyika, while the southern African fauna was retrieved as sister to the Angolan taxa. Within the East African clade we observed two independent transoceanic dispersal events, one to the Seychelles Archipelago and a second to Madagascar, while we observe a single transoceanic dispersal event from West Africa to São Tomé. The ancestral area estimation suggested a West African/East African ancestral range for the family with multiple dispersal events between southern Africa and East Africa, and between East Africa and Central Africa The taxonomic implications of our results are discussed in light of the widespread paraphyly evident among a number of genera.

Complete mitochondrial genome sequence of the Eastern gorilla (Gorilla beringei) and implications for african ape biogeography

The Western and Eastern species of gorillas (Gorilla gorilla and Gorilla beringei) began diverging in the mid-Pleistocene, but in a complex pattern with ongoing gene flow following their initial split. We sequenced the complete mitochondrial genomes of 1 Eastern and 1 Western gorilla to provide the most accurate date for their mitochondrial divergence, and to analyze patterns of nucleotide substitutions. The most recent common ancestor of these genomes existed about 1.9 million years ago, slightly more recent than that of chimpanzee and bonobo. We in turn use this date as a calibration to reanalyze sequences from the Eastern lowland and mountain gorilla subspecies to estimate their mitochondrial divergence at approximately 380000 years ago. These dates help frame a hypothesis whereby populations became isolated nearly 2 million years ago with restricted maternal gene flow, followed by ongoing male migration until the recent past. This process of divergence with prolonged hybridization occurred against the backdrop of the African Pleistocene, characterized by intense fluctuations in temperature and aridity, while at the same time experiencing tectonic uplifting and consequent shifts in the drainage of major river systems. Interestingly, this same pattern of introgression following divergence and discrepancies between mitochondrial and nuclear loci is seen in fossil hominins from Eurasia, suggesting that such processes may be common in hominids and that living gorillas may provide a useful model for understanding isolation and migration in our extinct relatives.