Molecular studies have demonstrated a deep lineage split between the two gorilla species, as well as divisions within these taxa;estimates place this divergence in the mid-Pleistocene, with gene flow continuing until approximately 80,000 years ago. Here, we present analyses of skeletal data indicating the presence of substantial recent gene flow among gorillas at all taxonomic levels: between populations, subspecies, and species. Complementary analyses of DNA sequence variation suggest that low-level migration occurred primarily in a westerly-to-easterly direction. In western gorillas, the locations of hybrid phenotypes map closely to expectations based on population refugia and riverine barrier hypotheses, supporting the presence of significant vicariancedriven structuring and occasional admixture within this taxon. In eastern lowland gorillas, the high frequency of hybrid phenotypes is surprising, suggesting that this region represents a zone of introgression between eastern gorillas and migrants from the west, and underscoring the conservation priority of this critically endangered group. These results highlight the complex nature of evolutionary divergence in this genus, indicate that historical gene flow has played a major role in structuring gorilla diversity, and demonstrate that our understanding of the evolutionary processes responsible for shaping biodiversity can benefit immensely from consideration of morphological and molecular data in conjunction. Historical and environmental characteristics of a region influence both the origin and maintenance of its biodiversity at multiple levels (Avise 2000;Vellend and Geber 2005). For example, the repeated glacial cycles of the Pleistocene had a significant impact on the biogeographic patterns of many organisms (Avise et al. 1998;Hewitt 2004a); habitat fragmentation during glacial maxima led to isolation and, in many cases, local extinction as populations retreated and their ranges became subdivided into geographically distinct glacial refugia. During these protracted periods of isolation, populations differentiated to varying degrees; under contrasting selection regimes and together with the increased effects of genetic drift, character divergence under allopatry emerged in many lineages (Endler 1977;Avise 2000). Following climate amelioration and periods of population range expansion, zones of secondary contact are likely to have formed between previously refugial populations, with admixture often resulting in increased variation, the production of character clines, and the emergence of novel phenotypic traits (Endler 1977;Barton and Hewitt 1985;Arnold 1997).To date, a large number of studies using methods of phylogeographic inference have provided clear evidence across myriad organisms for resumed gene flow among once isolated (Pleistocene) refugial populations (see examples in reviews by Avise et al. 1998;Taberlet et al. 1998;Hewitt 2000Hewitt , 2004b, and highlight research on hybrid zones as central to our understanding of historic subdivisions and their...