Genetic adaptations of the plateau zokor in high-elevation burrows

Synchronized Expansion and Contraction of Olfactory, Vomeronasal, and Taste Receptor Gene Families in Hystricomorph Rodents

Chemical senses, including olfaction, pheromones, and taste, are crucial for the survival of most animals. There has long been a debate about whether different types of senses might influence each other. For instance, primates with a strong sense of vision are thought to have weakened olfactory abilities, although the oversimplified trade-off theory is now being questioned. It is uncertain whether such interactions between different chemical senses occur during evolution. To address this question, we examined four receptor gene families related to olfaction, pheromones, and taste: olfactory receptor (OR), vomeronasal receptor type 1 and type 2 (V1R and V2R), and bitter taste receptor (T2R) genes in Hystricomorpha, which is morphologically and ecologically the most diverse group of rodents. We also sequenced and assembled the genome of the grasscutter, Thryonomys swinderianus. By examining 16 available genome assemblies alongside the grasscutter genome, we identified orthologous gene groups among hystricomorph rodents for these gene families to separate the gene gain and loss events in each phylogenetic branch of the Hystricomorpha evolutionary tree. Our analysis revealed that the expansion or contraction of the four gene families occurred synchronously, indicating that when one chemical sense develops or deteriorates, the others follow suit. The results also showed that V1R/V2R genes underwent the fastest evolution, followed by OR genes, and T2R genes were the most evolutionarily stable. This variation likely reflects the difference in ligands of V1R/V2Rs, ORs, and T2Rs: species-specific pheromones, environment-based scents, and toxic substances common to many animals, respectively.

Coping with extremes: Alternations in diet, gut microbiota, and hepatic metabolic functions in a highland passerine

In wild animals, diet and gut microbiota interactions are critical moderators of metabolic functions and are highly contingent on habitat conditions. Challenged by the extreme conditions of high-altitude environments, the strategies implemented by highland animals to adjust their diet and gut microbial composition and modulate their metabolic substrates remain largely unexplored. By employing a typical human commensal species, the Eurasian tree sparrow (Passer montanus, ETS), as a model species, we studied the differences in diet, digestive tract morphology and enzyme activity, gut microbiota, and metabolic energy profiling between highland (the Qinghai-Tibet Plateau, QTP; 3230 m) and lowland (Shijiazhuang, Hebei; 80 m) populations. Our results showed that highland ETSs had enlarged digestive organs and longer small intestinal villi, while no differences in key digestive enzyme activities were observed between the two populations. The 18S rRNA sequencing results revealed that the dietary composition of highland ETSs were more animal-based and less plant-based than those of the lowland ones. Furthermore, 16S rRNA sequencing results suggested that the intestinal microbial communities were structurally segregated between populations. PICRUSt metagenome predictions further indicated that the expression patterns of microbial genes involved in material and energy metabolism, immune system and infection, and xenobiotic biodegradation were strikingly different between the two populations. Analysis of liver metabolomics revealed significant metabolic differences between highland and lowland ETSs in terms of substrate utilization, as well as distinct sex-specific alterations in glycerophospholipids. Furthermore, the interplay between diet, liver metabolism, and gut microbiota suggests a dietary shift resulting in corresponding changes in gut microbiota and metabolic functions. Our findings indicate that highland ETSs have evolved to optimize digestion and absorption, rely on more protein-rich foods, and possess gut microbiota tailored to their dietary composition, likely adaptive physiological and ecological strategies adopted to cope with extreme highland environments.

Mosaic environment-driven evolution of the deep-sea mussel Gigantidas platifrons bacterial endosymbiont

BACKGROUND

The within-species diversity of symbiotic bacteria represents an important genetic resource for their environmental adaptation, especially for horizontally transmitted endosymbionts. Although strain-level intraspecies variation has recently been detected in many deep-sea endosymbionts, their ecological role in environmental adaptation, their genome evolution pattern under heterogeneous geochemical environments, and the underlying molecular forces remain unclear.

RESULTS

Here, we conducted a fine-scale metagenomic analysis of the deep-sea mussel Gigantidas platifrons bacterial endosymbiont collected from distinct habitats: hydrothermal vent and methane seep. Endosymbiont genomes were assembled using a pipeline that distinguishes within-species variation and revealed highly heterogeneous compositions in mussels from different habitats. Phylogenetic analysis separated the assemblies into three distinct environment-linked clades. Their functional differentiation follows a mosaic evolutionary pattern. Core genes, essential for central metabolic function and symbiosis, were conserved across all clades. Clade-specific genes associated with heavy metal resistance, pH homeostasis, and nitrate utilization exhibited signals of accelerated evolution. Notably, transposable elements and plasmids contributed to the genetic reshuffling of the symbiont genomes and likely accelerated adaptive evolution through pseudogenization and the introduction of new genes.

CONCLUSIONS

The current study uncovers the environment-driven evolution of deep-sea symbionts mediated by mobile genetic elements. Its findings highlight a potentially common and critical role of within-species diversity in animal-microbiome symbioses. Video Abstract.

Divergent contributions of coding and noncoding sequences to initial high-altitude adaptation in passerine birds endemic to the Qinghai-Tibet Plateau

Early events in the evolution of an ancestral lineage can shape the adaptive patterns of descendant species, but the evolutionary mechanisms driving initial adaptation from an ancestor remain largely unexplored. High-altitude adaptations have been extensively explored from the viewpoint of protein-coding genes; however, the contribution of noncoding regions remains relatively neglected. Here, we integrate genomic and transcriptomic data to investigate adaptive evolution in the ancestor of three high-altitude snowfinch species endemic to the Qinghai-Tibet Plateau. Our genome-wide scan for adaptation in the snowfinch ancestor identifies strong adaptation signals in functions of development and metabolism for the coding genes, but in functions of the nervous system development for noncoding regions. This pattern is exclusive to the snowfinch ancestor compared to a control ancestral lineage subject to weak selection. Changes in noncoding regions in the snowfinch ancestor, especially those nearest to coding genes, may be disproportionately associated with the differential expression of genes in the brain tissue compared to other tissues. Extensive gene expression in the brain tissue can be further altered via genetic regulatory networks of transcription factors harbouring potential accelerated regulatory regions (e.g., the development-related transcription factor YEATS4). Altogether, our study provides new evidence concerning how coding and noncoding sequences work through decoupled pathways in initial adaptation to the selective pressure of high-altitude environments. The analysis highlights the idea that noncoding sequences may be promising elements in facilitating the rapid evolution and adaptation to high altitudes.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

