Traumatic brain injuries (TBI) are diverse with heterogeneous injury pathologies, which creates challenges for the clinical treatment and prevention of secondary pathologies such as post-traumatic epilepsy and subsequent dementias. To develop pharmacological strategies that treat TBI and prevent complications, animal models must capture the spectrum of TBI severity to better understand pathophysiological events that occur during and after injury. To address such issues, we improved upon our recent larval zebrafish TBI paradigm emphasizing titrating to different injury levels. We observed coordination between an increase in injury level and clinically relevant injury phenotypes including post-traumatic seizures (PTS) and tau aggregation. This preclinical TBI model is simple to implement, allows dosing of injury levels to model diverse pathologies, and can be scaled to medium- or high-throughput screening.

Traumatic brain injuries (TBI) lead to diverse pathophysiologies such as tissue damage, blood flow abnormalities, metabolic disturbances, oxidative stress, and inflammation (Werner and Engelhard, 2007). TBIs impact tens of millions of people worldwide, with estimates varying greatly due to differences in TBI classification and underreporting of non-hospitalized TBI patients (Langlois et al., 2006; Faul and Coronado, 2015). TBI has devastating effects on children, domestic abuse victims, military personnel and athletes. TBI can result in tragic and long-term consequences such as disability, behavioural and emotional changes, seizures, and an increased risk of developing dementia/neurodegenerative diseases (Whiteneck et al., 2016; Baguley et al., 2006; Perry et al., 2016; Annegers et al., 1998; Gavett et al., 2010; Gardner et al., 2014). We recently established a simple, accessible method to deliver blast TBI to larval zebrafish (Alyenbaawi et al., 2021). The method employs dropping a weight on a closed syringe containing zebrafish in their typical aqueous media, to create blast pressure waves in the syringe (Fig. 1) (Alyenbaawi et al., 2021). This TBI method offers various experimental and bioethical advantages, and we seek to improve upon and share those methods here.

Fig. 1.

Larval zebrafish model of TBI and the modifiable aspects investigated herein. The larval zebrafish TBI model provides an elegantly simple system to create a blast pressure wave and induce TBI. A weight is dropped onto the plunger of a syringe housing dozens of zebrafish larvae in media. The plunger depresses when struck by the falling weight, and this transduces a pressure wave through the liquid media and thus through the zebrafish larvae. In the present study, we modified the following aspects of the larval zebrafish TBI model to increase the blast pressure wave created and allow the user to administer a variety of injury doses: (A) the mass of the weight that was dropped; (B) the height from which the weight was dropped; (C) the size (volume) of the syringe that housed the group of zebrafish larvae. (D) The rigidity of the syringe mount, where a syringe clamp holder had arm movement which dissipates pressure upon impact versus a foam block syringe holder which halted syringe movement upon weight impact. Parameter changes are colour coded throughout the manuscript to match the numbered colours of this figure.

Fig. 1.

Larval zebrafish model of TBI and the modifiable aspects investigated herein. The larval zebrafish TBI model provides an elegantly simple system to create a blast pressure wave and induce TBI. A weight is dropped onto the plunger of a syringe housing dozens of zebrafish larvae in media. The plunger depresses when struck by the falling weight, and this transduces a pressure wave through the liquid media and thus through the zebrafish larvae. In the present study, we modified the following aspects of the larval zebrafish TBI model to increase the blast pressure wave created and allow the user to administer a variety of injury doses: (A) the mass of the weight that was dropped; (B) the height from which the weight was dropped; (C) the size (volume) of the syringe that housed the group of zebrafish larvae. (D) The rigidity of the syringe mount, where a syringe clamp holder had arm movement which dissipates pressure upon impact versus a foam block syringe holder which halted syringe movement upon weight impact. Parameter changes are colour coded throughout the manuscript to match the numbered colours of this figure.

Preclinical TBI models across various species and modes of injury have been productive in advancing our understanding and treatment of neurotrauma. Each model has benefits and challenges. The most prominent research uses rodent TBI models via fluid percussion injury, controlled cortical impact injury, weight-drop impact injury, or blast injury (Dixon et al., 1987, 1991; Marmarou et al., 1994; Cernak et al., 1996). Large animal micropig TBI models provide valuable information about injury events, such as biomechanical aspects of injury that cannot be recreated in small animal models (Ross et al., 1994; Smith et al., 1997; Browne et al., 2011). Drosophila TBI models allow cell-specific genetic manipulations ideal for dissecting injury pathways (Katzenberger et al., 2013; Buhlman et al., 2021). While these animal models are valuable, they are accompanied by various limitations that warrant consideration. For example, mammalian models are not readily accessible to optogenetic measures or manipulations, especially during physical injury, and this has hampered an appreciation of injury-relevant neural activity, neurotransmitters and receptors. Neurotrauma models in Drosophila offer increasing innovations and high throughput, but the drug screening is laborious and the central nervous system (CNS) biology of Drosophila and TBI patients is disparate (e.g. roles for glia in synapse homeostasis, glymphatic system for managing proteostasis, etc.) and requires care in translating the outcomes.

Several adult zebrafish TBI models exist such as stabbed lesion brain injury, closed head pulsed ultrasound injury, blunt force TBI, and laser irradiation TBI (Kishimoto et al., 2012; McCutcheon et al., 2017). In larval zebrafish, TBI has been modelled by stabbed lesion brain injury, linear deceleration induced mild TBI, and the blast pressure TBI model discussed in this study (Alyenbaawi et al., 2021; Gan et al., 2020; Beppi et al., 2022; Hentig et al., 2021; Tikhonova et al., 2022). Adult and larval zebrafish models both have unique strengths and both support translation to adult humans, but here we focus on the benefits of larval zebrafish models.

Zebrafish larvae are a powerful animal model of human disease because they have complex physiological systems akin to humans, making them advantageous compared to in vitro or invertebrate models. The intricate but accessible CNS of zebrafish larvae enables models of neurological diseases such as epilepsy, Parkinson's disease, and tauopathies (Alyenbaawi et al., 2021; Gawel et al., 2020; Feng et al., 2014; Lopez et al., 2017; Paquet et al., 2009; Van Laar et al., 2020; Frøyset et al., 2018). Conserved proteins and biochemistry (Howe et al., 2013; Patton et al., 2021) enable modelling of neurodegeneration and dissection of pathomechanisms (Bai and Burton, 2011; Wu et al., 2016a; Bretaud et al., 2007; Leimer et al., 1999). Larvae can be obtained in high quantities (hundreds per breeding tank) meaning thousands of samples can be run in a variety of assays each week. Zebrafish larvae are small, being several millimetres in length (and growing rapidly), making them amenable to 96-well plate formats for analysis. Zebrafish larvae are amenable to high throughput assays (Rihel and Schier, 2012; Kithcart and MacRae, 2017), which is a powerful untapped opportunity for TBI therapeutic research with the potential of testing small molecule libraries, CRISPR screening of genetic targets, and manipulation of molecular pathways to dissect pathophysiological events. High throughput larval zebrafish models are ideal for discovery of seizure/epilepsy therapeutics, which represent a major target in the form of post-traumatic seizures (PTS) and post-traumatic-epilepsy (PTE). For example, clemizole is in phase II clinical trials for Dravet's syndrome, a genetic epilepsy, which was discovered through assessing a larval zebrafish disease model with high throughput assays (Baraban et al., 2013). Clemizole and other zebrafish-based drug discoveries highlight the translational capabilities and practicality of larval zebrafish models for the discovery of clinical therapeutics (Patton et al., 2021; Goldshtein et al., 2020; Patten et al., 2017; DuVal et al., 2019; Shcheglovitov and Peterson, 2021; Chapela et al., 2019; MacRae and Peterson, 2015; Dash and Patnaik, 2023). Neural activity and other events can also be modelled during injurious impact using larval zebrafish, due to their transparent nature and innovation with calcium imaging paradigms like CaMPARI (Alyenbaawi et al., 2021; Fosque et al., 2015; Kanyo et al., 2020a,b, 2021).

Limitations of larval zebrafish for neurotrauma research include the paucity of robust assays of cognition (learning and memory), a lower husbandry temperature (28°C) compared to human physiology, and a relatively modest toolkit for cell- or tissue-specific genetic modifications (compared to Drosophila or mouse). These are offset by additional benefits such as diurnal habit (contrasting nocturnal rodents and associated logistics of investigating sleep homeostasis, behavioural metrics and glymphatic clearance, etc.). Zebrafish also provide important bioethical considerations within the principles of the three Rs: replacement, reduction and refinement (Russell and Burch, 1959). In our model, zebrafish larvae can serve to replace adult and mammalian animals that would otherwise need to be exposed to neurotraumatic interventions.

We have previously demonstrated that our TBI model produces TBI pathologies by examining dozens of individual larvae after TBI, with individual variability and robust statistics available (Alyenbaawi et al., 2021). These pathologies can be investigated in short timeframes compared to rodent models. We found that neural activity increased during TBI using optogenetics and that PTS and reduced blood flow occurred within hours post-TBI. Neuronal cell death and progressive tau pathology were apparent days after TBI. Excitingly, post-TBI therapeutics were tested using our model and we discovered treating post-TBI neurovolatility (aberrant neural activity such as clinical or non-clinical seizures) (Locskai et al., 2024) with anti-epileptic drugs reduced dementia pathology, whereas increasing neurovolatility with convulsants worsened pathology (Alyenbaawi et al., 2021). Our finding that therapeutically applying the anti-convulsant retigabine prevents TBI-mediated dementia pathology (Alyenbaawi et al., 2021) was validated by others using blast TBI in rodents (Vigil et al., 2023). Surprisingly, our data revealed that one convulsant was unexpectedly protective due to off-target effects which would not have been discoverable in vitro (Alyenbaawi et al., 2021).

Here, we seek to expand upon the model and make the methods openly accessible. Our design principles for the optimization of the larval zebrafish TBI model include: (1) Retaining the method's elegant simplicity, in that the materials required are inexpensive and accessible in most labs, and there are no intellectual property constraints. Any zebrafish researcher can access this TBI model. (2) TBI levels are optimized to trade-off between maximizing injury intensity and larval survival. (3) Various injury intensities can be implemented with reliable consistency to provide a low-cost dose-response paradigm. (4) Animal ethics are paramount, where larval zebrafish fulfil Replacement principles (Russell and Burch, 1959) for adult rodent TBI models; moreover, we present methods allowing researchers to refine (Russell and Burch, 1959) the method using an inexpensive pressure transducer rather than trial-and-error on animals.

