Differentiation underground: Range-wide multilocus genetic structure of the silvery mole-rat does not support current taxonomy based on mitochondrial sequences

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.

Functional and morphological divergence in the forelimb musculoskeletal system of scratch-digging subterranean mammals (Rodentia: Bathyergidae)

Whether the forelimb-digging apparatus of tooth-digging subterranean mammals has similar levels of specialization as compared to scratch-diggers is still unknown. We assessed the scapular morphology and forelimb musculature of all four solitary African mole rats (Bathyergidae): two scratch-diggers, Bathyergus suillus and Bathyergus janetta, and two chisel-tooth diggers, Heliophobius argenteocinereus and Georychus capensis. Remarkable differences were detected: Bathyergus have more robust neck, shoulder, and forearm muscles as compared to the other genera. Some muscles in Bathyergus were also fused and often showing wider attachment areas to bones, which correlate well with its more robust and larger scapula, and its wider and medially oriented olecranon. This suggests that shoulder, elbow, and wrist work in synergy in Bathyergus for generating greater out-forces and that the scapula and proximal ulna play fundamental roles as pivots to maximize and accommodate specialized muscles for better (i) glenohumeral and scapular stabilization, (ii) powerful shoulder flexion, (iii) extension of the elbow and (iv) flexion of the manus and digits. Moreover, although all bathyergids showed a similar set of muscles, Heliophobius lacked the m. tensor fasciae antebrachii (aiding with elbow extension and humeral retraction), and Heliophobius and Georychus lacked the m. articularis humeri (aiding with humeral adduction), indicating deeper morphogenetic differences among digging groups and suggesting a relatively less specialized scratch-digging ability. Nevertheless, Heliophobius and Bathyergus shared some similar adaptations allowing scratch-digging. Our results provide new information about the morphological divergence within this family associated with the specialization to distinct functions and digging behaviors, thus contributing to understand the mosaic of adaptations emerging in phylogenetically and ecologically closer subterranean taxa. This and previous anatomical studies on the Bathyergidae will provide researchers with a substantial basis on the form and function of the musculoskeletal system for future kinematic investigations of digging behavior, as well as to define potential indicators of scratch-digging ability.

Utility of cytochrome c oxidase I for the deciphering of unstable phylogeny and taxonomy of gorals, genus Nemorhaedus Hamilton Smith, 1827 (Bovidae, Ovibovina)

Gorals represent ungulate mammals of the Palearctic and Indo-Malayan realms that face habitat destruction and intense hunting pressure. Their classification has been the subject of various (mainly genetic) assessments in the last decade, but some results are conflicting, hampering some conservation-based decisions. Genetic sampling of gorals has increased considerably in recent years, at least for mitochondrial (mt) DNA. Results based on two mt genes (cytochrome b and the D-loop) are currently available. Still, the utility of cytochrome oxidase subunit I remains unanalysed, even though it belongs among the gene markers that enable a correct species identification in mammals. This study examines phylogenetic relationships and species delimitation in gorals using all currently available cytochrome oxidase subunit I sequences, including the not yet analysed goral population from Pakistan. Our results of various phylogenetic approaches, such as maximum parsimony, likelihood and Bayesian inference, and exploration of species boundaries via species delimitation support the validity of six species of goral, namely N.baileyi, N.caudatus, N.cranbrooki, N.evansi, N.goral, and N.griseus. This result accords well with results based on other mt genes, especially the cytochrome b from the highly exhaustive data sampling. Our study also summarises common sources of errors in the assessment of goral phylogeny and taxonomy and highlights future priorities in understanding goral diversification.

Cryptic diversity of Crocidura shrews in the savannahs of Eastern and Southern Africa

