Effect of nanobubbles (oxygen, ozone) on the Pacific white shrimp (Penaeus vannamei), Vibrio parahaemolyticus and water quality under lab conditions

Nguyen Huu Nghia1,*, Nguyen Thi Nguyen1, Phan Trong Binh1, Le Thi May1, Tong Tran Huy1, Pham Thai Giang1, Sophie St-Hilaire2, Phan Thi Van1
Author Information & Copyright
1Research Institute for Aquaculture No.1, Dinh Bang, Tu Son, Bac Ninh 16000, Vietnam
2Department of Infectious Diseases and Public Health, Jockey Club College of Veterinary Medicine and Life Sciences, City University of Hong Kong, Kowloon Tong 518057, Hong Kong
*Corresponding author: Nguyen Huu Nghia, Research Institute for Aquaculture No.1, Dinh Bang, Tu Son, Bac Ninh 16000, Vietnam, Tel: +84-2438785751,

Copyright © 2022 The Korean Society of Fisheries and Aquatic Science. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: Nov 16, 2021; Revised: Mar 07, 2022; Accepted: Jul 25, 2022

Published Online: Aug 31, 2022


This study assessed the effects of oxygen and ozone nanobubbles on gill morphology, weight gain, and mortality of Pacific white shrimp (Penaeus vannamei), as well as the level of Vibrio parahaemolyticus and water quality of shrimp culture tanks under lab conditions. Two experiments were carried out with oxygen macrobubble, ozone macrobubble, oxygen nanobubble, ozone nanobubble, and control treatments (air-stone macrobubble). Experiments were done in triplicate in 100 L tanks with 15‰ saline water, and 20 shrimp per tank. Tanks in Experiment 1 were not inoculated with bacteria; tanks in Experiment 2 were inoculated with V. parahaemolyticus at a concentration of 106 CFU/mL. The results revealed that short treatments with ozone nanobubbles had minimal impact on shrimp gills, mortality, and growth rates, reduced V. parahaemolyticus concentration in water compared to the other groups, and improve water quality. These laboratory results indicate that ozone nanobubble treatment may be useful for controlling V. parahaemolyticus. More work is needed to find the best protocol to apply the technology on a commercial scale.

Keywords: Ozone nanobubble; Penaeus vannamei; Gill morphology; Vibrio parahaemolyticus; Water quality


Aquaculture has experienced rapid growth in the past two decades. World aquaculture production increased from 35 million tons in 1997 to 114.5 million tons in 2018 (FAO, 2020). This pattern of growth has also occurred in Vietnam but faster with a nearly 100-fold increase in production from 492 thousand tons in 1997 to 4.8 million tons in 2021, with an export value of 8.89 billion USD. Production of farmed shrimp in 2021 was estimated at 930 tons of which 632 tons was whiteleg shrimp (Penaeus vannamei), 268 tons was tiger shrimp (Penaeus monodon), and the remainder was other shrimp species (TCTS, 2021). Due to the rapid growth of the aquaculture industry in Vietnam there have been issues with environmental degradation and disease outbreaks. One of the diseases causing severe losses to the shrimp industry is acute hepatopancreatic necrosis disease (AHPND), caused by the gram-negative bacterium Vibrio parahaemolyticus. It is estimated that losses from AHPND in the Mekong Delta region between 2011 and 2015 were approximately 26 million USD for whiteleg shrimp and 11 million USD for black tiger shrimp (FAO, 2018). There is currently no effective control measure for this disease as antibiotics are not effective. The industry therefore needs to find alternative measures to prevent and control this pathogen.

Nanobubble technology produces small gas bubbles with a diameter of less than 200 nm (Agrawal et al., 2011; Chaplin, 2021) which can exist in the water column for weeks (Azevedo et al., 2016; Parmar & Majumder, 2013). This phenomenon is associated with the physical characteristic of these small bubbles (Weijs & Lohse, 2013) and high zeta potential value which keeps them from aggregating into larger bubbles (Kirby, 2010; Sjogreen et al., 2018). Nanobubble technology has a lot of promising effective application in aquaculture because of its ability to increase and maintain oxygen content in water and control bacteria, a possible measure to reduce the use of antibiotics.

Oxygen nanobubbles optimize the dissolved oxygen in aquatic systems (Galang et al., 2019). High levels of oxygen may speed up metabolism, promote higher food intake, and increase growth of farmed shrimp (Ebina et al., 2013). Ozone nanobubble has recently been shown to effectively disinfect water in laboratory experiments and small mesocosm experiments (Imaizumi et al., 2018; Jhunkeaw et al., 2021; Nghia et al., 2021; Seki et al., 2017). The fundamental disinfection properties of ozone nanobubbles are that they attract negatively charged particles and create oxygen free radicals which disrupt the permeability of cell membranes (Gurung et al., 2016; Ikehata & Li, 2018; Temesgen et al., 2017). In water, ozone is converted to oxygen and therefore it increases dissolved oxygen levels while disinfecting water (Nano Bubble Technologies, 2020). In marine recirculating aquaculture system moderate ozonation seems to promote biological nitrification. This is likely mediated by the increase in oxygen promoting nitrogen converting bacterial populations. The limiting factor for the use of ozone in salt water is the toxicity to the fish (Schroeder et al., 2015). Effective ozone treatment depends on concentration, exposure time, pathogen load and organic matter content. The higher the ozone concentration, the higher the risk to famers and farmed species (Gonçalves & Gagnon, 2011).

