Hearing in blind subterranean Zambian mole-rats (Cryptomys sp.): collective behavioural audiogram in a highly social rodent

Thresholds for pure tone detection were examined in the common mole-rat, Cryptomys sp. (Bathyergidae, Rodentia) using a positive reinforcement procedure. To bypass the problems connected with testing isolated individuals of this extremely social species, a collective behavioural audiogram was determined for a family group of seven mole-rats. Within the tested frequency range of 225 to 18 kHz, the lowest thresholds (as low as 7.5 dB SPL, on average 24 dB SPL) were found at 800 Hz, the upper limit of hearing (at the level of 60 dB SPL) was at 18 kHz. The behavioural audiogram combines the results of previous studies on hearing in this species. It resembles the distortion threshold curve but differs from neurophysiological data as far as the high frequency cutoff is concerned. On the other hand, the region of the best hearing sensitivity is narrow in behavioural audiogram and neurophysiological curves but rather broad in the distortion threshold curve. In general, the behavioural audiogram of Cryptomys is in many aspects comparable with the available audiograms of other subterranean rodents. Accepted: 18 February 1997     

Heterochronic quantitative microevolution: Dental divergence in aboriginal Americans

Studies of tooth crown morphology alone have proven valuable in defining human population differentiation. We test the hypothesis that quantitative comparisons of more complex whole tooth structure may prove informative in understanding human diversity. Three disparate populations in Native American genetic history were compared: Kodiak Island Western Eskimos, Peruvian Inca Amerindians, and Southeast Asians. Enamel depth (an increasing gradient extended from Southeast Asians to the Inca) and root parameters were the most discriminating variables. The observed microevolution appears to result from variation in timing of known X-linked, Y-linked, and autosomal genes that affect either ameloblast or odontoblast differentiation. The dental traits were sexually dimorphic, the effect being more pronounced in aboriginal Americans, with male teeth having robust roots and thin enamel compared to female. Southeast Asians were isometrically related. The prominence of sexual dimorphism and the importance of sex-linked genes in the determination of the dental phenotypes suggest that sexual selection was one evolutionary force acting on early Asian populations. Subsequently, the selection appears to have been relaxed in Southeast Asians. Observed divergence of tooth shape among the populations, i.e., differences in the appropriation process of tooth primordia, was mainly the consequence of genetic drift modulating heterochronic regulators of homeotic genes. © 1996 Wiley-Liss, Inc.     

Latitudinal divergence in a widespread amphibian: Contrasting patterns of neutral and adaptive genomic variation

Stochastic effects from demographic processes and selection are expected to shape the distribution of genetic variation in spatially heterogeneous environments. As the amount of genetic variation is central for long‐term persistence of populations, understanding how these processes affect variation over large‐scale geographical gradients is pivotal. We investigated the distribution of neutral and putatively adaptive genetic variation, and reconstructed demographic history in the moor frog (Rana arvalis) using 136 individuals from 15 populations along a 1,700‐km latitudinal gradient from northern Germany to northern Sweden. Using double digest restriction‐site associated DNA sequencing we obtained 27,590 single nucleotide polymorphisms (SNPs), and identified differentiation outliers and SNPs associated with growing season length. The populations grouped into a southern and a northern cluster, representing two phylogeographical lineages from different post‐glacial colonization routes. Hybrid index estimation and demographic model selection showed strong support for a southern and northern lineage and evidence of gene flow between regions located on each side of a contact zone. However, patterns of past gene flow over the contact zone differed between neutral and putatively adaptive SNPs. While neutral nucleotide diversity was higher along the southern than the northern part of the gradient, nucleotide diversity in differentiation outliers showed the opposite pattern, suggesting differences in the relative strength of selection and drift along the gradient. Variation associated with growing season length decreased with latitude along the southern part of the gradient, but not along the northern part where variation was lower, suggesting stronger climate‐mediated selection in the north. Outlier SNPs included loci involved in immunity and developmental processes.     

RAPD divergence caused by microsite edaphic selection in wild barley

