Chromosomal speciation and adaptive radiation of mole rats in Asia Minor correlated with increased ecological stress

Evolutionary history and systematics of European blind mole rats (Rodentia: Spalacidae: Nannospalax): Multilocus phylogeny and species delimitation in a puzzling group

Species delimitation is a powerful approach to assist taxonomic decisions in challenging taxa where species boundaries are hard to establish. European taxa of the blind mole rats (genus Nannospalax) display small morphological differences and complex chromosomal evolution at a shallow evolutionary divergence level. Previous analyses led to the recognition of 25 'forms' in their distribution area. We provide a comprehensive framework to improve knowledge on the evolutionary history and revise the taxonomy of European blind mole rats based on samples from all but three of the 25 forms. We sequenced two nuclear-encoded genetic regions and the whole mitochondrial cytochrome b gene for phylogenetic tree reconstructions using concatenation and coalescence-based species-tree estimations. The phylogenetic analyses confirmed that Aegean N. insularis belongs to N. superspecies xanthodon, and that it represents the second known species of this superspecies in Europe. Mainland taxa reached Europe from Asia Minor in two colonisation events corresponding to two superspecies-level taxa: N. superspecies monticola (taxon established herewith) reached Europe c. 2.1 million years ago (Mya) and was followed by N. superspecies leucodon (re-defined herewith) c. 1.5 Mya. Species delimitation allowed the clarification of the taxonomic contents of the above superspecies. N. superspecies monticola contains three species geographically confined to the western periphery of the distribution of blind mole rats, whereas N. superspecies leucodon is more speciose with six species and several additional subspecies. The observed geographic pattern hints at a robust peripatric speciation process and rapid chromosomal evolution. The present treatment is thus regarded as the minimum taxonomic content of each lineage, which can be further refined based on other sources of information such as karyological traits, crossbreeding experiments, etc. The species delimitation models also allowed the recognition of a hitherto unnamed blind mole rat taxon from Albania, described here as a new subspecies.

Chromosome Asynapsis Is the Main Cause of Male Sterility in the Interspecies Hybrids of East Asian Voles (Alexandromys, Rodentia, Arvicolinae)

Closely related mammalian species often have differences in chromosome number and morphology, but there is still a debate about how these differences relate to reproductive isolation. To study the role of chromosome rearrangements in speciation, we used the gray voles in the Alexandromys genus as a model. These voles have a high level of chromosome polymorphism and substantial karyotypic divergence. We investigated testis histology and meiotic chromosome behavior in the captive-bred colonies of Alexandromys maximowiczii, Alexandromys mujanensis, two chromosome races of Alexandromys evoronensis, and their interracial and interspecies hybrids, to explore the relationship between karyotypic differences and male hybrid sterility. We found that the seminiferous tubules of the males of the parental species and the interracial hybrids, which were simple heterozygotes for one or more chromosome rearrangements, contained germ cells at all stages of spermatogenesis, indicating their potential fertility. Their meiotic cells displayed orderly chromosome synapsis and recombination. In contrast, all interspecies male hybrids, which were complex heterozygotes for a series of chromosome rearrangements, showed signs of complete sterility. Their spermatogenesis was mainly arrested at the zygotene- or pachytene-like stages due to the formation of complex multivalent chains, which caused extended chromosome asynapsis. The asynapsis led to the silencing of unsynapsed chromatin. We suggest that chromosome asynapsis is the main cause of meiotic arrest and male sterility in the interspecies hybrids of East Asian voles.

The bowfin genome illuminates the developmental evolution of ray-finned fishes

The bowfin (Amia calva) is a ray-finned fish that possesses a unique suite of ancestral and derived phenotypes, which are key to understanding vertebrate evolution. The phylogenetic position of bowfin as a representative of neopterygian fishes, its archetypical body plan and its unduplicated and slowly evolving genome make bowfin a central species for the genomic exploration of ray-finned fishes. Here we present a chromosome-level genome assembly for bowfin that enables gene-order analyses, settling long-debated neopterygian phylogenetic relationships. We examine chromatin accessibility and gene expression through bowfin development to investigate the evolution of immune, scale, respiratory and fin skeletal systems and identify hundreds of gene-regulatory loci conserved across vertebrates. These resources connect developmental evolution among bony fishes, further highlighting the bowfin's importance for illuminating vertebrate biology and diversity in the genomic era.

Phylogeny and biogeography of the African Bathyergidae: a review of patterns and processes

Background

We review genealogical relationships, biogeographic patterns and broad historical drivers of speciation within the Bathyergidae, a group of endemic African rodents, as well as identify key taxa which need further research.

Methods

