Genome evolution of blind subterranean mole rats: Adaptive peripatric versus sympatric speciation

Climate warming has divergent physiological impacts on sympatric lizards

Climate warming is expected to affect the vulnerability of sympatric species differentially due to their divergent traits, but the underlying physiological mechanisms of those impacts are poorly understood. We conducted field warming experiments (present climate vs. warm climate) using open-top chambers to determine the effects of climate warming on active body temperature, oxidative damage, immune competence, growth and survival in two sympatric desert-dwelling lizards, Eremias multiocellata and Eremias argus from May 2019 to September 2020. Our climate warming treatment did not affect survival of the two species, but it did increase active body temperatures and growth rate in E. multiocellata compared to E. argus. Climate warming also induced greater oxidative damage (higher malondialdehyde content and catalase activity) in E. multiocellata, but not in E. argus. Further, climate warming increased immune competence in E. multiocellata, but decreased immune competence in E. argus, with regards to white blood cell counts, bacteria killing ability and relative expression of immunoglobulin M. Our results suggest that climate warming enhances body temperature, and thereby oxidative stress, immune competence and growth in E. multiocellata, but decreases immune competence of E. argus, perhaps as a cost of thermoregulation to maintain body temperatures under climate warming. The divergent physiological effects of climate warming on sympatric species may have profound ecological consequences if it eventually leads to changes in reproductive activities, population dynamics and community structure. Our study highlights the importance of considering interspecific differences in physiological traits when we evaluate the impact of climate warming on organisms, even for those closely-related species coexisting within the same geographical area.

Genomic Insights into High-Altitude Adaptation: A Comparative Analysis of Roscoea alpina and R. purpurea in the Himalayas

Environmental stress at high altitudes drives the development of distinct adaptive mechanisms in plants. However, studies exploring the genetic adaptive mechanisms of high-altitude plant species are scarce. In the present study, we explored the high-altitude adaptive mechanisms of plants in the Himalayas through whole-genome resequencing. We studied two widespread members of the Himalayan endemic alpine genus Roscoea (Zingiberaceae): R. alpina (a selfing species) and R. purpurea (an outcrossing species). These species are distributed widely in the Himalayas with distinct non-overlapping altitude distributions; R. alpina is distributed at higher elevations, and R. purpurea occurs at lower elevations. Compared to R. purpurea, R. alpina exhibited higher levels of linkage disequilibrium, Tajima's D, and inbreeding coefficient, as well as lower recombination rates and genetic diversity. Approximately 96.3% of the genes in the reference genome underwent significant genetic divergence (FST ≥ 0.25). We reported 58 completely divergent genes (FST = 1), of which only 17 genes were annotated with specific functions. The functions of these genes were primarily related to adapting to the specific characteristics of high-altitude environments. Our findings provide novel insights into how evolutionary innovations promote the adaptation of mountain alpine species to high altitudes and harsh habitats.

Investigating the evolutionary dynamics of diploid number variation in Ctenomys (Ctenomyidae, Rodentia)

Contrary to predictions from classical hybrid sterility models of chromosomal speciation, some organisms display high rates of karyotype variation. Ctenomys are the current mammals with the greatest interspecific and intraspecific chromosomal variation. A large number of species have been studied cytogenetically. The diploid numbers of chromosomes range from 2n = 10 to 2n = 70. Here, we analyzed karyotype evolution in Ctenomys using comparative phylogenetic methods. We found a strong phylogenetic signal with chromosome number. This refutes the chromosomal megaevolution model, which proposes the independent accumulation of multiple chromosomal rearrangements in each closely related species. We found that Brownian motion (BM) described the observed characteristic changes more thoroughly than the Ornstein-Uhlenbeck and Early-Burst models. This suggests that the evolution of chromosome numbers occurs by a random walk along phylogenetic clades. However, our data indicate that the BM model alone does not fully characterize the chromosomal evolution of Ctenomys.

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

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

Conservation status of the blind mole rat populations in Hungary (Rodentia: Spalacinae: Nannospalax) revisited

Regular reviews of long-term research and conservation programs are useful sources of information for future directions in science and for the assessment of current conservation status of taxa. In this paper, we compiled all available data from the last 10 years related to Hungarian blind mole rat populations and assessed this information according to the following main themes: systematics, distribution and threats, and conservation actions. Based on the most recent information, national and global risk assessments are provided for the three species of Nannospalax (N. montanosyrmiensis, N. hungaricus and N. syrmiensis), currently accepted as part of the Hungarian fauna.

Historical demography and climatic niches of the Natal multimammate mouse (Mastomys natalensis) in the Zambezian region