Random amplified polymorphic DNA polymerase chain reaction (RAPDPCR) was used to assess genetic diversity in four subpopulations (86 individuals) of wild barley, Hordeum spontaneum, sampled from Tabigha microsite near the Sea of Galilee, Israel. The microsite consists of two 100 m transects that are topographically separated by 100 m, each equally subdivided into 50 m of basalt and terra rossa soil types. Despite the same macroclimate characterizing the area around the Sea of Galilee, the microsite offers two edaphically different microhabitats, with basalt being a more ecologically heterogeneous and broader-niche than the relatively drier but more homogeneous and narrow-niche terra rossa. Analysis of 118 putative loci revealed significant (P<0.05) genetic differentiation in polymorphism (P0.05) between the two soils across the transects with P being higher in the more heterogeneous basalt (mean P0.05 = 0.902), than in terra rossa (mean P0.05 = 0.820). Gene diversity (He) was higher in basalt (mean He=0.371), than in terra rossa (mean He=0.259). Furthermore, unique alleles were confined to one soil type, either in one or both transects. Rare alleles were observed more frequently in terra rossa than basalt, and in transect II only. Gametic phase disequilibria showed a larger multilocus association of alleles in basalt than terra rossa, and in transect I than II. Spearman rank correlation (r_s) revealed a strong association between specific loci and soil types, and transects. Also, analysis of multilocus organization revealed soil-specific multilocus-genotypes. Therefore, our results suggest an edaphically differentiated genetic structure, which corroborates the niche width-variation hypothesis, and can be explained, in part, by natural selection. This pattern of RAPD diversity is in agreement with allozyme and hordein protein diversities in the same subpopulations studied previously. This revised version was published online in July 2006 with corrections to the Cover Date.     

Body temperature patterns and rhythmicity in free-ranging subterranean Damaraland mole-rats, Fukomys damarensis

Body temperature (T(b)) is an important physiological component that affects endotherms from the cellular to whole organism level, but measurements of T(b) in the field have been noticeably skewed towards heterothermic species and seasonal comparisons are largely lacking. Thus, we investigated patterns of T(b) patterns in a homeothermic, free-ranging small mammal, the Damaraland mole-rat (Fukomys damarensis) during both the summer and winter. Variation in T(b) was significantly greater during winter than summer, and greater among males than females. Interestingly, body mass had only a small effect on variation in T(b) and there was no consistent pattern relating ambient temperature to variation in T(b). Generally speaking, it appears that variation in T(b) patterns varies between seasons in much the same way as in heterothermic species, just to a lesser degree. Both cosinor analysis and Fast Fourier Transform analysis revealed substantial individual variation in T(b) rhythms, even within a single colony. Some individuals had no T(b) rhythms, while others appeared to exhibit multiple rhythms. These data corroborate previous laboratory work showing multiplicity of rhythms in mole-rats and suggest the variation seen in the laboratory is a true indicator of the variation seen in the wild.     

Parasites in social subterranean Zambian mole-rats (Cryptomys spp., Bathyergidae, Rodentia)

Three out of 18 examined common mole-rats ( Cryptomys sp., Xkaryotype 2n = 58) from Zambia were infested with cestodes ( Inermicapsifer nmdagascariensis and an undetermined species) and a nematode ( Protospirura muricola ). Four out of 14 examined Zambian giant mole-rats ( Cryptomys mechowi ) hosted cestodes ( Raillietina (R) sp. and an undetermined species) and nematodes ( Protospirura muricola, Capillaria sp.). Helminths and circumstances of findings are briefly described. No ectoparasites were found. The influences of subterranean and social way of life and feeding habits upon occurrence of parasites in mole-rats are discussed. Burrow conditions and social behaviour seem not to favour infestation by parasites.     

Acoustic communication underground: vocalization characteristics in subterranean social mole-rats (Cryptomys sp., Bathyergidae)

In captive adult Zambian mole-rats 14 different sounds (13 true vocalizations) have been recorded during different behavioural contexts. The sound analysis revealed that all sounds occurred in a low and middle frequency range with main energy below 10 kHz. The majority of calls contained components of 1.6–2 kHz, 0.63–0.8 kHz, and/or 5–6.3 kHz. The vocalization range thus matched well the hearing range as established in other studies. The frequency content of courtship calls in two species of Zambian Cryptomys was compared with that in naked mole-rats (Heterocephalus glaber) and blind mole-rats (Spalax ehrenbergi) as described in the literature. The frequency range of maximum sound energy is negatively correlated with the body weight and coincides with the frequencies of best hearing in the respective species. In general, the vocalization range in subterranean mammals is shifted towards low frequencies which are best propagated in underground burrows. Accepted: 16 September 1996     

Phylogeographical patterns of genetic divergence and speciation in African mole-rats (Family: Bathyergidae)

