Short Evolutionary Timescales Research Articles

The Rapid Evolution of De Novo Proteins in Structure and Complex.

Recent studies in the rice genome-wide have established that de novo genes, evolving from non-coding sequences, enhance protein diversity through a stepwise process. However, the pattern and rate of their evolution in protein structure over time remain unclear. Here, we addressed these issues within a surprisingly short evolutionary timescale (<1 million years for 97% of Oryza de novo genes) with comparative approaches to gene duplicates. We found that de novo genes evolve faster than gene duplicates in the intrinsic disordered regions (IDRs, such as random coils), secondary structure elements (such as α-helix and β-strand), hydrophobicity, and molecular recognition features (MoRFs). In de novo proteins, specifically, we observed an 8-14% decay in random coils and IDR lengths and a 2.3-6.5% increase in structured elements, hydrophobicity, and MoRFs, per million years on average. These patterns of structural evolution align with changes in amino acid composition over time as well. We also revealed higher positive charges but smaller molecular weights for de novo proteins than duplicates. Tertiary structure predictions showed that most de novo proteins, though not typically well-folded on their own, readily form low-energy and compact complexes with other proteins facilitated by extensive residue contacts and conformational flexibility, suggesting a faster-binding scenario in de novo proteins to promote interaction. These analyses illuminate a rapid evolution of protein structure in de novo genes in rice genomes, originating from non-coding sequences, highlighting their quick transformation into active, protein complex-forming components within a remarkably short evolutionary timeframe.

Rapid evolutionary change in trait correlations of single proteins

Many organismal traits are genetically determined and covary in evolving populations. The resulting trait correlations can either help or hinder evolvability – the ability to bring forth new and adaptive phenotypes. The evolution of evolvability requires that trait correlations themselves must be able to evolve, but we know little about this ability. To learn more about it, we here study two evolvable systems, a yellow fluorescent protein and the antibiotic resistance protein VIM-2 metallo beta-lactamase. We consider two traits in the fluorescent protein, namely the ability to emit yellow and green light, and three traits in our enzyme, namely the resistance against ampicillin, cefotaxime, and meropenem. We show that correlations between these traits can evolve rapidly through both mutation and selection on short evolutionary time scales. In addition, we show that these correlations are driven by a protein’s ability to fold, because single mutations that alter foldability can dramatically change trait correlations. Since foldability is important for most proteins and their traits, mutations affecting protein folding may alter trait correlations mediated by many other proteins. Thus, mutations that affect protein foldability may also help shape the correlations of complex traits that are affected by hundreds of proteins.

Evolution of Transcript Abundance is Influenced by Indels in Protein Low Complexity Regions.

Protein Protein low complexity regions (LCRs) are compositionally biased amino acid sequences, many of which have significant evolutionary impacts on the proteins which contain them. They are mutationally unstable experiencing higher rates of indels and substitutions than higher complexity regions. LCRs also impact the expression of their proteins, likely through multiple effects along the path from gene transcription, through translation, and eventual protein degradation. It has been observed that proteins which contain LCRs are associated with elevated transcript abundance (TAb), despite having lower protein abundance. We have gathered and integrated human data to investigate the co-evolution of TAb and LCRs through ancestral reconstructions and model inference using an approximate Bayesian calculation based method. We observe that on short evolutionary timescales TAb evolution is significantly impacted by changes in LCR length, with insertions driving TAb down. But in contrast, the observed data is best explained by indel rates in LCRs which are unaffected by shifts in TAb. Our work demonstrates a coupling between LCR and TAb evolution, and the utility of incorporating multiple responses into evolutionary analyses.

Pervasive mimicry in flight behavior among aposematic butterflies

Flight was a key innovation in the adaptive radiation of insects. However, it is a complex trait influenced by a large number of interacting biotic and abiotic factors, making it difficult to unravel the evolutionary drivers. We investigate flight patterns in neotropical heliconiine butterflies, well known for mimicry of their aposematic wing color patterns. We quantify the flight patterns (wing beat frequency and wing angles) of 351 individuals representing 29 heliconiine and 9 ithomiine species belonging to ten color pattern mimicry groupings. For wing beat frequency and up wing angles, we show that heliconiine species group by color pattern mimicry affiliation. Convergence of down wing angles to mimicry groupings is less pronounced, indicating that distinct components of flight are under different selection pressures and constraints. The flight characteristics of the Tiger mimicry group are particularly divergent due to convergence with distantly related ithomiine species. Predator-driven selection for mimicry also explained variation in flight among subspecies, indicating that this convergence can occur over relatively short evolutionary timescales. Our results suggest that the flight convergence is driven by aposematic signaling rather than shared habitat between comimics. We demonstrate that behavioral mimicry can occur between lineages that have separated over evolutionary timescales ranging from <0.5 to 70 My.

