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.

Molecular systematics and biogeographic history of the African climbing-mouse complex (Dendromus)

Climbing mice in the genus Dendromus (sensu lato) are widely distributed in Africa, south of the Saharan Desert. The 17 currently recognized species in the genus range from widespread taxa to single-mountain endemics, and there is considerable variation across species with respect to habitats occupied. These habitats range from arid grasslands and savannahs to sub-alpine and alpine vegetation. Using the most comprehensive geographic and genetic survey to date and after reviewing many type specimens, we assess the systematics and biogeography of Dendromus. Given the structure of our molecular phylogenetic hypotheses, in which we recover six major clades, we propose the recognition of three genera within the Dendromus group (sensu lato): in addition to Dendromus (26 lineages), we suggest the retention of Megadendromus (monotypic) and the resurrection of the genus Poemys (six lineages). From our model-based molecular phylogenetic results and morphological comparisons, we suggest that six formerly synonymized taxa should be resurrected, and we highlight 14 previously undescribed lineages. We also constructed time-calibrations on our phylogeny, and performed ancestral area reconstructions using BioGeoBEARS. Based on fossil evidence, Dendromus appears to have had a widespread African distribution dating back to the Late Miocene (8-10 Ma), and our basal ancestral area reconstruction (Ethiopians Highlands + Eastern African Mountains + Zambezian region) supports this. Divergence of the six major clades we recover (Poemys, Megadendromus and four within Dendromus) occurred prior to or at the Miocene-Pliocene boundary 5.3 Ma. Biogeographically, Megadendromus is restricted to the Ethiopian Highlands. The ancestral area for Poemys is reconstructed as the Zambezian region, with species distributions ranging from South Africa to Western Africa. The ancestral area for Dendromus is reconstructed as the Ethiopian Highlands, with the ancestral areas of the four major clades being reconstructed as Ethiopian Highlands, Albertine Rift, South Africa or Western Africa. None of the four Dendromus clades are reciprocally monophyletic with respect to distributional area.

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.

Mechanisms of radiation in a bat group from the genus Pipistrellus inferred by phylogeography, demography and population genetics

Here, we present a study of the Pipistrellus pipistrellus species complex, a highly diversified bat group with a radiation centre in the Mediterranean biodiversity hotspot. The study sample comprised 583 animals from 118 localities representatively covering the bats' range in the western Palearctic. We used fast-evolving markers (the mitochondrial D-loop sequence and 11 nuclear microsatellites) to describe the phylogeography, demography and population structure of this model taxon and address details of its diversification. The overall pattern within this group includes a mosaic of phylogenetically basal, often morphologically distant, relatively small and mostly allopatric demes in the Mediterranean Basin, as well as two sympatric sibling species in the large continental part of the range. The southern populations exhibit constant size, whereas northern populations show a demographic trend of growth associated with range expansion during the Pleistocene climate oscillations. There is evidence of isolation by distance and female philopatry in P. pipistrellus sensu stricto. Although the northern populations are reproductively isolated, we detected introgression events among several Mediterranean lineages. This pattern implies incomplete establishment of reproductive isolating mechanisms in these populations as well as the existence of a past reinforcement stage in the continental siblings. The occurrence of reticulations in the radiation centre among morphologically and ecologically derived relict demes suggests that adaptive unequal gene exchange within hybridizing populations could play a role in speciation and adaptive radiation within this group.

Genetic structure and contrasting selection pattern at two major histocompatibility complex genes in wild house mouse populations

The mammalian major histocompatibility complex (MHC) is a tightly linked cluster of immune genes, and is often thought of as inherited as a unit. This has led to the hope that studying a single MHC gene will reveal patterns of evolution representative of the MHC as a whole. In this study we analyse a 1000-km transect of MHC variation traversing the European house mouse hybrid zone to compare signals of selection and patterns of diversification at two closely linked MHC class II genes, H-2Aa and H-2Eb. We show that although they are 0.01 cM apart (that is, recombination is expected only once in 10 000 meioses), disparate evolutionary patterns were detected. H-2Aa shows higher allelic polymorphism, faster allelic turnover due to higher mutation rates, stronger positive selection at antigen-binding sites and higher population structuring than H-2Eb. H-2Eb alleles are maintained in the gene pool for longer, including over separation of the subspecies, some H-2Eb alleles are positively and others negatively selected and some of the alleles are not expressed. We conclude that studies on MHC genes in wild-living vertebrates can give substantially different results depending on the MHC gene examined and that the level of polymorphism in a related species is a poor criterion for gene choice.

