1. Introduction
Cholera is a profuse diarrheal disease spread via the fecal-oral route through contaminated food and water and is endemic in communities around the world without potable water and sanitary facilities [
1]. Cholera patients pass a characteristic “rice-water stool” and, without rehydration treatment, are at risk of severe dehydration and death [
2]. The World Health Organization estimates an annual 1-4 million cholera cases and upwards of 100,000 deaths internationally, with the most vulnerable demographic being children younger than 5 years old [
3,
4]. The global risk of cholera is characterized as very high, due to the increasing number of outbreaks and lack of vaccine availability; prevalence is expected to worsen with climate change [
5,
6].
Vibrio cholerae, the etiological agent, is a motile Gram-negative bacillus found in aquatic environments [
7,
8].
V. cholerae isolates are classified into over 200 serogroups based on the somatic O-antigen. Some serogroups may cause mild to severe diarrheal symptoms, but only serogroups O1 and O139 are pandemic [
1,
8]. The O1 serogroup is further divided into two biotypes, classical and El Tor. The toxin-coregulated pilus (TCP) and cholera toxin (CT) are the two major virulence factors unique to these biotypes. These allow
V. cholerae to colonize the upper small intestinal epithelium in humans and induce severe diarrheal symptoms, respectively [
1,
9,
10,
11,
12]. Classical strains were the cause of the first six cholera pandemics beginning in 1817. However, beginning in 1961 El Tor became the dominant if not primary cause of the current seventh pandemic [
13]. In comparison to El Tor strains, classical strains cause more severe diarrheal symptoms. However, El Tor strains are thought to be more persistent in colonization, as infection duration lasts longer [
9,
14]. Currently, some circulating El Tor strains, known as atypical or variant El Tor, have acquired some of the classical biotype virulence traits and cause more severe symptoms [
15,
16,
17].
The human gut microbiome contains the majority of commensal bacteria in the body, that are crucial to digestion, nutrient absorption, metabolism, and the immune system [
18,
19]. Firmicutes and Bacteroidetes are normally the dominant phyla of the human gut microbiome, though individual variation occurs due to genetics, age, sex, diet, and lifestyle, among other factors [
20,
21,
22]. The mucus layer of the gut epithelium is essential to maintain a diverse microbiome, as it provides a surface for the resident gut microbiota and contains proteins for the bacteria to metabolize [
23,
24,
25]. During cholera infection,
V. cholerae colonizes the upper small intestine [
11,
26]. Here, CT induces physical efflux of ions and water, and the mucus shedding in tandem drastically decreases gut microbiome diversity. At this stage,
V. cholerae are the most abundant bacteria in the now-altered human gut, and resident gut microbiota fluctuate in response to colonization [
27,
28].
Growing evidence suggests that dysbiosis (a significant change in the microbiome composition) of the gut is implicated in neurological conditions such as anxiety, depression, autism, and Parkinson’s [
29,
30,
31]. Additionally, presence of a microbiome is necessary for normal neurological development [
32,
33]. Because
V. cholerae colonizes and perturbs the gut microbiome, it may be possible that infection affects neurological activity in a similar manner.
Current mammalian animal models used to study
V. cholerae are useful for understanding pathogenesis, but require invasive procedures, use of infant animals, or antibiotic-treated adult animals to enable
V. cholerae to colonize the gut [
34,
35,
36,
37,
38,
39]. Larval and adult zebrafish have been established as natural host models for
V. cholerae, as infection is done via immersion, and zebrafish exhibit diarrheal symptoms similar to humans [
40,
41,
42]. The zebrafish gut microbiome shifts over time among individuals and, without the use of antibiotics or invasive procedures required of other animal models, zebrafish have a mature gut microbiome before and during
V. cholerae colonization [
40,
43,
44,
45]. Adult zebrafish will clear the
V. cholerae infection on their own, and the duration varies based on the biotype strain used: classical strains are cleared within 72 hours post-infection (hpi), and El Tor strains can colonize beyond 144 hpi [
40]. Unpublished data from our lab suggests El Tor persists for up to two weeks. Zebrafish larvae can also be infected via immersion at 5 days post-fertilization (dpf), when the gut opens and uptake of exogenous food begins [
40,
46]. The larval gut microbiome is primarily colonized by bacteria of the phylum Pseudomonadota [
43,
47].
Previous studies have found that microbial colonization is required in germ-free zebrafish larvae for normal swim behavior and neurobehavioral development [
48,
49]. One study found that axenic zebrafish larvae infected with
V. cholerae at 1 dpf exhibited a decrease in hypoactivity at 10 dpf, while conventionally colonized zebrafish larvae infected with
V. cholerae had no change in locomotor activity [
49]. However, this was tested with a dark/light paradigm, whereas a second assay testing for a separate stimulus response would be beneficial to examining gut microbiome crosstalk with different regions of the larval brain and nervous system. An alternating light/dark assay tests the visual system, and a single-tap trial would test locomotor function separate from the visual system. Here, we report
V. cholerae El Tor C6706 colonization of the larval zebrafish gut has no effect on the response to a visual stimulus, but a dampening effect on the response to an acoustic stimulus.
2. Materials and Methods
2.1. Zebrafish Husbandry
Wild-type AB larvae were used, and larvae were fasted at least 12 h before behavioral trials. During experimentation, larvae were placed in autoclaved water from the Aquaneering aquatic housing system filtered by reverse osmosis and maintained at pH 7.0 to 7.5. Tank water was conditioned with Instant Ocean salt (Aquarium Systems, Mentor, OH) to a conductivity of 600 to 700 μS. Fish were kept in a glass-front incubator at 28 °C on a timed 14 hour light, 10 hour dark cycle. All animal protocols were approved by the Wayne State University IACUC.
2.2. Vibrio cholerae Infection
V. cholerae infection of zebrafish has been previously described [
40,
41]. Briefly, El Tor C6706 was incubated with aeration in 30 mL of Luria Broth (LB) agar with 100 µg/mL streptomycin at 37ºC for 12-14 hours. Cells were centrifuged at 8000g, washed and resuspended with 1x phosphate buffer saline (PBS) to a concentration of 10
9 CFU/mL by measuring at OD600 nm. Serial dilutions were plated onto LB agar with 100 µg/mL streptomycin and 100 µg/mL X-galactosidase (X-gal) for confirmation. In 6-well plates (Fisher Scientific, Pittsburgh, PA) the
V. cholerae dilution was pipetted into 5 mL of autoclaved water from the Aquaneering system to an average concentration for both behavioral trials, each repeated in triplicate, of 1.16x10
7 CFU per well. For the uninfected group in separate 6-well plates, the same amount of sterile 1x PBS was used.
Approximately 10-15 larvae at 5 dpf were then placed in each well with the inoculum (uninfected n=88, C6706 n=111). After 6 hpi, all larvae from both groups were removed from wells and separately placed in new wells with 5 mL of new autoclaved water. This rinsing process was repeated three times to wash off the initial inoculum as thoroughly as possible. The larvae were then placed in beakers with 80 mL of new autoclaved water. At 24 hpi larvae were fed larval food (ARTEMAC, Aquafauna Bio-Marine). A water change, rinsing, and feeding were repeated daily following infection. Additionally, 100 µL of water from each experimental beaker was plated directly onto LB Agar with 100 µg/mL streptomycin and 100 µg/mL X-gal daily, to confirm continued V. cholerae infection of C6706 group, and sterility of uninfected group. At 9 dpf, individual larvae were placed in flat bottom 24-well plates (diameter 1.65 inch wells) (Falcon, MA) to acclimate for 24 hours. At 10 dpf, or 5 days post-infection (dpi), all 24-well plates were closed with a lid and sealed in parafilm around the edges, and light-dark or single-tap stimulus trials were performed.
2.3. Behavioral Assays
All behavioral assays were performed as previously described [
50]. Briefly, the larvae in 24-well plates sealed with parafilm were placed in a DanioVision Observation Chamber (Noldus Information Technology, Wageningen, Netherlands) and tracked using the EthoVision XT13 software and a Basler Gen1 Camera (Basler acA1300-60). Camera resolution was set at 1280×960 and the frame rate was set at 25. The DanioVision Temperature Control Unit maintained the larvae at a constant temperature of 28.0±0.5 °C via a steady flow of water to the chamber. All behavioral trials were performed from 1-4 PM. Both behavioral assays were repeated in triplicate.
2.4. Light/Dark Trial
Larvae acclimated in the Noldus in darkness for 12 minutes. After acclimation, recording began and lasted for 24 minutes. A uniform light was emitted at 10,500 lux from below the stage. During all 4 alternating cycles of 3 minutes light and 3 minutes dark, distance moved (cm) and velocity (cm/sec) were tracked (uninfected n=25, C6706 n=60). Average distance and average velocity were calculated in Excel from 30 second time bins and analyzed on GraphPad using Mann-Whitney test with p<0.05 as cutoff for significance.
2.5. Single-Tap Trial
Larvae acclimated in the Noldus in darkness for 12 minutes. After acclimation, recording began and lasted for 2 minutes. Larvae were exposed to 1 minute of no stimulation, followed by a single-tap at the highest intensity setting, and ending with 1 minute of no stimulation. Distance moved (cm) was tracked (uninfected n=63, C6706 n=51). Average distance was calculated in Excel from 1 second time bins and analyzed on GraphPad using Mann-Whitney test with p<0.05 as cutoff for significance.
2.6. Zebrafish Euthanization and Homogenization
All larvae were euthanized in a lethal dose of MS-222 (300 mg/L Tris buffer, pH 7.0) for 30 minutes and homogenized with a pellet pestle (Fisher Scientific, Pittsburgh, PA). Before infection at 5 dpf, 15 larvae were homogenized together in 100 µL of 1x PBS. At 10 dpf, 5 uninfected larvae and 5 C6706-infected larvae were homogenized together in 100 µL of 1x PBS.
2.7. DNA Isolation and Sequencing
DNA from 100 µL of larval zebrafish’s combined homogenate was isolated using the DNEasy Powersoil Pro kit (Qiagen, Germantown, MD) per the manufacturer’s instructions. Two extractions using only sterile 1x PBS and no gut samples were also processed as kit controls. The V4 region of the 16s rRNA gene was targeted and amplified using 515F and 806R primers. These primers have yielded successful results in previous zebrafish microbiome studies, and the V4 region has demonstrated high reproducibility [
45,
51]. Samples were submitted to Michigan State University for Illumina MiSeq Sequencing using previously established methods [
52,
53].
All raw data were processed with R package ‘dada2tools’, available at
https://github.com/jp589/dada2tools, to efficiently correct Illumina amplicon errors without generating operational taxonomic units. Instead, amplicon sequence variants (ASVs) are derived based on 100% sequence similarity. Modifications to an online MiSeq protocol (
https://benjjneb.github.cio/dada2/tutorial.html) included allowing truncation lengths of 245 bp and 210 bp and a maximum number of errors of 2 bp for forward reads and 7 bp for reverse reads. Sequences were then classified into taxa using the silva_nr99_v138.1_train_set database with a minimum bootstrap value of 80% [
54]. classified as mitochondrial, chloroplast, or not classified at phylum level were removed. Based on the bacterial profiles of two blank DNA extraction kit samples, 1 ASV, an unclassified
Corynebacterium, was further removed from the dataset (
https://github.com/jp589/dada2tools).
65 ASVs were detected in the larval gut microbiome. The samples of the gut microbiome of larvae at time zero yielded far fewer 16S rRNA gene sequencing reads (6772) than those of either uninfected control (252,894) or C6706 (181,224) larvae at five days post-infection. Thus, relative abundance data were used for descriptive comparisons among these three groups of samples.
4. Discussion
The larval zebrafish brain has approximately 100,000 neurons [
55]. Identifying a dampened motor response in C6706-colonized larvae to an acoustic stimulus, but not to visual stimuli, could indicate that the gut microbiome has crosstalk with different regions and neurons of the larval brain. An open-source Zbrain atlas of the larval zebrafish brain has been established, opening the door for neuronal activity mapping [
55]. Acoustic stimuli were found to activate the ears and lateral line, which directly connects to the octavolateralis nucleus (ON), as well as particularly strong neuron activations in the torus semicircularis, thalamus, cerebellum, and remaining hindbrain [
55,
56,
57,
58]. Visual stimulus in the form of a 10 second light flash activated the retinal projections and diencephalic areas of the larval brain [
55]. In relation to our data, this could indicate a more direct pathway between the gut microbiome and the larval brain regions associated with a response to an acoustic stimulus. More specifically, a
V. cholerae El Tor C6706-colonized gut microbiome is either indirectly or directly weakening crosstalk due to the presence of C6706 or absence of resident gut microbiota. Larval locomotion studies are translational in that activation of reticulospinal neurons in the brain stem and the vestibulospinal tract are conserved in vertebrates [
59]. Structures of the larval zebrafish brain have evolutionarily conserved homologous functions to other vertebrates’ [
60,
61].
Our study aimed to identify any behavioral response from two separate stimuli. Screening for other larval behaviors may yield results highlighting other regions of the brain affected by a gut microbiome shift. This includes multi-tap assays to measure habituation, circadian rhythm assays, and prey-capture to measure decision-making [
59,
62,
63,
64,
65]. Whole-brain imaging would be another direction towards larval neurology during
V. cholerae colonization. Behavioral studies are also possible for adult zebrafish colonized with
V. cholerae. While adult zebrafish exhibit much more complex, continuous behaviors, on a broader scale, translational relevance is equally complex [
66]. For example, adult zebrafish social phenotypes were parallel to social interaction cues observed in humans, such as head direction and physical distance [
67]. A link with
V. cholerae colonization to behavior in larval and adult zebrafish could translate to environmental observations, where
V. cholerae and zebrafish potentially interact.
The Comamonadaceae family,
Rheinheimera coerulea, and Flectobacillus species have been isolated from freshwater environments [
68,
69,
70].
Pseudomonas has been established as a part of the zebrafish core gut microbiome [
47]. Stephens et al. has shown that unclassified Comamonadaceae took up 97.5% of all larval intestines [
43].
Brevundimonas kwangchunensis is the one ASV that has not been thoroughly described, although the
Brevundimonas genus has been found in soil and water samples [
71,
72]. It is interesting that
Vibrio was not one of the most abundant ASVs from the C6706-infected group, though the few larval homogenates plated from this time did yield X-gal blue CFUs indicative of
V. cholerae (not shown). It could be possible that the specific group of larval homogenates submitted for sequencing were poorly infected and colonized with
V. cholerae in comparison to other larvae. The gut microbiome abundances and number of ASVs significantly changing from 5 dpf to 10 dpf prove that gut microbiome diversity can quickly shift at this early developmental stage [
43]. Additionally, perhaps plating homogenates on LB without streptomycin would provide confirmation of other, more abundant bacterial species’ growth. Since these behavioral analyses were the first to consider
V. cholerae, future trials can include more larvae for sampling to provide a more comprehensive model of the gut microbiome.
In terms of V. cholerae colonization, while V. cholerae El Tor C6706 was tested, it is possible that classical and non-O1 biotype strains induce different behavioral responses. Future studies could determine if interactions between the larval gut and brain play a role specific to the more persistent colonization of V. cholerae El Tor strains.
Author Contributions
I.C., R.T. and J.H.W.; methodology, X.L. and S.B.; software, X.L. and S.B.; validation, I.C..; formal analysis, I.C.; investigation, I.C.; resources, X.L. and R.T.; data curation, I.C., K.R.T. and J.P..; writing—original draft preparation, I.C.; writing—review and editing, I.C., J.P. and J.H.W.; visualization, I.C.; supervision, J.H.W.; project administration, R.T. and J.H.W.; funding acquisition, J.H.W. All authors have read and agreed to the published version of the manuscript.