We sourced comparable cytochrome b sequence data (comparable data available for all members for the Family) and geographic information for all six genera of the African subterranean rodent. This information was combined into the most comprehensive and geographically representative evolutionary study for the Bathyergidae to date.

Results

Species richness within the Bathyergidae appears to be underestimated, with undescribed taxa in five of the six genera. Biogeographic patterns suggest large historical distributions, which were repeatedly fragmented by major landscape changes (especially rifting, uplift and drainage evolution) since the Miocene. Aside from vicariant events, other factors (ecological specialization, population-level responses and climatic change) may have been instrumental in driving divergences in the Bathyergidae. As such, adaptive differences may exist among both populations and species across their discrete ranges, driving independent evolutionary trajectories among taxa. In addition, highly fragmented distributions of divergent (and often relict) lineages indicates the possibility of narrow endemics restricted to diminishing suitable habitats. From this, it is clear that a systematic revision of the Bathyergidae is necessary; such a revision should include comprehensive sampling of all putative taxa, the addition of genomic information to assess adaptive differences, as well as ecological information.

Habitat and Burrow System Characteristics of the Blind Mole Rat Spalax galili in an Area of Supposed Sympatric Speciation

A costly search for food in subterranean rodents resulted in various adaptations improving their foraging success under given ecological conditions. In Spalax ehrenbergi superspecies, adaptations to local ecological conditions can promote speciation, which was recently supposed to occur even in sympatry at sites where two soil types of contrasting characteristics abut each other. Quantitative description of ecological conditions in such a site has been, nevertheless, missing. We measured characteristics of food supply and soil within 16 home ranges of blind mole rats Spalax galili in an area subdivided into two parts formed by basaltic soil and pale rendzina. We also mapped nine complete mole rat burrow systems to compare burrowing patterns between the soil types. Basaltic soil had a higher food supply and was harder than rendzina even under higher moisture content and lower bulk density. Population density of mole rats was five-times lower in rendzina, possibly due to the lower food supply and higher cover of Sarcopoterium shrubs which seem to be avoided by mole rats. A combination of food supply and soil parameters probably influences burrowing patterns resulting in shorter and more complex burrow systems in basaltic soil.

Adaptation of pelage color and pigment variations in Israeli subterranean blind mole rats, Spalax ehrenbergi [corrected]

Background

Concealing coloration in rodents is well established. However, only a few studies examined how soil color, pelage color, hair-melanin content, and genetics (i.e., the causal chain) synergize to configure it. This study investigates the causal chain of dorsal coloration in Israeli subterranean blind mole rats, Spalax ehrenbergi.

Methods

We examined pelage coloration of 128 adult animals from 11 populations belonging to four species of Spalax ehrenbergi superspecies (Spalax galili, Spalax golani, Spalax carmeli, and Spalax judaei) and the corresponding coloration of soil samples from the collection sites using a digital colorimeter. Additionally, we quantified hair-melanin contents of 67 animals using HPLC and sequenced the MC1R gene in 68 individuals from all four mole rat species.

Results

Due to high variability of soil colors, the correlation between soil and pelage color coordinates was weak and significant only between soil hue and pelage lightness. Multiple stepwise forward regression revealed that soil lightness was significantly associated with all pelage color variables. Pelage color lightness among the four species increased with the higher southward aridity in accordance to Gloger's rule (darker in humid habitats and lighter in arid habitats). Darker and lighter pelage colors are associated with darker basalt and terra rossa, and lighter rendzina soils, respectively. Despite soil lightness varying significantly, pelage lightness and eumelanin converged among populations living in similar soil types. Partial sequencing of the MC1R gene identified three allelic variants, two of which were predominant in northern species (S. galili and S. golani), and the third was exclusive to southern species (S. carmeli and S. judaei), which might have caused the differences found in pheomelanin/eumelanin ratio.

Conclusion/Significance

Darker dorsal pelage in darker basalt and terra rossa soils in the north and lighter pelage in rendzina and loess soils in the south reflect the combined results of crypsis and thermoregulatory function following Gloger's rule.

Possible incipient sympatric ecological speciation in blind mole rats (Spalax)

