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.

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.

Increased Adaptive Variation Despite Reduced Overall Genetic Diversity in a Rapidly Adapting Invader.

Invasive species often evolve rapidly following introduction despite genetic bottlenecks that may result from small numbers of founders; however, some invasions may not fit this "genetic paradox". The invasive cane toad (Rhinella marina) displays high phenotypic variation across its introduced Australian range. Here, we used three genome-wide datasets to characterize their population structure and genetic diversity. We found that toads form three genetic clusters: 1) native range toads, 2) toads from the source population in Hawaii and long-established areas near introduction sites in Australia, and 3) toads from more recently established northern Australian sites. Although we find an overall reduction in genetic diversity following introduction, we do not see this reduction in loci putatively under selection, suggesting that genetic diversity may have been maintained at ecologically relevant traits, or that mutation rates were high enough to maintain adaptive potential. Nonetheless, toads encounter novel environmental challenges in Australia, and the transition between genetic clusters occurs at a point along the invasion transect where temperature rises and rainfall decreases. We identify environmentally associated loci known to be involved in resistance to heat and dehydration. This study highlights that natural selection occurs rapidly and plays a vital role in shaping the structure of invasive populations.

Light pollution increases West Nile virus competence of a ubiquitous passerine reservoir species.

Among the many anthropogenic changes that impact humans and wildlife, one of the most pervasive but least understood is light pollution. Although detrimental physiological and behavioural effects resulting from exposure to light at night are widely appreciated, the impacts of light pollution on infectious disease risk have not been studied. Here, we demonstrate that artificial light at night (ALAN) extends the infectious-to-vector period of the house sparrow (Passer domesticus), an urban-dwelling avian reservoir host of West Nile virus (WNV). Sparrows exposed to ALAN maintained transmissible viral titres for 2 days longer than controls but did not experience greater WNV-induced mortality during this window. Transcriptionally, ALAN altered the expression of gene regulatory networks including key hubs (OASL, PLBD1 and TRAP1) and effector genes known to affect WNV dissemination (SOCS). Despite mounting anti-viral immune responses earlier, transcriptomic signatures indicated that ALAN-exposed individuals probably experienced pathogen-induced damage and immunopathology, potentially due to evasion of immune effectors. A simple mathematical modelling exercise indicated that ALAN-induced increases of host infectious-to-vector period could increase WNV outbreak potential by approximately 41%. ALAN probably affects other host and vector traits relevant to transmission, and additional research is needed to advise the management of zoonotic diseases in light-polluted areas.

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.

Viral Discovery in the Invasive Australian Cane Toad (Rhinella marina) Using Metatranscriptomic and Genomic Approaches.

Cane toads are a notorious invasive species, inhabiting over 1.2 million km2 of Australia and threatening native biodiversity. The release of pathogenic cane toad viruses is one possible biocontrol strategy yet is currently hindered by the poorly described cane toad virome. Metatranscriptomic analysis of 16 cane toad livers revealed the presence of a novel and full-length picornavirus, Rhimavirus A (RhiV-A), a member of a reptile- and amphibian-specific cluster of the Picornaviridae basal to the Kobuvirus-like group. In the combined liver transcriptome, we also identified a complete genome sequence of a distinct epsilonretrovirus, Rhinella marina endogenous retrovirus (RMERV). The recently sequenced cane toad genome contains 8 complete RMERV proviruses as well as 21 additional truncated insertions. The oldest full-length RMERV provirus was estimated to have inserted 1.9 million years ago (MYA). To screen for these viral sequences in additional toads, we analyzed publicly available transcriptomes from six diverse Australian locations. RhiV-A transcripts were identified in toads sampled from three locations across 1,000 km of Australia, stretching to the current Western Australia (WA) invasion front, while RMERV transcripts were observed at all six sites. Finally, we scanned the cane toad genome for nonretroviral endogenous viral elements, finding three sequences related to small DNA viruses in the family Circoviridae This shows ancestral circoviral infection with subsequent genomic integration. The identification of these current and past viral infections enriches our knowledge of the cane toad virome, an understanding of which will facilitate future work on infection and disease in this important invasive species.IMPORTANCE Cane toads are poisonous amphibians that were introduced to Australia in 1935 for insect control. Since then, their population has increased dramatically, and they now threaten many native Australian species. One potential method to control the population is to release a cane toad virus with high mortality rates, yet few cane toad viruses have been characterized. This study samples cane toads from different Australian locations and uses an RNA sequencing and computational approach to find new viruses. We report novel complete picornavirus and retrovirus sequences that were genetically similar to viruses infecting frogs, reptiles, and fish. Using data generated in other studies, we show that these viral sequences are present in cane toads from distinct Australian locations. Three sequences related to circoviruses were also found in the toad genome. The identification of new viral sequences will aid future studies that investigate their prevalence and potential as agents for biocontrol.