Crocidura (Eulipotyphla, Soricidae) is the most species-rich genus among mammals, with high cryptic diversity and complicated taxonomy. The hirta-flavescens group of Crocidura represents the most abundant and widespread shrews in savannahs of eastern and southern Africa, making them a suitable phylogeographical model for assessing the role of paleoclimatic changes on current biodiversity in open African habitats. We present the first comprehensive study on the phylogeography, evolutionary history, geographical distribution, systematics, and taxonomy of the group, using the integration of mitochondrial, genome-wide (ddRAD sequencing), morphological and morphometrical data collected from specimens over most of the known geographic distribution. Our genomic data confirmed the monophyly of this group and its sister relationship with the olivieri group of Crocidura. There is a substantial genetic variation within the hirta-flavescens group, with three highly supported clades showing parapatric distribution and which can be distinguished morphologically: C. hirta, distributed in both the Zambezian and Somali-Masai bioregions, C. flavescens, known from South Africa and south-western Zambia, and C. cf. flavescens, which is known to occur only in central and western Tanzania. Morphometric data revealed relatively minor differences between C. hirta and C. cf. flavescens, but they differ in the colouration of the pelage. Diversification of the hirta-flavescens group has most likely happened during phases of grassland expansion and contraction during Plio-Pleistocene climatic cycles. Eastern African Rift system, rivers, and the distinctiveness of Zambezian and Somali-Masai bioregions seem to have also shaped the pattern of their diversity, which is very similar to sympatric rodent species living in open habitats. Finally, we review the group's taxonomy and propose to revalidate C. bloyeti, currently a synonym of C. hirta, including the specimens treated as C. cf. flavescens.

Biogeographical Importance of the Livingstone Mountains in Southern Tanzania: Comparative Genetic Structure of Small Non-volant Mammals

The Livingstone Mountains (LM; also known as the Kipengere Range) found in south-western Tanzania at the northern end of Lake Nyasa are an important region for understanding the biogeography of Eastern Africa. The two branches of the East African Rift Valley meet here and the mountains might represent stepping stones for colonization and migration between different parts of the Eastern Afromontane Biodiversity Hotspot (especially the link between the Eastern Arc Mountains, EAM, and the Southern Rift Mountains, SRM), as well as an efficient barrier to gene flow for taxa living in drier savannahs in lower elevations. Here we combine new mitochondrial sequence data from 610 recently sampled rodents and shrews with available georeferenced genetic data (3538 specimens) from southern Tanzania, northern Malawi/Zambia and northern Mozambique and compare the spatial genetic structure among different taxa. There is no universal phylogeographic pattern in taxa preferring humid montane habitats. For some of them, the Makambako Gap acts as a barrier between the SRM and the EAM, but other taxa can bridge this gap. Barriers within the EAM (frequently) and within the SRM (sometimes) appear more important. The Rukwa rift between the SRM and the ARM is an important barrier that perhaps can only be crossed by taxa that are not that strictly tied to humid montane environments. For mammals living in lower-elevation savannah-like habitats, the LM can act as a strict barrier to gene flow, and together with the Ufipa Plateau, Lake Nyasa and the EAM create a very similar phylogeographic pattern with three recognizable genetic groups in most savannah-dwellers. The Livingstone Mountains thus appear to be one of the most important biogeographic crossroads in Eastern Africa.

Evoked auditory potentials from African mole-rats and coruros reveal disparity in subterranean rodent hearing

Hearing in subterranean rodents exhibits numerous peculiarities, including low sensitivity and restriction to a narrow range of comparatively low frequencies. Past studies provided two conflicting hypotheses explaining how these derived traits evolved: structural degeneration and adaptive specialization. To further elucidate this issue, we recorded auditory brainstem responses from three species of social subterranean rodents that differ in the degree of specialization to the underground habitat: the naked mole-rat (Heterocephalus glaber) and the Mashona mole-rat (Fukomys darlingi), which represent the ancient lineage of African mole-rats (Bathyergidae), and the coruro (Spalacopus cyanus), a South American rodent (Octodontidae) that adopted a subterranean lifestyle in more recent geological time. Additionally, we measured call amplitudes of social vocalizations to study auditory vocal coupling. We found elevated auditory thresholds and severe hearing range restrictions in the African mole-rats, with hearing in naked mole-rats tending to be more sensitive than in Mashona mole-rats, in which hearing notably deteriorated with increasing age. In contrast, hearing in coruros was similar to that of epigeic rodents, with its range extending into ultrasonic frequencies. However, as in the mole-rats, the coruros' region of best hearing was located at low frequencies close to 1 kHz. We argue that the auditory sensitivity of African mole-rats, although remarkably poor, has been underestimated by recent studies, whereas data on coruros conform to previous results. Considering the available evidence, we propose to be open to both degenerative and adaptive interpretations of hearing physiology in subterranean mammals, as each may provide convincing explanations for specific auditory traits observed.

Species limits and phylogeographic structure in two genera of solitary African mole-rats Georychus and Heliophobius

