Genetic diversity of the largest African mole-rat genus, Bathyergus. One, two or four species?

Recent advances in sequencing technology and phylogenetic methods allow us to solve puzzling taxonomic questions using detailed analyses of genetic diversity of populations and gene flow between them. The genus of solitary-living dune mole-rat, Bathyergus, is quite unique among six genera of African mole-rats. The animals are by far the largest and the only scratch digging mole-rat genus possessing a skull less adapted to digging, grooved upper incisors, and more surface locomotor activity. Most authors recognize two species of dune mole-rats, B. suillus and B. janetta, but according to others, the genus is monotypic. In addition, recent molecular studies have revealed cryptic genetic diversity and suggested the existence of up to four species. In our study, we used mitochondrial and genome-wide nuclear data collected throughout the distribution of the genus to investigate the number of species. In agreement with previous studies, we found Bathyergus to be differentiated into several distinct lineages, but we also found evidence for a degree of gene flow between some of them. Furthermore, we confirmed that B. janetta is nested within B. suillus, making the latter paraphyletic and we documented an instance of local mitochondrial introgression between these two nominal species. Phylogeographic structure of the genus was found to be very shallow. Although traditionally dated to the Miocene, we found the first split within the genus to be much younger estimated to 0.82 Ma before present. Genealogical distinctiveness of some lineages was very low, and the coancestry matrix showed extensive sharing of closely related haplotypes throughout the genus. Accordingly, Infomap clustering on the matrix showed all populations to form a single cluster. Overall, our study tends to support the existence of only one species of Bathyergus namely, B. suillus. Environmental niche modelling confirmed its dependence on sandy soils and the preference for soils with relatively high carbon content. Bayesian skyline plots indicate recent population decline in the janetta lineage, probably related to global environmental change.