Sympatric speciation has been controversial since it was first proposed as a mode of speciation. Subterranean blind mole rats (Spalacidae) are considered to speciate allopatrically or peripatrically. Here, we report a possible incipient sympatric adaptive ecological speciation in Spalax galili (2n = 52). The study microsite (0.04 km(2)) is sharply subdivided geologically, edaphically, and ecologically into abutting barrier-free ecologies divergent in rock, soil, and vegetation types. The Pleistocene Alma basalt abuts the Cretaceous Senonian Kerem Ben Zimra chalk. Only 28% of 112 plant species were shared between the soils. We examined mitochondrial DNA in the control region and ATP6 in 28 mole rats from basalt and in 14 from chalk habitats. We also sequenced the complete mtDNA (16,423 bp) of four animals, two from each soil type. Remarkably, the frequency of all major haplotype clusters (HC) was highly soil-biased. HCI and HCII are chalk biased. HC-III was abundant in basalt (36%) but absent in chalk; HC-IV was prevalent in basalt (46.5%) but was low (20%) in chalk. Up to 40% of the mtDNA diversity was edaphically dependent, suggesting constrained gene flow. We identified a homologous recombinant mtDNA in the basalt/chalk studied area. Phenotypically significant divergences differentiate the two populations, inhabiting different soils, in adaptive oxygen consumption and in the amount of outside-nest activity. This identification of a possible incipient sympatric adaptive ecological speciation caused by natural selection indirectly refutes the allopatric alternative. Sympatric ecological speciation may be more prevalent in nature because of abundant and sharply abutting divergent ecologies.

Stress, adaptation, and speciation in the evolution of the blind mole rat, Spalax, in Israel

Environmental stress played a major role in the evolution of the blind mole rat superspecies Spalax ehrenbergi, affecting its adaptive evolution and ecological speciation underground. Spalax is safeguarded all of its life underground from aboveground climatic fluctuations and predators. However, it encounters multiple stresses in its underground burrows including darkness, energetics, hypoxia, hypercapnia, food scarcity, and pathogenicity. Consequently, it evolved adaptive genomic, proteomic, and phenomic complexes to cope with those stresses. Here I describe some of these adaptive complexes, and their theoretical and applied perspectives. Spalax mosaic molecular and organismal evolution involves reductions or regressions coupled with expansions or progressions caused by evolutionary tinkering and natural genetic engineering. Speciation of Spalax in Israel occurred in the Pleistocene, during the last 2.00-2.35 Mya, generating four species associated intimately with four climatic regimes with increasing aridity stress southwards and eastwards representing an ecological speciational adaptive trend: (Spalax golani, 2n=54→S. galili, 2n=52→S. carmeli, 2n=58→S. judaei, 2n=60). Darwinian ecological speciation occurred gradually with relatively little genetic change by Robertsonian chromosomal and genic mutations. Spalax genome sequencing has just been completed. It involves multiple adaptive complexes to life underground and is an evolutionary model to a few hundred underground mammals. It involves great promise in the future for medicine, space flight, and deep-sea diving.

Karyotypic and molecular polymorphisms in Ctenomys torquatus (Rodentia: Ctenomyidae): taxonomic considerations

The rodent genus Ctenomys (tuco-tucos) comprises more than 60 described species, and shows extraordinary inter- and intraspecific karyotypic variation. The most widely distributed species of Ctenomys in Brazil is C. torquatus. Although several cytogenetic studies have been done, the karyotypic variability of this species is still poorly known. In this paper we report two new diploid numbers for C. torquatus: 2n = 40 and 2n = 42, both showing AN = 72. The new distribution limits of C. torquatus here reported include localities in the southern, central and western parts of Rio Grande do Sul (RS) State in southern Brazil. The phylogenetic relationship between C. torquatus from Alegrete, RS, and Ctenomys sp. from Corrientes, Argentina, is described by means of mtDNA cytochrome b analysis. Although both entities share similar karyotypes and sperm morphology, these two species are not phylogenetically close.

Chromosomal phylogeny and evolution of the African mole-rats (Bathyergidae)

The subterranean African mole-rats (Family Bathyergidae) show considerable variation in their diploid numbers, but there is limited understanding of the events that shaped the extant karyotypes. Here we investigate chromosomal evolution in specimens representative of six genera and an outgroup species, the cane rat Thryonomys swinderianus, using flow-sorted painting probes isolated from the naked mole-rat, Heterocephalus glaber (2n = 60). A chromosomal phylogeny based on the cladistic analysis of adjacent syntenies detected by cross-species chromosome painting was not consistent with that obtained using DNA sequences due, in large part, to the conserved nature of the Bathyergus, Georychus and Cryptomys karyotypes. In marked contrast, the Fukomys and Heliophobius species showed extensive chromosome reshuffling, permitting their analysis by a computational approach that has conventionally been employed in comparative genomic studies for retrieving phylogenetic information based on DNA sequence or gene order data. Using the multiple genome rearrangements (MGR) algorithm and chromosomal rearrangement data detected among F. damarensis, F. darlingi, F. mechowi and the sister taxa B. janetta and Heliophobius argenteocinereus, cytogenetic support for the monophyly of Fukomys and a sister association for F. darlingi + F. damarensis was retrieved, mirroring the published sequence-based topology. We show that F. damarensis, a lineage adapted to arid and climatically unpredictable environments in Southern Africa, is characterized by a large number of fissions the fixation of which has probably been favoured by environmental factors and/or its particular eusocial structure.

