Skip to main content

Biological characterization of compounds from Rhinella schneideri poison that act on the complement system

Abstract

Background

The skin secretions of toads of the family Bufonidae contain biogenic amines, alkaloids, steroids (bufotoxins), bufodienolides (bufogenin), peptides and proteins. The poison of Rhinella schneideri, formerly classified as Bufo paracnemis, presents components that act on different biological systems, including the complement system. The aim of this study was to isolate and examine the activity of Rhinella schneideri poison (RsP) components on the complement system.

Methods

The components active on the complement system were purified in three chromatographic steps, using a combination of cation-exchange, anion-exchange and gel filtration chromatography. The resulting fractions were analyzed by SDS-PAGE and screened for their activity in the hemolytic assay of the classical/lectin complement pathways. Fractions active on the complement system were also assessed for their ability to generate C3 fragments evaluated by two dimensional immunoelectrophoresis assay, C3a and C5a by neutrophil chemotaxis assay and SC5b-9 complex by ELISA assay.

Results

The fractionation protocol was able to isolate the component S5 from the RsP, as demonstrated by SDS-PAGE and the RP-FPLC profile. S5 is a protein of about 6000 Da, while S2 presents components of higher molecular mass (40,000 to 50,000 Da). Fractions S2 and S5 attenuated the hemolytic activity of the classical/lectin pathways after preincubation with normal human serum. Both components stimulated complement-dependent neutrophil chemotaxis and the production of C3 fragments, as shown by two-dimensional immunoelectrophoresis. S2 showed a higher capacity to generate the SC5b-9 complex than the other fractions. This action was observed after the exposure of normal human serum to the fractions.

Conclusions

This is the first study to examine the activity of RsP components on the complement system. Fractions S2 and S5 reduced the complement hemolytic activity, stimulated complement-dependent neutrophil chemotaxis and stimulated the production of C3 fragments, indicating that they were able to activate the complement cascade. Furthermore, fraction S2 was also able to generate the SC5b-9 complex. These components may be useful tools for studying dysfunction of the complement cascade.

Background

The family Bufonidae, with more than 590 species distributed among 50 genera, is one of the largest Anuran families [1]. The genus Rhinella is composed of 88 species, of which 36 are found in Brazil [1]. Rhinella schneideri, previously known as Bufo paracnemis, is the species most commonly encountered in Brazil [2, 3].

Amphibian skin secretions contain a large number of biologically active compounds that are involved in the regulation of physiological functions of the skin, as well as in defense mechanisms against predators and microorganisms [4]. The skin glands produce mucus, peptides, biogenic amines, steroids, and alkaloids. Pharmacologically, these substances may be neurotoxic, cardiotoxic, hemotoxic or myotoxic, and can provoke anesthetic, hypotensive and/or hypertensive effects [5, 6].

Dried poison from the skin glands of the Chinese toad (Bufo bufo gargarizans cantor) has been used as a therapeutic agent in traditional Chinese medicine, as well as in other Asian countries [79]. Isolated components from toad glands have been used to treat several types of cancer [1015]. A previous report described the influence of Rhinella schneideri poison (RsP) on the lytic activity of the complement system [16].

The complement system (CS) is one of the main defense mechanisms of vertebrates and encompasses over 30 proteins, some of which circulate in the plasma as precursors. Depending on the stimulus, complement activation occurs by classical, alternative or lectin pathways (CP, AP and LP, respectively), leading to a cascade of component interactions and the generation of products that can exert biological activities such as anaphylaxis, chemotaxis, opsonization, immune complex solubilization and participation in the immune response. After recognition, a series of serine proteases is activated, culminating in formation of the “membrane attack complex” (MAC) within the membrane that leads to lysis or cell activation. Two important mediators of the inflammatory reaction, C3a and C5a, are produced as a consequence of CS activation. However, inappropriate activation can result in substantial injury. To prevent undesired complement activation, inhibitors acting at different stages of the activation pathways are used. Despite the large number of inhibitory compounds identified so far, there is still a need for selective complement system modulators [1719].

Since the poisonous secretion of the parotoid gland of the R. schneideri toad presents anticomplement activity, this work aimed to purify the active components and to investigate their effects on the complement system [16].

Methods

Poison

The poison was collected by applying pressure to parotoid glands of Rhinella schneideri toads, immediately desiccated under vacuum and stored at –20 °C until usage. Prior to the assays, the poison or toxin solutions were filtered through sterilizing membranes (Merck-Millipore, Germany – cellulose ester filters: 0.45 μm and 0.22 μm, respectively).

Experimental animals

An adult male sheep from the animal facility of the University of São Paulo in Ribeirão Preto was kept in accordance with the ethical guidelines established by the Brazilian College of Animal Experimentation (COBEA). All experiments were approved and conducted in accordance with the ethical principles in animal experimentation adopted by the Ethics Commission for the Use of Animals (CEUA), Campus of Ribeirão Preto, USP (protocol no 05.1.637.53.6).

Fractionation of R. schneideri poison

The soluble material from the desiccated poison (500 mg) was clarified by filtration through membranes (0.45 μm and then 0.22 μm, Merck-Millipore, Germany). The material was chromatographed at 4 °C on a 2.5 × 63.0 cm column of CM-cellulose-52 (Whatman, USA), which was equilibrated and initially eluted with 300 mL of 0.05 M NH4HCO3 buffer, pH 7.8, when a convex concentration gradient was started from 0.05 to 1.00 M NH4HCO3 buffer. Fractions of 3.0 mL were collected. Absorbance at 280 nm and buffer concentration profiles were then traced as previously described [20].