How voles adapt to subterranean lifestyle: Insights from RNA-seq

Life under the earth surface is highly challenging and associated with a number of morphological, physiological and behavioral modifications. Subterranean niche protects animals from predators, fluctuations in environmental parameters, but is characterized by high levels of carbon dioxide and low levels of oxygen and implies high energy requirements associated with burrowing. Moreover, it lacks most of the sensory inputs available above ground. The current study describes results from RNA-seq analysis of four subterranean voles from subfamily Arvicolinae: Prometheomys schaposchnikowi, Ellobius lutescens, Terricola subterraneus, and Lasiopodomys mandarinus. Original RNA-seq data were obtained for eight species, for nine species, SRA data were downloaded from the NCBI SRA database. Additionally assembled transcriptomes of Mynomes ochrogaster and Cricetulus griseus were included in the analysis. We searched for the selection signatures and parallel amino acid substitutions in a total of 19 species. Even within this limited data set, we found significant changes of dN/dS ratio by free-ratio model analysis for subterranean Arvicolinae. Parallel substitutions were detected in genes RAD23B and PYCR2. These genes are associated with DNA repair processes and response to oxidative stress. Similar substitutions were discovered in the RAD23 genes for highly specialized subterranean Heterocephalus glaber and Fukomys damarensis. The most pronounced signatures of adaptive evolution related to subterranean niche within species of Arvicolinae subfamily were detected for Ellobius lutescens. Our results suggest that genomic adaptations can occur very quickly so far as the amount of selection signatures was found to be compliant with the degree of specialization to the subterranean niche and independent from the evolutionary age of the taxon. We found that the number of genomic signatures of selection does not depend on the age of the taxon, but is positively correlated with the degree of specialization to the subterranean niche.

Gut microbiota facilitates adaptation of the plateau zokor (Myospalax baileyi) to the plateau living environment

Gut microbiota not only helps the hosts to perform many key physiological functions such as food digestion, energy harvesting and immune regulation, but also influences host ecology and facilitates adaptation of the host to extreme environments. Plateau zokors epitomize successful physiological adaptation to their living environment in the face of the harsh environment characterized by low temperature, low pressure and hypoxia in the Tibetan plateau region and high concentrations of CO₂ in their burrows. Therefore, here we used a metagenomic sequencing approach to explore how gut microbiota contributed to the adaptive evolution of the plateau zokor on the Qinghai-Tibet Plateau. Our metagenomic results show that the gut microbiota of plateau zokors on the Tibetan plateau is not only enriched in a large number of species related to energy metabolism and production of short-chain fatty acids (SCFAs), but also significantly enriched the KO terms that involve carbohydrate uptake pathways, which well address energy uptake in plateau zokors while also reducing inflammatory responses due to low pressure, hypoxia and high CO₂ concentrations. There was also a significant enrichment of tripeptidyl-peptidase II (TPPII) associated with antigen processing, apoptosis, DNA damage repair and cell division, which may facilitate the immune response and tissue damage repair in plateau zokors under extreme conditions. These results suggest that these gut microbiota and their metabolites together contribute to the physiological adaptation of plateau zokors, providing new insights into the contribution of the microbiome to the evolution of mammalian adaptation.

Divergence in interspecific and intersubspecific gene expression between two closely related horseshoe bats (Rhinolophus)

Closely related species have been used as representative systems to investigate the genetic mechanisms involved in the early stages of species differentiation. Previous studies have indicated that variation in gene expression might be a sensitive indicator of initial species divergence, although the role of expression divergence, and especially that associated with phenotypic variation remained relatively undefined. For three organs (cochlea, brain, and liver) from two closely related bat species (Rhinolophus siamensis and R. episcopus), the interspecific and intersubspecific gene expression profiles were compared using transcriptomics in this study. Striking organ specificity of expression was observed, and expression profiles exhibited similarities between cochlea and brain tissues. Numerous differentially expressed genes (DEGs) were identified for each organ in the interspecific comparison (cochlea/brain/liver: 1,069/647/692) and intersubspecific comparison (608/528/368). Functional enrichment analysis indicated vital variation in expression related to the immune system, ion activities, neuronal function, and multisensory system regulation in both comparisons. DEGs relevant to the variation in echolocation calls (RF) were found, and some of them were involved in the pivotal patterns of expression variation. The regulation of immune, ion channel, neural activity, and sophisticated sensory functions at the expression level might be key mechanisms in the early species divergence of bats, and the expression variation related to acoustical signal could have played a crucial part. This study expands our knowledge of gene expression and patterns of variation for three key organs to echolocation at both the interspecific and intersubspecific levels. Further, the framework described here provides insight into the genetic basis of phenotypic variation during the incipient stage of species differentiation.

Dynamic changes in the home range of the subterranean rodent Myospalax baileyi

As ecosystem engineers, subterranean rodents excavate and inhabit burrow systems. However, the changes in their use of underground space are poorly recorded. There is conflicting evidence about whether the burrow systems of subterranean rodents, once established, are relatively stable as a result of the high energy costs of digging. We monitored the size of the home ranges of the plateau zokor (Myospalax baileyi) during different stages of its life cycle to show whether mating behavior and the characteristics of its habitat influence the size and location of its home range. We used radio-tracking to quantify the changes in, and overlap of, the home range of M. baileyi during a one-year period. The average size of the home ranges of male zokors was 6.5 times larger than that of female zokors during the mating season. The males expanded their burrows to overlap with multiple females to increase their chances of mating. However, there was no overlap between estrus females or males, perhaps to reduce the number of encounters and unnecessary fights. The home ranges of male and female zokors were similar in size after courtship and the home ranges of single zokors overlapped with those of several neighbors. Most individuals remained territorial and excluded intraspecific interactions from their home ranges. The location of female zokors was stable throughout the year, but half of the males changed the location of their nests and established completely new home ranges in the non-breeding season, mainly in October. The use of space by M. baileyi is flexible in response to a need for physical contact during the mating season and food resources. The home ranges of subterranean plateau zokors are dynamic and the home ranges of male zokors can change within one breeding cycle.

