Comparative and population genomics approaches reveal the basis of adaptation to deserts in a small rodent

Selective sweep analysis of the adaptability of the Yarkand hare (Lepus yarkandensis) to hot arid environments using SLAF-seq

The Yarkand hare (Lepus yarkandensis) inhabits arid desert areas and is endemic to China. It has evolved various adaptations to survive in hot arid environments, including stress responses, the ability to maintain water homeostasis and heat tolerance. Here, we performed a selective sweep analysis to identify the candidate genes for adaptation to hot arid environments in the Yarkand hare. A total of 397 237 single-nucleotide polymorphisms were obtained from 80 Yarkand hares, which inhabit hot arid environments, and 36 Tolai hares (Lepus tolai), which inhabit environments with a mild climate, via specific-locus amplified fragment sequencing. We identified several candidate genes that were associated with the heat stress response (HSPE1), oxidative stress response (SLC23A and GLRX2), immune response (IL1R1 and IRG1), central nervous system development (FGF13, THOC2, FMR1 and MECP2) and regulation of water homeostasis (CDK1) according to fixation index values and θπ ratios in the selective sweep analysis, and six of these genes (GLRX2, IRG1, FGF13, FMR1, MECP2 and CDK1) are newly discovered genes. To the best of our knowledge, this is the first study to identify candidate genes for adaptation to hot arid environments in the Yarkand hare. The results of this study enhance our understanding of the adaptation of the Yarkand hare to hot arid environments and will aid future studies aiming to functionally verify these candidate genes.

Impacts of dietary fat on multi tissue gene expression in the desert-adapted cactus mouse

Understanding the relationship between dietary fat and physiological responses is crucial in species adapted to arid environments where water scarcity is common. In this study, we present a comprehensive exploration of gene expression across five tissues (kidney, liver, lung, gastrointestinal tract, and hypothalamus) and 19 phenotypic measurements, investigating the effects of dietary fat in the desert-adapted cactus mouse ( Peromyscus eremicus ). We show impacts on immune function, circadian gene regulation, and mitochondrial function for mice fed a lower-fat diet compared to mice fed a higher-fat diet. In arid environments with severe water scarcity, even subtle changes in organismal health and water balance can affect physical performance, potentially impacting survival and reproductive success. The study sheds light on the complex interplay between diet, physiological processes, and environmental adaptation, providing valuable insights into the multifaceted impacts of dietary choices on organismal well-being and adaptation strategies in arid habitats.

scalepopgen: Bioinformatic Workflow Resources Implemented in Nextflow for Comprehensive Population Genomic Analyses

Population genomic analyses such as inference of population structure and identifying signatures of selection usually involve the application of a plethora of tools. The installation of tools and their dependencies, data transformation, or series of data preprocessing in a particular order sometimes makes the analyses challenging. While the usage of container-based technologies has significantly resolved the problems associated with the installation of tools and their dependencies, population genomic analyses requiring multistep pipelines or complex data transformation can greatly be facilitated by the application of workflow management systems such as Nextflow and Snakemake. Here, we present scalepopgen, a collection of fully automated workflows that can carry out widely used population genomic analyses on the biallelic single nucleotide polymorphism data stored in either variant calling format files or the plink-generated binary files. scalepopgen is developed in Nextflow and can be run locally or on high-performance computing systems using either Conda, Singularity, or Docker. The automated workflow includes procedures such as (i) filtering of individuals and genotypes; (ii) principal component analysis, admixture with identifying optimal K-values; (iii) running TreeMix analysis with or without bootstrapping and migration edges, followed by identification of an optimal number of migration edges; (iv) implementing single-population and pair-wise population comparison-based procedures to identify genomic signatures of selection. The pipeline uses various open-source tools; additionally, several Python and R scripts are also provided to collect and visualize the results. The tool is freely available at https://github.com/Popgen48/scalepopgen.

When the tap runs dry: The multi-tissue gene expression and physiological responses of water deprived Peromyscus eremicus

The harsh and dry conditions of desert environments have resulted in genomic adaptations, allowing for desert organisms to withstand prolonged drought, extreme temperatures, and limited food resources. Here, we present a comprehensive exploration of gene expression across five tissues (kidney, liver, lung, gastrointestinal tract, and hypothalamus) and 19 phenotypic measurements to explore the whole-organism physiological and genomic response to water deprivation in the desert-adapted cactus mouse (Peromyscus eremicus). The findings encompass the identification of differentially expressed genes and correlative analysis between phenotypes and gene expression patterns across multiple tissues. Specifically, we found robust activation of the vasopressin renin-angiotensin-aldosterone system (RAAS) pathways, whose primary function is to manage water and solute balance. Animals reduce food intake during water deprivation, and upregulation of PCK1 highlights the adaptive response to reduced oral intake via its actions aimed at maintained serum glucose levels. Even with such responses to maintain water balance, hemoconcentration still occurred, prompting a protective downregulation of genes responsible for the production of clotting factors while simultaneously enhancing angiogenesis which is thought to maintains tissue perfusion. In this study, we elucidate the complex mechanisms involved in water balance in the desert-adapted cactus mouse, P. eremicus. By prioritizing a comprehensive analysis of whole-organism physiology and multi-tissue gene expression in a simulated desert environment, we describe the complex and successful response of regulatory processes.