Although nanobubble can increase and maintain oxygen in water, gas overdose can cause hyperoxia. The effect of air supersaturation on aquatic animals has documented in many studies. Hyperoxia caused a 67% increase in total interlamellar cell mass (Tzaneva et al., 2011). Boyd & Fast (1992) reported that fish and shrimp mortality increased when oxygen exceeded 20 mg/L and or the saturation rate exceeded 300%. This condition can result in gas bubble disease (gas bubble trauma), which has been described in a wide variety of fishes and invertebrates. Weitkamp & Katz (1980) reported that bubbles were observed in the gills of dying fish as well as between fin rays and under scales when dissolved oxygen levels are over 300% saturation and numerous marine fish died when dissolved oxygen concentrations were above 250% of saturation in Galveston Bay. Anecdotal reports from shrimp industries in Asia suggests oxygen nanobubble technology is effective at increasing dissolved oxygen in pond shrimp culture. However, despite some knowledge of the physical and chemical properties of nanobubbles, much remains unknown about their effect on aquatic animals (Atkinson et al., 2019).

In this study, we investigated the effects of nanobubbles of oxygen and ozone gases on gill morphology, weight gain and mortality of Pacific white shrimp (P. vannamei) as well as V. parahaemolyticus load and water quality of shrimp culture tanks under laboratory conditions.

Materials and Methods

Laboratory set up

Experiments were carried out from May to September 2020 in a laboratory at the Centre for Environment and Disease Monitoring in Aquaculture (CEDMA), Research Institute for Aquaculture No. 1, Bac Ninh, Vietnam (RIA1). Trials were conducted in 140 L cylindrical composite tanks with 100 L of 15‰ saline water, which is optimal for V. parahaemolyticus (WHO & FAO, 2011). Saline water was prepared by adding 1.5 kg Blue Treasures Sea salt (Qingdao Sea-Salt Aquarium Technology, Qingdao, China) to 98.5 L fresh water.

A nanobubble generator model aQua+75MO (AquaPro Solutions, Singapore, Singapore) with a 1 HP power (0.75 kW) motor and a water flow output rate of 1,000 L/h was used to generate nanobubbles. The machine produced nanobubbles with an average size of 168.9 ± 73.8 nm, and a concentration of 1.04 × 109 ± 2.6 × 108 particles/mL (AquaPro Solutions, 2019). To create oxygen nanobubbles, we used an oxygen concentrator (model Yuwell 7F-10, Yuwell, Danyang, China; 10 L/min max flowrate, ≥ 95.5% oxygen concentration, flow rate at 4 L/min) provided oxygen for nanobubble generator. To create ozone nanobubbles we used the above mentioned oxygen concentrator provided oxygen to ozone generator (model OM-Z10, OzoneMaxx, Hanoi, Vietnam; output flow rate of 15 L/min, produce 10 grams of ozone per hour) to produce ozone, and using the ozone for nanobubble generator as described by Nghia et al. (2021).

Experimental design

Two separate experiments were designed to investigate the effects of oxygen nanobubbles and ozone nanobubbles on shrimp’s gill morphology, weight gain, and mortality. These experiments were also designed to study the effects of the gases nanobubbles to V. parahaemolyticus concentration in water, and physicochemical parameters of water in the experiment tanks.

Experiment 1 (without V. parahaemolyticus) consisted of five treatments: oxygen macrobubble (O2), ozone macrobubble (O3), oxygen nanobubbles (O2-NB), ozone nanobubbles (O3-NB), and negative control which was given air macrobubble (NCTRL). Each group was replicated three times. The O2 and O2-NB tanks were treated for 4 minutes once per day, while O3 and O3-NB tanks were treated for 2 minutes once per day. Optimal treatment times were set based on previous research results (Nghia et al., 2021) and a pilot trial prior to this experiment. There were no bacteria added to the systems in Experiment 1 (only in Experiment 2), but we monitored bacterial count in the water to determine if different gases could prevent bacterial proliferation and improve water quality. The treatments were continued for 21 days.

Experiment 2 (with V. parahaemolyticus) consisted of six treatments: O2, O3, O2-NB, O3-NB, NCTRL no V. parahaemolyticus exposure, only air macrobubble, and positive control (PCTRL) exposed to V. parahaemolyticus and air macrobubble. Each treatment was replicated three times. The O2, O2-NB, O3-NB tanks were treated for 10 minutes once per day, while the O3 tanks were treated for 2 minutes once per day. Except for the NCTRL, all tanks in Experiment 2 were inoculated with V. parahaemolyticus at a concentration 106 CFU/mL. The reason we set 10 minutes treatments for all treatments except the O3 instead of the shorter exposure used in Experiment 1 because we needed more time to disinfect against the bacteria. O3 macrobubble treatment were set 2 minutes because macrobubbles has more effect to shrimp than nanobubbles, and a longer duration would have resulted in mortality of the shrimp (data not shown). Every 48 hours, 50 mL of Nutrient Broth culture medium was added to each tank to stimulate and maintain the optimum condition for bacteria growth. The treatments were continued for 21 days.