Functional genomics analysis reveals the evolutionary adaptation and demographic history of pygmy lorises

Lorises are a group of globally threatened strepsirrhine primates that exhibit many unusual physiological and behavioral features, including a low metabolic rate, slow movement, and hibernation. Here, we assembled a chromosome-level genome sequence of the pygmy loris (Xanthonycticebus pygmaeus) and resequenced whole genomes from 50 pygmy lorises and 6 Bengal slow lorises (Nycticebus bengalensis). We found that many gene families involved in detoxification have been specifically expanded in the pygmy loris, including the GSTA gene family, with many newly derived copies functioning specifically in the liver. We detected many genes displaying evolutionary convergence between pygmy loris and koala, including PITRM1. Significant decreases in PITRM1 enzymatic activity in these two species may have contributed to their characteristic low rate of metabolism. We also detected many evolutionarily convergent genes and positively selected genes in the pygmy loris that are involved in muscle development. Functional assays demonstrated the decreased ability of one positively selected gene, MYOF, to up-regulate the fast-type muscle fiber, consistent with the lower proportion of fast-twitch muscle fibers in the pygmy loris. The protein product of another positively selected gene in the pygmy loris, PER2, exhibited weaker binding to the key circadian core protein CRY, a finding that may be related to this species' unusual circadian rhythm. Finally, population genomics analysis revealed that these two extant loris species, which coexist in the same habitat, have exhibited an inverse relationship in terms of their demography over the past 1 million years, implying strong interspecies competition after speciation.

The Complex and Well-Developed Morphological and Histological Structures of the Gastrointestinal Tract of the Plateau Zokor Improve Its Digestive Adaptability to High-Fiber Foods

The morphological and histological traits of the gastrointestinal tract (GIT) enable the animal to perform some specific functions that enhance the species’ adaptability to environments. The plateau zokor (Eospalax baileyi) is a subterranean rodent that mainly forages on plant roots in the Qinghai-Tibet Plateau, but little is known about the mechanism by which the plateau zokor digests roots that have high fiber contents. In this study, we used comparative anatomy methods to compare the morphological and histological traits of the GIT of both the plateau zokor and the plateau pika (Ochotona curzoniae), a small, fossorial lagomorph that forages aboveground plant parts, in order to clarify the traits of the plateau zokor’s GIT and to understand its adaptations to high-fiber foods. The results showed that the foods which plateau zokors eat have a higher fiber content than those which the plateau pikas eat. The plateau zokor has a double-chambered and hemi-glandular stomach (the tubular glands are only in the gastric corpus II, and the gastric fundus is keratinized), whereas the plateau pika has a simple, wholly glandular stomach. The gross morphological indicators (organ index and relative length) of the GIT were significantly lower in the plateau zokor than they were in the plateau pika (p < 0.001). However, the thickness of the gastric corpus II mucosal layer and the gastric fundus muscle layer are significantly higher in the plateau zokor than they are in the plateau pika (p < 0.001), and the thickness of each layer of intestinal tissue is higher in the plateau zokor than it is in the plateau pika. Additionally, the small intestinal villi also are higher and wider in the plateau zokor than they are in the plateau pika. Our results suggest that instead of adapting to digest the high-fiber diet by expanding the size of the GIT, the plateau zokor has evolved a complex stomach and a well-developed gastrointestinal histological structure, and that these specialized GIT structures are consistent with an optimal energy-economy evolutionary adaptation strategy.

Vitamin D3 Metabolic Enzymes in Plateau Zokor (Myospalax baileyi) and Plateau Pika (Ochotona curzoniae): Expression and Response to Hypoxia

Vitamin D3 (D3) is produced endogenously from 7-dehydrocholesterol by irradiation and is an important secosteroid for the absorption of calcium and phosphate. Lithocholic acid (LCA) increases intestinal paracellular calcium absorption in a vitamin D receptor-dependent manner in vitamin D-deficient rats. The plateau zokor (Myospalax baileyi), a strictly subterranean species, and plateau pika are endemic to the Qinghai-Tibet Plateau. To verify whether the zokors were deficient in D3 and reveal the effects of hypoxia on D3 metabolism in the zokors and pikas, we measured the levels of 25(OH)D3, calcium, and LCA, and quantified the expression levels of D3 metabolism-related genes. The results showed an undetectable serum level of 25(OH)D3 and a significantly higher concentration of LCA in the serum of plateau zokor, but its calcium concentration was within the normal range compared with that of plateau pika and Sprague-Dawley rats. With increasing altitude, the serum 25(OH)D3 levels in plateau pika decreased significantly, and the mRNA and protein levels of CYP2R1 (in the liver) and CYP27B1 (in the kidney) in plateau pika decreased significantly. Our results indicate that plateau zokors were deficient in D3 and abundant in LCA, which might be a substitution of D3 in the zokor. Furthermore, hypoxia suppresses the metabolism of D3 by down-regulating the expression of CYP2R1 and CYP27B1 in plateau pika.

Response to Different Oxygen Partial Pressures and Evolution Analysis of Apoptosis-Related Genes in Plateau Zokor (Myospalax baileyi)

