Structural and inhibitor sensitivity analysis of influenza B-like viral neuraminidases derived from Asiatic toad and spiny eel.

Influenza virus neuraminidase (NA) is an important target for antiviral development because it plays a crucial role in releasing newly assembled viruses. Two unique influenza-like virus genomes were recently reported in the Wuhan Asiatic toad and Wuhan spiny eel. Their NA genes appear to be highly divergent from all known influenza NAs, raising key questions as to whether the Asiatic toad influenza-like virus NA (tNA) and spiny eel NA (eNA) have canonical NA activities and structures and whether they show sensitivity to NA inhibitors (NAIs). Here, we found that both tNA and eNA have neuraminidase activities. A detailed structural analysis revealed that tNA and eNA present similar overall structures to currently known NAs, with a conserved calcium binding site. Inhibition assays indicated that tNA is resistant to NAIs, while eNA is still sensitive to NAIs. E119 is conserved in canonical NAs. The P119E substitution in tNA can restore sensitivity to NAIs, and, in contrast, the E119P substitution in eNA decreased its sensitivity to NAIs. The structures of NA-inhibitor complexes further provide a detailed insight into NA-inhibitor interactions at the atomic level. Moreover, tNA and eNA have unique N-glycosylation sites compared with canonical NAs. Collectively, the structural features, NA activities, and sensitivities to NAIs suggest that fish- and amphibian-derived influenza-like viruses may circulate in these vertebrates. More attention should be paid to these influenza-like viruses because their NA molecules may play roles in the emergence of NAI resistance.

A cane toad (Rhinella marina) N-methyltransferase converts primary indolethylamines to tertiary psychedelic amines.

Psychedelic indolethylamines have emerged as potential medicines to treat several psychiatric pathologies. Natural sources of these compounds include 'magic mushrooms' (Psilocybe spp.), plants used to prepare ayahuasca, and toads. The skin and parotid glands of certain toads accumulate a variety of specialized metabolites including toxic guanidine alkaloids, lipophilic alkaloids, poisonous steroids, and hallucinogenic indolethylamines such as DMT, 5-methoxy-DMT, and bufotenin. The occurrence of psychedelics has contributed to the ceremonial use of toads, particularly among Mesoamerican peoples. Yet, the biosynthesis of psychedelic alkaloids has not been elucidated. Herein, we report a novel indolethylamine N-methyltransferase (RmNMT) from cane toad (Rhinella marina). The RmNMT sequence was used to identify a related NMT from the common toad, Bufo bufo. Close homologs from various frog species were inactive, suggesting a role for psychedelic indolethylamine biosynthesis in toads. Enzyme kinetic analyses and comparison with functionally similar enzymes showed that recombinant RmNMT was an effective catalyst and not product inhibited. The substrate promiscuity of RmNMT enabled the bioproduction of a variety of substituted indolethylamines at levels sufficient for purification, pharmacological screening, and metabolic stability assays. Since the therapeutic potential of psychedelics has been linked to activity at serotonergic receptors, we evaluated binding of derivatives at 5-HT1A and 5-HT2A receptors. Primary amines exhibited enhanced affinity at the 5-HT1A receptor compared with tertiary amines. With the exception of 6-substituted derivatives, N,N-dimethylation also protected against catabolism by liver microsomes.

Isolation, Identification and Chemical Modification of Bufadienolides from Bufo melanostictus Schneider and Their Cytotoxic Activities against Prostate Cancer Cells.

The traditional Chinese medicine toad venom (Venenum bufonis) has been extensively used to treat various diseases, including cancers, in China and other Southeast Asian countries. The major constituents of toad venom, e.g., bufadienolides and alkaloids, exhibit broad-spectrum pharmacological effects in cancers. Herein, two new bufadienolides (1 and 2), along with eleven known compounds (3-13) were successfully isolated from Bufo melanostictus Schneider. Their structures were elucidated by extensive spectroscopic data and X-ray diffraction analysis. Furthermore, four lactam derivatives were synthesized through the transformation of bufadienolides lactones. The inhibitory effects of these compounds against human prostate cancer cell lines PC-3 and DU145 were evaluated. The outcomes indicated a notable trend, with a substantial subset displaying nanomolar range IC50 values against PC-3 and DU145 cells, underscoring their pronounced cytotoxicity. Moreover, a noteworthy distinction surfaces, wherein lactones consistently outperformed their lactam counterparts, further validating their heightened potency for the treatment of prostate cancer. This study contributes significant preclinical evidence substantiating the therapeutic viability of bufadienolides and toad venom as intervention strategies for prostate cancer.

