- Department of Life Sciences, Changzhi University, Changzhi, China
Due to its high efficacy as a wide-spectrum disinfectant and its potential for the degradation of pollutants and pesticides, ozone has broad application prospects in agricultural production. In this study, micro/nano bubble technology was applied to achieve a saturation state of bubble nutrient solution, including micro-nano oxygen (O2 group) and micro-nano ozone (O3 group) bubble nutrient solutions. The effects of these solutions on lettuce physiological indices as well as changes in the microbial community within the rhizosphere substrate were studied. The application of micro/nano (O2 and O3) bubble nutrient solutions to substrate-cultured lettuce plants increased the amount of dissolved oxygen in the nutrient solution, increased the lettuce yield, and elevated the net photosynthetic rate, conductance of H2O and intercellular carbon dioxide concentration of lettuce plants. Diversity analysis of the rhizosphere microbial community revealed that both the abundance and diversity of bacterial and fungal communities in the substrate increased after plant cultivation and decreased following treatment with micro/nanobubble nutrient solutions. RDA results showed that the microbial community in the S group was positively associated with EC, that in the CK and O2 groups exhibited a positive correlation with SC, and that in the O3 group displayed a positive correlation with CAT and POD. Overall, the implementation of micro/nanobubble generation technology in soilless substrates can effectively increase the lettuce growth and yield, and O3 had a more pronounced effect on lettuce yield and quality and the microbial community structure in the substrate than O2. Our study would provide a reference and theoretical basis for developing sustainable and green technology for promoting lettuce production and can be a promising alternative to conventional methods for improving crop yields.
1 Introduction
Currently, soil disinfection has become widely used in agriculture, which can reduce the occurrence of pests and diseases at the source and improve vegetable quality and yield (Gullino et al., 2022). The use of conventional chemical pesticides and disinfectants generally leads to detrimental effects (Rani et al., 2021). Improper utilization of pesticides can also adversely affect the quality of water, soil and air (Tudi et al., 2021). Consequently, uncomplicated and effective disinfection techniques used for decontaminating less contaminated soil is gaining traction. Recently, ozone has been considered an environmentally friendly method extensively applied in various fields, including shelf life extension of products, deodorization, food, water and air purification, industrial wastewater treatment (Li et al., 2009), preservation of vegetables and fruits (Miller et al., 2013), and rhizosphere soil sterilization (Díaz-López et al., 2022).
Ozone is an allotropic form of oxygen with a strong oxidation potential, and it can decompose into oxygen without harmful residues (Tomiyasu et al., 1985). Compared with alternative methods of soil disinfection, ozone sterilization exhibits rapid disinfection speed and robust sterilization efficacy without inducing pathogenic microorganism resistance (Botondi et al., 2021; Díaz-López et al., 2022). Ozone can target proteins, unsaturated lipids, respiratory enzymes and peptidoglycans in cell membranes, enzymes and nucleic acids in the cytoplasm, and proteins and peptidoglycans in sporangia and viral capsids (Khadre et al., 2006). Due to its rapid decomposition and antibacterial properties, ozone treatment can improve crop quality and antioxidant capacity (Piechowiak et al., 2022) and greatly increase the stress resistance, photosynthetic capacity and yield of tomatoes (Xu et al., 2021). It can also inactivate Escherichia coli and mycelia in dry figs, thereby reduce aflatoxin levels (Zorlugenç et al., 2008).
Soil cultivation is the predominant method for crop production and requires regular soil aeration and irrigation; however, numerous pathogenic microorganisms also proliferate in the soil, competing for nutrients with crops and consequently impeding crop growth and defense mechanisms (Das et al., 2022). Soilless cultivation has gradually become as a pivotal component of facility horticulture and the primary mode of factory production for crops owing to its inherent advantages, such as fertilizer conservation, water preservation, enhanced yield potential, superior quality output, and impeccable hygiene standards (Gruda, 2022). Soilless substrate can replace natural soil to anchor crop roots, and supply essential nutrients through drip or trickle irrigation to ensure optimal growth and development (Asaduzzaman et al., 2015; Fields et al., 2021). Moreover, soilless substrate cultivation not only prevents soil disturbance and reduces land pressure and degradation, but also promotes crop intensification through year-round production (Gebreegziher, 2023). Compared to traditional soil cultivation, the composition and concentration of nutrient solution and growth media can be precisely controlled in soilless cultivation, and the yield and quality of horticultural crops can also be greatly improved (Xu et al., 1995; Elia et al., 1999). It enables precise control of trace element concentrations, ensuring continuous exposure of roots to fortified nutrient solutions without soil interference. Furthermore, soilless cultivation can maximize the absorption, transport and accumulation of nutrients in the edible part of plant roots (Rouphael and Kyriacou, 2018). The utilization of soilless culture also facilitates the mitigation of soil-borne diseases and the optimization of plant nutrition (Gebreegziher, 2023).
Micro/nanobubble generation technology, developed at the end of the last century, has garnered increasing attention as a clean and efficient method for water treatment (Khuntia et al., 2012). Microbubbles are referred as bubbles with a diameter less than 100 microns and situated between micron- and nanoscale dimensions (Chang et al., 2024). They can generate through various mechanisms, such as rotational shear, pressure dissolution, electrochemistry, micropore pressure and mixed jet flow (Díaz-López et al., 2022). These bubbles possessed unique physical and chemical properties, including slow bubble rise in solution, prolonged residence time, enhanced solubility, expanded specific surface area, accelerated gas−liquid mass transfer rate (Wu et al., 2019), elevated interface point, and spontaneous generation of free radicals (Li et al., 2009). Recently, ozone micro/nanobubbles have emerged as a promising disinfection technology, facilitating prolonged aqueous-phase ozone reactivity. This innovative approach holds potential for preserving the freshness and quality of fruits and vegetables while safeguarding their overall nutritional value (Botondi et al., 2021). Moreover, it can also increase the ozone saturation concentration in water, leading to its decomposition into oxygen for crop growth and development, thereby mitigating hypoxia-related issues in crops (Fan et al., 2020).
Lettuce (Lactuca sativa L.) is one of the most commonly consumed leafy vegetables worldwide and widely used for soilless cultivation (Nerlich and Dannehl, 2021). Currently, hypoxia is primarily mitigated in crop soilless cultivation via water circulation by a pump to ensure continuous nutrient solution flow and enhance dissolved oxygen levels through the synergistic effects of water flow and air contact. However, the efficiency in large-scale hydroponic systems is relatively low, which poses challenges to the normal growth of plants. Therefore, it is necessary to develop new methods for irrigating crops and improving the oxygen content in water. The integration of soilless culture substrate and micro/nanobubble generation technology enables efficient and rapid dissolution of ozone into the nutrient solution, thereby augmenting the dissolved oxygen levels in the solution. In this study, lettuce plants were cultivated in a soilless substrate and irrigated with Yamazaki nutrient solution, micro/nano-ozone nutrient solution or micro/nano-oxygen nutrient solution (CK, O2 and O3). The effects on the physiological indices and rhizosphere microbial community diversity of lettuce were analyzed, providing a reference and theoretical basis for developing sustainable and green technology and can be a promising alternative to conventional methods for improving crop yields.
2 Materials and methods
2.1 Preparation of three types of nutrient solutions
(1) Japan Yamazaki lettuce nutrient solution (CK): The mother liquid of the nutrient solution was prepared according to the formula listed in Supplementary Table S1, mixed in a ratio of liquid A:B:C of 1:1:0.1 and diluted 100 times before use.
(2) Micro/nano oxygen bubble nutrient solution (O2): Air was compressed with an oxygen concentrator 8F-5AW device (Jiangsu Yuyue Medical Equipment Co., Ltd.), resulting in oxygen generation. Subsequently, the oxygen was conveyed through an intake pipe to a micro/nanobubble generator MF-5000 device (Shanghai Xingheng Technology Co., Ltd.) where it was allowed to interact with the CK nutrient solution for 20 minutes, forming a micro/nanooxygen bubble nutrient solution.
(3) Micro/nano-ozone bubble nutrient solution (O3): Air was compressed in the oxygen generator, resulting in oxygen generation. Subsequently, the produced oxygen was fed into the ozone generator VMUS-4 (Weihua Dike Beijing Science and Technology Co., Ltd.) for conversion into ozone. The concentration of gaseous-phase ozone could be determined using an online UV-300H ultraviolet ozone detector (Weihua Dike Beijing Science and Technology Co., Ltd.) and the OZONE DESTRUCT 30168-02 system (Weihua Dike Beijing Science and Technology Co., Ltd.). The generated ozone was then conveyed to the micro/nanobubble generator MF-5000 (Shanghai Xingheng Technology Co., Ltd.) through an intake pipe and allowed to interact with the CK nutrient solution for 20 minutes, thereby forming a micro/nano-ozone bubble nutrient solution. The ozone concentration in the water was measured using a GreenPrima PM8200CL ozone analyzer (GreenPrima Instrument Co., Ltd., Shanghai), which consists of a Bsens650 double platinum electrode sensor, a PM8200CL controller, and a BAF615 constant current trough.
2.2 Lettuce treatment and physical index measurements
Seeds of the Italian lettuce variety were obtained from China Vegetable Seed Technology Co., Ltd. (Beijing) and cultivated in a 4×8 seedling tray at 25 °C within a controlled greenhouse environment. Once the third leaf emerged, the seedlings were transplanted into Danish peat soil (PINDSTRUP 0-10 mm) substrate. The 36 lettuce plants with identical growth statuses were randomly divided into three groups. The plants were irrigated with 300 mL of micro/nano-ozone bubble nutrient solution (O3, ~30 mg/L), micro/nano-oxygen bubble nutrient solution (O2, ~30 mg/L), or Japan Yamazaki lettuce nutrient solution (CK) every two days. During the treatment, the temperature was maintained at 25 °C, and the pH of the nutrient solution was maintained at 6.2. The lettuce plants were harvested after 55 days, and the following physiological indicators were assessed. (1) The net photosynthetic rate, transpiration rate, conductance to H2O and intercellular carbon dioxide concentration of the lettuce plants were measured by an LI-6400XT portable photosynthesometer (Beijing Ecotek Technology Co., Ltd.) following the manufacturer’s instructions. (2) The rapid induction kinetic curve of chlorophyll fluorescence (OJIP curve) in lettuce leaves was assessed using a Handy-PE portable plant efficiency analyzer (Hansha Scientific Instruments Ltd.) to determine the chlorophyll fluorescence parameters. (3) The lettuce fresh weight was determined using a balance, followed by oven drying at 105 °C for 30 minutes. Subsequently, the dried plants were further subjected to drying at 80 °C for 48 hours before their dry weight was measured. (4) The nitrate content was determined according to the methods described by Wang et al. (2017a). The ascorbic acid content was determined using the 2,6-dichloro-isophenol method, while root activity was assessed using 2, 3, 5-triphenyltetrazolium chloride method (TTC) TTC method. (5) The plant height and spreading width of the lettuce plants were measured by a ruler, and the number of leaves was recorded.
2.3 Analysis of substrate physical and chemical properties
The rhizosphere substrate samples were collected after lettuce harvesting. A five-point sampling method was employed, and the substrates were gathered within a radius of 0-15 cm surrounding the root system (Zhao et al., 2024). Part of the samples were air-dried and used to determine physicochemical properties and enzyme activities. The remaining samples were immediately cryopreserved with liquid nitrogen and maintained at a temperature of -80 °C for amplicon sequencing. The samples included S, representing the original substrate; CK, representing the control group with nutrient solution irrigation; O2, representing the micro/nano oxygen bubble nutrient solution irrigation; and O3, representing the micro/nano ozone bubble nutrient solution irrigation. Each treatment consisted of 3 biological replicates.
Catalase (CAT) activity was measured with a Soil CAT Activity Assay Kit via spectrophotometry at a wavelength of 240 nm. The peroxidase (POD) activity was assessed by employing a soil-specific POD assay kit (Soil Peroxidase (S-POD) Activity Assay Kit, Spectrophotometer, Solarbio) using colorimetric analysis at a visible wavelength of 430 nm. The determination of urease activity (UE) was conducted using the Soil Urease (S-UE) Activity Assay Kit at a visible wavelength of 630 nm with a spectrophotometer. The sucrase (SC) activity was determined using a spectrophotometer with the Soil Saccharase (S-SC) Activity Assay Kit at a visible wavelength of 540 nm. The pH was determined using a pH meter (ST3100, Ohaus Instrument (Changzhou) Co., Ltd.). The soluble salt concentration (EC) was determined by a portable EC meter (TZS-ECW-G, Zhejiang Topu Yunnong Technology Co., Ltd.). The significant differences in all detection indicators were determined using ANOVA in R (v4.3.1), while pairwise t tests were used to assess the significant differences between pairwise samples.
2.4 DNA extraction, PCR amplification and sequencing
Total genomic DNA for each sample was extracted using the OMEGA Soil DNA Kit (M5635-02) (Omega Bio Tek, Norcross, GA, USA) following the manufacturer’s instructions and stored at -20 °C prior to further analysis. The quantity of the obtained DNA was measured using a NanoDrop NC2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and agarose gel electrophoresis.
PCR amplification of the V3-V4 region of the bacterial 16S rRNA gene was performed using the forward primer 338F (5’-ACTCCTACGGGAGGCAGCA-3’) and reverse primer 806R (5’-GGACTACHVGGGTWTCTAAT-3’). The fungal primers ITS5F (5’-GGAAGTAAAAGTCGTAACAAGG-3’) and ITS1R (5’-GCTGCGTTCTTCATCGATGC-3’) were used to amplify the V1 region of the ITS gene. Sample-specific 7-bp barcodes were incorporated into the primers for multiplex sequencing. The PCR mixture contained 5 μL of reaction buffer (5×) and high GC buffer (5×), 0.25 μL of Q5 high-fidelity DNA polymerase (5 U/μL), 2 μL (10 mM) of dNTPs, 1 μL (10 µM) of each forward and reverse primer, 2 μL of DNA template, and 8.75 μL of ddH2O. Thermal cycling consisted of initial denaturation at 98 °C for 30 s, followed by 25~30 cycles of denaturation at 98 °C for 15 s, annealing at 50 °C for 30 s, and extension at 72 °C for 30 s, with a final extension of 5 min at 72 °C. PCR amplicons were purified with an Axygen Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) and quantified using a Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA) on a microplate reader (BioTek, FLx800). After the individual quantification step, amplicons were pooled in equal amounts, and paired-end 2×250 bp sequencing was performed using the Illumina NovaSeq platform with a NovaSeq 6000 SP Reagent Kit (500 cycles) at Shanghai Personal Biotechnology Co., Ltd. (Shanghai, China).
2.5 Quality control and analysis of sequence data
The QIIME2 (v2019.4) protocol was used with slight modifications according to the official tutorials (https://docs.qiime2.org/2019.4/tutorials/) to conduct data analysis (Bolyen et al., 2018). After removing the primers using the cutadapt plugin and filtering, denoising, merging, and removing chimeras from the raw sequences using the DADA2 plugin, the obtained high-quality sequences that exhibited 100% similarity were subsequently clustered into amplicon sequence variants (ASVs) (Callahan et al., 2016). For the annotation of bacterial 16S rDNA and ITS, we used the Greengenes (v13.8) (DeSantis et al., 2006) and UNITE databases (v8.0) (Kõljalg et al., 2013) to assign taxonomic information. ASVs with abundances less than 0.001% of the total samples were removed, and the abundance matrix with rare ASVs removed was used for subsequent analyses. We used the QIIME feature-table rarefy function for data normalized for differences in sequencing depth, and the leveling depth was set to 95% of the minimum sample sequence size. We conducted a comparative analysis of the bacterial and fungal proportions and percentages across various taxonomic levels (domain, kingdom, phylum, class, order, family and genus) within each sample to further investigate variations in the composition of the substrate microbiota.
2.6 Species diversity analysis
The α diversity indices, including Chao1, Shannon (Shannon, 1948), observed species, Simpson, Pielou’s evenness, and Goods’ coverage, were calculated using QIIME2 and visualized as box plots to compare the species richness and evenness of bacteria and fungi. Using R and QIIME2 software, β-diversity analysis was performed using the Bray-Curtis metric to investigate changes in microbial community structure between samples, which can be visualized by PCoA. The unweighted pair-group method with arithmetic means (UPGMA) hierarchical clustering of samples according to Euclidean distance based on species composition profiles at the genus level was used. To compare the memberships and structures of communities between samples, heatmaps were generated with the top 20 ASVs using Mothur.
2.7 Correlation analysis between physicochemical properties and microbial communities
Based on the physicochemical parameters and the relative abundances of the bacterial and fungal communities, redundancy analysis (RDA) was conducted using Canoco 4.5 software.
2.8 Analysis of key microbiota
To identify key bacteria and fungi within different groups, we employed linear discriminant analysis effect size (LEfSe) to analyze characteristic strains using the default parameters (Segata et al., 2011). The Kruskal-Wallis test was applied with a significance threshold of 0.05 for variance analysis. Similarly, the Wilcoxon rank sum test was conducted with a significance threshold of 0.05. Moreover, discriminating features were determined based on a logarithmic LDA score threshold of 3.0.
3 Results
3.1 Effects of different nutrient solutions on plant height, degree of development and yield in lettuce
After irrigation with O2, O3 or CK nutrient solution, we measured the lettuce leaf count, plant height, and individual plant weight (Figure 1). Compared with the CK group, the O3 treatment significantly increased the number of lettuce leaves and promoted both the growth and development of the plants. Although the O2 group also exhibited some promotion effects on lettuce plants, these effects did not reach statistical significance (Figures 1B–D, p > 0.05). The application of micro/nanobubble nutrient solution for irrigating lettuce significantly increased the wight per plant (Figures 1E–G, p <0.01), with a more pronounced promoting effect observed in the O3 group than in the O2 group. Therefore, irrigation with micro/nanobubble nutrient solutions, especially those containing ozone, can effectively enhance lettuce production.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1393905/fpls-15-1393905-HTML/image_m/fpls-15-1393905-g001.jpg)
Figure 1 Effects of different irrigation patterns on plant height, spreading width and weight per lettuce plant. (A) Graph of the lettuce plants harvested from the three different groups. (B–G) Comparisons of leaf number, plant height, spreading width, root weight, leaf weight and whole plant weight between the three groups. The significance of differences was determined using ANOVA with pairwise t tests. Different letters on different columns indicate a significant difference between these two treatments (p < 0.05).
3.2 Effects of different nutrient solutions on root activity and ascorbic acid and nitrate contents of lettuce plants
There was a significant increase in the root vitality of the lettuce plants in the O2 and O3 groups compared to that in the CK group (p < 0.05); the root viability coefficient of the O3 group was significantly greater than that of the O2 group (Table 1). The ascorbic acid content of the lettuce plants was not significantly affected by the nutrient solutions. Compared with that in the CK group, the application of the micro/nanobubble nutrient solution resulted in a significant reduction in the nitrate content (p < 0.01), and the O3 group exhibited significantly lower nitrate levels than the O2 group. Therefore, during substrate cultivation, the application of micro/nanobubble generation technology to introduce oxygen or ozone into nutrient solution could effectively enhance root vitality and significantly mitigate nitrate accumulation in lettuce plants.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1393905/fpls-15-1393905-HTML/image_m/fpls-15-1393905-t001.jpg)
Table 1 Effects of microbubble technology in nutrient solution on the root activity, ascorbic acid content and nitrate content of lettuce plants.
3.3 Effects of different nutrient solutions on the photosynthetic parameters of the lettuce plants
Compared with that of the CK group, the net photosynthetic rate of the O2 group exhibited a slight increase but lacked statistical significance; however, the net photosynthetic rate of the O3 group significantly improved (Table 2). Moreover, both the O2 and O3 groups exhibited significantly increased transpiration rates, intercellular carbon dioxide (CO2) concentrations, and conductance to H2O (p < 0.05). Furthermore, the net photosynthetic rate, transpiration rate, intercellular carbon dioxide concentration, and conductance to H2O in the O3 group were significantly greater than those in the O2 group (p < 0.05). These results suggested that incorporating micro/nanobubble generation technology into soilless substrate cultivation processes was beneficial for enhancing leaf conductance to H2O, increasing the intercellular CO2 concentration, and substantially improving the lettuce transpiration rate.
3.4 Effect of micro/nanobubble nutrient solutions on the chlorophyll fluorescence parameters of lettuce plants
Compared to those of the CK group, the Sm, Vj, DIo/CSo and DIo/CSm values of the O3 and O2 groups decreased significantly; the potential photochemical activity (Fv/Fo), maximum photochemical efficiency (Fv/Fm), ABS/RC, TRo/RC, REo/RC, ETo/RC and ETo/CSo increased, respectively (Figures 2A, B). Our results indicated that irrigating lettuce with nutrient solution after micro/nano treatment enhances the efficiency of chlorophyll absorption and light energy transfer, especially the O3 group.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1393905/fpls-15-1393905-HTML/image_m/fpls-15-1393905-g002.jpg)
Figure 2 Effects of different micro/nanobubble nutrient solutions on the chlorophyll fluorescence parameters of lettuce leaves. (A, B) The electron transport and reaction center activity of the PSII receptor side of lettuce leaves. (C) OJIP curves of lettuce leaves in the CK, O2 and O3 groups. (D) Fluorescence O-P logarithmic curve. (E) OJIP curves of the fluorescence kinetics of the O2 and O3 groups after global (O~P) double standardization with the CK group. (F) Logarithmic OJIP curves of the fluorescence kinetics of the O2 and O3 groups after global (O~P) double standardization with the CK group.
As shown in Figure 2C, compared with that of the CK group, the initial fluorescence intensity (Fo) of the O3 and O2 groups decreased by 29% and 37%, respectively; the O3 group exhibited a 13% increase in the maximum fluorescence intensity (Fm), while the O2 group showed a 9% decrease. Additionally, when the PSII reaction centers were completely closed, the O3 group demonstrated a greater fluorescence yield than the O2 group did. Compared to the CK group, the fluorescence intensity of the J-spot was attenuated, while that of the I-spot and P-spot were enhanced. O3 irrigation resulted in a reduction in the fluorescence intensity of the J-spot, I-spot and P-spot. After standardizing the OJIP curve (Figure 2D), it was observed that the fluorescence intensity of the J and P points in the O3 and O2 groups was significantly greater than that the CK group. The results suggested that O3 and O2 groups can enhance electron transfer in both the PSII reaction center and acceptor side. Notably, significant changes were observed at points J (2 ms) and I (30 ms) (Figures 2E, F). Overall, the effects of O2 and O3 treatment on the OJIP curve of lettuce leaves were comparable, while the impact of O3 on the oxygen emission complex (OEC) of lettuce leaves exceeded that of O2.
3.5 Analysis of physical and chemical properties of the substrate
Figure 3 shows that the indices of pH, CAT activity, POD activity, UE activity, and SC activity increased across all the three groups irrigating with nutrient solution (O3, O2 and CK) compared to those of the S group. Additionally, a decrease in EC was observed compared to that in the S group. Notably, the EC of the O3 group exhibited a slight increase compared to that of both the CK and O2 groups (Figure 3A). Compared to those of the CK group, both the O3 and O2 groups exhibited a decrease in pH (Figure 3B). POD enzyme activity significantly increased, with the greatest increase observed in the O3 group (Figure 3C). CAT activity increased more in the O2 group than in the O3 group (Figure 3D). SC and UE activities were lower in both the O2 and O3 groups than the CK group, with a particularly significant decrease observed in SC activity within the O3 group (Figures 3E, F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1393905/fpls-15-1393905-HTML/image_m/fpls-15-1393905-g003.jpg)
Figure 3 Analysis of the physical and chemical properties of the substrate. Values of EC (A), pH (B), POD (C), CAT (D), SC (E) and UE (F) between the three groups. The significant differences were determined using ANOVA with pairwise t tests.
3.6 Microbial sequencing analysis
Totally, 12 libraries were sequenced, which obtained 1,663,178 and 907,842 raw bacterial and fungal sequences. After quality control (QC) and data filtering, 1,500,641 and 828,149 clean reads were obtained for analysis of bacterial and fungal communities (Supplementary Table S2). In the bacterial community, the numbers of ASVs in the S, CK, O2 and O3 groups were 6,625, 8,565, 8,373 and 7,293, respectively. The number of common ASVs across all four groups was 636, accounting for 9.6%, 7.43%, 7.60% and 8.72% of the total number of ASVs identified in each group, respectively (Figure 4A). The fungal community contained 195 ASVs in the S group, while those in the CK, O2 and O3 groups numbered 789, 738 and 685, respectively. There were 58 shared ASVs among these four groups, accounting for 29.74%, 7.35%, 7.86% and 8.47% of the total number of ASVs identified in each corresponding group, respectively (Figure 4B). Overall, the ASV increased in the substrate after lettuce cultivation and irrigation, but there was a significant difference in composition among the three groups. The total number of ASVs in both groups treated with O2 and O3 was lower than that in the CK group, indicating a reduced survival ability for microorganisms in these environments. Especially, the O3 group had a more noticeable impact on the total number of ASVs.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1393905/fpls-15-1393905-HTML/image_m/fpls-15-1393905-g004.jpg)
Figure 4 ASV comparison of microbial communities and the abundance profile of the top 20 genera between different groups. (A, B) Vann comparisons of ASVs identified among the bacterial and fungal communities. (C, D) The abundance profiles of the top 20 genera among four samples, and UPGMA clustering of the samples was conducted according to the Euclidean distance of the species composition data.
Clustering analysis based on the genus-level abundance of different classes in the microbial communities revealed that Acidothermus, Dokdonella, Dongia, Bryobacter, Candidatus_Solibacter, Pseudolabrys and 67-14 exhibited relatively greater abundances in the S group. Conversely, Subgroup 6, Candidatus_Kaiserbacteria, Sphingomonas, Devosia, Saccharimonadales, Mucilaginibacter, BIrii41, Haliangium, A4b, Bauldia, SWB02, Terrimonas and Opitutus displayed lower abundances than they did in the other three groups (CK, O2 and O3) (Figure 4C). The species compositions of the CK and O2 groups exhibited a relatively high degree of similarity. In terms of the fungal community, significant differences were observed in the distribution of the dominant bacteria, with Chrysosporium, Clitopilus, Mortierella, and Lecanicillium being relatively abundant in the S group. Conversely, Chrysosporium, Plectosphaerella, Olpidium, Epicoccum, Paramyrothecium, Capronia, Alternaria and Scytalidium were relatively highly abundant in the O2 group. Monilinia, Cladosporium, Mycosphaerella, Thermoascus, Penicillium, Malassezia and Phialemonium were relatively abundant in the O3 group. Additionally, the relative abundances of Alternaria, Scytalidium, Conlarium, and Trichoderma were greater in the CK group (Figure 4D). These results demonstrated that the relative abundance of the dominant bacteria was influenced by treatment with the micro/nanobubble nutrient solutions.
3.7 Alpha diversity of the microbial community
Compared with those of the CK group, the α diversity of S group exhibited lower Chao1, Simpson, Pielou’s evenness, Shannon, and observed species indices in the bacterial community (Figure 5A). Conversely, the CK group displayed higher values than the O2 and O3 groups. These findings suggested that lettuce cultivation enhanced bacterial diversity and that micro/nanobubble nutrient solution decreased this diversity. However, no significant difference was detected in terms of the α diversity (Figure 5B, p > 0.05). That the rarefaction curves tended to be flat indicated that the sequencing results effectively reflected the diversity present in the sample (Supplementary Figures S1A, B). The bacterial abundance curve exhibited a relatively gentle slope (Supplementary Figure S1C), indicating small differences in abundance between the ASVs and high uniformity in community composition. In contrast, the fungal community displayed a steep line (Supplementary Figure S1D), suggesting low uniformity in community composition. Hence, both the abundance and diversity of bacterial and fungal communities in the substrate increased after plant cultivation and decreased following treatment with micro/nanobubble nutrient solutions.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1393905/fpls-15-1393905-HTML/image_m/fpls-15-1393905-g005.jpg)
Figure 5 The α diversity of the microbial community among the different groups. (A) The α diversity of the bacterial community. (B) The α diversity of the fungal community, including Chao1, Good’s coverage, Simpson, Pielou’s evenness, Shannon and observed species indices. * indicate significant differences among the different groups (Dunn’s test, p < 0.05).
3.8 Analysis of the microbial community structure composition
To investigate the impact of micro/nanobubble nutrient solutions on the dominant bacterial community structure and fungi, we analyzed and compared the species composition of each group. Figure 6A showed the classification tree of the bacterial communities in the four substrate samples. There were mostly four phyla that were the most abundant, namely, Proteobacteria, Actinobacteria, Acidobacteria, and Chloroflexi. Among them, Alphaproteobacteria, Actinobacteria, Acidobacteriia, and Ktedonobacteria belong to these phyla. Moreover, two classes, Bacteroidia and Verrucomicrobiae, also showed high abundance. Figure 6B showed the classification tree of the fungal community, which mostly included two phyla, Ascomycota and Basidiomycota. Among them, Eurotiomycetes, Sordariomycetes, Dothideomycetes and Agaricomycetes had the highest abundances.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1393905/fpls-15-1393905-HTML/image_m/fpls-15-1393905-g006.jpg)
Figure 6 Classification tree and taxonomic composition of the microbial community structures among the different groups. (A, B) The classification tree of the detected bacteria and fungi. (C, D) The taxonomic composition of bacterial and fungal community structures at the phylum level.
Proteobacteria consistently dominated the bacterial communities across the three treatment groups (Figure 6C, CK, O2 and O3); however, variations in relative abundance were observed. The phyla Proteobacteria, Actinobacteria, Acidobacteria, Bacteroidetes, and Chloroflexi collectively accounted for 88% of the bacterial composition among all groups. Conversely, Verrucomicrobia, Patescibacteria, Gemmatimonadetes, Planctomycetes, and Firmicutes exhibited relative abundances below 5%. Compared to that in the CK group, the relative abundance of Acidobacteria decreased in both the O2 and O3 groups (Figure 6C). In the O2 group, there was a decrease in the proportion of Actinobacteria, while the relative abundances of Proteobacteria and Bacteroidetes increased. In the O3 group, there was a decreasing trend in the relative abundances of Proteobacteria and Bacteroidetes, while the proportions of Actinobacteria and Patescibacteria increased.
The composition of the fungal community was relatively simple compared to that of the bacterial community (Figure 6D). In terms of relative abundance, the top 8 phyla of samples included Basidiomycota, Ascomycota, Mortierellomycota, Chytridiomycota, Olpidiomycota, Mucoromycota, Aphelidiomycota, and Rozellomycota. Compared with those in the S group, the relative abundance of Basidiomycota and Mortierellomycota decreased, and that of Ascomycota increased in the substrate (CK, O2 and O3). Compared with than in the CK group, the relative abundance of Mortierellomycota in the O2 and O3 groups significantly decreased. Moreover, there was a decreasing trend in the relative abundance ratio of Basidiomycota and Chytridiomycota, while an increase was noted in the relative abundance ratio of Olpidiomycota. Compared to those in the CK group, the relative abundances of Ascomycota and Mucoromycota in the O2 group decreased. Conversely, in the O3 group, there was an increase in the relative abundance of Ascomycota, Mucoromycota, and Aphelidiomycota.
Principal Coordinate Analysis (PCoA) is a conventional unconstrained ranking technique that can be employed to depict variations or associations in the internal composition of microbial communities. As shown in Figure 7A, principal component 1 (PCo1) accounted for 50.6% of the bacterial community diversity within the samples, while PCo2 explained an additional 18.1%. When considering fungal community diversity (Figure 7B), PCo1 and PCo2 jointly contributed to approximately 61.4% of the observed variation. Figures 7C, D show that all three biological replicates from each treatment, except the O2 treatment, were clustered into one cluster at the genus level (the same was also true at the other levels), indicating that the micro/nanobubble nutrient solutions changed the structural distribution of bacteria and fungi in the substrate. The topological structure of the hierarchical clustering of bacteria and fungi exhibited remarkable similarity, with both the CK and O2 groups closely clustered together. The S group was positioned at the base of the topology, indicating that lettuce cultivation most significantly influenced substrate microorganisms. Intergroup difference analysis revealed that the distance among the three sampling points in the CK group was the smallest, and that the distance between CK and O2 or CK and O3 was much greater (Figures 7E, F). This further confirmed that the effects of the intergroup treatment outweighed the effects of the sampling deviations. Totally, these results indicated that the microbial and fungal communities in the rhizosphere substrate changed after treatment with two micro/nanobubble nutrient solutions and that the community composition was more influenced by O3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1393905/fpls-15-1393905-HTML/image_m/fpls-15-1393905-g007.jpg)
Figure 7 β diversity analysis of the microbial community structures. (A, B) PCoA plots of bacterial and fungal communities. (C, D) The panel on the left is a hierarchical clustering diagram in which samples are clustered according to their similarity of bacterial and fungal communities. The panel on the right is a stacked bar chart of the top 10 most abundant genera. (E, F) Analysis of the distance from O2 and O3 to CK based on their microbial abundance profiles.
3.9 Principal coordinate analysis and redundancy analysis
To further investigate the impact of environmental factors on the distribution of substrate microbial community components, such as CAT, POD, UE, SC, pH and EC, RDA was employed to explore the relationships between the composition and structure of the substrate microbial community and environmental variables. As shown in Figure 8A, it shows that the significant differences of the bacterial community structure between the S, CK, O2 and O3 groups were correlated with the environmental conditions. RDA1 accounted for 62.51% of the interpretive variance, while RDA2 explained an additional 12.58%. Together, these two components contributed to a cumulative interpretive degree of 75.09%, indicating a high level of explanatory power. The RDA of fungal communities revealed that the incorporation of environmental factors and other variables led to more pronounced dissimilarities in the community structure composition between the S, CK, O2 and O3 groups. The first component accounted for 73.91% and the second component 4.52% (Figure 8B). The CAT, POD, UE, SC and pH were found to be important environmental factors affecting both the bacterial and fungal community. For the bacterial community (Figure 8A), the S group was positively correlated with EC; the CK group positively correlated with pH, UE and SC, the O2 group exhibited a positive correlation with SC; and the O3 group displayed a positive correlation with CAT and POD. For the fungal communities (Figure 8B), the S group was positively correlated with EC; the CK group showed a positive correlation with SC; the O2 group exhibited a positive correlation with pH, UE and SC; and the O3 group displayed a positive correlation with SC, CAT, and POD.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1393905/fpls-15-1393905-HTML/image_m/fpls-15-1393905-g008.jpg)
Figure 8 RDA plots of microbial communities with environmental variables. (A) Bacterial communities. (B) Fungal communities. POD, CAT, UE, SC and EC.
3.10 LEfSe difference analysis
LEfSe analysis was employed to identify biomarkers in microbial communities exhibiting significant differences among the groups. For the bacterial communities, 38, 14, 25, and 21 biomarkers were characterized for the S, CK, O2 and O3 groups, respectively (Figure 9A). Actinobacteria exhibited the greatest difference in abundance in the S group, whereas Terrimonas and Deltaproteobacteria were identified as the most abundant bacterial taxa in the CK and O2 groups, respectively. Devosia emerged as the predominant taxon in the O3 group. In the fungal community, there were 11, 6, 6, and 18 biomarkers in S, CK, O2 and O3, respectively (Figure 9B). Among these groups, Clitopilus was identified as the most significant fungus in the S group, Fusarium exhibited the highest abundance in the CK group, Plectosphaerella showed significant enrichment in the O2 group, and Sordariomycetes emerged as the most significant fungus in the O3 group.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1393905/fpls-15-1393905-HTML/image_m/fpls-15-1393905-g009.jpg)
Figure 9 LEfSe analysis of microbial sequences. (A) Bacterial 16S rDNA sequences. (B) Fungal ITS rDNA sequences. The concentric circles depict the hierarchical classification levels ranging from phylum to genus. Each small circle at a specific level represents a taxonomic classification, with its size directly proportional to the relative abundance of that taxon. Hollow nodes represent taxa with no significant differences between groups, while colored nodes indicate significant differences between groups S, CK, O2 and O3, with higher abundances in samples represented by these colors. Letters were used to identify taxa that exhibited significant differences between groups. The legend on the right provides the names of the species denoted by letters. The size of each node corresponds to the average relative abundance of the taxon.
4 Discussion
Recently, the use of ozone water in agricultural production has significantly increased, including its application in the preservation of agricultural products (Shi et al., 2023), enhancement of fruit quality (Silva Neto et al., 2019), sterilization and disinfection, and disease prevention and control (Trinetta et al., 2011). By incorporating micro/nanobubble generation technology into agricultural practices, it not only conserves agricultural water but also enhances crop quality and yield (Liu et al., 2019). Due to the historical significance and economic value as a vital vegetable crop, employing micro/nanobubble generation technology alongside soilless substrate culture has immense potential for optimizing lettuce cultivation systems and increasing productivity.
4.1 Effects of micro/nano-ozone bubble nutrient solution on the growth, quality, and photosynthetic parameters of lettuce
Improving the lettuce quality and yield is a major concern in lettuce cultivation. Yang et al. (2017) reported that the use of glutamine as a nitrogen source in hydroponic cultivation can increase the accumulation of glycosylated flavonoids, ascorbic acid, and most amino acids, and improve the nutritional value of lettuces. Foliar application of selenium solution for lettuce plants can enhance their nitrogen utilization efficiency (Piñero et al., 2022). It has also been shown that irrigation with oxygenated water can increase crop yield (Zhou et al., 2019). In this study, compared with the CK group, the O2 and O3 treatments increased the lettuce yield (Figure 1G) and improved the root system vitality (Table 1). This treatment may facilitate the promotion of aerobic respiration in plant roots, thereby augment nutrient transfer (Sey et al., 2009; Liu et al., 2019), which can stimulate the growth of lettuce plants.
Ascorbic acid is abundant in fresh vegetables and fruits (Zheng et al., 2022), with lettuce being considered a moderate dietary source of this essential nutrient (Yang et al., 2021). In this study, the content of ascorbic acid in the O3 group was greater than the CK group. Yao et al. (2021) demonstrated that the utilization of micro/nano oxygen aeration technology can enhance the oxidation efficiency of ammonia and effectively eliminate NH4+ through improved oxygen transfer efficiency. Similar to Yao et al. (2021), the nitrate content in the lettuce leaves of both O2 and O3 groups decreased compared to the CK group, especially the O2 group (Table 1).
It has demonstrated that the introduction of oxygen can enhance crop yield, water use efficiency, and the leaf photosynthetic rate (Zhou et al., 2019). The net photosynthetic rate, transpiration rate, intercellular carbon dioxide concentration, and stomatal development in the O3 group were significantly greater than those in the O2 group (Table 2). When a nutrient solution treated with micro/nanobubble generation technology is applied for lettuce irrigation, it can improve the net photosynthetic rate, and promote the efficient absorption and transfer of light energy, thereby promote lettuce growth and development. In conclusion, the micro/nanobubble nutrient solution not only promotes lettuce growth but also enhances its quality, especially the O3 treatment.
4.2 Effect of micro/nano-ozone bubble nutrient solution on microbial community diversity in the substrate
The exceptional disinfectant and antibacterial properties of ozone stem from its potent oxidizing effect. It has demonstrated bactericidal activity against both gram-positive and gram-negative bacteria, in addition to its ability to effectively deactivate yeast and mold spores (Wang et al., 2021). Ozone micro/nanobubbles can effectively remove the antibiotics ciprofloxacin and levofloxacin from water (Babaee et al., 2023). Dissolved ozone exhibits greater antibacterial efficacy than ozone in its gaseous state and can effectively control bacteria, fungi, viruses, protozoa, and spores by impacting cellular membranes and other cellular constituents (Premjit et al., 2022; Shelake et al., 2022).
Ozone water irrigation may induce denaturation and deactivation of harmful substances within the plant rhizosphere, thereby enhancing plant growth. There was a significant increase in both the abundance and diversity of bacterial and fungal communities present in the substrate (Figure 5). The species number and diversity of bacteria and fungi in the substrates of the O2 groups were slightly lower than the CK group, while those in the O3 group were significantly lower. This may be due to ozone’s potent oxidation and antibacterial activity. Microbial community analysis revealed similar dominant phyla of bacteria and fungi among the sampled substrates; however, their relative abundances varied. Actinobacterium is pivotal in the decomposition of cellulose and various chemical compounds, thereby significantly contributing to the carbon cycle (Lewin et al., 2016). The relative abundance of Actinobacteria in both the CK and O2 groups exhibited low, whereas that in the O3 and S groups showed high. Bacteroidetes are the predominant bacteria in soil (Larsbrink and McKee, 2020) and can secrete a diverse array of carbohydrate-active enzymes. Compared to those in the S group, the relative abundances of Bacteroidetes in the CK, O2 and O3 groups increased. The abundance of Bacteroidetes in the O3 group was lower than that in the O2 group but was consistent with that in the CK group (Figure 6C). Previous studies have shown that bacteria of the Acidobacteria phylum are oligotrophic (Wang et al., 2017b) and that the abundance of Acidobacteria in soil is often greater under lower organic carbon conditions (Fierer et al., 2007). The relative abundance of Acidobacteria in the substrates of groups O2 and O3 decreased. Proteobacteria are copiotrophically attributed and favored by nutrient-rich conditions with high carbon content (Fierer et al., 2007; Newton and McMahon, 2010). In all the groups, more than 40% of the bacteria belonged to Proteobacteria. Chloroflexi is an anaerobic bacterium that undergoes parthenogenetic reproduction and primarily engages in carbon and nitrogen fixation (Logue and Lindström, 2010). The relative abundance of Chloroflexi remained consistent across all groups. The findings suggested that the introduction of oxygenated irrigation did not impact the population of beneficial bacteria but did positively influence nutrient accumulation in the substrate.
The relative abundance of Basidiomycota decreased in the CK, O2 and O3 groups compared to that in the S group, while there was a significant decrease in the relative abundance of Mortierellomycota and an increase in that of Ascomycota. The Basidiomycota phylum can be broadly categorized into saprophytic, symbiotic, and parasitic/pathogenic taxa based on their nutritional strategies. Saprophytic fungi can decompose organic matter such as wood, leaf litter or feces, while symbiotic fungi establish mutualistic associations with other organisms. Therefore, these fungi are pivotal in ecosystem functioning. Pathogenic or parasitic fungi can infect plants and other fungi (Schmidt-Dannert, 2016). The Mortierella fungus is a prevalent soil fungus that is vital in soil ecology and crop health. It facilitates the absorption of mineral elements by plant roots while inhibits pathogens (Miao et al., 2016; Li et al., 2018). A subset of saprophytic bacteria within the Ascomycota phylum is crucial in decomposing recalcitrant organic matter present in soil, thereby facilitating nutrient cycling. However, certain members of Ascomycota can also induce plant diseases such as root rot, stem rot, fruit rot and branch blight (Beimforde et al., 2014; Zhao et al., 2019).
The concentration of ozone in nutrient solutions may also impact crop growth and yield. Studies have shown that high levels of ozone can damage crops and vegetation (Rai et al., 2015) and can adversely affect plant growth, active ingredients and yield (Lee et al., 2022). For example, high concentrations of ozone can reduce the biomass of potatoes, cause leaf damage, and decrease the total number of potato tubers. Increased ozone concentrations can also impact the diversity of root-associated bacteria in crops (Feng et al., 2015). For instance, O3 has inhibitory effects on the microbial community in wheat rhizosphere soil (Li et al., 2012) and indirectly affects the composition of microbial communities in poplar rhizosphere soil negatively (Wang et al., 2024). An increase in O3 concentration can result in a decline in the biomass of fungi and gram-negative bacteria (Díaz-López et al., 2022), thereby perturbing the microbial community composition, particularly by reducing nitrifying bacteria and enhancing the relative abundance of fungal groups capable of denitrification. Elevated O3 not only impacts the structure of soil microbial communities but also has potential implications for their functional capabilities (Chen et al., 2019). In this study, a saturated oxygen and ozone micro/nanobubble nutrient solution was used. However, there is no clear categorization of the dissolved oxygen mass concentration in micro/nanobubble nutrient solutions, and further investigation is needed to determine the optimal gas concentration promoting plant growth and optimizing microbial community of substrate. We aimed to develop an irrigation frequency that is more suitable for meeting the requirements of plant growth while also being cost-effective, scientifically sound, and sustainable.
5 Conclusion
In summary, the application of micro/nanoscale (O2 and O3) bubble nutrient solution can promote photosynthesis and increase the yield of substrate-cultured lettuce plants. The majority physical and chemical indexes of the substrate also exhibited an increasing trend, while SC and UE activities decreased. The change in the rhizosphere microbial community with O2 and O3 treatment revealed that both the abundance and diversity of bacterial and fungal communities in the substrate increased after plant cultivation and decreased following treatment with micro/nanobubble nutrient solutions. The O3 group had a more pronounced impact on the microbial community structure in the substrate than the O2 group. Our study provides a method that combines micro/nanobubble generation technology with a soilless substrate to cultivate lettuce plants, which can promote robust seedling growth and substantially increase lettuce yield.
Data availability statement
The raw sequencing data have been deposited in the NCBI Sequence Read Archive (SRA) under accession number PRJNA1078579.
Author contributions
QZ: Data curation, Funding acquisition, Investigation, Methodology, Validation, Writing – original draft. JD: Formal analysis, Funding acquisition, Validation, Writing – original draft, Writing – review & editing. SL: Data curation, Investigation, Writing – original draft. WL: Formal analysis, Writing – original draft. AL: Conceptualization, Visualization, Writing – review & editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the fund for Shanxi “1331 Project” Key Subjects Construction (1331KSC), Technological Innovation Programs of Higher Education Institutions in Shanxi (2021L523, 2023L332).
Acknowledgments
We thank Ms. Liying Meng (Shanxi University) for her technical support for uploading raw data to NCBI database. We also thank AJE for language editing and Shanghai Personal Biotechnology Co., Ltd for sequencing service.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2024.1393905/full#supplementary-material
Supplementary Figure 1 | Rarefaction curves and ASV rank-abundance curves of different substrate samples. (A) Rarefaction curves of bacterial communities. (B) Rarefaction curves of fungal communities. (C) ASV rank-abundance curves of the bacterial communities. (D) ASV rank-abundance curves of fungal communities. The horizontal coordinate is the sequencing depth, and the vertical coordinate is the median value of the α diversity index calculated 10 times, which is presented in a box plot.
Supplementary Table 1 | The mother liquor formula of Japanese Yamazaki nutrient solution for lettuce.
Supplementary Table 2 | Sequencing information for different samples.
References
Asaduzzaman, M., Saifullah, M., Mollick, A.K.M.S.R., Hossain, M. M., Halim, G. M. A., Asao, T. (2015). “Influence of Soilless Culture Substrate on Improvement of Yield and Produce Quality of Horticultural Crops,” in Soilless Culture - Use of Substrates for the Production of Quality Horticultural Crops. (London, UNITED KINGDOM: IntechOpen Limited), 1–32. doi: 10.5772/58679
Babaee, Y., Saghravani, S. F., Feizy, J. (2023). Predicting ciprofloxacin and levofloxacin decomposition utilizing ozone micro-nano bubbles through the central composite design method. Environ. Technol., 1–13. doi: 10.1080/09593330.2023.2260123
Beimforde, C., Feldberg, K., Nylinder, S., Rikkinen, J., Tuovila, H., Dörfelt, H., et al. (2014). Estimating the Phanerozoic history of the Ascomycota lineages: Combining fossil and molecular data. Mol. Phylogenet. Evol. 78, 386–398. doi: 10.1016/j.ympev.2014.04.024
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C., Al-Ghalith, G. A., et al. (2018). QIIME 2: reproducible, interactive, scalable, and extensible microbiome data science. PeerJ Preprints 6, e27295v27292. doi: 10.7287/peerj.preprints.27295v2
Botondi, R., Barone, M., Grasso, C. (2021). A review into the effectiveness of ozone technology for improving the safety and preserving the quality of fresh-cut fruits and vegetables. Foods 10, 748. doi: 10.3390/foods10040748
Callahan, B. J., Mcmurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., Holmes, S. P. (2016). DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Chang, J., Zhang, H., Liu, F., Cui, K. (2024). Exploring the mechanism of micro-nano bubble water in enhancing the mechanical properties of sulfoaluminate cement-based materials. Construction Building Materials 411, 134400. doi: 10.1016/j.conbuildmat.2023.134400
Chen, Z., Maltz, M. R., Cao, J., Yu, H., Shang, H., Aronson, E. (2019). Elevated O3 alters soil bacterial and fungal communities and the dynamics of carbon and nitrogen. Sci. Total Environ. 677, 272–280. doi: 10.1016/j.scitotenv.2019.04.310
Das, P. P., Singh, K. R. B., Nagpure, G., Mansoori, A., Singh, R. P., Ghazi, I. A., et al. (2022). Plant-soil-microbes: A tripartite interaction for nutrient acquisition and better plant growth for sustainable agricultural practices. Environ. Res. 214, 113821. doi: 10.1016/j.envres.2022.113821
DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72, 5069–5072. doi: 10.1128/AEM.03006-05
Díaz-López, M., Siles, J. A., Ros, C., Bastida, F., Nicolás, E. (2022). The effects of ozone treatments on the agro-physiological parameters of tomato plants and the soil microbial community. Sci. Total Environ. 812, 151429. doi: 10.1016/j.scitotenv.2021.151429
Elia, A., Serio, F., Gonnella, M., Santamaria, P. (1999). Growing nitrate free endive in soilless systems. Acta Hortic. 481, 267–272. doi: 10.17660/ActaHortic.1999.481.28
Fan, W., An, W.-G., Huo, M.-X., Yang, W., Zhu, S.-Y., Lin, S.-S. (2020). Solubilization and stabilization for prolonged reactivity of ozone using micro-nano bubbles and ozone-saturated solvent: A promising enhancement for ozonation. Separation Purification Technol. 238, 116484. doi: 10.1016/j.seppur.2019.116484
Feng, Y., Yu, Y., Tang, H., Zu, Q., Zhu, J., Lin, X. (2015). The contrasting responses of soil microorganisms in two rice cultivars to elevated ground-level ozone. Environ. pollut. 197, 195–202. doi: 10.1016/j.envpol.2014.11.032
Fields, J. S., Owen, J., Lamm, A., Altland, J. E., Jackson, B. E., Zheng, Y., et al. (2021). Soilless substrate science: a North American needs assessment to steer soilless substrate research into the future. Acta Hortic. 1317, 313–318. doi: 10.17660/ActaHortic.2021.1317.36
Fierer, N., Bradford, M. A., Jackson, R. B. (2007). Toward an ecological classification of soil bacteria. Ecology 88, 1354–1364. doi: 10.1890/05-1839
Gebreegziher, W. G. (2023). Soilless culture technology to transform vegetable farming, reduce land pressure and degradation in drylands. Cogent Food Agric. 9, 2265106. doi: 10.1080/23311932.2023.2265106
Gruda, N. S. (2022). Advances in soilless culture and growing media in today’s horticulture—An editorial. Agronomy 12, 2773. doi: 10.3390/agronomy12112773
Gullino, M. L., Garibaldi, A., Gamliel, A., Katan, J. (2022). Soil disinfestation: from soil treatment to soil and plant health. Plant Dis. 106, 1541–1554. doi: 10.1094/PDIS-09-21-2023-FE
Khadre, M. A., Yousef, A. E., Kim, J. G. (2006). Microbiological aspects of ozone applications in food: A review. J. Food Sci. 66, 1242–1252. doi: 10.1111/j.1365-2621.2001.tb15196.x
Khuntia, S., Majumder, S. K., Ghosh, P. (2012). Microbubble-aided water and wastewater purification: a review. Rev. Chem. Eng. 28, 191–221. doi: 10.1515/revce-2012-0007
Kõljalg, U., Nilsson, R. H., Abarenkov, K., Tedersoo, L., Taylor, A. F. S., Bahram, M., et al. (2013). Towards a unified paradigm for sequence-based identification of fungi. Mol. Ecol. 22, 5271–5277. doi: 10.1111/mec.12481
Larsbrink, J., McKee, L. S. (2020). Bacteroidetes bacteria in the soil: Glycan acquisition, enzyme secretion, and gliding motility. Adv. Appl. Microbiol. 110, 63–98. doi: 10.1016/bs.aambs.2019.11.001
Lee, J. K., Kwak, M. J., Jeong, S. G., Woo, S. Y. (2022). Individual and interactive effects of elevated ozone and temperature on plant responses. Horticulturae 8, 211. doi: 10.3390/horticulturae8030211
Lewin, G. R., Carlos, C., Chevrette, M. G., Horn, H. A., Mcdonald, B. R., Stankey, R. J., et al. (2016). Evolution and ecology of actinobacteria and their bioenergy applications. Annu. Rev. Microbiol. 70, 235–254. doi: 10.1146/annurev-micro-102215-095748
Li, F., Chen, L., Redmile-Gordon, M., Zhang, J., Zhang, C., Ning, Q., et al. (2018). Mortierella elongata’s roles in organic agriculture and crop growth promotion in a mineral soil. Land Degradation Dev. 29, 1642–1651. doi: 10.1002/ldr.2965
Li, P., Takahashi, M., Chiba, K. (2009). Degradation of phenol by the collapse of microbubbles. Chemosphere 75, 1371–1375. doi: 10.1016/j.chemosphere.2009.03.031
Li, X., Deng, Y., Li, Q., Lu, C., Wang, J., Zhang, H., et al. (2012). Shifts of functional gene representation in wheat rhizosphere microbial communities under elevated ozone. ISME J. 7, 660–671. doi: 10.1038/ismej.2012.120
Liu, Y., Zhou, Y., Wang, T., Pan, J., Zhou, B., Muhammad, T., et al. (2019). Micro-nano bubble water oxygation: Synergistically improving irrigation water use efficiency, crop yield and quality. J. Clean. Prod. 222, 835–843. doi: 10.1016/j.jclepro.2019.02.208
Logue, J. B., Lindström, E. S. (2010). Species sorting affects bacterioplankton community composition as determined by 16S rDNA and 16S rRNA fingerprints. ISME J. 4, 729–738. doi: 10.1038/ismej.2009.156
Miao, C.-P., Mi, Q.-L., Qiao, X.-G., Zheng, Y.-K., Chen, Y.-W., Xu, L.-H., et al. (2016). Rhizospheric fungi of Panax notoginseng: diversity and antagonism to host phytopathogens. J. Ginseng Res. 40, 127–134. doi: 10.1016/j.jgr.2015.06.004
Miller, F. A., Silva, C. L. M., Brandão, T. R. S. (2013). A review on ozone-based treatments for fruit and vegetables preservation. Food Eng. Rev. 5, 77–106. doi: 10.1007/s12393-013-9064-5
Nerlich, A., Dannehl, D. (2021). Soilless cultivation: dynamically changing chemical properties and physical conditions of organic substrates influence the plant phenotype of lettuce. Front. Plant Sci. 11, 601455. doi: 10.3389/fpls.2020.601455
Newton, R. J., McMahon, K. D. (2010). Seasonal differences in bacterial community composition following nutrient additions in a eutrophic lake. Environ. Microbiol. 13, 887–899. doi: 10.1111/j.1462-2920.2010.02387.x
Piechowiak, T., Grzelak-Błaszczyk, K., Sójka, M., Balawejder, M. (2022). One-time ozone treatment improves the postharvest quality and antioxidant activity of Actinidia arguta fruit. Phytochemistry 203, 113393. doi: 10.1016/j.phytochem.2022.113393
Piñero, M. C., Otálora, G., Collado-González, J., López-Marín, J., Del Amor, F. M. (2022). Effects triggered by foliar selenium application on growth, enzyme activities, mineral nutrients and carbohydrates in lettuce under an aquaculture system. Plant Physiol. Biochem. 180, 1–8. doi: 10.1016/j.plaphy.2022.03.028
Premjit, Y., Sruthi, N. U., Pandiselvam, R., Kothakota, A. (2022). Aqueous ozone: Chemistry, physiochemical properties, microbial inactivation, factors influencing antimicrobial effectiveness, and application in food. Compr. Rev. Food Sci. Food Saf. 21, 1054–1085. doi: 10.1111/1541-4337.12886
Rai, R., Agrawal, M., Kumar Choudhary, K., Agrawal, S. B., Emberson, L., Büker, P. (2015). Application of ethylene diurea (EDU) in assessing the response of a tropical soybean cultivar to ambient O3: Nitrogen metabolism, antioxidants, reproductive development and yield. Ecotoxicol. Environ. Saf. 112, 29–38. doi: 10.1016/j.ecoenv.2014.10.031
Rani, L., Thapa, K., Kanojia, N., Sharma, N., Singh, S., Grewal, A. S., et al. (2021). An extensive review on the consequences of chemical pesticides on human health and environment. J. Clean. Prod. 283, 124657. doi: 10.1016/j.jclepro.2020.124657
Rouphael, Y., Kyriacou, M. C. (2018). Enhancing quality of fresh vegetables through salinity eustress and biofortification applications facilitated by soilless cultivation. Front. Plant Sci. 22, 1254. doi: 10.3389/fpls.2018.01254
Schmidt-Dannert, C. (2016). Biocatalytic portfolio of basidiomycota. Curr. Opin. Chem. Biol. 31, 40–49. doi: 10.1016/j.cbpa.2016.01.002
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12, R60. doi: 10.1186/gb-2011-12-6-r60
Sey, B. K., Manceur, A. M., Whalen, J. K., Gregorich, E. G., Rochette, P. (2009). Root-derived respiration and nitrous oxide production as affected by crop phenology and nitrogen fertilization. Plant Soil 326, 369–379. doi: 10.1007/s11104-009-0018-x
Shannon, C. E. (1948). A mathematical theory of communication. Bell Syst. Tech. J. 27, 623–656. doi: 10.1002/bltj.1948.27.issue-4
Shelake, P. S., Mohapatra, D., Tripathi, M. K., Giri, S. K., Kate, A., Kumar, M. (2022). Inactivation of Aspergillus Niger and Erwinia carotovora in onion (Allium cepa L.) bulbs subjected to pulsed ozone treatment. Postharvest Biol. Technol. 192, 111998. doi: 10.1016/j.postharvbio.2022.111998
Shi, J., Cai, H., Qin, Z., Li, X., Yuan, S., Yue, X., et al. (2023). Ozone micro-nano bubble water preserves the quality of postharvest parsley. Food Res. Int. 170, 113020. doi: 10.1016/j.foodres.2023.113020
Silva Neto, O. P. D., Pinto, E. V. D. S., Ootani, M. A., Silva Junior, J. L. D., Lima, J.L.D.S.B., Sousa, A. E. D. D. (2019). Ozone slows down anthracnose and increases shelf life of papaya fruits. Rev. Bras. Fruticultura 41, e–439. doi: 10.1590/0100-29452019439
Tomiyasu, H., Fukutomi, H., Gordon, G. (1985). Kinetics and mechanism of ozone decomposition in basic aqueous solution. Inorganic Chem. 24, 2962–2966. doi: 10.1021/ic00213a018
Trinetta, V., Vaidya, N., Linton, R., Morgan, M. (2011). A comparative study on the effectiveness of chlorine dioxide gas, ozone gas and e-beam irradiation treatments for inactivation of pathogens inoculated onto tomato, cantaloupe and lettuce seeds. Int. J. Food Microbiol. 146, 203–206. doi: 10.1016/j.ijfoodmicro.2011.02.014
Tudi, M., Daniel Ruan, H., Wang, L., Lyu, J., Sadler, R., Connell, D., et al. (2021). Agriculture development, pesticide application and its impact on the environment. Int. J. Environ. Res. Public Health 18, 1112. doi: 10.3390/ijerph18031112
Wang, J., Song, Y., Ma, T., Raza, W., Li, J., Howland, J. G., et al. (2017b). Impacts of inorganic and organic fertilization treatments on bacterial and fungal communities in a paddy soil. Appl. Soil Ecol. 112, 42–50. doi: 10.1016/j.apsoil.2017.01.005
Wang, Q., Yang, Q., Zhang, M., Ma, J., Qu, L. (2024). Effects of ozone stress on rhizosphere soil of poplar seedlings. Forests 15, 205. doi: 10.3390/f15010205
Wang, Q. H., Yu, L. J., Liu, Y., Lin, L., Lu, R. G., Zhu, J. P., et al (2017a). Methods for the detection and determination of nitrite and nitrate: A review. Talanta 165, 709–720. doi: 10.1016/j.talanta.2016.12.044
Wang, T., Yun, J., Zhang, Y., Bi, Y., Zhao, F., Niu, Y. (2021). Effects of ozone fumigation combined with nano-film packaging on the postharvest storage quality and antioxidant capacity of button mushrooms (Agaricus bisporus). Postharvest Biol. Technol. 176, 111501. doi: 10.1016/j.postharvbio.2021.111501
Wu, C., Li, P., Xia, S., Wang, S., Wang, Y., Hu, J., et al. (2019). The role of interface in microbubble ozonation of aromatic compounds. Chemosphere 220, 1067–1074. doi: 10.1016/j.chemosphere.2018.12.174
Xu, H.-L., Gauthier, L., Gosselin, A. (1995). Effects of fertigation management on growth and photosynthesis of tomato plants grown in peat, rockwool and NFT. Scientia Hortic. 63, 11–20. doi: 10.1016/0304-4238(95)00791-Q
Xu, J.-P., Yu, Y.-C., Zhang, T., Ma, Q., Yang, H.-B. (2021). Effects of ozone water irrigation and spraying on physiological characteristics and gene expression of tomato seedlings. Horticulture Res. 8, 180. doi: 10.1038/s41438-021-00618-8
Yang, X., Feng, L., Zhao, L., Liu, X., Hassani, D., Huang, D. (2017). Effect of glycine nitrogen on lettuce growth under soilless culture: a metabolomics approach to identify the main changes occurred in plant primary and secondary metabolism. J. Sci. Food Agric. 98, 467–477. doi: 10.1002/jsfa.8482
Yang, X., Gil, M. I., Yang, Q., Tomás-Barberán, F. A. (2021). Bioactive compounds in lettuce: Highlighting the benefits to human health and impacts of preharvest and postharvest practices. Compr. Rev. Food Sci. Food Saf. 21, 4–45. doi: 10.1111/1541-4337.12877
Yao, G.-J., Ren, J.-Q., Zhou, F., Liu, Y.-D., Li, W. (2021). Micro-nano aeration is a promising alternative for achieving high-rate partial nitrification. Sci. Total Environ. 795, 148899. doi: 10.1016/j.scitotenv.2021.148899
Zhao, Q., Dong, J., Yan, Z., Xu, L., Liu, A. (2024). Effect of cucumber continuous monocropping on traditional Chinese medicine residue through analysis of physicochemical characteristics and microbial diversity. Agronomy 14, 709. doi: 10.3390/agronomy14040709
Zhao, S., Qiu, S., Xu, X., Ciampitti, I. A., Zhang, S., He, P. (2019). Change in straw decomposition rate and soil microbial community composition after straw addition in different long-term fertilization soils. Appl. Soil Ecol. 138, 123–133. doi: 10.1016/j.apsoil.2019.02.018
Zheng, X., Yuan, Y., Huang, B., Hu, X., Tang, Y., Xu, X., et al. (2022). Control of fruit softening and Ascorbic acid accumulation by manipulation of SlIMP3 in tomato. Plant Biotechnol. J. 20, 1213–1225. doi: 10.1111/pbi.13804
Zhou, Y., Zhou, B., Xu, F., Muhammad, T., Li, Y. (2019). Appropriate dissolved oxygen concentration and application stage of micro-nano bubble water oxygation in greenhouse crop plantation. Agric. Water Manage. 223, 105713. doi: 10.1016/j.agwat.2019.105713
Keywords: lettuce, micro/nano-oxygen nutrient solution, micro/nano-ozone bubble nutrient solution, physiological index, microbial community diversity
Citation: Zhao Q, Dong J, Li S, Lei W and Liu A (2024) Effects of micro/nano-ozone bubble nutrient solutions on growth promotion and rhizosphere microbial community diversity in soilless cultivated lettuces. Front. Plant Sci. 15:1393905. doi: 10.3389/fpls.2024.1393905
Received: 29 February 2024; Accepted: 27 March 2024;
Published: 11 April 2024.
Edited by:
Zhen Jiao, Zhengzhou University, ChinaReviewed by:
Hui Yang, Hunan Academy of Agricultural Sciences (CAAS), ChinaZhigang Liu, Jiangsu University, China
Xuanzhen Li, Henan Agricultural University, China
Copyright © 2024 Zhao, Dong, Li, Lei and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ake Liu, akeliu@126.com
†These authors have contributed equally to this work