The resulting pools, designated C1 to C7, were then lyophilized until salt-free. The fraction C1 showed the lowest percentage of hemolysis. Therefore, C1 was submitted to the next fractionation step. The soluble material from fraction C1 (56.6 mg in 5 mL of 0.05 M Tris–HCl, pH 7.8, centrifuged at 15,700 × g, at 4 °C, for 10 min) was applied on a 1.0 × 10.0 cm DEAE-Sepharose column at room temperature, previously equilibrated with 0.05 M Tris–HCl, pH 7.8 (buffer A). Elution was performed with a linear gradient of buffer B (0.05 M Tris–HCl supplemented with 1.0 M NaCl, pH 7.8), at a flow rate of 0.5 mL/min. Absorbance was monitored at 280 nm. The chromatography was performed in an Äkta™ Prime system (GE Healthcare, Sweden) and the resulting pools, denominated D1 to D4, were lyophilized. Fraction D3 showed the highest activity on the complement system and was submitted to a molecular filtration on a Sephacryl S-200 column (1.6 cm × 60 cm) at room temperature, previously equilibrated with PBS (phosphate buffered saline), pH 7.4, at a flow rate of 0.4 mL/min. The absorbance was monitored at 254 nm. The resulting pools, designated S1 to S5 were grouped according to their respective absorbance peaks and stored at –20 °C.

Fractions S2 and S5, which showed activity on the complement system, were submitted to a reversed phase FPLC using a C2C18 column (0.46 × 10 cm, Amersham Biosciences, Sweden). The column was equilibrated with 0.1 % (V/V) trifluoroacetic acid (TFA, solution A); and the components were eluted by a step concentration gradient from 0 to 100 % of solution B (80 % acetonitrile, 0.1 % trifluoroacetic acid, V/V), at a flow rate of 0.5 mL/min, at room temperature. The absorbance (λ = 214 nm) was registered by the Äkta™ Prime system (GE Healthcare, Sweden).

Polyacrylamide gel electrophoresis

Sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) was run as described by Laemmli [21]. The gel was stained with Silver Staining Kit Protein (Pharmacia Biotech, Sweden) or Coomassie Blue R-350. Conditions of voltage and amperage (maximum values: 90 V, 40 mA and 15 W) were controlled by an EPS 3500 XL Electrophoresis Power Supply (Pharmacia Biotech, Sweden).

Solutions

Cells were washed in PBS, pH 7.4, and complement fixation diluent (CFD) containing 0.1 % gelatin (gel) was used for hemolytic assays of CP/LP activity as described by Harrison and Lachmann [22]. Modified Alsever’s solution [23] was used as an anticoagulant for sheep blood storage.

Normal human serum (NHS) and erythrocytes

Human blood was obtained from healthy donors (approval certificate by the Research Ethics Committee – CAAE, protocol n° 0022.0.212.000–08). Blood samples were collected from healthy volunteers of both sexes (aged 20 to 30 years) without anticoagulant and allowed to clot for one hour at room temperature, after which they were centrifuged at 556 × g, for ten minutes at 4 °C, and the NHS obtained was stored at –70 °C.

Healthy adult male sheep were bled by jugular vein puncture; the blood was collected in two volumes of Alsever’s modified solution, stored at 4 °C and utilized for 15 days as a source of erythrocytes for CP/LP hemolytic assays. Sheep blood was centrifuged (556 × g, 15 min, 4 °C), after which the plasma and buffy coat were discarded. Red cells were washed twice in PBS, suspended in CFD/Gel and mixed with an appropriate volume of anti-sheep erythrocyte antibody. This erythrocyte-antibody suspension was maintained at 4 °C for 15 min and its absorbance at 700 nm was adjusted to 0.70–0.80.

Hemolytic complement assay

NHS was diluted in CFD/Gel at a ratio of 1:20, V/V. Fractions (100 μL in PBS, S1 – A280 ~ 0.35; S2 – A280 ~ 0.20; S3 – A280 ~ 0.17; S4 – A280 ~ 0.10 and S5 – A280 ~ 0.16) obtained from D3 molecular filtration (Sephacryl S-200) were incubated with CFD/Gel solution (12.5 μL) and diluted serum (1:20; 37.5 μL) for one hour at 37 °C. After the incubation period, erythrocyte-antibody suspension (100 μL) was added to the samples and a new incubation was performed for 30 min at 37 °C. At the end of the incubation, cold PBS (250 μL) was added to the samples, which were then centrifuged at 556 × g for ten minutes. The percentage of lysis was determined by the absorbance at 412 nm, using as 100 % lysis control the suspension of lysed erythrocytes in water, and as 0 % lysis control the cells incubated in CFD/Gel. The positive control was prepared under the same reaction conditions except that the fraction volume was replaced by PBS (100 μL). This assay was employed to monitor the activity of the fractions on the complement system during the purification process.

Human neutrophils suspension

Human blood from healthy donors was mixed with modified Alsever’s solution (V/V) and centrifuged at 978 × g for ten minutes. Neutrophils were isolated by the gelatin method, as described by Paula et al. [24] with modifications. Briefly, after blood centrifugation, the plasma and buffy coat were discarded, and the cell pellet was suspended in two volumes of 2.5 % gelatin solution prepared in 0.15 M NaCl. This suspension was incubated for 15 min at 37 °C. After incubation, the upper neutrophil-rich layer was collected, diluted in 30 mL of 0.15 M NaCl solution and centrifuged at 757 × g for ten minutes at room temperature. The cell pellet was suspended in 20 mL of 0.83 % NH4Cl solution, pH 7.8, and incubated for five minutes at 37 °C, in order to lyse remaining erythrocytes. After incubation, the supernatant was discarded and the suspension was centrifuged at 757 × g for ten minutes at room temperature. The cell pellet was washed in 30 mL of 0.15 NaCl solution and centrifuged at 757 × g for ten minutes. The supernatant was discarded and 1 mL of the neutrophil suspension was suspended in 1 mL of Hank’s solution containing 0.1 % gelatin. Cells were diluted (1:10) in Turk solution and counted in a Neubauer Chamber. The neutrophil purity of 80–90 % with viability higher than 95 % was accomplished by trypan blue exclusion test. One neutrophil suspension was standardized to contain 1.2 x 106 cells mL−1 and used in the neutrophil chemotaxis assay.

Neutrophil chemotaxis assay