The Natal multimammate mouse (Mastomys natalensis) is the most widespread rodent species in sub-Saharan Africa, often studied as an agricultural pest and reservoir of viruses. Its mitochondrial (Mt) phylogeny revealed six major lineages parapatrically distributed across open habitats of sub-Saharan Africa. In this study we used 1949 sequences of the mitochondrial cytochrome b gene to elaborate on distribution and evolutionary history of three Mt lineages inhabiting the open habitats of the Zambezian region (corresponding roughly to the African savannas south of the Equator). We describe in more detail contact zones between the lineages—their location and extent of co-occurrence within localities—and infer past population trends. The estimates are interpreted in the light of climatic niche models. The lineages underwent reduction in effective population size during the last glacial, but they spread widely after that: two of them after the last glacial maximum and the last one in mid-Holocene. The centers of expansion, i.e., possible long-term savanna refugia, were estimated to lie close to the Eastern Arc Mountains and lakes of the Great African Rift, geomorphological structures likely to have had long-term influence on geographical distribution of the lineages. Environmental niche modeling shows climate could also affect the broad scale distribution of the lineages but is unlikely to explain the narrow width of the contact zones. The intraspecific Mt differentiation of M. natalensis echoes phylogeographic patterns observed in multiple co-distributed mammal species, which suggests the mammal communities in the region are shaped by the same long-term processes.

Out of the Qinghai-Tibet plateau: Genomic biogeography of the alpine monospecific genus Megadenia (Biscutelleae, Brassicaceae)

Numerous high-elevation alpine plants of the Qinghai-Tibet Plateau (QTP) also have disjunct distribution in adjacent low-altitude mountains. The out-of-QTP versus into-the-QTP hypothesis of alpine plants provide strong evidence for the highly disputed assumption of the massive ice sheet developed in the central plateau during the Last Glacial Maximum (LGM). In this study, we sequenced the genomes of most known populations of Megadenia, a monospecific alpine genus of Brassicaceae distributed primarily in the QTP, though rarely found in adjacent low-elevation mountains of north China and Russia (NC-R). All sequenced samples clustered into four geographic genetic groups: one pair was in the QTP and another was in NC-R. The latter pair is nested within the former, and these findings support the out-of-QTP hypothesis. Dating the four genetic groups and niche distribution suggested that Megadenia migrated out of the QTP to adjacent regions during the LGM. The NC-R group showed a decrease in the effective population sizes. In addition, the genes with high genetic divergences in the QTP group were mainly involved in habitat adaptations during low-altitude colonization. These findings reject the hypothesis of development massive ice sheets, and support glacial survival of alpine plants within, as well as further migration out of, the QTP.

Sympatric or micro-allopatric speciation in a glacial lake? Genomic islands support neither

Apparent cases of sympatric speciation may actually be due to micro-allopatric or micro-parapatric speciation. One way to distinguish between these models is to examine the existence and nature of genomic islands of divergence, wherein divergent DNA segments are interspersed with low-divergence segments. Such islands should be rare or absent under micro-allopatric speciation but common in cases of speciation with gene flow. Sympatric divergence of endemic fishes is known from isolated saline, crater, postglacial, and ancient lakes. Two morphologically distinct cyprinid fishes, Gymnocypris eckloni scoliostomus (GS) and G. eckloni eckloni (GE), in a small glacial lake on the Qinghai-Tibet Plateau, Lake Sunmcuo, match the biogeographic criteria of sympatric speciation. In this study, we examined genome-wide variation in 46 individuals from these two groups. The divergence time between the GS and GE lineages was estimated to be 20-60 Kya. We identified 54 large genomic islands (≥100 kb) of speciation, which accounted for 89.4% of the total length of all genomic islands. These islands harboured divergent genes related to olfactory receptors and olfaction signals that may play important roles in food selection and assortative mating in fishes. Although the genomic islands clearly indicated speciation with gene flow and rejected micro-allopatric speciation, they were too large to support the hypothesis of sympatric speciation. Theoretical and recent empirical studies suggested that continual gene flow in sympatry should give rise to many small genomic islands (as small as a few kilobases in size). Thus, the observed pattern is consistent with the extensive evidence on parapatric speciation, in which adjacent habitats facilitate divergent selection but also permit gene flow during speciation. We suggest that many, if not most, of the reported cases of sympatric speciation are likely to be micro-parapatric speciation.

Sympatric Speciation in Mole Rats and Wild Barley and Their Genome Repeatome Evolution: A Commentary

