An association between differential expression and genetic divergence in the Patagonian olive mouse (Abrothrix olivacea)

Phylogenomics and species delimitation of an abundant and little-studied Amazonian forest spiny rat

Species delimitation studies based on integrating different datasets such as genomic, morphometric, and cytogenetics data are rare in studies focused on Neotropical rodents. As a consequence, the evolutionary history of most of these genera remains poorly understood. Proechimys is a highly diverse and widely distributed genus of Neotropical spiny rats with unique traits like multiple sympatry, micro-habitat segregation, and fuzzy species limits. Here, we applied RAD-Seq to infer the phylogenetic relationships, estimate the species boundaries, and estimate the divergence times for Proechimys, one of the most common and least studied small mammals in the Amazon. We tested whether inferred lineages in the phylogenetic trees could be considered distinct species based on the genomic dataset and morphometric data. Analyses revealed the genus is not monophyletic, with Proechimys hoplomyoides sister to a group of Hoplomys gymnurus + all other Proechimys species, contesting the generic status of Hoplomys. There are five main clades in Proechimys stricto sensu (excluding H. gymnurus and P. hoplomyoides). Species delimitation analyses supported 25 species within the genus Proechimys. The five main clades in Proechimys stricto sensu also showed similar ages for their origins, and two rapid diversification events were identified in the Early Pliocene and in the Early Pleistocene. Most cases of sympatry in Proechimys occur among species from the different main clades, and although Proechimys is an inhabitant of the Amazon, three species occupied the Cerrado biome during the Pleistocene. We could associate available nominal taxon, cytogenetics information, and DNA sequences in Genbank to most of the 25 species we hypothesized from our delimitation analyses. Based on our analyses, we estimate that eight forms represent putative new species that need a taxonomic revision.

Climate influences the genetic structure and niche differentiation among populations of the olive field mouse Abrothrix olivacea (Cricetidae: Abrotrichini)

Even when environmental variation over time and space is commonly considered as an important driver of population divergence, few evaluations of intraspecific genetic variation explicitly assess whether observed structure has been caused by or is correlated with landscape heterogeneity. Several phylogeographic studies have characterized the mitochondrial diversity of Abrothrix olivacea, but none has incorporated landscape genetics analyses and ecological niche modeling, leaving a gap in the understanding of the species evolutionary history. Here, these aspects were addressed based on 186 single nucleotide polymorphisms, extracted from sequences of 801 bp of Cytb gene, gathered from 416 individuals collected at 103 localities in Argentina and Chile. Employing multivariate statistical analyses (gPCA, Mantel and Partial Mantel Tests, Procrustes Analysis, and RDA), associations between genetic differences and geographic and climatic distances were evaluated. Presence data was employed to estimate the potential geographic distribution of this species during historical and contemporary climatic scenarios, and to address differences among the climatic niches of their main mitochondrial lineages. The significant influence of landscape features in structuring mitochondrial variability was evidenced at different spatial scales, as well as the role of past climatic dynamics in driving geographic range shifts, mostly associated to Quaternary glaciations. Overall, these results suggest that throughout geographic range gene flow is unevenly influenced by climatic dissimilarity and the geographic distancing, and that studied lineages do not exhibit distributional signals of climatic niche conservatism. Additionally, genetic differentiation occurred by more complex evolutionary processes than mere disruption of gene flow or drift.

A South American Mouse Morbillivirus Provides Insight into a Clade of Rodent-Borne Morbilliviruses

Morbilliviruses are negative-sense single-stranded monosegmented RNA viruses in the family Paramyxoviridae (order Mononegavirales). Morbilliviruses infect diverse mammals including humans, dogs, cats, small ruminants, seals, and cetaceans, which serve as natural hosts. Here, I report the identification and characterization of novel viruses detected in public RNAseq datasets of South American long-haired and olive field mice. The divergent viruses dubbed Ratón oliváceo morbillivirus (RoMV) detected in renal samples from mice collected from Chile and Argentina are characterized by an unusually large genome including long intergenic regions and the presence of an accessory protein between the F and H genes redounding in a genome architecture consisting in 3'-N-P/V/C-M-F-hp-H-L-5'. Structural and functional annotation, genetic distance, and evolutionary insights suggest that RoMV is a member of a novel species within genus Morbillivirus tentatively named as South American mouse morbillivirus. Phylogenetic analysis suggests that this mouse morbillivirus is closely related to and clusters into a monophyletic group of novel rodent-borne morbilliviruses. This subclade of divergent viruses expands the host range, redefines the genomic organization and provides insights on the evolutionary history of genus Morbillivirus.

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.

Transcriptomic signatures of temperature adaptation in the eastern oyster Crassostrea virginica

