Brain transcriptome analysis reveals gene expression differences associated with dispersal behaviour between range-front and range-core populations of invasive cane toads in Australia.

Understanding the mechanisms allowing invasive species to adapt to novel environments is a challenge in invasion biology. Many invaders demonstrate rapid evolution of behavioural traits involved in range expansion such as locomotor activity, exploration and risk-taking. However, the molecular mechanisms that underpin these changes are poorly understood. In 86 years, invasive cane toads (Rhinella marina) in Australia have drastically expanded their geographic range westward from coastal Queensland to Western Australia. During their range expansion, toads have undergone extensive phenotypic changes, particularly in behaviours that enhance the toads' dispersal ability. Common-garden experiments have shown that some changes in behavioural traits related to dispersal are heritable. At the molecular level, it is currently unknown whether these changes in dispersal-related behaviour are underlain by small or large differences in gene expression, nor is known the biological function of genes showing differential expression. Here, we used RNA-seq to gain a better understanding of the molecular mechanisms underlying dispersal-related behavioural changes. We compared the brain transcriptomes of toads from the Hawai'ian source population, as well as three distinct populations from across the Australian invasive range. We found markedly different gene expression profiles between the source population and Australian toads. By contrast, toads from across the Australian invasive range had very similar transcriptomic profiles. Yet, key genes with functions putatively related to dispersal behaviour showed differential expression between populations located at each end of the invasive range. These genes could play an important role in the behavioural changes characteristic of range expansion in Australian cane toads.

Captivity induces large and population-dependent brain transcriptomic changes in wild-caught cane toads (Rhinella marina).

Gene expression levels are key molecular phenotypes at the interplay between genotype and environment. Mounting evidence suggests that short-term changes in environmental conditions, such as those encountered in captivity, can substantially affect gene expression levels. Yet, the exact magnitude of this effect, how general it is, and whether it results in parallel changes across populations are not well understood. Here, we take advantage of the well-studied cane toad, Rhinella marina, to examine the effect of short-term captivity on brain gene expression levels, and determine whether effects of captivity differ between long-colonized and vanguard populations of the cane toad's Australian invasion range. We compared the transcriptomes of wild-caught toads immediately assayed with those from toads captured from the same populations but maintained in captivity for seven months. We found large differences in gene expression levels between captive and wild-caught toads from the same population, with an over-representation of processes related to behaviour and the response to stress. Captivity had a much larger effect on both gene expression levels and gene expression variability in toads from vanguard populations compared to toads from long-colonized areas, potentially indicating an increased plasticity in toads at the leading edge of the invasion. Overall, our findings indicate that short-term captivity can induce large and population-specific transcriptomic changes, which has significant implications for studies comparing phenotypic traits of wild-caught organisms from different populations that have been held in captivity.

Integrative Analyses of Transcriptomes to Explore Common Molecular Effects of Antipsychotic Drugs.

There is little understanding of the underlying molecular mechanism(s) involved in the clinical efficacy of antipsychotics for schizophrenia. This study integrated schizophrenia-associated transcriptional perturbations with antipsychotic-induced gene expression profiles to detect potentially relevant therapeutic targets shared by multiple antipsychotics. Human neuronal-like cells (NT2-N) were treated for 24 h with one of the following antipsychotic drugs: amisulpride, aripiprazole, clozapine, risperidone, or vehicle controls. Drug-induced gene expression patterns were compared to schizophrenia-associated transcriptional data in post-mortem brain tissues. Genes regulated by each of four antipsychotic drugs in the reverse direction to schizophrenia were identified as potential therapeutic-relevant genes. A total of 886 genes were reversely expressed between at least one drug treatment (versus vehicle) and schizophrenia (versus healthy control), in which 218 genes were commonly regulated by all four antipsychotic drugs. The most enriched biological pathways include Wnt signaling and action potential regulation. The protein-protein interaction (PPI) networks found two main clusters having schizophrenia expression quantitative trait loci (eQTL) genes such as PDCD10, ANK2, and AKT3, suggesting further investigation on these genes as potential novel treatment targets.