Climate Change Alters Sexual Signaling in a Desert-Adapted Frog.

AbstractClimate change is altering species' habitats, phenology, and behavior. Although sexual behaviors impact population persistence and fitness, climate change's effects on sexual signals are understudied. Climate change can directly alter temperature-dependent sexual signals, cause changes in body size or condition that affect signal production, or alter the selective landscape of sexual signals. We tested whether temperature-dependent mating calls of Mexican spadefoot toads (Spea multiplicata) had changed in concert with climate in the southwestern United States across 22 years. We document increasing air temperatures, decreasing rainfall, and changing seasonal patterns of temperature and rainfall in the spadefoots' habitat. Despite increasing air temperatures, spadefoots' ephemeral breeding ponds have been getting colder at most elevations, and male calls have been slowing as a result. However, temperature-standardized call characters have become faster, and male condition has increased, possibly due to changes in the selective environment. Thus, climate change might generate rapid, complex changes in sexual signals with important evolutionary consequences.

"In vivo biosynthesis of N,N-dimethyltryptamine, 5-MeO-N,N-dimethyltryptamine, and bufotenine in E.coli".

N,N-dimethyltryptamine (DMT), 5-methoxy-N,N-dimethyltryptamine (5-MeO-DMT) and 5-hydroxy-N,N-dimethyltryptamine (bufotenine) are psychedelic tryptamines found naturally in both plants and animals and have shown clinical potential to help treat mental disorders, such as anxiety and depression. Advances in both metabolic and genetic engineering make it possible to engineer microbes as cell factories to produce DMT and its aforementioned derivatives to meet demand for ongoing clinical study. Here, we present the development of a biosynthetic production pathway for DMT, 5-MeO-DMT, and bufotenine in the model microbe Escherichia coli. Through the application of genetic optimization techniques and process optimization in benchtop fermenters, the in vivo production of DMT in E. coli was observed. DMT production with tryptophan supplementation reached maximum titers of 74.7 ± 10.5 mg/L under fed batch conditions in a 2-L bioreactor. Additionally, we show the first reported case of de novo production of DMT (from glucose) in E. coli at a maximum titer of 14.0 mg/L and report the first example of microbial 5-MeO-DMT and bufotenine production in vivo. This work provides a starting point for further genetic and fermentation optimization studies with the goal to increase methylated tryptamine production metrics to industrially competitive levels.

Easternmost distribution of Bufo bufo (Linnaeus, 1758) in Türkiye: implications for the putative contact zone between B. bufo and B. verrucosissimus.

The geographic range of a species is crucial for obtaining information on the exact distribution of the species. The geographic data are important for delimiting distinct species or exploring the degree of differentiation among different populations of a species. The local details of species boundaries facilitate the study of the importance of phylogeographic background, secondary contacts, and hybrid zones, along with the relations between the species and its extrinsic environmental factors. In the present study, the range boundaries of Bufo bufo and Bufo verrucosissimus in the north-eastern region of Türkiye were delineated using an integrative taxonomic approach that utilized a combination of molecular and morphological data. According to the mtDNA results of the present study, B. bufo inhabits a single distribution from Iyidere town to Çayeli town in Rize, while B. verrucosissimus is distributed from Savsat town of Artvin to Ardesen town in Rize. In addition, the two species coexist in Pazar, Hemsin, and Çamlihemsin towns in Rize. The demographic analyses indicated a distinct population expansion for the B. verrucosissimus species after the Last Glacial Maximum, while the same did not occur for B. bufo. The univariate and multivariate statistical analyses conducted for the morphological data of the two species corroborated the presence of a putative contact zone between B. bufo and B. verrucosissimus. In summary, the present study resolved the non-distinct geographic boundaries between B. bufo and B. verrucosissimus species and also revealed the easternmost distribution of B. bufo in Türkiye. In addition, important evidence on the putative contact zone between the two species was indicated using an integrative taxonomic approach.