When the tap runs dry: the physiological effects of acute experimental dehydration in Peromyscus eremicus

Desert organisms have evolved physiological, biochemical and genomic mechanisms to survive the extreme aridity of desert environments. Studying desert-adapted species provides a unique opportunity to investigate the survival strategies employed by organisms in some of the harshest habitats on Earth. Two of the primary challenges faced in desert environments are maintaining water balance and thermoregulation. We collected data in a simulated desert environment and a captive colony of cactus mice (Peromyscus eremicus) and used lab-based experiments with real time physiological measurements; energy expenditure, water loss rate and respiratory exchange rate, to characterize the response to water deprivation. Mice without access to water had significantly lower energy expenditures and in turn, reduced water loss compared to mice with access to water after the first 24 h of the experiment. Additionally, we observed significant mass loss that is probably due to dehydration-associated anorexia a response to limit fluid loss by reducing waste and the solute load as well as allowing water reabsorption from the kidneys and gastrointestinal tract. Finally, we observed body temperature correlated with sex, with males without access to water maintaining body temperature when compared with hydrated males, whereas body temperature decreased for females without access to water, suggesting daily metabolic depression in females.

Local adaptation and future climate vulnerability in a wild rodent

As climate change continues, species pushed outside their physiological tolerance limits must adapt or face extinction. When change is rapid, adaptation will largely harness ancestral variation, making the availability and characteristics of that variation of critical importance. Here, we used whole-genome sequencing and genetic-environment association analyses to identify adaptive variation and its significance in the context of future climates in a small Palearctic mammal, the bank vole (Clethrionomys glareolus). We found that peripheral populations of bank vole in Britain are already at the extreme bounds of potential genetic adaptation and may require an influx of adaptive variation in order to respond. Analyses of adaptive loci suggest regional differences in climate variables select for variants that influence patterns of population adaptive resilience, including genes associated with antioxidant defense, and support a pattern of thermal/hypoxic cross-adaptation. Our findings indicate that understanding potential shifts in genomic composition in response to climate change may be key to predicting species' fate under future climates.

Mobilisation of jerboa kidney gene networks during dehydration and opportunistic rehydration

Desert animals have evolved systems that enable them to thrive under dry conditions. Focusing on the kidney, we have investigated the transcriptomic adaptations that enable a desert rodent, the Lesser Egyptian Jerboa (Jaculus jaculus), to withstand water deprivation and opportunistic rehydration. Analysis of the whole kidney transcriptome showed many differentially expressed genes in the Jerboa kidney, 6.4% of genes following dehydration and an even greater number (36.2%) following rehydration compared to control. Genes correlated with the rehydration condition included many ribosomal protein coding genes suggesting a concerted effort to accelerate protein synthesis when water is made available. We identify an increase in TGF-beta signaling antagonists in dehydration (e.g., GREM2). We also describe expression of multiple aquaporin and solute carrier transporters mapped to specific nephron segments. The desert adapted renal transcriptome presented here is a valuable resource to expand our understanding of osmoregulation beyond that derived from model organisms.

When the tap runs dry: The physiological effects of acute experimental dehydration in Peromyscus eremicus

Desert organisms have evolved physiological, biochemical, and genomic mechanisms to survive the extreme aridity of desert environments. Studying desert-adapted species provides a unique opportunity to investigate the survival strategies employed by organisms in some of the harshest habitats on Earth. Two of the primary challenges faced in desert environments are maintaining water balance and thermoregulation. We collected data in a simulated desert environment and a captive colony of cactus mice (Peromyscus eremicus) and used lab-based experiments with real time physiological measurements to characterize the response to water-deprivation. Mice without access to water had significantly lower energy expenditures and in turn, reduced water loss compared to mice with access to water after the first 24 hours of the experiment. Additionally, we observed significant weight loss likely related to dehydration-associated anorexia a response to limit fluid loss by reducing waste and the solute load as well as allowing water reabsorption from the kidneys and gastrointestinal tract. Finally, we observed body temperature correlated with sex, with males without access to water maintaining body temperature when compared to hydrated males while body temperature decreased for females without access to water compared to hydrated, suggesting daily torpor in females.

High-Quality Chromosome-Level Genome Assembly of the Corsac Fox (Vulpes corsac) Reveals Adaptation to Semiarid and Harsh Environments