Effects of Psychotropic Drugs on Ribosomal Genes and Protein Synthesis.

Altered protein synthesis has been implicated in the pathophysiology of several neuropsychiatric disorders, particularly schizophrenia. Ribosomes are the machinery responsible for protein synthesis. However, there remains little information on whether current psychotropic drugs affect ribosomes and contribute to their therapeutic effects. We treated human neuronal-like (NT2-N) cells with amisulpride (10 µM), aripiprazole (0.1 µM), clozapine (10 µM), lamotrigine (50 µM), lithium (2.5 mM), quetiapine (50 µM), risperidone (0.1 µM), valproate (0.5 mM) or vehicle control for 24 h. Transcriptomic and gene set enrichment analysis (GSEA) identified that the ribosomal pathway was altered by these drugs. We found that three of the eight drugs tested significantly decreased ribosomal gene expression, whilst one increased it. Most changes were observed in the components of cytosolic ribosomes and not mitochondrial ribosomes. Protein synthesis assays revealed that aripiprazole, clozapine and lithium all decreased protein synthesis. Several currently prescribed psychotropic drugs seem to impact ribosomal gene expression and protein synthesis. This suggests the possibility of using protein synthesis inhibitors as novel therapeutic agents for neuropsychiatric disorders.

Signatures of selection in a recent invasion reveal adaptive divergence in a highly vagile invasive species.

A detailed understanding of population genetics in invasive populations helps us to identify drivers of successful alien introductions. Here, we investigate putative signals of selection in Australian populations of invasive common starlings, Sturnus vulgaris, and seek to understand how these have been influenced by introduction history. We used reduced representation sequencing to determine population structure, and identify Single Nucleotide Polymorphisms (SNPs) that are putatively under selection. We found that since their introduction into Australia, starling populations have become genetically differentiated despite the potential for high levels of dispersal, and that starlings have responded to selective pressures imposed by a wide range of environmental conditions across their geographic range. Isolation by distance appears to have played a strong role in determining genetic substructure across the starling's Australian range. Analyses of candidate SNPs that are putatively under selection indicated that aridity, precipitation and temperature may be important factors driving adaptive variation across the starling's invasive range in Australia. However, we also noted that the historic introduction regime may leave footprints on sites flagged as being under adaptive selection, and encourage critical interpretation of selection analyses in non-native populations.

Transcriptional Modulation of the Hippo Signaling Pathway by Drugs Used to Treat Bipolar Disorder and Schizophrenia.

Recent reports suggest a link between positive regulation of the Hippo pathway with bipolar disorder (BD), and the Hippo pathway is known to interact with multiple other signaling pathways previously associated with BD and other psychiatric disorders. In this study, neuronal-like NT2 cells were treated with amisulpride (10 µM), aripiprazole (0.1 µM), clozapine (10 µM), lamotrigine (50 µM), lithium (2.5 mM), quetiapine (50 µM), risperidone (0.1 µM), valproate (0.5 mM), or vehicle control for 24 h. Genome-wide mRNA expression was quantified and analyzed using gene set enrichment analysis (GSEA), with genes belonging to Hippo, Wnt, Notch, TGF- ß, and Hedgehog retrieved from the KEGG database. Five of the eight drugs downregulated the genes of the Hippo pathway and modulated several genes involved in the interacting pathways. We speculate that the regulation of these genes, especially by aripiprazole, clozapine, and quetiapine, results in a reduction of MAPK and NFκB pro-inflammatory signaling through modulation of Hippo, Wnt, and TGF-ß pathways. We also employed connectivity map analysis to identify compounds that act on these pathways in a similar manner to the known psychiatric drugs. Thirty-six compounds were identified. The presence of antidepressants and antipsychotics validates our approach and reveals possible new targets for drug repurposing.