The plateau zokor (Myospalax baileyi) is a native species of the Qinghai–Tibet Plateau that spends its entire life underground in sealed burrows with hypoxic conditions. The present study aimed to assess the sequence characteristics of apoptosis-related genes and the response to different oxygen partial pressures (pO₂) in plateau zokor and Sprague-Dawley rats. The sequences of the p53-induced protein with a death domain (Pidd), p53-upregulated modulator of apoptosis (Puma), insulin-like growth factor binding protein 3 (Igfbp3), and apoptosis protease-activating factor 1 (Apaf1) were evaluated concerning homology and convergent evolution sites, and their mRNA levels were evaluated in different tissues under 14.13 (3,300 m) and 16.12 kPa (2,260 m) pO₂ conditions. Our results showed that, (1) the sequences of the apoptosis-related genes in plateau zokor were highly similar to those of Nannospalax galili, followed by Rattus norvegicus; (2). Pidd, Puma, Igfbp3, and Apaf1 of plateau zokor were found to have five, one, two, and five convergent sites in functional domains with N. galili, respectively. Lastly (3), under low pO₂, the expression of Pidd and Puma was downregulated in the lung of plateau zokors. In turn, Igfbp3 and Apaf1 were upregulated in the liver and lung, and Puma was upregulated in the skeletal muscle of plateau zokor under low pO₂. In Sprague-Dawley rats, low pO₂ downregulated Puma and Apaf1 expression in the liver and downregulated Igfbp3 and Puma in the lung and skeletal muscle separately. In contrast, low pO₂ upregulated Pidd expression in the liver and skeletal muscle of Sprague-Dawley rats. Overall, the expression patterns of Apaf1, Igfbp3, and Puma showed the opposite pattern in the liver, lung, and skeletal muscle, respectively, of plateau zokor as compared with Sprague-Dawley rats. In conclusion, for the long-time adaptation to hypoxic environments, Pidd, Puma, Igfbp3, and Apaf1 of plateau zokor underwent convergent evolution, which we believe may have led to upregulation of their levels under low oxygen partial pressures to induce apoptosis, so as to suppress tumorigenesis under hypoxic environments in plateau zokor.

A New Homotetramer Hemoglobin in the Pulmonary Surfactant of Plateau Zokors (Myospalax Baileyi)

The plateau zokor (Myospalax baileyi) is a native species to the Qinghai-Tibetan Plateau, inhabiting hypoxia and hypercapnia sealed subterranean burrows that pose several unique physiological challenges. In this study, we observed a novel heme-containing protein in the pulmonary surfactant (PS) of plateau zokor, identified the encoding gene of the protein, predicted its origination and structure, verified its expression in alveolar epithelial cells, and determined the protein’s affinity to oxygen and its effect on the oxygen-dissolving capability in the PS of plateau zokors. The protein is an unusual homotetramer hemoglobin consisting of four γ-like subunits, and the subunit is encoded by a paralog gene of γ, that is γ-like. The divergence time of γ-like from γ is estimated by the molecular clock to be about 2.45 Mya. The generation of γ-like in plateau zokors might well relate to long-time stress of the high land hypoxia. Unlike γ, the γ-like has a hypoxia response element (HRE) and a lung tissue-specific enhancer in its upstream region, and it is expressed specifically in lung tissues and up-regulated by hypoxia. The protein is named as γ₄-like which is expressed specifically in Alveolar epithelial type II (ATII) cells and secreted into the alveolar cavities through the osmiophilic multilamellar body (LBs). The γ₄-like has a higher affinity to oxygen, and that increases significantly oxygen-dissolving capability in the PS of plateau zokors by its oxygenation function, which might be beneficial for the plateau zokors to obtain oxygen from the severe hypoxia environments by facilitating oxygen diffusion from alveoli to blood.

Gene losses may contribute to subterranean adaptations in naked mole-rat and blind mole-rat

Background

Naked mole-rats (Heterocephalus glaber, NMRs) and blind mole-rats (Spalax galili, BMRs) are representative subterranean rodents that have evolved many extraordinary traits, including hypoxia tolerance, longevity, and cancer resistance. Although multiple candidate loci responsible for these traits have been uncovered by genomic studies, many of them are limited to functional changes to amino acid sequence and little is known about the contributions of other genetic events. To address this issue, we focused on gene losses (unitary pseudogenes) and systematically analyzed gene losses in NMRs and BMRs, aiming to elucidate the potential roles of pseudogenes in their adaptation to subterranean lifestyle.

Results

We obtained the pseudogene repertoires in NMRs and BMRs, as well as their respective aboveground relatives, guinea pigs and rats, on a genome-wide scale. As a result, 167, 139, 341, and 112 pseudogenes were identified in NMRs, BMRs, guinea pigs, and rats, respectively. Functional enrichment analysis identified 4 shared and 2 species-specific enriched functional groups (EFGs) in subterranean lineages. Notably, the pseudogenes in these EFGs might be associated with either regressive (e.g., visual system) or adaptive (e.g., altered DNA damage response) traits. In addition, several pseudogenes including TNNI3K and PDE5A might be associated with specific cardiac features observed in subterranean lineages. Interestingly, we observed 20 convergent gene losses in NMRs and BMRs. Given that the functional investigations of these genes are generally scarce, we provided functional evidence that independent loss of TRIM17 in NMRs and BMRs might be beneficial for neuronal survival under hypoxia, supporting the positive role of eliminating TRIM17 function in hypoxia adaptation. Our results also suggested that pseudogenes, together with positively selected genes, reinforced subterranean adaptations cooperatively.

Conclusions

Our study provides new insights into the molecular underpinnings of subterranean adaptations and highlights the importance of gene losses in mammalian evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01243-0.

Long-read genome sequencing provides molecular insights into scavenging and societal complexity in spotted hyena Crocuta crocuta

The spotted hyena (Crocuta crocuta) is a large and unique terrestrial carnivore. It is a particularly fascinating species due to its distinct phenotypic traits, especially its complex social structure and scavenging lifestyle, with associated high dietary exposure to microbial pathogens. However, the underlying molecular mechanisms related to these phenotypes remain elusive. Here, we sequenced and assembled a high-quality long-read genome of the spotted hyena, with a contig N50 length of ∼13.75 Mb. Based on comparative genomics, immunoglobulin family members (e.g., IGKV4-1) showed significant adaptive duplications in the spotted hyena and striped hyena. Furthermore, immune-related genes (e.g., CD8A, LAG3, and TLR3) experienced species-specific positive selection in the spotted hyena lineage. These results suggest that immune tolerance between the spotted hyena and closely related striped hyena has undergone adaptive divergence to cope with prolonged dietary exposure to microbial pathogens from scavenging. Furthermore, we provided the potential genetic insights underlying social complexity, hinting at social behavior and cognition. Specifically, the RECNE-associated genes (e.g., UGP2 and ACTR2) in the spotted hyena genome are involved in regulation of social communication. Taken together, our genomic analyses provide molecular insights into the scavenging lifestyle and societal complexity of spotted hyenas.