All tanks in both experiments stocked with 20 shrimp (40–45 days post-larva) from a local farm and supplied with an aquarium air-stone operated continuously throughout the trials. All animals were acclimated in a composite tank (100 L) at 15‰ salinity for 1 week prior to the study. Shrimp were fed 3 times a day for a total of 3% of their body weight and examined daily. Waste at the bottom of the tanks were collected daily, water was added to the tanks to compensate for water lost during cleaning. Moribund animals (i.e., shrimps that were weak, had no reaction, had no food in their gut, and sank to the bottom of the tank) were collected for sampling immediately. Water quality was monitored daily for 21 days.

Effects of nanobubbles on the Pacific white shrimp (Penaeus vannamei)
Gill morphology

Shrimp gill samples were prepared for histology following Lightner (1996). Sections (5 µm thickness) were stained with hematoxylin-Mayer and eosin. Preparations were observed under microscope (Nikon E200 LED, China) at 400× magnification. The degree of gill pathology was graded into 5 main levels: normal filament cells (L0), filament cells adhesion (L1), cell’s nucleus adhesion (L2), gill arch atrophy (L3), and gill arch caseation (L4). Each major level was further divided into 4 sub-levels (score) (e.g., L0.1 to L0.4) to indicate the severity of the gill damage. The number of analyzed samples (shrimp) from each treatment depended on the number of moribund caught shrimp and remaining live shrimp at the end of experiment. To obtain an overall score for each gill sample we multiplied the main level (e.g., 1–5) by the sub-level (e.g., 1–4) and these numbers were summed up for each shrimp, the calculation is expressed by the following formula:

Total score = y = 1 y = 4 x = 0 x = 4 ( x + 1 ) × y i .

Where x is the main level and y is sub-level.

These scores were then recoded to normal (total < 2) and then subsequently categorized into four groups based on quartiles of the scored data, i.e., light, medium, severe, and very severe pathology with approximately equal sample size. In order to have the same scale in both trials, we based the selection of cut-off points on trial 1 for both trials.

Weight gain

Before starting the experiment, shrimp were weighed individually using an analytical balance (PA214, Ohaus, Parsippany, NJ, USA) after removing water from the shrimp with a soft tissue paper. At the end of the experiment, live shrimp in the experimental tanks were weighed as described above. We calculated the average weight of the shrimp in each tank at the start and at the end of the study. We compared the change in average weight across tanks in different treatment groups.


The number of dead shrimps was recorded according to the experimental tank. Shrimp activity and status were checked two times a day during the experiments. Shrimp that sank to the bottom of the tank were collected immediately, euthanized, and fixed with Davison solution for histopathology analysis. A bacterial culture of the hepatopancreas on thiosulfate citrate bile sucrose agar (TCBS) media was done to test for V. parahaemolyticus to determine if this was the cause of death prior to fixation.

Effects of nanobubbles on Vibrio parahaemolyticus
Vibrio parahaemolyticus preparation and infection

The AHPND strain of V. parahaemolyticus was collected from the bacteria storage bank in CEDMA, RIA1 and was prepared as follows: 1) V. parahaemolyticus bacteria were thawed from ultra-low temperature freezer –80°C and defrosted at room temperature; 2) the bacteria were cultured in selective medium (TCBS) in an incubator at 29°C. After 18–24 hours, single colonies were put into sterilized Erlenmeyer flask containing nutrient broth with 2% NaCl medium and shaken at 29°C for 18 hours; and 3) the bacterial concentration was optically estimated at OD600 nm with an Eppendorf Bio Spectrometer (Eppendorf, Hamburg, Germany) and enumerated by dilution and plating. Bacteria was introduced into experiment tank at the concentration 106 colony forming units per mL of water (CFU/mL).

Vibrio parahaemolyticus count in the water

Bacterial counts were measured at 0, 6, 12, 24 hours at the beginning of the experiments and thereafter once a day. The tank’s water was stirred before sampling. Bacterial concentration in water was quantified following Buller (2004). Collected water samples from experimental tanks were diluted 10-fold with 2% sodium chloride solution. A volume of 100 µL of dilution solution was inoculated on a TCBS agar plate and a glass was used to evenly spread the bacteria on the agar surface. Cultured TCBS agar plates were incubated in the incubator at 29.0°C. After 24 hours, colonies were counted, and the bacterial concentration was calculated using the following formula:

X = A V × K .

Where: X is bacterial concentration in 1 mL sample (CFU/mL), A is the number of colonies growing on agar plates, V is the volume of water introduced into the culture (for example, if 100 µL is added, V = 0,1), K is the dilution factor (for example, at a concentration of 10-1, K = 10). The minimal detection limit in our lab was 14 CFU/mL.

Vibrio parahaemolyticus detection in shrimp

The hepatopancreas of each shrimp was dissected and DNA was extracted as described in MOSTE (2019). AP3 primer pairs with sequences: 5’-ATGAGTAACAATATAAAAACATGAAAC-3’ and 3’ GTGGTAATAGATTGTACAGAA-5’ was used to test for the presence of V. parahaemolyticus by PCR using gel electrophoreses to visualize results.

Water quality parameters assessment

Water quality parameters including temperature, pH, oxygen (DO), alkalinity (Alk), oxidation reduction potential (ORP), nitrite (NO2), hydrogen sulfide (H2S), ammonia (NH3), ammonium (NH4+) were measured at 0, 6, 12, 24 hours, and thereafter once a day. Tank water was stirred before sampling. Temperature, DO, pH, ORP were measured by a Pro1020 Dissolved Oxygen and pH or ORP Instrument (YSI, Yellow Springs, OH, USA). Alk, NO2, NO3, H2S, NH3, NH4 were analyzed by the SMEWW 2302B: 2011, SMEWW 4500-NO2 B: 2011, 4500-NO3 B: 2011, SMEWW 4500-S2- D: 2011, SMEWW 4500-NH3 F: 2011 method, respectively (Baird & Bridgewater, 2017).