Transcriptional Effects of Psychoactive Drugs on Genes Involved in Neurogenesis.

Although neurogenesis is affected in several psychiatric diseases, the effects and mechanisms of action of psychoactive drugs on neurogenesis remain unknown and/or controversial. This study aims to evaluate the effects of psychoactive drugs on the expression of genes involved in neurogenesis. Neuronal-like cells (NT2-N) were treated with amisulpride (10 µM), aripiprazole (0.1 µM), clozapine (10 µM), lamotrigine (50 µM), lithium (2.5 mM), quetiapine (50 µM), risperidone (0.1 µM), or valproate (0.5 mM) for 24 h. Genome wide mRNA expression was quantified and analysed using gene set enrichment analysis, with the neurogenesis gene set retrieved from the Gene Ontology database and the Mammalian Adult Neurogenesis Gene Ontology (MANGO) database. Transcription factors that are more likely to regulate these genes were investigated to better understand the biological processes driving neurogenesis. Targeted metabolomics were performed using gas chromatography-mass spectrometry. Six of the eight drugs decreased the expression of genes involved in neurogenesis in both databases. This suggests that acute treatment with these psychoactive drugs negatively regulates the expression of genes involved in neurogenesis in vitro. SOX2 and three of its target genes (CCND1, BMP4, and DKK1) were also decreased after treatment with quetiapine. This can, at least in part, explain the mechanisms by which these drugs decrease neurogenesis at a transcriptional level in vitro. These results were supported by the finding of increased metabolite markers of mature neurons following treatment with most of the drugs tested, suggesting increased proportions of mature relative to immature neurons consistent with reduced neurogenesis.

Understanding sequencing data as compositions: an outlook and review.

Motivation: Although seldom acknowledged explicitly, count data generated by sequencing platforms exist as compositions for which the abundance of each component (e.g. gene or transcript) is only coherently interpretable relative to other components within that sample. This property arises from the assay technology itself, whereby the number of counts recorded for each sample is constrained by an arbitrary total sum (i.e. library size). Consequently, sequencing data, as compositional data, exist in a non-Euclidean space that, without normalization or transformation, renders invalid many conventional analyses, including distance measures, correlation coefficients and multivariate statistical models. Results: The purpose of this review is to summarize the principles of compositional data analysis (CoDA), provide evidence for why sequencing data are compositional, discuss compositionally valid methods available for analyzing sequencing data, and highlight future directions with regard to this field of study. Supplementary information: Supplementary data are available at Bioinformatics online.

Benchmarking differential expression analysis tools for RNA-Seq: normalization-based vs. log-ratio transformation-based methods.

BACKGROUND: Count data generated by next-generation sequencing assays do not measure absolute transcript abundances. Instead, the data are constrained to an arbitrary "library size" by the sequencing depth of the assay, and typically must be normalized prior to statistical analysis. The constrained nature of these data means one could alternatively use a log-ratio transformation in lieu of normalization, as often done when testing for differential abundance (DA) of operational taxonomic units (OTUs) in 16S rRNA data. Therefore, we benchmark how well the ALDEx2 package, a transformation-based DA tool, detects differential expression in high-throughput RNA-sequencing data (RNA-Seq), compared to conventional RNA-Seq methods such as edgeR and DESeq2. RESULTS: To evaluate the performance of log-ratio transformation-based tools, we apply the ALDEx2 package to two simulated, and two real, RNA-Seq data sets. One of the latter was previously used to benchmark dozens of conventional RNA-Seq differential expression methods, enabling us to directly compare transformation-based approaches. We show that ALDEx2, widely used in meta-genomics research, identifies differentially expressed genes (and transcripts) from RNA-Seq data with high precision and, given sufficient sample sizes, high recall too (regardless of the alignment and quantification procedure used). Although we show that the choice in log-ratio transformation can affect performance, ALDEx2 has high precision (i.e., few false positives) across all transformations. Finally, we present a novel, iterative log-ratio transformation (now implemented in ALDEx2) that further improves performance in simulations. CONCLUSIONS: Our results suggest that log-ratio transformation-based methods can work to measure differential expression from RNA-Seq data, provided that certain assumptions are met. Moreover, these methods have very high precision (i.e., few false positives) in simulations and perform well on real data too. With previously demonstrated applicability to 16S rRNA data, ALDEx2 can thus serve as a single tool for data from multiple sequencing modalities.