The Corsac fox (Vulpes corsac) is a species of fox distributed in the arid prairie regions of Central and Northern Asia, with distinct adaptations to dry environments. Here, we applied Oxford-Nanopore sequencing and a chromosome structure capture technique to assemble the first Corsac fox genome, which was then assembled into chromosome fragments. The genome assembly has a total length of 2.2 Gb with a contig N50 of 41.62 Mb and a scaffold N50 of 132.2 Mb over 18 pseudo-chromosomal scaffolds. The genome contained approximately 32.67% of repeat sequences. A total of 20,511 protein-coding genes were predicted, of which 88.9% were functionally annotated. Phylogenetic analyses indicated a close relation to the Red fox (Vulpes vulpes) with an estimated divergence time of ~3.7 million years ago (MYA). We performed separate enrichment analyses of species-unique genes, the expanded and contracted gene families, and positively selected genes. The results suggest an enrichment of pathways related to protein synthesis and response and an evolutionary mechanism by which cells respond to protein denaturation in response to heat stress. The enrichment of pathways related to lipid and glucose metabolism, potentially preventing stress from dehydration, and positive selection of genes related to vision, as well as stress responses in harsh environments, may reveal adaptive evolutionary mechanisms in the Corsac fox under harsh drought conditions. Additional detection of positive selection for genes associated with gustatory receptors may reveal a unique desert diet strategy for the species. This high-quality genome provides a valuable resource for studying mammalian drought adaptation and evolution in the genus Vulpes.

Similar adaptative mechanism but divergent demographic history of four sympatric desert rodents in Eurasian inland

Phenotypes associated with metabolism and water retention are thought to be key to the adaptation of desert species. However, knowledge on the genetic changes and selective regimes on the similar and divergent ways to desert adaptation in sympatric and phylogenetically close desert organisms remains limited. Here, we generate a chromosome level genome assembly for Northern three-toed jerboa (Dipus sagitta) and three other high-quality genome assemblies for Siberian jerboa (Orientallactaga sibirica), Midday jird (Meriones meridianus), and Desert hamster (Phodopus roborovskii). Genomic analyses unveil that desert adaptation of the four species mainly result from similar metabolic pathways, such as arachidonic acid metabolism, thermogenesis, oxidative phosphorylation, insulin related pathway, DNA repair and protein synthesis and degradation. However, the specific evolved genes in the same adaptative molecular pathway often differ in the four species. We also reveal similar niche selection but different demographic histories and sensitivity to climate changes, which may be related to the diversified genomic adaptative features. In addition, our study suggests that nocturnal rodents have evolved some specific adaptative mechanism to desert environments compared to large desert animals. Our genomic resources will provide an important foundation for further research on desert genetic adaptations.

Shared Patterns of Gene Expression and Protein Evolution Associated with Adaptation to Desert Environments in Rodents

Desert specialization has arisen multiple times across rodents and is often associated with a suite of convergent phenotypes, including modification of the kidneys to mitigate water loss. However, the extent to which phenotypic convergence in desert rodents is mirrored at the molecular level is unknown. Here, we sequenced kidney mRNA and assembled transcriptomes for three pairs of rodent species to search for shared differences in gene expression and amino acid sequence associated with adaptation to deserts. We conducted phylogenetically independent comparisons between a desert specialist and a non-desert relative in three families representing ∼70 million years of evolution. Overall, patterns of gene expression faithfully recapitulated the phylogeny of these six taxa providing a strong evolutionary signal in levels of mRNA abundance. We also found that 8.6% of all genes showed shared patterns of expression divergence between desert and non-desert taxa, much of which likely reflects convergent evolution, and representing more than expected by chance under a model of independent gene evolution. In addition to these shared changes, we observed many species-pair-specific changes in gene expression indicating that instances of adaptation to deserts include a combination of unique and shared changes. Patterns of protein evolution revealed a small number of genes showing evidence of positive selection, the majority of which did not show shared changes in gene expression. Overall, our results suggest that convergent changes in gene regulation play an important role in the complex trait of desert adaptation in rodents.

An eco-evo-devo genetic network model of stress response

The capacity of plants to resist abiotic stresses is of great importance to agricultural, ecological and environmental sustainability, but little is known about its genetic underpinnings. Existing genetic tools can identify individual genetic variants mediating biochemical, physiological, and cellular defenses, but fail to chart an overall genetic atlas behind stress resistance. We view stress response as an eco-evo-devo process by which plants adaptively respond to stress through complex interactions of developmental canalization, phenotypic plasticity, and phenotypic integration. As such, we define and quantify stress response as the developmental change of adaptive traits from stress-free to stress-exposed environments. We integrate composite functional mapping and evolutionary game theory to reconstruct omnigenic, information-flow interaction networks for stress response. Using desert-adapted Euphrates poplar as an example, we infer salt resistance-related genome-wide interactome networks and trace the roadmap of how each SNP acts and interacts with any other possible SNPs to mediate salt resistance. We characterize the previously unknown regulatory mechanisms driving trait variation; i.e. the significance of a SNP may be due to the promotion of positive regulators, whereas the insignificance of a SNP may result from the inhibition of negative regulators. The regulator-regulatee interactions detected are not only experimentally validated by two complementary experiments, but also biologically interpreted by their encoded protein-protein interactions. Our eco-evo-devo model of genetic interactome networks provides an approach to interrogate the genetic architecture of stress response and informs precise gene editing for improving plants' capacity to live in stress environments.

