A priori and a posteriori approaches for finding genes of evolutionary interest in non-model species: osmoregulatory genes in the kidney transcriptome of the desert rodent Dipodomys spectabilis (banner-tailed kangaroo rat)

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.

Remotely sensed environmental measurements detect decoupled processes driving population dynamics at contrasting scales

The increasing availability of satellite imagery has supported a rapid expansion in forward-looking studies seeking to track and predict how climate change will influence wild population dynamics. However, these data can also be used in retrospect to provide additional context for historical data in the absence of contemporaneous environmental measurements. We used 167 Landsat-5 Thematic Mapper (TM) images spanning 13 years to identify environmental drivers of fitness and population size in a well-characterized population of banner-tailed kangaroo rats (Dipodomys spectabilis) in the southwestern United States. We found evidence of two decoupled processes that may be driving population dynamics in opposing directions over distinct time frames. Specifically, increasing mean surface temperature corresponded to increased individual fitness, where fitness is defined as the number of offspring produced by a single individual. This result contrasts with our findings for population size, where increasing surface temperature led to decreased numbers of active mounds. These relationships between surface temperature and (i) individual fitness and (ii) population size would not have been identified in the absence of remotely sensed data, indicating that such information can be used to test existing hypotheses and generate new ecological predictions regarding fitness at multiple spatial scales and degrees of sampling effort. To our knowledge, this study is the first to directly link remotely sensed environmental data to individual fitness in a nearly exhaustively sampled population, opening a new avenue for incorporating remote sensing data into eco-evolutionary studies.

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.

Identification of selection signatures in Iranian dromedary and Bactrian camels using whole genome sequencing data

The Old World camels play an important role as one of the main food sources in large parts of Asia and Africa. Natural selection combined with artificial selection by human has affected parts of the domestic animal genome for adapting them to their habitats and meeting human needs. Here, we used whole genome sequencing data of 34 camels (including 14 dromedaries and 20 Bactrian camels) to identify the genomic signature of selection in the Iranian dromedary (ID) and Bactrian camels (IB). To detect the mentioned regions, we used two methods including population differentiation index (Fst) and cross-population extended haplotype homozygosity (XP-EHH) with 50 kb sliding window and 25 kb step size. Based on gene ontology analysis on the candidate genes identified for IB camels, we found GO terms associated with lung development, nervous system development, immune system and behavior. Also, we identified several genes related to body thermoregulation (ZNF516), meat quality (ANK1 and HSPA13), and high-altitude adaptation (OPA1) for IB camels. In the list of detected candidate genes under selection in ID camels, the genes related to energy metabolism (BDH1), reproduction (DLG1, IMMP2L and FRASI), long-term memory (GRIA1), kidney (SLC12A1), lung development (EMILIN2 and FBN1) and immunity (SOCS2, JAK1, NRROS and SENP1) were found. Our findings, along with further studies in this field, will strengthen our knowledge about the effect of selection on the camelid genome under different geographical, climatic and even cultural conditions.

High-quality reference genome for an arid-adapted mammal, the banner-tailed kangaroo rat (Dipodomys spectabilis)

Kangaroo rats in the genus Dipodomys are found in a variety of habitat types in western North America, including deserts, arid and semiarid grasslands, and scrublands. Many Dipodomys species are experiencing strong population declines due to increasing habitat fragmentation, with two species listed as federally endangered in the United States. The precarious state of many Dipodomys populations, including those occupying extreme environments, make species of this genus valuable subjects for studying the impacts of habitat degradation and fragmentation on population genomic patterns and for characterizing the genomic bases of adaptation to harsh conditions. To facilitate exploration of such questions, we assembled and annotated a reference genome for the banner-tailed kangaroo rat (Dipodomys spectabilis) using PacBio HiFi sequencing reads, providing a more contiguous genomic resource than two previously assembled Dipodomys genomes. Using the HiFi data for D. spectabilis and publicly available sequencing data for two other Dipodomys species (Dipodomys ordii and Dipodomys stephensi), we demonstrate the utility of this new assembly for studies of congeners by conducting inference of historic effective population sizes (N_e) and linking these patterns to the species’ current extinction risk statuses. The genome assembly presented here will serve as a valuable resource for population and conservation genomic studies of Dipodomys species, comparative genomic research within mammals and rodents, and investigations into genomic adaptation to extreme environments and changing landscapes.

Inferring the genetic responses to acute drought stress across an ecological gradient

BACKGROUND