African mole-rats are subterranean Hystricomorph rodents, distributed widely throughout sub-Saharan Africa, and displaying a range of social and reproductive strategies from solitary dwelling to the 'insect-like' sociality of the naked mole-rat, Heterocephalus glaber. Both molecular systematic studies of Rodentia and the fossil record of bathyergids indicate an ancient origin for the family. This study uses an extensive molecular phylogeny and mitochondrial cytochrome b and 12s rRNA molecular clocks to examine in detail the divergence times, and patterns of speciation of the five extant genera in the context of rift valley formation in Africa. Based on a value of 40-48 million years ago (Myr) for the basal divergence of the family (Heterocephalus), we estimate divergence times of 32-40 Myr for Heliophobius, 20-26 Myr for Georychus/Bathyergus and 12-17 Myr for Cryptomys, the most speciose genus. While early divergences may have been independent of rifting, patterns of distribution of later lineages may have been influenced directly by physical barriers imposed by the formation of the Kenya and Western Rift, and indirectly by accompanying climatic and vegetative changes. Rates of chromosomal evolution and speciation appear to vary markedly within the family. In particular, the genus Cryptomys appears to have undergone an extensive radiation and shows the widest geographical distribution. Of the two distinct clades within this genus, one exhibits considerable karyotypic variation while the other does not, despite comparatively high levels of sequence divergence between some taxa. These different patterns of speciation observed both within the family and within the genus Cryptomys may have been a result of environmental changes associated with rifting.     

The genomic consequences of adaptive divergence and reproductive isolation between species of manakins

The processes of adaptation and speciation are expected to shape genomic variation within and between diverging species. Here we analyze genomic heterogeneity of genetic differentiation and introgression in a hybrid zone between two bird species (Manacus candei and M. vitellinus) using 59 100 SNPs, a whole genome assembly, and Bayesian models. Measures of genetic differentiation () and introgression (genomic cline center [α] and rate [β]) were highly heterogeneous among loci. We identified thousands of loci with elevated parameter estimates, some of which are likely to be associated with variation in fitness in Manacus populations. To analyze the genomic organization of differentiation and introgression, we mapped SNPs onto a draft assembly of the M. vitellinus genome. Estimates of , α, and β were autocorrelated at very short physical distances (< 100 bp), but much less so beyond this. In addition, average statistical associations (linkage disequilibrium) between SNPs were generally low and were not higher in admixed populations than in populations of the parental species. Although they did not occur with a constant probability across the genome, loci with elevated , α, and β were not strongly co-localized in the genome. Contrary to verbal models that predict clustering of loci involved in adaptation and isolation in discrete genomic regions, these results are consistent with the hypothesis that genetic regions involved in adaptive divergence and reproductive isolation are scattered throughout the genome. We also found that many loci were characterized by both exceptional genetic differentiation and introgression, consistent with the hypothesis that loci involved in isolation are also often characterized by a history of divergent selection. However, the concordance between isolation and differentiation was only partial, indicating a complex architecture and history of loci involved in isolation.     

History repeats itself: genomic divergence in copepods

Press stop, erase everything from now till some arbitrary time in the past and start recording life as it evolves once again. Would you see the same tape of life playing itself over and over, or would a different story unfold every time? The late Steven Jay Gould called this experiment replaying the tape of life and argued that any replay of the tape would lead evolution down a pathway radically different from the road actually taken (Gould 1989). This thought experiment has puzzled evolutionary biologists for a long time: how repeatable are evolutionary events? And if history does indeed repeat itself, what are the factors that may help us predict the path taken? A powerful means to address these questions at a small evolutionary scale is to study closely related populations that have evolved independently, under similar environmental conditions. This is precisely what Pereira et al. ( 2016 ) set out to do using marine copepods Tigriopus californicus, and present their results in this issue of Molecular Ecology. They show that evolution can be repeatable and even partly predictable, at least at the molecular level. As expected from theory, patterns of divergence were shaped by natural selection. At the same time, strong genetic drift due to small population sizes also constrained evolution down a similar evolutionary road, and probably contributed to repeatable patterns of genomic divergence.     

Gene flow,linked selection,and divergent sorting of ancient polymorphism shape genomic divergence landscape in a group of edaphic specialists