The chemotaxis assay was performed using a modified version of Boyden’s technique [25], in which 120 μL of NHS with 50 μL of CFD/Gel and 50 μL of each fraction, S1 (A280 ~ 0.35), S2 (A280 ~ 0.20), S3 (A280 ~ 0.17), S4 (A280 ~ 0.10) and S5 (A280 ~ 0.16), obtained from molecular filtration of D3, were placed in the lower migration chamber and covered with a filter of 13 mm diameter and 3 μm pore (SSWP 01300, Merck-Millipore, Germany). The upper compartment of the chamber was filled with 300 μL of a suspension of human neutrophils (1.2 × 106 cells mL−1). Next, all chambers were closed and incubated at 37 °C for 60 min in a humid atmosphere. After incubation, the filters were removed from the chambers, fixed in propanol, stained with Harris hematoxylin, dehydrated in isopropanol and cleared with xylene. Each filter was placed between a slide and a coverslip with Entellan (Merck KGaA, Germany). A mixture of NHS (120 μL) with CFD/Gel (100 μL) and zymosan (75 μL, 1 mg/mL) was used as positive control, and NHS (120 μL) with CFD/Gel (100 μL) as negative control.

Neutrophil migration within the filter was determined under light microscopy by the leading-front method, measuring in micrometers the greatest distance crossed by three cells per field [26]. At least ten fields were examined at 100× magnification for each filter.

Two-dimensional immunoelectrophoresis (2D-IEP)

For this analysis, 50 μL of fractions S2 (A280 ~ 0.2) and S5 (A280 ~ 0.16) were preincubated in a water bath with 100 μL of NHS 1:2 by 60 min at 37 °C. Immunoelectrophoresis was performed according to the method of Clark and Freeman [27], using glass plates (5.5 × 7.5 × 0.2 cm) and 1.3 % agarose in buffer (0.025 M Tris–HCl, 0.027 M glycine, 0.02 M sodium barbital, 0.01 M EDTA, pH 8.8). In the first dimension, the positive control (31.25 μL of zymosan plus 100 μL of 1:2 NHS), the negative control (100 μL of 1:2 NHS plus 50 μL of PBS) and fraction S2 and S5 (50 μL of each fraction plus 100 μL 1:2 NHS) were electrophoresed for four hours, at 140 V and 5 mA/plate. For the second dimension, the plates were completed with 1.3 % agarose (5 mL) containing 1 % anti-human C3 antibody (Calbiochem/Merck, Germany) and electrophoresed for 14 h, at 10 W and 5 mA/plate. The plates were dried at room temperature, stained with Ponceau 0.5 % and destained with 10 % acetic acid.

Evaluation of the capacity to generate the SC5b-9 complex

The capacity of fractions (S1 to S5) to generate the SC5b-9 complex was evaluated by enzyme-linked immunosorbent assay (ELISA, Quidel SC5b-9 Complement® kit, USA) after exposure of the NHS to 50 μL of each fraction [28].

Statistical analysis

The results were expressed as the mean ± SEM. The groups were compared statistically by ANOVA followed by the Tukey post-hoc test. All data were analyzed via Prism™ v.5 (GraphPad Inc., USA).

Results

Fractionation of R. schneideri poison

The components from Rhinella schneideri poison with activity on the CS were obtained by three chromatographic steps: cation-exchange, anion-exchange and molecular exclusion. The chromatographic profile of soluble poison on CM-cellulose-52 (cation-exchange) showed seven different fractions, denominated C1 to C7 (Fig. 1a). Fraction C1 presents the highest inhibition of hemolytic complement activity, as previously demonstrated by our group [29]. RsP and fraction C1 were assayed by SDS-PAGE (Fig. 1b), where C1 appeared as a complex mixture of proteins. Therefore, it was submitted to the next fractionation step on a DEAE-Sepharose column (Fig. 1c).

Fig. 1
figure 1

Fractionation of Rhinella schneideri poison (RsP). a Chromatographic profile of RsP on CM-cellulose-52. The column was equilibrated with 0.05 M ammonium bicarbonate, pH 7.8. The sample (extract from 500 mg) was applied at a flow rate of 20 drops/min; and adsorbed /components were eluted using a convex concentration gradient of NH4HCO3 (0.05 to 1.0 M, pH 7.8). Fractions (3.0 mL/tube) were collected at 4 °C. b SDS-PAGE using 13.5 % separation gel. Lanes 1, 2 and 3 – ultra- mass markers, respectively. Lane 4 – fraction C1; lanes 5 and 6 – RsP. c Chromatographic low (Sigma-Aldrich, USA), low-high (GE Healthcare, Sweden) and high (GE Healthcare, Sweden)molecular profile of fraction C1 on DEAE-Sepharose. The column was equilibrated with 0.05 M Tris–HCl, pH 7.8 (buffer A). The sample (56.6 mg of C1) was applied at a flow rate of 0.5 mL/min; and adsorbed components were eluted using a linear gradient from 0–1 M NaCl in equilibrating buffer (buffer B). Elution with 100 % buffer B was achieved after 150 mL. d SDS-PAGE using 13.5 % separation gel. Lane 1 – RsP; lane 2 – fraction C1; lane 3 – fraction D3; lane 4 – ultra-low molecular mass markers (Sigma-Aldrich, USA). e Chromatographic profile of fraction D3 on Sephacryl S-200. The column, equilibrated with PBS, pH 7.4, was eluted with this same buffer (flow rate: 0.4 mL/min) and fractions of 1 mL were collected. In (a) and (c), the elution profiles were monitored at 280 nm, while in (e) the profile was monitored at 254 nm. f SDS-PAGE using 13.5 % separation gel. Lanes 1 and 2 – fraction S5; lanes 3 and 5 – ultra-low-molecular-mass markers (Sigma-Aldrich, USA); lanes 4 and 6 – low-molecular-mass markers (GE Healthcare, Sweden); lane 7 – fraction S2

Among the four fractions (D1, D2, D3 and D4) obtained from the rechromatography of fraction C1, the fraction D3 presented the highest activity on CS [29]. Unfortunately, it was composed of low- and high-molecular-mass components, highlighting a protein with an approximate molecular mass of 6 kDa, observed in the SDS-PAGE (Fig. 1d). In order to isolate some component presenting action on complement system, fraction D3 was submitted to a gel filtration on a Sephacryl-S200 column (Fig. 1e). Five fractions, designated S1 to S5, were obtained; and the active fractions S2 and S5 were analyzed by SDS-PAGE (Fig. 1f). Fraction S2 was composed of high-molecular-mass proteins (40,000 to 50,000 Da), while S5 was a protein of about 6 kDa. The recoveries of chromatographic fractions with activity on CS are shown in Table 1.

