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.

Ladies in armor: A micro-computed tomographic study of skin calcification in European toads (genus Bufo)

Micro-computed tomography is a powerful tool toward the detailed reconstruction of internal and external morphology, in particular for ossified and other dense tissues. Here, we document and compare the level of calcification in the skin of the head and the parotoids (the external skin glands) in males and females of common and spined toads, Bufo bufo and B. spinosus. In some anurans, including Bufo species, a specific acellular calcified tissue layer within the dermis has been documented (the Eberth-Katschenko, or EK-layer). By a combination of micro-computed tomography and classical histology, we detected additional calcium deposits located in the dermal layer stratum spongiosum, positioned above the EK-layer. We showed that the level of calcification and the presence of additional calcium deposits are size and sex related, increasing in the order B. bufo males, B. spinosus males, B. bufo females to B. spinosus females. The last of these groups is the least variable. Bufo spinosus females have dense calcium deposits in the parotoids and the dorsal and ventral skin. Three-dimensional volume renderings and cross-sectional slices obtained by micro-CT scanning indicate that this approach is a promising technique for further studies on bufonid skin anatomy and geographic variation in skin calcification.

The genome sequence of the common toad, Bufo bufo (Linnaeus, 1758)

We present a genome assembly from an individual male Bufo bufo (the common toad; Chordata; Amphibia; Anura; Bufonidae). The genome sequence is 5.04 gigabases in span. The majority of the assembly (99.1%) is scaffolded into 11 chromosomal pseudomolecules. Gene annotation of this assembly by the NCBI Eukaryotic Genome Annotation Pipeline has identified 21,517 protein coding genes.

Emergence of a bufonid herpesvirus in a population of the common toad Bufo bufo in Germany

Bufonid herpesvirus 1 (BfHV1) was initially described in 2014 from cases of mortalities and dermatitis in Swiss populations of the common toad Bufo bufo. We identified a closely related herpesvirus strain in a German common toad population affected by an ongoing epidemic of multifocal proliferative to ulcerative skin disease since 2018.

Hybridization and introgression between toads with different sex chromosome systems

The growing interest in the lability of sex determination in non‐model vertebrates such as amphibians and fishes has revealed high rates of sex chromosome turnovers among closely related species of the same clade. Can such lineages hybridize and admix with different sex‐determining systems, or could the changes have precipitated their speciation? We addressed these questions in incipient species of toads (Bufonidae), where we identified a heterogametic transition and characterized their hybrid zone with genome‐wide markers (RADseq). Adult and sibship data confirmed that the common toad B. bufo is female heterogametic (ZW), while its sister species the spined toad B. spinosus is male heterogametic (XY). Analysis of a fine scale transect across their parapatric ranges in southeastern France unveiled a narrow tension zone (∼10 km), with asymmetric mitochondrial and nuclear admixture over hundreds of kilometers southward and northward, respectively. The geographic extent of introgression is consistent with an expansion of B. spinosus across B. bufo’s former ranges in Mediterranean France, as also suggested by species distribution models. However, widespread cyto‐nuclear discordances (B. spinosus backrosses carrying B. bufo mtDNA) run against predictions from the dominance effects of Haldane's rule, perhaps because Y and W heterogametologs are not degenerated. Common and spined toads can thus successfully cross‐breed despite fundamental differences in their sex determination mechanisms, but remain partially separated by reproductive barriers. Whether and how the interactions of their XY and ZW genes contribute to these barriers shall provide novel insights on the debated role of labile sex chromosomes in speciation.

Relationships Between Chemical Defenses of Common Toad (Bufo bufo) Tadpoles and Bacterial Community Structure of their Natural Aquatic Habitat