Statistical analysis

Gill scores were analyzed using ordinal logistic regression to compare categories (normal, light, medium, severe, and very severe changes) among treatments and to estimate probabilities for occurrence of these categories. Change in weight of shrimp was compared between treatments. An interaction between the two factors (treatment and time) was used in the linear regression. Total mortality at the end of the experiments was compared among treatments using a logistic regression model. In Experiment 2, survival curves were compared among treatments using the Cox proportional hazards model. We also depicted mortality curves using Kaplan-Meier failure function.

Log transformed bacterial counts (x) Log10(x + 1), and the water quality parameters were compared between treatments and over time using linear regression. We used an interaction term between time and treatment in the model to control for the repeated measures. For bacteria counts, the negative control was excluded from the analysis. Model assumptions were assessed using residual plots.

Analyses were carried out using XLSTAT statistical and data analysis solution (XLSTAT, 2020). The p-values < 0.05 were considered significant.


Effects of nanobubbles on the Pacific white shrimp (Penaeus vannamei)
Shrimp gill morphology

In Experiment 1, all shrimp from the control group scored as normal, i.e., no gill change observed. The highest gill damage level L4.2 occurred in 18 shrimp in the O3 treatment, 7 shrimp in the O2 treatment, and 2 shrimp in O2-NB treatment. No L4.2 gill scores were observed in the O3-NB treatment (Table 1). The relative occurrence of the various composite gill change scores did not differ significantly between O2 and O2-NB treatment and had considerably fewer animals with severe gill damage than the O3 treatment (p < 0.001). Differences among the O3 treatment and the two O2 treatments were not significant (Fig. 1).

Table 1. Summary occurrences of the various pathology indicators by treatment of the E1 and E2, negative control is excluded
Experiment E1 E2 E1 E2 E1 E2 E1 E2 E1 E2
Treatment NCTRL PCTRL O2 O2 O2-NB O2-NB O3 O3 O3-NB O3-NB
No. of shrimp NA 31 58 37 54 22 49 31 53 38
L0.1 NA 0 0 0 0 0 0 0 0 0
L0.2 NA 6 0 0 0 3 4 2 18 13
L0.3 NA 8 12 4 9 5 20 12 14 19
L0.4 NA 6 22 12 20 10 24 7 22 3
L1.1 NA 17 31 27 27 10 30 25 31 35
L1.2 NA 0 20 0 17 2 5 4 5 3
L1.3 NA 0 1 0 4 0 11 0 5 0
L1.4 NA 0 0 0 3 0 0 0 0 0
L2.1 NA 10 44 8 26 11 35 13 32 19
L2.2 NA 2 8 2 18 6 2 0 1 0
L2.3 NA 0 1 0 1 0 0 0 1 0
L2.4 NA 0 0 0 0 0 0 0 0 0
L3.1 NA 5 10 4 7 6 0 5 2 1
L3.2 NA 1 2 3 3 0 0 0 0 0
L3.3 NA 0 0 0 0 0 0 0 0 0
L3.4 NA 0 0 0 0 0 0 0 0 0
L4.1 NA 11 18 7 16 8 0 11 3 6
L4.2 NA 2 7 8 2 1 18 2 0 1
L4.3 NA 8 0 14 0 2 0 7 0 2
L4.4 NA 1 0 3 0 0 0 0 0 0

E1, Experiment 1; E2, Experiment 2; NCTRL, negative control which was given air macrobubble; PCTRL, positive control which was given air macrobubble; O2, oxygen macrobubble; O2-NB, oxygen nanobubble; O3, ozone macrobubble; O3-NB, ozone nanobubble; NA, not available.

Download Excel Table
Fig. 1. Predicted probabilities for occurrences of the various categories of composite scores of gill change of the treatments: NCTRL, PCTRL, O2, O2-NB, O3, O3-NB of E1 and E2. E1, Experiment 1; E2, Experiment 2; NCTRL, negative control which was given air macrobubble; PCTRL, positive control which was given air macrobubble; O2, oxygen macrobubble; O2-NB, oxygen nanobubble; O3, ozone macrobubble; O3-NB, ozone nanobubble.
Download Original Figure

In Experiment 2, all treatments except the negative control had gill damage. There was a few shrimp (i.e., one in the positive control and four in the O2 treatment) had very high gill damage scores of (Table 1). All shrimp from the negative control and a few shrimp from the other treatments had normal scores. The two O2 treatments differed significantly (p < 0.001) from the O2-NB which had fewer animals with lower damage scores compared to the O2 treatment. The two O3 treatments also differed significantly (p < 0.01) in which O3-NB groups had only a few shrimp with very high gill damage scores. Relative occurrences of different gill scores were not different between the positive control, O2-NB, and O3 treatments (Fig. 1).

Shrimp weight gain

In Experiment 1, the size of shrimp did not vary significantly among treatments at the start of the experiment. The average weight gain in the control group at the end of the study was 0.7871 ± 0.044 g. The weight gain in the O3 treatment group was significantly lower than in the other four groups (p < 0.05). Differences in weight gain among the other four treatments did not differ significantly.