Proteins from toad's parotoid macroglands: do they play a role in gland functioning and chemical defence?

BACKGROUND: Parotoid gland secretion of bufonid toads is a rich source of toxic molecules that are used against predators, parasites and pathogens. Bufadienolides and biogenic amines are the principal compounds responsible for toxicity of parotoid secretion. Many toxicological and pharmacological analyses of parotoid secretions have been performed, but little is known about the processes related to poison production and secretion. Therefore, our aim was to investigate protein content in parotoids of the common toad, Bufo bufo, to understand the processes that regulate synthesis and excretion of toxins as well as functioning of parotoid macroglands. RESULTS: Applying a proteomic approach we identified 162 proteins in the extract from toad's parotoids that were classified into 11 categories of biological functions. One-third (34.6%) of the identified molecules, including acyl-CoA-binding protein, actin, catalase, calmodulin, and enolases, were involved in cell metabolism. We found many proteins related to cell division and cell cycle regulation (12.0%; e.g. histone and tubulin), cell structure maintenance (8.4%; e.g. thymosin beta-4, tubulin), intra- and extracellular transport (8.4%), cell aging and apoptosis (7.3%; e.g. catalase and pyruvate kinase) as well as immune (7.0%; e.g. interleukin-24 and UV excision repair protein) and stress (6.3%; including heat shock proteins, peroxiredoxin-6 and superoxide dismutase) response. We also identified two proteins, phosphomevalonate kinase and isopentenyl-diphosphate delta-isomerase 1, that are involved in synthesis of cholesterol which is a precursor for bufadienolides biosynthesis. STRING protein-protein interaction network predicted for identified proteins showed that most proteins are related to metabolic processes, particularly glycolysis, stress response and DNA repair and replication. The results of GO enrichment and KEGG analyses are also consistent with these findings. CONCLUSION: This finding indicates that cholesterol may be synthesized in parotoids, and not only in the liver from which is then transferred through the bloodstream to the parotoid macroglands. Presence of proteins that regulate cell cycle, cell division, aging and apoptosis may indicate a high epithelial cell turnover in parotoids. Proteins protecting skin cells from DNA damage may help to minimize the harmful effects of UV radiation. Thus, our work extends our knowledge with new and important functions of parotoids, major glands involved in the bufonid chemical defence.

Increasingly cautious sampling, not the black colouration of unpalatable prey, is used by fish in avoidance learning.

The efficiency of aposematic colouration of prey is based on the innate bias or facilitation of avoidance learning of predators. In many toxic bufonids, larvae are uniformly black, which is considered a warning signal. We compared fish predation on normal (black) and 'transient albino' (greyish) common toad Bufo bufo tadpoles that did not differ in toxicity or activity. In a two-stage experiment, each fish was presented with tadpoles of one colour in the first trial and the other colour in a subsequent trial. While tadpoles sampled by fish were typically not ingested, some died from injuries. The attack rate did not differ between tadpole phenotypes nor trials, irrespective of which phenotype was the first exposed to the fish. However, during the second trial, the sampled tadpoles, independent of colouration, were mouthed by fish for shorter periods and tadpole mortality decreased. The duration of mouthing also declined with an increasing number of attacks during subsequent trials. We conclude that in single-species prey populations, black tadpole colouration is not a warning signal as it does not accelerate predator learning about prey unprofitability. Our results indicate that with growing experience, predators sample potentially toxic prey more cautiously. This may explain why natural selection does not eliminate aposematic morphs even if predators continuously sample conspicuous prey.

Environmental and Anthropogenic Factors Shape the Skin Bacterial Communities of a Semi-Arid Amphibian Species.