Many organisms synthesize secondary metabolites against natural enemies. However, to which environmental factors the production of these metabolites is adjusted to is poorly investigated in animals, especially so in vertebrates. Bufadienolides are steroidal compounds that are present in a wide range of plants and animals and, if present in large quantities, can provide protection against natural enemies, such as pathogens. In a correlative study involving 16 natural populations we investigated how variation in bufadienolide content of larval common toads (Bufo bufo) is associated with the bacterial community structure of their aquatic environment. We also evaluated pond size, macrovegetation cover, and the abundance of predators, conspecifics and other larval amphibians. We measured toxin content of tadpoles using HPLC-MS and determined the number of bufadienolide compounds (NBC) and the total quantity of bufadienolides (TBQ). AICc-based model selection revealed strong relationships of NBC and TBQ with bacterial community structure of the aquatic habitat as well as with the presence of conspecific tadpoles. The observed relationships may have arisen due to adaptation to local bacterial communities, phenotypic plasticity, differential biotransformation of toxin compounds by different bacterial communities, or a combination of these processes. Bacterial groups that contribute to among-population variation in toxin content remain to be pinpointed, but our study suggesting that toxin production may be influenced by the bacterial community of the environment represents an important step towards understanding the ecological and evolutionary processes leading to microbiota-mediated variation in skin toxin profiles of aquatic vertebrates.

Testing an hypothesis of hybrid zone movement for toads in France

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.

Effect-based environmental monitoring for thyroid disruption in Swedish amphibian tadpoles

It is well-known that the metamorphosis process in amphibians is dependent on thyroid hormones. Laboratory studies have shown that several environmental contaminants can affect the function of thyroid hormones leading to alterations in the amphibian metamorphosis. The basic idea of the present study was to elucidate if the amphibian metamorphosis might be a useful tool as biomarker for effect-based environmental monitoring, examining wild tadpoles for potential thyroid hormone disruption. A laboratory test was performed to identify the responses from exposure to 6-propylthiouracil (PTU), which has a well-known mechanism on the thyroid system, on Swedish tadpoles from the Rana genus. This was followed by an environmental monitoring study where tadpoles of Rana arvalis, R. temporaria, and Bufo bufo were sampled from various sites in Sweden. Morphological data such as body weight, histopathological measurements of the thyroid glands, and environmental parameters were recorded. The results revealed that Rana tadpoles respond similar as other amphibians to PTU exposure, with interrupted development and increased size relative to the developmental stage. Data on some wild tadpoles showed similar features as the PTU exposed, such as high body weight, thus suggesting potential thyroid disrupting effects. However, histological evaluation of thyroid glands and pesticide analyses of the water revealed no clear evidence of chemical interactions. To a minor degree, the changes in body weight may be explained by natural circumstances such as pH, forest cover, and temperature. The present study cannot fully explain whether the high body weights recorded in some tadpoles have natural or chemical explanations. However, the study reveals that it is clearly achievable to catch tadpoles in suitable stages for the use in this type of biomonitoring and that the use of these biomarkers for assessment of thyroid disruption seems to be highly relevant.

The toad fly Lucilia bufonivora: its evolutionary status and molecular identification

The blow fly genus Lucilia is composed largely of saprophages and facultative myasis agents, including the economically important species Lucilia cuprina (Wiedemann) (Diptera: Calliphoridae) and Lucilia sericata (Meigen). Only one species is generally recognized as an obligate agent of myiasis, Lucilia bufonivora Moniez, and this is an obligate parasite of toads. Lucilia silvarum (Meigen), a sister species, behaves mainly as a carrion breeder; however, it has also been reported as a facultative parasite of amphibians. Morphologically, these species are almost identical, and historically this has led to misidentification, taxonomic ambiguity and a paucity of studies of L. bufonivora. In this study, dipterous larvae were analysed from toad myiasis cases from the U.K., The Netherlands and Switzerland, together with adult specimens of fly species implicated in amphibian parasitism: L. bufonivora, L. silvarum and Lucilia elongata Shannon (from North America). Partial sequences of two genes, cox1 and ef1α, were amplified. Seven additional blow fly species were analysed as outgroups. Bayesian inference trees of cox1, ef1α and a combined-gene dataset were constructed. All larvae isolated from toads were identified as L. bufonivora and no specimens of L. silvarum were implicated in amphibian myiasis. This study confirms L. silvarum and L. bufonivora as distinct sister species and provides unambiguous molecular identification of L. bufonivora.

