Mitochondrial DNA variation of the common hippopotamus: evidence for a recent population expansion

Distribution range and human-hippopotamus (Hippopotamus amphibious) conflict in the Lake Tana biosphere reserve, Ethiopia

BACKGROUND

Hippopotamus (Hippopotamus amphibious) is a semi-aquatic mammal that is considered the third largest living land mammal still alive on earth. The main objective of this study was to determine the distribution range and human-hippopotamus conflicts in the Lake Tana Biosphere Reserve, Ethiopia. Data were collected from June 2022 to January 2023. Total count methods were used to estimate the hippopotamus population, and Arc GIS was used to determine the distribution range and habitat suitability of hippopotamus. Structured questionnaire and focus group discussions were used to collect data on human-hippopotamus conflict. Data were analysed using One-way ANOVA, independent t-test, and chi-sqaure test.

RESULTS

The mean population estimation of hippopotamus in in the study area showed 252.67 ± 15.9 with a density of 0.072 individuals per km2. The suitable slope analysis of hippopotamus along the the Lake Tana Biosphere Reserve's shoreline revealed that 4.7% of the shoreline was most suitable, 48.4% moderately suitable, and 46.9% not suitable. Variations in the hippopotamus' habitat suitability within the study area have been associated with a number of factors including availability of water, food, elevation, and vegetation cover. Crop raiding and overgrazing were the primary causes for conflict.

CONCLUSION

In order to conserve hippopotamus in the in the Lake Tana Biosphere Reserve, appropriate conservation measures including habitat protection and restoration, and community engagement and education need to be developed. Moreover, further research on the ecological aspects of the Lake's ecosystem is required to ensure the conservation of hippopotamus.

Genetic diversity and demographic history of the largest remaining migratory population of brindled wildebeest (Connochaetes taurinus taurinus) in southern Africa

The blue wildebeest (Connochaetes taurinus) is a keystone species in the savannahs of southern Africa, where it maintains shortgrass plains and serves as an important prey source for large carnivores. Despite being the second-largest migratory wildebeest population, the brindled wildebeest (C. t. taurinus) of the Greater Liuwa Ecosystem (GLE) of western Zambia have remained largely unstudied, until recently. While studies have increased understanding of recent demography, migration, and population limiting factors, the level of genetic diversity, patterns of gene flow, and long-term demographic history of brindled wildebeest in the GLE remains unknown. Most genetic studies of wildebeest have focused on small, heavily-managed populations, rather than large, migratory populations of high conservation significance. We used restriction-site associated DNA sequencing (RAD-seq) to assess genetic diversity, population structure, and demographic history of brindled wildebeest in the GLE. Using SNPs from 1,730 loci genotyped across 75 individuals, we found moderate levels of genetic diversity in GLE brindled wildebeest (He =  0.210), very low levels of inbreeding (FIS =  0.033), and an effective population size of about one tenth the estimated population size. No genetic population structure was evident within the GLE. Analyses of the site frequency spectrum found signatures of expansion during the Middle Pleistocene followed by population decline in the Late Pleistocene and early Holocene, a pattern previously observed in other African ungulates. These results will supplement field studies in developing effective conservation plans for wildebeest as they face continued and increasing threats of habitat loss, poaching, and other human impacts across their remaining range.

Analysis of the mtDNA D-loop Region Casts New Light on Philippine Red Junglefowl Phylogeny and Relationships to Other Junglefowl Species in Asia

Red junglefowl (RJF) is considered the ancestor of domestic chickens. However, the possible maternal origin, genetic diversity, and subspecies classification of the Philippine (PH) RJF remains uncertain. In this study, the complete mitochondrial DNA (mtDNA) D-loop sequence of 55 PH RJFs collected from the mountainous areas of Occidental Mindoro, Palawan, Agusan del Norte, Capiz, Leyte, Iloilo, and Guimaras were analyzed and compared with chicken reference sequences. Phylogenetic analysis revealed multiple maternal origins of the PH RJFs based on haplogroups D, E, and Y classification. This was supported by PH RJFs and RJFs from other Asian countries sharing a clade. A median-joining network also revealed the haplotype sharing of the PH RJFs and Indonesian RJF, demonstrating common maternal ancestry. High haplotype and nucleotide diversity were also observed at all sampling sites. Analysis of molecular variance indicated that the principal molecular variance existed within populations (81.23%) rather than among populations (18.77%). A population neutrality test and Bayesian skyline plot (BSP) analysis elucidated the RJF maternal effective population size expansion in the Philippines that possibly started approximately 2,800-3,000 years ago. The co-existence of Gallus gallus bankiva and Gallus gallus gallus in the Philippines was also verified. The haplotype sharing of the current RJF samples with commercial chickens suggested the need to formulate conservation programs that would protect the RJFs in the Philippines.