Tissue-specific and cis-regulatory changes underlie parallel, adaptive gene expression evolution in house mice.

Changes in gene regulation have long been appreciated as a driving force of adaptive evolution, however the relative contributions of cis- and trans-acting changes to gene regulation over short evolutionary timescales remain unclear. Instances of recent, parallel phenotypic evolution provide an opportunity to assess whether parallel patterns are seen at the level of gene expression, and to assess the relative contribution of cis- and trans- changes to gene regulation in the early stages of divergence. Here, we studied gene expression in liver and brown adipose tissue in two wild-derived strains of house mice that independently adapted to cold, northern environments, and we compared them to a strain of house mice from a warm, tropical environment. To investigate gene regulatory evolution, we studied expression in parents and allele-specific expression in F1 hybrids of crosses between warm-adapted and cold-adapted strains. First, we found that the different cold-adapted mice showed both unique and shared changes in expression, but that the proportion of shared changes (i.e. parallelism) was greater than expected by chance. Second, we discovered that expression evolution occurred largely at tissue-specific and cis-regulated genes, and that these genes were over-represented in parallel cases of evolution. Finally, we integrated the expression data with scans for selection in natural populations and found substantial parallelism in the two northern populations for genes under selection. Furthermore, selection outliers were associated with cis-regulated genes more than expected by chance; cis-regulated genes under selection influenced phenotypes such as body size, immune functioning, and activity level. These results demonstrate that parallel patterns of gene expression in mice that have independently adapted to cold environments are driven largely by tissue-specific and cis-regulatory changes, providing insight into the mechanisms of adaptive gene regulatory evolution at the earliest stages of divergence.

Bacteria use a catabolic patchwork pathway of apparently recent origin for degradation of the synthetic buffer compound TRIS.

The synthetic buffer compound TRIS (2-amino-2-(hydroxymethyl)propane-1,3-diol) is used in countless applications, and no detailed information on its degradation has been published so far. Herein, we describe the discovery of a complete bacterial degradation pathway for TRIS. By serendipity, a Pseudomonas strain was isolated from sewage sludge that was able to grow with TRIS as only carbon and nitrogen source. Genome and transcriptome analyses revealed two adjacent gene clusters embedded in a mobile genetic element on a conjugative plasmid to be involved in TRIS degradation. Heterologous gene expression revealed cluster I to encode a TRIS uptake protein, a TRIS alcohol dehydrogenase, and a TRIS aldehyde dehydrogenase, catalyzing the oxidation of TRIS into 2-hydroxymethylserine. Gene cluster II encodes a methylserine hydroxymethyltransferase (mSHMT) and a d-serine dehydratase that plausibly catalyze the conversion of 2-hydroxymethylserine into pyruvate. Conjugational plasmid transfer into Pseudomonas putida KT2440 enabled this strain to grow with TRIS and with 2-hydromethylserine, demonstrating that the complete TRIS degradation pathway can be transmitted by horizontal gene transfer. Subsequent enrichments from wastewater purification systems led to the isolation of further TRIS-degrading bacteria from the Pseudomonas and Shinella genera carrying highly similar TRIS degradation gene clusters. Our data indicate that TRIS degradation evolved recently via gene recruitment and enzyme adaptation from multiple independent metabolic pathways, and database searches suggest that the TRIS degradation pathway is now globally distributed. Overall, our study illustrates how engineered environments can enhance the emergence of new microbial metabolic pathways in short evolutionary time scales.

Gene expression and alternative splicing contribute to adaptive divergence of ecotypes.