How do xerophytic species thrive in environments that experience extreme annual drought? Although critical to the survival of many species, the genetic responses to drought stress in many non-model organisms has yet to be explored. We investigated this question in Mentzelia section Bartonia (Loasaceae), which occurs throughout western North America, including arid lands. To better understand the genetic responses to drought stress among species that occur in different habitats, the gene expression levels of three species from Mentzelia were compared across a precipitation gradient. Two de novo reference transcriptomes were generated and annotated. Leaf and root tissues were collected from control and drought shocked plants and compared to one another for differential expression. A target-gene approach was also implemented to better understand how drought-related genes from model and crop species function in non-model systems.

RESULTS

When comparing the drought-shock treatment plants to their respective control plants, we identified 165 differentially expressed clusters across all three species. Differentially expressed genes including those associated with water movement, photosynthesis, and delayed senescence. The transcriptome profiling approach was coupled with a target genes approach that measured expression of 90 genes associated with drought tolerance in model organisms. Comparing differentially expressed genes with a ≥ 2 log-fold value between species and tissue types showed significant differences in drought response. In pairwise comparisons, species that occurred in drier environments differentially expressed greater genes in leaves when drought shocked than those from wetter environments, but expression in the roots mostly produced opposite results.

CONCLUSIONS

Arid-adapted species mount greater genetic responses compared to the mesophytic species, which has likely evolved in response to consistent annual drought exposure across generations. Drought responses also depended on organ type. Xerophytes, for example, mounted a larger response in leaves to downregulate photosynthesis and senescence, while mobilizing carbon and regulating water in the roots. The complexity of drought responses in Mentzelia suggest that whole organism responses need to be considered when studying drought and, in particular, the physiological mechanisms in which plants regulate water, carbon, cell death, metabolism, and secondary metabolites.

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.

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.

Evolution of pathogen response genes associated with increased disease susceptibility during adaptation to an extreme drought in a Brassica rapa plant population

Background

Pathogens are key components in natural and agricultural plant systems. There is evidence of evolutionary changes in disease susceptibility as a consequence of climate change, but we know little about the underlying genetic basis of this evolution. To address this, we took advantage of a historical seed collection of a Brassica rapa population, which we previously demonstrated evolved an increase in disease susceptibility to a necrotrophic fungal pathogen following a drought.

Results

Previously, we combined a resurrection experiment with genome-wide sequencing of 124 pooled ancestral and descendant plants. Here, using these previously generated sequence data (Franks et al. in Mol Ecol 25(15):3622–3631, 2016), we show that well-characterized necrotrophic fungal pathogen response (NFPR) genes have evolved, as indicated by changes in allele frequency, between ancestors and descendants, with several of them identified as extreme F_ST outliers. The jasmonic acid (JA) signaling pathway in particular seems to underlie the evolution of disease susceptibility, in addition to its well characterized role in plastic disease response. We identify a list of 260 genes that are both NFPR genes and are differentially expressed in response to drought, based on publicly available data. We present evidence that five of these genes evolved between ancestors and descendants, suggesting that the drought acted as the evolutionary driver, and that the accompanying increase in disease susceptibility may have been a consequence of genetic pleiotropy.

Conclusions

Our study provides evidence that for this population, standing variation in NFPR genes is affected by natural selection related to climate change. Our results reveal potentially important candidates that may underlie trait evolution in both crops and natural systems. Additionally, this trade-off between adaptation to biotic and abiotic stresses is an example of how climate change can have diverse and unexpected consequences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01789-7.

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.

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.

De novo comparative transcriptome analysis of a rare cicada, with identification of candidate genes related to adaptation to a novel host plant and drier habitats

Background

Although the importance of host plant chemistry in plant–insect interactions is widely recognized, our understanding about the genetic basis underlying the relationship between changes in midgut proteins and adaptation of plant-feeding insects to novel host plants and habitats is very limited. To address this knowledge gap, the transcriptional profiles of midguts among three populations of the cicada Subpsaltria yangi Chen were compared. Among which, the Hancheng (HC) and Fengxiang (FX) populations occurring in the Loess Plateau feed on Ziziphus jujuba Mill. var. spinosa (Bunge) Hu ex H. F. Chow, while the population occurring in a much drier habitat in the Helan (HL) Mountains is locally specialized on a chemically divergent plant, Ephedra lepidosperma C. Y. Cheng.

Results