The large geographic distribution of the eastern oyster, Crassostrea virginica, makes it an ideal species to test how populations have adapted to latitudinal gradients in temperature. Despite inhabiting distinct thermal regimes, populations of C. virginica near the species' southern and northern geographic range show no population differences in their physiological response to temperature. In this study, we used comparative transcriptomics to understand how oysters from either end of the species' range maintain enantiostasis across three acclimation temperatures (10, 20, and 30°C). With this approach, we identified genes that were differentially expressed in response to temperature between individuals of C. virginica collected from New Brunswick, Canada and Louisiana, USA. We observed a core set of genes whose expression responded to temperature in both populations, but also an even larger set of genes with expression patterns that were unique to each population. Intriguingly, the genes with population-specific responses to temperature had elevated F_ST and Ka/Ks ratios compared to the genome-wide average. In contrast, genes showing only a response to temperature were found to only have elevated F_ST values suggesting that divergent F_ST may be due to selection on linked regulatory regions rather than positive selection on protein coding regions. Taken together, our results suggest that, despite coarse-scale physiological similarities, natural selection has shaped divergent gene expression responses to temperature in geographically separated populations of this broadly eurythermal marine invertebrate.

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.

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.

Genetic Diversity and Demographic History of the Shaggy Soft-Haired Mouse Abrothrix hirta (Cricetidae; Abrotrichini)

Genetic information on species can inform decision making regarding conservation of biodiversity since the response of organisms to changing environments depend, in part, on their genetic makeup. Territories of central-southern Chile and Argentina have undergone a varying degree of impact during the Quaternary, where the response of local fauna and flora was rather species-specific. Here, we focus on the sigmodontine rodent Abrothrix hirta, distributed from 35° S in Chile and Argentina to northern Tierra del Fuego. Based on 119,226 transcriptome-derived SNP loci from 46 individuals of A. hirta, we described the geographic distribution of the genetic diversity of this species using a maximum likelihood tree, principal component and admixture analyses. We also addressed the demographic history of the main intraspecific lineages of A. hirta using GADMA. We found that A. hirta exhibited four allopatric intraspecific lineages. Three main genetic groups were identified by a Principal Component Analysis and by Ancestry analysis. The demographic history of A. hirta was characterized by recent population stability for populations at the northernmost part of the range, while southern populations experienced a recent population expansion.

Comparative transcriptome analysis of the newly discovered insect vector of the pine wood nematode in China, revealing putative genes related to host plant adaptation

BACKGROUND

In many insect species, the larvae/nymphs are unable to disperse far from the oviposition site selected by adults. The Sakhalin pine sawyer Monochamus saltuarius (Gebler) is the newly discovered insect vector of the pine wood nematode (Bursaphelenchus xylophilus) in China. Adult M. saltuarius prefers to oviposit on the host plant Pinus koraiensis, rather than P. tabuliformis. However, the genetic basis of adaptation of the larvae of M. saltuarius with weaken dispersal ability to host environments selected by the adult is not well understood.

RESULTS

In this study, the free amino and fatty acid composition and content of the host plants of M. saltuarius larvae, i.e., P. koraiensis and P. tabuliformis were investigated. Compared with P. koraiensis, P. tabuliformis had a substantially higher content of various free amino acids, while the opposite trend was detected for fatty acid content. The transcriptional profiles of larval populations feeding on P. koraiensis and P. tabuliformis were compared using PacBio Sequel II sequencing combined with Illumina sequencing. The results showed that genes relating to digestion, fatty acid synthesis, detoxification, oxidation-reduction, and stress response, as well as nutrients and energy sensing ability, were differentially expressed, possibly reflecting adaptive changes of M. saltuarius in response to different host diets. Additionally, genes coding for cuticle structure were differentially expressed, indicating that cuticle may be a potential target for plant defense. Differential regulation of genes related to the antibacterial and immune response were also observed, suggesting that larvae of M. saltuarius may have evolved adaptations to cope with bacterial challenges in their host environments.

CONCLUSIONS

The present study provides comprehensive transcriptome resource of M. saltuarius relating to host plant adaptation. Results from this study help to illustrate the fundamental relationship between transcriptional plasticity and adaptation mechanisms of insect herbivores to host plants.

Genetic variation of the Chilean endemic long-haired mouse Abrothrix longipilis (Rodentia, Supramyomorpha, Cricetidae) in a geographical and environmental context