Regulation of gene expression is a critical link between genotype and phenotype explaining substantial heritable variation within species. However, we are only beginning to understand the ways that specific gene regulatory mechanisms contribute to adaptive divergence of populations. In plants, the post-transcriptional regulatory mechanism of alternative splicing (AS) plays an important role in both development and abiotic stress response, making it a compelling potential target of natural selection. AS allows organisms to generate multiple different transcripts/proteins from a single gene and thus may provide a source of evolutionary novelty. Here, we examine whether variation in alternative splicing and gene expression levels might contribute to adaptation and incipient speciation of dune-adapted prairie sunflowers in Great Sand Dunes National Park, Colorado, USA. We conducted a common garden experiment to assess transcriptomic variation among ecotypes and analyzed differential expression, differential splicing, and gene coexpression. We show that individual genes are strongly differentiated for both transcript level and alternative isoform proportions, even when grown in a common environment, and that gene coexpression networks are disrupted between ecotypes. Furthermore, we examined how genome-wide patterns of sequence divergence correspond to divergence in transcript levels and isoform proportions and find evidence for both cis and trans-regulation. Together, our results emphasize that alternative splicing has been an underappreciated mechanism providing source material for natural selection at short evolutionary time scales.

Quantification and modeling of turnover dynamics of de novo transcripts in Drosophila melanogaster.

Most of the transcribed eukaryotic genomes are composed of non-coding transcripts. Among these transcripts, some are newly transcribed when compared to outgroups and are referred to as de novo transcripts. De novo transcripts have been shown to play a major role in genomic innovations. However, little is known about the rates at which de novo transcripts are gained and lost in individuals of the same species. Here, we address this gap and estimate the de novo transcript turnover rate with an evolutionary model. We use DNA long reads and RNA short reads from seven geographically remote samples of inbred individuals of Drosophila melanogaster to detect de novo transcripts that are gained on a short evolutionary time scale. Overall, each sampled individual contains around 2500 unspliced de novo transcripts, with most of them being sample specific. We estimate that around 0.15 transcripts are gained per year, and that each gained transcript is lost at a rate around 5× 10-5 per year. This high turnover of transcripts suggests frequent exploration of new genomic sequences within species. These rate estimates are essential to comprehend the process and timescale of de novo gene birth.

Concurrent invasions of European starlings in Australia and North America reveal population‐specific differentiation in shared genomic regions

AbstractA species' success during the invasion of new areas hinges on an interplay between the demographic processes common to invasions and the specific ecological context of the novel environment. Evolutionary genetic studies of invasive species can investigate how genetic bottlenecks and ecological conditions shape genetic variation in invasions, and our study pairs two invasive populations that are hypothesized to be from the same source population to compare how each population evolved during and after introduction. Invasive European starlings (Sturnus vulgaris) established populations in both Australia and North America in the 19th century. Here, we compare whole‐genome sequences among native and independently introduced European starling populations to determine how demographic processes interact with rapid evolution to generate similar genetic patterns in these recent and replicated invasions. Demographic models indicate that both invasive populations experienced genetic bottlenecks as expected based on invasion history, and we find that specific genomic regions have differentiated even on this short evolutionary timescale. Despite genetic bottlenecks, we suggest that genetic drift alone cannot explain differentiation in at least two of these regions. The demographic boom intrinsic to many invasions as well as potential inversions may have led to high population‐specific differentiation, although the patterns of genetic variation are also consistent with the hypothesis that this infamous and highly mobile invader adapted to novel selection (e.g., extrinsic factors). We use targeted sampling of replicated invasions to identify and evaluate support for multiple, interacting evolutionary mechanisms that lead to differentiation during the invasion process.

How the switch to hummingbird pollination has greatly contributed to our understanding of evolutionary processes.

The evolutionary switch to hummingbird pollination exemplifies complex adaptation, requiring evolutionary change in multiple component traits. Despite this complexity, diverse lineages have converged on hummingbird-adapted flowers on a relatively short evolutionary timescale. Here, I review how features of the genetic basis of adaptation contribute to this remarkable evolutionary lability. Large-effect substitutions, large mutational targets for adaptation, adaptive introgression, and concentrated architecture all contribute to the origin and maintenance of hummingbird-adapted flowers. The genetic features of adaptation are likely shaped by the ecological and geographic context of the switch to hummingbird pollination, with implications for future evolutionary trajectories.

Exploring the evolution and function of Canoe's intrinsically disordered region in linking cell-cell junctions to the cytoskeleton during embryonic morphogenesis.