The Drosophila melanogaster Genetic Reference Panel.

A major challenge of biology is understanding the relationship between molecular genetic variation and variation in quantitative traits, including fitness. This relationship determines our ability to predict phenotypes from genotypes and to understand how evolutionary forces shape variation within and between species. Previous efforts to dissect the genotype-phenotype map were based on incomplete genotypic information. Here, we describe the Drosophila melanogaster Genetic Reference Panel (DGRP), a community resource for analysis of population genomics and quantitative traits. The DGRP consists of fully sequenced inbred lines derived from a natural population. Population genomic analyses reveal reduced polymorphism in centromeric autosomal regions and the X chromosome, evidence for positive and negative selection, and rapid evolution of the X chromosome. Many variants in novel genes, most at low frequency, are associated with quantitative traits and explain a large fraction of the phenotypic variance. The DGRP facilitates genotype-phenotype mapping using the power of Drosophila genetics.

Multiple dispersal vectors drive range expansion in an invasive marine species.

The establishment and subsequent spread of invasive species is widely recognized as one of the most threatening processes contributing to global biodiversity loss. This is especially true for marine and estuarine ecosystems, which have experienced significant increases in the number of invasive species with the increase in global maritime trade. Understanding the rate and mechanisms of range expansion is therefore of significant interest to ecologists and conservation managers alike. Using a combination of population genetic surveys, environmental DNA (eDNA) plankton sampling and hydrodynamic modelling, we examined the patterns of introduction of the predatory Northern Pacific seastar (Asterias amurensis) and pathways of secondary spread within southeast Australia. Genetic surveys across the invasive range reveal some genetic divergence between the two main invasive regions and no evidence of ongoing gene flow, a pattern that is consistent with the establishment of the second invasive region via a human-mediated translocation event. In contrast, hydrodynamic modelling combined with eDNA plankton sampling demonstrated that the establishment of range expansion populations within a region is consistent with natural larval dispersal and recruitment. Our results suggest that both anthropogenic and natural dispersal vectors have played an important role in the range expansion of this species in Australia. The multiple modes of spread combined with high levels of fecundity and a long larval duration in A. amurensis suggests it is likely to continue its range expansion and significantly impact Australian marine ecosystems.

A genetic perspective on rapid evolution in cane toads (Rhinella marina).

The process of biological invasion exposes a species to novel pressures, in terms of both the environments it encounters and the evolutionary consequences of range expansion. Several invaders have been shown to exhibit rapid evolutionary changes in response to those pressures, thus providing robust opportunities to clarify the processes at work during rapid phenotypic transitions. The accelerating pace of invasion of cane toads (Rhinella marina) in tropical Australia during its 80-year history has been well characterized at the phenotypic level, including common-garden experiments that demonstrate heritability of several dispersal-relevant traits. Individuals from the invasion front (and their progeny) show distinctive changes in morphology, physiology and behaviour that, in combination, result in far more rapid dispersal than is true of conspecifics from long-colonized areas. The extensive body of work on cane toad ecology enables us to place into context studies of the genetic basis of these traits. Our analyses of differential gene expression from toads from both ends of this invasion-history transect reveal substantial upregulation of many genes, notably those involved in metabolism and cellular repair. Clearly, then, the dramatically rapid phenotypic evolution of cane toads in Australia has been accompanied by substantial shifts in gene expression, suggesting that this system is well suited to investigating the genetic underpinnings of invasiveness.