Evolutionary trade-offs between heat and cold tolerance limit responses to fluctuating climates

The evolutionary potential of species to cope with short-term temperature fluctuations during reproduction is critical to predicting responses to future climate change. Despite this, vertebrate research has focused on reproduction under high or low temperatures in relatively stable temperate climates. Here, we characterize the genetic basis of reproductive thermal tolerance to temperature fluctuations in the ostrich, which lives in variable environments in tropical and subtropical Africa. Both heat and cold tolerance were under selection and heritable, indicating the potential for evolutionary responses to mean temperature change. However, we found evidence for a negative, genetic correlation between heat and cold tolerance that should limit the potential for adaptation to fluctuating temperatures. Genetic constraints between heat and cold tolerance appear a crucial, yet underappreciated, factor influencing responses to climate change.

Genomic landscape of a relict fir-associated fungus reveals rapid convergent adaptation towards endophytism

Comparative and pan-genomic analyses of the endophytic fungus Pezicula neosporulosa (Helotiales, Ascomycota) from needles of the relict fir, Abies beshanzuensis, showed expansions of carbohydrate metabolism and secondary metabolite biosynthetic genes characteristic for unrelated plant-beneficial helotialean, such as dark septate endophytes and ericoid mycorrhizal fungi. The current species within the relatively young Pliocene genus Pezicula are predominantly saprotrophic, while P. neosporulosa lacks such features. To understand the genomic background of this putatively convergent evolution, we performed population analyses of 77 P. neosporulosa isolates. This revealed a mosaic structure of a dozen non-recombining and highly genetically polymorphic subpopulations with a unique mating system structure. We found that one idiomorph of a probably duplicated mat1-2 gene was found in putatively heterothallic isolates, while the other co-occurred with mat1-1 locus suggesting homothallic reproduction for these strains. Moreover, 24 and 81 genes implicated in plant cell-wall degradation and secondary metabolite biosynthesis, respectively, showed signatures of the balancing selection. These findings highlight the evolutionary pattern of the two gene families for allowing the fungus a rapid adaptation towards endophytism and facilitating diverse symbiotic interactions.

Chromosome size affects sequence divergence between species through the interplay of recombination and selection

The structure of the genome shapes the distribution of genetic diversity and sequence divergence. To investigate how the relationship between chromosome size and recombination rate affects sequence divergence between species, we combined empirical analyses and evolutionary simulations. We estimated pairwise sequence divergence among 15 species from three different mammalian clades-Peromyscus rodents, Mus mice, and great apes-from chromosome-level genome assemblies. We found a strong significant negative correlation between chromosome size and sequence divergence in all species comparisons within the Peromyscus and great apes clades but not the Mus clade, suggesting that the dramatic chromosomal rearrangements among Mus species may have masked the ancestral genomic landscape of divergence in many comparisons. Our evolutionary simulations showed that the main factor determining differences in divergence among chromosomes of different sizes is the interplay of recombination rate and selection, with greater variation in larger populations than in smaller ones. In ancestral populations, shorter chromosomes harbor greater nucleotide diversity. As ancestral populations diverge, diversity present at the onset of the split contributes to greater sequence divergence in shorter chromosomes among daughter species. The combination of empirical data and evolutionary simulations revealed that chromosomal rearrangements, demography, and divergence times may also affect the relationship between chromosome size and divergence, thus deepening our understanding of the role of genome structure in the evolution of species divergence.

Multiple types of genomic variation contribute to adaptive traits in the mustelid subfamily Guloninae

Species of the mustelid subfamily Guloninae inhabit diverse habitats on multiple continents, and occupy a variety of ecological niches. They differ in feeding ecologies, reproductive strategies and morphological adaptations. To identify candidate loci associated with adaptations to their respective environments, we generated a de novo assembly of the tayra (Eira barbara), the earliest diverging species in the subfamily, and compared this with the genomes available for the wolverine (Gulo gulo) and the sable (Martes zibellina). Our comparative genomic analyses included searching for signs of positive selection, examining changes in gene family sizes, as well as searching for species-specific structural variants (SVs). Among candidate loci associated with phenotypic traits, we observed many related to diet, body condition and reproduction. For example, for the tayra, which has an atypical gulonine reproductive strategy of aseasonal breeding, we observe species-specific changes in many pregnancy-related genes. For the wolverine, a circumpolar hypercarnivore that must cope with seasonal food scarcity, we observed many changes in genes associated with diet and body condition. All types of genomic variation examined (single nucleotide polymorphisms, gene family expansions, structural variants) contributed substantially to the identification of candidate loci. This strongly argues for consideration of variation other than single nucleotide polymorphisms in comparative genomics studies aiming to identify loci of adaptive significance.