Based on comparative analysis, 1826 (HL vs HC) differentially expressed genes (DEGs) and 723 DEGs (HL vs FX) were identified between the populations utilizing different host plants, including 20, 36, 2, 5 and 2 genes related to digestion, detoxification, oxidation-reduction, stress response and water-deprivation response, respectively, and 35 genes presumably associated with osmoregulation. However, only 183 DEGs were identified between the HC and FX populations, including two genes related to detoxification, two genes related to stress response, and one gene presumably associated with osmoregulation. These results suggest that the weakest expression differences were between the populations utilizing the same host plant and occurring in the closest habitats, which may help explain the metabolic mechanism of adaptation in S. yangi populations to novel host plants and new niches.

Conclusions

The observed differences in gene expression among S. yangi populations are consistent with the hypothesis that the host plant shift and habitat adaptation in the HL population was facilitated by differential regulation of genes related to digestion, detoxification, oxidation-reduction, stress response, water-deprivation response and osmoregulation. The results may inform future studies on the molecular mechanisms underlying the relationship between changes in midgut proteins and adaptation of herbivorous insects to novel host plants and new niches.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5547-y) contains supplementary material, which is available to authorized users.

Physiological and biochemical changes associated with acute experimental dehydration in the desert adapted mouse, Peromyscus eremicus

Characterizing traits critical for adaptation to a given environment is an important first step in understanding how phenotypes evolve. How animals adapt to the extreme heat and aridity commonplace to deserts is an exceptionally interesting example of these processes, and has been the focus of study for decades. In contrast to those studies, where experiments are conducted on either wild animals or captive animals held in non-desert conditions, the study described here leverages a unique environmental chamber that replicates desert conditions for captive Peromyscus eremicus (cactus mouse). Here, we establish baseline values for daily water intake and for serum electrolytes, as well as the response of these variables to acute experimental dehydration. In brief, P eremicus daily water intake is very low. Its serum electrolytes are distinct from many previously studied animals, and its response to acute dehydration is profound, though not suggestive of renal impairment, which is atypical of mammals.

Characterizing the reproductive transcriptomic correlates of acute dehydration in males in the desert-adapted rodent, Peromyscus eremicus

BACKGROUND

The understanding of genomic and physiological mechanisms related to how organisms living in extreme environments survive and reproduce is an outstanding question facing evolutionary and organismal biologists. One interesting example of adaptation is related to the survival of mammals in deserts, where extreme water limitation is common. Research on desert rodent adaptations has focused predominantly on adaptations related to surviving dehydration, while potential reproductive physiology adaptations for acute and chronic dehydration have been relatively neglected. This study aims to explore the reproductive consequences of acute dehydration by utilizing RNAseq data in the desert-specialized cactus mouse (Peromyscus eremicus).

RESULTS

We exposed 22 male cactus mice to either acute dehydration or control (fully hydrated) treatment conditions, quasimapped testes-derived reads to a cactus mouse testes transcriptome, and then evaluated patterns of differential transcript and gene expression. Following statistical evaluation with multiple analytical pipelines, nine genes were consistently differentially expressed between the hydrated and dehydrated mice. We hypothesized that male cactus mice would exhibit minimal reproductive responses to dehydration; therefore, this low number of differentially expressed genes between treatments aligns with current perceptions of this species' extreme desert specialization. However, these differentially expressed genes include Insulin-like 3 (Insl3), a regulator of male fertility and testes descent, as well as the solute carriers Slc45a3 and Slc38a5, which are membrane transport proteins that may facilitate osmoregulation.

CONCLUSIONS

These results suggest that in male cactus mice, acute dehydration may be linked to reproductive modulation via Insl3, but not through gene expression differences in the subset of other a priori tested reproductive hormones. Although water availability is a reproductive cue in desert-rodents exposed to chronic drought, potential reproductive modification via Insl3 in response to acute water-limitation is a result which is unexpected in an animal capable of surviving and successfully reproducing year-round without available external water sources. Indeed, this work highlights the critical need for integrative research that examines every facet of organismal adaptation, particularly in light of global climate change, which is predicted, amongst other things, to increase climate variability, thereby exposing desert animals more frequently to the acute drought conditions explored here.

Severe acute dehydration in a desert rodent elicits a transcriptional response that effectively prevents kidney injury

Animals living in desert environments are forced to survive despite severe heat, intense solar radiation, and both acute and chronic dehydration. These animals have evolved phenotypes that effectively address these environmental stressors. To begin to understand the ways in which the desert-adapted rodent Peromyscus eremicus survives, reproductively mature adults were subjected to 72 h of water deprivation, during which they lost, on average, 23% of their body weight. The animals reacted via a series of changes in the kidney, which included modulating expression of genes responsible for reducing the rate of transcription and maintaining water and salt balance. Extracellular matrix turnover appeared to be decreased, and apoptosis was limited. In contrast to the canonical human response, serum creatinine and other biomarkers of kidney injury were not elevated, suggesting that changes in gene expression related to acute dehydration may effectively prohibit widespread kidney damage in the cactus mouse.