Population genomics of the Wolbachia endosymbiont in Drosophila melanogaster.

Wolbachia are maternally inherited symbiotic bacteria, commonly found in arthropods, which are able to manipulate the reproduction of their host in order to maximise their transmission. The evolutionary history of endosymbionts like Wolbachia can be revealed by integrating information on infection status in natural populations with patterns of sequence variation in Wolbachia and host mitochondrial genomes. Here we use whole-genome resequencing data from 290 lines of Drosophila melanogaster from North America, Europe, and Africa to predict Wolbachia infection status, estimate relative cytoplasmic genome copy number, and reconstruct Wolbachia and mitochondrial genome sequences. Overall, 63% of Drosophila strains were predicted to be infected with Wolbachia by our in silico analysis pipeline, which shows 99% concordance with infection status determined by diagnostic PCR. Complete Wolbachia and mitochondrial genomes show congruent phylogenies, consistent with strict vertical transmission through the maternal cytoplasm and imperfect transmission of Wolbachia. Bayesian phylogenetic analysis reveals that the most recent common ancestor of all Wolbachia and mitochondrial genomes in D. melanogaster dates to around 8,000 years ago. We find evidence for a recent global replacement of ancestral Wolbachia and mtDNA lineages, but our data suggest that the derived wMel lineage arose several thousand years ago, not in the 20th century as previously proposed. Our data also provide evidence that this global replacement event is incomplete and is likely to be one of several similar incomplete replacement events that have occurred since the out-of-Africa migration that allowed D. melanogaster to colonize worldwide habitats. This study provides a complete genomic analysis of the evolutionary mode and temporal dynamics of the D. melanogaster-Wolbachia symbiosis, as well as important resources for further analyses of the impact of Wolbachia on host biology.

Is developmental plasticity triggered by DNA methylation changes in the invasive cane toad (Rhinella marina)?

Many organisms can adjust their development according to environmental conditions, including the presence of conspecifics. Although this developmental plasticity is common in amphibians, its underlying molecular mechanisms remain largely unknown. Exposure during development to either 'cannibal cues' from older conspecifics, or 'alarm cues' from injured conspecifics, causes reduced growth and survival in cane toad (Rhinella marina) tadpoles. Epigenetic modifications, such as changes in DNA methylation patterns, are a plausible mechanism underlying these developmental plastic responses. Here we tested this hypothesis, and asked whether cannibal cues and alarm cues trigger the same DNA methylation changes in developing cane toads. We found that exposure to both cannibal cues and alarm cues was associated with local changes in DNA methylation patterns. These DNA methylation changes affected genes putatively involved in developmental processes, but in different genomic regions for different conspecific-derived cues. Genetic background explains most of the epigenetic variation among individuals. Overall, the molecular mechanisms triggered by exposure to cannibal cues seem to differ from those triggered by alarm cues. Studies linking epigenetic modifications to transcriptional activity are needed to clarify the proximate mechanisms that regulate developmental plasticity in cane toads.

Epithelial de-differentiation triggered by co-ordinate epigenetic inactivation of the EHF and CDX1 transcription factors drives colorectal cancer progression.

Colorectal cancers (CRCs) often display histological features indicative of aberrant differentiation but the molecular underpinnings of this trait and whether it directly drives disease progression is unclear. Here, we identify co-ordinate epigenetic inactivation of two epithelial-specific transcription factors, EHF and CDX1, as a mechanism driving differentiation loss in CRCs. Re-expression of EHF and CDX1 in poorly-differentiated CRC cells induced extensive chromatin remodelling, transcriptional re-programming, and differentiation along the enterocytic lineage, leading to reduced growth and metastasis. Strikingly, EHF and CDX1 were also able to reprogramme non-colonic epithelial cells to express colonic differentiation markers. By contrast, inactivation of EHF and CDX1 in well-differentiated CRC cells triggered tumour de-differentiation. Mechanistically, we demonstrate that EHF physically interacts with CDX1 via its PNT domain, and that these transcription factors co-operatively drive transcription of the colonic differentiation marker, VIL1. Compound genetic deletion of Ehf and Cdx1 in the mouse colon disrupted normal colonic differentiation and significantly enhanced colorectal tumour progression. These findings thus reveal a novel mechanism driving epithelial de-differentiation and tumour progression in CRC.

