Karyotype and genetic evolution in speciation of subterranean mole rats of the genus Spalax in Turkey

Karyotype (2fl) and allozymc diversity at 37 gene loci were determined in 69 subterranean mole rats in Turkey belonging to the two superspecies: the ancestor Spalax leucodon ( n = 55; 20 populations) and the descendant S. ehrenbergi ( n = 14; lour populations). We identified remarkable variation of diploid chromosome numbers in the S. leucodon superspecies: 2 n = 38. 40. 50, 54, 60 and 62; and in the S. ehrenbergi superspecies: 2 n = 52, 56 and 58. Genetic diversity indices were low on average in both S. leucodon and S. ehrenbergi superspecies: Allele diversity, A = 1.081 and 1.074; polymorphism, P-5%= 0.077 and 0.068; heterozygosity, H = 0.038 and 0.027; and gene diversity, H, = 0.038 and 0.034, respectively. H ranged from 0 in mesic or semimesic regions to 0.088 in arid Anatolia. We consider the populations with different diploid chromosome numbers, 2 n , as good biological species. Karyotypic diversity may mark extensive ecological speciation. Nei's genetic distances, D (average 0.174, range 0.002 0.422) and ecogeographical criteria suggest that almost each population may represent a different biological species, but critical future testing is necessary to support this claim. Karyotypes and allozymes are nonrandomly distributed across Turkey, displaying remarkable correlations with climatic and biolir factors. Both In and H are significantly correlated with aridity stress (2 n /rainfall, r , =–0.74; P < 0.001), and in our region also with climatic unpredictability. These results support the niche-width genetic variation hypothesis in space and time. Climatic selection in Turkey appears to be a major architect of karyotype and genetic (allozymc) diversity and divergence in mole rat evolution, in both speciation and adaptation.     

Mitochondrial DNA polymorphisms in subterranean mole-rats of the Spalax ehrenbergi superspecies in Israel, and its peripheral isolates

Patterns of mitochondrial DNA (mtDNA) variation were examined in 133 mole-rats constituting all four chromosomal species (2n = 52, 2n = 54, 2n = 58, and 2n = 60) of the Spalax ehrenbergi superspecies in Israel, as well as the peripheral isolates of 2n = 60. In the main range of the complex, a total of 28 mtDNA haplotypes were found in 64 mole-rats, with most haplotypes being unique to either a single chromosomal species or population. mtDNA divergence increased from low to high diploid number in a north-to-south direction in Israel. Overall levels of mtDNA diversity were unexpectedly the highest in the 2n = 60, the youngest species of the complex. The mtDNA haplotypes can be separated into two major groups, 2n = 52-54 and 2n = 58-60, and a phylogenetic analysis for each group revealed evidence of a few haplotypes not sorted by diploid number. The overall patterns of mtDNA divergence seen within and among the four chromosomal species are consistent with the parapatric mode of speciation as suggested from previous studies of allozyme and DNA hybridization. In a separate data set the patterns of mtDNA variation were examined across the main geographic range and across peripheral semi-isolates and isolates of the 2n = 60 chromosomal species. Fifteen haplotypes were found in 69 mole-rats. High levels of mtDNA diversity characterized the main range, semi-isolated, and even some desert isolated populations. The peripheral isolates contain much mtDNA diversity, including novel haplotypes.      

Variations of molar morphology in the Spalax ehrenbergi superspecies: adaptive and phylogenetic significance

This study compares and contrasts variations of enamel fold pattern on the crowns of the molars of subterranean mole rats ( Spalax ) belonging to two superspecies, S. ehrenbergi and S. leucodon , and involving about 20 chromosomal species, distributed parapatrically and ranging on different soil types. The sample studied involved 397 skulls, subdivided as follows: 280 of S. ehrenbergi from Israel, representing the four chromosomal species 2n = 52, 54, 58 and 60; 59 of S. ehrenbergi from Egypt, Lebanon, Syria, Jordan and Turkey; 58 of S. leucodon from Turkey.
In making the comparisons it was necessary to take into account the marked changes of molar pattern that take place as a result of wear.
There is a reasonably good agreement between differences of molar pattern and genetic and immunological distances. The species in Israel and Egypt form a compact group, clearly differentiated from Turkish leucodon; ehrenbergi from Turkey and Syria are intermediate, in accordance with their geographical distribution.
In Israel, differences associated with soil type are at about the same level as differences between species. Local populations show a higher level of diversity, and in a number of cases neighbouring populations differed significantly, sometimes over very short distances. Population differences within species are generally associated with differences of soil and vegetation, indicating the influence of natural selection. There are also differences across species boundaries, implying a phylogenetic (heritage) effect. The nature of the adaptive link between molar pattern and soil is unknown: possible factors are contamination of the food by soil, and toughness of the food which necessitates changes of enamel fold length to improve grinding efficiency.     