African mole-rats (Bathyergidae) are an intensively studied family of subterranean rodents including three highly social and three solitary genera. Although their phylogenetic interrelations are clear, genetic diversity and the number of species within each genus is much less certain. Among the solitary genera, Heliophobius and Georychus were for a long time considered as monotypic, but molecular studies demonstrated strong phylogeographic structure within each genus and proposed that they represent complexes of cryptic species. The present study re-evaluates their internal genetic/phylogenetic structure using a combination of methodological approaches. We generated datasets of one mitochondrial and six specifically selected nuclear markers as well as of a large number of double digest restriction site associated (ddRAD) loci and then applied species delimitation analyses based on the multispecies coalescent model or clustering on co-ancestry matrices. The population structure was largely congruent across all analyses, but the methods differed in their resolution scale when determining distinct gene pools. While the multispecies coalescent model distinguished five Georychus and between eleven to thirteen Heliophobius gene pools in both Sanger sequenced and ddRAD loci, two clustering algorithms revealed significantly finer or coarser structure in ddRAD based co-ancestry matrices. Tens of clusters were distinguished by fineRADstructure and one (in Georychus) or two clusters (in Heliophobius) by Infomap. The divergence dating of the bathyergid phylogeny estimated that diversification within both genera coincided with the onset of the Pleistocene and was likely driven by repeated large-scale climatic changes. Based on this updated genetic evidence, we suggest recognizing one species of Georychus and two species of Heliophobius, corresponding to a northern and southern major lineage, separated by the Eastern Arc Mountains. Yet, the final taxonomic revision should await integrated evidence stemming from e.g.. morphological, ecological, or behavioral datasets.

Historical biogeography, systematics, and integrative taxonomy of the non-Ethiopian speckled pelage brush-furred rats (Lophuromys flavopunctatus group)

BACKGROUND

The speckled-pelage brush-furred rats (Lophuromys flavopunctatus group) have been difficult to define given conflicting genetic, morphological, and distributional records that combine to obscure meaningful accounts of its taxonomic diversity and evolution. In this study, we inferred the systematics, phylogeography, and evolutionary history of the L. flavopunctatus group using maximum likelihood and Bayesian phylogenetic inference, divergence times, historical biogeographic reconstruction, and morphometric discriminant tests. We compiled comprehensive datasets of three loci (two mitochondrial [mtDNA] and one nuclear) and two morphometric datasets (linear and geometric) from across the known range of the genus Lophuromys.

RESULTS

The mtDNA phylogeny supported the division of the genus Lophuromys into three primary groups with nearly equidistant pairwise differentiation: one group corresponding to the subgenus Kivumys (Kivumys group) and two groups corresponding to the subgenus Lophuromys (L. sikapusi group and L. flavopunctatus group). The L. flavopunctatus group comprised the speckled-pelage brush-furred Lophuromys endemic to Ethiopia (Ethiopian L. flavopunctatus members [ETHFLAVO]) and the non-Ethiopian ones (non-Ethiopian L. flavopunctatus members [NONETHFLAVO]) in deeply nested relationships. There were distinctly geographically structured mtDNA clades among the NONETHFLAVO, which were incongruous with the nuclear tree where several clades were unresolved. The morphometric datasets did not systematically assign samples to meaningful taxonomic units or agree with the mtDNA clades. The divergence dating and ancestral range reconstructions showed the NONETHFLAVO colonized the current ranges over two independent dispersal events out of Ethiopia in the early Pleistocene.

CONCLUSION

The phylogenetic associations and divergence times of the L. flavopunctatus group support the hypothesis that paleoclimatic impacts and ecosystem refugia during the Pleistocene impacted the evolutionary radiation of these rodents. The overlap in craniodental variation between distinct mtDNA clades among the NONETHFLAVO suggests unraveling underlying ecomorphological drivers is key to reconciling taxonomically informative morphological characters. The genus Lophuromys requires a taxonomic reassessment based on extensive genomic evidence to elucidate the patterns and impacts of genetic isolation at clade contact zones.

Multilocus phylogeny of African striped grass mice (Lemniscomys): Stripe pattern only partly reflects evolutionary relationships