Challenges of devising a milk recipe in a hand-reared hippopotamus (Hippopotamus amphibius)

In January 2017, a Nile hippopotamus (Hippopotamus amphibious) was born approximately six weeks premature at the Cincinnati Zoo & Botanical Garden. Due to the calf's weakened condition and lack of interest from the dam, management at the zoo made the decision to hand-rear the calf. Limited published information was available on hand-rearing this species of hippopotamus (hippo). To devise a nutritionally appropriate recipe, milk samples were acquired from the dam on Days 1, 3, 8, and 9 postpartum, and assayed for sugar, protein, fat, mineral, and water content using standard methods validated for multiple species of mammals at the Smithsonian National Zoo's Nutrition Science Laboratory. The sugar content stayed relatively constant (mean = 4.5%; range: 4.3%-4.7%). The fat consistently increased from 0.48% to 4.24% (mean = 2%). Excluding Day 9, the protein content gradually decreased from 9.56% to 6.39% (mean = 8%). The dry matter (DM) ranged from 14.38% to 16.72% (i.e., water content of 85.62%-83.28%), with the sum of the solids (sugar, protein, fat, and ash) averaging 98.5% of measured DM. Fat content was lower than expected but within the range of other artiodactyls. Between Days 1 and 8, the trend of decreasing protein and increasing fat was consistent with a change from colostrum to mature milk. The sharp increase in fat and protein with a decrease in sugar on Day 9 may indicate the beginning of the involution of the mammary gland due to lack of nursing stimulus. Utilizing this information, the Cincinnati Zoo was able to formulate a successful milk replacement recipe that allowed the calf to be raised through weaning to maturity.

The origin of sheep settlement in Western Mediterranean

The arrival of Neolithic culture in North Africa, especially domestic animals has been essentially documented from archaeological records. As the data relative to sheep are scarce, we studied the genetic relationship between Moroccan sheep breeds and Mediterranean ones using the sequencing of 628 bp of the mitochondrial DNA control region in 193 Moroccan individuals, belonging to six breeds, and 652 sequences from other breeds in Europe and Middle East. Through Network analysis and an original phylogenetically derived method, the connection proportions of each Moroccan breed to foreign ones were estimated, highlighting the strong links between Moroccan and Iberian breeds. The first founders of Moroccan sheep population were issued at 79% from Iberia and 21% from a territory between Middle East and Africa. Their calculated expansion times were respectively 7,100 and 8,600 years B.P. This suggests that Neolithization was introduced by a double influence, from Iberia and from another route, maybe Oriental or Sub-Saharan. The consequence of the environmental changes encountered by founders from Iberia was tested using different neutrality tests. There are significant selection signatures at the level of Moroccan and European breeds settled in elevated altitudes, and an erosion of nucleotide diversity in Moroccan breeds living in arid areas.

Genetic diversity and population structure of Terapon jarbua (Forskål, 1775) (Teleostei, Terapontidae) in Malaysian waters

A background study is important for the conservation and stock management of a species. Terapon jarbua is a coastal Indo-Pacific species, sourced for human consumption. This study examined 134 samples from the central west and east coasts of Peninsular (West) Malaysia and East Malaysia. A 1446-bp concatenated dataset of mtDNA COI and Cyt b sequences was used in this study and 83 haplotypes were identified, of which 79 are unique haplotypes and four are shared haplotypes. Populations of T. jarbua in Malaysia are genetically heterogenous as shown by the high level of haplotype diversity ranging from 0.9167-0.9952, low nucleotide diversity ranging from 0.0288-0.3434, and high FST values (within population genetic variation). Population genetic structuring is not distinct as shown by the shared haplotypes between geographic populations and mixtures of haplotypes from different populations within the same genetic cluster. The gene flow patterns and population structuring observed among these regions are likely attributed to geographical distance, past historical events, allopatric speciation, dispersal ability and water currents. For instance, the mixture of haplotypes revealed an extraordinary migration ability of T. jarbua (>1200 km) via ancient river connectivity. The negative overall value of the neutrality test and a non-significant mismatch distribution are consistent with demographic expansion(s) in the past. The median-joining network concurred with the maximum likelihood haplotype tree with three major clades resolved. The scarcity of information on this species is an obstacle for future management and conservation purposes. Hence, this study aims to contribute information on the population structure, genetic diversity, and historical demography of T. jarbua in Malaysia.