Characterization of a male reproductive transcriptome for Peromyscus eremicus (Cactus mouse)

Rodents of the genus Peromyscus have become increasingly utilized models for investigations into adaptive biology. This genus is particularly powerful for research linking genetics with adaptive physiology or behaviors, and recent research has capitalized on the unique opportunities afforded by the ecological diversity of these rodents. Well characterized genomic and transcriptomic data is intrinsic to explorations of the genetic architecture responsible for ecological adaptations. Therefore, this study characterizes the transcriptome of three male reproductive tissues (testes, epididymis and vas deferens) of Peromyscus eremicus (Cactus mouse), a desert specialist. The transcriptome assembly process was optimized in order to produce a high quality and substantially complete annotated transcriptome. This composite transcriptome was generated to characterize the expressed transcripts in the male reproductive tract of P. eremicus, which will serve as a crucial resource for future research investigating our hypothesis that the male Cactus mouse possesses an adaptive reproductive phenotype to mitigate water-loss from ejaculate. This study reports genes under positive selection in the male Cactus mouse reproductive transcriptome relative to transcriptomes from Peromyscus maniculatus (deer mouse) and Mus musculus. Thus, this study expands upon existing genetic research in this species, and we provide a high quality transcriptome to enable further explorations of our proposed hypothesis for male Cactus mouse reproductive adaptations to minimize seminal fluid loss.

Oxytocin and arginine vasopressin receptor evolution: implications for adaptive novelties in placental mammals

Oxytocin receptor (OXTR) and arginine vasopressin receptors (AVPR1a, AVPR1b, and AVPR2) are paralogous genes that emerged through duplication events; along the evolutionary timeline, owing to speciation, numerous orthologues emerged as well. In order to elucidate the evolutionary forces that shaped these four genes in placental mammals and to reveal specific aspects of their protein structures, 35 species were selected. Specifically, we investigated their molecular evolutionary history and intrinsic protein disorder content, and identified the presence of short linear interaction motifs. OXTR seems to be under evolutionary constraint in placental mammals, whereas AVPR1a, AVPR1b, and AVPR2 exhibit higher evolutionary rates, suggesting that they have been under relaxed or experienced positive selection. In addition, we describe here, for the first time, that the OXTR, AVPR1a, AVPR1b, and AVPR2 mammalian orthologues preserve their disorder content, while this condition varies among the paralogues. Finally, our results reveal the presence of short linear interaction motifs, indicating possible functional adaptations related to physiological and/or behavioral taxa-specific traits.

Characterization of the kidney transcriptome of the long-haired mouse Abrothrix hirta (Rodentia, Sigmodontinae) and comparison with that of the olive mouse A. olivacea

To understand how small mammals cope with the challenge of water homeostasis is a matter of interest for ecologists and evolutionary biologists. Here we take a step towards the understanding of the transcriptomic functional response of kidney using as a model the long-haired mouse (Abrothrix hirta) a species that distributes across Patagonian steppes and Austral temperate rainforests in Argentina and Chile. Specifically, we i) characterize the renal transcriptome of A. hirta, and ii) compare it with that-already available-of the co-generic and co-distributed A. olivacea. Renal mRNA transcripts from 16 specimens of A. hirta from natural populations were analyzed. Over 500 million Illumina paired-end reads were assembled de novo under two approaches, an individual assembly for each specimen, and a single in-silico normalized joint assembly including all reads from all specimens. The total number of annotated genes was similar with both strategies: an average of 14,956 in individual assemblies and 14,410 in the joint assembly. Overall, 15,463 distinct genes express in the kidney of A. hirta. Transcriptomes of A. hirta and A. olivacea were similar in terms of gene abundance and composition: 95.6% of the genes of A. hirta were also found in A. olivacea making their functional profiles also similar. However, differences in the transcriptome of these two species were observed in the set of highly expressed genes, in terms of private genes for each species and the functional profiles of highly expressed genes. As part of the novel transcriptome characterization, we provide distinct gene lists with their functional annotation that would constitute the basis for further research on these or any other species of the subfamily Sigmodontinae, which includes about 400 living species distributed from Tierra del Fuego to southern United States.

Transcriptomic characterization of the immunogenetic repertoires of heteromyid rodents

Background