Ribosomal DNA (rDNA) spacer polymorphism in mole rats

Genetic differentiation of ribosomal DNA (rDNA) nontranscribed-spacer (NTS) polymorphism was analyzed in 50 individuals from 13 populations among the four chromosomal species (2n = 52, 2n = 54, 2n = 58, and 2n = 60) of subterranean mole rats of the Spalax ehrenbergi complex in Israel. Southern blot analysis with a mouse rDNA probe and two restriction enzymes, EcoRI and BamHI, revealed various sizes of major restriction fragments. The assumption that this variation is due to length polymorphism of NTS DNA was supported by the construction of restriction-site maps. On the basis of the EcoRI and BamHi fragment lengths, we could characterize the major types of NTS rDNA repeating units in each individual. Each member of the central population in the four chromosomal groups of mole rats has a characteristic combination of the NTS types, suggesting that the karyotype groups were genetically diverged. Some near-hybrid-zone populations reflect similarities with the rDNA spectra of a neighbor chromosomal group. This might have resulted from gene flow across the hybrid zones.      

MtDNA polymorphisms: evolutionary significance in adaptation and speciation of subterranean mole rats

Patterns of mtDNA diversity in subterranean mole rats of the Spalax ehrenbergi superspecies (2n = 52, 54, 58 and 60) were previously studied in the main ranges of the four chromosomal species, and specifically in the 2n = 60 species and its peripheral steppe semi-isolates and desert isolates. In the present study we correlated mtDNA diversity indices, nucleon diversity, h , and nucleotide divergence, π, with physical (climatic), biotic (parasites) and biological (genetical, morphological, physiological and behavioural) factors, showing that mtDNA diversity is structured ecogeographically and biologically. The following significant correlations of mtDNA diversity were indicated with: (i) climatic heterogeneity and unpredictability; (ii) levels of ecto- and endoparasites; and (iii) biological diversities, primarily with physiological diversity associated with the energy budget. Small steppe semi-isolates and desert isolates harbour high levels of mtDNA haplotype diversity, some novel, which may be a prerequisite for future speciation events. We conclude that the ecogeographical and biological correlates, as well as the maintenance of mtDNA polymorphisms in small isolated populations, strongly suggest that mtDNA diversity is not neutral. Diversifying natural selection appears to be an important differentiating factor of mtDNA diversity in the twin evolutionary processes of adaptive radiation and active speciation. We suggest critical experiments to substantiate our conclusions and highlight the contribution of mtDNA diversity to fitness, i.e. to the biological function of mtDNA diversity in the evolutionary process.     

The lens protein alpha-B-crystallin of the blind subterranean mole-rat: high homology with sighted mammals

Blind subterranean mole rats, Spalax ehrenbergi, retain a subcutaneous, degenerated eye, which is visually non-functional but which does function in circadian entrainment. Crystallins, members of the small heat shock protein family, constitute approximately 90% of the water-soluble proteins of the transparent eye lens and are crucial for its optical properties, but they are also expressed in other tissues. In our attempt to understand the role of the eye in the blind mole-rat, we now describe the cloning, sequencing, and expression of the cDNA of alpha-B-Crystallin from two species of Spalax (S. galili and S. Judaei, with diploid chromosome numbers 2n=52 and 60, respectively). Spalax alpha- B-Crystallin is highly conserved. It is expressed in many tissues of Spalax, among them Spalax eye. The sequence of the cDNA of alpha-B-Crystallin in the eye and in the heart of Spalax is identical. Further studies are essential to clarify the role of this gene in the lens of an atrophied eye of a visually blind mammal.     