Indication of spatially random occurrence of Chlamydia-like organisms in Bufo bufo tadpoles from ponds located in the Geneva metropolitan area

Occurrence of bacteria belonging to the order Chlamydiales was investigated for the first time in common toad (Bufo bufo) tadpole populations collected from 41 ponds in the Geneva metropolitan area, Switzerland. A Chlamydiales-specific real-time PCR was used to detect and amplify the Chlamydiales 16S ribosomal RNA–encoding gene from the tails of 375 tadpoles. We found the studied amphibian populations to host Chlamydia-like organisms (CLOs) attributable to the genera Similichlamydia, Neochlamydia, Protochlamydia and Parachlamydia (all belonging to the family Parachlamydiaceae), Simkania (family Simkaniaceae) and Estrella (family Criblamydiaceae); additionally, DNA from the genus Thermoanaerobacter (family Thermoanaerobacteriaceae) was detected. Global autocorrelation analysis did not reveal a spatial structure in the observed CLOs occurrence rates, and association tests involving land cover characteristics did not evidence any clear effect on CLOs occurrence rates in B. bufo. Although preliminary, these results suggest a random and ubiquitous distribution of CLOs in the environment, which would support the biogeographical expectation ‘everything is everywhere’ for the concerned microorganisms.

A matter of proportion? Associational effects in larval anuran communities under fish predation

In Batesian mimicry, a species lacking defences against predators benefits from mimicking the aposematic signal of a defended species, while the model may incur the costs of reduced defensive efficacy. Similar reciprocal indirect effects may emerge even when the signal is not mimicked; termed associational effects, such interactions are well known in plants sharing herbivores but have received little attention in animal studies. We investigated associational interactions in a system where unequally defended prey (chemically defended Bufo bufo and undefended Rana temporaria tadpoles), sharing general morphology but not an aposematic signal, were exposed to predation by the carp Cyprinus carpio along a gradient of relative prey abundance. In the absence of fish, the assemblage composition had no effect on the survival of Rana, while that of Bufo decreased with increasing abundance of Rana. Fish reduced the survival of tadpoles from both species. However, increased relative abundance of Bufo in the community led to enhanced survival in both Bufo and Rana. Increasing relative proportions of heterospecifics reduced metamorph mass only in Bufo, indicating greater sensitivity to interspecific competition compared to Rana; the effect was reduced in the presence of fish. Our results show that undefended non-mimetic prey enjoy reduced predation with increasing relative abundance of chemically defended prey, which in turn suffer greater mortality with an increasing proportion of the undefended species. Associational resistance/susceptibility, driven by current assemblage composition, not by selection for resemblance, can shape the dynamics of mixed communities of defended and undefended prey in the absence of mimicry.

Toxic activity and protein identification from the parotoid gland secretion of the common toad Bufo bufo

Anuran toxins released from the skin glands are involved in defence against predators and microorganisms. Secretion from parotoid macroglands of bufonid toads is a rich source of bioactive compounds with the cytotoxic, cardiotoxic and hemolytic activity. Bufadienolides are considered the most toxic components of the toad poison, whereas the protein properties are largely unknown. In the present work, we analysed the cardio-, myo-, and neurotropic activity of extract and the selected proteins from Bufo bufo parotoids in in vitro physiological bioassays carried out on two standard model organisms: beetles and frogs. Our results demonstrate a strong cardioactivity of B. bufo gland extract. The toad poison stimulates (by 16%) the contractility of the insect heart and displays the cardioinhibitory effect on the frog heartbeat frequency (a 27% decrease), coupled with an irreversible cardiac arrest. The gland extract also exhibits significant myotropic properties (a 10% decrease in the muscle contraction force), whereas its neuroactivity remains low (a 4% decrease in the nerve conduction velocity). Among identified peptides present in the B. bufo parotoid extract are serine proteases, muscle creatine kinase, phospholipid hydroperoxide glutathione peroxidase, cytotoxic T-lymphocyte protein, etc. Some proteins contribute to the cardioinhibitory effect. Certain compounds display the paralytic (myo- and neurotropic) properties. As the toad gland extract exhibits a strong cardiotoxic activity, we conclude that the poison is a potent agent capable of slaying a predator. Our results also provide the guides for the use of toad poison-peptides in therapeutics and new drug development.