Interpreting the formation of genomic variation landscape, especially genomic regions with elevated differentiation (i.e. islands), is fundamental to a better understanding of the genomic consequences of adaptation and speciation. Edaphic islands provide excellent systems for understanding the interplay of gene flow and selection in driving population divergence and speciation. However, discerning the relative contribution of these factors that modify patterns of genomic variation remains difficult. We analysed 132 genomes from five recently divergent species in Primulina genus, with four species distributed in Karst limestone habitats and the fifth one growing in Danxia habitats. We demonstrated that both gene flow and linked selection have contributed to genome-wide variation landscape, where genomic regions with elevated differentiation (i.e., islands) were largely derived by divergent sorting of ancient polymorphism. Specifically, we identified several lineage-specific genomic islands that might have facilitated adaptation of P. suichuanensis to Danxia habitats. Our study is amongst the first cases disentangling evolutionary processes that shape genomic variation of plant specialists, and demonstrates the important role of ancient polymorphism in the formation of genomic islands that potentially mediate adaptation and speciation of endemic plants in special soil habitats.     

Circadian genes in a blind subterranean mammal III: molecular cloning and circadian regulation of cryptochrome genes in the blind subterranean mole rat, Spalax ehrenbergi superspecies

The blind subterranean mole rat superspecies Spalax ehrenbergi is an extreme example of mammalian adaptation to life underground. Though this rodent is totally visually blind, harboring a drastically degenerated subcutaneous rudimentary eye, its daily activity rhythm is entrainable to LD cycles. This indicates that it confers light information to the clock, as has been previously shown by the authors in behavioral studies as well as by molecular analyses of its Clock/MOP3 and its three Per genes. The Cryptochrome (Cry) genes found in animals and plants act both as photoreceptors and as essential components of the negative feedback mechanism of the biological clock. To further understand the circadian system of this unique mammal, the authors cloned and characterized the open reading frame of Spalax Cry1 and Cry2. The Spalax CRY1 protein is significantly closer to the human homolog than to the mice one, in contrast to the evolutionary expectations. They have found two isoforms of Cry2 in Spalax, which differ in their 5' end of the open reading frame and defined their expression in Spalax populations. They found a large and significant excess of heterozygotes of sCry2 (sCry2L/S genotype). Both sCry1 and sCry2 mRNAs were found in the SCN, the eye, the harderian gland, as well as in a wide range of peripheral tissues. Their expression pattern under different LD conditions has also been analyzed. As was already shown for other circadian genes, despite being blind and living in darkness, the Cry genes of Spalax behave in a similar, though not identical, pattern as in sighted animals. Once again, the results indicate that the uniquely hypertrophied harderian gland of Spalax plays a key role in its circadian system.     

Mammalian glycosyltransferases: genomic organization and protein structure.

In recent years, several glycosyltransferase genes and cDNAs have been cloned and characterized. Although the glycosyltransferases seem to share the same general architecture, there is only little sequence similarity between the various enzymes. Moreover, a comparison of the organization of the genes shows that there is no common pattern of intron-exon structure. In addition, there seems to be little or no correlation between glycosyltransferase exons and protein domains. Taken together, these observations suggest that many of the glycosyltransferase genes evolved independently. So far, only two glycosyltransferase gene families have been described. These families may have evolved by exon-shuffling, or by gene duplication and subsequent divergence. For specific glycosyltransferases, mechanisms such as alternative splicing and alternative promoter usage play a role in the production of multiple protein isoforms from a single gene. These isoenzymes may differ in their enzymatic properties or cellular localization.     

Identifying genomic hotspots of differentiation and candidate genes involved in the adaptive divergence of pea aphid host races

Identifying the genomic bases of adaptation to novel environments is a long‐term objective in evolutionary biology. Because genetic differentiation is expected to increase between locally adapted populations at the genes targeted by selection, scanning the genome for elevated levels of differentiation is a first step towards deciphering the genomic architecture underlying adaptive divergence. The pea aphid Acyrthosiphon pisum is a model of choice to address this question, as it forms a large complex of plant‐specialized races and cryptic species, resulting from recent adaptive radiation. Here, we characterized genomewide polymorphisms in three pea aphid races specialized on alfalfa, clover and pea crops, respectively, which we sequenced in pools (poolseq). Using a model‐based approach that explicitly accounts for selection, we identified 392 genomic hotspots of differentiation spanning 47.3 Mb and 2,484 genes (respectively, 9.12% of the genome size and 8.10% of its genes). Most of these highly differentiated regions were located on the autosomes, and overall differentiation was weaker on the X chromosome. Within these hotspots, high levels of absolute divergence between races suggest that these regions experienced less gene flow than the rest of the genome, most likely by contributing to reproductive isolation. Moreover, population‐specific analyses showed evidence of selection in every host race, depending on the hotspot considered. These hotspots were significantly enriched for candidate gene categories that control host–plant selection and use. These genes encode 48 salivary proteins, 14 gustatory receptors, 10 odorant receptors, five P450 cytochromes and one chemosensory protein, which represent promising candidates for the genetic basis of host–plant specialization and ecological isolation in the pea aphid complex. Altogether, our findings open new research directions towards functional studies, for validating the role of these genes on adaptive phenotypes.     