Chromosome-Level Reference Genome Assembly for the American Pika (Ochotona princeps)

The American pika (Ochotona princeps) is an alpine lagomorph found throughout western North America. Primarily inhabiting talus slopes at higher elevations (>2000 m), American pikas are well adapted to cold, montane environments. Warming climates on both historical and contemporary scales have contributed to population declines in American pikas, positioning them as a focal mammalian species for investigating the ecological effects of climate change. To support and expand ongoing research efforts, here, we present a highly contiguous and annotated reference genome assembly for the American pika (OchPri4.0). This assembly was produced using Dovetail de novo proximity ligation methods and annotated through the NCBI Eukaryotic Genome Annotation pipeline. The resulting assembly was chromosome- scale, with a total length of 2.23 Gb across 9350 scaffolds and a scaffold N50 of 75.8 Mb. The vast majority (>97%) of the total assembly length was found within 36 large scaffolds; 33 of these scaffolds correlated to whole autosomes, while the X chromosome was covered by 3 large scaffolds. Additionally, we identified 17 enriched gene ontology terms among American pika-specific genes putatively related to adaptation to high-elevation environments. This high-quality genome assembly will serve as a springboard for exploring the evolutionary underpinnings of behavioral, ecological, and taxonomic diversification in pikas as well as broader-scale eco-evolutionary questions pertaining to cold-adapted species in general.

Searching for signatures of positive selection in cytochrome b gene associated with subterranean lifestyle in fast-evolving arvicolines (Arvicolinae, Cricetidae, Rodentia)

Background

Mitochondrial genes encode proteins involved in oxidative phosphorylation. Variations in lifestyle and ecological niche can be directly reflected in metabolic performance. Subterranean rodents represent a good model for testing hypotheses on adaptive evolution driven by important ecological shifts. Voles and lemmings of the subfamily Arvicolinae (Rodentia: Cricetidae) provide a good example for studies of adaptive radiation. This is the youngest group within the order Rodentia showing the fastest rates of diversification, including the transition to the subterranean lifestyle in several phylogenetically independent lineages.

Results

We evaluated the signatures of selection in the mitochondrial cytochrome b (cytB) gene in 62 Arvicolinae species characterized by either subterranean or surface-dwelling lifestyle by assessing amino acid sequence variation, exploring the functional consequences of the observed variation in the tertiary protein structure, and estimating selection pressure. Our analysis revealed that: (1) three of the convergent amino acid substitutions were found among phylogenetically distant subterranean species and (2) these substitutions may have an influence on the protein complex structure, (3) cytB showed an increased ω and evidence of relaxed selection in subterranean lineages, relative to non-subterranean, and (4) eight protein domains possess increased nonsynonymous substitutions ratio in subterranean species.

Conclusions

Our study provides insights into the adaptive evolution of the cytochrome b gene in the Arvicolinae subfamily and its potential implications in the molecular mechanism of adaptation. We present a framework for future characterizations of the impact of specific mutations on the function, physiology, and interactions of the mtDNA-encoded proteins involved in oxidative phosphorylation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01819-4.

Adaptive evolution of β-globin gene in subterranean in South America octodontid rodents

The convergent evolution of subterranean rodents is an excellent model to study how natural selection operates and the genetic bases of these adaptations, but the study on the different taxa has been very uneven and still insufficient. In the octodontoid caviomorph rodent superfamily there are two independent lineages where they have recently evolved into totally underground lifestyles: the genera Ctenomys (tuco-tucos) and Spalacopus (coruro). The underground habitat is characterized by an hypoxic and hypercapnic atmosphere, thus gas exchange is one of the most important challenges for these animals. The invasion of the underground niche could have modified the selective regimes of proteins involved in the respiration and transport of O₂ of these rodents, positively selecting mutations of higher affinity for O₂. Here we examine the sequence variation in the beta globin gene in these two lineages, within a robust phylogenetic context. Using different approaches (classical and Bayesian maximum likelihood (PAML/Datamonkey) and alternatives methods (TreeSAAP)) we found at least three sites with evidence of positive selection in underground lineages, especially the basal branch that leads to the Octodontidae family and the branch that leads to the coruro, suggesting some adaptive changes to the underground life. We also found a convergence with another underground rodent, which cannot be identified by the above methods.

Cross-Species Insights Into Genomic Adaptations to Hypoxia

Over millions of years, vertebrate species populated vast environments spanning the globe. Among the most challenging habitats encountered were those with limited availability of oxygen, yet many animal and human populations inhabit and perform life cycle functions and/or daily activities in varying degrees of hypoxia today. Of particular interest are species that inhabit high-altitude niches, which experience chronic hypobaric hypoxia throughout their lives. Physiological and molecular aspects of adaptation to hypoxia have long been the focus of high-altitude populations and, within the past decade, genomic information has become increasingly accessible. These data provide an opportunity to search for common genetic signatures of selection across uniquely informative populations and thereby augment our understanding of the mechanisms underlying adaptations to hypoxia. In this review, we synthesize the available genomic findings across hypoxia-tolerant species to provide a comprehensive view of putatively hypoxia-adaptive genes and pathways. In many cases, adaptive signatures across species converge on the same genetic pathways or on genes themselves [i.e., the hypoxia inducible factor (HIF) pathway). However, specific variants thought to underlie function are distinct between species and populations, and, in most cases, the precise functional role of these genomic differences remains unknown. Efforts to standardize these findings and explore relationships between genotype and phenotype will provide important clues into the evolutionary and mechanistic bases of physiological adaptations to environmental hypoxia.

Comparative transcriptomic analysis of fireflies (Coleoptera: Lampyridae) to explore the molecular adaptations to fresh water