The amphibian skin microbiome is important in maintaining host health, but is vulnerable to perturbation from changes in biotic and abiotic conditions. Anthropogenic habitat disturbance and emerging infectious diseases are both potential disrupters of the skin microbiome, in addition to being major drivers of amphibian decline globally. We investigated how host environment (hydrology, habitat disturbance), pathogen presence, and host biology (life stage) impact the skin microbiome of wild Dhofar toads (Duttaphrynus dhufarensis) in Oman. We detected ranavirus (but not Batrachochytrium dendrobatidis) across all sampling sites, constituting the first report of this pathogen in Oman, with reduced prevalence in disturbed sites. We show that skin microbiome beta diversity is driven by host life stage, water source, and habitat disturbance, but not ranavirus infection. Finally, although trends in bacterial diversity and differential abundance were evident in disturbed versus undisturbed sites, bacterial co-occurrence patterns determined through network analyses revealed high site specificity. Our results therefore provide support for amphibian skin microbiome diversity and taxa abundance being associated with habitat disturbance, with bacterial co-occurrence (and likely broader aspects of microbial community ecology) being largely site specific.

DO Toxic Invasive Prey Become a Toxin Source for Native Consumers?

Toxic organisms can become food that potentially harms consumers. When these organisms become invasive species, the harm often turns to a serious threat that disrupts native ecosystems. On the other hand, there are consumers that can exploit toxic organisms for food and sequester intact toxins from them for the consumers' own chemical defense. Therefore, it can be expected that toxic invasive prey can become a toxin source for native consumers. Here, we focused on the relationship between toads, which are one of the major toxic invasive organisms and possess bufadienolides (BDs), and Rhabdophis snakes, which sequester BDs from toads. On Sado Island, Japan, R. tigrinus is native, but no toads had inhabited this island until Bufo japonicus formosus was introduced as a domestic invasive species in 1963 and 1964. At present, invasive toads are distributed only in the southwestern part of the island. We collected a total of 25 and 24 R. tigrinus from areas allopatric and sympatric with toads, respectively. Then, we investigated the possession of BDs and the BD profile of these snakes. We found that only R. tigrinus sympatric with toads possessed BDs, whereas all snakes allopatric with toads lacked BDs. Based on the characteristics of the BD profile, the toxin source was identified as B. j. formosus. Our findings show that a new case of impact caused by toxic invasive species, i.e., "toxin supply to native consumers from invasive prey", could occur.

Toxicity and effects on anuran tadpole metamorphosis of the anthranilic diamide insecticides chlorantraniliprole and cyantraniliprole.

The present study examined the acute and chronic toxicity attributed to commercial formulations of the anthranilic diamide insecticides chlorantraniliprole (CHLO) and cyantraniliprole (CYAN) on the neotropical amphibian species Rhinella arenarum, Rhinella fernandezae and Scinax granulatus. The median lethal concentrations obtained after 96 hr exposure (96 hr-LC50) were generally greater than 100 mg/L, except for stage 25 S. Granulatus, which were the most sensitive animals tested with a 96 hr-LC50 value of 46.78 mg/L. In subchronic exposures of R. arenarum, the 21day-LC50 were 151.4 mg/L for CHLO and >160 mg/L for CYAN, the weight gain of the tadpoles during this period not being markedly affected in both cases. Finally, when tadpoles of R. arenarum were exposed to CHLO throughout the metamorphic process, an inverted U-shaped non-monotonic dose-response relationship was observed between exposure concentrations and both % of individuals transiting between stage 39 and 42 and the time required to accomplish this. Data obtained raise the hypothesis of an effect of CHLO on the hypothalamic-pituitary-thyroid (HPT) axis, either directly or through an interaction with the stress-hormone system, as metamorphic progression from stage 39 to S42 occurs under the strict control of thyroid hormones. These observations are important as the anthranilic diamide insecticides are not currently known as endocrine disruptors. Further investigations are needed to clarify the pathways leading to these effects and examine whether environmentally-relevant aquatic concentrations of anthranilic diamides might be impacting amphibian populations in the wild.

Neoplasia in an assurance population of Puerto Rican crested toads (Peltophryne lemur).

