Variations in ribosomal DNA and mitochondrial DNA among chromosomal species of subterranean mole rats

Restriction site variations in nuclear ribosomal DNA (rDNA) spacers and mitochondrial DNA (mtDNA) were examined in several populations of mole rats with variable numbers of chromosomes, which represented the two superspecies Spalax leucodon (2n = 38, 54, or 62) and Spalax ehrenbergi (2n = 52, 54, 58, or 60). Sequence divergence of rDNA spacers between the members of the superspecies was approximately 8%, while the variation within each superspecies was 4% on average. The intrasuperspecies differentiation of rDNA spacers was generally associated with changes in the diploid number of chromosomes. However, substantial divergence (approximately 1.5%) was also detected among populations with the same diploid number of chromosomes. The sequence divergence of mtDNA among haplotypes of S. ehrenbergi was 10% or higher and among haplotypes of S. leucodon it was approximately 12%. By contrast, the range of sequence divergence between superspecies was 7.4%-12%. The large divergence of mtDNAs within each superspecies of mole rats may be explained by the preservation of ancient mtDNA polymorphisms.      

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.     

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).     

The effect of protein and salt loading on urinary concentrating ability in four chromosomal species of Spalax ehrenbergi

Urinary concentrating ability was tested under protein and salt load in the four chromosomal species of subterranean mole-rats ( Spalax ehrenbergi superspecies) found in Israel. Protein stress induced by a diet of soyabeans supplemented by agar gel, demonstrated a significant increase in urinary osmolarity (UO). In the species living in the driest and warmest region (2n = 60). UO (1423·101 mOsmol/kg) was significantly (P<0·05) higher than in the other three species (2n = 52. 1172·31 mOsmol/kg; 2n = 54, 1160·116 mOsmol kg; and 2n = 58, 1216·145 mOsmol/kg). Upon salt loading this diet with 0·3 mol NaCl, UO increased significantly. However, when the salt load was increased to 0·45 mol NaCl, UO decreased significantly in all but one species (2n = 60) which maintained UO at 1522·65 mOsmol/kg. A decline in UO was attributed to diuresis resulting from a significant increase in urine and sodium excretion. The kidney, of only the xeric ranging species (2n = 60), demonstrated the ability to produce a hyperosmotic urine, in spite of the high salt load. These results might explain the restricted distribution of S. ehrenbergi: the only species (2n = 60) found in an environment rich in succulents and halophyte plants (steppe). This species appears to push speciation and adaptive radiation to the southern limit of its superspecies.     

Seismic communication in the blind subterranean mole-rat: patterns of head thumping and of their detection in the Spalax ehrenbergi superspecies in Israel

Patterns of seismic communication and somatosensory responses were studied in four chromosomal species (2n= 52, 54, 58, 60) of blind subterranean mole-rats, belonging to the Spalax ehrenbergi superspecies in Israel. The patterns involved head thumping signals and surface-recorded cerebral potentials evoked by vibrational signals. The patterns correspond to the differences in territory size and population density which characterize the four species. It is suggested that this correspondence is related to the adaptive significance of the seismic communication.     

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.     

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.     

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.     

Brain size diversity in adaptation and speciation of subterranean mole rats

We have tested brain size diversity and encephalization in the actively speciating subterranean mole rats of the Spalax ehrenbergi superspecies in Israel. Our sample involved 171 individuals comprising 12 populations and 4 chromosomal species (2n = 52, 54, 58 and 60) distributed parapatrically from the northern Mediterranean region southward (2n = 52, 54→+58→60) into increasingly more arid and unpredictable climatic regimes, approaching the Negev Desert. Our results indicate that relative brain size and encephalization are highest in 2n = 60 as compared with 2n = 52, 54 and 58. We hypothesize that this pattern is adaptive and molded by natural selection. Brain evolution and higher encephalization in the S. ehrenbergi complex appears to be associated with increasing ecological stresses of aridity and climatic unpredictability.     

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.     

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.     

Adaptive differentiation of body size in speciating mole rats