When populations evolve under disparate environmental conditions, they experience different selective pressures that shape patterns of sequence evolution and gene expression. These may be manifested in genetic and phenotypic differences such as a diverse immunogenetic repertoire in species from tropical latitudes that have greater and/or different parasite burdens than more temperate species. To test this idea, we compared the transcriptomes of one tropical species (Heteromys desmarestianus) and two species from temperate latitudes (Dipodomys spectabilis and Chaetodipus baileyi) from the Heteromyidae. We did so in a search for positive selection on sequences and/or differential expression, while controlling for phylogenetic history in our choice of species.

Results

We identified 127,812 contigs and annotated 34,878 of these, identifying immune genes associated with interleukins, cytokines, and the production of mast cells. We identified 632 genes that were upregulated in H. desmarestianus (8.7% of genes tested) and 492 (6.7%) that were downregulated. Gene ontology terms including “immune response” were associated with 31 (4.9%) of the 632 upregulated genes. We found preliminary evidence for positive selection on three genes (Palmitoyltransferase ZDHHC5 Ubiquitin-conjugating enzyme E2 N, Krueppel-like factor 10, and Spindle and kinetochore-associated protein 1) along the H. desmarestianus lineage.

Conclusions

Overall our findings pinpoint genes in species from disparate environments that are on different evolutionary trajectories in terms of expression levels and/or nucleotide sequence. Our data indicate there are significant differences in the expression of genes among the spleen transcriptomes of these species and that a number of these differentially expressed genes do not show the same pattern of differential expression in another tissue type. This points to the possibility of expression differences between these species specific to the spleen transcriptome.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-929) contains supplementary material, which is available to authorized users.

Natural selection and the genetic basis of osmoregulation in heteromyid rodents as revealed by RNA-seq

One adaptation of ecological and evolutionary interest is the extraordinary ability of desert rodents to retain water during waste production. Much is known regarding the unique kidney physiology of kangaroo rats (Dipodomys spp.) and their ability to retain water during waste production, yet the genetic basis of these physiological adaptations is relatively unknown. Herein, we utilized RNA-seq data to conduct a comparative study to identify osmoregulatory genes expressed in heteromyid rodents. We sequenced kidney tissue from two temperate desert species (Dipodomys spectabilis and Chaetodipus baileyi) from two separate subfamilies of the Heteromyidae and compared these transcriptomes to a tropical mesic species (Heteromys desmarestianus) from a third subfamily. The evolutionary history of these subfamilies provided a robust phylogenetic control that allowed us to separate shared evolutionary history from convergence. Using two methods to detect differential expression (DE), we identified 1890 genes that showed consistent patterns of DE between the arid and mesic species. A three-species reciprocal BLAST analysis revealed 3511 sets of putative orthologues that, upon comparison to known Mus musculus sequences, revealed 323 annotated and full-length genic regions. Selection tests displayed evidence of positive selection (dn/ds > 1) on six genes in the two desert species and remained significant for one of these genes after correction for multiple testing. Thus, our data suggest that both the coding sequence and expression of genes have been shaped by natural selection to provide the genetic architecture for efficient osmoregulation in desert-adapted heteromyid rodents.

Microsatellite analyses across three diverse vertebrate transcriptomes (Acipenser fulvescens, Ambystoma tigrinum, and Dipodomys spectabilis)

Historically, many population genetics studies have utilized microsatellite markers sampled at random from the genome and presumed to be selectively neutral. Recent studies, however, have shown that microsatellites can occur in transcribed regions, where they are more likely to be under selection. In this study, we mined microsatellites from transcriptomes generated by 454-pyrosequencing for three vertebrate species: lake sturgeon (Acipenser fulvescens), tiger salamander (Ambystoma tigrinum), and kangaroo rat (Dipodomys spectabilis). We evaluated (i) the occurrence of microsatellites across species; (ii) whether particular gene ontology terms were over-represented in genes that contained microsatellites; (iii) whether repeat motifs were located in untranslated regions or coding sequences of genes; and (iv) in silico polymorphism. Microsatellites were less common in tiger salamanders than in either lake sturgeon or kangaroo rats. Across libraries, trinucleotides were found more frequently than any other motif type, presumably because they do not cause frameshift mutations. By evaluating variation across reads assembled to a given contig, we were able to identify repeat motifs likely to be polymorphic. Our study represents one of the first comparative data sets on the distribution of vertebrate microsatellites within expressed genes. Our results reinforce the idea that microsatellites do not always occur in noncoding DNA, but commonly occur in expressed genes.