We describe a screening system’s workflow and data analysis for evaluating chemical compound toxicity based on the zebrafish embryo vibration startle response. The system records the movements of zebrafish embryos upon exposure to a vibration stimulus and allows for an integrated evaluation of general toxicity/lethality and neuromuscular toxicity.
We developed a simple screening system for the evaluation of neuromuscular and general toxicity in zebrafish embryos. The modular system consists of electrodynamic transducers above which tissue culture dishes with embryos can be placed. Multiple such loudspeaker-tissue culture dish pairs can be combined. Vibrational stimuli generated by the electrodynamic transducers induce a characteristic startle and escape response in the embryos. A belt-driven linear drive sequentially positions a camera above each loudspeaker to record the movement of the embryos. In this way, alterations to the startle response due to lethality or neuromuscular toxicity of chemical compounds can be visualized and quantified. We present an example of the workflow for chemical compound screening using this system, including the preparation of embryos and treatment solutions, operation of the recording system, and data analysis to calculate benchmark concentration values of compounds active in the assay. The modular assembly based on commercially available simple components makes this system both economical and flexibly adaptable to the needs of particular laboratory setups and screening purposes.
In recent years, zebrafish have become highly popular model organisms for the evaluation of chemical compound effects, encompassing research areas from drug development to environmental toxicology1. As vertebrates, zebrafish share many aspects of their genetic makeup and overall physiology with humans2,3. Therefore, results obtained in this model often are directly relevant to human health. Several drug candidates currently in clinical trials have been identified in compound screens using zebrafish4.
Toxicity assessment is one major application where tests using zebrafish embryonic stages are of interest. Various Organisation for Economic Co-operation and Development (OECD) test guidelines exist for the use of zebrafish in environmental toxicity testing5,6. The small size and rapid development of zebrafish embryos make them highly suitable for screening approaches on a medium to high throughput scale1,3,4. Toxicological endpoints targeted by such screens include embryonic malformations and lethality7, endocrine disruption8, organ toxicity9, and behavioral assessments indicative of neural toxicity10,11. The behavioral assays are possible because zebrafish embryos show various types of locomotor responses to different stimuli depending on their stage of development. For example, 1 day post fertilization (dpf) embryos show spontaneous tail coiling12 and respond to a sequence of light pulses with a typical sequence of movements, the so-called photomotor response (PMR)10. After hatching, typically occurring around 48-72 hours post fertilization (hpf), the freely swimming eleutheroembryos13 gradually develop startle and escape responses to vibrational stimuli starting around 4 dpf14. These responses are characterized by a distinctive bend to the direction opposite the direction of the stimulus (the so-called C-bend or C-start), which is followed by a smaller counter bend and swimming behavior14,15,16,17. Notably, embryonic behaviors are governed by neural circuits using various neurotransmitter systems, allowing for probing chemical compound effects targeting these systems. For example, the PMR assay revealed the effects of compounds interfering with cholinergic, adrenergic, and dopaminergic signaling10, while the startle response involves cholinergic, glutamatergic, and glycinergic neurons16,18. Furthermore, compounds that damage the muscles or the neuro-muscular interface will also affect these behaviors, as will compounds toxic to the inner ear/lateral line hair cells19,20. Observing zebrafish locomotor behavior in response to a stimulus is thus a suitable means to assess not only neurotoxicity but equally ototoxicity and myotoxicity. Scoring locomotor behavior also serves as a proxy for general toxicity/lethality assessment since dead embryos do not move. Thereby, embryonic locomotion behaviors represent an integrative readout for a first-tier toxicity screening approach, which indicates lethal and neuromuscular compound effects in one setup. Given that the eleutheroembryos are already capable of metabolizing compounds, the approach may also detect the effects of metabolic transformation products7,21,22. Importantly, zebrafish embryos are not considered as protected life stage under some animal protection legislations until the stage of free feeding after 120 hpf13. Therefore, they are regarded as an alternative to animal toxicity testing.
Figure 1: Vibration startle response system setup. (A) Overview of the system. Plates with embryos exposed to the test compounds are placed on the electrodynamic transducer array ("loudspeakers"). The camera is sequentially moved by the belt-driven linear drive into the recording position above the targeted transducer. (B) Detailed view of the transducer/loudspeaker with tissue culture dish inserted on top. The plates are illuminated from below by an LED light sheet at 4000-5000 lux. An LED light next to the speaker lights up while the stimulus is given. (C) Still image of video recorded by the camera upon stimulation of the embryos. (D) Screenshot of the configuration file. (E) Screenshot of the recording software interface. Please click here to view a larger version of this figure.
Here, we describe a testing protocol for the evaluation of compound effects on the vibration startle response using an in-house build simple testing device based on vibration stimuli generated by electrodynamic transducers coupled with an automated video recording of several freely moving embryos in a tissue culture dish23. The system is modular and allows for sequential recording from several tissue culture dishes in parallel. In the setup currently used, five electrodynamic transducers provide a vibrational stimulus (500 Hz, duration 1 ms) to tissue culture dishes containing 20 embryos placed on top of them (Figure 1). The plates are illuminated from below at 4000-5000 lux with LED light sheets. An LED light next to each transducer indicates periods of stimulus application, and an oscilloscope indicates waveforms and frequency of the applied stimulus (for details, see Ref. 23). The behavior of the embryos is recorded by a high-speed camera (Table of Materials) at 1000 frames per second (fps), which is moved above the targeted speaker by a belt-driven linear drive. This recording speed is required to reliably resolve the startle response. The system provides a low-cost, individually adaptable alternative to current commercial systems. The precise workflow detailed below is currently performed in the framework of the Precision Toxicology initiative24 in order to determine suitable exposure conditions for OMICS data acquisition from zebrafish embryos treated with a selected set of toxicants.
All zebrafish husbandry and handling were performed in compliance with the German animal protection standards and approved by the Government of Baden-Württemberg, Regierungspräsidium Karlsruhe, Germany (Aktenzeichen 35-9185.64/BH KIT).
1. Preparing stock solutions of chemicals to be tested
2. Collecting and raising zebrafish embryos
3. Preparing working dilution of chemicals to be tested
4. Exposing embryos to the chemicals to be tested
5. Performing the vibration startle assay
NOTE: Analyse the plates in the same order as recorded in step 4.5. Each run should include an E3 control plate.
6. Data analysis
Figure 2A shows the percentage of immotile embryos in 48 clutches of untreated wild-type embryos (AB2O2 strain). On average, 14.33% of untreated wild-type embryos do not react to the vibration stimulus. In 4 clutches, the percentage of immotile larvae reached 50%, but 75% of the clutches had a percentage of immotile larvae below 20%.
Figure 2B,C show an example of a typical calculation of a benchmark concentration/dose (BMC/BMD29,30) for compound effects on motility with the vibration startle assay workflow, as currently performed within the PrecisionTox consortium24. BMC10, BMC25, and BMC50 values correspond to the concentrations at which 10%, 25%, and 50% of embryos show immotility levels higher than the background, respectively. Only embryos that are completely immotile are included in this calculation, not those that still show partial responses, such as only a C-bend without subsequent escape swimming or only tail movements (Figure 2B). The embryos were exposed to 8 concentrations of the sodium channel inhibitor tricaine methanesulfonate, which is frequently used for fish anesthesia31. The data indicate a background level of around 25% immotility in response to the vibration stimulus. Starting at 1% tricaine, motility is reduced and then ceases above 2.5%. The KNIME workflow calculates the BMC50 as 164.9 µM, which corresponds to 1.07% tricaine and an immotility level of 75% (Figure 2C). The small 95% confidence intervals (indicated by the grey shades in the curve) indicate robust reproducibility of the motility values in this assay.
Figure 2D shows an example of a suboptimal assay run, the data of which should not be used for BMC calculations. Five E3 treated control groups with different embryos derived from the same clutch are shown (AB2O2 [wild-type strain] replicate 1-5). Only the first group shows a near normal response, showing around 25% immotility that is consistent with literature values32 and those obtained in the assay described here, as shown in Figure 2A, while all other groups show reduced and/or incomplete behavioral responses (e.g., showing only a C-bend not followed by swimming activity, or motility without a clear C-bend at the beginning). Such a response may occur when embryos do not develop properly and are in an immature state due to developmental delay, which impacts the robustness of the startle response14,33.
Figure 2: Example of a typical result, including benchmark dose calculation. (A) Percentage of non-responsive embryos after the sound pulse for wild-type untreated larvae for 48 clutches (n = 10 per clutch). The mean (14.33%) and standard deviation (±16.19%) are indicated in red. (B) Evaluation of startle response behavior of embryos (n=20 per condition) treated with the indicated concentration of tricaine in E3 medium or with E3 alone as a control. Behavior is classified according to the color scheme and cartoons indicated to the right of the graph, with each embryo assigned to only one of the following classes: "immotile": embryo does not show any movement; "tail moving": embryo shows tail movement, but neither C-bend nor swimming behavior; "motile": embryo shows swimming movement, but no C-bend in response to the vibrational stimulus; "C-bend only": embryo shows C-bend, but not escape swimming; "C-bend + motile": embryo shows typical C-bend behavior followed by escape swimming (the typical full startle response). The different behaviors are shown as a percentage of the total number of embryos for each treatment. (C) BMC calculation graph generated by the KNIME workflow, indicating the percentage of "immotile" embryos for each treatment concentration. Blue, red, and black lines indicate the BMC10, BMC25, and BMC50 values, i.e., the concentrations at which 10%, 25%, and 50% of embryos show immotility levels higher than the background, respectively. (D) Example of a discarded assay run. Five E3-treated control runs with different AB2O2 wild-type embryos derived from the same clutch are shown (replicate 1-5). Only replicate 1 shows a nearly normal response, while embryos of the remaining runs do not show the typical C-bend + escape swimming response. Please click here to view a larger version of this figure.
Table 1: Statistical parameters to estimate the goodness of fit and thresholds to accept determined BMC values. Please click here to download this Table.
Table 2: Properties of a selection of vibrational startle response assay systems. Please click here to download this Table.
Supplementary File 1: Excel template for configuration file. Please click here to download this File.
Supplementary File 2: KNIME input template with an example data set. Please click here to download this File.
Supplementary File 3: KNIME output file example. Please click here to download this File.
We present the workflow and data analysis for chemical compound evaluation using a custom-built zebrafish embryo vibration startle assay setup. The workflow generates robust data that allow the calculation of typical parameters specifying compound toxicity, such as benchmark concentration/dose (BMC/BMD). The modularity of the setup allows adaptation to different needs for throughput and space requirements. As the system is made from low-cost basal components, following a relatively simple setup, it provides a cheap alternative to existing commercial systems, which are generally designed for several assay types at once, rely on proprietary software, and remain relatively costly.
Both these commercial systems and other custom-made systems allow for the assessment of single embryos or larvae in multiwell plates (e.g., 12-well34, 16-well32,35, 24-well20,33,36, 48-well37, 96-well38,39,40,41,42 and even 384-well [as 4×96 well]43), but the spatial restriction in the wells makes the analysis of some data parameters of the escape response (e.g., distance traveled) more challenging. Furthermore, in some of these setups, imaging is restricted to a subset of the wells of the plate, reducing the throughput36,39. Imaging embryos in dishes allows for better assessment of escape response parameters and enables recording the behavior of several embryos at once (up to 30 in a 6 cm dish, for example). Usually, dish-based imaging is limited to one dish per run44,45,46,47,48 (exceptions perform imaging in parallel on 6 dishes with one larva each49 or on 4 larvae in 2 split dishes50), a drawback that can be solved by parallel designs such as in our case. We have summarized some characteristics of the system used in this study and other commercial and custom-made solutions in Table 220,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50,
51,52.
One advantage of the method is a readout capturing both lethality and behavioral changes, which can increase the performance of toxicity assessments. For example, while the zebrafish fish embryo acute toxicity test (FET)5 has been shown to predict toxicity in the adult fish acute toxicity test53 quite well, its prediction accuracy was improved by including behavioral readouts54. The reason for this is the weak mortality induced by neuroactive compounds seen in fish embryos, probably due to the lack of respiratory failure syndrome causing enhanced toxicity in juvenile or adult fish. Neuroactivity can, however, be identified by assessment of behavior. Furthermore, behavioral readouts can also capture myotoxic and ototoxic effects as well as other, more subtle toxic effects on physiology, which are sublethal yet influence the behavioral performance of the organism.
When conducting the assay, it is critical to ensure proper handling of compounds as well as using a homogenously developing batch of zebrafish embryos. Thus, using glass vials for compound storage should minimize the decline in concentrations of chemicals, particularly hydrophobic compounds, due to absorbance to plastic material. In the case of compounds of high absorptive potential to "plastic" polystyrene, glass plates can also be used for the incubation. Cleaning of the eggs in the tissue culture dishes used for collection and removal of dead embryos is a critical step to ensure standard development. Normal speed of development is important, as developmental delays may affect the maturity of neural networks underlying the assessed behavior14,33. Also, to enable comparison of compound effects, eggs should be derived from the same strain since different strains have been reported to present different behavioral profiles38,55,56,57. During exposure, it is important to incubate the embryos in a humidified chamber in order to avoid excessive evaporation of the E3 medium, which would alter the concentrations tested.
E3 controls should be incorporated into each run in order to determine the baseline response level of the particular batch of embryos used in the test series. Typically, we run one plate of controls along each set of 5 measurements. As illustrated in Figure 2D, this approach also allows for the detection of batches with suboptimal responses due to delayed development or for other reasons, such as genetic background effects. In case of an unexpected lack of response to the stimulus, also watch out for potential transducer failure. Typically, the startle responses show a sigmoidal concentration-response behavior that allows for curve fitting using a log-logistic model. However, in rare cases with biphasic responses, other models may have to be employed, such as Gaussian or Cedergreen models. They are available within the R packages drc and bdm27,28.
The lack of response to the vibrational stimulus may indicate simply the death of the embryos or severely impaired life functions due to general cytotoxicity, but might also reflect more specific toxicity targeting neural circuits of stimulus perception, integration, and locomotor output. Other possible compound effects are interference with the neuromuscular interface or with muscle structure and function. To distinguish between these possibilities, further assays are necessary. For example, the structural integrity of the muscles can be assessed with a birefringency assay58,59, and transgenic lines are available to assess perturbance of muscular and neural function60,61. However, the recorded video data already allow for a more detailed analysis of the morphology and the behavioral response of the embryos that can provide first additional information. Is only the C-bend impaired, or all motility? Are still remnants of neuromuscular activity present, as indicated by weak or trembling tail movements? Do such altered behaviors go along with changes in morphology, such as edema or increased body curvature? Additionally, parameters such as latency time until the C-bend or the distance travelled during the escape response can be evaluated (see, for example, Ref. 44).
The screening protocol described here allows for rapid and robust compound toxicity evaluations, with the added value of specifically detecting non-lethal neurotoxic, ototoxic, and myotoxic compounds. The provided analysis workflow is easy to implement and provides a robust readout. Modifications of the stimulus protocols used in the vibration startle assay have been used to address compound effects on more complex aspects of startle behavior as well, such as prepulse inhibition (PPI)39,44 and habituation32,33, and could be adapted to the electrodynamic transducer-based stimulus setup used in this study.
A main application of startle-response-based screening systems is the assessment of compound effects in chemical screens, which is of relevance both for human toxicity evaluation and drug development1,4,62. At the same time, by testing the early life stages of an aquatic organism, the results obtained have direct relevance for ecotoxicological risk assessment63,64. In addition, startle response systems can be used for behavioral phenotyping in genetic screens65,66,67,68,69. Our easily implementable and adaptable system provides an affordable setup to smaller laboratories intending to conduct their own specific screening projects in these various domains of application.
The authors have nothing to disclose.
We thankfully acknowledge the excellent technical assistance of the support staff at the IBCS-BIP fish facility and screening center. This work has received funding from the European Union's Horizon 2020 research and innovation program under Grant Agreement No 965406 (PrecisionTox). This output reflects only the authors' view, and the European Union cannot be held responsible for any use that may be made of the information contained therein.
Fine test sieves, Brass frame, pore size 250 μm | Sigma-Aldrich | Z289744-1EA | Or comparable material |
High-speed camera | XIMEA | MQ013MG-ON USB 3 | |
Laboratory Bottles, Narrow Neck, with Screw Cap | VWR | 215-3261 | Reference number for 50 mL, available up to 20 L. Or comparable material |
Pipette tip, working volume: 10 µL | SARSTEDT | 70.3010.210 | Or comparable material |
Pipette tip, working volume: 1000 µL | SARSTEDT | 70.3050.100 | Or comparable material |
Pipette tip, working volume: 20 µL | SARSTEDT | 70.3020.210 | Or comparable material |
Pipette tip, working volume: 200 µL | SARSTEDT | 70.3030.100 | Or comparable material |
Serological pipette 10 mL | SARSTEDT | 86.1254.001 | Or comparable material |
Serological pipette 25 mL | SARSTEDT | 86.1685.001 | Or comparable material |
Serological pipette 5 mL | SARSTEDT | 86.1253.001 | Or comparable material |
Tissue culture dish 60,0 mm/15,0 mm vented (Polystyrene) | Greiner bio-one | 628102 | Or comparable material |
Tissue culture dish 100, suspension (Polystyrene) | SARSTEDT | 83.3902.500 | Or comparable material |
Transfer pipette 6 mL | SARSTEDT | 86.1175 | Or comparable material |
Tube 15 mL 120 mm x 17 mm PP | SARSTEDT | 62.554.502 | Or comparable material |
Tube 50 mL 114mm x 28 mm PP | SARSTEDT | 62.5472.54 | Or comparable material |