Shrinkage of body size in rodents as a strategy of populations under anthropogenic conditions (results of 50 years of study of rodent populations)

Based on results of 50 years of continuous observations of populations of the pine vole (Microtus subterraneus de Selys-Longchamps, 1836) and yellow-necked wood mouse (Sylvaemus flavicollis Melchior, 1834), the phenomenon of reducing body size of individuals was revealed, particularly a significant reduction in body dimensions and body mass. The research was conducted at the Kaniv Nature Reserve covering various stages of the local ecosystem's existence under anthropogenic impact. Four cycles of density dynamics of two populations are chosen for comparison of their characteristics. The first three cycles correspond to qualitatively different periods in the existence of the protected ecosystem and populations of the studied species, whereas the last cycle reflects the current situation. It has been established that the tendency to decrease body dimensions is also characteristic for other rodent species, its manifestations are observed for 35 years, and the scale gradually increases. By fatness indexes, individuals of the pine vole on average lose 23.3% and yellow-necked wood mice lose 16.7% of the former values. It was shown that the level of reduce in exterior parameters is always greater in females than in males, and in the group of breeding adults it can reach 32%. It has been found that the phenomenon also applies to juveniles, whose fatness decreases by an average of 21%. The phenomenon occurs against the background of violations of various aspects of population dynamics, which allows it to be associated with anthropogenic changes in the environment. It is suggested that shrinking can be realized by various mechanisms. First, as a result of mortality, the largest individuals and reproductive females with their greatest energy needs fall out of the population, and, second, the growth and weight gain of young animals is slower. As a result, the decrease in external parameters of individuals reduces their specific energy needs and allows them to better survive under adverse conditions. From these points of view, reducing body dimensions can be considered as a specific population strategy to maintain ecological balance.

Disentangling environmental drivers of circadian metabolism in desert-adapted mice

Metabolism is a complex phenotype shaped by natural environmental rhythms, as well as behavioral, morphological and physiological adaptations. Metabolism has been historically studied under constant environmental conditions, but new methods of continuous metabolic phenotyping now offer a window into organismal responses to dynamic environments, and enable identification of abiotic controls and the timing of physiological responses relative to environmental change. We used indirect calorimetry to characterize metabolic phenotypes of the desert-adapted cactus mouse (Peromyscus eremicus) in response to variable environmental conditions that mimic their native environment versus those recorded under constant warm and constant cool conditions, with a constant photoperiod and full access to resources. We found significant sexual dimorphism, with males being more prone to dehydration than females. Under circadian environmental variation, most metabolic shifts occurred prior to physical environmental change and the timing was disrupted under both constant treatments. The ratio of CO2 produced to O2 consumed (the respiratory quotient) reached greater than 1.0 only during the light phase under diurnally variable conditions, a pattern that strongly suggests that lipogenesis contributes to the production of energy and endogenous water. Our results are consistent with historical descriptions of circadian torpor in this species (torpid by day, active by night), but reject the hypothesis that torpor is initiated by food restriction or negative water balance.

Sex-biased gene expression and recent sex chromosome turnover

Cichlids are well known for their propensity to radiate generating arrays of morphologically and ecologically diverse species in short evolutionary time. Following this rapid evolutionary pace, cichlids show high rates of sex chromosome turnover. We here studied the evolution of sex-biased gene (SBG) expression in 14 recently diverged taxa of the Lake Tanganyika Tropheini cichlids, which show different XY sex chromosomes. Across species, sex chromosome sequence divergence predates divergence in expression between the sexes. Only one sex chromosome, the oldest, showed signs of demasculinization in gene expression and potentially contribution to the resolution of sexual conflict. SBGs in general showed high rates of turnovers and evolved mostly under drift. Sexual selection did not shape the rapid evolutionary changes of SBGs. Male-biased genes evolved faster than female-biased genes, which seem to be under more phylogenetic constraint. We found a relationship between the degree of sex bias and sequence evolution driven by sequence differences among the sexes. Consistent with other species, strong sex bias towards sex-limited expression contributes to resolving sexual conflict in cichlids. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Local adaptation of bitter taste and ecological speciation in a wild mammal

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.

Multiomic analysis of the Arabian camel (Camelus dromedarius) kidney reveals a role for cholesterol in water conservation