Table 1 Recovery of the chromatographic components obtained during the fractionation procedure

The active fractions S2 and S5 were submitted to a reversed phase FPLC using a C2C18 column (Fig. 2). S5 showed higher purity than S2, which presented a chromatographic profile with two major peaks, S2.1 and S2.2.

Fig. 2
figure 2

Reversed-phase FPLC of fractions S2 and S5. The C2C18 column was equilibrated with 0.1 % (V/V) trifluoroacetic acid (TFA, solution A). Adsorbed proteins were eluted using a concentration gradient from 0 to 100 % of solution B (80 % acetonitrile in 0.1 % TFA, V/V). Fractions of 0.5 mL/tube were collected at a flow rate of 0.5 mL/min

Hemolytic complement assay

All fractions obtained from the last chromatographic procedure were subjected to the hemolytic assay of the classical/lectin pathway, in which 100 μL volumes of each fraction – S1 (A280 ~ 0.35), S2 (A280 ~ 0.2), S3 (A280 ~ 0.17), S4 (Ab280 ~ 0.1) and S5 (A280 ~ 0.16) – were used. The hemolytic activities of the classical/lectin pathways observed in the presence of all fractions were significantly lower than the positive control, particularly in the presence of fractions S2 and S5 (Fig. 3).

Fig. 3
figure 3

Effect of fractions S1 - S5 on the classical/lectin hemolytic pathways of complement activation. The fractions (100 μL in PBS, S1 – A280 ~ 0.35; S2 – A280 ~ 0.20; S3 – A280 ~ 0.17; S4 – A280 ~ 0.10 and S5 –A280 ~ 0.16)) were incubated for one hour at 37 °C with normal human serum diluted 1:20 (37.5 μL) and CFD/Gel solution (12.5 μL). The positive control was run under the same conditions but in the absence of fractions. The absorbances of supernatants from cells incubated in CFD/Gel buffer (0 % lysis) and lysed in water (100 % lysis) were employed to calculate the percentage of lysis induced by NHS in the absence (positive control) or presence of the fractions (tests). The columns represent means ± SEM of an experiment performed in duplicate. ** p < 0.01, *** p < 0.001 and **** p < 0.0001 compared to the positive control

Two-dimensional immunoelectrophoresis

The 2D-IEP profile of the positive control shows two protein peaks, corresponding to C3 and C3b, indicating the ability of zymosan to activate the complement system, leading to partial cleavage of C3 (Fig. 4a). A symmetrical peak was observed in the negative control, corresponding to intact C3. The 2D-IEP profiles obtained in the presence of S2 and S5 (Fig. 4b and c, respectively) also showed two peaks, similar to the positive control, corresponding to C3 and C3b, indicating that S2 and S5 were able to activate the complement system. The background of these 2D-IEP profiles (Fig. 4a, b and c) were removed in order to highlight the presence of one peak for negative control and two peaks for positive control, S2 and S5 assays (Fig. 4d).

Fig. 4
figure 4

Immunoelectrophoretic analysis of C3 in NHS incubated with fractions S2 and S5. a Positive control (C+) with zymosan (31.25 μL, 1 mg/mL) and NHS (100 μL, 1:2), and negative control (C-) with PBS (50 μL) and NHS (100 μL, 1:2). b Fraction S2 (50 μL, A280 ~ 0.2) with SHN (100 μL, 1:2). c Fraction S5 (50 μL, A280 ~ 0.16) with NHS (100 μL, 1:2). All mixtures were incubated for 60 min in a water bath at 37 °C. The plates were dried at room temperature, stained with 0.5 % Ponceau and bleached with 10 % acetic acid. Electrophoretic conditions: first dimension – four hours, 140 V at 15 mA and 10 W; second dimension – 14 h at 15 mA and 10 W. d Figures were manipulated to remove the background highlighting the presence of one peak for C- and two peaks for C+, S2 and S5 assays

Neutrophil chemotaxis assay

A significant (p < 0.001) increase in neutrophil migration was observed by preincubation of poison components S2 and S5 with NHS (Fig. 5). These results indicate that S2 and S5 were able to induce activation of the complement system, leading to the formation of chemotactic factors.

Fig. 5
figure 5

Neutrophil chemotaxis induced by normal human serum (NHS) preincubated with fractions S1 - S5. The fractions were preincubated with NHS for 60 min at 37 °C. The positive control consisted of 120 μL of NHS with 100 μL of CFD/Gel buffer and 75 μL of zymosan (1 mg/mL), while the negative control was 120 μL of NHS with 100 μL of CFD/Gel buffer. Neutrophil migration was assessed by the leading-front technique, in which at least ten microscopic fields were analyzed per filter at 100× magnification. The columns represent the means ± SEM of one experiment done in triplicate. *** p < 0.001 compared to the negative control

Evaluation of the capacity to generate the SC5b-9 complex

The concentrations of SC5b-9 complex produced after exposure to NHS with fractions S1, S2, S3, S4 and S5, as well as zymosan (positive control) were determined by enzyme-linked immunosorbent assay (Fig. 6). S2 showed significant capacity to generate the SC5b-9 complex compared to the negative control (p < 0.01) and § compared to S4 (p < 0.05).

Fig. 6
figure 6

Formation of the SC5b-9 complex. NHS was incubated for 60 min with PBS (negative control), zymosan (positive control; 1 mg/mL) and fractions S1 to S5. The assay was done using a commercial kit (Quidel SC5b-9 Complement® kit, USA). The columns represent the means ± SEM of one experiment done in duplicate. ** p < 0.01 and **** p < 0.0001 compared to the negative control, and § p < 0.05 compared to S4

Discussion