The theories of sympatric speciation (SS) and coding and noncoding (cd and ncd =repeatome)  genome function are still contentious. Studies on SS in our two new models, "Evolution Canyon" and "Evolution Plateau", in Israel, divergent microclimatically and geologically-edaphically, respectively, indicated that in ecologically divergent microsites SS is a common speciation model across life from bacteria to mammals. Genomically, the intergenic ncd repeatome was and is still regarded by many biologists as "selfish," "junk," and non-functional. In contrast, it is considered by the encyclopedia of DNA elements discovery as biochemically functional and regulatory, and the transposable elements were considered earlier by Barbara McClintock as "controlling elements" of genes. Remarkably, it is found that repeated elements can statistically identify significantly, the five species of subterranean mole rats of Spalax ehrenbergi superspecies adapted to increasingly arid climatic trend southward in Israel. Moreover, it is first discovered in the SS studies in two distant taxa, subterranean mole rats and wild barley, and later also in spiny mice in Israel and subterranean zokors in China, that the noncoding repeatome is genomically mirroring the image of the protein-coding genome in divergent ecologies. It is shown that this mirroring image is statistically significant both within and between the ecologically divergent taxa supporting the hypothesis that much of the repeatome might be regulatory and selected as the protein-coding genome by the same ecological stresses.

Host diet shapes functionally differentiated gut microbiomes in sympatric speciation of blind mole rats in Upper Galilee, Israel

The gut microbiome is important for host nutrient metabolism and ecological adaptation. However, how the gut microbiome is affected by host phylogeny, ecology and diet during sympatric speciation remain unclear. Here, we compare and contrast the gut microbiome of two sympatric blind mole rat species and correlate them with their corresponding host phylogeny, ecology soil metagenomes, and diet to determine how these factors may influence their gut microbiome. Our results indicate that within the host microbiome there is no significant difference in community composition, but the functions between the two sympatric species populations vary significantly. No significant correlations were found between the gut microbiome differentiation and their corresponding ecological soil metagenomes and host phylogeny. Functional enrichment analysis suggests that the host diets may account for the functional divergence of the gut microbiome. Our results will help us understand how the gut microbiome changes with corresponding ecological dietary factors in sympatric speciation of blind subterranean mole rats.

Sympatric genetic divergence between early- and late-season weedy rice populations

Sympatric genetic divergence is the most appealing and controversial pattern in the theory of ecological speciation. Examples that support sympatric genetic divergence in plant species are extremely rare. Solid evidence of sympatric genetic divergence will provide deep insights for revealing the underlying mechanisms of ecological speciation. We analysed the total genomic DNA sequences of 120 weedy rice (WR; Oryza sativa f. spontanea) plants, representing three WR population pairs separately from three early- and late-season rice fields, in comparison with those of the co-occurring rice cultivars and other rice materials. We detected substantial genetic divergence within the pairs of the sympatric early- and late-season WR populations, although genetic divergence was unevenly distributed across the genomes. Restricted gene flow was determined between the sympatric WR populations, resulting in their distinct genetic structures. We also detected relatively low genetic diversity that was likely to be associated with stronger selection in early-season WR populations. Our findings provide strong evidence for sympatric genetic divergence between the WR populations in the same fields but in different seasons. We conclude that temporal isolation plays an important role in creating genetic divergence between sympatric populations/species in plants.

An Association Study of DNA Methylation and Gene Expression in Angelman Syndrome: A Bioinformatics Approach

Angelman syndrome (AS) is a neurodevelopmental disorder caused by the loss of function of the E3-ligase UBE3A. Despite multiple studies, AS pathophysiology is still obscure and has mostly been explored in rodent models of the disease. In recent years, a growing body of studies has utilized omics datasets in the attempt to focus research regarding the pathophysiology of AS. Here, for the first time, we utilized a multi-omics approach at the epigenomic level and the transcriptome level, for human-derived neurons. Using publicly available datasets for DNA methylation and gene expression, we found genome regions in proximity to gene promoters and intersecting with gene-body regions that were differentially methylated and differentially expressed in AS. We found that overall, the genome in AS postmortem brain tissue was hypo-methylated compared to healthy controls. We also found more upregulated genes than downregulated genes in AS. Many of these dysregulated genes in neurons obtained from AS patients are known to be critical for neuronal development and synaptic functioning. Taken together, our results suggest a list of dysregulated genes that may be involved in AS development and its pathological features. Moreover, these genes might also have a role in neurodevelopmental disorders similar to AS.

Incipient Sympatric Speciation and Evolution of Soil Bacteria Revealed by Metagenomic and Structured Non-Coding RNAs Analysis

Simple Summary