Age-dependent changes in sensitivity to a pesticide in tadpoles of the common toad (Bufo bufo)

The worldwide en masse application of pesticides and the frequently reported malign effects on several non-target organisms underpin the importance of ecotoxicological research on these anthropogenic pollutants. Previous studies showed that sensitivity to herbicides can vary widely depending on additional stress factors, on the species and even on the population investigated. However, there is little information about how sensitivity changes during ontogeny, and how the duration of exposure is linked to the magnitude of malign effects, even though this knowledge would be important for the interpretation of toxicity test results and for formulating recommendations regarding the timing of pesticide application. We exposed tadpoles of the common toad (Bufo bufo) to three concentrations (0, 2 and 4mg a.e./L) of a glyphosate-based herbicide during the 1st, 2nd, 3rd, 4th, or 5th period of larval development or during the entire experiment, and measured survival, time until metamorphosis and body mass at metamorphosis to estimate fitness-consequences. Younger tadpoles were more sensitive to the herbicide in all measured traits than older ones, and this age-dependence was especially pronounced at the high herbicide concentration. Furthermore, tadpoles exposed to the herbicide during the entire experiment developed slower than tadpoles exposed only early on, but we did not observe a similar effect either on body mass or survival. The observed age-dependence of sensitivity to herbicides draws attention to the fact that results of toxicity tests obtained for one age-class are not necessarily generalizable across ontogeny. Also, the age of test animals has to be considered when planning ecotoxicological studies and interpreting their results. Finally, taking into account the temporal breeding habits of local amphibians when planning pesticide application would be highly favourable: if tadpoles would not get exposed to the herbicide during their most sensitive early development, they would sustain less anthropogenic damage from our efforts of controlling weeds.

Screening breeding sites of the common toad (Bufo bufo) in England and Wales for evidence of endocrine disrupting activity

Anuran amphibians are often present in agricultural landscapes and may therefore be exposed to chemicals in surface waters used for breeding. We used passive accumulation devices (SPMD and POCIS) to sample contaminants from nine breeding sites of the Common toad (Bufo bufo) across England and Wales, measuring endocrine activity of the extracts in a recombinant yeast androgen screen (YAS) and yeast estrogen screen (YES) and an in vitro vitellogenin induction screen in primary culture of Xenopus laevis hepatocytes. We also assessed hatching, growth, survival, and development in caged larvae in situ, and sampled metamorphs for gonadal histopathology. None of the SPMD extracts exhibited estrogen receptor or androgen receptor agonist activity, while POCIS extracts from two sites in west-central England exhibited concentration-dependent androgenic activity in the YAS. Three sites exhibited significant estrogenic activity in both the YES and the Xenopus hepatocyte. Hatching rates varied widely among sites, but there was no consistent correlation between hatching rate and intensity of agricultural activity, predicted concentrations of agrochemicals, or endocrine activity measured in YES/YAS assays. While a small number of intersex individuals were observed, their incidence could not be associated with predicted pesticide exposure or endocrine activitity measured in the in vitro screens. There were no significant differences in sex ratio, as determined by gonadal histomorphology among the study sites, and no significant correlation was observed between proportion of males and predicted exposure to agrochemicals. However, a negative correlation did become apparent in later sampling periods between proportion of males and estrogenic activity of the POCIS sample, as measured in the YES. Our results suggest that larval and adult amphibians may be exposed to endocrine disrupting chemicals in breeding ponds, albeit at low concentrations, and that chemical contaminants other than plant protection products may contribute to endocrine activity of surface waters in the agricultural landscape.