Research involving animal-derived substances that act on the complement system has been well documented in the literature. Spider (Loxosceles), snake (Elapidae, Crotalidae and Viperidae), honeybee, wasp and scorpion venoms have demonstrated the ability to activate the CS [3035]. This activation can be initiated by cleavage of a specific component or by interaction with other CS components resulting in the formation of the “membrane attack complex” [32]. Tityus serrulatus venom activates the CS, leading to factor B and C3 cleavage, reduction of serum lytic activity and generation of complement chemotactic factors [30]. Assis et al. [16] showed that fractionation of the poisonous secretion of B. marinus paracnemis Lutz (currently named Rhinella schneideri), by dialysis and chromatography on QAE-Sephadex, yielded a fraction with anticomplementary activity when incubated with human serum. This effect was evaluated by measuring the kinetics of lytic activity on sensitized sheep red blood cells (classical pathway) and unsensitized rabbit cells (alternative pathway). A study of the skin secretion of six species of common toads in China revealed that only Bombina maximum poison showed direct hemolytic activity at a dose of 20 μg/mL [36].

This study describes the effects of two components from RsP that interfere with CP/LP of the complement system. To the best of our knowledge, only one study reported the interaction of RsP with the CS until now [16]. The ability of this poison to induce serum-related leukocyte recruitment was evaluated as an indicator of complement activation and consequent generation of complement chemotactic factors.

The procedure used in this study to fractionate the active compound from RsP was relatively simple, involving only three chromatographic steps, a cationic and an anionic chromatography followed by a gel filtration. The SDS-PAGE results show that fraction C1, active on CS, is composed of high- and low-molecular-mass proteins. The major protein of fraction D3 has a molecular mass of about 60,000 Da and corresponds to the isolated S5 active component. On the other hand, the fraction S2 is composed mainly of high-molecular-weight proteins (40,000–50,000 Da). The chromatographic profile of the RP-FPLC S5 confirms the high purity of this component.

The complement activation occurs along classical, alternative or lectin pathways leading to a cascade of component interactions and generation of products that present such biological activities as anaphylaxis, chemotaxis, opsonization, immune complex solubilization, participation in the immune response and other activities [1719]. Two important mediators of the inflammatory reaction, C3a and C5a, are produced as a consequence of CS activation [1719].

The hemolytic complement assay was used to ensure the functional integrity of the whole pathways (classical or alternative) with the terminal pathway. The results obtained showed that all fractions (S1 – A280 ~ 0.35, S2 – A280 ~ 0.2, S3 – A280 ~ 0.17, S4 – Ab280 ~ 0.1 and S5 – A280 ~ 0.16) induced significant reductions in hemolytic activity of classical/lectin pathways, but smaller hemolysis values were obtained in the presence of S2 and S5 (p < 0.0001). The solutions of the fractions (S1-S5) used in the hemolytic assay showed different absorptions at 280 nm. Contrary to what might be expected in this context, our objective was to perform only a qualitative analysis of the effect of the fractions on the CS. The fraction solutions with the highest possible concentration were used, considering their solubility as well as the amount obtained from each fraction (S3 and S4 are present in low proportion in the fraction D3 – Fig. 1e).

We have chosen to use the absorbance at 280 nm and the volume as a means to quantify the samples, because toad poison is composed of proteic (protein and peptides) and non-proteic (mucus, biogenic amines, steroids and alkaloids) compounds. The non-proteic compounds interfere with many protein quantification assays, invalidating the measurement of the sample. According to Marongio [29], the protein concentration, determined by the biuret method, of a dispersion of 5 mg/mL of the R. schneideri poison was only 1.32 mg/mL, corresponding to 26 % of the total weight of the poison.

The reduction in CP/LP lytic activity induced by S2 and S5 suggests an activation of the complement cascade during the preincubation phase (NHS + fractions) and subsequent inactivation (unstable components). The CS activation preceding the addition of red blood cells would reduce serum concentrations of complement components, thereby leading to decreased NHS lytic activity during the hemolytic assay. Similar results were observed in studies of snake venoms from the genera Bothrops (B. jararaca, B. moojeni and B. cotiara), Micrurus (M. ibiboboca and M. spixii) and Naja (N. naja, N. melanoleuca and N. nigricollis) [32, 37, 38].

The presence of an inhibitor of CS in RsP is also possible, since protease inhibitors have been identified on the skin of Anura species [3941]. Several compounds can modify or interact with CS by activating or inhibiting it [16, 3038, 42]. Peptides synthesized from phage-displayed peptide libraries based on C1q binding are able to inhibit the hemolytic activity of the classical complement pathway [43]. Another peptide from phage-displayed libraries, the compstatin peptide, a 13-amino-acid cyclic peptide, binds to a β-chain of C3 and inhibits activation of both the classical and alternative pathways [44, 45].

The immunoelectrophoresis assay showed that cleavage of C3 in serum incubated with S2 and S5 (Fig. 4b and c, respectively), was similar to that induced by incubation of NHS with zymosan (positive control, Fig. 4a), corroborating the hypothesis that poison components induce activation of the CS. Bertazzi et al. [30] showed that Tityus serrulatus venom was also able to alter C3 immunoelectrophoresis migration after incubation with NHS.

The chemotaxis assay served as an indicator of the activation of CS and consequent generation of neutrophil chemoattractant factors. S2 and S5 increased neutrophil migration by interacting with CS components, leading to subsequent cleavage of C3 and C5, which produced the active fragments C3a and C5a (anaphylatoxins) during the preincubation phase (60 min at 37 °C) of NHS with fractions. These results were similar to that presented by zymosan (positive control) and confirm that S2 and S5 are able to activate the complement system. A similar effect was observed in a prior study on Tityus serrulatus venom [30]. BaP1, a 24 kDa metalloprotease isolated from Bothrops asper venom, induced neutrophil chemotaxis that was mediated by agents derived from activation of the complement system [37, 46].

The assay performed to evaluate the capacity of the RsP components to induce SC5b-9 complex formation showed that only S2 was able to induce a significant increase (p <0.01) in the SC5b-9 concentration, compared to negative control (Fig. 6). This assay was conducted to better clarify the action of fractions on the complement system and provided an additional indication of the terminal complement system activation induced by S2.

Primary among the effects of SC5b-9 complex is tissue injury through cell lysis or stimulation of pro-inflammatory mediators [47]. It is known that more than 80 % of C5a and SC5b-9 is generated by activation of the mannose-binding lectin or classical pathway [48, 49]. High levels of activation and generation of SC5b-9 complex are related to several pathological states, including lupus erythematosus and rheumatoid arthritis [47].