Disease monitoring of amphibian assurance populations is an important buffer against ongoing global extinctions. This study documents a high incidence of neoplasia in a zoo-managed assurance population of Puerto Rican crested toads (Peltophryne lemur; PRCTs). Over 5 years, neoplasia was diagnosed in 17/49 (35%) submitted adult PRCTs and was the cause of death or euthanasia in 13/17 (72%). Most toads were male (16/17; 94%) and 6 to 11-years-old (average 8.1 years). Notably, seven toads (41%) had multiple neoplasms. Of the 29 neoplasms identified, 17 (59%) were cutaneous or subcutaneous. The most common neoplasms included mast cell tumors (MCTs; 8/29; 28%), histiocytic sarcomas (6/29; 21%), lymphoma/leukemia (4/29; 14%), and squamous cell carcinomas (3/29; 10%). Distant metastases were documented in 6/8 (75%) toads with MCTs. Causes for neoplasia in this population were not determined though may include genetic or environmental factors. Continued investigations of managed endangered amphibians will help elucidate mechanisms of carcinogenesis.

Identification of Protein Quality Markers in Toad Venom from Bufo gargarizans.

Toad venom is a traditional Chinese medicine with high medicinal value. The existing quality evaluation standards of toad venom have obvious limitations because of the lack of research on proteins. Thus, it is necessary to screen suitable quality markers and establish appropriate quality evaluation methods for toad venom proteins to guarantee their safety and efficacy in clinical applications. SDS-PAGE, HPLC, and cytotoxicity assays were used to analyze differences in protein components of toad venom from different areas. Functional proteins were screened as potential quality markers by proteomic and bioinformatic analyses. The protein components and small molecular components of toad venom were not correlated in content. Additionally, the protein component had strong cytotoxicity. Proteomics analysis showed that 13 antimicrobial proteins, four anti-inflammatory and analgesic proteins, and 20 antitumor proteins were differentially expressed extracellular proteins. A candidate list of functional proteins was coded as potential quality markers. Moreover, Lysozyme C-1, which has antimicrobial activity, and Neuropeptide B (NPB), which has anti-inflammatory and analgesic activity, were identified as potential quality markers for toad venom proteins. Quality markers can be used as the basis of quality studies of toad venom proteins and help to construct and improve safe, scientific, and comprehensive quality evaluation methods.

Application of a novel multicomponent nanoemulsion to tumor therapy Based on the theory of "unification of drugs and excipients".

Toad skin has many pharmacological activities and bufadienolides are regarded as its main anti-tumor components. The poor water solubility, high toxicity, rapid elimination and less selectivity in vivo of bufadienolides limit the application of toad skin. Based on the "unification of drugs and excipients" theory, the toad skin extracts (TSE) and Brucea javanica oil (BJO) nanoemulsions (NEs) were designed to solve the aforementioned problems. BJO as the main oil phase was not only used to prepare the NEs, but played a synergistic therapeutic role combined with TSE. TSE-BJO NEs showed 155 nm particle size, entrapment efficiency of >95% and good stability. TSE-BJO NEs demonstrated superior anti-tumor activity compared with the TSE or BJO NEs alone. The mechanism of TSE-BJO NEs to enhance the antineoplastic efficacy involved several pathways, such as inhibiting cell proliferation, inducing tumor cell apoptosis >40% and arresting cell cycle at G2/M. TSE-BJO NEs could co-deliver drugs into the target cells efficiently and exhibit satisfying synergism. Besides, TSE-BJO NEs facilitated the long circulation of bufadienolides contributing to the high accumulation of drugs at tumor sites and the improvement of anti-tumor efficacy. The study achieves the combinative administration of the toxic TSE and BJO with high efficacy and safety.

Disentangling the developmental origins of a novel phenotype: enhancement versus reversal of environmentally induced gene expression.

Increasing evidence suggests that many novel traits might have originated via plasticity-led evolution (PLE). Yet, little is known of the developmental processes that underpin PLE, especially in its early stages. One such process is 'phenotypic accommodation', which occurs when, in response to a change in the environment, an organism experiences adjustments across variable parts of its phenotype that improve its fitness. Here, we asked if environmentally induced changes in gene expression are enhanced or reversed during phenotypic accommodation of a novel, complex phenotype in spadefoot toad tadpoles (Spea multiplicata). More genes than expected were affected by both the environment and phenotypic accommodation in the liver and brain. However, although phenotypic accommodation primarily reversed environmentally induced changes in gene expression in liver tissue, it enhanced these changes in brain tissue. Thus, depending on the tissue, phenotypic accommodation may either minimize functional disruption via reversal of gene expression patterns or promote novelty via enhancement of existing expression patterns. Our study thereby provides insights into the developmental origins of a novel phenotype and the incipient stages of PLE.

