Mole rats are exclusively subterranean and highly specialized rodents. Their long lifespans, remarkable anti-cancer mechanisms, and various distinctive adaptive features make them a useful research model. Moreover, opposing convergence of morphological traits, they have developed extremely high karyotype variability. Thus, 74 chromosomal forms have been described so far and new ones are being revealed continuously. These evolved during the process of rapid radiation and occur in different biogeographical regions. During research into their reproductive biology we have already provided substantial evidence for species-level separation of these taxa. Here, we review diverse chromosomal forms of the lesser blind mole rat, Mediterranean Nannospalax leucodon, distributed in South-eastern Europe, their karyotype records, biogeography, origin, and phylogeny from our extensive research. In the light of new data from molecular genetic studies, we question some former valuations and propose a cryptospecies rank for seven reproductively isolated chromosomal forms with sympatric and parapatric distribution and clear ecogeographical discrepances in their habitats, as well as new experimental and theoretical methods for understanding the courses of speciation of these unique fossorial mammals.

As a complex skeletal organ consisting of 2 functional and developmental units (ascending ramus and alveolar region), the mandible represents a well-established model in morphological integration studies. The concept of morphological integration assumes that developmentally or functionally related traits are more correlated than others and hence evolve together. We compared the level and pattern of mandibular morphological integration between groups of adult yellow-necked field mice ( Apodemus flavicollis ), with and without B chromosomes (Bs) in a population from Mt. Avala, Serbia. Bs are dispensable supernumerary chromosomes characterized by irregular and non–Mendelian modes of inheritance. The level of morphological integration was higher in animals with Bs. One of the 2 regions of the mandible tested (alveolar region) was significantly more affected by the presence of Bs than the other, with an increase in intensity of integration of 41.61% versus 15.86%. The hypothesis of morphological integration, which postulates disunion of the mandible into 2 distinct functional and developmental modules, was confirmed in animals with Bs. Bs probably have a function because they affect mandible phenotype (although the mechanism is unknown), increase variability within populations, and could lead to selective advantage.

We explored the cryptic speciation of the Nannospalax leucodon species complex, characterised by intense karyotype evolution and reduced phenotypic variability that has produced different lineages, out of which 25 are described as chromosomal forms (CFs), so many cryptic species remain unnoticed. Although some of them should be classified as threatened, they lack the official nomenclature necessary to be involved in conservation strategies. Reproductive isolation between seven CFs has previously been demonstrated. To investigate the amount and dynamics of genetic discrepancy that follows chromosomal changes, infer speciation levels, and obtain phylogenetic patterns, we analysed mitochondrial 16S rRNA and MT-CYTB nucleotide polymorphism among 17 CFs—the highest number studied so far. Phylogenetic trees delineated 11 CFs as separate clades. Evolutionary divergence values overlapped with acknowledged higher taxonomic categories, or sometimes exceeded them. The fact that CFs with higher 2n are evolutionary older corresponds to the fusion hypothesis of Nannospalax karyotype evolution. To participate in conservation strategies, N. leucodon classification should follow the biological species concept, and proposed cryptic species should be formally named, despite a lack of classical morphometric discrepancy. We draw attention towards the syrmiensis and montanosyrmiensis CFs, estimated to be endangered/critically endangered, and emphasise the need for detailed monitoring and population survey for other cryptic species.

Phylogenetic relationships among three species of the genus Apodemus (A. flavicollis,A. sylvaticus and A. agrarius) from seven different localities in Serbia were estimated using arbitrarily primed PCR (AP-PCR) profiling, which allowed clear distinction between the sibling species A. flavicollis and A. sylvaticus. Eight oligonucleotide primers out of the 22 screened, gained a total of 397 amplified bands. The average pairwise genetic distances within Apodemus species ranged from 0.51 between A. flavicollis and A. sylvaticus, to 0.74 between A. sylvaticus and A. agrarius. The large genetic distances obtained among populations of A. flavicollis (0.13) indicate that this part of the Balkan Peninsula was probably a refuge during the last glaciation. Furthermore, differences obtained between populations south and north from the river Danube suggest that this species survived the Last Glacial Maximum in at least two separate refuges. Although our data for molecular markers are from a limited number of localities, extended analyses of C-band distribution on chromosomes from previously studied localities south of the rivers Danube and Sava, indicate that distribution of A. sylvaticus is limited to northern Serbia. Therefore, it is possible that this part of the Balkan Peninsula was not a refuge for this species during the last glaciation.