In Experiment 2, there were significant differences among the groups at the start of the experiment; the ozone group had lower mean weight than the other groups (p = 0.007). The average weight gain in the negative control group was 0.5977 ± 0.06 g and differences among treatments were not significant (Fig. 2).

Fig. 2. Weight gain of shrimp by treatments: NCTRL, PCTRL, O2, O2-NB, O3, O3-NB of E1 and E2. E1, Experiment 1; E2, Experiment 2; NCTRL, negative control which was given air macrobubble; PCTRL, positive control which was given air macrobubble; O2, oxygen macrobubble; O2-NB, oxygen nanobubble; O3, ozone macrobubble; O3-NB, ozone nanobubble.
Download Original Figure
Shrimp mortality

In Experiment 1, the odds of shrimp dying in the O3 treatment group was 4.67 times higher than in the control treatment (p < 0.05), while differences among the other treatments were not significant. There was not a significant difference in the odds of shrimp dying in the O2-NB group or in the O3-NB groups compared to shrimp in the control group (odds ratio [OR] = 1.27, p = 0.73; OR=1.85, p = 0.35).

In Experiment 2, the mortality rate in the negative control was 0.05 (5% of shrimps had died). All treatments had a higher mortality at the end of the study than the negative control. The odds of dying in the positive control was 28.5 times that of the negative control (p < 0.001). Similar mortality rates 62.4 (p < 0.001), 32.8 (p < 0.001) and 19 (p < 0.001) were detected in the O2 treatment, O2-NB treatment and the O3 treatment, respectively. The odds of dying in the O3-NB treatment was only 3.35 times that of the negative control group and this was not statistically significant (p = 0.081). Differences among PCTRL, O2 and O2-NB treatments were not statistically significant. Mortality was higher (p < 0.001) in the O3 treatment than in the O3-NB treatment (Fig. 3). Mortality curves differed significantly among treatments (p < 0.001). Mortality occurred during the first few days of the experiment, which was low after then.

Fig. 3. Mortalities rate of different treatments at the end of experiment by treatments: NCTRL, PCTRL, O2, O2-NB, O3, O3-NB of E1 and E2. E1, Experiment 1; E2, Experiment 2; NCTRL, negative control which was given air macrobubble; PCTRL, positive control which was given air macrobubble; O2, oxygen macrobubble; O2-NB, oxygen nanobubble; O3, ozone macrobubble; O3-NB, ozone nanobubble.
Download Original Figure
Effects of nanobubbles on Vibrio parahaemolyticus
Vibrio parahaemolyticus concentration in water

Changes in bacteria counts [as Log10(count + 1)] over time was modelled. Bacteria counts declined in a curvilinear manner with Log10(time in hours) and therefore we used polynomial regression to analyze the differences in patterns. Bacteria counts declined quickly in all tanks but particularly so in the two O3 treatment groups (Fig. 4). No bacteria were observed in the negative control. Bacterial counts at important milestones such as 6, 120, 240, and 480 hours of the experiment were compared using marginal means. By 6 hours post treatment bacteria counts were markedly lower in the two ozone treatments than in the positive control and the two oxygen treatments (p < 0.001). The same pattern was also observed at the other time points although bacteria density had gone down in all treatments.

Fig. 4. Scatterplot of Vibrio parahaemolyticus counts (points) and predicted counts (lines) based on polynomial regression. The dotted line is the initial concentration, and the dashed line indicate the detection limit. PCTRL, O2, O2-NB, O3, O3-NB of Experiment 1 and Experiment 2. PCTRL, positive control which was given air macrobubble; O2, oxygen macrobubble; O2-NB, oxygen nanobubble; O3, ozone macrobubble; O3-NB, ozone nanobubble.
Download Original Figure
Vibrio parahaemolyticus detection in shrimp

Additional, PCR test of dead and live shrimp suggest V. parahaemolyticus was found in 29 of 38 shrimp (76.3%) in the positive control. The odds of death associated with V. parahaemolyticus (ratio of positive and negative probability) was 0.5862 (n = 26), 0.6466 (n = 37), 0.3103 (n = 22) and 0.1194 (n = 18) in the O2, O2-NB, O3 and O3-NB treatments, respectively.

Effects of nanobubbles on water quality

Temperature and pH in both experiments changed in a similar pattern for all treatments. Temperature in the NB treatment groups was significantly higher than the other treatments (p < 0.001). The pH level in the O2 treatments was significantly higher than in the O3 treatments (p < 0.001). In both experiments, ORP values of the O3 treatments were significantly higher than in other treatment groups (p < 0.001); and DO content of the O3 and O2 treatment tanks were significantly higher than the NCTRL and PCTRL treatment tanks (p < 0.001). O3 was measured only in O3 and O3-NB treatments. O3 values of O3-NB treatment were significantly lower than O3 treatment tanks (p < 0.001) in Experiment 1 (Table 2).