One central question for cell and developmental biologists is defining how epithelial cells can change shape and move during embryonic development without tearing tissues apart. This requires robust yet dynamic connections of cells to one another, via the cell-cell adherens junction, and of junctions to the actin and myosin cytoskeleton, which generates force. The last decade revealed that these connections involve a multivalent network of proteins, rather than a simple linear pathway. We focus on Drosophila Canoe, homolog of mammalian Afadin, as a model for defining the underlying mechanisms. Canoe and Afadin are complex, multidomain proteins that share multiple domains with defined and undefined binding partners. Both also share a long carboxy-terminal intrinsically disordered region (IDR), whose function is less well defined. IDRs are found in many proteins assembled into large multiprotein complexes. We have combined bioinformatic analysis and the use of a series of canoe mutants with early stop codons to explore the evolution and function of the IDR. Our bioinformatic analysis reveals that the IDRs of Canoe and Afadin differ dramatically in sequence and sequence properties. When we looked over shorter evolutionary time scales, we identified multiple conserved motifs. Some of these are predicted by AlphaFold to be alpha-helical, and two correspond to known protein interaction sites for alpha-catenin and F-actin. We next identified the lesions in a series of eighteen canoe mutants, which have early stop codons across the entire protein coding sequence. Analysis of their phenotypes are consistent with the idea that the IDR, including the conserved motifs in the IDR, are critical for protein function. These data provide the foundation for further analysis of IDR function.

Individual identity information persists in learned calls of introduced parrot populations.

Animals can actively encode different types of identity information in learned communication signals, such as group membership or individual identity. The social environments in which animals interact may favor different types of information, but whether identity information conveyed in learned signals is robust or responsive to social disruption over short evolutionary timescales is not well understood. We inferred the type of identity information that was most salient in vocal signals by combining computational tools, including supervised machine learning, with a conceptual framework of "hierarchical mapping", or patterns of relative acoustic convergence across social scales. We used populations of a vocal learning species as a natural experiment to test whether the type of identity information emphasized in learned vocalizations changed in populations that experienced the social disruption of introduction into new parts of the world. We compared the social scales with the most salient identity information among native and introduced range monk parakeet (Myiopsitta monachus) calls recorded in Uruguay and the United States, respectively. We also evaluated whether the identity information emphasized in introduced range calls changed over time. To place our findings in an evolutionary context, we compared our results with another parrot species that exhibits well-established and distinctive regional vocal dialects that are consistent with signaling group identity. We found that both native and introduced range monk parakeet calls displayed the strongest convergence at the individual scale and minimal convergence within sites. We did not identify changes in the strength of acoustic convergence within sites over time in the introduced range calls. These results indicate that the individual identity information in learned vocalizations did not change over short evolutionary timescales in populations that experienced the social disruption of introduction. Our findings point to exciting new research directions about the robustness or responsiveness of communication systems over different evolutionary timescales.

Dynamics and stage-specificity of between-population gene expression divergence in the Drosophila melanogaster larval fat body

Gene expression variation is pervasive across all levels of organismal organization, including development. Few studies, however, have examined variation in developmental transcriptional dynamics among populations, or how it contributes to phenotypic divergence. Indeed, the evolution of gene expression dynamics when both the evolutionary and temporal timescale are comparatively short remains relatively uncharacterized. Here, we examined coding and non-coding gene expression in the fat body of an ancestral African and a derived European Drosophila melanogaster population across three developmental stages spanning ten hours of larval development. Between populations, expression divergence was largely stage-specific. We detected higher expression variation during the late wandering stage, which may be a general feature of this stage. During this stage, we also detected higher and more extensive lncRNA expression in Europe, suggesting that lncRNA expression may be more important in derived populations. Interestingly, the temporal breadth of protein-coding and lncRNA expression became more restricted in the derived population. Taken together with the signatures of potential local adaptation that we detected at the sequence level in 9-25% of candidate genes (those showing evidence of expression divergence between populations), this finding suggests that gene expression becomes more developmental stage-specific during adaptation to new environments. We further used RNAi to identify several candidate genes that likely contribute to known phenotypic divergence between these populations. Our results shed light on the evolution and dynamics of expression variation over short developmental and evolutionary timescales, and how this variation contributes to population and phenotypic divergence.

The Evolution of Mitochondrial Genomes between Two Cymbidium Sister Species: Dozens of Circular Chromosomes and the Maintenance and Deterioration of Genome Synteny.