Rapid ecological specialization despite constant population sizes

BACKGROUND

The bushbuck, Tragelaphus scriptus, is a widespread and ecologically diverse ungulate species complex within the spiral-horned antelopes. This species was recently found to consist of two genetically divergent but monophyletic lineages, which are paraphyletic at mitochondrial (mt)DNA owing to an ancient interspecific hybridization event. The Scriptus lineage (T. s. scriptus) inhabits the north-western half of the African continent while Sylvaticus (T. s. sylvaticus) is found in the south-eastern half. Here we test hypotheses of historical demography and adaptation in bushbuck using a higher-resolution framework, with four nuclear (MGF, PRKCI, SPTBN, and THY) and three new mitochondrial markers (cytochrome b, 12S rRNA, and 16S rRNA).

METHODS

Genealogies were reconstructed for the mitochondrial and nuclear data sets, with the latter dated using fossil calibration points. We also inferred the demographic history of Scriptus and Sylvaticus using coalescent-based methods. To obtain an overview of the origins and ancestral colonisation routes of ancestral bushbuck sequences across geographic space, we conducted discrete Bayesian phylogeographic and statistical dispersal-vicariance analyses on our nuclear DNA data set.

RESULTS

Both nuclear DNA and mtDNA support previous findings of two genetically divergent Sylvaticus and Scriptus lineages. The three mtDNA loci confirmed 15 of the previously defined haplogroups, including those with convergent phenotypes. However, the nuclear tree showed less phylogenetic resolution at the more derived parts of the genealogy, possibly due to incomplete lineage sorting of the slower evolving nuclear genome. The only exception to this was the montane Menelik's bushbuck (Sylvaticus) of the Ethiopian highlands, which formed a monophyletic group at three of four nuclear DNA loci. We dated the coalescence of the two lineages to a common ancestor ∼2.54 million years ago. Both marker sets revealed similar demographic histories of constant population size over time. We show that the bushbuck likely originated in East Africa, with Scriptus dispersing to colonise suitable habitats west of the African Rift and Sylvaticus radiating from east of the Rift into southern Africa via a series of mainly vicariance events.

DISCUSSION

Despite lower levels of genetic structure at nuclear loci, we confirmed the independent evolution of the Menelik's bushbuck relative to the phenotypically similar montane bushbuck in East Africa, adding further weight to previous suggestions of convergent evolution within the bushbuck complex. Perhaps the most surprising result of our analysis was that both Scriptus and Sylvaticus populations remained relatively constant throughout the Pleistocene, which is remarkable given that this was a period of major climatic and tectonic change in Africa, and responsible for driving the evolution of much of the continent's extant large mammalian diversity.

Genetic diversity and biogeography of T. officinale inferred from multi locus sequence typing approach

Taraxacum officinale (Asteraceae) is widely distributed weedy plant used as a traditional medicinal herb. The population genetics and historical biogeography of this plant have remained relatively unexplored. This study explores phylogeny, population genetics and ancestral reconstructions adopting multi locus sequence typing (MLST) approach. MLST sequences dataset was generated from genomics and chloroplast DNA sequences obtained from 31 T. officinale haplotypes located in 16 different countries. Phylogenetic analysis distributed these haplotypes in well differentiated geographic clades. The study suggested a close relationship between Europe and adjacent Asian countries. Populations of these regions predominantly formed common haplogroups, showed considerable level of gene flow and evidence for recombination events across European and Asian population. Biogeographical inferences obtained by applying statistical dispersal-vicariance analysis (S-DIVA) and Bayesian binary MCMC (BBM) analysis showed that T. officinale was putatively originated in Europe. Molecular clock analysis based on ITS dataset suggested that the divergence between Europe and East Asian populations can be dated to 1.07 Mya with subsequent dispersal and vicariance events. Among different spatial process long distance seed dispersal mediated by wind had potentially assisted the population expansion of T. officinale.

Genetic diversity, phylogeographic structure and effect of selection at the mitochondrial hypervariable region of Nigerian chicken populations