Summary We report the body weight of 1,653 subterranean mole rats comprising 12 populations and 4 chromosomal species (2n=52, 54, 58 and 60) of the Spalax ehrenbergi superspecies in Israel. The sample was collected from 1970 to 1985 and includes all captured animals with a minimal age of 10 months. The results indicated the following. (a) Body weight of males was significantly higher than that of females. (b) There is a southward latitudinal gradient in body size. Northern animals living in cooler and more productive mesic environments are larger than southern animals living in warmer and less productive xeric environments. (c) The interspecific differences for each sex are statistically significant. (d) Body size is negatively correlated with temperature variables, and positively correlated with plant cover (reflecting productivity or food resources) and rainy days. (e) The best predictors of body size, explaining up to 87% of the variation in size included various combinations of temperature variables and plant cover.We conclude that in both adaptation and speciation natural selection is a major agent of differentiation of body size in accordance with multiple factors, primarily temperature and food resources operating on the energetics balance.     

Adaptive variation in structure and function of kidneys of speciating subterranean mole rats

Summary We report on kidney structure and function in subterranean mammals of four chromosomal species (2n=52, 54, 58 and 60) belonging to the Spalax ehrenbergi superspecies, in relation to their speciation and adaptive radiation from mesic (2n=52) to xeric (2n=60) environments in Israel. Structural variables measured involved: (1) Relative Medullary Thickness, (RMT); (2) Relative Kidney Weight. (RKW); and (3) Percentage of Kidney out of Body Weight (PKW). Functional variables involved: (i) Urine Solid Concentration, (USC); and (ii) Urine Osmotic Concentration (UOC). The results for chromosomal species 2n=52, 54, 58 and 60 indicated nonsignificant increase southward for RMT, but displayed significant increase along the same transect for RKW, PKW, and USC. The UOC was significantly lower in mesic 2n=52 as compared to the other three species when experimental animals were fed in the laboratory on regular carrot food. However, protein stress food (soybean) and salt stress of 0.45 mol NaCl, caused significant, three and a half fold increase of UOC in 2n=52, 54 and 58; but four and a half fold increase in 2n=60, significantly higher than in the other three species. We conclude that both structurally and functionally, the kidneys differentiated adaptively during the Pleistocene evolution of S. ehrenbergi in Israel, in accordance with aridity stress and halophyte food resources towards the desert. Nevertheless, Spalax generally shows clear upper limits in kidney structural and functional capacities, preventing it from colonizing the true desert, south of the 100 mm isohyete.     

Hypoxic survival differs between two mole rat species (Spalax ehrenbergi) of humid and arid habitats.

1. Two chromosomal species, 2n = 52 and 2n = 60 of the mole rat superspecies (Spalax ehrenbergi), occupy humid (2n = 52) and arid (2n = 60) habitats in Israel. 2. Gas conductivity of the soil of the 2n = 52 mole rat is lower than that of the 2n = 60 mole rat, and the 2n = 52 mole rat is better adapted to hypoxia. 3. The hypothesis that the 2n = 52 mole rat can survive to a lower pO2 than the 2n = 60 mole rat was tested. 4. Terminal pO2 (Torr) of four females 2n = 52 was lower, 18.0 +/- 2.9 (SD), than the terminal pO2 of five females 2n = 60, 28.2 +/- 5.1 (SD). 5. The hypoxic survival of the 2n = 52 mole rat as compared to that of the 2n = 60 mole rat correlates with other physiological traits: breathing and heart frequencies, blood hemoglobin and tissue gas tensions.     

Sexual selection and natural selection in body size differentiation of subterranean mole rats

We exlored indirectly, the operation of sexual selection in subterranean mole rats of the Spalax ebrenbergi superspecies in Israel comprising four chromosomal species, 2n = 52, 54, 58 and 60. We reanalzed two previously available data sets of 1. body size differentiation (Nevo et al. 1986a) and 2. the intensity of “Total Aggression” in mole rats (Nevo et al. 1986b). We correlated the mean size difference between the two sexes, in each of the 12 populations of the chromosomal species, with the mean level of agression, and with climatic factors, both displaying significant correlations. The results indicated that for 2n = 52, 54, 58 and 60, the population averae difference in body weight between the sexes decreases southward as follows: 37.7g (30.8 % of females body weight), 39.3g(29.0%) 26.3g(22.8%) and 20.3g (19.3%), respectively. We interpret the higher body size diherential ketween the sexes in the north as due to sexual selection.     

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%).     