Identification of the first Toll-like receptor gene in passerine birds: TLR4 orthologue in zebra finch (Taeniopygia guttata)

Toll-like receptors (TLRs) are the basic components of the vertebrate pathogen recognition system. Despite uniform general structure, remarkable variability in domain composition can be found in individual TLRs among species. Knowledge of interspecific differences is of particular importance to our understanding of selective pressures on TLRs. Currently, most TLRs are characterized only in a limited number of model species, including domestic chicken as a universal avian model. Here, we describe structure and expression pattern of TLR4 in zebra finch, a widely used passerine model species. The tgTlr4 gene consists of three exons (204, 167 and 3033-3043 bp) that are transcribed into messenger RNA with a relatively long 3'-untranslated region (788 bp). Predicted protein is composed of 842 amino acids (aas) forming extracellular domain with nine leucine-rich repeat (LRR) motives flanked at the carboxy-terminal end by leucine-rich repeat carboxy-terminal domain, transmembrane domain and cytoplasmic toll/interleukin-1 receptor domain. The overall structure is similar to other known TLR4 molecules with 32%-49% aa identity to various mammals and 74% to chicken. Although the position of most of the domains in zebra finch TLR4 resembles their position in chicken, there is one extra LRR at the aa position 207-229 in tgTLR4 and one LRR known in chTLR4 is missing. The gene is highly expressed in the bone marrow and in the spleen, intermediately in the gut and low expression was found in the liver and lungs. For the first time in birds, expression of tgTLR4 in peritoneal macrophages was found to be enhanced by the Escherichia coli lipopolysaccharide treatment.

Comparative phylogeography of two sibling species of forest-dwelling rodent (Praomys rostratus and P. tullbergi) in West Africa: different reactions to past forest fragmentation

Two sibling species of the rodent genus Praomys occur in West African forests: P. tullbergi and P. rostratus. By sampling across their geographical ranges (459 individuals from 77 localities), we test the hypothesis that climatic oscillations during the Quaternary made an impact on the observed pattern of cytochrome b sequence variation. We show that, although these two species have parapatric geographical distributions, their phylogeographical histories are dissimilar, which could be related to their distinct ecological requirements. Since the arid phases of the Pleistocene were characterized by isolated forest patches, and intervening wetter periods by forest expansion, these changes in forest cover may be the common mechanism responsible for the observed phylogeographical patterns in both of these species. For example, in both species, most clades had either allopatric or parapatric geographical distributions; however, genetic diversity was much lower in P. tullbergi than in P. rostratus. The genetic pattern of P. tullbergi fits the refuge hypothesis, indicating that a very small number of populations survived in distinct forest blocks during the arid phases, then expanded again with forest recovery. In contrast, a number of populations of P. rostratus appear to have survived during the dry periods in more fragmented forest habitats, with varying levels of gene flow between these patches depending on climatic conditions and forest extent. In addition, historical variations of the West African hydrographic network could also have contributed to the pattern of genetic differentiation observed in both species.

Multiple parasites mediate balancing selection at two MHC class II genes in the fossorial water vole: insights from multivariate analyses and population genetics

We investigated the factors mediating selection acting on two MHC class II genes (DQA and DRB) in water vole (Arvicola scherman) natural populations in the French Jura Mountains. Population genetics showed significant homogeneity in allelic frequencies at the DQA1 locus as opposed to neutral markers (nine microsatellites), indicating balancing selection acting on this gene. Moreover, almost exhaustive screening for parasites, including gastrointestinal helminths, brain coccidia and antibodies against viruses responsible for zoonoses, was carried out. We applied a co-inertia approach to the genetic and parasitological data sets to avoid statistical problems related to multiple testing. Two alleles, Arte-DRB-11 and Arte-DRB-15, displayed antagonistic associations with the nematode Trichuris arvicolae, revealing the potential parasite-mediated selection acting on DRB locus. Selection mechanisms acting on the two MHC class II genes thus appeared different. Moreover, overdominance as balancing selection mechanism was showed highly unlikely in this system.

Density-related changes in selection pattern for major histocompatibility complex genes in fluctuating populations of voles

