Sensory systems are attractive evolutionary models to address how organisms adapt to local environments that can cause ecological speciation. However, tests of these evolutionary models have focused on visual, auditory, and olfactory senses. Here, we show local adaptation of bitter taste receptor genes in two neighboring populations of a wild mammal-the blind mole rat Spalax galili-that show ecological speciation in divergent soil environments. We found that basalt-type bitter receptors showed higher response intensity and sensitivity compared with chalk-type ones using both genetic and cell-based functional analyses. Such functional changes could help animals adapted to basalt soil select plants with less bitterness from diverse local foods, whereas a weaker reception to bitter taste may allow consumption of a greater range of plants for animals inhabiting chalk soil with a scarcity of food supply. Our study shows divergent selection on food resources through local adaptation of bitter receptors, and suggests that taste plays an important yet underappreciated role in speciation.

• Klíčová slova
• bitter taste, ecological speciation, functional assay, local adaptation,
• MeSH
• chuť genetika   MeSH
• fyziologická adaptace genetika   MeSH
• savci MeSH
• Spalax * genetika   MeSH
• vznik druhů (genetika) MeSH
• zvířata MeSH
• Check Tag
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Research Support, N.I.H., Extramural MeSH

Speciation mechanisms remain controversial. Two speciation models occur in Israeli subterranean mole rats, genus Spalax: a regional speciation cline southward of four peripatric climatic chromosomal species and a local, geologic-edaphic, genic, and sympatric speciation. Here we highlight their genome evolution. The five species were separated into five genetic clusters by single nucleotide polymorphisms, copy number variations (CNVs), repeatome, and methylome in sympatry. The regional interspecific divergence correspond to Pleistocene climatic cycles. Climate warmings caused chromosomal speciation. Triple effective population size, Ne , declines match glacial cold cycles. Adaptive genes evolved under positive selection to underground stresses and to divergent climates, involving interspecies reproductive isolation. Genomic islands evolved mainly due to adaptive evolution involving ancient polymorphisms. Repeatome, including both CNV and LINE1 repetitive elements, separated the five species. Methylation in sympatry identified geologically chalk-basalt species that differentially affect thermoregulation, hypoxia, DNA repair, P53, and other pathways. Genome adaptive evolution highlights climatic and geologic-edaphic stress evolution and the two speciation models, peripatric and sympatric.

• Klíčová slova
• genomic sequencing, methylation, repeatome, speciation models, subterranean rodents,
• MeSH
• biologická adaptace MeSH
• biologická evoluce * MeSH
• epigeneze genetická MeSH
• genetická variace MeSH
• genom MeSH
• jednonukleotidový polymorfismus MeSH
• molekulární evoluce MeSH
• populační genetika MeSH
• reprodukční izolace MeSH
• Spalax genetika   fyziologie   MeSH
• sympatrie * MeSH
• tok genů MeSH
• variabilita počtu kopií segmentů DNA MeSH
• vazebná nerovnováha MeSH
• zvířata MeSH
• Check Tag
• mužské pohlaví MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Geografické názvy
• Izrael MeSH

Spiny mice, Acomys cahirinus, colonized Israel 30,000 y ago from dry tropical Africa and inhabited rocky habitats across Israel. Earlier, we had shown by mtDNA that A. cahirinus incipiently sympatrically speciates at Evolution Canyon I (EC I) in Mount Carmel, Israel because of microclimatic interslope divergence. The EC I microsite consists of a dry and hot savannoid "African" slope (AS) and an abutting humid and cool-forested "European" slope (ES). Here, we substantiate incipient SS in A. cahirinus at EC I based on the entire transcriptome, showing that multiple slope-specific adaptive complexes across the transcriptome result in two divergent clusters. Tajima's D distribution of the abutting Acomys interslope populations shows that the ES population is under stronger positive selection, whereas the AS population is under balancing selection, harboring higher genetic polymorphisms. Considerable sites of the two populations were differentiated with a coefficient of FST = 0.25-0.75. Remarkably, 24 and 37 putatively adaptively selected genes were detected in the AS and ES populations, respectively. The AS genes involved DNA repair, growth arrest, neural cell differentiation, and heat-shock proteins adapting to the local AS stresses of high solar radiation, drought, and high temperature. In contrast, the ES genes involved high ATP associated with energetics stress. The sharp ecological interslope divergence led to strong slope-specific selection overruling the interslope gene flow. Earlier tests suggested slope-specific mate choice. Habitat interslope-adaptive selection across the transcriptome and mate choice substantiate sympatric speciation (SS), suggesting its prevalence at EC I and commonality in nature.

• Klíčová slova
• RNA-seq, adaptive ecological speciation, microclimate, natural selection,
• MeSH
• biologická evoluce MeSH
• jednonukleotidový polymorfismus MeSH
• mikroklima MeSH
• mozek metabolismus   MeSH
• Murinae genetika   MeSH
• tok genů MeSH
• transkriptom MeSH
• vznik druhů (genetika) * MeSH
• zvířata MeSH
• Check Tag
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• Geografické názvy
• Izrael MeSH

Incipient sympatric speciation in blind mole rat, Spalax galili, in Israel, caused by sharp ecological divergence of abutting chalk-basalt ecologies, has been proposed previously based on mitochondrial and whole-genome nuclear DNA. Here, we present new evidence, including transcriptome, DNA editing, microRNA, and codon usage, substantiating earlier evidence for adaptive divergence in the abutting chalk and basalt populations. Genetic divergence, based on the previous and new evidence, is ongoing despite restricted gene flow between the two populations. The principal component analysis, neighbor-joining tree, and genetic structure analysis of the transcriptome clearly show the clustered divergent two mole rat populations. Gene-expression level analysis indicates that the population transcriptome divergence is displayed not only by soil divergence but also by sex. Gene ontology enrichment of the differentially expressed genes from the two abutting soil populations highlights reproductive isolation. Alternative splicing variation of the two abutting soil populations displays two distinct splicing patterns. L-shaped FST distribution indicates that the two populations have undergone divergence with gene flow. Transcriptome divergent genes highlight neurogenetics and nutrition characterizing the chalk population, and energetics, metabolism, musculature, and sensory perception characterizing the abutting basalt population. Remarkably, microRNAs also display divergence between the two populations. The GC content is significantly higher in chalk than in basalt, and stress-response genes mostly prefer nonoptimal codons. The multiple lines of evidence of ecological-genomic and genetic divergence highlight that natural selection overrules the gene flow between the two abutting populations, substantiating the sharp ecological chalk-basalt divergence driving sympatric speciation.

• Klíčová slova
• DNA editing, ecological adaptive speciation, microRNA regulation, natural selection, nonoptimal codon usage,
• MeSH
• ekosystém MeSH
• mikro RNA metabolismus   MeSH
• půda MeSH
• silikáty MeSH
• Spalax genetika   metabolismus   MeSH
• sympatrie * MeSH
• tok genů MeSH
• transkriptom * MeSH
• uhličitan vápenatý MeSH
• vznik druhů (genetika) * MeSH
• zvířata MeSH
• Check Tag
• mužské pohlaví MeSH
• ženské pohlaví MeSH
• zvířata MeSH
• Publikační typ
• časopisecké články MeSH
• práce podpořená grantem MeSH
• srovnávací studie MeSH
• Názvy látek
• basalt MeSH Prohlížeč
• mikro RNA MeSH
• půda MeSH
• silikáty MeSH
• uhličitan vápenatý MeSH