Mhc Class II DNA Polymorphisms within and between Chromosomal Species of the Spalax Ehrenbergi Superspecies in Israel

Restriction fragment length polymorphisms (RFLPs) of two major histocompatibility class II genes (P alpha 1 and Q beta) were studied in 13 populations of four chromosomal species (2n = 52, 54, 58 and 60) of the mole rat, Spalax ehrenbergi superspecies in Israel. A substantial frequency of allelic fragments was found in both genes for all populations, including a desert isolate. In the P alpha 1 gene, one allelic fragment is a result of a deletion mutation which is diagnostic of the 2n = 52 chromosomal species. All other ten allelic variants are the result of point mutations. All mutations are located in a short region flanking the 3' end of the gene. Based on Mhc polymorphisms we confirm earlier evidence that gene flow does not occur between the older chromosomal species (2n = 52, 54, 58), and that reproductive isolation decreases, progressively from the oldest to the youngest species (2n = 60).     

Total body water and adaptive water turnover rate in four chromosomal species of subterranean mole rats of the Spalax ehrenbergi superspecies in Israel

Total body water (TBW) and water turnover rate (WTR) were measured in 24 subterranean mole rats comprising four populations, each belonging to a different chromosomal species of the Spalax ehrenbergi superspecies in Israel. The four species range in different climates: humid-cool (2n = 52); semiarid-cool (2n = 54); humid-warm (2n = 58) and arid-warm (2n = 60). TBW, as a percentage of body weight, measured by tritiated water (HTO), was 72.4%±4–7 in 2n = 52, significantly ( P < 005) higher than the similar estimates 61.7%± 7.2, and 59.4%± 5–3, for 2n = 60 and 58, respectively. A comparison of HTO space, as a percentage of TBW, closely approximated TBW, ranging from 97% to 108%. WTR was high, 218.1 and 230.9 ml/kg^0-75/day in the mesic populations of 2n = 58 and 52, respectively. By contrast, WTR estimates were significantly lower , ( P < 0.001), 150.2 and 148.9 ml/kg^0-75/day in the xeric populations of 2n = 54 and 60, respectively. The biological half-life time, T_1/2, was similar and faster, 32.7 and 27.9 hours in the mesic populations of 2n = 52 and 58, as compared with slower, 47.9 and 40.8 hours in the xeric populations of 2n = 54 and 60, respectively. Urine osmolarity in the most xeric northern Negev steppic population of 2n = 60 (737 ± 45 mmol/kg) was significantly ( P < 0.001) higher than in the other species.
We conclude that adaptive radiation in the Spalax ehrenbergi superspecies involves speciation in semiarid (2n = 54) and arid (2n = 60) climates by physiological adaptations of kidney water conservation, along with multiple morphological, physiological and behavioural syndromes of climatic adaptations to increasing aridity (Nevo, 1986).     

DNA sequence variation in the mitochondrial control region of subterranean mole rats,Spalax ehrenbergi superspecies,in Israel

The complete mitochondrial control region was sequenced for 60 individuals representing different populations for each of the four species of the subterranean mole rat Spalax ehrenbergi superspecies in Israel: Spalax galili (2n = 52), S. golani (2n = 54), S. carmeli (2n = 58), and S. judaei (2n = 60). The control region of all species and populations is very similar both in length (979 to 983 bp) and in base composition. As in agreement with previous surveys on mitochondrial control regions on mammals, the mole rat control region can be divided into a central domain and two flanking domains, ETAS (extended termination associated sequences) and CSB (conserved sequence blocks). Along with the common conserved blocks found in these domains (ETAS1, ETAS2, CSB1, CSB2, and CSB3), we have also detected in all individuals an ETAS1-like and a CSB1-like element, both in the ETAS domain. The most conserved region was the central domain, followed by the CSB and ETAS domains, showing important differences in the four species analyzed. Phylogenetic analysis supported the existence of two clades. One clade contained individuals belonging to Spalax galili (2n = 52) and S. golani (2n = 54), separated in two different branches depending on the species. The other clade contained individuals belonging to S. carmeli (2n = 58) and S. judaei (2n = 60) mixed together, suggesting a more recent event of speciation. Within species we have observed a southward trend of increasing variability. These results have been explained as a consequence of the adaptation of the species to ecological factors such as aridity and temperature stresses.     