Improving amphibian genomic resources: a multitissue reference transcriptome of an iconic invader.

Background: Cane toads (Rhinella marina) are an iconic invasive species introduced to 4 continents and well utilized for studies of rapid evolution in introduced environments. Despite the long introduction history of this species, its profound ecological impacts, and its utility for demonstrating evolutionary principles, genetic information is sparse. Here we produce a de novo transcriptome spanning multiple tissues and life stages to enable investigation of the genetic basis of previously identified rapid phenotypic change over the introduced range. Findings: Using approximately 1.9 billion reads from developing tadpoles and 6 adult tissue-specific cDNA libraries, as well as a transcriptome assembly pipeline encompassing 100 separate de novo assemblies, we constructed 62 202 transcripts, of which we functionally annotated ∼50%. Our transcriptome assembly exhibits 90% full-length completeness of the Benchmarking Universal Single-Copy Orthologs data set. Robust assembly metrics and comparisons with several available anuran transcriptomes and genomes indicate that our cane toad assembly is one of the most complete anuran genomic resources available. Conclusions: This comprehensive anuran transcriptome will provide a valuable resource for investigation of genes under selection during invasion in cane toads, but will also greatly expand our general knowledge of anuran genomes, which are underrepresented in the literature. The data set is publically available in NCBI and GigaDB to serve as a resource for other researchers.

Effects of invasion history on physiological responses to immune system activation in invasive Australian cane toads.

The cane toad (Rhinella marina) has undergone rapid evolution during its invasion of tropical Australia. Toads from invasion front populations (in Western Australia) have been reported to exhibit a stronger baseline phagocytic immune response than do conspecifics from range core populations (in Queensland). To explore this difference, we injected wild-caught toads from both areas with the experimental antigen lipopolysaccharide (LPS, to mimic bacterial infection) and measured whole-blood phagocytosis. Because the hypothalamic-pituitary-adrenal axis is stimulated by infection (and may influence immune responses), we measured glucocorticoid response through urinary corticosterone levels. Relative to injection of a control (phosphate-buffered saline), LPS injection increased both phagocytosis and the proportion of neutrophils in the blood. However, responses were similar in toads from both populations. This null result may reflect the ubiquity of bacterial risks across the toad's invaded range; utilization of this immune pathway may not have altered during the process of invasion. LPS injection also induced a reduction in urinary corticosterone levels, perhaps as a result of chronic stress.

The ribosome-associated complex antagonizes prion formation in yeast.

The number of known fungal proteins capable of switching between alternative stable conformations is steadily increasing, suggesting that a prion-like mechanism may be broadly utilized as a means to propagate altered cellular states. To gain insight into the mechanisms by which cells regulate prion formation and toxicity we examined the role of the yeast ribosome-associated complex (RAC) in modulating both the formation of the [PSI(+)] prion - an alternative conformer of Sup35 protein - and the toxicity of aggregation-prone polypeptides. The Hsp40 RAC chaperone Zuo1 anchors the RAC to ribosomes and stimulates the ATPase activity of the Hsp70 chaperone Ssb. We found that cells lacking Zuo1 are sensitive to over-expression of some aggregation-prone proteins, including the Sup35 prion domain, suggesting that co-translational protein misfolding increases in Δzuo1 strains. Consistent with this finding, Δzuo1 cells exhibit higher frequencies of spontaneous and induced prion formation. Cells expressing mutant forms of Zuo1 lacking either a C-terminal charged region required for ribosome association, or the J-domain responsible for Ssb ATPase stimulation, exhibit similarly high frequencies of prion formation. Our findings are consistent with a role for the RAC in chaperoning nascent Sup35 to regulate folding of the N-terminal prion domain as it emerges from the ribosome.