Tunable Toxicity of Bufadienolides is Regulated through a Configuration Inversion Catalyzed by a Short-Chain Dehydrogenase/Reductase.

Bufadienolides are toxic components widely found in amphibious toads that exhibit a wide range of biological activities. Guided by UPLC-QTOF-MS analysis, several 3-epi-bufadienolides with unique structures were isolated from the bile of the Asiatic toad, Bufo gargarizans. However, the enzymatic machinery of this epimerization in toads and its significance in chemical ecology remains poorly understood. Herein, we firstly compared the toxicities of two typical bufadienolides, bufalin (featuring a 14ß-hydroxyl) and resibufogenin (containing a 14, 15-epoxy group), with their corresponding 3-epi isomers in a zebrafish model. The results of the toxicology assays showed that the ratio of maximum non-toxic concentrations of these two pairs of compounds are 256 and 96 times, respectively, thereby indicating that 3-hydroxyl epimerization leads to a significant decrease in toxicity. Aiming to investigate the biotransformation of 3-epi bufadienolides in toads, we applied liver lysate to transform bufalin and found that it could stereoselectively catalyze the conversion of bufalin into its 3α-hydroxyl epimer. Following this, we cloned and characterized a short-chain dehydrogenase/reductase, HSE-1, from the toad liver cDNA library and verified its 3(ß→α)-hydroxysteroid epimerization activity. To the best of our knowledge, this is the first hydroxyl epimerase identified from amphibians that regulates the toxicity of animal-derived natural products.

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.

Population genomics of a predatory mammal reveals patterns of decline and impacts of exposure to toxic toads.

Mammal declines across northern Australia are one of the major biodiversity loss events occurring globally. There has been no regional assessment of the implications of these species declines for genomic diversity. To address this, we conducted a species-wide assessment of genomic diversity in the northern quoll (Dasyurus hallucatus), an Endangered marsupial carnivore. We used next generation sequencing methods to genotype 10,191 single nucleotide polymorphisms (SNPs) in 352 individuals from across a 3220-km length of the continent, investigating patterns of population genomic structure and diversity, and identifying loci showing signals of putative selection. We found strong heterogeneity in the distribution of genomic diversity across the continent, characterized by (i) biogeographical barriers driving hierarchical population structure through long-term isolation, and (ii) severe reductions in diversity resulting from population declines, exacerbated by the spread of introduced toxic cane toads (Rhinella marina). These results warn of a large ongoing loss of genomic diversity and associated adaptive capacity as mammals decline across northern Australia. Encouragingly, populations of the northern quoll established on toad-free islands by translocations appear to have maintained most of the initial genomic diversity after 16 years. By mapping patterns of genomic diversity within and among populations, and investigating these patterns in the context of population declines, we can provide conservation managers with data critical to informed decision-making. This includes the identification of populations that are candidates for genetic management, the importance of remnant island and insurance/translocated populations for the conservation of genetic diversity, and the characterization of putative evolutionarily significant units.

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.

Evidence that catecholaminergic systems mediate dynamic colour change during explosive breeding events in toads.

Many animals communicate by rapidly (within minutes or seconds) changing their body coloration; however, we know little about the physiology of this behaviour. Here we study how catecholaminergic hormones regulate rapid colour change in explosive breeding toads (Duttaphrynus melanostictus), where large groups of males gather and quickly change their colour from brown to bright yellow during reproduction. We find that both epinephrine (EP) and/or norepinephrine (NE) cause the toads' skin to become yellow in minutes, even in the absence of social and environmental cues associated with explosive breeding. We hypothesize that natural selection drives the evolution of rapid colour change by co-opting the functional effects of catecholaminergic action. If so, then hormones involved in 'fight or flight' responses may mechanistically facilitate the emergence of dynamic visual signals that mediate communication in a sexual context.