Microsatellite diversity in populations of blind subterranean mole rats (Spalax ehrenbergi superspecies) in Israel: speciation and adaptation

The genetic diversity, structure, and relationships of 12 populations of blind subterranean mole rats belonging to four species of the Spalax ehrenbergi superspecies were examined. Diversity levels and divergence of 19 nuclear microsatellite markers were analysed. Results indicate that the 12 populations fall into two major clusters representing the species pairs Spalax galili (2 n  = 52) and S. golani (2 n  = 54) vs. S. carmeli (2 n  = 58) and S. judaei (2 n  = 60). There is very little gene flow between the species pairs, except for one population (Afiq) of S. carmeli (2 n  = 58), which is located near the hybrid zone between S. golani (2 n  = 54) and S. carmeli (2 n  = 58). Higher values of the mean for alleles per locus, gene diversity, and number of unique alleles were observed for the southern species pair, S. carmeli and S. judaei, than for the northern species pair, S. galili and S. golani . Microsatellite diversity increases southward toward the Negev Desert and is positively correlated with aridity stress. Climatic factors of water availability coupled with temperature account for 76.7% of microsatellite variance. Natural selection appears to adaptively determine microsatellite evolution in Spalax regionally across Israel.  © 2004 The Linnean Society of London , Biological Journal of the Linnean Society , 2004, 83 , 229–241.     

Genetic polymorphism of urine deoxyribonuclease I isomerases of subterranean mole rats,Spalax ehrenbergi superspecies,in Israel: Ecogeographical patterns and correlates

Genetic polymorphism of urine deoxyribonuclease I (DNase I) of mole rats was analyzed by isoelectric focusing in a thin-layer polyacrylamide gel (IEF-PAGE). One hundred and three subterranean mole rats, comprising 13 populations belonging to the four chromosomal species (2n = 52, 54, 58, 60) of the actively speciating Spalax ehrenbergi superspecies in Israel, were tested. The following results were indicated. (i) Spalax DNase I consisted of 6-12 major isozymes. (ii) Four phenotypes (numbers in parentheses) were 1 (92), 1-2 (5), 1-3 (4), and 2 (1). The decreasing order of genetic diversity, He, in the four species was 0.37, 0.13, 0.10, and 0.0 for 2n = 58, 52, 54, and 60, respectively. (iii) Spearman rank correlations and multiple regression analyses indicated associations of allele frequencies and genetic diversity with climatic and vegetation factors. We concluded that (a) climatic selection, either directly or indirectly through plant (i.e., food resources) diversity, plays an important role in DNase genetic differentiation and (b) no gene flow and introgression occur between the recent derivative of speciation (2n = 60) and its ancestor (2n = 58), suggesting the operation of reproductive isolation between both species despite natural hybridization.     

Oxygen and carbon dioxide fluctuations in burrows of subterranean blind mole rats indicate tolerance to hypoxic-hypercapnic stresses

The composition of oxygen (O2), carbon dioxide (CO2), and soil humidity in the underground burrows from three species of the Israeli subterranean mole rat Spalax ehrenbergi superspecies were studied in their natural habitat. Two geographically close populations of each species from contrasting soil types were probed. Maximal CO2 levels (6.1%) and minimal O2 levels (7.2%) were recorded in northern Israel in the breeding mounds of S. carmeli in a flooded, poor drained field of heavy clay soil with very high volumetric water content. The patterns of gas fluctuations during the measurement period among the different Spalax species studied were similar. The more significant differentiation in gas levels was not among species, but between neighboring populations inhabiting heavy soils or light soils: O2 was lower and CO2 was higher in the heavy soils (clay and basaltic) compared to the relatively light soils (terra rossa and rendzina). The extreme values of gas concentration, which occurred during the rainy season, seemed to fluctuate with partial flooding of the tunnels, animal digging activity, and over-crowded breeding mounds inhabited by a nursing female and her offspring. The gas composition and soil water content in neighboring sites with different soil types indicated large differences in the levels of hypoxic-hypercapnic stress in different populations of the same species. A growing number of genes associated with hypoxic stress have been shown to exhibit structural and functional differences between the subterranean Spalax and the above-ground rat (Rattus norvegicus), probably reflecting the molecular adaptations that Spalax went through during 40 million years of evolution to survive efficiently in the severe fluctuations in gas composition in the underground habitat.     