Proteolytic activity evaluation of the fractions S2 and S5 was performed using chromogenic substrate for alpha-chymotrypsin (Sigma-S7388, N-Succinyl-Ala-Ala-Pro-Phe p-nitroanilide, Sigma-Aldrich, USA) and for coagulation proteases [Sigma-T6140, N-(p-Tosyl)-Gly-Pro-Lys 4-nitroanilide acetate salt, Sigma-Aldrich, USA]. Additionally, these fractions were subjected to assays to evaluate inhibitory activity against trypsin and chymotrypsin proteases. Neither fraction showed proteolytic or inhibitory activity (data not shown), indicating that their actions on the CS are not by proteolysis or inhibition of the complement cascade proteases.

Several approaches are being proposed for the development of new pharmacological agents directed at diseases in which the CS is active [47, 5053]. Cobra venom factor (CVF) is a non-toxic venom compound with functional and structural characteristics very similar to the complement component C3 [53, 54]. The development of the humanized version of CVF is a promising therapeutic agent for many pathologies [50, 52, 53].

Conclusion

In summary, our results indicate that RsP presents components, especially S2 and S5, that are able to activate the complement cascade, as evidenced by the decreased serum lytic activity, production of C3 fragments, increased leukocyte migration and SC5b-9 generation. Based on these findings, the RsP may be considered a rich source of substances that can be used as molecular tools to study dysfunction of the CS, since they are able to modulate the activity of this system.

Ethics committee approval

All experiments were approved and conducted in accordance with the ethical principles in animal experimentation adopted by the Ethics Commission for the Use of Animals (CEUA), Campus of Ribeirão Preto, USP (protocol no 05.1.637.53.6). The use of human blood was approved by the Research Ethics Committee of the School of Pharmaceutical Sciences of Ribeirão Preto, University of São Paulo (USP) under protocol n° 0022.0.212.000–08.

Abbreviations

AP:

Alternative pathway

CFD:

Complement fixation diluents

CP/LP:

Classical/lectin pathway

CVF:

Cobra venom factor

FPLC:

Fast protein liquid chromatography

Gel:

Gelatin

CS:

Complement system

2D-IEP:

Two-dimensional immunoelectrophoresis

MAC or SC5b-9:

Membrane attack complex

NHS:

Normal human serum

PBS:

Phosphate-buffered saline

RsP:

Rhinella schneideri poison

SDS-PAGE:

Sodium dodecyl sulfate polyacrylamide gel electrophoresis

TFA:

Trifluoroacetic acid