Effects of psychoactive drugs on cellular bioenergetic pathways.

OBJECTIVES: To investigate the actions of lithium, valproate, lamotrigine and quetiapine on bioenergetic pathways in cultured NT2-N neuronal-like cells and C8-B4 microglial cells. METHODS: NT2-N and C8-B4 cells were cultured and treated with lithium (2.5 mM), valproate (0.5 mM), quetiapine (0.05 mM) or lamotrigine (0.05 mM) for 24 hours. Gene expression and the mitochondrial bioenergetic profile were measured in both cell lines. RESULTS: In NT2-N cells, valproate increased oxidative phosphorylation (OXPHOS) gene expression, mitochondrial uncoupling and maximal respiratory capacity, while quetiapine decreased OXPHOS gene expression and respiration linked to ATP turnover, as well as decreasing the expression of genes in the citric acid cycle. Lamotrigine decreased OXPHOS gene expression but had no effect on respiration, while lithium reduced the expression of genes in the citric acid cycle. In C8-B4 cells, valproate and lithium increased OXPHOS gene expression, and valproate increased basal respiratory rate and maximal and spare respiratory capacities. In contrast, quetiapine significantly reduced basal respiratory rate and maximal and spare respiratory capacities. CONCLUSIONS: Overall our data suggest that some drugs used to treat neuropsychiatric and affective disorders have actions on a range of cellular bioenergetic processes, which could impact their effects in patients.

Biological Mechanism(s) Underpinning the Association between Antipsychotic Drugs and Weight Gain.

Weight gain and consequent metabolic alterations are common side-effects of many antipsychotic drugs. Interestingly, several studies have suggested that improvement in symptoms and adverse metabolic effects are correlated. We used next generation sequencing data from NT-2 (human neuronal) cells treated with aripiprazole, amisulpride, risperidone, quetiapine, clozapine, or vehicle control, and compared with the Pillinger P-score (ranked from 0 to 1, indicating greater increase in weight gain and related metabolic parameters) to identify the genes most associated with the drugs' propensity to cause weight gain. The top 500 genes ranked for their correlation with the drugs' propensity to cause weight gain were subjected to pathway analysis using DAVID (NIH). We further investigated transcription factors (TFs) that are more likely to regulate the genes involved in these processes using the prediction tool of key TFs from TRRUST. The results suggest an enrichment for genes involved in lipid biosynthesis and metabolism, which are of interest for mechanisms underpinning weight-gain. The list of genes involved in the lipid pathways that correlated with weight gain was enriched for genes transcriptionally regulated by SREBF1 and SREBF2. Furthermore, quetiapine significantly increased the expression of SREBF1 and SREBF2 in NT-2 cells. Our results suggest that the effects of these antipsychotic drugs on lipid metabolism may be mediated, at least in part, via regulation of SREBF1/SREBF2 expression, with evidence of a direct effect of quetiapine on the expression of SREBF1/2. The effects of antipsychotic drugs on lipid metabolism may influence white matter structure (therapeutic effect) and the risk of weight gain, lipid disturbances, and, consequently, metabolic syndrome (adverse effects). Understanding the different molecular effects of these drugs could inform a personalized medicine approach in treating patients with schizophrenia.

Discovery of Novel Viruses Associated With the Invasive Cane Toad (Rhinella marina) in Its Native and Introduced Ranges.