The microevolutionary dynamics of soil bacteria under microclimatic differences are largely unexplored in contrast to our improving knowledge of their vast diversity. In this study, we performed a comparative metagenomic analysis of two sharply divergent rocks and soil types at the Evolution Plateau (EP) in eastern Upper Galilee, Israel. We have identified the significant differences in bacterial taxonomic diversity, functions, and patterns of RNA-based gene regulation between the bacteria from two different soil types. Furthermore, we have identified several species with a significant genetic divergence of the same species between the two soil types, highlighting the soil bacteria’s incipient sympatric speciation.

Abstract

Soil bacteria respond rapidly to changes in new environmental conditions. For adaptation to the new environment, they could mutate their genome, which impacts the alternation of the functional and regulatory landscape. Sometimes, these genetic and ecological changes may drive the bacterial evolution and sympatric speciation. Although sympatric speciation has been controversial since Darwin suggested it in 1859, there are several strong theoretical or empirical evidences to support it. Sympatric speciation associated with soil bacteria remains largely unexplored. Here, we provide potential evidence of sympatric speciation of soil bacteria by comparison of metagenomics from two sharply contrasting abutting divergence rock and soil types (Senonian chalk and its rendzina soil, and abutting Pleistocene basalt rock and basalt soil). We identified several bacterial species with significant genetic differences in the same species between the two soil types and ecologies. We show that the bacterial community composition has significantly diverged between the two soils; correspondingly, their functions were differentiated in order to adapt to the local ecological stresses. The ecologies, such as water availability and pH value, shaped the adaptation and speciation of soil bacteria revealed by the clear-cut genetic divergence. Furthermore, by a novel analysis scheme of riboswitches, we highlight significant differences in structured non-coding RNAs between the soil bacteria from two divergence soil types, which could be an important driver for functional adaptation. Our study provides new insight into the evolutionary divergence and incipient sympatric speciation of soil bacteria under microclimatic ecological differences.

Genomic insights into zokors' phylogeny and speciation in China

SignificancePhylogenies are the basis of many ecological and evolutionary studies. However, zokor phylogeny and speciation patterns are heavily debated. This study disentangled the phylogeny and speciation of zokors genomically. Six species of the Eospalax were separated into high-altitude E. baileyi and E. smithi and the rest four low-altitude species by recent Qinghai-Tibet Plateau uplift 3.6 million y ago. E. rothschildi and E. smithi speciated south of the Qinling-Huaihe Line, where refuges were supplied during glaciation. Introgression and incomplete lineage sorting led to the complex phylogeny of zokors. Genomic islands were formed due to ancient polymorphisms and divergence hitchhiking. This study concluded that climatic, geological, and tectonic events shaped the phylogeny and speciation of zokors in China.

Adaptation of mammals to hypoxia

Oxygen plays a pivotal role in the metabolism and activities of mammals. However, oxygen is restricted in some environments-subterranean burrow systems or habitats at high altitude or deep in the ocean-and this could exert hypoxic stresses such as oxidative damage on organisms living in these environments. In order to cope with these stresses, organisms have evolved specific strategies to adapt to hypoxia, including changes in physiology, gene expression regulation, and genetic mutations. Here, we review how mammals have adapted to the three high-altitude plateaus of the world, the limited oxygen dissolved in deep water habitats, and underground tunnels, with the aim of better understanding the adaptation of mammals to hypoxia.

Comparative Genomics Analysis of Repetitive Elements in Ten Gymnosperm Species: "Dark Repeatome" and Its Abundance in Conifer and Gnetum Species

Repetitive elements (RE) and transposons (TE) can comprise up to 80% of some plant genomes and may be essential for regulating their evolution and adaptation. The “repeatome” information is often unavailable in assembled genomes because genomic areas of repeats are challenging to assemble and are often missing from final assembly. However, raw genomic sequencing data contain rich information about RE/TEs. Here, raw genomic NGS reads of 10 gymnosperm species were studied for the content and abundance patterns of their “repeatome”. We utilized a combination of alignment on databases of repetitive elements and de novo assembly of highly repetitive sequences from genomic sequencing reads to characterize and calculate the abundance of known and putative repetitive elements in the genomes of 10 conifer plants: Pinus taeda, Pinus sylvestris, Pinus sibirica, Picea glauca, Picea abies, Abies sibirica, Larix sibirica, Juniperus communis, Taxus baccata, and Gnetum gnemon. We found that genome abundances of known and newly discovered putative repeats are specific to phylogenetically close groups of species and match biological taxa. The grouping of species based on abundances of known repeats closely matches the grouping based on abundances of newly discovered putative repeats (kChains) and matches the known taxonomic relations.