The Arabian camel (Camelus dromedarius) is the most important livestock animal in arid and semi-arid regions and provides basic necessities to millions of people. In the current context of climate change, there is renewed interest in the mechanisms that enable camelids to survive in arid conditions. Recent investigations described genomic signatures revealing evolutionary adaptations to desert environments. We now present a comprehensive catalogue of the transcriptomes and proteomes of the dromedary kidney and describe how gene expression is modulated as a consequence of chronic dehydration and acute rehydration. Our analyses suggested an enrichment of the cholesterol biosynthetic process and an overrepresentation of categories related to ion transport. Thus, we further validated differentially expressed genes with known roles in water conservation which are affected by changes in cholesterol levels. Our datasets suggest that suppression of cholesterol biosynthesis may facilitate water retention in the kidney by indirectly facilitating the AQP2-mediated water reabsorption.

Chromosome-Level Assembly of the Atlantic Silverside Genome Reveals Extreme Levels of Sequence Diversity and Structural Genetic Variation

The levels and distribution of standing genetic variation in a genome can provide a wealth of insights about the adaptive potential, demographic history, and genome structure of a population or species. As structural variants are increasingly associated with traits important for adaptation and speciation, investigating both sequence and structural variation is essential for wholly tapping this potential. Using a combination of shotgun sequencing, 10x Genomics linked reads and proximity-ligation data (Chicago and Hi-C), we produced and annotated a chromosome-level genome assembly for the Atlantic silverside (Menidia menidia)-an established ecological model for studying the phenotypic effects of natural and artificial selection-and examined patterns of genomic variation across two individuals sampled from different populations with divergent local adaptations. Levels of diversity varied substantially across each chromosome, consistently being highly elevated near the ends (presumably near telomeric regions) and dipping to near zero around putative centromeres. Overall, our estimate of the genome-wide average heterozygosity in the Atlantic silverside is among the highest reported for a fish, or any vertebrate (1.32-1.76% depending on inference method and sample). Furthermore, we also found extreme levels of structural variation, affecting ∼23% of the total genome sequence, including multiple large inversions (> 1 Mb and up to 12.6 Mb) associated with previously identified haploblocks showing strong differentiation between locally adapted populations. These extreme levels of standing genetic variation are likely associated with large effective population sizes and may help explain the remarkable adaptive divergence among populations of the Atlantic silverside.

Gene expression plasticity and desert adaptation in house mice

Understanding how organisms adapt to new environments is a key problem in evolution, yet it remains unclear whether phenotypic plasticity generally facilitates or hinders this process. Here we studied evolved and plastic responses to water-stress in lab-born descendants of wild house mice (Mus musculus domesticus) collected from desert and non-desert environments and measured gene expression and organismal phenotypes under control and water-stressed conditions. After many generations in the lab, desert mice consumed significantly less water than mice from other localities, indicating that this difference has a genetic basis. Under water-stress, desert mice maintained more weight than non-desert mice, and exhibited differences in blood chemistry related to osmoregulatory function. Gene expression in the kidney revealed evolved differences between mice from different environments as well as plastic responses between hydrated and dehydrated mice. Desert mice showed reduced expression plasticity under water-stress compared to non-desert mice. Importantly, non-desert mice under water-stress generally showed shifts toward desert-like expression, consistent with adaptive plasticity. Finally, we identify several co-expression modules linked to phenotypes of interest. These findings provide evidence for local adaptation after a recent invasion and suggest that adaptive plasticity may have facilitated colonization of the desert environment.

Integrating genetics, physiology and morphology to study desert adaptation in a lizard species

Integration of multiple approaches is key to understand the evolutionary processes of local adaptation and speciation. Reptiles have successfully colonized desert environments, that is, extreme and arid conditions that constitute a strong selective pressure on organisms. Here, we studied genomic, physiological and morphological variations of the lizard Liolaemus fuscus to detect adaptations to the Atacama Desert. By comparing populations of L. fuscus inhabiting the Atacama Desert with populations from the Mediterranean forests from central Chile, we aimed at characterizing features related to desert adaptation. We combined ddRAD sequencing with physiological (evaporative water loss, metabolic rate and selected temperature) and morphological (linear and geometric morphometrics) measurements. We integrated the genomic and phenotypic data using redundancy analyses. Results showed strong genetic divergence, along with a high number of fixed loci between desert and forest populations. Analyses detected 110 fixed and 30 outlier loci located within genes, from which 43 were in coding regions, and 12 presented non-synonymous mutations. The candidate genes were associated with cellular membrane and development. Desert lizards presented lower evaporative water loss than those from the forest. Morphological data showed that desert lizards had smaller body size, different allometry, larger eyeballs and more dorsoventrally compressed heads. Our results suggest incipient speciation between desert and forest populations. The adaptive signal must be cautiously interpreted since genetic drift could also contribute to the divergence pattern. Nonetheless, we propose water and resource availability, and changes in habitat structure, as the most relevant challenges for desert reptiles. This study provides insights of the mechanisms that allow speciation as well as desert adaptation in reptiles at multiple levels, and highlights the benefit of integrating independent evidence.