In this study, the maternal genetic diversity, phylogenetic relationship and effect of natural selection on indigenous chickens from Nigeria were assessed. A total of 397-bp fragment of the mitochondrial DNA (mtDNA) D-loop region of 171 indigenous chickens from four populations of Nigeria and four commercial egg line strains (two Anak titan, one Giriraja and one Yaffa) as out-groups were analysed. Thirty-one haplotypes (28 from Nigerian chickens and three from commercial strains) and 34 polymorphic sites were identified. The mean haplotypic and nucleotide diversity were found to be 0.39 ± 0.05 and 0.02 ± 0.02, respectively. Majority of Nigerian chicken haplotypes observed were grouped into haplogroup D which originated from Indian subcontinent, suggesting a single maternal lineage. Genetic variation within and between populations accounted for 97.30 and 2.70% of the total genetic variation, respectively, which is in agreement with a recent and maternal founding effect. High number (4) of negatively selected sites observed based on single likelihood ancestral counting (SLAC) model indicated that the sampled Nigerian chicken populations were undergoing purifying selection. This study concluded that there was relatively high genetic diversity and differentiation, thus, this information will probably paveway for further evaluation studies, preservation and improvement of Nigerian chickens as genetic resources towards ensuring food security.

Genetic models reveal historical patterns of sea lamprey population fluctuations within Lake Champlain

The origin of sea lamprey (Petromyzon marinus) in Lake Champlain has been heavily debated over the past decade. Given the lack of historical documentation, two competing hypotheses have emerged in the literature. First, it has been argued that the relatively recent population size increase and concomitant rise in wounding rates on prey populations are indicative of an invasive population that entered the lake through the Champlain Canal. Second, recent genetic evidence suggests a post-glacial colonization at the end of the Pleistocene, approximately 11,000 years ago. One limitation to resolving the origin of sea lamprey in Lake Champlain is a lack of historical and current measures of population size. In this study, the issue of population size was explicitly addressed using nuclear (nDNA) and mitochondrial DNA (mtDNA) markers to estimate historical demography with genetic models. Haplotype network analysis, mismatch analysis, and summary statistics based on mtDNA noncoding sequences for NCI (479 bp) and NCII (173 bp) all indicate a recent population expansion. Coalescent models based on mtDNA and nDNA identified two potential demographic events: a population decline followed by a very recent population expansion. The decline in effective population size may correlate with land-use and fishing pressure changes post-European settlement, while the recent expansion may be associated with the implementation of the salmonid stocking program in the 1970s. These results are most consistent with the hypothesis that sea lamprey are native to Lake Champlain; however, the credibility intervals around parameter estimates demonstrate that there is uncertainty regarding the magnitude and timing of past demographic events.

Genetic consequences of population expansions and contractions in the common hippopotamus (Hippopotamus amphibius) since the Late Pleistocene

Over the past two decades, an increasing amount of phylogeographic work has substantially improved our understanding of African biogeography, in particular the role played by Pleistocene pluvial-drought cycles on terrestrial vertebrates. However, still little is known on the evolutionary history of semi-aquatic animals, which faced tremendous challenges imposed by unpredictable availability of water resources. In this study, we investigate the Late Pleistocene history of the common hippopotamus (Hippopotamus amphibius), using mitochondrial and nuclear DNA sequence variation and range-wide sampling. We documented a global demographic and spatial expansion approximately 0.1-0.3 Myr ago, most likely associated with an episode of massive drainage overflow. These events presumably enabled a historical continent-wide gene flow among hippopotamus populations, and hence, no clear continental-scale genetic structuring remains. Nevertheless, present-day hippopotamus populations are genetically disconnected, probably as a result of the mid-Holocene aridification and contemporary anthropogenic pressures. This unique pattern contrasts with the biogeographic paradigms established for savannah-adapted ungulate mammals and should be further investigated in other water-associated taxa. Our study has important consequences for the conservation of the hippo, an emblematic but threatened species that requires specific protection to curtail its long-term decline.

Genetic variation and phylogenetic relationships of Indian buffaloes of uttar pradesh

India possesses a total buffalo population of 105 million out of which 26.1% inhabit Uttar Pradesh. The buffalo of Uttar Pradesh are described as nondescript or local buffaloes. Currently, there is no report about the genetic diversity, phylogenetic relationship and matrilineal genetic structure of these buffaloes. To determine the origin and genetic diversity of UP buffaloes, we sequenced and analysed the mitochondrial DNA D-loop sequences in 259 samples from entire Uttar Pradesh. One hundred nine haplotypes were identified in UP buffaloes that were defined by 96 polymorphic sites. We implemented neutrality tests to assess signatures of recent historical demographic events like Tajima’s D test and Fu’s Fs test. The phylogenetic studies revealed that there was no geographic differentiation and UP buffaloes had a single maternal lineage while buffaloes of Eastern UP were distinctive from rest of the UP buffaloes.