We have empirically adjusted components of our TBI method to characterize the pressure dynamics and applied this information to optimize injury for the study of moderate-to-severe TBI. The complex heterogenous nature of TBI is a primary factor for the lack of successful translation of therapeutics in clinical settings (Morales et al., 2005; Marklund and Hillered, 2011). By optimizing our parameters, different TBI intensities can be modelled in our system to further study heterogenous injury pathomechanisms and potential therapeutics. A contemporaneous publication provides step-by-step instructions with set-up and trouble-shooting advice (Gill et al., 2024), and together these writings are intended to allow easy implementation of zebrafish TBI research in any institution. Various TBI model parameters were modified (see Fig. 1) including the mass of the dropped weight, the height the weight was dropped from, the size of the syringe housing the larvae, the rigidity of the syringe holder, and the number of times the weight was dropped per injury group. We chose the ethically favourable route of replacing animals with a custom pressure sensor to optimize pressure levels. Once optimized mechanically, we then validated the outcomes against the biological effects of various injury severity levels, including by measuring PTS, stimuli response, and tau pathology, representing both behavioural and neuroanatomical outcomes of neurotrauma.

Characterising pressure for the optimisation of injury

Each of the four adjustable parameters (denoted as ‘A’ thru ‘D’ in Fig. 1) was combined to provide a broad range of TBI pressures that span three orders of magnitude (Fig. 2). The weight drops created a pressure waveform consisting of multiple pressure waves as the weight bounced (Fig. 2A,B), so we first focused on the maximal pressure of each parameter change individually. Changes to the mass of the weight dropped (as depicted in Fig. 1A) had the largest impact on maximal pressure ranging from a 6- to 21-fold increase across all TBI setups, whereas the drop height of the weight (Fig. 1B) had the second largest effect increasing maximal pressure up to >3-fold (Fig. 2C-F, Tables S1-S4).

Fig. 2.

Maximal pressure created in alternate configurations of the larval zebrafish TBI model. To optimise TBI forces, four equipment parameters were adjusted in concert: (1) a 10 ml or 20 ml syringe was used, and (2) the syringe was held by either a clamp or a rigid foam block. Next, (3) weights of various masses were dropped (4) from various heights (108 cm, 54 cm, 27 cm). Weights ranging from 65-300 g were dropped onto the syringe, and the pressure levels generated for each configuration were measured using Arduino IDE software. Pressure values displayed are the mean values of the maximal pressure reading achieved from three independent weight drop recordings. (A,B) Example time course of pressure dynamics and average pressure generated per pressure wave in the 108 cm/20 ml syringe/foam block TBI setup. These data are duplicated and expanded in Figs. S2-S5 to highlight that additional pressure waves are generated due to the weight bouncing on impact. (C-F) Maximal pressure measurements for each indicated parameter combination. Maximal pressure can be manipulated across three orders of magnitude through changing various TBI assay settings. Alterations to the mass of the weight dropped had the largest impact on maximal pressure, followed by changes to the drop height of the weight. Using a foam block to hold the syringe and decreasing the syringe size both modestly increased maximal pressure.

Fig. 2.

Maximal pressure created in alternate configurations of the larval zebrafish TBI model. To optimise TBI forces, four equipment parameters were adjusted in concert: (1) a 10 ml or 20 ml syringe was used, and (2) the syringe was held by either a clamp or a rigid foam block. Next, (3) weights of various masses were dropped (4) from various heights (108 cm, 54 cm, 27 cm). Weights ranging from 65-300 g were dropped onto the syringe, and the pressure levels generated for each configuration were measured using Arduino IDE software. Pressure values displayed are the mean values of the maximal pressure reading achieved from three independent weight drop recordings. (A,B) Example time course of pressure dynamics and average pressure generated per pressure wave in the 108 cm/20 ml syringe/foam block TBI setup. These data are duplicated and expanded in Figs. S2-S5 to highlight that additional pressure waves are generated due to the weight bouncing on impact. (C-F) Maximal pressure measurements for each indicated parameter combination. Maximal pressure can be manipulated across three orders of magnitude through changing various TBI assay settings. Alterations to the mass of the weight dropped had the largest impact on maximal pressure, followed by changes to the drop height of the weight. Using a foam block to hold the syringe and decreasing the syringe size both modestly increased maximal pressure.

Decreasing syringe size when using the clamp syringe holder (Fig. 1C) decreased maximal pressure up to 0.7-fold for lower drop weight masses but modestly increased (2-fold) maximal pressure at higher drop weight masses (Fig. 2C,E, Table S1 and S2). Decreasing syringe size when using the foam block syringe holder resulted in a 1.62- to 3.05-fold increase in maximal pressure (Fig. 2D,F, Table S3 and S4). Increasing syringe holder rigidity by using the foam block syringe holder instead of the clamp (Fig. 1D) increased maximal pressure 1.29- to 2.98-fold when using a 10 ml syringe (Fig. 2C,D, Table S1 and S3). Comparatively, the rigid foam block syringe holder decreased maximal pressure up to 0.59-fold for 30 g (all heights) and 65 g at 27 cm height but increased pressure up to 2.10-fold otherwise when using a 20 ml syringe (Fig. 2E,F, Table S2 and S4). Overall, maximal pressure was achieved by using the 10 ml syringe, rigid foam block holder, and 300 g weight dropped from 108 cm.

Consistency of pressure values across parameter changes was compared using the relative standard deviation (RSD) of the maximal pressure (Fig. S1). RSD decreased as the mass of the weight increased, except for the 65 g weight which displayed a lower RSD versus the 30 g and 100 g weight (Fig. S1). The lower RSD values for 65 g weights may be because it had a diameter similar to the higher weights, whereas the 35 g and 100 g weights were smaller. The diameter of the weight being mismatched to the diameter of the guide tube used may introduce variability when dropped, perhaps by being more prone to falling at different angles.

Although maximal pressure (discussed above) is a useful metric to compare pressure between TBI setups, secondary pressure waves generated add complexity to the total injurious force created (Fig. 2A,B). We measured the average pressure per wave generated after impact to investigate how parameter changes impact secondary pressure waves. Overall, an increased mass of the weight dropped and decreased syringe size increased the number of additional pressure waves generated (Figs S2-S5, Tables S1-S4).

These secondary pressure waves add complexity when comparing similar maximal pressure measurements (as calculated in Fig. 2) of different TBI setups. For example, 200 g dropped on the 10 ml syringe/foam block generated a primary and secondary pressure wave 210% and 75.4% of the pressure of the initial wave generated by 100 g in the same setup, respectively (Fig. S4B, Table S3). The injurious impact of 200 g would be more than double the impact of 100 g, followed by a second impact on par with the primary impact from 100 g, meaning the pressure experienced by larvae is much greater than that suggested by maximal pressure or the primary waves alone. We suggest that secondary pressure waves should be considered when cross-comparing the injury phenotypes measured in each TBI model setup at similar pressure levels.

Post-traumatic seizures (PTS) increase with injury severity, but severe injury confounds behavioural tracking of seizures

To assess PTS, we used a well-established behavioural assay that quantifies larval movement, using the software EthovisionXT, as a proxy of seizure intensity following TBI. Our goal was to measure how changing our TBI model parameters impacts the pathobiological outputs of neurotrauma injury and how these pathobiological changes are related to pressure measurements.

First, we measured the locomotor activity of larvae injured in the 20 ml syringe and clamp TBI setup, which we have previously used to detect seizures (Alyenbaawi et al., 2021). Overall, TBI increased locomotor activity as the height of the drop weight was increased up to more than 3-fold, but activity decreased when injury became too high (Fig. 3A). Larvae injured with the 100 g and 200 g weight showed hyperactivity versus no TBI larvae as weight drop height increased to 54 cm and 108 cm (Fig. 3B,C, P<0.05-0.0001; effect sizes show doubling or tripling of activity levels and varying with mass of the weight and height). Larvae injured by the 300 g showed hyperactivity at 27 cm and 54 cm weight drop heights, but activity decreased for the 108 cm injury group (Fig. 3D, P<0.05).

Fig. 3.