Limited Evidence for Parallel Evolution Among Desert-Adapted Peromyscus Deer Mice

Warming climate and increasing desertification urge the identification of genes involved in heat and dehydration tolerance to better inform and target biodiversity conservation efforts. Comparisons among extant desert-adapted species can highlight parallel or convergent patterns of genome evolution through the identification of shared signatures of selection. We generate a chromosome-level genome assembly for the canyon mouse (Peromyscus crinitus) and test for a signature of parallel evolution by comparing signatures of selective sweeps across population-level genomic resequencing data from another congeneric desert specialist (Peromyscus eremicus) and a widely distributed habitat generalist (Peromyscus maniculatus), that may be locally adapted to arid conditions. We identify few shared candidate loci involved in desert adaptation and do not find support for a shared pattern of parallel evolution. Instead, we hypothesize divergent molecular mechanisms of desert adaptation among deer mice, potentially tied to species-specific historical demography, which may limit or enhance adaptation. We identify a number of candidate loci experiencing selective sweeps in the P. crinitus genome that are implicated in osmoregulation (Trypsin, Prostasin) and metabolic tuning (Kallikrein, eIF2-alpha kinase GCN2, APPL1/2), which may be important for accommodating hot and dry environmental conditions.

Molecular Evolution of Ecological Specialisation: Genomic Insights from the Diversification of Murine Rodents

Adaptive radiations are characterized by the diversification and ecological differentiation of species, and replicated cases of this process provide natural experiments for understanding the repeatability and pace of molecular evolution. During adaptive radiation, genes related to ecological specialization may be subject to recurrent positive directional selection. However, it is not clear to what extent patterns of lineage-specific ecological specialization (including phenotypic convergence) are correlated with shared signatures of molecular evolution. To test this, we sequenced whole exomes from a phylogenetically dispersed sample of 38 murine rodent species, a group characterized by multiple, nested adaptive radiations comprising extensive ecological and phenotypic diversity. We found that genes associated with immunity, reproduction, diet, digestion, and taste have been subject to pervasive positive selection during the diversification of murine rodents. We also found a significant correlation between genome-wide positive selection and dietary specialization, with a higher proportion of positively selected codon sites in derived dietary forms (i.e., carnivores and herbivores) than in ancestral forms (i.e., omnivores). Despite striking convergent evolution of skull morphology and dentition in two distantly related worm-eating specialists, we did not detect more genes with shared signatures of positive or relaxed selection than in a nonconvergent species comparison. Although a small number of the genes we detected can be incidentally linked to craniofacial morphology or diet, protein-coding regions are unlikely to be the primary genetic basis of this complex convergent phenotype. Our results suggest a link between positive selection and derived ecological phenotypes, and highlight specific genes and general functional categories that may have played an integral role in the extensive and rapid diversification of murine rodents.

Life in Deserts: The Genetic Basis of Mammalian Desert Adaptation

Deserts are among the harshest environments on Earth. The multiple ages of different deserts and their global distribution provide a unique opportunity to study repeated adaptation at different timescales. Here, we summarize recent genomic research on the genetic mechanisms underlying desert adaptations in mammals. Several studies on different desert mammals show large overlap in functional classes of genes and pathways, consistent with the complexity and variety of phenotypes associated with desert adaptation to water and food scarcity and extreme temperatures. However, studies of desert adaptation are also challenged by a lack of accurate genotype-phenotype-environment maps. We encourage development of systems that facilitate functional analyses, but also acknowledge the need for more studies on a wider variety of desert mammals.

Compromised Function of the Pancreatic Transcription Factor PDX1 in a Lineage of Desert Rodents

Gerbils are a subfamily of rodents living in arid regions of Asia and Africa. Recent studies have shown that several gerbil species have unusual amino acid changes in the PDX1 protein, a homeodomain transcription factor essential for pancreatic development and β-cell function. These changes were linked to strong GC-bias in the genome that may be caused by GC-biased gene conversion, and it has been hypothesized that this caused accumulation of deleterious changes. Here we use two approaches to examine if the unusual changes are adaptive or deleterious. First, we compare PDX1 protein sequences between 38 rodents to test for association with habitat. We show the PDX1 homeodomain is almost totally conserved in rodents, apart from gerbils, regardless of habitat. Second, we use ectopic gene overexpression and gene editing in cell culture to compare functional properties of PDX1 proteins. We show that the divergent gerbil PDX1 protein inefficiently binds an insulin gene promoter and ineffectively regulates insulin expression in response to high glucose in rat cells. The protein has, however, retained the ability to regulate some other β-cell genes. We suggest that during the evolution of gerbils, the selection-blind process of biased gene conversion pushed fixation of mutations adversely affecting function of a normally conserved homeodomain protein. We argue these changes were not entirely adaptive and may be associated with metabolic disorders in gerbil species on high carbohydrate diets. This unusual pattern of molecular evolution could have had a constraining effect on habitat and diet choice in the gerbil lineage.