Pan-African genetic structure in the African buffalo (Syncerus caffer): investigating intraspecific divergence

The African buffalo (Syncerus caffer) exhibits extreme morphological variability, which has led to controversies about the validity and taxonomic status of the various recognized subspecies. The present study aims to clarify these by inferring the pan-African spatial distribution of genetic diversity, using a comprehensive set of mitochondrial D-loop sequences from across the entire range of the species. All analyses converged on the existence of two distinct lineages, corresponding to a group encompassing West and Central African populations and a group encompassing East and Southern African populations. The former is currently assigned to two to three subspecies (S. c. nanus, S. c. brachyceros, S. c. aequinoctialis) and the latter to a separate subspecies (S. c. caffer). Forty-two per cent of the total amount of genetic diversity is explained by the between-lineage component, with one to seventeen female migrants per generation inferred as consistent with the isolation-with-migration model. The two lineages diverged between 145 000 to 449 000 years ago, with strong indications for a population expansion in both lineages, as revealed by coalescent-based analyses, summary statistics and a star-like topology of the haplotype network for the S. c. caffer lineage. A Bayesian analysis identified the most probable historical migration routes, with the Cape buffalo undertaking successive colonization events from Eastern toward Southern Africa. Furthermore, our analyses indicate that, in the West-Central African lineage, the forest ecophenotype may be a derived form of the savanna ecophenotype and not vice versa, as has previously been proposed. The African buffalo most likely expanded and diverged in the late to middle Pleistocene from an ancestral population located around the current-day Central African Republic, adapting morphologically to colonize new habitats, hence developing the variety of ecophenotypes observed today.

Comparative phylogeography of African savannah ungulates

The savannah biome of sub-Saharan Africa harbours the highest diversity of ungulates (hoofed mammals) on Earth. In this review, we compile population genetic data from 19 codistributed ungulate taxa of the savannah biome and find striking concordance in the phylogeographic structuring of species. Data from across taxa reveal distinct regional lineages, which reflect the survival and divergence of populations in isolated savannah refugia during the climatic oscillations of the Pleistocene. Data from taxa across trophic levels suggest distinct savannah refugia were present in West, East, Southern and South-West Africa. Furthermore, differing Pleistocene evolutionary biogeographic scenarios are proposed for East and Southern Africa, supported by palaeoclimatic data and the fossil record. Environmental instability in East Africa facilitated several spatial and temporal refugia and is reflected in the high inter- and intraspecific diversity of the region. In contrast, phylogeographic data suggest a stable, long-standing savannah refuge in the south.

Species delimitation and phylogeography of Aphonopelma hentzi (Araneae, Mygalomorphae, Theraphosidae): cryptic diversity in North American tarantulas

BACKGROUND

The primary objective of this study is to reconstruct the phylogeny of the hentzi species group and sister species in the North American tarantula genus, Aphonopelma, using a set of mitochondrial DNA markers that include the animal "barcoding gene". An mtDNA genealogy is used to consider questions regarding species boundary delimitation and to evaluate timing of divergence to infer historical biogeographic events that played a role in shaping the present-day diversity and distribution. We aimed to identify potential refugial locations, directionality of range expansion, and test whether A. hentzi post-glacial expansion fit a predicted time frame.

METHODS AND FINDINGS

A Bayesian phylogenetic approach was used to analyze a 2051 base pair (bp) mtDNA data matrix comprising aligned fragments of the gene regions CO1 (1165 bp) and ND1-16S (886 bp). Multiple species delimitation techniques (DNA tree-based methods, a "barcode gap" using percent of pairwise sequence divergence (uncorrected p-distances), and the GMYC method) consistently recognized a number of divergent and genealogically exclusive groups.

CONCLUSIONS

The use of numerous species delimitation methods, in concert, provide an effective approach to dissecting species boundaries in this spider group; as well they seem to provide strong evidence for a number of nominal, previously undiscovered, and cryptic species. Our data also indicate that Pleistocene habitat fragmentation and subsequent range expansion events may have shaped contemporary phylogeographic patterns of Aphonopelma diversity in the southwestern United States, particularly for the A. hentzi species group. These findings indicate that future species delimitation approaches need to be analyzed in context of a number of factors, such as the sampling distribution, loci used, biogeographic history, breadth of morphological variation, ecological factors, and behavioral data, to make truly integrative decisions about what constitutes an evolutionary lineage recognized as a "species".