Molecular Evolution of Cytochrome b of Subterranean Mole Rats, Spalax ehrenbergi Superspecies, in Israel

We describe the molecular evolution of cytochrome b of blind subterranean mole rats. We examined 12 individuals for nucleotide differences in the region of 402 base pairs of mitochondrial cytochrome b. Each individual represents a different population from the entire ecological and speciational range of the four chromosomal species in Israel (2n= 52, 54, 58, and 60) belonging to the Spalax ehrenbergi superspecies. Our results indicate the following. (i) There are seven first-position transitional differences, compared to 34 variable third positions, with no observed second-position substitutions. (ii) A maximum of four amino acids differences occurs across the range. (iii) Within-species diversity increases southward. Only 1 autoapomorphic substitution characterizes either 2n= 52 or 2n= 54, but 6–11 substitutions characterize 2n= 58, and 9–13 substitutions characterize 2n= 60. (iv) Both parsimony and maximum-likelihood trees suggest two monophyletic groups: (a) 2n= 52 and 54, and (b) 2n= 58 and 60, as identified earlier by other protein and DNA markers. (v) Mitochondrial cytochrome b heterogeneity is significantly correlated with climatic factors (rainfall) and biotic factors (body size and allozymes). We hypothesize that two selective regimes direct cytochrome b evolution in the S. ehrenbergi superspecies: (i) purifying selection in the flooded, mesic, hypoxic northern range of 2n= 52 and 54 and (ii) diversifying selection in the climatically spatiotemporal, xeric, and variable southern range of 2n= 58 and 60. Thus, the molecular evolution of mitochondrial cytochrome b in S. ehrenbergi is explicable by opposite selective stresses across the range of S. ehrenbergi in Israel, associated with the ecological adaptive radiation of the complex. Received: 23 October 1998 / Accepted: 2 May 1999     

Adaptive variation of pelage color within and between species of the subterranean mole rat (Spalax ehrenbergi) in Israel

Summary Color notations of dorsal pelage were analyzed in 451 adult subterranean mole rats, comprising 19 populations and 4 chromosomal species (2n=52, 54, 58 and 60) of the Spalax ehrenbergi superspecies in Israel. In addition, the color notations of soil samples from the collection sites were measured. In an attempt to evaluate the degree of correlation between pelage color of mole rats and the local soil color, each color was compared both macro- and microgeographically. The macrogeographic comparisons were among populations of the 4 chromosomal species, and the microgeographic comparisons were among pairs of geographically neighboring subpopulations which occur under similar climatic conditions but differ in soil type and color. The results indicate that 1. no differences in pelage color notation exist between the sexes; 2. pelage of mole rats is generally characterized by an over-all grayish color, but its shades vary regionally and locally in accordance with varying soil color; 3. a positive correlation was found between pelage and local soil colors in the macrogeographic analysis and this correlation was verified in a microgeographic comparison; 4. specimens from populations of the northern chromosomal species (2n=52, 54) mainly inhabit terra rossa and basalt soils, respectively, with reddish brown tones, and they tend to have a more reddish shade, whereas animals from the central species (2n=58), which inhabits mainly dark soils including alluvial clays, terra rossa and brown rendzina and the souther species 2n=60, which occurs mainly in light soil types such as pale rendzina, sandy loams and loess, tend to have a more yellowish shade. The general appearance of specimens from 2n=58 was dark gray while that of 2n=60 specimens was light gray; 5. mole rats living in xeric environments (particularly 2n=60, but also 2n=54) are lighter than those which live in mesic environments (2n=52, 58). The results support the conclusion that pelage color of strictly subterranean mole rats is subject to selection pressures of overground predation of disharmonious types with their background soil color. The possibility that the pelage color variation and patterns also contributes to better thermoregulation while mole rats are above ground is discussed.     

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.