Autosomal rearrangements in Graomys griseoflavus (Rodentia): a model of non-random Robertsonian divergence

The south American rodent Graomys griseoflavus exhibits a remarkable chromosome polymorphism as a consequence of four Robertsonian fusions. Focusing on the genetic analysis of the taxon, genome organization of all karyomorphs was studied at chromosome and molecular organization level. Cytogenetic (G, NOR and Re banding) and molecular (satellite and mitochondrial DNAs) events accompanying chromosome divergence allowed tracing a phylogenetic relationship among all karyomorphs. Available data led to propose that chromosome evolution of G. griseoflavus occurred in a non-random sequence of centric fusions, supporting the hypothesis of single origin for Robertsonian karyomorphs.

Pericentromeric organization at the fusion point of mouse Robertsonian translocation chromosomes

In mammals, Robertsonian (Rb) translocation (the joining of two telo/acrocentric chromosomes at their centromere to form a metacentric) is the most effective process in chromosomal evolution leading to speciation; its occurrence also affects human health (through the induction of trisomies) and the fertility of farm animals. To understand the mechanism of Rb translocation, we used the house mouse as a model system and studied the organization of pericentromeric satellite DNAs (satDNA) of telocentrics and Rb chromosomes, both minor and major satDNA. The chromosome-orientation fluorescence in situ hybridization (CO-FISH) technique was used to analyze the major satDNA. To detect the very small amount of minor satDNA, a procedure was developed that combines CO-FISH with primed in situ labeling and conventional FISH and is five times more sensitive than the CO-FISH procedure alone. It was found that both the major and the minor satDNA tandem repeats are oriented head-to-tail in telocentric and Rb chromosomes, and their polarity is always the same relative to the centromere. We suggest that all tandemly repetitive satDNAs in a species probably are locked into such a symmetry constraint as a universal consequence of chromosomal evolution. Rb translocation breakpoints were found localized within the minor satDNA of telocentrics, and these sequences contributed symmetrically to the formation of the centromeric region of the Rb chromosomes. These results are important for an understanding of the geometry of Rb translocations and suggest the study of DNA orientation as a new tool for investigating these rearrangements.

A new derived and highly polymorphic chromosomal race of Liolaemus monticola (Iguanidae) from the 'Norte Chico' of Chile

A multiple Robertsonian fission chromosomal race of the Liolaemus monticola complex in Chile is described and is shown to be the most derived and the most complex among the Liolaemus examined thus far. The 29 karyotyped lizards analysed from the locality of Mina Hierro Viejo, Petorca, Provincia de Valparaiso, Chile, exhibited a diploid chromosomal number ranging from 42 to 44, and several polymorphisms. The polymorphisms included: a pair 1 fission; a pair 2 fission plus a pericentric inversion in one of the fission products, which moved the NOR and satellite from the tip of the long arm of the metacentric 2 to the short arm of the fission product; a fission in pair 3; a polymorphism for an enlarged chromosome pair 6; and a polymorphism for a pericentric inversion in pair 7. This population is fixed for a fission of chromosome pair 4. A total of 76% of the lizards analysed were polymorphic for one or more pairs of chromosomes. We have compared these data with other Liolaemus monticola chromosomal races and calculated the Hardy-Weinberg ratios for the polymorphic chromosome pairs in this Multiple-Fission race. Karyotypic differences between the Northern (2n = 38-40) and the Multiple-Fission (2n = 42-44) races were attributed mainly to Robertsonian fissions, an enlarged chromosome and pericentric inversions involving the macrochromosomes and one microchromosome pair.

Selection versus random drift: long-term polymorphism persistence in small populations (evidence and modelling)

Our data on a subterranean mammal, Spalax ehrenbergi , and other evidence, indicate that appreciable polymorphism can be preserved in small isolated populations consisting of several dozens of, or a hundred, individuals. Current theoretical models predict fast gene fixation in small panmictic populations without selection, mutation, or gene inflow. Using simple multilocus models, we demonstrate here that moderate stabilizing selection (with stable or fluctuating optimum) for traits controlled by additive genes could oppose random fixation in such isolates during thousands of generations. We also show that in selection-free models polymorphism persists only for a few hundred generations even under high mutation rates. Our multi-chromosome models challenge the hitchhiking hypothesis of polymorphism maintenance for many neutral loci due to close linkage with few selected loci.