Population trend and distribution of the Vulnerable common hippopotamus Hippopotamus amphibius in the Mara Region of Kenya

The common hippopotamus Hippopotamus amphibius can significantly influence the dynamics of ecosystems and engender serious conflicts with people but, in Kenya, one of the species strongholds, it has been little studied or monitored. We surveyed the hippopotamus population in the Masai Mara National Reserve and the adjoining pastoral ranches in 2006 using foot counts along 155.3 km of the main rivers. We counted 4,170 hippopotamuses in 171 schools. Comparisons with earlier surveys suggest that this population increased by 169.6% between 1971 and 1980 within the reserve and, although it did not increase within the reserve during 1980–2006, it increased by 359.4% outside the reserve during this period against a background of deteriorating habitat conditions. The overall density in 2006 was 26.9 hippopotamuses km-1 of river, equivalent to a biomass of 26,677 kg km-1 of river. The ratio of calves to 100 adults was 9:100 inside the reserve, 10:100 outside the reserve and 6:100 along tributaries of the Mara River, implying that the population is either increasing or that its spatial distribution is being compressed because of range contraction. The apparent increase in the hippopotamus population contrasts with marked contemporaneous declines in the populations of most other large mammalian herbivore species in the Reserve. We discuss possible reasons underlying the increase in the hippopotamus population.

Population genetic structure of savannah elephants in Kenya: conservation and management implications

We investigated population genetic structure and regional differentiation among African savannah elephants in Kenya using mitochondrial and microsatellite markers. We observed mitochondrial DNA (mtDNA) nucleotide diversity of 1.68% and microsatellite variation in terms of average number of alleles, expected and observed heterozygosities in the total study population of 10.20, 0.75, and 0.69, respectively. Hierarchical analysis of molecular variance of mtDNA variation revealed significant differentiation among the 3 geographical regions studied (F(CT) = 0.264; P < 0.05) and a relatively lower differentiation among populations within regions (F(SC) = 0.218; P < 0.0001). Microsatellite variation significantly differentiated among populations within regions (F(SC) = 0.019; P < 0.0001) but not at the regional levels (F(CT) = 0.000; P > 0.500). We attribute the high differentiation at the mitochondrial genome to the matrilineal social structure of elephant populations, female natal philopatry, and probably ancient vicariance. Lack of significant regional differentiation at the nuclear loci vis-a-vis strong differences at mtDNA loci between regions is likely the effect of subsequent homogenization through male-mediated gene flow. Our results depicting 3 broad regional mtDNA groups and the observed population genetic differentiation as well as connectivity patterns should be incorporated in the planning of future management activities such as translocations.

Hybridization between subspecies of waterbuck (Kobus ellipsiprymnus) in zones of overlap with limited introgression

Two subspecies of waterbuck (Kobus ellipsiprymnus), common (Kobus ellipsiprymnus ellipsiprymnus) and defassa (Kobus ellipsiprymnus defassa), are recognized based on differences in rump pattern, coat colour and geographical distribution. These forms are parapatrically distributed with an area of range overlap in East Africa, where phenotypically intermediate populations occur. Variation in 478 bp of the mitochondrial DNA control region and 14 polymorphic microsatellite loci were used to describe the genetic structure and phylogeographical pattern of the species, and to assess if the intermediate populations are the results of hybridization. In total, 186 individuals from 11 localities were analysed. A high degree of genetic differentiation was found between subspecies, although this was most evident from the microsatellite data. Hybridization was suggested in the phenotypically and geographically intermediate Nairobi NP population in Kenya. A neighbour-joining (NJ) tree based on microsatellite population genetic distances grouped Nairobi between the common and defassa populations, and a Bayesian analysis clearly showed introgression. Individuals sampled in Samburu NP, Kenya, had a common waterbuck phenotype, but introgression was suggested by both markers. Although a high degree of maternal defassa input was indicated from the sequence data, the Samburu population grouped with the common waterbuck in the microsatellite population genetic distance tree, with high support. Analyses of linkage disequilibrium and maximum-likelihood estimates of genetic drift suggested that admixture between subspecies is a recent event. The fact that introgression is limited between subspecies could be caused by chromosomal differences, hindering gene flow between common and defassa waterbuck.