Plant mitochondrial genomes (mitogenomes) exhibit fluid genome architectures, which could lead to the rapid erosion of genome synteny over a short evolutionary time scale. Among the species-rich orchid family, the leafy Cymbidium lancifolium and leafless Cymbidium macrorhizon are sister species with remarkable differences in morphology and nutritional physiology. Although our understanding of the evolution of mitochondria is incomplete, these sister taxa are ideal for examining this subject. In this study, the complete mitogenomes of C. lancifolium and C. macrorhizon, totaling 704,244 bp and 650,751 bp, respectively, were assembled. In the 2 mitogenomes, 38 protein-coding genes, 18 cis- and 6 trans-spliced introns, and approximately 611 Kb of homologous sequences are identical; overall, they have 99.4% genome-wide similarity. Slight variations in the mitogenomes of C. lancifolium and C. macrorhizon in repeat content (21.0 Kb and 21.6 Kb, respectively) and mitochondrial DNA of plastid origin (MIPT; 38.2 Kb and 37.5 Kb, respectively) were observed. The mitogenome architectures of C. lancifolium and C. macrorhizon are complex and comprise 23 and 22 mini-circular chromosomes, respectively. Pairwise comparisons indicate that the two mitogenomes are largely syntenic, and the disparity in chromosome numbers is likely due to repeat-mediated rearrangements among different chromosomes. Notably, approximately 93.2 Kb C. lancifolium mitochondrial sequences lack any homology in the C. macrorhizon mitogenome, indicating frequent DNA gains and losses, which accounts mainly for the size variation. Our findings provide unique insights into mitogenome evolution in leafy and leafless plants of sister species and shed light on mitogenome dynamics during the transition from mixotrophy to mycoheterotrophy.

Sexually antagonistic selection maintains genetic variance when sexual dimorphism evolves.

Genetic variance (VG) in fitness related traits is often unexpectedly high, evoking the question how VG can be maintained in the face of selection. Sexually antagonistic (SA) selection favouring alternative alleles in the sexes is common and predicted to maintain VG, while directional selection should erode it. Both SA and sex-limited directional selection can lead to sex-specific adaptations but how each affect VG when sexual dimorphism evolves remain experimentally untested. Using replicated artificial selection on the seed beetle Callosobruchus maculatus body size we recently demonstrated an increase in size dimorphism under SA and male-limited (ML) selection by 50% and 32%, respectively. Here we test their consequences on genetic variation. We show that SA selection maintained significantly more ancestral, autosomal additive genetic variance than ML selection, while both eroded sex-linked additive variation equally. Ancestral female-specific dominance variance was completely lost under ML, while SA selection consistently sustained it. Further, both forms of selection preserved a high genetic correlation between the sexes (rm,f). These results demonstrate the potential for sexual antagonism to maintain more genetic variance while fuelling sex-specific adaptation in a short evolutionary time scale, and are in line with predicted importance of sex-specific dominance reducing sexual conflict over alternative alleles.

Evolution of female song and duetting in the chaffinch (Fringilla) species complex

Female song is ancestral to songbirds and shows considerable phylogenetic signal, but its presence also appears to be labile and correlated with life‐history and ecology. While previous studies have examined the evolution of female song across species‐rich families, here we studied female song in island populations of a recently diverged species‐complex, the chaffinches (genus Fringilla). We show this behaviour has recently evolved in these populations, and probably on two independent occasions. In F. teydea, female song is performed in loose duets with males, while in F. coelebs, female song is produced solo. Populations with singing females showed year‐round territoriality and were found in regions with low seasonality – both factors previously connected with high rates of female singing. To determine the relative saliency of female songs to conspecifics, and the degree to which they instigate territorial defence behaviours, we performed a series of speaker playback experiments. In F. c. canariensis, female song could induce comparable responses to male song, whilst duetting F. teydea playback produced similar responses to solo songs, and therefore may relate to within‐pair communication instead. Our results suggest female song can be a highly labile trait that can evolve over short evolutionary timescales.

Evolutionary conservation and diversification of auditory neural circuits that process courtship songs in Drosophila

Acoustic communication signals diversify even on short evolutionary time scales. To understand how the auditory system underlying acoustic communication could evolve, we conducted a systematic comparison of the early stages of the auditory neural circuit involved in song information processing between closely-related fruit-fly species. Male Drosophila melanogaster and D. simulans produce different sound signals during mating rituals, known as courtship songs. Female flies from these species selectively increase their receptivity when they hear songs with conspecific temporal patterns. Here, we firstly confirmed interspecific differences in temporal pattern preferences; D. simulans preferred pulse songs with longer intervals than D. melanogaster. Primary and secondary song-relay neurons, JO neurons and AMMC-B1 neurons, shared similar morphology and neurotransmitters between species. The temporal pattern preferences of AMMC-B1 neurons were also relatively similar between species, with slight but significant differences in their band-pass properties. Although the shift direction of the response property matched that of the behavior, these differences are not large enough to explain behavioral differences in song preferences. This study enhances our understanding of the conservation and diversification of the architecture of the early-stage neural circuit which processes acoustic communication signals.