Comparative genomic analysis of sifakas (Propithecus) reveals selection for folivory and high heterozygosity despite endangered status

Sifakas (genus Propithecus) are critically endangered, large-bodied diurnal lemurs that eat leaf-based diets and show corresponding anatomical and microbial adaptations to folivory. We report on the genome assembly of Coquerel's sifaka (P. coquereli) and the resequenced genomes of Verreaux's (P. verreauxi), the golden-crowned (P. tattersalli), and the diademed (P. diadema) sifakas. We find high heterozygosity in all sifakas compared with other primates and endangered mammals. Demographic reconstructions nevertheless suggest declines in effective population size beginning before human arrival on Madagascar. Comparative genomic analyses indicate pervasive accelerated evolution in the ancestral sifaka lineage affecting genes in several complementary pathways relevant to folivory, including nutrient absorption and xenobiotic and fatty acid metabolism. Sifakas show convergent evolution at the level of the pathway, gene family, gene, and amino acid substitution with other folivores. Although sifakas have relatively generalized diets, the physiological challenges of habitual folivory likely led to strong selection.

A linked-read approach to museomics: Higher quality de novo genome assemblies from degraded tissues

High-throughput sequencing technologies are a proposed solution for accessing the molecular data in historical specimens. However, degraded DNA combined with the computational demands of short-read assemblies has posed significant laboratory and bioinformatics challenges for de novo genome assembly. Linked-read or "synthetic long-read" sequencing technologies, such as 10× Genomics, may provide a cost-effective alternative solution to assemble higher quality de novo genomes from degraded tissue samples. Here, we compare assembly quality (e.g., genome contiguity and completeness, presence of orthogroups) between four new deer mouse (Peromyscus spp.) genomes assembled using linked-read technology and four published genomes assembled from a single shotgun library. At a similar price-point, these approaches produce vastly different assemblies, with linked-read assemblies having overall higher contiguity and completeness, measured by larger N50 values and greater number of genes assembled, respectively. As a proof-of-concept, we used annotated genes from the four Peromyscus linked-read assemblies and eight additional rodent taxa to generate a phylogeny, which reconstructed the expected relationships among species with 100% support. Although not without caveats, our results suggest that linked-read sequencing approaches are a viable option to build de novo genomes from degraded tissues, which may prove particularly valuable for taxa that are extinct, rare or difficult to collect.

A Roadmap for Understanding the Evolutionary Significance of Structural Genomic Variation

Structural genomic variants (SVs) are ubiquitous and play a major role in adaptation and speciation. Yet, comparative and population genomics have focused predominantly on gene duplications and large-effect inversions. The lack of a common framework for studying all SVs is hampering progress towards a more systematic assessment of their evolutionary significance. Here we (i) review how different types of SVs affect ecological and evolutionary processes; (ii) suggest unifying definitions and recommendations for future studies; and (iii) provide a roadmap for the integration of SVs in ecoevolutionary studies. In doing so, we lay the foundation for population genomics, theoretical, and experimental approaches to understand how the full spectrum of SVs impacts ecological and evolutionary processes.

Comparative and population genomics approaches reveal the basis of adaptation to deserts in a small rodent

Organisms that live in deserts offer the opportunity to investigate how species adapt to environmental conditions that are lethal to most plants and animals. In the hot deserts of North America, high temperatures and lack of water are conspicuous challenges for organisms living there. The cactus mouse (Peromyscus eremicus) displays several adaptations to these conditions, including low metabolic rate, heat tolerance, and the ability to maintain homeostasis under extreme dehydration. To investigate the genomic basis of desert adaptation in cactus mice, we built a chromosome‐level genome assembly and resequenced 26 additional cactus mouse genomes from two locations in southern California (USA). Using these data, we integrated comparative, population, and functional genomic approaches. We identified 16 gene families exhibiting significant contractions or expansions in the cactus mouse compared to 17 other Myodontine rodent genomes, and found 232 sites across the genome associated with selective sweeps. Functional annotations of candidate gene families and selective sweeps revealed a pervasive signature of selection at genes involved in the synthesis and degradation of proteins, consistent with the evolution of cellular mechanisms to cope with protein denaturation caused by thermal and hyperosmotic stress. Other strong candidate genes included receptors for bitter taste, suggesting a dietary shift towards chemically defended desert plants and insects, and a growth factor involved in lipid metabolism, potentially involved in prevention of dehydration. Understanding how species adapted to deserts will provide an important foundation for predicting future evolutionary responses to increasing temperatures, droughts and desertification in the cactus mouse and other species.