Quaternary climate and associated vegetational changes affected the fauna of the Chilean Mediterranean ecosystem. Here we studied the genetic variation of the long-haired mouse, Abrothrix longipilis, a sigmodontine rodent endemic to this area. Within an environmentally explicit context, we examined the geographic distribution of the genetic diversity and demographic history of the species based on sequences of the mitochondrial Cytochrome-b gene of 50 individuals from 13 localities and a large panel of single nucleotide polymorphisms of 17 individuals from 6 localities. The gene genealogy of A. longipilis revealed three intraspecific lineages that are allopatric and latitudinally segregated (northern, central, and southern lineages) with an estimated crown age for the whole species clade of 552.3 kyr B.P. A principal component analysis based on 336,596 SNP loci is in line with the information given by the the mitochondrial gene genealogy. Along its complete distributional range, A. longipilis showed patterns of isolation by distance and also isolation by environment. The general pattern of historical demography showed stability for most intraspecific lineages of A. longipilis. Northern and central lineages showed signals of historical demographic stability, while the southern lineage showed contrasting signals. In agreement with this, the niche models performed showed that in the northern range of A. longipilis, areas of high suitability for this species increased towards the present time; areas of central range would have remained relatively stable, while southern areas would have experienced more change through time. In summary, our study shows three distinct allopatric lineages of A. longipilis, each showing slightly different demographic history.

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.

The transcriptional correlates of divergent electric organ discharges in Paramormyrops electric fish

Background

Understanding the genomic basis of phenotypic diversity can be greatly facilitated by examining adaptive radiations with hypervariable traits. In this study, we focus on a rapidly diverged species group of mormyrid electric fish in the genus Paramormyrops, which are characterized by extensive phenotypic variation in electric organ discharges (EODs). The main components of EOD diversity are waveform duration, complexity and polarity. Using an RNA-sequencing based approach, we sought to identify gene expression correlates for each of these EOD waveform features by comparing 11 specimens of Paramormyrops that exhibit variation in these features.

Results

Patterns of gene expression among Paramormyrops are highly correlated, and 3274 genes (16%) were differentially expressed. Using our most restrictive criteria, we detected 145–183 differentially expressed genes correlated with each EOD feature, with little overlap between them. The predicted functions of several of these genes are related to extracellular matrix, cation homeostasis, lipid metabolism, and cytoskeletal and sarcomeric proteins. These genes are of significant interest given the known morphological differences between electric organs that underlie differences in the EOD waveform features studied.

Conclusions

In this study, we identified plausible candidate genes that may contribute to phenotypic differences in EOD waveforms among a rapidly diverged group of mormyrid electric fish. These genes may be important targets of selection in the evolution of species-specific differences in mate-recognition signals.

Environment predicts repeated body size shifts in a recent radiation of Australian mammals

Closely related species that occur across steep environmental gradients often display clear body size differences, and examining this pattern is crucial to understanding how environmental variation shapes diversity. Australian endemic rodents in the Pseudomys Division (Muridae: Murinae) have repeatedly colonized the arid, monsoon, and mesic biomes over the last 5 million years. Using occurrence records, body mass data, and Bayesian phylogenetic models, we test whether body mass of 31 species in the Pseudomys Division can be predicted by their biome association. We also model the effect of eight environmental variables on body mass. Despite high phylogenetic signal in body mass evolution across the phylogeny, we find that mass predictably increases in the mesic biome and decreases in arid and monsoon biomes. As per Bergmann's rule, temperature is strongly correlated with body mass, as well as several other variables. Our results highlight two important findings. First, body size in Australian rodents has tracked with climate through the Pleistocene, likely due to several environmental variables rather than a single factor. Second, support for both Brownian motion and predictable change at different taxonomic levels in the Pseudomys Division phylogeny demonstrates how the level at which we test hypotheses can alter interpretation of evolutionary processes.

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.

Contrasting patterns of divergence at the regulatory and sequence level in European Daphnia galeata natural populations

Understanding the genetic basis of local adaptation has long been a focus of evolutionary biology. Recently, there has been increased interest in deciphering the evolutionary role of Daphnia's plasticity and the molecular mechanisms of local adaptation. Using transcriptome data, we assessed the differences in gene expression profiles and sequences in four European Daphnia galeata populations. In total, ~33% of 32,903 transcripts were differentially expressed between populations. Among 10,280 differentially expressed transcripts, 5,209 transcripts deviated from neutral expectations and their population-specific expression pattern is likely the result of local adaptation processes. Furthermore, a SNP analysis allowed inferring population structure and distribution of genetic variation. The population divergence at the sequence level was comparatively higher than the gene expression level by several orders of magnitude consistent with strong founder effects and lack of gene flow between populations. Using sequence homology, the candidate transcripts were annotated using a comparative genomics approach. Additionally, we also performed a weighted gene co-expression analysis to identify population-specific regulatory patterns of transcripts in D. galeata. Thus, we identified candidate transcriptomic regions for local adaptation in this key species of aquatic ecosystems in the absence of any laboratory-induced stressor.