Mammalian gene evolution: Nucleotide sequence divergence between mouse and rat

As a paradigm of mammalian gene evolution, the nature and extent of DNA sequence divergence between homologous protein-coding genes from mouse and rat have been investigated. The data set examined includes 363 genes totalling 411 kilobases, making this by far the largest comparison conducted between a single pair of species. Mouse and rat genes are on average 93.4% identical in nucleotide sequence and 93.9% identical in amino acid sequence. Individual genes vary substantially in the extent of nonsynonymous nucleotide substitution, as expected from protein evolution studies; here the variation is characterized. The extent of synonymous (or silent) substitution also varies considerably among genes, though the coefficient of variation is about four times smaller than for nonsynonymous substitutions. A small number of genes mapped to the X-chromosome have a slower rate of molecular evolution than average, as predicted if molecular evolution is male-driven. Base composition at silent sites varies from 33% to 95% G + C in different genes; mouse and rat homologues differ on average by only 1.7% in silent-site G + C, but it is shown that this is not necessarily due to any selective constraint on their base composition. Synonymous substitution rates and silent site base composition appear to be related (genes at intermediate G + C have on average higher rates), but the relationship is not as strong as in our earlier analyses. Rates of synonymous and nonsynonymous substitution are correlated, apparently because of an excess of substitutions involving adjacent pairs of nucleotides. Several factors suggest that synonymous codon usage in rodent genes is not subject to selection.     

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.      

A population genomic analysis of species boundaries: neutral processes, adaptive divergence and introgression between two hybridizing plant species

Little is known about the nature of species boundaries between closely related plant species and about the extent of introgression as a consequence of hybridization upon secondary contact. To address these topics we analyzed genome-wide differentiation between two closely related Silene species, Silene latifolia and S. dioica , and assessed the strength of introgression in sympatry. More than 300 AFLP markers were genotyped in three allopatric and three sympatric populations of each species. Outlier analyses were performed separately for sympatric and allopatric populations. Both positive and negative outlier loci were found, indicating that divergent and balancing selection, respectively, have shaped patterns of divergence between the two species. Sympatric populations of the two species were found to be less differentiated genetically than allopatric populations, indicating that hybridization has led to gene introgression. We conclude that differentiation between S. latifolia and S. dioica has been shaped by a combination of introgression and selection. These results challenge the view that species differentiation is a genome-wide phenomenon, and instead support the idea that genomes can be porous and that species differentiation has a genic basis.     

Exploring the correlations between sequence evolution rate and phenotypic divergence across the Mammalian tree provides insights into adaptive evolution

Sequence evolution behaves in a relatively consistent manner, leading to one of the fundamental paradigms in biology, the existence of a ??molecular clock??. The molecular clock can be distilled to the concept of accumulation of substitutions, through time yielding a stable rate from which we can estimate lineage divergence. Over the last 50?years, evolutionary biologists have obtained an in-depth understanding of this clock??s nuances. It has been fine-tuned by taking into account the vast heterogeneity in rates across lineages and genes, leading to ??relaxed?? molecular clock methods for timetree reconstruction. Sequence rate varies with life history traits including body size, generation time and metabolic rate, and we review recent studies on this topic. However, few studies have explicitly examined correlates between molecular evolution and morphological evolution. The patterns observed across diverse lineages suggest that rates of molecular and morphological evolution are largely decoupled. We discuss how identifying the molecular mechanisms behind rapid functional radiations are central to understanding evolution. The vast functional divergence within mammalian lineages that have relatively ??slow?? sequence evolution refutes the hypotheses that pulses in diversification yielding major phenotypic change are the result of steady accumulation of substitutions. Patterns rather suggest phenotypic divergence is likely caused by regulatory alterations mediated through mechanisms such as insertions/deletions in functional regions. These can rapidly arise and sweep to fixation faster than predicted from a lineage??s sequence neutral substitution rate, enabling species to leapfrog between phenotypic ??islands??. We suggest research directions that could illuminate mechanisms behind the functional diversity we see today.