Table 2. Summary of physico-chemical factors i.e., mean (min–max) of two experiments
Treat NCTRL PCTRL O2 O2-NB O3 O3-NB p 1 p 2 p 3
Experiment 1
 Temp (°C) 28.8 (26.6–30.0) 28.6 (26.9–29.8) 29.0 (27.3–30.4) 28.7 (27.1–29.8) 28.8 (27.2–29.9) < 0.001 0.0515 < 0.001
 pH 8.32 (7.99–8.59) 8.25 (7.91–8.71) 8.23 (7.87–8.48) 8.28 (8.01–8.59) 8.26 (7.95–8.47) < 0.001 < 0.01 < 0.001
 ORP (mg/L) 227 (158–576) 254 (220–348) 307 (258–391) 560 (258–845) 449 (257–765) < 0.01 1.0000 < 0.001
 DO (mg/L) 5.3 (4.3–6.3) 15.0 (5.2–19.4) 16.8 (4.1–22.6) 12.2 (5.3–18.6) 13.8 (5.0–18.5) 1.0000 0.9898 < 0.001
 O3 (mg/L) 0.234 (0–0.40) 0.027 (0–0.10) < 0.001 < 0.001 < 0.01
Experiment 2
 Temp (°C) 28.8 (27.9–30.1) 28.9 (28.1–29.8) 28.9 (28.1–29.8) 29.8 (28.1–30.8) 29.0 (28.1–30.7) 29.9 (28.2–30.8) < 0.001 < 0.01 < 0.001
 pH 8.07 (7.89–8.17) 8.06 (7.91–8.19) 8.09 (7.99–8.21) 8.05 (7.88–8.15) 8.07 (7.94–8.19) 8.04 (7.86–8.15) < 0.001 0.5364 < 0.001
 ORP (mV) 237 (170–283) 234 (176–289) 258 (202–326) 270 (231–299) 269 (222–343) 277 (29–319) < 0.001 0.9534 < 0.001
 DO (mg/L) 6.2 (4.5–11.0) 6.0 (4.6–9.9) 23.2 (5.9–37.7) 20.9 (5.2–36.7) 11.9 (5.9–20.7) 21.3 (5.0–31.4) < 0.001 0.9672 < 0.001
 O3 (mg/L) 0.05 (0–0.39) 0.07 (0.00–0.43) < 0.001 0.9058 < 0.01

p1 is p-value of comparison of parameter values variation in the control treatment over time, p2 is p-value of comparison of treatments at time = 0, p3 is p-value of comparison over time between treatments.

NCTRL, negative control which was given air macrobubble; PCTRL, positive control which was given air macrobubble; O2, oxygen macrobubble; O2-NB, oxygen nanobubble; O3, ozone macrobubble; O3-NB, ozone nanobubble; ORP, oxidation reduction potential; DO, oxygen.

Download Excel Table


Our findings suggested that treatments with nanobubbles could result in mild to moderate gill pathology. In our first experiment, the majority of shrimp with gill changes were classified with mild changes, and shrimp were active and had normal growth. Our findings were consistent with other studies reporting supersaturation with oxygen adversely affecting the gill structure of shrimp (Boyd & Fast, 1992; Weitkamp & Katz, 1980) and fish (Tzaneva et al., 2011). A recent study by Jhunkeaw et al. (2021) also did not find mortality and or cell damage in tilapia (Oreochromis niloticus) 48 hours post exposure to a 10 minutes treatment with O3-NB. However, perhaps more frequent shorter exposures to O3-NB would eventually irritate the gills of fish.

The morphological changes observed in our study were more severe in Experiment 2 than in Experiment 1 likely because shrimp were exposed to nanobubbles for a longer period, and they were exposed to bacteria. Experiment 2 suggests the mild damage to gills associated with exposure to O3-NB do not outweigh the benefits that this treatment has on disinfecting water during an outbreak of Vibriosis.

Although it is possible that poor water quality during our treatments could account for some of the gill pathology as nitrogen compounds can result in epithelial oedema, haemolymph infiltration, multifocal melanization and necrosis in shrimp (Fregoso et al., 2017), the values for NH3, NO2, and NO3 in our tanks were within normal range for shrimp.

The result of both experiments suggested that repeated O3 treatments (macrobubble) reduced growth rate of shrimp during our 21 days experiment. However, using O3 in nanobubble form did not affect the growth rate of shrimp. This observation paralleled the gill pathology and the mortality results. In effect, when O3 is delivered using nanobubble it appears to be less toxic and irritating to shrimp than when it was delivered with an aerator (macrobubble). Several water quality parameters differences were observed during the treatment with O3-NB and O2-NB compared to other treatments and these may have contributed to the weight gain differences. For example, the use of nanobubble in tanks increases water temperature by up to 2°C which could have accounted for the weight gain result in our study, however, the changes in weight were not large and were not significant except for the reduction in growth in the O3 macrobubble group compared to the other treatments. The use of O3 treatment without nanobubble appears to have a negative impact on growth of shrimp, so care should be taken when applying this gas using an aerator on commercial farms.

In the absence of bacteria (Experiment 1) the O3 macrobubble treatment had the highest mortality. This finding was supported by Gonçalves & Gagnon (2011) who found that higher ozone concentrations are a risk to farmed species causing gross tissue damage and mortality. However, when O3 was applied in nanobubble form, we reduced the mortality. This implies that O3-NB treatment might be useful for controlling bacteria and increasing the DO without negative impact on shrimp.

When V. parahaemolyticus exponentially increased in the tanks, O3-NB treatment had a disinfection effect, which resulted in significantly lower mortality in shrimp (Experiment 2). Use of O3-NB may cause some gill changes but appears to be beneficial for reducing mortality in the presence of a pathogen. Given the results of Experiment 2 this benefit appears to be due to the reduction in bacteria. It might also be possible to the control bacteria levels without having to treat shrimp every day with ozone nanobubbles, which may further reduce the negative impact on culture species. The result of the study suggested that short frequent O3-NB treatments may be a solution to controlling bacteria, maintain DO content, and limit mortality associated with V. parahaemolyticus in shrimp culture.