Seizure-like locomotor activity of 6 dpf larvae following TBI increases at moderate injury levels but decreases at more severe injury levels. Following TBI injury with various parameters, seizure-like activity was assessed by placing the larvae in a 96-well plate recorded using EthovisionXT. Each data point represents an individual zebrafish larva. (A) Exemplar activity traces of individual zebrafish larvae within wells of a 96-well plate. The red depicts the movement pattern of the larvae during one minute of activity and orange boxes denote groups that displayed significantly increased mean activity. (B-D) Mean locomotor activity of larval zebrafish after dropping 100 g (B), 200 g (C), or 300 g (D) from various heights onto a 20 ml syringe held by a clamp compared to a no TBI control. Each experimental group had two to four replicates. Each dot represents behaviour of an individual larva. No TBI control data (N=60) is shared between each graph to prevent bias from splitting the individual control replicates amongst each comparative graph (B) 100 g N=21, 23, 12 for the 27 cm, 54 cm, and 108 cm drop height groups; (C) 200 g N=21, 17, and 12 for the 27 cm, 54 cm, and 108 cm drop height groups; (D) 300 g N=24, 18, and 12 for the 27 cm, 54 cm, and 108 cm drop height groups. (E,F) Linear regression analysis of larval activity as a function of maximal pressure and the average pressure of the primary pressure wave, respectively, for each weight (maximal pressure R2: 0.99 P<0.01 all weights; average pressure of primary pressure wave R2: 0.99 (100 and 200 g) and 0.90 (300 g), P<0.05 (100 g and 200 g). (G) Higher pressure foam block holder induced a decrease in zebrafish locomotor activity after dropping 100 g (N= 48), 200 g (N=48), or 300 g (N=48) from 108 cm onto a 10mL syringe held by a foam block compared to a no TBI control (N=48) (four experimental replicates). (H) Zebrafish locomotor activity as a function of maximal pressure experienced during injury shows a range of 90-300 kPa where behavioural seizures are detected (green) and trend with increased pressure (line), whereas higher pressure induces decreased locomotor activity. A one-way ANOVA with Dunnett's multiple comparisons of means was used to test for statistical significance for behavioural experiments (*=P<0.05, **=P<0.01).

Fig. 3.

Seizure-like locomotor activity of 6 dpf larvae following TBI increases at moderate injury levels but decreases at more severe injury levels. Following TBI injury with various parameters, seizure-like activity was assessed by placing the larvae in a 96-well plate recorded using EthovisionXT. Each data point represents an individual zebrafish larva. (A) Exemplar activity traces of individual zebrafish larvae within wells of a 96-well plate. The red depicts the movement pattern of the larvae during one minute of activity and orange boxes denote groups that displayed significantly increased mean activity. (B-D) Mean locomotor activity of larval zebrafish after dropping 100 g (B), 200 g (C), or 300 g (D) from various heights onto a 20 ml syringe held by a clamp compared to a no TBI control. Each experimental group had two to four replicates. Each dot represents behaviour of an individual larva. No TBI control data (N=60) is shared between each graph to prevent bias from splitting the individual control replicates amongst each comparative graph (B) 100 g N=21, 23, 12 for the 27 cm, 54 cm, and 108 cm drop height groups; (C) 200 g N=21, 17, and 12 for the 27 cm, 54 cm, and 108 cm drop height groups; (D) 300 g N=24, 18, and 12 for the 27 cm, 54 cm, and 108 cm drop height groups. (E,F) Linear regression analysis of larval activity as a function of maximal pressure and the average pressure of the primary pressure wave, respectively, for each weight (maximal pressure R2: 0.99 P<0.01 all weights; average pressure of primary pressure wave R2: 0.99 (100 and 200 g) and 0.90 (300 g), P<0.05 (100 g and 200 g). (G) Higher pressure foam block holder induced a decrease in zebrafish locomotor activity after dropping 100 g (N= 48), 200 g (N=48), or 300 g (N=48) from 108 cm onto a 10mL syringe held by a foam block compared to a no TBI control (N=48) (four experimental replicates). (H) Zebrafish locomotor activity as a function of maximal pressure experienced during injury shows a range of 90-300 kPa where behavioural seizures are detected (green) and trend with increased pressure (line), whereas higher pressure induces decreased locomotor activity. A one-way ANOVA with Dunnett's multiple comparisons of means was used to test for statistical significance for behavioural experiments (*=P<0.05, **=P<0.01).

The 200 g weight dropped from 54 cm versus 300 g dropped from 27 cm created similar pressure dynamics (Table S2, Fig. S3D,F) and both resulted in similar magnitude post-injury hyperactivity. 200 g dropped from 108 cm and 300 g dropped from 54 cm also produced similar pressure dynamics (Table S2, Fig. S3B,D), but the magnitude of hyperactivity for 300 g/54 cm injured larvae was 43.2% lower. Similarly, the pressure dynamics of 100 g dropped from 108 cm and 200 g dropped from 27 cm were comparable (Table S2, Fig. S3B,F), but only the 100 g/108 cm injured fish displayed seizure phenotypes.

Slight differences in pressure-wave dynamics between TBI settings could contribute to larval injury. For example, the number of pressure waves differs for 100 g dropped from 108 cm and 200 g dropped from 27 cm (Fig. S3B,F). Additionally, TBI settings that caused a lower seizure phenotype compared to settings of similar pressure all had a shorter duration of time between the first and second pressure wave (average of 73-86 ms), whereas TBI settings that had similar pressure values and seizure phenotypes only differed by 17 ms on average. These data suggest that small differences between the pressure dynamics may alter phenotypes. Furthermore, the pressure values created by 200 g from 27 cm and 100 g from 108 cm are similar, but the gravitational potential energy of the former is 50% of the latter. Mechanical properties such as syringe and syringe holder movement etc. could potentially impact the injury and resulting phenotypes detected.

Overall, our larval TBI model was consistent across comparable groups and showed an upward trend in seizure activity as the drop height was increased for a given weight. Increases to maximal pressure and average primary wave pressure significantly correlate with the magnitude of detected seizures (Fig. 3E,H,F, P<0.05; R2=0.99). These results suggest that the highest drop height is best for producing injury phenotypes with the largest dynamic range. When comparing TBI assay changes modifying one parameter at a time (e.g. drop height or weight) may minimise variables impacting pressure dynamics and produce the most consistent injury phenotype ‘dose response’.

Locomotor activity decreased at the highest injury level (Fig. 3D), consistent with reduced locomotion TBI symptoms reported in human and animal models (McCutcheon et al., 2017; Beppi et al., 2022; Zhao et al., 2012). In one study consistent with our results, moderately injured rodents had higher locomotor activity, but severely injured rodents had decreased locomotor activity (Chen et al., 2023). Additionally, intense larval zebrafish stage III seizures cause convulsive events and loss of posture (LOP) instead of regular swimming patterns, reducing detectable activity. To test if further increases in pressure and injury resulted in decreased locomotor activity, we administered TBI through our 10mL syringe/foam block setup (Tables S1 and S3). All weight groups had significantly reduced activity (Fig. 3G, P<0.01-0.001). Seizure phenotypes were detectable for injuries of approximately 90-300 kPa maximal pressure (Fig. 3H, green), with additional factors aside from maximal pressure influencing the intensity of seizures as discussed above. injuries above 300 kPa maximal pressure reduced activity (Fig. 3H). These results support that when injury is increased above a certain level larval zebrafish display a reduced locomotor activity phenotype which may confound locomotor-based seizure detection.

To investigate PTS at higher injury levels, we manually quantified seizures. Larvae were injured in the 20 ml syringe clamp or foam block TBI setup (Fig. 4A). Stage I seizures were measured as hyperactive movements with minor convulsions and LOP, while stage II/III seizures were measured as distinct whirlpool motions and major convulsive episodes with distinct bouts of LOP (Baraban et al., 2005; Bandara et al., 2020). The proportion of seizing larvae increased as the mass of the weight used for injury increased for the clamp syringe holder group (Fig. 4A). Larvae injured using the foam block syringe holder displayed a similar proportion of seizing larvae across all weight groups, but more high-intensity seizures were detected (Fig. 4A). Logistic regression analysis indicated the odds of a larva displaying seizure behaviour increased 0.5% for each 1 kPa increase in maximal pressure or 64.7% per 100 kPa (OR per kPa: 1.005, Table 1). Larvae in the clamp/200 g versus foam block/100 g injury groups experienced similar maximal pressures but displayed different phenotypes (Fig. 4A,B), further highlighting that the number of pressure waves and timing between weight bounces may alter the injury phenotype.

Fig. 4.

Post-traumatic seizures and unresponsive zebrafish behaviours increase with more severe injury. To assess decreased locomotor activity following higher injury and pressure intensities, seizures and loss of consciousness were measured manually. (A) Zebrafish larvae were recorded after injury using either the foam block or clamp syringe holder. Larvae were manually scored (blinded) for the presence of seizure behaviour using three categories: no seizure, stage I seizure, or stage II/III seizure. The proportion of larvae displaying seizure trended with increased maximal pressure levels (N=48 no TBI, N=24,100 g clamp, N=24,200 g clamp, N=24,300 g clamp, N=24,100 g foam block, N=19,200 g foam block, and N=22,300 g foam block injured larvae. Each experimental group was repeated twice). (B) A threshold of 1 min or more of inactivity was set to quantify larval inactivity after injury in the videos used for manual seizure scoring. The total duration of these inactive bouts was quantified. Inactivity was significantly higher than no TBI controls for the highest pressure clamp injury group (300 g) and all foam block injury groups (P<0.05-0.001). (C) Following TBI larval zebrafish display a stunned phenotype where they do not respond to stimuli in the form of a fin poke, whereas non-injured fish respond vigorously. The proportion of larvae that exhibit this stunned phenotype increased for higher injury intensity groups versus lower injury groups (purple=highest intensity, N=20; blue=middle intensity, N=13; green=lowest intensity, N=11. Each group had two replicates). Across experiments, dead larvae were removed from quantification. Larvae were verified as dead by a lack of heartbeat in the fin poke response assay and zero movements while floating during video recording. A one-way ANOVA with Dunnett's multiple comparisons of means was used to test for statistical significance between injured larva inactivity and no TBI inactivity (*P<0.05, ** P<0.01, ***P<0.001).

Fig. 4.

Post-traumatic seizures and unresponsive zebrafish behaviours increase with more severe injury. To assess decreased locomotor activity following higher injury and pressure intensities, seizures and loss of consciousness were measured manually. (A) Zebrafish larvae were recorded after injury using either the foam block or clamp syringe holder. Larvae were manually scored (blinded) for the presence of seizure behaviour using three categories: no seizure, stage I seizure, or stage II/III seizure. The proportion of larvae displaying seizure trended with increased maximal pressure levels (N=48 no TBI, N=24,100 g clamp, N=24,200 g clamp, N=24,300 g clamp, N=24,100 g foam block, N=19,200 g foam block, and N=22,300 g foam block injured larvae. Each experimental group was repeated twice). (B) A threshold of 1 min or more of inactivity was set to quantify larval inactivity after injury in the videos used for manual seizure scoring. The total duration of these inactive bouts was quantified. Inactivity was significantly higher than no TBI controls for the highest pressure clamp injury group (300 g) and all foam block injury groups (P<0.05-0.001). (C) Following TBI larval zebrafish display a stunned phenotype where they do not respond to stimuli in the form of a fin poke, whereas non-injured fish respond vigorously. The proportion of larvae that exhibit this stunned phenotype increased for higher injury intensity groups versus lower injury groups (purple=highest intensity, N=20; blue=middle intensity, N=13; green=lowest intensity, N=11. Each group had two replicates). Across experiments, dead larvae were removed from quantification. Larvae were verified as dead by a lack of heartbeat in the fin poke response assay and zero movements while floating during video recording. A one-way ANOVA with Dunnett's multiple comparisons of means was used to test for statistical significance between injured larva inactivity and no TBI inactivity (*P<0.05, ** P<0.01, ***P<0.001).

Table 1.

Logistic regression analysis considering how changes in TBI parameters impact the probability of larval zebrafish displaying TBI phenotypes after injury

TBI phenotypeParameter changeOdds-ratio (β1)95% confidence intervalP valueModel fit (Area under ROC curve)
Seizures Maximal pressure 1.005 (per kPa) 1.003-1.006 <0.0001 0.77 
Stimuli response (30 min post-TBI) Maximal pressure 0.9978 (per kPa) 0.9963-0.9992 <0.01 0.79 
Tau pathology Number of weight drops per injury group 1.099 (per weight drop increase) 1.013-1.194 <0.05 0.81 
Tau pathology Maximal pressure 1.004 (per kPa) 1.003-1.006 <0.0001 0.81 
TBI phenotypeParameter changeOdds-ratio (β1)95% confidence intervalP valueModel fit (Area under ROC curve)
Seizures Maximal pressure 1.005 (per kPa) 1.003-1.006 <0.0001 0.77 
Stimuli response (30 min post-TBI) Maximal pressure 0.9978 (per kPa) 0.9963-0.9992 <0.01 0.79 
Tau pathology Number of weight drops per injury group 1.099 (per weight drop increase) 1.013-1.194 <0.05 0.81 
Tau pathology Maximal pressure 1.004 (per kPa) 1.003-1.006 <0.0001 0.81 

Larvae with more intense seizures appeared to have long bouts of LOP and inactivity, so we next quantified cumulative inactive time across injury groups. Larvae injured in the clamp syringe holder had a non-significant increase in inactive time when injured with 100 g, inactive time comparable to controls when injured with 200 g, and a significant increase in inactive time when injured with 300 g (Fig. 4B, P<0.001). All larval groups injured in the foam block syringe holder displayed a significant increase in inactive time versus controls (Fig. 4B, P<0.05-0.001). These results suggest that as injury intensity increases larval activity decreases even when seizures are apparent, further supporting our behavioural seizure results (Fig. 3D,G).

To further characterize inactivity post-TBI we measured larval fin poke response, an assay used to assess larval loss of consciousness after anesthetization (Du et al., 2018). The average stimuli response rate of larvae 30 min post-TBI decreased as the injury intensity increased (Fig. 4C). Logistic regression analysis of the larval response rate 30 min post-TBI indicated a 0.22% decrease in stimuli response odds per 1 kPa increase in maximal pressure, or 19% per 100 kPa (OR per 1kPa: 0.9978, Table 1). These results suggest that more larvae develop a stunned or unresponsive injury phenotype as injury increases.

These data show that our larval zebrafish TBI assay can be optimized for the detection of PTS at multiple intensities. Altering weight and height parameters impacts behavioural seizure activity between 90-300 kPa, but additional factors influence the magnitude of seizures. To reduce variability across TBI assay settings we recommend changing one parameter at a time. In our experiments, height was associated with increased seizure activity when weight was held constant, and a weight drop height of 108 cm induced the strongest phenotype. We recommend a weight drop height of 108 cm and adjustments to the weight if injury dosing is needed. Since increased injury suppresses activity, other assays like manual seizure scoring, optogenetic assays like CaMPARI, which we have validated during TBI (Alyenbaawi et al., 2021), or EEG may be needed for proper quantification.

Tau aggregation increases with injury levels

Longer-term consequences of TBI include CTE, and a risk for Alzheimer's disease (AD). Thus, we next sought to identify the TBI parameters that produce a useful dynamic range for assessing tauopathy, using a genetically encoded Tau-GFP biosensor. This biosensor approach has been fruitful for assessing tauopathy in vitro (Wu et al., 2016b; Sanders et al., 2014) and was recently deployed in larval zebrafish (Alyenbaawi et al., 2021). Tauopathy is scored in the spinal cord because it is easier to visualize and accurately quantify compared to tauopathy in the brain, and previous work showed the tau abundance in the two structures vary in good coordination (Alyenbaawi et al., 2021). To determine whether increasing injury produces more tauopathy, various parameters were compared, while injury dose was further modulated by increasing the number of times the weight was dropped (Fig. 5). Larval survival sharply declined in the highest-pressure TBI assay so survival was tracked from injury to quantification to highlight the tradeoff of increasing weight drops. Larval survival when injured in the clamp syringe holder was consistent with uninjured controls but starkly decreased when the number of weight drops was increased in the 10 ml foam block setup (Fig. 5E).

Fig. 5.

Tau aggregation increases with injury levels. Tau biosensor larvae were used to quantify how increasing injury intensity impacts the amount of detectable tau pathology. These transgenic larvae express full-length human tau and the human tau 4R domain linked to GFP, allowing for the visualization of tau aggregates as GFP+ puncta. Data represents individual larvae and bars represent +/- SEM and each experimental group was replicated at least twice. (A) Representative images of larval zebrafish spines, where quantifiable GFP+ Tau puncta increase as the severity of injury and number of weight drops increase. Image i depicts the full body of a larva for orientation and scale (scale bar: 0.5 mm). White arrows denote GFP+ Tau puncta. (B) Quantification of GFP+ tau aggregates in larvae injured in a 20mL syringe held by a clamp (N=41, 34, 14, 14, 26, and 24 for the 0, 3, 5, 7, 9 and 11 drops respectively). (C) Quantification of GFP+ tau aggregates in larvae injured in a 10mL syringe held by a clamp (N=27, 12, 13, 12, 19, and 18 for the 0, 3, 5, 7, 9 and 11 drops respectively). (D) Quantification of GFP+ tau aggregates in larvae injured in a 10mL syringe held by a foam block (N=13, 12, 13, 10, 5, and 4 for the 0, 3, 5, 7, 9 and 11 drops respectively). (E) Percent of zebrafish larvae that survived to 7 dpf after TBI using each TBI setup (N=equivalent to B-D). (F) Percent of zebrafish larvae with one or more tau GFP+ puncta in each TBI model setup. (G) Linear regression of mean spinal cord tau puncta (data from Fig. 5B-D) as a function of maximal pressure for larval injury groups that received differing numbers of weight drops. Every weight drop group except the three weight drop group significantly trended with maximal pressure. The five and seven weight drop groups had the highest R2.

Fig. 5.

Tau aggregation increases with injury levels. Tau biosensor larvae were used to quantify how increasing injury intensity impacts the amount of detectable tau pathology. These transgenic larvae express full-length human tau and the human tau 4R domain linked to GFP, allowing for the visualization of tau aggregates as GFP+ puncta. Data represents individual larvae and bars represent +/- SEM and each experimental group was replicated at least twice. (A) Representative images of larval zebrafish spines, where quantifiable GFP+ Tau puncta increase as the severity of injury and number of weight drops increase. Image i depicts the full body of a larva for orientation and scale (scale bar: 0.5 mm). White arrows denote GFP+ Tau puncta. (B) Quantification of GFP+ tau aggregates in larvae injured in a 20mL syringe held by a clamp (N=41, 34, 14, 14, 26, and 24 for the 0, 3, 5, 7, 9 and 11 drops respectively). (C) Quantification of GFP+ tau aggregates in larvae injured in a 10mL syringe held by a clamp (N=27, 12, 13, 12, 19, and 18 for the 0, 3, 5, 7, 9 and 11 drops respectively). (D) Quantification of GFP+ tau aggregates in larvae injured in a 10mL syringe held by a foam block (N=13, 12, 13, 10, 5, and 4 for the 0, 3, 5, 7, 9 and 11 drops respectively). (E) Percent of zebrafish larvae that survived to 7 dpf after TBI using each TBI setup (N=equivalent to B-D). (F) Percent of zebrafish larvae with one or more tau GFP+ puncta in each TBI model setup. (G) Linear regression of mean spinal cord tau puncta (data from Fig. 5B-D) as a function of maximal pressure for larval injury groups that received differing numbers of weight drops. Every weight drop group except the three weight drop group significantly trended with maximal pressure. The five and seven weight drop groups had the highest R2.

The lowest maximal pressure 20 ml syringe held by clamp setup resulted in the least amount of tauopathy, with mean tau puncta counts and the percentage of larvae with tau remaining stagnant across weight drop groups (Fig. 5A(iv), B, P<0.01-0.001 versus control, 5F (green)). Larvae injured in the higher pressure 10 ml syringe/clamp setup had a 1.7-to 4.31-fold increase in tauopathy versus the 20 ml syringe/clamp group (Fig. 5A(v), C). Five or more weight drop events increased detectable tau and the percentage of larvae with tau versus the three-weight-drop group [Fig. 5C, P<0.05-0.01 versus control, Fig. 5F (pink)]. The highest pressure 10 ml syringe/foam block TBI setup resulted in a 2.83−9.96- and 2.31-3.71-fold increase in tauopathy versus the 20 ml and 10 ml syringe held by clamp groups, respectively (Fig. 5A(iii), A(vi), D). The mean tau count and percentage of larvae with tau trended with increased weight drops but decreased at the highest weight drop levels [Fig. 5D, P<0.001 versus control, Fig. 5F (brown)]. The decrease in the percentage of tau GFP+ larvae, mean puncta counts and high variability for nine- and 11-weight-drop groups is most likely because of high lethality (Fig. 5E). Increasing weight-drop events can steeply increases mortality at higher pressure levels, therefore slightly lower injury levels may be best to optimize detectable tauopathy.

Multiple logistic regression analysis of these data predicted that each weight drop increased the odds of larvae having tauopathy by 9.9% and each 1 kPa increase in maximal pressure increased the odds by 0.4%, or 49% per 100 kPa (OR per hit: 1.099, OR per 1kPa: 1.004, Table 1). Linear regression showed tauopathy in the five and seven weight drop groups scaled with maximal pressure the most (Fig. 5G), suggesting five to seven weight-drop events may be optimal, especially when considering optimal survival rates (Fig. 5E). These results show that higher injury intensities and more injury events increase tauopathy in larval zebrafish, much like what is reported in humans, further validating the relevance of our assay.

In this study, we characterized the pressure dynamics in our larval TBI model to optimize model settings for the study of moderate-severe TBI. Modified TBI parameters like syringe size, drop weight mass, drop weight height, and syringe holder rigidity, all alter the blast pressure and injury produced (Fig. 1). TBI pathophysiology is heterogenous between and within injury severity groups, yet current TBI classifications used in research treat TBI severity groups as a singular disease processes (Rosenbaum and Lipton, 2012). For example, most mild TBI patients often recover completely after injury but 30% of patients suffer from persistent debilitating pathology (Alexander, 1995), which underlines differences in pathophysiological progression and outcome for TBI patients. Studying multiple injury severities of TBI is important for the successful translation of preclinical research to clinically significant findings (Xiong et al., 2013). Our data shows that altering parameters of our TBI assay reliably changed the blast pressure created which in turn correlated with increased injury levels as shown by the markers we tested.

Our study addresses gaps in available TBI models in two primary ways. (1) Our model can be used as a high-throughput vertebrate assay with hundreds to thousands of zebrafish larvae undergoing TBI in a day. The number of individual specimens that can be studied increases the amount of heterogenous TBI pathology one would encounter in comparison to other animal models of TBI. Additionally, the high throughput nature of our TBI assay makes it an excellent system to easily test large numbers of potential clinically relevant pharmaceuticals. (2) Our model has bioethical advantages functioning as a clinically relevant replacement (Russell and Burch, 1959) alternative to higher sentience vertebrate models. Reduction strategies are also available with our model where pressure measurements can be used for pre-experiment optimization without the need for trial-and-error use of animals.

Novel insights of the larval zebrafish TBI model

Our results display a strong relationship between impact intensity and multiple injury phenotypes. The ability to reliably create a range of injury phenotypes through altering injury intensity strengthens our model's ability for studying heterogenous TBI outcomes. For example, our model could be used to establish a variety of injury levels and post-TBI phenotypes followed by transcriptomic or proteomic analysis of thousands of individual larvae to better understand pleiotropic mechanisms of TBI pathophysiology and outcomes.

We have also revealed that this TBI method introduces a non-linearity that results from changes in weight bounce, weight fall, the wobble of the syringe stand etc., suggesting that small forces can have an impact on injury. Careful experimental and assay design is needed when creating such TBI assays to ensure the similarity of injury levels. Furthermore, studies into the mechanical events behind injury induction may provide new information about how variability in injury events can impact the biological manifestation of injury phenotypes. We additionally found a biphasic relationship between injury intensity and locomotor seizure phenotypes. These results suggest that high-throughput examination of behavioural seizures require artificial intelligence (AI) to detect severe seizures as injury increases, meaning that further refinement of available tracking systems can expand the high-throughput capabilities of PTS research in larval zebrafish.

Optimal TBI settings for the study of moderate/severe TBI

The optimal dosing of injury within an experiment is attained by using the foam block and 10mL syringe with a large drop height. These settings combine to provide the broadest range of TBI pressures (spanning 33-1105 kPa maximal pressure) and injury severity to the larvae by altering the mass of the weight dropped. Additionally, we suggest adjusting the guide tube diameter to match the diameter of the weight being used to reduce trial variability. If larval injury needs to be increased further than the maximum injury level obtained using the highest weight in the TBI setup, the number of weight drops per injury group can be increased to exacerbate injury, but researchers should be cognizant of the steep decline in larval survival (Russell and Burch, 1959) (Fig. 5E). If the goal of the researcher is to establish an injury ‘dose-response’, maximal pressure measurements are an appropriate metric when modulating pressure through changing an individual parameter within the same TBI setup. When multiple TBI settings are changed at once the pressure dynamics created within the system can change in multiple ways (propensity for second pressure waves, time between weight bounces on the syringe plunger etc.) meaning the maximal pressure values alone may not be the best comparison.

Characteristics of moderate-severe TBI in the larval zebrafish TBI model

TBI is categorized as mild, moderate, and severe based on criteria such as structural imaging, loss-of-consciousness (LOC), amnesia, and Glasgow Coma Scale score (van Baalen et al., 2003; Group MoCmW, 2009). Some of these criteria cannot be measured in larval zebrafish but other metrics can be used to understand TBI severity at different pressure levels. In our injury groups, we saw that larvae experienced a stunned phenotype, seizures, and dementia pathology.

After injury, larvae are prone to falling on their side for long periods (up to hours for intense injury) and exhibit a suppressed response to stimuli, suggesting compromised neuronal function. Such injury phenotypes may be akin to a larval zebrafish LOC-like state that increases with injury severity, similar to how LOC increases with injury level in humans (Group MoCmW, 2009). Further studies are needed to characterize this stunned phenotype as injury increases, but at this stage, the stunned phenotype is important to consider when designing experiments (especially behavioural) and could differentiate TBI severity in larval zebrafish.

PTS occurs in 5-7% of patients hospitalized with TBI, with injury severity being a primary risk factor, increasing to 11% and 35-50% for severe non-penetrating and severe penetrating TBI respectively (Yablon, 1993; Annegers and Coan, 2000; Teasell et al., 2007). Our larval zebrafish showed no seizure activity for low-intensity injuries, increasing injury further caused low-grade seizure responses like hyperactivity and whirlpool swimming, and in the most severely injured larvae, intense convulsive seizures were detected. The number of larvae experiencing seizure phenotypes also increased with the injury severity. Locomotor activity decreased with higher injury intensities, as reported in other TBI studies (McCutcheon et al., 2017; Zhao et al., 2012), leading to decreased seizure detection via behavioural software. A combination of manual seizure scoring, calcium imaging (Alyenbaawi et al., 2021), and EEG assays may be of high value when detecting PTS after severe injuries in larval zebrafish.

PTS models are essential for developing pharmaceutical interventions for PTS and post-traumatic epilepsy (PTE). Acute PTS (seizure(s) <7 days post-TBI) are effectively treated with anti-seizure drugs and can often be prevented with early prophylactic intervention, whereas PTE (1 or more seizure >7 days post-TBI) cannot be prevented through prophylactic anti-seizure medication and 30-40% of epilepsy patients experience seizures that are drug resistant (Rao and Parko, 2015). We have previously demonstrated that anti-seizure medications are successful in treating/preventing pathology following TBI in our larval zebrafish model (Alyenbaawi et al., 2021). Our zebrafish model's similarity to the clinical etiology of TBI pathophysiology and high throughput potential make our paradigm an enticing option for dissecting the mechanisms of pharmaceutical-resistant PTS/PTE and testing potential anti-epileptogenic drugs.

TBI is an important factor of dementia as it is the primary cause of chronic traumatic encephalopathy (CTE) and a major risk factor AD (Plassman et al., 2000; McKee et al., 2009; Fleminger et al., 2003; Ramos-Cejudo et al., 2018; Schaffert et al., 2018). Repetitive mild TBI leads to tau pathology over years-decades in patients with CTE, but single moderate-severe TBI has also been linked to increased tau pathology (Johnson et al., 2012; Walker et al., 2022). Using our ‘tau biosensor larvae’, which express full-length human tau and Tau4R-GFP (Alyenbaawi et al., 2021), we measured tau pathology after varying TBI intensities to optimize detectable tauopathy. Low-pressure injuries resulted in little tauopathy, and high-pressure level injuries resulted in an abundance of tauopathy, displaying a trend that depended on injury severity. Our previous work quantifying tauopathy in this fashion (focused on the spinal cord as a convenient proxy for tauopathy in the brain) allowed us to identify and characterise various drugs and mechanisms regarding post-traumatic epilepsy (Alyenbaawi et al., 2021). Additionally, here we found that applying additional injury events increased the detectable tauopathy, which is true in humans, further strengthening the validity of our model.

Optimising the detection of tauopathy in our TBI model is important because of the potential to dissect pathomechanisms using our transgenic biosensor larvae. We have previously shown that PTS are a druggable link to tauopathy using these tau biosensor larvae in our TBI model, where blocking PTS with pharmaceuticals reduced dementia pathology (Alyenbaawi et al., 2021). Our recent review details current research connecting TBI to PTS and seizures to dementia, but a prominent gap exists in the literature since very few studies have directly tested the connection between TBI, PTS, and dementia (Locskai et al., 2024). The refinements discussed here in TBI assay settings to detect tauopathy and PTS provides us and others with easier access to a model that can further define the mechanisms between TBI, PTS, and dementias, providing further insight into potential therapeutics.

Novel TBI models that reflect the heterogenous nature of TBI are imperative to understand the pathomechanisms that shape long-term clinical outcomes like post-concussion syndrome, PTS/PTE, and dementias. The larval zebrafish TBI method can be accurately optimized to reflect clinically relevant pathology and phenotypes of TBI. The stunned unresponsive phenotype, PTS, and tauopathy all increase with injury severity. This TBI method provides advantages that complement available TBI methods such as its high-throughput capabilities where potential therapeutics can be tested towards an array of different pathological events. Large numbers of larvae can be injured at once (hundreds to thousands in a day) increasing the scale of individual specimen injury that can be examined. These aspects of the TBI model make it an excellent candidate for the preclinical stage of drug discovery. Lastly, the larval zebrafish TBI model has bioethical Replacement and Reduction advantages. Since larval zebrafish are less sentient than other vertebrae models and initial injury settings can be optimized via an inexpensive pressure transducer without larvae present, reducing the number of animals injured.

Animal ethics

Zebrafish were raised and maintained following protocol AUP00000077 approved by the Animal Care and Use Committee: Biosciences at the University of Alberta, operating under the guidelines of the Canadian Council of Animal Care. The fish were raised and maintained within the University of Alberta fish facility under a 14:10 light: dark cycle at 28°C as previously described (Westerfield, 2000).

Larval zebrafish TBI model

In the present study, we refined our previously described larval zebrafish blast pressure TBI model (Alyenbaawi et al., 2021) (Fig. 1). Detailed step-by-step instructions and setup considerations have recently been made available (Gill et al., 2024). Larvae were placed into a 10mL or 20mL syringe (Becton Dickinson #309604, #302830) with 1 ml of E3 media, which was then closed with a Luer-Lok stopper valve attachment (Cole-Parmer #UZ-30600-00). The syringe was placed in either: (1) a three-prong clamp centred on a support rod and stand (e.g. Thermo Fisher Scientific # S41710) or (2) a more rigid foam block with a hole drilled into the centre matching the syringe size, fashioned from a standard yoga foam block (23×15×7.6 cm dimensions). Once mounted, the syringe plunger was centred underneath a guide tube made from rolled paper with a circumference slightly bigger than the largest weight used (40 mm diameter for a standard 38 mm calibration weight; tube heights measured 27 cm, 54 cm, or 108 cm). A weight (e.g. scale calibration weight) was dropped into the guide tube, creating a pressure wave upon impact. Between weight drops, the syringe stopper was opened to remove air bubbles and reset larvae back into the syringe barrel, and the time between drops was <10 min.

Pressure detection setup

To measure pressure in the syringe, we used an Arduino Uno Rev3 microcontroller board (Arduino #A000066) and an automotive fuel line pressure transducer (AUTEX GSND-0556629788) as described previously (Gill et al., 2024). The Arduino board was connected in a circuit to the pressure transducer and a photoresistor (eBoot EBOOT-RESISTOR-05) using a breadboard (Haraqi ESH-PB-01). The photoresistor was positioned opposite from a light source at the top of the guide tube. When the weight was dropped into the guide tube the plane of light would break, altering the photoresistor voltage and triggering the recording. The pressure transducer then measured for one second after the photoresistor voltage dropped (alternatively the light was kept off if longer recordings were required). The pressure transducer output was collected in the Arduino software (Arduino IDE 2.0.3) with a setting of 2,000,000 baud. The output was converted to pressure (kPa) with the formulas [voltage=(5×output)/1023] where output is the values given by the pressure transducer, [PSI=(voltage – 0.5)×(37.5)], and [kPa=PSI×6.895]. For pressure wave measurements, the first pressure spike reading from baseline until the first measurement back at baseline was considered a single pressure wave.

To validate the transducer calibration, static weights of varying mass were rested on the syringe plunger, and readings of the transducer were compared to expected values. Expected pressure values of each weight at rest are calculated using Pascal's principle: pressure=weight / cross-sectional area of the syringe barrel. Caution may be warranted that our calibration of the pressure used a static weight, and information regarding the pressure transducers’ dynamic performance is not provided herein.

Behavioural detection of post-traumatic (PTS) seizures

To measure PTS, we used a locomotor assay previously described (Alyenbaawi et al., 2021; Kanyo et al., 2020a,b, 2021). In short, larval zebrafish were subjected to TBI at 6 days post fertilization (dpf) using the methods described. Larvae were used at 6 dpf because their CNS is well developed by that age, the larvae are robust to manipulation, and it was convenient to our workflow, but other ages could be useful for other workers. Larvae were then placed into a 96-well plate in 200 µl of E3 media and recorded for 30 min using the behavioural software EthoVisonXT (11.5) 45 min after TBI.

Manual scoring of seizures

Seizures were assessed manually to account for the confounding effect of larval inactivity at higher injury intensities. Larvae were administered TBI of varying levels and then recorded for 2 h using EthovisionXT as described above. A high activity bout threshold was set at 64% of the highest maximal activity movement detected in the uninjured group. This threshold detected abnormally high bouts of activity and when they occurred. The time points of these high-activity bouts were then scored (by blinded observers) for the presence of seizure characteristics. Larval zebrafish seizures are defined by three stages, where stage I seizures are increased activity, stage II seizures are whirlpool motions, and stage III seizures are clonic seizure bouts with periods of loss of posture (LOP) and sinking (Baraban et al., 2005; Bandara et al., 2020). Here, stage I seizures were scored as distinct hyperactive movements with slight convulsive behaviour and minor LOP (LOP <1 s/larvae still moving). Stage II/III seizures were scored as whirlpool movements, clonic seizure events, and LOP (fully on side with inactivity >1 s). Larvae that floated in their well and did not display any movement were considered dead and removed from the data.

Larval zebrafish fin poke response

To further assess larval inactivity after TBI, the stimuli response of larvae was assessed post-TBI using a fin poke response assay. This assay is used to measure loss of consciousness after anesthetization in zebrafish (Du et al., 2018). Larvae at 6 dpf were injured with various TBI intensities and placed in a 10 ml Petri dish by group. The larvae were poked once with a sharp metal point on the posterior fin and observed for response every 5 min for 30 min (blinded). Responsive larvae immediately swim away, whereas unresponsive larvae do not respond. The heartbeats of the larvae were assessed after quantification to ensure no dead larvae were included in the data.

Tau-GFP quantification

To measure the formation of GFP+ tau puncta we used Tg(eno2:Hsa.MAPT_Q244-E372−EGFP)ua3171 larvae (ZFIN ID: ZDB-ALT-211005-6) crossed with Tg(eno2:hsa.MAPT-ires-egfp)Pt406 larvae (ZFIN ID: ZDB-ALT-080122–6), which express full-length 4R human tau and the 4R domain of human tau linked to GFP, both under the eno2 promoter (Alyenbaawi et al., 2021; Bai et al., 2007). These fish have been previously validated to develop tau GFP+ aggregates after TBI (Alyenbaawi et al., 2021). Larvae subjected to TBI (as described above) at 3 dpf were anesthetized at 7 dpf with MS-222 and imaged under a fluorescent dissecting scope (Leica M165 FC). GFP+ puncta were counted from a lateral view (blinded).

Statistics

All statistical analyses were done using GraphPad Prism Software (version 10.0.0). Mean and maximum pressures were calculated +/- standard error from the mean in units of kilopascals (kPa). Relative standard deviations were calculated taking the proportion of the SD to the mean, multiplied by 100. One-way ordinary ANOVA tests were used for statistical analysis of EthovisionXT seizure data followed by Dunnett's multiple comparison of means. Linear relationships were analysed using simple linear regression, whereas binary relationships were analysed with multiple or simple logistic regression.

The authors appreciate the support of Science Animal Support Services at the University of Alberta, Canada.

Author contributions

Conceptualization: L.F.L., A.H.B., H.A., E.A.B., W.A.; Methodology: L.F.L.; Validation: L.F.L.; Formal analysis: L.F.L., S.A.W.T.; Investigation: L.F.L., T.G., S.A.W.T.; Data curation: L.F.L.; Writing - original draft: L.F.L., T.G.; Writing - review & editing: L.F.L., T.G., S.A.W.T., A.H.B., H.A., E.A.B., W.A.; Visualization: L.F.L.; Supervision: E.A.B., W.A.; Project administration: W.A.; Funding acquisition: W.A.

Funding

L.F.L. is supported by studentships from the Canadian Institutes of Health Research (CIHR) and the Alberta Synergies in Alzheimer's and Related Disorders (SynAD) program which is funded by the Alzheimer Society of Alberta and Northwest Territories through their Hope for Tomorrow program and the University Hospital Foundation. H.A. appreciates supports from the Deanship of Scientific Research at Majmaah University. T.G.’s work for this manuscript was supported in part by a studentship from Natural Sciences and Engineering Research Council of Canada (NSERC). The authors appreciate the support to WTA for operating costs from donors who prefer to remain anonymous. EAB was supported by the US Department of Veterans Affairs (I01-BX003168). Open Access funding provided by University of Alberta. Deposited in PMC for immediate release.

Data availability

All relevant data can be found within the article and its supplementary information.

Alexander
,
M. P.
(
1995
).
Mild traumatic brain injury: pathophysiology, natural history, and clinical management
.
Neurology
45
,
1253
-
1260
.
Alyenbaawi
,
H.
,
Kanyo
,
R.
,
Locskai
,
L. F.
,
Kamali-Jamil
,
R.
,
Duval
,
M. G.
,
Bai
,
Q.
,
Wille
,
H.
,
Burton
,
E. A.
and
Allison
,
W. T.
(
2021
).
Seizures are a druggable mechanistic link between TBI and subsequent tauopathy
.
Elife
10
.
Annegers
,
J. F.
and
Coan
,
S. P.
(
2000
).
The risks of epilepsy after traumatic brain injury
.
Seizure
9
,
453
-
457
.
Annegers
,
J. F.
,
Hauser
,
W. A.
,
Coan
,
S. P.
and
Rocca
,
W. A.
(
1998
).
A population-based study of seizures after traumatic brain injuries
.
N. Engl. J. Med.
338
,
20
-
24
.
Baguley
,
I. J.
,
Cooper
,
J.
and
Felmingham
,
K.
(
2006
).
Aggressive behavior following traumatic brain injury: how common is common?
J. Head Trauma Rehabil.
21
,
45
-
56
.
Bai
,
Q.
and
Burton
,
E. A.
(
2011
).
Zebrafish models of Tauopathy
.
Biochim. Biophys. Acta
1812
,
353
-
363
.
Bai
,
Q.
,
Garver
,
J. A.
,
Hukriede
,
N. A.
and
Burton
,
E. A.
(
2007
).
Generation of a transgenic zebrafish model of Tauopathy using a novel promoter element derived from the zebrafish eno2 gene
.
Nucleic Acids Res.
35
,
6501
-
6516
.
Bandara
,
S. B.
,
Carty
,
D. R.
,
Singh
,
V.
,
Harvey
,
D. J.
,
Vasylieva
,
N.
,
Pressly
,
B.
,
Wulff
,
H.
and
Lein
,
P. J.
(
2020
).
Susceptibility of larval zebrafish to the seizurogenic activity of GABA type A receptor antagonists
.
Neurotoxicology
76
,
220
-
234
.
Baraban
,
S. C.
,
Taylor
,
M. R.
,
Castro
,
P. A.
and
Baier
,
H.
(
2005
).
Pentylenetetrazole induced changes in zebrafish behavior, neural activity and c-fos expression
.
Neuroscience
131
,
759
-
768
.
Baraban
,
S. C.
,
Dinday
,
M. T.
and
Hortopan
,
G. A.
(
2013
).
Drug screening in Scn1a zebrafish mutant identifies clemizole as a potential Dravet syndrome treatment
.
Nat. Commun.
4
,
2410
.
Beppi
,
C.
,
Penner
,
M.
,
Straumann
,
D.
and
Bögli
,
S. Y.
(
2022
).
A non-invasive biomechanical model of mild TBI in larval zebrafish
.
PLoS One
17
,
e0268901
.
Bretaud
,
S.
,
Allen
,
C.
,
Ingham
,
P. W.
and
Bandmann
,
O.
(
2007
).
p53-dependent neuronal cell death in a DJ-1-deficient zebrafish model of Parkinson's disease
.
J. Neurochem.
100
,
1626
-
1635
.
Browne
,
K. D.
,
Chen
,
X. H.
,
Meaney
,
D. F.
and
Smith
,
D. H.
(
2011
).
Mild traumatic brain injury and diffuse axonal injury in swine
.
J. Neurotrauma
28
,
1747
-
1755
.
Buhlman
,
L. M.
,
Krishna
,
G.
,
Jones
,
T. B.
and
Thomas
,
T. C.
(
2021
).
Drosophila as a model to explore secondary injury cascades after traumatic brain injury
.
Biomed. Pharmacother.
142
,
112079
.
Cernak
,
I.
,
Savic
,
J.
,
MalicevicZunic
,
G.
,
Radosevic
,
P.
,
Ivanovic
,
I.
and
Davidovic
,
L.
(
1996
).
Involvement of the central nervous system in the general response to pulmonary blast injury
.
J. Trauma
40
Suppl. 3
,
S100
-
S104
.
Chapela
,
D.
,
Sousa
,
S.
,
Martins
,
I.
,
Cristóvão
,
A. M.
,
Pinto
,
P.
,
Corte-Real
,
S.
and
Saúde
,
L.
(
2019
).
A zebrafish drug screening platform boosts the discovery of novel therapeutics for spinal cord injury in mammals
.
Sci. Rep.
9
,
10475
.
Chen
,
Z. P.
,
Wang
,
S.
,
Zhao
,
X.
,
Fang
,
W.
,
Wang
,
Z.
,
Ye
,
H.
,
Wang
,
M.-J.
,
Ke
,
L.
,
Huang
,
T.
,
Lv
,
P.
et al.
(
2023
).
Lipid-accumulated reactive astrocytes promote disease progression in epilepsy
.
Nat. Neurosci.
26
,
542
-
554
.
Dash
,
S. N.
and
Patnaik
,
L.
(
2023
).
Flight for fish in drug discovery: a review of zebrafish-based screening of molecules
.
Biol. Lett.
19
,
20220541
.
Dixon
,
C. E.
,
Lyeth
,
B. G.
,
Povlishock
,
J. T.
,
Findling
,
R. L.
,
Hamm
,
R. J.
,
Marmarou
,
A.
,
Young
,
H. F.
and
Hayes
,
R. L.
(
1987
).
A fluid percussion model of experimental brain injury in the rat
.
J. Neurosurg.
67
,
110
-
119
.
Dixon
,
C. E.
,
Clifton
,
G. L.
,
Lighthall
,
J. W.
,
Yaghmai
,
A. A.
and
Hayes
,
R. L.
(
1991
).
A controlled cortical impact model of traumatic brain injury in the rat
.
J. Neurosci. Methods
39
,
253
-
262
.
Du
,
W. J.
,
Zhang
,
R. W.
,
Li
,
J.
,
Zhang
,
B.-B.
,
Peng
,
X.-L.
,
Cao
,
S.
,
Yuan
,
J.
,
Yuan
,
C.-D.
,
Yu
,
T.
and
Du
,
J.-
(
2018
).
The locus coeruleus modulates intravenous general anesthesia of zebrafish via a cooperative mechanism
.
Cell Rep.
24
,
3146
-
3155.e3
.
Duval
,
M. G.
,
Hinge
,
V. K.
,
Snyder
,
N.
,
Kanyo
,
R.
,
Bratvold
,
J.
,
Pokrishevsky
,
E.
,
Cashman
,
N. R.
,
Blinov
,
N.
,
Kovalenko
,
A.
and
Allison
,
W. T.
(
2019
).
Tryptophan 32 mediates SOD1 toxicity in a in vivo motor neuron model of ALS and is a promising target for small molecule therapeutics
.
Neurobiol. Dis.
124
,
297
-
310
.
Faul
,
M.
and
Coronado
,
V.
(
2015
).
Epidemiology of traumatic brain injury
.
Handb. Clin. Neurol.
127
,
3
-
13
.
Feng
,
C. W.
,
Wen
,
Z. H.
,
Huang
,
S. Y.
,
Hung
,
H.-C.
,
Chen
,
C.-H.
,
Yang
,
S.-N.
,
Chen
,
N.-F.
,
Wang
,
H.-M.
,
Hsiao
,
C.-D.
and
Chen
,
W.-F.
(
2014
).
Effects of 6-hydroxydopamine exposure on motor activity and biochemical expression in zebrafish (Danio rerio) larvae
.
Zebrafish
11
,
227
-
239
.
Fleminger
,
S.
,
Oliver
,
D. L.
and
Lovestone
,
S.
(
2003
).
Head injury as a risk factor for Alzheimer's disease: the evidence 10 years on; a partial replication
.
J. Neurol. Neurosurg. Psychiatry
74
,
857
-
862
.
Fosque
,
B. F.
,
Sun
,
Y.
,
Dana
,
H.
,
Yang
,
C.-T.
,
Ohyama
,
T.
,
Tadross
,
M. R.
,
Patel
,
R.
,
Zlatic
,
M.
,
Kim
,
D. S.
,
Ahrens
,
M. B.
et al.
(
2015
).
Neural circuits. Labeling of active neural circuits in vivo with designed calcium integrators
.
Science
347
,
755
-
760
.
Frøyset
,
A. K.
,
Edson
,
A. J.
,
Gharbi
,
N.
,
Khan
,
E. A.
,
Dondorp
,
D.
,
Bai
,
Q.
,
Tiraboschi
,
E.
,
Suster
,
M. L.
,
Connolly
,
J. B.
,
Burton
,
E. A.
et al.
(
2018
).
Astroglial DJ-1 over-expression up-regulates proteins involved in redox regulation and is neuroprotective in vivo
.
Redox Biol.
16
,
237
-
247
.
Gan
,
D.
,
Wu
,
S.
,
Chen
,
B.
and
Zhang
,
J.
(
2020
).
Application of the Zebrafish traumatic brain injury model in assessing cerebral inflammation
.
Zebrafish
17
,
73
-
82
.
Gardner
,
R. C.
,
Burke
,
J. F.
,
Nettiksimmons
,
J.
,
Kaup
,
A.
,
Barnes
,
D. E.
and
Yaffe
,
K.
(
2014
).
Dementia risk after traumatic brain injury vs nonbrain trauma: the role of age and severity
.
JAMA Neurol
71
,
1490
-
1497
.
Gavett
,
B. E.
,
Stern
,
R. A.
,
Cantu
,
R. C.
,
Nowinski
,
C. J.
and
Mckee
,
A. C.
(
2010
).
Mild traumatic brain injury: a risk factor for neurodegeneration
.
Alzheimers Res. Ther.
2
,
18
.
Gawel
,
K.
,
Langlois
,
M.
,
Martins
,
T.
,
Van Der Ent
,
W.
,
Tiraboschi
,
E.
,
Jacmin
,
M.
,
Crawford
,
A. D.
and
Esguerra
,
C. V.
(
2020
).
Seizing the moment: Zebrafish epilepsy models
.
Neurosci. Biobehav. Rev.
116
,
1
-
20
.
Gill
,
T.
,
Locskai
,
L. F.
,
Burton
,
A. H.
et al.
. (
2024
).
Delivering traumatic brain injury to larval Zebrafish
. In
Zebrafish: Methods and Protocols
, Vol.
2707
(ed.
J. F.
Amatruda
), pp.
3
-
22
.
New York
:
Springer Nature
.
Goldshtein
,
H.
,
Muhire
,
A.
,
Petel Légaré
,
V.
,
Pushett
,
A.
,
Rotkopf
,
R.
,
Shefner
,
J. M.
,
Peterson
,
R. T.
,
Armstrong
,
G. A. B.
and
Russek-Blum
,
N.
(
2020
).
Efficacy of Ciprofloxacin/Celecoxib combination in zebrafish models of amyotrophic lateral sclerosis
.
Ann. Clin. Transl. Neurol.
7
,
1883
-
1897
.
Group MoCmW
(
2009
).
VA/DoD Clinical practice guideline for management of concussion/mild traumatic brain injury
.
J. Rehabil. Res. Dev.
46
,
CP1
-
C68
.
Hentig
,
J.
,
Cloghessy
,
K.
,
Lahne
,
M.
,
Jung
,
Y. J.
,
Petersen
,
R. A.
,
Morris
,
A. C.
and
Hyde
,
D. R.
(
2021
).
Zebrafish blunt-force TBI induces heterogenous injury pathologies that mimic human TBI and responds with sonic hedgehog-dependent cell proliferation across the neuroaxis
.
Biomedicines
9
.
Howe
,
K.
,
Clark
,
M. D.
,
Torroja
,
C. F.
,
Torrance
,
J.
,
Berthelot
,
C.
,
Muffato
,
M.
,
Collins
,
J. E.
,
Humphray
,
S.
,
McLaren
,
K.
,
Matthews
,
L.
et al.
(
2013
).
The zebrafish reference genome sequence and its relationship to the human genome
.
Nature
496
,
498
-
503
.
Johnson
,
V. E.
,
Stewart
,
W.
and
Smith
,
D. H.
(
2012
).
Widespread τ and amyloid-β pathology many years after a single traumatic brain injury in humans
.
Brain Pathol.
22
,
142
-
149
.
Kanyo
,
R.
,
Leighton
,
P. L. A.
,
Neil
,
G. J.
,
Locskai
,
L. F.
and
Allison
,
W. T.
(
2020a
).
Amyloid-β precursor protein mutant zebrafish exhibit seizure susceptibility that depends on prion protein
.
Exp. Neurol.
328
.
Kanyo
,
R.
,
Wang
,
C. K.
,
Locskai
,
L. F.
,
Li
,
J.
,
Allison
,
W. T.
and
Kurata
,
H. T.
(
2020b
).
Functional and behavioral signatures of Kv7 activator drug subtypes
.
Epilepsia
61
,
1678
-
1690
.
Kanyo
,
R.
,
Amin
,
M. R.
,
Locskai
,
L. F.
,
Bouvier
,
D. D.
,
Olthuis
,
A. M.
,
Allison
,
W. T.
and
Ali
,
D. W.
(
2021
).
Medium-throughput zebrafish optogenetic platform identifies deficits in subsequent neural activity following brief early exposure to cannabidiol and Δ
.
Sci. Rep.
11
,
11515
.
Katzenberger
,
R. J.
,
Loewen
,
C. A.
,
Wassarman
,
D. R.
,
Petersen
,
A. J.
,
Ganetzky
,
B.
and
Wassarman
,
D. A.
(
2013
).
A Drosophila model of closed head traumatic brain injury
.
Proc. Natl. Acad. Sci. USA
110
,
E4152
-
E4159
.
Kishimoto
,
N.
,
Shimizu
,
K.
and
Sawamoto
,
K.
(
2012
).
Neuronal regeneration in a zebrafish model of adult brain injury
.
Dis. Model. Mech.
5
,
200
-
209
.
Kithcart
,
A.
and
MacRae
,
C. A.
(
2017
).
Using Zebrafish for high-throughput screening of novel cardiovascular drugs
.
JACC Basic Transl. Sci.
2
,
1
-
12
.
Langlois
,
J. A.
,
Rutland-Brown
,
W.
and
Wald
,
M. M.
(
2006
).
The epidemiology and impact of traumatic brain injury: a brief overview
.
J. Head Trauma Rehabil.
21
,
375
-
378
.
Leimer
,
U.
,
Lun
,
K.
,
Romig
,
H.
,
Walter
,
J.
,
Grünberg
,
J.
,
Brand
,
M.
and
Haass
,
C.
(
1999
).
Zebrafish (Danio rerio) presenilin promotes aberrant amyloid beta-peptide production and requires a critical aspartate residue for its function in amyloidogenesis
.
Biochemistry
38
,
13602
-
13609
.
Locskai
,
L. F.
,
Alyenbaawi
,
H.
and
Allison
,
W. T.
(
2024
).
Antiepileptic drugs as potential dementia prophylactics following traumatic brain injury
.
Annu. Rev. Pharmacol. Toxicol.
64
,
577
-
598
.
Lopez
,
A.
,
Lee
,
S. E.
,
Wojta
,
K.
,
Ramos
,
E. M.
,
Klein
,
E.
,
Chen
,
J.
,
Boxer
,
A. L.
,
Gorno-Tempini
,
M. L.
,
Geschwind
,
D. H.
,
Schlotawa
,
L.
et al.
(
2017
).
A152T tau allele causes neurodegeneration that can be ameliorated in a zebrafish model by autophagy induction
.
Brain
140
,
1128
-
1146
.
MacRae
,
C. A.
and
Peterson
,
R. T.
(
2015
).
Zebrafish as tools for drug discovery
.
Nat. Rev. Drug Discov.
14
,
721
-
731
.
Marklund
,
N.
and
Hillered
,
L.
(
2011
).
Animal modelling of traumatic brain injury in preclinical drug development: where do we go from here?
Br. J. Pharmacol.
164
,
1207
-
1229
.
Marmarou
,
A.
,
Foda
,
M. A.
,
Van Den Brink
,
W.
,
Campbell
,
J.
,
Kita
,
H.
and
Demetriadou
,
K.
(
1994
).
A new model of diffuse brain injury in rats. Part I: Pathophysiology and biomechanics
.
J. Neurosurg.
80
,
291
-
300
.
McCutcheon
,
V.
,
Park
,
E.
,
Liu
,
E.
,
Sobhebidari
,
P.
,
Tavakkoli
,
J.
,
Wen
,
X.-Y.
and
Baker
,
A. J.
(
2017
).
A novel model of traumatic brain injury in adult Zebrafish demonstrates response to injury and treatment comparable with mammalian models
.
J. Neurotrauma
34
,
1382
-
1393
.
McKee
,
A. C.
,
Cantu
,
R. C.
,
Nowinski
,
C. J.
,
Hedley-Whyte
,
E. T.
,
Gavett
,
B. E.
,
Budson
,
A. E.
,
Santini
,
V. E.
,
Lee
,
H.-S.
,
Kubilus
,
C. A.
and
Stern
,
R. A.
(
2009
).
Chronic traumatic encephalopathy in athletes: progressive tauopathy after repetitive head injury
.
J. Neuropathol. Exp. Neurol.
68
,
709
-
735
.
Morales
,
D. M.
,
Marklund
,
N.
,
Lebold
,
D.
,
Thompson
,
H. J.
,
Pitkanen
,
A.
,
Maxwell
,
W. L.
,
Longhi
,
L.
,
Laurer
,
H.
,
Maegele
,
M.
,
Neugebauer
,
E.
et al.
(
2005
).
Experimental models of traumatic brain injury: do we really need to build a better mousetrap?
Neuroscience
136
,
971
-
989
.
Paquet
,
D.
,
Bhat
,
R.
,
Sydow
,
A.
,
Mandelkow
,
E.-M.
,
Berg
,
S.
,
Hellberg
,
S.
,
Fälting
,
J.
,
Distel
,
M.
,
Köster
,
R. W.
,
Schmid
,
B.
et al.
(
2009
).
A zebrafish model of tauopathy allows in vivo imaging of neuronal cell death and drug evaluation
.
J. Clin. Invest.
119
,
1382
-
1395
.
Patten
,
S. A.
,
Aggad
,
D.
,
Martinez
,
J.
,
Tremblay
,
E.
,
Petrillo
,
J.
,
Armstrong
,
G. A. B.
,
La Fontaine
,
A.
,
Maios
,
C.
,
Liao
,
M.
,
Ciura
,
S.
et al.
(
2017
).
Neuroleptics as therapeutic compounds stabilizing neuromuscular transmission in amyotrophic lateral sclerosis
.
JCI Insight
2
,
e97152
.
Patton
,
E. E.
,
Zon
,
L. I.
and
Langenau
,
D. M.
(
2021
).
Zebrafish disease models in drug discovery: from preclinical modelling to clinical trials
.
Nat. Rev. Drug Discov.
20
,
611
-
628
.
Perry
,
D. C.
,
Sturm
,
V. E.
,
Peterson
,
M. J.
,
Pieper
,
C. F.
,
Bullock
,
T.
,
Boeve
,
B. F.
,
Miller
,
B. L.
,
Guskiewicz
,
K. M.
,
Berger
,
M. S.
,
Kramer
,
J. H.
et al.
(
2016
).
Association of traumatic brain injury with subsequent neurological and psychiatric disease: a meta-analysis
.
J. Neurosurg.
124
,
511
-
526
.
Plassman
,
B. L.
,
Havlik
,
R. J.
,
Steffens
,
D. C.
,
Helms
,
M. J.
,
Newman
,
T. N.
,
Drosdick
,
D.
,
Phillips
,
C.
,
Gau
,
B. A.
,
Welsh-Bohmer
,
K. A.
,
Burke
,
J. R.
et al.
(
2000
).
Documented head injury in early adulthood and risk of Alzheimer's disease and other dementias
.
Neurology
55
,
1158
-
1166
.
Ramos-Cejudo
,
J.
,
Wisniewski
,
T.
,
Marmar
,
C.
,
Zetterberg
,
H.
,
Blennow
,
K.
,
De Leon
,
M. J.
and
Fossati
,
S.
(
2018
).
Traumatic brain injury and Alzheimer's disease: the cerebrovascular link
.
EBioMedicine
28
,
21
-
30
.
Rao
,
V. R.
and
Parko
,
K. L.
(
2015
).
Clinical approach to posttraumatic epilepsy
.
Semin. Neurol.
35
,
57
-
63
.
Rihel
,
J.
and
Schier
,
A. F.
(
2012
).
Behavioral screening for neuroactive drugs in zebrafish
.
Dev. Neurobiol.
72
,
373
-
385
.
Rosenbaum
,
S. B.
and
Lipton
,
M. L.
(
2012
).
Embracing chaos: the scope and importance of clinical and pathological heterogeneity in mTBI
.
Brain Imaging Behav.
6
,
255
-
282
.
Ross
,
D. T.
,
Meaney
,
D. F.
,
Sabol
,
M. K.
,
Smith
,
D. H.
and
Gennarelli
,
T. A.
(
1994
).
Distribution of forebrain diffuse axonal injury following inertial closed head injury in miniature swine
.
Exp. Neurol.
126
,
291
-
299
.
Russell
,
W. M. S.
and
Burch
,
R. L
. (
1959
).
The Principles of Humane Experimental Technique
.
Methuen, London
.
Sanders
,
D. W.
,
Kaufman
,
S. K.
,
Devos
,
S. L.
,
Sharma
,
A. M.
,
Mirbaha
,
H.
,
Li
,
A.
,
Barker
,
S. J.
,
Foley
,
A. C.
,
Thorpe
,
J. R.
,
Serpell
,
L. C.
et al.
(
2014
).
Distinct tau prion strains propagate in cells and mice and define different tauopathies
.
Neuron
82
,
1271
-
1288
.
Schaffert
,
J.
,
Lobue
,
C.
,
Didehbani
,
N.
,
Chiang
,
H.-S.
,
Didehbani
,
N.
,
Lacritz
,
L.
,
Rossetti
,
H.
,
Dieppa
,
M.
,
Hart
,
J.
and
Cullum
,
C. M.
(
2018
).
Traumatic brain injury history is associated with an earlier age of dementia onset in autopsy-confirmed Alzheimer's disease
.
Neuropsychology
32
,
410
-
416
.
Shcheglovitov
,
A.
and
Peterson
,
R. T.
(
2021
).
Screening platforms for genetic epilepsies-zebrafish, iPSC-derived neurons, and organoids
.
Neurotherapeutics
18
,
1478
-
1489
.
Smith
,
D. H.
,
Chen
,
X. H.
,
Xu
,
B. N.
,
McIntosh
,
T. K.
,
Gennarelli
,
T. A.
and
Meaney
,
D. F.
(
1997
).
Characterization of diffuse axonal pathology and selective hippocampal damage following inertial brain trauma in the pig
.
J. Neuropathol. Exp. Neurol.
56
,
822
-
834
.
Teasell
,
R.
,
Bayona
,
N.
,
Lippert
,
C.
,
Villamere
,
J.
and
Hellings
,
C.
(
2007
).
Post-traumatic seizure disorder following acquired brain injury
.
Brain Inj.
21
,
201
-
214
.
Tikhonova
,
M. A.
,
Maslov
,
N. A.
,
Bashirzade
,
A. A.
,
Nehoroshev
,
E. V.
,
Babchenko
,
V. Y.
,
Chizhova
,
N. D.
,
Tsibulskaya
,
E. O.
,
Akopyan
,
A. A.
,
Markova
,
E. V.
,
Yang
,
Y.-L.
et al.
(
2022
).
A novel laser-based zebrafish model for studying traumatic brain injury and its molecular targets
.
Pharmaceutics
14
. (
Van Baalen
,
B.
,
Odding
,
E.
,
Maas
,
A. I.
,
Ribbers
,
G. M.
,
Bergen
,
M. P.
and
Stam
,
H. J.
(
2003
).
Traumatic brain injury: classification of initial severity and determination of functional outcome
.
Disabil. Rehabil.
25
,
9
-
18
.
Van Laar
,
V. S.
,
Chen
,
J.
,
Zharikov
,
A. D.
,
Bai
,
Q.
,
Di Maio
,
R.
,
Dukes
,
A. A.
,
Hastings
,
T. G.
,
Watkins
,
S. C.
,
Greenamyre
,
J. T.
,
St Croix
,
C. M.
et al.
(
2020
).
α-Synuclein amplifies cytoplasmic peroxide flux and oxidative stress provoked by mitochondrial inhibitors in CNS dopaminergic neurons in vivo
.
Redox Biol.
37
,
101695
.
Vigil
,
A. F.
,
Belchior
,
H.
,
Bugay
,
V.
,
Bazaldua
,
I. I.
,
Stoja
,
A.
,
Dantas
,
D. C.
,
Chung
,
S. H.
,
Farmer
,
A.
,
Bozdemir
,
E.
,
Holstein
,
D. M.
et al.
(
2023
).
Acute treatment with the M-channel (Kv7, KCNQ) opener retigabine reduces the long- term effects of repetitive blast traumatic brain injuries
.
Neurotherapeutics
20
,
853
-
869
.
Walker
,
A.
,
Chapin
,
B.
,
Abisambra
,
J.
and
Dekosky
,
S. T.
(
2022
).
Association between single moderate to severe traumatic brain injury and long-term tauopathy in humans and preclinical animal models: a systematic narrative review of the literature
.
Acta Neuropathol. Commun.
10
,
13
.
Werner
,
C.
and
Engelhard
,
K.
(
2007
).
Pathophysiology of traumatic brain injury
.
Br. J. Anaesth.
99
,
4
-
9
.
Westerfield
,
M.
(
2000
).
The Zebrafish Book: A Guide for the Laboratory Use of Zebrafish (Danio rerio)
, 4th edn
Eugene
:
University of Oregon press
.
Whiteneck
,
G. G.
,
Cuthbert
,
J. P.
,
Corrigan
,
J. D.
and
Bogner
,
J. A.
(
2016
).
Prevalence of self-reported lifetime history of traumatic brain injury and associated disability: a statewide population-based survey
.
J. Head Trauma Rehabil.
31
,
E55
-
E62
.
Wu
,
B. K.
,
Yuan
,
R. Y.
,
Lien
,
H. W.
,
Hung
,
C.-C.
,
Hwang
,
P.-P.
,
Chen
,
R. P.-Y.
,
Chang
,
C.-C.
,
Liao
,
Y.-F.
and
Huang
,
C.-J.
(
2016a
).
Multiple signaling factors and drugs alleviate neuronal death induced by expression of human and zebrafish tau proteins in vivo
.
J. Biomed. Sci.
23
,
25
.
Wu
,
J. W.
,
Hussaini
,
S. A.
,
Bastille
,
I. M.
,
Rodriguez
,
G. A.
,
Mrejeru
,
A.
,
Rilett
,
K.
,
Sanders
,
D. W.
,
Cook
,
C.
,
Fu
,
H.
,
Boonen
,
R. A. C. M.
et al.
(
2016b
).
Neuronal activity enhances tau propagation and tau pathology in vivo
.
Nat. Neurosci.
19
,
1085
-
1092
.
Xiong
,
Y.
,
Mahmood
,
A.
and
Chopp
,
M.
(
2013
).
Animal models of traumatic brain injury
.
Nat. Rev. Neurosci.
14
,
128
-
142
.
Yablon
,
S. A.
(
1993
).
Posttraumatic seizures
.
Arch. Phys. Med. Rehabil.
74
,
983
-
1001
.
Zhao
,
Z.
,
Loane
,
D. J.
,
Murray
,
M. G.
,
Stoica
,
B. A.
and
Faden
,
A. I.
(
2012
).
Comparing the predictive value of multiple cognitive, affective, and motor tasks after rodent traumatic brain injury
.
J. Neurotrauma
29
,
2475
-
2489
.

Competing interests

The authors declare no competing or financial interests.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

Supplementary information