Behavioural responses of an Australian colubrid snake (Dendrelaphis punctulatus) to a novel toxic prey item (the Cane Toad Rhinella marina)

Whole snake genomes from eighteen families of snakes (Serpentes: Caenophidia) and their applications to systematics

We present genome assemblies for 18 snake species representing 18 families (Serpentes: Caenophidia): Acrochordus granulatus, Aparallactus werneri, Boaedon fuliginosus, Calamaria suluensis, Cerberus rynchops, Grayia smithii, Imantodes cenchoa, Mimophis mahfalensis, Oxyrhabdium leporinum, Pareas carinatus, Psammodynastes pulverulentus, Pseudoxenodon macrops, Pseudoxyrhopus heterurus, Sibynophis collaris, Stegonotus admiraltiensis, Toxicocalamus goodenoughensis, Trimeresurus albolabris, and Tropidonophis doriae. From these new genome assemblies, we extracted thousands of loci commonly used in systematic and phylogenomic studies on snakes, including target-capture datasets composed of ultraconserved elements (UCEs) and anchored hybrid enriched loci (AHEs), as well as traditional Sanger loci. Phylogenies inferred from the two target-capture loci datasets were identical with each other and strongly congruent with previously published snake phylogenies. To show the additional utility of these non-model genomes for investigative evolutionary research, we mined the genome assemblies of two New Guinea island endemics in our dataset (S. admiraltiensis and T. doriae) for the ATP1a3 gene, a thoroughly researched indicator of resistance to toad toxin ingestion by squamates. We find that both these snakes possess the genotype for toad toxin resistance despite their endemism to New Guinea, a region absent of any toads until the human-mediated introduction of Cane Toads in the 1930s. These species possess identical substitutions that suggest the same bufotoxin resistance as their Australian congenerics (Stegonotus australis and Tropidonophis mairii) which forage on invasive Cane Toads. Herein, we show the utility of short-read high-coverage genomes, as well as improving the deficit of available squamate genomes with associated voucher specimens.

Variation in size and shape of toxin glands among cane toads from native-range and invasive populations

If optimal investment in anti-predator defences depends on predation risk, invading new regions (and thus, encountering different predators) may favour shifts in that investment. Cane toads offer an ideal system to test this prediction: expensive anti-predator toxins are stored mainly in parotoid glands whose dimensions are easy to measure, and toad invasions have changed the suites of predators they encounter. Although plasticity may influence parotoid morphology, comparisons between parents and progeny revealed that gland dimensions were highly heritable. That heritability supports the plausibility of an evolved basis to variation in gland dimensions. Measurements of 3779 adult toads show that females have larger glands than males, invasive populations have larger glands than in the native-range, and that parotoid sexual size dimorphism varies strongly among invaded areas. Geographic variation in parotoid morphology may be driven by predation risk to both adult toads and offspring (provisioned with toxins by their mother), with toxins allocated to eggs exacerbating the risk of cannibalism but reducing the risk of interspecific predation. Investment into chemical defences has evolved rapidly during the cane toad’s international diaspora, consistent with the hypothesis that organisms flexibly adjust resource allocation to anti-predator tactics in response to novel challenges.