References

  1. AmphibiaWeb. Information on amphibian biology and conservation. Species Numbers, 2000. http://www.amphibiaweb.org/amphibian/speciesnums.html. Accessed 13 Nov 2014.

  2. Frost D. Amphibian species of the world 6.0, an online reference. http://research.amnh.org/herpetology/amphibia/index.html. Accessed 10 Nov 2014.

  3. Maciel NM, Collevatti RG, Colli GR, Schwartz EF. Late Miocene diversification and phylogenetic relationships of the huge toads in the Rhinella marina (Linnaeus, 1758) species group (Anura: Bufonidae). Mol Phylogenet Evol. 2010;57(2):787–97.

    Article  PubMed  Google Scholar 

  4. Clarke BT. The natural history of amphibian skin secretions, their normal functioning and potential medical applications. Biol Rev Camb Philos Soc. 1997;72(3):365–79.

    Article  CAS  PubMed  Google Scholar 

  5. Sakate M. Lucas de Oliveira PC. Toad envenoming in dogs: effects and treatment. J Venom Anim Toxins. 2000;6(1):52–62.

    Article  Google Scholar 

  6. Toledo RC, Jared C. Cutaneous granular glands and amphibian venoms. Comp Biochem Phys A. 1995;111(1):1–29.

    Article  Google Scholar 

  7. Lee S, Lee Y, Choi YJ, Han KS, Chung HW. Cyto-/genotoxic effects of the ethanol extract of Chan Su, a traditional chinese medicine, in human cancer cell lines. J Ethnopharmacol. 2014;152(2):372–6.

    Article  PubMed  Google Scholar 

  8. Chow L, Johnson M, Wells A, Dasgupta A. Effect of the traditional Chinese medicines Chan Su, Lu-Shen-Wan, Dan Shen, and Asian ginseng on serum digoxin measurement by Tina-quant (Roche) and Synchron LX system (Beckman) digoxin immunoassays. J Clin Lab Anal. 2003;17(1):22–7.

    Article  CAS  PubMed  Google Scholar 

  9. Chen KK, Koraríková A. Pharmacology and toxicology of toad venom. J Pharm Sci. 1967;56(12):1535–41.

    Article  CAS  PubMed  Google Scholar 

  10. Gomes A, Giri B, Alam A, Mukherjee S, Bhattacharjee P, Gomes A. Anticancer activity of a low immunogenic protein toxin (BMP1) from Indian toad (Bufo melanostictus, Schneider) skin extract. Toxicon. 2011;58(1):85–92.

    Article  CAS  PubMed  Google Scholar 

  11. Zhang L, Nakaya K, Yoshida T, Kuroiwa Y. Induction by bufalin of differentiation of human leukemia cells HL60, U937, and ML1 toward macrophage/monocyte-like cells and its potent synergistic effect on the differentiation of human leukemia cells in combination with other inducers. Cancer Res. 1992;52(17):4634–41.

    CAS  PubMed  Google Scholar 

  12. Zuo-Qing L. Traditional Chinese medicine for primary liver cancer. World J Gastroenterol. 1998;4(4):360–4.

    Google Scholar 

  13. Yeh JY, Huang WJ, Kan SF, Wang PS. Effects of bufalin and cinobufagin on the proliferation of androgen dependent and independent prostate cancer cells. Prostate. 2003;54(2):112–24.

    Article  CAS  PubMed  Google Scholar 

  14. Shen S, Zhang Y, Wang Z, Liu R, Gong X. Bufalin induces the interplay between apoptosis and autophagy in glioma cells through endoplasmic reticulum stress. Int J Biol Sci. 2014;10(2):212–24.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  15. Jiang L, Zhao MN, Liu TY, Wu XS, Weng H, Ding Q, et al. Bufalin induces cell cycle arrest and apoptosis in gallbladder carcinoma cells. Tumour Biol. 2014;35(11):10931–41.

    Article  CAS  PubMed  Google Scholar 

  16. Assis AI, Barbosa JE, Carvalho IF. Anticomplementary fraction from the poisonous secretion of paratoid gland of the toad (Bufo marinus paracnemis Lutz). Experientia. 1985;41(7):940–2.

    Article  PubMed  Google Scholar 

  17. Janeway CA, Travers P, Walport M, Shlomchik MJ. O sistema complemento e a imunidade inata. In: Janeway Junior CA, Shlomchik MJ, Travers P, Walport M, editors. Imunobiologia: o sistema imune na saúde e na doença. 5th ed. Porto Alegre: Artmed; 2002. p. 63–83.

    Google Scholar 

  18. Sarma JV, Ward PA. The complement system. Cell Tissue Res. 2011;343(1):227–35.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  19. Ricklin D, Hajishengallis G, Yang K, Lambris JD. Complement: a key system for immune surveillance and homeostasis. Nat Immunol. 2010;11(9):785–97.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  20. Arantes EC, Prado WA, Sampaio SV, Giglio JR. A simplified procedure for the fractionation of Tityus serrulatus venom: isolation and partial characterization of TsTX-IV, a new neurotoxin. Toxicon. 1989;27(8):907–16.

    Article  CAS  PubMed  Google Scholar 

  21. Laemmli UK. Cleavage of structural protein during the assembly of the head of bacteriophage T4. Nature. 1970;227(5259):680–5.

    Article  CAS  PubMed  Google Scholar 

  22. Harrison RA, Lachmann PJ. Complement technology. In: Weir DM, editor. Handbook of experimental immunology. Oxford, UK: Blackwell Scientific Publications; 1986. p. 39.1–39.49.

    Google Scholar 

  23. Hoffmann LG, Mayer MM. Immune hemolysis and complement fixation. In: Williams CA, Chase MW, editors. Methods in immunology and immunochemistry. New York: Academic; 1977. p. 140.

    Google Scholar 

  24. Paula FS, Kabeya LM, Kanashiro A, de Figueiredo AS, Azzolini AE, Uyemura AS, et al. Modulation of human neutrophil oxidative metabolism and degranulation by extract of Tamarindus indica L. fruit pulp. Food Chem Toxicol. 2009;47(1):163–70.

    Article  CAS  PubMed  Google Scholar 

  25. Boyden S. The chemotactic effect of mixtures of antibody and antigen on polymorphonuclear leukocytes. J Exp Med. 1962;115:453–66.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  26. Zigmond SH, Hirsch JG. Leukocyte locomotion and chemotaxis: new methods for evaluation, and demonstration of a cell-derived chemotactic factor. J Exp Med. 1973;137(2):387–410.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  27. Clarke HG, Freeman T. Quantitative immunoelectrophoresis of human serum proteins. Clin Sci. 1968;35(2):403–13.

    CAS  PubMed  Google Scholar 

  28. Deppisch R, Schmitt V, Bommer J, Hänsch GM, Ritz E, Rauterberg EW. Fluid phase generation of terminal complement complex as a novel index of bioincompatibility. Kidney Int. 1990;37(2):696–706.

    Article  CAS  PubMed  Google Scholar 

  29. Marongio AFQ. Isolamento e caracterização parcial do componente do veneno de Bufo paracnemis com ação sobre o sistema complemento. Estudos in vitro e atividade proteolítica. MSs Thesis. São Paulo: Faculdade de Ciências Farmacêuticas de Ribeirão Preto – USP; 2006.

    Google Scholar 

  30. Bertazzi DT, de Assis-Pandochi AI, Talhaferro VL. Caleiro Seixas Azzolini AE, Pereira Crott LS, Arantes EC. Activation of the complement system and leukocyte recruitment by Tityus serrulatus scorpion venom. Int Immunopharmacol. 2005;5(6):1077–84.

    Article  CAS  PubMed  Google Scholar 

  31. Bertazzi DT, de Assis-Pandochi AI, Azzolini AE, Talhaferro VL, Lazzarini M, Arantes EC. Effect of Tityus serrulatus scorpion venom and its major toxin, TsTX-I, on the complement system in vivo. Toxicon. 2003;41(4):501–8.

    Article  CAS  PubMed  Google Scholar 

  32. Farsky SH, Gonçalves LR, Gutíerrez JM, Correa AP, Rucavado A, Gasque P, et al. Bothrops asper snake venom and its metalloproteinase BaP-1 activate the complement system. Role in leucocyte recruitment. Mediators Inflamm. 2000;9(5):213–21.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  33. Silva WDD, Tambourgi DV, Campos ACMR, Magnoli F, Petricevich VL, Kipnis TL. Complement activation by animal venoms. Toxin Rev. 1995;14(3):375–400.

    Article  Google Scholar 

  34. Gebel HM, Finke JH, Elgert KD, Campbell BJ, Barret JT. Inactivation of complement by Loxosceles reclusa spider venom. Am J Trop Med Hyg. 1979;28(4):756–62.

    CAS  PubMed  Google Scholar 

  35. Tanaka GD, Pidde-Queiroz G, Furtado MFD, van den Berg C, Tambourgi DV. Micrurus snake venoms activate human complement system and generate anaphylatoxins. BMC Immunol. 2012;13:4.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Lai R, Zhao Y, Yang D-M, Zha H, Lee W, Zhang Y. Comparative study of the biological activities of the skin secretions from six common Chinese amphibians. Zool Res. 2000;23(2):113–9.

    Google Scholar 

  37. Farsky SH, Walber J, Costa-Cruz M, Cury Y, Teixeira CF. Leukocyte response induced by Bothrops jararaca crude venom: in vivo and in vitro studies. Toxicon. 1997;35(2):185–93.

    Article  CAS  PubMed  Google Scholar 

  38. Tambourgi DV, Dos Santos MC, Furtado MFD, de Freitas MCW, da Silva WD, Kipnis TL. Pro-inflammatory activities in elapid snake venoms. Br J Pharmacol. 1994;112(3):723–7.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  39. Chen T, Shaw C. Identification and molecular cloning of novel trypsin inhibitor analogs from the dermal venom of the Oriental fire-bellied toad (Bombina orientalis) and the European yellow-bellied toad (Bombina variegata). Peptides. 2003;24(6):873–80.

    Article  CAS  PubMed  Google Scholar 

  40. Conlon JM, Kim JB. A protease inhibitor of the Kunitz family from skin secretions of the tomato frog, Dyscophus guineti (Microhylidae). Biochem Biophys Res Commun. 2000;279(3):961–4.

    Article  CAS  PubMed  Google Scholar 

  41. Mignogna G, Pascarella S, Wechselberger C, Hinterleitner C, Mollay C, Amiconi G, et al. BSTI, a trypsin inhibitor from skin secretions of Bombina bombina related to protease inhibitors of nematodes. Protein Sci. 1996;5(2):357–62.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  42. Chakrabarty AK, Saha K, Sen P, Sharma KK, Agarwal SK. Effects of antibacterial agents on the complement system. Immunopharmacology. 1981;3(4):281–7.

    Article  CAS  PubMed  Google Scholar 

  43. Roos A, Nauta AJ, Broers D, Faber-Krol MC, Trouw LA, Drijfhout JW, et al. Specific inhibition of the classical complement pathway by C1q-binding peptides. J Immunol. 2001;167(12):7052–9.

    Article  CAS  PubMed  Google Scholar 

  44. Furlong S, Dutta A, Coath M, Gormley J, Hubbs S, Lloyd D, et al. C3 activation is inhibited by analogs of compstatin but not by serine protease inhibitors or peptidyl alpha-ketoheterocycles. Immunopharmacology. 2000;48(2):199–212.

    Article  CAS  PubMed  Google Scholar 

  45. Soulika AM, Holland MC, Sfyroera G, Sahu A, Lambris JD. Compstatin inhibits complement activation by binding to the beta-chain of complement factor 3. Mol Immunol. 2006;43(12):2023–9.

    Article  CAS  PubMed  Google Scholar 

  46. Gutiérrez JM, Romero M, Diaz C, Borkow G, Ovadia M. Isolation and characterization of a metalloproteinase with weak hemorrhagic activity from the venom of the snake Bothrops asper (terciopelo). Toxicon. 1995;33(1):19–29.

    Article  PubMed  Google Scholar 

  47. Rus V, Malide D, Bolosiu HD, Parasca I, Dutu AL. Levels of SC5b--9 complement complex in plasma and synovial fluid of patients with rheumatic disease. Med Interne. 1990;28(4):305–10.

    CAS  PubMed  Google Scholar 

  48. Harboe M, Ulvund G, Vien L, Fung M, Mollnes TE. The quantitative role of alternative pathway amplification in classical pathway induced terminal complement activation. Clin Exp Immunol. 2004;138(3):439–46.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  49. Harboe M, Garred P, Karlstrom E, Lindstad JK, Stahl GL, Mollnes TE. The down-stream effects of mannan-induced lectin complement pathway activation depend quantitatively on alternative pathway amplification. Mol Immunol. 2009;47(2–3):373–80.

    Article  CAS  PubMed  Google Scholar 

  50. Vogel CW, Finnegan PW, Fritzinger DC. Humanized cobra venom factor: structure, activity, and therapeutic efficacy in preclinical disease models. Mol Immunol. 2014;61(2):191–203.

    Article  CAS  PubMed  Google Scholar 

  51. Mao YF, Yu QH, Zheng XF, Liu K, Liang WQ, Wang YW, et al. Pre-treatment with Cobra venom factor alleviates acute lung injury induced by intestinal ischemia-reperfusion in rats. Eur Rev Med Pharmacol Sci. 2013;17(16):2207–17.

    CAS  PubMed  Google Scholar 

  52. Vogel CW, Fritzinger DC, Gorsuch WB, Stahl GL. Complement depletion with humanised cobra venom factor: efficacy in preclinical models of vascular diseases. Thromb Haemost. 2015;113(3):548–52.

    Article  PubMed  Google Scholar 

  53. Vogel CW, Fritzinger DC. Cobra venom factor: structure, function, and humanization for therapeutic complement depletion. Toxicon. 2010;56(7):1198–222.

    Article  CAS  PubMed  Google Scholar 

  54. Vogel CW, Fritzinger DC, Hew BE, Thorne M, Bammert H. Recombinant cobra venom factor. Mol Immunol. 2004;41(2–3):191–9.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