Host-pathogen interactions are of particular interest in studies of the interplay between population dynamics and natural selection. The major histocompatibility complex (MHC) genes of demographically fluctuating species are highly suitable markers for such studies, because they are involved in initiating the immune response against pathogens and display a high level of adaptive genetic variation. We investigated whether two MHC class II genes (DQA1, DRB) were subjected to contemporary selection during increases in the density of fossorial water vole (Arvicola terrestris) populations, by comparing the neutral genetic structure of seven populations with that estimated from MHC genes. Tests for heterozygosity excess indicated that DQA1 was subject to intense balancing selection. No such selection operated on neutral markers. This pattern of selection became more marked with increasing abundance. In the low-abundance phase, when populations were geographically isolated, both overall differentiation and isolation-by-distance were more marked for MHC genes than for neutral markers. Model-based simulations identified DQA1 as an outlier (i.e. under selection) in a single population, suggesting the action of local selection in fragmented populations. The differences between MHC and neutral markers gradually disappeared with increasing effective migration between sites. In the high-abundance year, DQA1 displayed significantly lower levels of overall differentiation than the neutral markers. This gene therefore displayed stronger homogenization than observed under drift and migration alone. The observed signs of selection were much weaker for DRB. Spatial and temporal fluctuations in parasite pressure and locus-specific selection are probably the most plausible mechanisms underlying the observed changes in selection pattern during the demographic cycle.

Duplication, balancing selection and trans-species evolution explain the high levels of polymorphism of the DQA MHC class II gene in voles (Arvicolinae)

Major histocompatibility complex (MHC) genes play important role in host-parasite interactions and parasites are crucial factors influencing the population dynamics of hosts. We described the structure and diversity of exon 2 of the MHC class II DQA gene in three species of voles (Arvicolinae) exhibiting regular multi-annual fluctuations of population density and analysed the processes leading to the observed MHC polymorphism. By using cloning-sequencing methodology and capillary electrophoresis-single strand conformation polymorphism, we described seven sequences in the water, eight in the common, and seven in the bank voles coming from an area of 70 km(2) around the Nozeroy canton in the Jura Mountains (Franche Comté, France). All exon 2 sequences translate to give unique amino acid sequences and positive selection was found to act very intensively on antigen binding sites. We documented the presence of recombination at vole DQA region but its importance in generating allelic polymorphism seems to be relatively limited. For the first time within rodents, we documented the duplication of the DQA gene in all three species with both copies being transcriptionally active. Phylogenetic analysis of allelic sequences revealed extensive trans-species polymorphism within the subfamily although no alleles were shared between species in our data set. We discuss possible role of parasites in forming the recent polymorphism pattern of the DQA locus in voles.

Sexual selection for male dominance reduces opportunities for female mate choice in the European bitterling (Rhodeus sericeus)

Sexual selection involves two main mechanisms: intrasexual competition for mates and intersexual mate choice. We experimentally separated intrasexual (male-male interference competition) and intersexual (female choice) components of sexual selection in a freshwater fish, the European bitterling (Rhodeus sericeus). We compared the roles of multiple morphological and behavioural traits in male success in both components of sexual competition, and their relation to male reproductive success, measured as paternity of offspring. Body size was important for both female choice and male-male competition, though females also preferred males that courted more vigorously. However, dominant males often monopolized females regardless of female preference. Subordinate males were not excluded from reproduction and sired some offspring, possibly through sneaked ejaculations. Male dominance and a greater intensity of carotenoid-based red colouration in their iris were the best predictors of male reproductive success. The extent of red iris colouration and parasite load did not have significant effects on female choice, male dominance or male reproductive success. No effect of parasite load on the expression of red eye colouration was detected, though this may have been due to low parasite prevalence in males overall. In conclusion, we showed that even though larger body size was favoured in both intersexual and intrasexual selection, male-male interference competition reduced opportunities for female choice. Females, despite being choosy, had limited control over the paternity of their offspring. Our study highlights the need for reliable measures of male reproductive success in studies of sexual selection.

Retinoid X receptor suppresses transformation by the v-myb oncogene

The v-myb oncogene of avian myeloblastosis virus causes acute monoblastic leukemia in vivo and transforms myelomonocytic cells in culture. Retinoids are potent regulators of proliferation and differentiation in various cell types, and they can initiate differentiation in certain types of leukemic cells. However, the BM2 v-myb-transformed chicken monoblastic cell line is resistant to retinoic acid treatment. We found that overexpression of the retinoid X receptor confers sensitivity of BM2 cells to retinoic acid, resulting in induction of growth arrest and terminal differentiation. In contrast, the frequency of apoptosis was not affected by the retinoid X receptor in this cell type. We also demonstrated that suppression of transformation by v-Myb results from the negative effect of retinoid X receptor on v-Myb transactivation function, similar to that previously described for the retinoic acid receptor. The retinoid X receptor-induced inhibition of transactivation by v-Myb seems to be enhanced by a cell type-specific factor(s), which is not required by retinoic acid receptor.