POPULATION STRUCTURE OF THE BOTTLENOSE DOLPHIN (TURSIOPS TRUNCATUS) AS DETERMINED BY RESTRICTION ENDONUCLEASE ANALYSIS OF MITOCHONDRIAL DNA

Abstract: Restriction fragment length polymorphisms of mitochondrial DNA (mtDNA) were used to test for population subdivision in the bottlenose dolphin (Tursiops truncatus). Atlantic and Pacific dolphin mtDNA samples exhibited distinctly different haplotypes (approximately 2.4% sequence divergence), indicating a lack of gene exchange. Within the Atlantic Ocean, mtDNA samples from the Gulf of Mexico and the Atlantic Coast were also found to be distinct, with a sequence divergence of approximately 0.6%. The Atlantic Coast–Gulf of Mexico dichotomy is consistent with patterns of genetic variation from other marine and coastal organisms from this region, and supports the hypothesized role of bio-geographic events in promoting the divergence of these and other forms. Regional differentiation was identified along the Atlantic Coast, whereas low sequence divergences among haplotypes and consistent haplotype frequencies across populations suggested considerable gene exchange among Gulf of Mexico populations. A highly divergent haplotype found in two individuals from two localities in the Gulf of Mexico is best explained by dispersal from either a distinct offshore Gulf stock or an unsampled Atlantic Coast stock. Additional samples are required to test for the existence of a distinct offshore race and, if it exists, to identify its distribution and contribution to population structure.     

Mole rat hemoglobin: primary structure and evolutionary aspects in a second karyotype of Spalax ehrenbergi, Rodentia, (2n = 52)

The hemoglobins of the four karyotypes of Spalax ehrenbergi (2n = 52, 54, 58, 60) did not show any differences in their electrophoretic pattern and in high performance liquid chromatography. The complete amino-acid sequence of mole rat hemoglobin (Spalax ehrenbergi), chromosome species 2n = 52, is presented. It was elucidated by automatic Edman degradation of the chains, the tryptic peptides, and the C-terminal peptide obtained by acid hydrolysis of the Asp-Pro bond in beta-chains. The alpha- and beta-chains are identical with those of the chromosome species 2n = 60. A comparison of the hemoglobins of mole rat, mouse, and other rodents shows homology but no indication of adaptation to subterranean life. In all probability alpha 11(A9)Arg and alpha 120(H3)Gly, unique in mole rat among all mammalian hemoglobins, are not involved in high oxygen affinity. The construction of a phylogenetic tree by the maximum parsimony method, based on hemoglobin sequences, made it possible to show that Rodentia originated as a monophyletic clade, and to find the phylogenetic relationship of Spalacidae to other Rodentia (Mus, Rattus, Ondatra, Mesocricetus, Citellus, and Cavia). Among all rodents the slowest rate of nucleotide replacements occurred in the lineage to Spalax (20%) and the fastest in the lineage to Cavia (59%).     

Ribosomal DNA non-transcribed spacer polymorphism in subterranean mole rats: genetic differentiation,environmental correlates and phylogenetic relationships

Summary An analysis is presented of genetic differentiation in the non-transcribed spacers of ribosomal DNA (NTS rDNA). Diversity, environmental correlates and the phylogenetic relationships are examined within and between species of the actively speciating subterranean mole rat, superspeciesSpalax ehrenbergi (2n=52, 54, 58, 60) in Israel. This analysis is based on a previous study of the geographic distribution of restriction fragment length polymorphisms of NTS rDNA. Here we present results indicating that NTS rDNA diversity exists mostly (66%) within populations, while 20% is between populations within species, and 14% between species. Multivariate discriminant analysis succeeded in separating 10 of the 13 populations (77%) into their correct chromosomal species, on the basis of the combination of three NTS rDNA repetypes. The phylogenetic relationships suggest that the complex involves two pairs of closely related species (2n=52–54 and 2n=58–60). NTS rDNA diversity, as well as the decrease southward in frequency of repetype C, are correlated with climatic factors of humidity and temperature. These data are discussed in terms of the evolutionary forces of migration and selection which may cause NTS rDNA differentiation. Climatic selection appears to be the major differentiating factor of NTS rDNA.     