The decline in the concentration of V. parahaemolyticus after O3 macrobubble and nanobubble treatments was consistent with our previous published study (Nghia et al., 2021). The application of O3 in nanobubble form seems slightly more efficient at killing bacteria than O3 delivered via macrobubble after 6 hours despite the fact the O3 was delivered for a shorter period than the O3-NB. A study by Ikehata & Li (2018) also found O3 and O3-NB had similar disinfection properties.

The disinfection properties observed in our treatment groups was attributed to high levels of ozone. Ozone levels were higher in the macrobubble group than in nanobubble group, which could explain the negative impact on the shrimp in this study. ORP was higher in the O3 treatments than the O2 treatments. ORP values of about 150–250 mV can promote nitrification, control bacterial growth, but above 350 mV starts killing bacteria (Cefas, 2010; Higgins, 2008). It should be noted that ORP is difficult to measure in water with high organic matter (Suslow, 2004) that reflects lower ORP values in Experiment 2. In water ozone and ORP values are positively correlated (Cefas, 2010; Nghi et al., 2018) so we can validate the accuracy of measurement by comparing their values at the same selected time points.

DO values massively increased during the O2 and O3 treatments, which was similar to what has been reported by Wang et al. (2018). The excessive levels of dissolved oxygen in our O2 treatment groups may have resulted in the gill pathology observed in this group of shrimp. The levels were well above the supersaturation level recommended for shrimp. The O3-NB treatment groups also had elevated DO but did not have severe gill pathology. O3 and O3-NB treatment both increased DO, ORP and decreased NO2 compared with all other treatments; and decreased NO3, NH3, and NH4 compared to the control. The results demonstrate the suitability of applying nanobubble to improve water quality for shrimp farming.


High numbers of shrimp with mild to moderate gill damage were observed in the O3, and O2 (macrobubble) treatments but these lesions were much less common in the O2-NB and O3-NB (nanobubble) treatment groups in both experiments. These results suggest shrimp gills were most affected by the O2 and O3 treatments when these gases were not delivered using nanobubble form. The delivery of these gases using nanobubbles appeared to reduce the negative effect of supersaturation. All treatments except the O3-NB treatment had a higher mortality at the end of the study than the negative control in the Experiment 2. Bacteria counts in water declined quickly in all tanks especially in the two O3 treatments suggested that besides increasing oxygen in the water, the use of O3-NB had the added benefit of reducing the bacterial load in the water and improved shrimp survival. Our study suggests that O3-NB can be applied to shrimp to effectively control V. parahaemolyticus, and hence, reduce AHPND. However, field studies are needed to confirm laboratory results and refine treatments for maximum effectiveness and use of O3-NB in intensive commercial shrimp farms.

Competing interests

No potential conflict of interest relevant to this article was reported.

Funding sources

This work was carried out with financial support (Grant No. 109050-002) from UK Government-Department of Health and Social Care (DHSC), Global AMR Innovation Fund (GAMRIF) and the International Development Research Centre, Ottawa, ON, Canada.


Not applicable.

Availability of data and materials

Upon reasonable request, the datasets of this study can be available from the corresponding author.

Ethics approval and consent to participate

This article does not require IRB/IACUC approval because there are no human and animal participants.



Agrawal A, Ng WJ, Liu Y. Principle and applications of microbubble and nanobubble technology for water treatment. Chemosphere. 2011; 84:1175-80


AquaPro Solutions. aQua+0.75M technical specification. Singapore: AquaPro Solutions. 2019.


Atkinson AJ, Apul OG, Schneider O, Garcia-Segura S, Westerhoff P. Nanobubble technologies offer opportunities to improve water treatment. Acc Chem Res. 2019; 52:1196-205


Azevedo A, Etchepare R, Calgaroto S, Rubio J. Aqueous dispersions of nanobubbles: generation, properties and features. Miner Eng. 2016; 94:29-37


Baird R, Bridgewater L. Standard methods for the examination of water and wastewater. 23rd ed Washington, DC: American Public Health Association. 2017.


Boyd CE, Fast AW. Pond monitoring and management.In In: Fast AW, Lester LJ, editors.editors Marine shrimp culture. Amsterdam: Elsevier. 1992; p p. 497-513


Buller NB. Bacteria from fish and other aquatic animals: a practical identification manual. Oxfordshire: CABI. 2004


Cefas. Ozone–potential application in depuration systems in the UK: cefas discussion document. Weymouth: Cefas Weymouth Laboratory. 2010.


Chaplin M. Nanobubbles (ultrafine bubbles) [Internet]. Water Structure and Science 2021cited 2021 Apr 19


Ebina K, Shi K, Hirao M, Hashimoto J, Kawato Y, Kaneshiro S, et al. Oxygen and air nanobubble water solution promote the growth of plants, fishes, and mice. PLOS ONE. 2013; 8e65339


Food and Agriculture Organization of the United Nations [FAO]. The state of world fisheries and aquaculture 2018–meeting the sustainable development goals. Rome: FAO. 2018.


Food and Agriculture Organization of the United Nations [FAO]. The state of world fisheries and aquaculture 2020. Sustainability in action. Rome: FAO. 2020.