Cane toads (Rhinella marina) are notoriously successful invaders: from 101 individuals brought to Australia in 1935, poisonous toads now cover an area >1.2 million km2 with adverse effects on native fauna. Despite extensive research on the role of macroparasites in cane toad invasion, viral research is lagging. We compared viral prevalence and diversity between toads in their native range (French Guiana, n=25) and two introduced ranges: Australia (n=151) and Hawai'i (n=10) with a metatranscriptomic and metagenomic approach combined with PCR screening. Australian toads almost exclusively harbor one of seven viruses detected globally. Rhimavirus-A (Picornaviridae) exhibited low genetic diversity and likely actively infected 9% of sampled Australian toads extending across ~2,000km of Northern Australia and up to the current invasion front. In native range cane toads, we identified multiple phylogenetically distinct viruses (Iridoviridae, Picornaviridae, Papillomaviridae, and Nackedna-like virus). None of the same viruses was detected in both ranges, suggesting that Australian cane toads have largely escaped the viral infection experienced by their native range counterparts. The novel native range viruses described here are potential biocontrol agents, as Australian toads likely lack prior immunological exposure to these viruses. Overall, our evidence suggests that there may be differences between viruses infecting cane toads in their native vs. introduced ranges, which lays the groundwork for further studies on how these viruses have influenced the toads' invasion history.

Load carrying during locomotion in the barnacle goose (Branta leucopsis): the effect of load placement and size.

Load carrying has been used to study the energetics and mechanics of locomotion in a range of taxa. Here we investigated the energetic and kinematic effects of trunk and limb loading in walking barnacle geese (Branta leucopsis). A directly proportional relationship between increasing back-load mass and metabolic rate was established, indicating that the barnacle goose can carry back loads (up to 20% of body mass) more economically than the majority of mammals. The increased cost of supporting and propelling the body during locomotion is likely to account for a major proportion of the extra metabolic cost. Sternal loads up to 15% of body mass were approximately twice as expensive to carry as back loads. Given the key role in dorso-ventral movement of the sternum during respiration we suggest that moving this extra mass may account for the elevated metabolic rate. Loading the distal limb with 5% extra mass incurred the greatest proportional rise in metabolism, and also caused increases in stride length, swing duration and stride frequency during locomotion. The increased work required to move the loaded limb may explain the high cost of walking.

Immune and environment-driven gene expression during invasion: An eco-immunological application of RNA-Seq.

Host-pathogen associations change rapidly during a biological invasion and are predicted to impose strong selection on immune function. It has been proposed that the invader may experience an abrupt reduction in pathogen-mediated selection ("enemy release"), thereby favoring decreased investment into "costly" immune responses. Across plants and animals, there is mixed support for this prediction. Pathogens are not the only form of selection imposed on invaders; differences in abiotic environmental conditions between native and introduced ranges are also expected to drive rapid evolution. Here, we use RNA-Seq to assess the expression patterns of immune and environmentally associated genes in the cane toad (Rhinella marina) across its invasive Australian range. Transcripts encoding mediators of costly immune responses (inflammation, cytotoxicity) showed a curvilinear relationship with invasion history, with highest expression in toads from oldest and newest colonized areas. This pattern is surprising given theoretical expectations of density dynamics in invasive species and may be because density influences both intraspecific competition and parasite transmission, generating conflicting effects on the strength of immune responses. Alternatively, this expression pattern may be the result of other evolutionary forces, such as spatial sorting and genetic drift, working simultaneously with natural selection. Our findings do not support predictions about immune function based on the enemy release hypothesis and suggest instead that the effects of enemy release are difficult to isolate in wild populations, especially in the absence of information regarding parasite and pathogen infection. Additionally, expression patterns of genes underlying putatively environmentally associated traits are consistent with previous genetic studies, providing further support that Australian cane toads have adapted to novel abiotic challenges.