Phylogeography of the Alpine salamander, Salamandra atra (Salamandridae) and the influence of the Pleistocene climatic oscillations on population divergence

Fifty individuals of the endemic Alpine salamander, Salamandra atra, representing 13 populations throughout the range of the two currently recognized subspecies, atra and aurorae, were examined for sequence variation in a large portion (1050 bp) of the mitochondrial cytochrome b gene. We revealed a large number of mitochondrial DNA (mtDNA) haplotypes (10). Interpopulation sequence divergence was very low, ranging from 0 to 3.1%. The relationships among haplotypes were poorly resolved. The divergence time estimate between several mtDNA haplotypes suggested a pre-Pleistocene differentiation approximately 3 million years ago. Moreover, the impact of the Pleistocene glaciations on the phylogeographical patterns appears to have been secondary, although a somewhat reduced genetic variability was found in populations living in areas that were directly affected by the glaciation.     

A switch from diurnal to nocturnal activity in S. ehrenbergi is accompanied by an uncoupling of light input and the circadian clock

The subterranean mole rat Spalax ehrenbergi superspecies represents an extreme example of adaptive visual and neuronal reorganization. Despite its total visual blindness, its daily activity rhythm is entrainable to light-dark cycles, indicating that it can confer light information to the clock. Although most individuals are active during the light phase under laboratory conditions (diurnal animals), some individuals switch their activity period to the night (nocturnal animals). Similar to other rodents, the Spalax circadian clock is driven by a set of clock genes, including the period (sPer) genes. In this work, we show that diurnal mole rats express the Per genes sPer1 and sPer2 with a peak during the light period. Light can synchronize sPer gene expression to an altered light-dark cycle and thereby reset the clock. In contrast, nocturnal Spalax express sPer2 in the dark period and sPer1 in a biphasic manner, with a light-dependent maximum during the day and a second light-independent maximum during the night. Although sPer1 expression remains light inducible, this is not sufficient to reset the molecular clockwork. Hence, the strict coupling of light, Per expression, and the circadian clock is lost. This indicates that Spalax can dissociate the light-driven resetting pathway from the central clock oscillator.     

Penial differentiation in speciation of subterranean mole rats Spalax ehrenbergi in Israel

Geographic and ontogenetic variation was examined in the glans penes and bacula of 92 specimens representing four biological species with the karyotypes 2n=52, 54, 58 and 60 of the Spalax ehrenbergi superspecies in 13 populations from Israel and one disjunct population from northern Egypt. Five measured and three calculated penial characters were analysed and employed in multivariate statistics. Comparison was made between penial morphological characters and craniometric biochemical and chromosomal data. The baculum is strongly species-specific and its length bears no significant relationships with head and body length. The results of these analyses suggest that different bacular morphologies may provide a basis for species-specific prezygotic isolating mechanisms, thus, together with other factors, contributing to reproductive isolation and speciation.     

Preliminary results on the genetic variability of mitochondrial DNA in the signal crayfish,Pacifastacus leniusculus Dana

This study presents a revised method for the extraction of mitochondrial genome (mtDNA) from crayfish species of aquacultural interest, Pacifastacus leniusculus. The mean size of the mitochondrial genome is approximately 18 000 bp. Restriction fragment length polymorphism (RFLP) analysis performed on the mtDNA has revealed extensive genetic variability within this species. The estimated percentage of nucleotide sequence divergence among the four haplotypes found for P. leniusculus ranges from 1.16 to 3.99%. These data suggest that RFLP analysis of mtDNA may provide greater resolution than protein electrophoresis for population identification among signal crayfish populations. This marker offers new avenues for aquaculture and for understanding the genetic structure and evolutionary history of crayfish populations.