Aquatic insects are well adapted to freshwater environments, but the molecular basis of these adaptations remains largely unknown. Most firefly species (Coleoptera: Lampyridae) are terrestrial, but the larvae of several species are aquatic. Here, larval and adult transcriptomes from Aquatica leii (freshwater) and Lychnuris praetexta (terrestrial) were generated to test whether the genes associated with metabolic efficiency and morphology have undergone adaptive evolution to fresh water. The aquatic fireflies had a significantly lower ratio of nonsynonymous to synonymous substitutions than the terrestrial insects, indicating a genomewide evolutionary constraint in the aquatic fireflies. We identified 341 fast-evolving genes and 116 positively selected genes in the aquatic fireflies. Of these, 76 genes exhibiting both fast evolution and positive selection were primarily involved in ATP production, energy metabolism and the hypoxia response. We identified 7,271 differentially expressed genes (DEGs) in A. leii (adults versus larvae) and 8,309 DEGs in L. praetexta (adults versus larvae). DEGs specific to the aquatic firefly (n = 1,445) were screened via interspecific comparisons (A. leii versus L. praetexta) and were significantly enriched for genes involved in metabolic efficiency (e.g., ATP production, hypoxia, and immune responses) and certain aspects of morphology (e.g., cuticle chitin, tracheal and compound eye morphology). These results indicate that sequence and expression-level changes in genes associated with both metabolic efficiency and morphological attributes related to the freshwater lifestyle contributed to freshwater adaptation in fireflies. This study provides new insights into the molecular mechanisms of aquatic adaptation in insects.

Geometric morphometric analysis of the plateau zokor (Eospalax baileyi) revealed significant effects of environmental factors on skull variations

The plateau zokor (Eospalax baileyi) is employed as an ideal model for examining the relationships between phenotypic and ecological adaptations to the underground conditions in which the skull morphology evolves to adapt to tunnel environment. We evaluated the influence of environmental factors (altitude, temperature, and precipitation) and geographical distance on the variations in skull morphology of a native subterranean rodent plateau zokor population. Thin-plate spline showed that the trend of morphological changes along the CV1 axis was as follows: the two zygomatic arch and the two postorbital processes moved down, the two mastoid processes and the tooth row moved upward, and the tympanic bulla grew longer. The changes along the CV2 axis were as follows: the nasal bone and the tooth row became longer, the distance between the two anterior tips of zygomatic arch lengthened, the infraorbital foramen became smaller, the whole posterior part of the skull became shorter, the zygomatic bone and the two posterior tips of zygomatic arch moved down, and the foramen magnum became bigger. Thus we found significant differences in the skull shape among the seven populations studied. Along with the reduction in the altitude and increase in the mean annual temperature and mean annual precipitation, the nasal bone became shorter, the distance between the two anterior tips of the zygomatic arch became shorter, the whole posterior part of the skull lengthened, the infraorbital foramen became smaller, the two mastoid processes moved upward, and the occipital bone moved down on the dorsal surface of the skull. On the ventral surface of the skull, with an increase in the altitude, mean annual temperature, and mean annual precipitation, the tympanic bulla became shorter, the tooth row moved down, and the foramen magnum became smaller. The morphological changes in the skull were significantly positively correlated with environmental factors. Finally, there was a significant positive correlation between the Procrustes distance matrix of the skull and the geographic distance matrix, which indicates that the evolution of the plateau zokor follows the distance isolation model, but it needs to be further explored from genetic perspectives.

Genome of the webworm Hyphantria cunea unveils genetic adaptations supporting its rapid invasion and spread

BACKGROUND

The fall webworm Hyphantria cunea is an invasive and polyphagous defoliator pest that feeds on nearly any type of deciduous tree worldwide. The silk web of H. cunea aids its aggregating behavior, provides thermal regulation and is regarded as one of causes for its rapid spread. In addition, both chemosensory and detoxification genes are vital for host adaptation in insects.

RESULTS

Here, a high-quality genome of H. cunea was obtained. Silk-web-related genes were identified from the genome, and successful silencing of the silk protein gene HcunFib-H resulted in a significant decrease in silk web shelter production. The CAFE analysis showed that some chemosensory and detoxification gene families, such as CSPs, CCEs, GSTs and UGTs, were expanded. A transcriptome analysis using the newly sequenced H. cunea genome showed that most chemosensory genes were specifically expressed in the antennae, while most detoxification genes were highly expressed during the feeding peak. Moreover, we found that many nutrient-related genes and one detoxification gene, HcunP450 (CYP306A1), were under significant positive selection, suggesting a crucial role of these genes in host adaptation in H. cunea. At the metagenomic level, several microbial communities in H. cunea gut and their metabolic pathways might be beneficial to H. cunea for nutrient metabolism and detoxification, and might also contribute to its host adaptation.

CONCLUSIONS

These findings explain the host and environmental adaptations of H. cunea at the genetic level and provide partial evidence for the cause of its rapid invasion and potential gene targets for innovative pest management strategies.

Suppressed LPS-mediated TLR4 signaling in the plateau zokor (Eospalax baileyi) compared to the bamboo rat (Rhizomys pruinosus) and rat (Rattus norvegicus)

Ecological immunology involves the study of the immune function of wildlife, which is seldom compared with that of model animals. Here, we evaluated and compared the level of the innate immune response in the plateau zokor (Eospalax baileyi), an indigenous underground rodent from the Tibetan Plateau, with that in the bamboo rat (Rhizomys pruinosus) and Sprague-Dawley (SD) rat (Rattus norvegicus). The spleen was observed by ordinary light and transmission electron microscopy, and the spleen index was calculated. After liposaccharide (LPS) challenge, the expression of Toll-like receptor 2 (TLR2), TLR4, and hypoxia-inducible factor 1α (HIF-1α) in the spleen was detected by Western blot analysis and immunofluorescence. The expression of nuclear factor-κB1 (NF-κB1) and mitogen-activated protein kinase 14 (MAPK14) in the spleen was detected by real-time quantitative polymerase chain reaction, and the levels of interleukin 6 (IL-6), tumor necrosis factor-α (TNF-α), and interferon-β (IFN-β) in the spleen were detected by enzyme-linked immunoassay. The spleen index of the plateau zokor was lower than that of the bamboo rat and SD rat. The expression of TLR4, NF-κB1, and MAPK14 and the levels of IL-6 and TNF-α in the spleen of the plateau zokor were lower than those of the bamboo rat and SD rat, while the expression of TLR2 and HIF-1α and the level of IFN-β were higher than those of the bamboo rat and SD rat. We speculate that suppression of the TLR4 signaling pathway in the plateau zokor is an adaptation to hypoxic tunnels that decreases antigenic risk and maintains immune homeostasis. Moreover, the spleen of the plateau zokor is reduced in size, reducing the innate immunity investment in the spleen. We also noted that high levels of HIF-1α in the spleen of the plateau zokor suppressed crosstalk between HIF-1α and TLR4, promoting the innate immune response.