Testing the success of palaeontological methods in the delimitation of clam shrimp (Crustacea, Branchiopoda) on extant species

AbstractFossil spinicaudatan taxonomy heavily relies on carapace features (size, shape, ornamentation) and palaeontologists have greatly refined methods to study and describe carapace variability. Whether carapace features alone are sufficient for distinguishing between species of a single genus has remained untested. In our study, we tested common palaeontological methods on 481 individuals of the extant Australian genus Ozestheria that have been previously assigned to ten species based on genetic analysis. All species are morphologically distinct based on geometric morphometrics (p ≤ 0.001), but they occupy overlapping regions in Ozestheria morphospace. Linear discriminant analysis of Fourier shape coefficients reaches a mean model performance of 93.8% correctly classified individuals over all possible 45 pairwise species comparisons. This can be further increased by combining the size and shape datasets. Nine of the ten examined species are clearly sexually dimorphic but male and female morphologies strongly overlap within species with little influence on model performance. Ornamentation is commonly species‐diagnostic; seven ornamentation types are distinguished of which six are species‐specific while one is shared by four species. A transformation of main ornamental features (e.g. from punctate to smooth) can occur among closely related species suggesting short evolutionary timescales. Our overall results support the taxonomic value of carapace features, which should also receive greater attention in the taxonomy of extant species. The extensive variation in carapace shape and ornamentation is noteworthy and several species would probably have been assigned to different genera or families if these had been fossils, bearing implications for the systematics of fossil Spinicaudata.

Neural networks: solving the chemistry of the interstellar medium

ABSTRACT Non-equilibrium chemistry is a key process in the study of the interstellar medium (ISM), in particular the formation of molecular clouds and thus stars. However, computationally, it is among the most difficult tasks to include in astrophysical simulations, because of the typically high (&amp;gt;40) number of reactions, the short evolutionary time-scales (about 104 times less than the ISM dynamical time), and the characteristic non-linearity and stiffness of the associated ordinary differential equations system (ODEs). In this proof of concept work, we show that Physics Informed Neural Networks (PINN) are a viable alternative to traditional ODE time integrators for stiff thermochemical systems, i.e. up to molecular hydrogen formation (9 species and 46 reactions). Testing different chemical networks in a wide range of densities (−2 &amp;lt; log n/cm−3 &amp;lt; 3) and temperatures (1 &amp;lt; log T/K &amp;lt; 5), we find that a basic architecture can give a comfortable convergence only for simplified chemical systems: to properly capture the sudden chemical and thermal variations, a Deep Galerkin Method is needed. Once trained (∼103 GPUhr), the PINN well reproduces the strong non-linear nature of the solutions (errors $\lesssim 10{{\ \rm per\ cent}}$) and can give speed-ups up to a factor of ∼200 with respect to traditional ODE solvers. Further, the latter have completion times that vary by about $\sim 30{{\ \rm per\ cent}}$ for different initial n and T, while the PINN method gives negligible variations. Both the speed-up and the potential improvement in load balancing imply that PINN-powered simulations are a very palatable way to solve complex chemical calculation in astrophysical and cosmological problems.

Polygenic plague resistance in the great gerbil uncovered by population sequencing.

Pathogens can elicit high selective pressure on hosts, potentially altering genetic diversity over short evolutionary timescales. Intraspecific variation in immune response is observable as variable survivability from specific infections. The great gerbil (Rhombomys opimus) is a rodent plague host with a heterogenic but highly resistant phenotype. Here, we investigate the genomic basis for plague-resistant phenotypes by exposing wild-caught great gerbils to plague (Yersinia pestis). Whole genome sequencing of 10 survivors and 10 moribund individuals revealed a subset of genomic regions showing elevated differentiation. Gene ontology analysis of candidate genes in these regions demonstrated enrichment of genes directly involved in immune functions, cellular metabolism and the regulation of apoptosis as well as pathways involved in transcription, translation, and gene regulation. Transcriptomic analysis revealed that the early activated great gerbil immune response to plague consisted of classical components of the innate immune system. Our approach combining challenge experiments with transcriptomics and population level sequencing, provides new insight into the genetic background of plague-resistance and confirms its complex nature, most likely involving multiple genes and pathways of both the immune system and regulation of basic cellular functions.