The authors are grateful to the National Council for Scientific and Technological Development (CNPq), the State of São Paulo Research Foundation (FAPESP – scholarship to FAPA, n. 2008/56327–0) and the Coordination for the Improvement of Higher Education Personnel (CAPES – scholarship to FAPA) and the Support Nucleus for Research on Animal Toxins (NAP-TOXAN-USP, grant n. 12–125432.1.3) for financial support. The authors also acknowledge the biologist Luiz Henrique Anzaloni Pedrosa for extracting the toad venom. Thanks are also due to the Center for the Study of Venoms and Venomous Animals (CEVAP) of UNESP for enabling the publication of this special collection (CNPq process 469660/2014–7).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Eliane C. Arantes.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

FAPA conducted the experiments, performed data analysis and wrote the manuscript. FPL, AECSA and JCP assisted in the experiments. KCFB assisted in the experiments and contributed to the writing of the manuscript. LSPC participated in the research design and contributed to the writing of the manuscript. ECA is the corresponding author and also participated in the research design and contributed to the writing of the manuscript. All authors read and approved the final manuscript.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Anjolette, F.A.P., Leite, F.P., Bordon, K.C.F. et al. Biological characterization of compounds from Rhinella schneideri poison that act on the complement system. J Venom Anim Toxins Incl Trop Dis 21, 25 (2015). https://doi.org/10.1186/s40409-015-0024-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40409-015-0024-9

Keywords