Murine rodents are one of the most evolutionary successful groups of extant mammals. They are also important for human as vectors and reservoirs of zoonoses and agricultural pests. Unfortunately, their fast and relatively recent diversification impedes our understanding of phylogenetic relationships and species limits of many murine taxa, including those with very conspicuous phenotype that has been frequently used for taxonomic purposes. One of such groups are the striped grass mice (genus Lemniscomys), distributed across sub-Saharan Africa in 11 currently recognized species. These are traditionally classified into three morphological groups according to different pelage colouration on the back: (a) L. barbarus group (three species) with several continuous pale longitudinal stripes; (b) L. striatus group (four species) with pale stripes diffused into short lines or dots; and (c) L. griselda group (four species) with a single mid-dorsal black stripe. Here we reconstructed the most comprehensive molecular phylogeny of the genus Lemniscomys to date, using the largest currently available multi-locus genetic dataset of all but two species. The results show four main lineages (=species complexes) with the distribution corresponding to the major biogeographical regions of Africa. Surprisingly, the four phylogenetic lineages are only in partial agreement with the morphological classification, suggesting that the single-stripe and/or multi-striped phenotypes evolved independently in multiple lineages. Divergence dating showed the split of Lemniscomys and Arvicanthis genera at the beginning of Pleistocene; most of subsequent speciation processes within Lemniscomys were affected by Pleistocene climate oscillations, with predominantly allopatric diversification in fragmented savanna biome. We propose taxonomic suggestions and directions for future research of this striking group of African rodents.

Evolutionary history of Pneumocystis fungi in their African rodent hosts

Pneumocystis is a genus of parasitic fungi infecting lung tissues in a wide range of mammal species, displaying a strong host specificity and patterns of co-speciation with their hosts. However, a recent study on Asiatic murids challenged these patterns reporting several Pneumocystis lineages/species shared by different host species or even genera in the Rattini and Murini tribes. Here we screened lung samples of 27 species of African rodents from five families for the presence of Pneumocystis DNA. Using reconstructed multi-locus phylogenies of both hosts and parasites, we tested the hypothesis of their co-evolution. We found that Pneumocystis is widespread in African rodents, detected in all but seven screened host species, with species-level prevalence ranging from 5.9 to 100%. Several host species carry pairs of highly divergent Pneumocystis lineages/species. The retrieved co-phylogenetic signal was highly significant (p = .0017). We found multiple co-speciations, sorting events and two host-shift events, which occurred between Murinae and Deomyinae hosts. Comparison of genetic distances suggests higher substitution rates for Pneumocystis relative to the rodent hosts on neutral loci and slower rates on selected ones. We discuss life-history traits and population dynamics factors which could explain the observed results.

Multiple radiations of spiny mice (Rodentia: Acomys) in dry open habitats of Afro-Arabia: evidence from a multi-locus phylogeny

Background

Spiny mice of the genus Acomys are distributed mainly in dry open habitats in Africa and the Middle East, and they are widely used as model taxa for various biological disciplines (e.g. ecology, physiology and evolutionary biology). Despite their importance, large distribution and abundance in local communities, the phylogeny and the species limits in the genus are poorly resolved, and this is especially true for sub-Saharan taxa. The main aims of this study are (1) to reconstruct phylogenetic relationships of Acomys based on the largest available multilocus dataset (700 genotyped individuals from 282 localities), (2) to identify the main biogeographical divides in the distribution of Acomys diversity in dry open habitats in Afro-Arabia, (3) to reconstruct the historical biogeography of the genus, and finally (4) to estimate the species richness of the genus by application of the phylogenetic species concept.

Results

The multilocus phylogeny based on four genetic markers shows presence of five major groups of Acomys called here subspinosus, spinosissimus, russatus, wilsoni and cahirinus groups. Three of these major groups (spinosissimus, wilsoni and cahirinus) are further sub-structured to phylogenetic lineages with predominantly parapatric distributions. Combination of alternative species delimitation methods suggests the existence of 26 molecular operational taxonomic units (MOTUs), potentially corresponding to separate species. The highest genetic diversity was found in Eastern Africa. The origin of the genus Acomys is dated to late Miocene (ca. 8.7 Ma), when the first split occurred between spiny mice of eastern (Somali-Masai) and south-eastern (Zambezian) savannas. Further diversification, mostly in Plio-Pleistocene, and the current distribution of Acomys were influenced by the interplay of global climatic factors (e.g., Messinian salinity crisis, intensification of Northern Hemisphere glaciation) with local geomorphology (mountain chains, aridity belts, water bodies). Combination of divergence dating, species distribution modelling and historical biogeography analysis suggests repeated “out-of-East-Africa” dispersal events into western Africa, the Mediterranean region and Arabia.

Conclusions

The genus Acomys is very suitable model for historical phylogeographic and biogeographic reconstructions of dry non-forested environments in Afro-Arabia. We provide the most thorough phylogenetic reconstruction of the genus and identify major factors that influenced its evolutionary history since the late Miocene. We also highlight the urgent need of integrative taxonomic revision of east African taxa.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1380-9) contains supplementary material, which is available to authorized users.