Convergent evolution in human and domesticate adaptation to high-altitude environments

Humans and their domestic animals have lived and thrived in high-altitude environments worldwide for thousands of years. These populations have developed a number of adaptations to survive in a hypoxic environment, and several genomic studies have been conducted to identify the genes that drive these adaptations. Here, we discuss the various adaptations and genetic variants that have been identified as adaptive in human and domestic animal populations and the ways in which convergent evolution has occurred as these populations have adapted to high-altitude environments. We found that human and domesticate populations have adapted to hypoxic environments in similar ways. Specific genes and biological pathways have been involved in high-altitude adaptation for multiple populations, although the specific variants differ between populations. Additionally, we found that the gene EPAS1 is often a target of selection in hypoxic environments and has been involved in multiple adaptive introgression events. High-altitude environments exert strong selective pressures, and human and animal populations have evolved in convergent ways to cope with a chronic lack of oxygen.

This article is part of the theme issue ‘Convergent evolution in the genomics era: new insights and directions'.

Limited Evidence for Parallel Molecular Adaptations Associated with the Subterranean Niche in Mammals: A Comparative Study of Three Superorders

Among mammals, several lineages have independently adapted to a subterranean niche and possess similar phenotypic traits for burrowing (e.g., cylindrical bodies, short limbs, and absent pinnae). Previous research on mole-rats has revealed molecular adaptations for coping with reduced oxygen, elevated carbon dioxide, and the absence of light. In contrast, almost nothing is known regarding molecular adaptations in other subterranean lineages (e.g., true moles and golden moles). Therefore, the extent to which the recurrent phenotypic adaptations of divergent subterranean taxa have arisen via parallel routes of molecular evolution remains untested. To address these issues, we analyzed ∼8,000 loci in 15 representative subterranean taxa of four independent transitions to an underground niche for signatures of positive selection and convergent amino acid substitutions. Complementary analyses were performed in nonsubterranean "control" taxa to assess the biological significance of results. We found comparable numbers of positively selected genes in each of the four subterranean groups; however, correspondence in terms of gene identity between gene sets was low. Furthermore, we did not detect evidence of more convergent amino acids among subterranean species pairs compared with levels found between nonsubterranean controls. Comparisons with nonsubterranean taxa also revealed loci either under positive selection or with convergent substitutions, with similar functional enrichment (e.g., cell adhesion, immune response, and coagulation). Given the limited indication that positive selection and convergence occurred in the same loci, we conclude that selection may have acted on different loci across subterranean mammal lineages to produce similar phenotypes.

Function-related Drivers of Skull Morphometric Variation and Sexual Size Dimorphism in a Subterranean Rodent, Plateau Zokor (Eospalax baileyi)

Sexual dimorphism is prevalent in most living organisms. The difference in size between sexes of a given species is generally known as sexual size dimorphism (SSD). The magnitude of the SSD is determined by Rensch's rule where size dimorphism increases with increasing body size when the male is the larger sex and decreases with increasing average body size when the female is the larger sex. The unique underground environment that zokors (Eospalax baileyi) live under in the severe habitat of the Qinghai‐Tibetan Plateau (QTP) could create SSD selection pressures that may or may not be supported by Rensch's rule, making this scientific question worthy of investigation. In this study, we investigated the individual variation between sexes in body size and SSD of plateau zokors using measurements of 19 morphological traits. We also investigated the evolutionary mechanisms underlying SSD in plateau zokors. Moreover, we applied Rensch's rule to all extant zokor species. Our results showed male‐biased SSD in plateau zokors: The body‐ and head‐related measurements were greater in males than in females. Linear regression analysis between body length, body weight, and carcass weight showed significant relationships with some traits such as skull length, lower incisor length, and tympanic bulla width, which might support our prediction that males have faster growth rates than females. Further, the SSD pattern corroborated the assumption of Rensch's rule in plateau zokors but not in the other zokor species. Our findings suggest that the natural underground habitat and behavioral differences between sexes can generate selection pressures on male traits and contribute to the evolution of SSD in plateau zokors.

Identification and evolution of olfactory genes in the small poplar longhorn beetle Saperda populnea

Saperda populnea is a serious pest of poplar and willow trees in the Palaearctic region, causing extensive damage to forests and the lumber industry. Until recently, there is no safe and effective chemical method to control this pest due to the lack of sufficient knowledge on the molecular basis of its olfactory genes, moreover, the evolutionary history of the olfactory gene family in subfamily Lamiinae is still fully unknown. Our RNA sequencing of the antennae of S. populnea identified 43 odorant binding proteins (OBPs), 15 chemosensory proteins (CSPs), 56 odorant receptors (ORs) and 24 inotropic receptors (IRs) in S. populnea. The RT-PCR results showed several genes were expressed in a sex specific manner, suggesting that these genes might play key role in their olfactory-sensing and sex-related behaviors. Further evolutionary studies were performed on these olfactory genes, overall comparison of the Ka/Ks values of orthologous genes in S. populnea and two other Lamiinae species showed three main conclusions: 1. olfactory genes have evolved more rapidly than the non-olfactory genes in the tested long horn beetles; 2. the IR gene family are under a strong purifying selection; 3. the OBPs of Monochamus alternatus evolved more rapidly than the other two species, which is speculated to be correlated with differentiation of selective pressure in different geographic origins. Detailed evolutionary studies on each olfactory genes showed that several OBPs and ORs are under significantly purifying/relaxed selective pressure, and several positive selection sites were also detected, modeling of SpopOR14 and SpopOBP4/5 showed that most of the positive selection sites were distributed at the N-terminus of SpopOR14, while the positive selection sites in SpopOBP4/5 were located in H-bond donors, results suggest that these sites are more likely to be linked with the selectivity of modeled olfactory genes. The research provided a better understanding of the molecular basis and evolutionary history of the olfactory genes in Lamiinae, through elaborating the mechanism whereby amino structural evolution affects specific variants in OBPs and ORs.