How complex is the Bufo bufo species group?

Species delineation remains one of the most challenging tasks in the study of biodiversity, mostly owing to the application of different species concepts, which results in contrasting taxonomic arrangements. This has important practical consequences, since species are basic units in fields like ecology and conservation biology. We here review molecular genetic evidence relevant to the systematics of toads in the Bufo bufo species group (Anura, Bufonidae). Two studies recently published in this journal (Recuero et al., MPE 62: 71-86 and García-Porta et al., MPE 63: 113-130) addressed this issue but reached opposing conclusions on the taxonomy of the group (four versus two species). In particular, allozyme data in the latter paper were interpreted as evidence for hybridization across species (between B. bufo-B. spinosus and B. bufo-B. verrucosissimus). We tested claims for hybridization through re-analysis of allozyme data for individuals instead of populations, to be able to distinguish between sympatry with and without admixture, and found no evidence of hybridization across taxa. We propose alternative explanations for the observed patterns that García-Porta et al. (2012) failed to consider. In the absence of unequivocal evidence for hybridization and introgression, we reject the proposal to downgrade Bufo spinosus and Bufo verrucosissimus to the subspecies level.

Types of supernumerary outgrowths produced after inverting the dorsoventral limb axis of the anuranBufo bufo

Contralateral grafts were performed on the larval limb buds of the anuranBufo bufo. The dorsoventral axis of 80 buds at stages IV or V was inverted. Ten tadpoles were used as controls. Fifty-two supernumerary structures developed, all of them in dorsal or ventral locations on the host stump. The majority (32 out of the 44 outgrowths with more than 3 toes) were normal limbs of stump handedness. However, the following abnormal structures were also observed: 2 double-posterior, 3 mixed-symmetric, and 7 undetermined cases. These results are in agreement with the predictions of a hierarchical polar coordinate model for epimorphic regeneration.

Seasonal changes in the migratory behaviour of the toad Bufo bufo: direction and magnitude of movements

The migratory behaviour of the toad Bufo bufo was studied from February 1985 to April 1986 in the submontane region of Bavaria, West Germany. Toads were fitted with a mechanical tracking device to record individual paths of migration. Three aspects of migratory behaviour were quantified: orientation in relation to the breeding site, straightness of path, and locomotory activity. The annual activity period began with migration from the hibernation sites to the breeding pond in April. The paths went straight towards the breeding pond independent of the distance (70-420 m). During the period of oviposition the preference for the breeding site direction vanished and toads moved away from the breeding pond, but in less straight paths than before. In summer migratory activity decreased considerably and was restricted to small areas, the home ranges, at distances of 55-1600 m from the natal breeding pond. The straightness of path was rather low, because toads often returned to their starting points. During rainy nights toads occasionally left their home ranges for extensive excursions. In autumn most toads again migrated towards the breeding pond, but paths were significantly less straight and direct than in spring. However, toads stopped before reaching the breeding pond and hibernated in holes or under the leaf layer. The mortality rate of tracked toads was about 45%. The relative influence of 17 environmental variables on locomotory activity was evaluated by principal component analysis and stepwise multiple regression. Temperature at night and rainfall variables accounted for significant amounts of variance, whereas temperature by day, air humidity, and atmospheric pressure showed no correlation. Activity decreased if temperature approached 0° C or after long periods without rainfall. Within a certain range of tolerance, however, the locomotory activity of the toads was widely independent of environmental factors, indicating that endogeneous factors are more important sources of variation in the migratory behaviour of these toads than commonly assumed.