Fregoso-López MG, Morales-Covarrubias MS, Franco-Nava MA, Ramírez-Rochín J, Fierro-Sañudo JF, Ponce-Palafox JT, et al. Histological alterations in gills of shrimp Litopenaeus vannamei in low-salinity waters under different stocking densities: potential relationship with nitrogen compounds. Aquac Res. 2017; 48:5854-63


Galang DP, Ashari AK, Sulmatiwi L, Mahasri G, Prayogo , Sari LA. The oxygen content and dissolved oxygen consumption level of white shrimp Litopenaeus vannamei in the nanobubble cultivation system. IOP Conf Ser Earth Environ Sci. 2019; 236:012014


Gonçalves AA, Gagnon GA. Ozone application in recirculating aquaculture system: an overview. Ozone Sci Eng. 2011; 33:345-67


Gurung A, Dahl O, Jansson K. The fundamental phenomena of nanobubbles and their behavior in wastewater treatment technologies. Geosyst Eng. 2016; 19:133-42


Higgins P. ORP management in wastewater as an indicator of process efficiency [Internet]. YSI 2013cited 2020 May 5


Ikehata K, Li Y. Ozone-based processes.In In: Ameta SC, Ameta R, editors.editors Advanced oxidation processes for wastewater treatment. London: Academic Press. 2018; p p. 115-34


Imaizumi K, Tinwongger S, Kondo H, Hirono I. Disinfection of an EMS/AHPND strain of Vibrio parahaemolyticus using ozone nanobubbles. J Fish Dis. 2018; 41:725-7


Jhunkeaw C, Khongcharoen N, Rungrueng N, Sangpo P, Panphut W, Thapinta A, et al. Ozone nanobubble treatment in freshwater effectively reduced pathogenic fish bacteria and is safe for Nile tilapia (Oreochromis niloticus). Aquaculture. 2021; 534:736286


Kirby BJ. Micro- and nanoscale fluid mechanics: transport in microfluidic devices. Cambridge: Cambridge University Press. 2010.


Lightner DV. A handbook of shrimp pathology and diagnostic procedures for diseases of cultured penaeid shrimp. Baton Rouge: World Aquaculture Society. 1996.


Ministry of Science and Technology [MOSTE]. National Standard TCVN 8710-19: 2019 on aquatic diseases - diagnostic procedure – part 19: acute hepatopancreatic necrosis disease in shrimp [Internet]. MARD 2019cited 2020 Feb 22


Nano Bubble Technologies. Aquaculture [Internet]. Nano bubble technologies 2020cited 2020 Feb 22


Nghi NH, Cuong LC, Dieu TV, Ngu T, Oanh DTY. Ozonation process and water disinfection. Vietnam J Chem. 2018; 56:717-20


Nghia NH, Van PT, Giang PT, Hanh NT, St-Hilaire S, Domingos JA. Control of Vibrio parahaemolyticus (AHPND strain) and improvement of water quality using nanobubble technology. Aquac Res. 2021; 52:2727-39


Parmar R, Majumder SK. Microbubble generation and microbubble-aided transport process intensification—a state-of-the-art report. Chem Eng Process. 2013; 64:79-97


Schroeder JP, Klatt SF, Schlachter M, Zablotski Y, Keuter S, Spieck E, et al. Impact of ozonation and residual ozone-produced oxidants on the nitrification performance of moving-bed biofilters from marine recirculating aquaculture systems. Aquac Eng. 2015; 65:27-36


Seki M, Ishikawa T, Terada H, Nashimoto M. Microbicidal effects of stored aqueous ozone solution generated by nano-bubble technology. In Vivo. 2017; 31:579-83


Sjogreen CA, Téllez DAL, Pérez JER, Hurtado PCP, Roa-Rojas J. Experimental study of nanobubbles in salt solutions. Rev Acad Colomb Cienc Exact Fis Nat. 2018; 42:41-8


Suslow TV. Oxidation-reduction potential (ORP) for water disinfection monitoring, control, and documentation. Oakland, CA: University of California. 2004


TCTS. Kết quả thực hiện năm 2021 và triển khai nhiệm vụ trọng tâm năm 2022 lĩnh vực thủy sản [Fisheries annual report-2021, Hanoi]. Hochiminh: TCTS. 2021.


Temesgen T, Bui TT, Han M, Kim TI, Park H. Micro and nanobubble technologies as a new horizon for water-treatment techniques: a review. Adv Colloid Interface Sci. 2017; 246:40-51


Tzaneva V, Bailey S, Perry SF. The interactive effects of hypoxemia, hyperoxia, and temperature on the gill morphology of goldfish (Carassius auratus). Am J Physiol Regul Integr Comp Physiol. 2011; 300:R1344-51


Wang L, Miao X, Ali J, Lyu T, Pan G. Quantification of oxygen nanobubbles in particulate matters and potential applications in remediation of anaerobic environment. ACS Omega. 2018; 3:10624-30


Weijs JH, Lohse D. Why surface nanobubbles live for hours. Phys Rev Lett. 2013; 110:054501


Weitkamp DE, Katz M. A review of dissolved gas supersaturation literature. Trans Am Fish Soc. 1980; 109:659-702


World Health Organization [WHO], Food and Agriculture Organization of the United Nations [FAO]. Risk assessment of Vibrio parahaemolyticus in seafood: interpretative summary and technical report. Geneva: WHO. 2011.


XLSTAT. The leading data analysis and statistical solution for Microsoft EXCEL® [Internet]. Addinsoftcited 2020 Feb 22