p53 gene cloning and response to hypoxia in the plateau zokor, Myospalax baileyi

The plateau zokor (Myospalax baileyi) is a specialized subterranean rodent that lives on the Qinghai-Tibet Plateau. The species has evolved a series of strategies to adapt to its hypoxic environment and hypercapnia. p53 is a tumour suppressor gene that plays a crucial role in the cellular response to hypoxia by inducing cell cycle arrest, cell apoptosis, DNA damage repair and angiogenesis. To investigate the sequence characteristics of p53 and the response to hypoxia in plateau zokor, we cloned the p53 coding DNA sequence, analysed it, and measured the expression level of p53 at different altitudes in plateau zokor and rats. Our results show that the coding DNA sequence is 1179 bp, consisting of 392 amino acid residues. Compared to human p53, the subterranean rodents have two mutation sites in common with the human hotspots in the DNA-binding domain. Compared to subterranean rodents, plateau zokor have a mutation at residue 309. In addition, subterranean rodents have two convergent sites at residues 78 and 84. The expression levels of p53 in plateau zokor tissues increase significantly from 2260 m to 3300 m, but there was no significant difference in rats at those altitudes. Our results suggest that subterranean rodents have two mutation sites in common with the human hotspots in the DNA-binding domain, the mutation of Gly309Asp is a unique mutation site of plateau zokor p53, and there are two convergent sites enhancing subterranean rodent adaptation to hypoxic conditions. In addition, p53 is sensitive to the oxygen concentration in plateau zokor, and hypoxia upregulates the levels of p53. Generally, plateau zokor use this strategy to adapt to a hypoxic environment.

Comparative transcriptomic analysis of Tibetan Gynaephora to explore the genetic basis of insect adaptation to divergent altitude environments

Adaptation of insects to different altitudes remain largely unknown, especially those endemic to the Tibetan Plateau (TP). Here, we generated the transcriptomes of Gynaephora menyuanensis and G. alpherakii, inhabiting different high altitudes on the TP, and used these and the previously available transcriptomic and genomic sequences from low-altitude insects to explore potential genetic basis for divergent high-altitude adaptation in Gynaephora. An analysis of 5,869 orthologous genes among Gynaephora and other three low-altitude insects uncovered that fast-evolving genes and positively selected genes (PSGs) in the two Gynaephora species were enriched in energy metabolism and hypoxia response categories (e.g. mitochondrion, oxidation-reduction process, and response to oxidative stress). Particularly, mTOR signaling pathway involving hypoxia was enriched by PSGs, indicating this well-known pathway in mammal hypoxia adaptation may be an important signaling system in Gynaephora. Furthermore, some PSGs were associated with response to hypoxia (e.g. cytochrome proteins), cold (e.g. dehydrogenase) and DNA repair (e.g. DNA repair proteins). Interestingly, several insect-specific genes that were associated with exoskeleton and cuticle development (e.g. chitinase and ecdysteroids) had experienced positive selection, suggesting the specific adaptive mechanisms in insects. This study is favourable for understanding the adaptive evolution of Gynaephora and even TP insects to divergent altitudes.

Genetic Adaptation of Schizothoracine Fish to the Phased Uplifting of the Qinghai-Tibetan Plateau

Many species of Schizothoracine, a subfamily of Cyprinidae, are highly endemic to the Qinghai–Tibetan Plateau (QTP). To characterize the adaptive changes associated with the Schizothoracine expansion at high altitudes, we sequenced tissue transcriptomes of two highland and two subhighland Schizothoracines and analyzed gene evolution patterns by comparing with lowland cyprinids. Phylogenetic tree reconstruction and divergence time estimation indicated that the common ancestor of Schizothoracine fish lived ∼32.7 million years ago (MYA), coinciding with the timing of the first phase of QTP uplifting. Both high- and subhigh-Schizothoracines demonstrated elevated dN/dS ratios in the protein-coding genes compared to lowland cyprinids, from which some biological processes implicated in altitude adaptation were commonly identified. On the other hand, the highland and subhighland lineages presented drastically divergent landscapes of positively selected genes (PSGs), enriched with very different gene ontology (GO) profiles, including those in “sensory organ morphogenesis,” “regulation of protein ubiquitination,” “blood circulation,” and “blood vessel development.” These results indicated different selection pressures imposed on the highland and subhighland lineages of the Schizothoracine subfamily, with a higher number of genes in the high-altitude species involved in adaptations such as sensory perception, blood circulation, and protein metabolism. Our study indicated divergent genetic adaptations in the aquatic species facing the phased uplifting of QTP.

Unraveling the message: insights into comparative genomics of the naked mole-rat

Animals have evolved to survive, and even thrive, in different environments. Genetic adaptations may have indirectly created phenotypes that also resulted in a longer lifespan. One example of this phenomenon is the preternaturally long-lived naked mole-rat. This strictly subterranean rodent tolerates hypoxia, hypercapnia, and soil-based toxins. Naked mole-rats also exhibit pronounced resistance to cancer and an attenuated decline of many physiological characteristics that often decline as mammals age. Elucidating mechanisms that give rise to their unique phenotypes will lead to better understanding of subterranean ecophysiology and biology of aging. Comparative genomics could be a useful tool in this regard. Since the publication of a naked mole-rat genome assembly in 2011, analyses of genomic and transcriptomic data have enabled a clearer understanding of mole-rat evolutionary history and suggested molecular pathways (e.g., NRF2-signaling activation and DNA damage repair mechanisms) that may explain the extraordinarily longevity and unique health traits of this species. However, careful scrutiny and re-analysis suggest that some identified features result from incorrect or imprecise annotation and assembly of the naked mole-rat genome: in addition, some of these conclusions (e.g., genes involved in cancer resistance and hairlessness) are rejected when the analysis includes additional, more closely related species. We describe how the combination of better study design, improved genomic sequencing techniques, and new bioinformatic and data analytical tools will improve comparative genomics and ultimately bridge the gap between traditional model and nonmodel organisms.