Two Food Waste By-Products Selectively Stimulate Beneficial Resident Citrus Host-Associated Microbes in a Zero-Runoff Indoor Plant Production System

The global production of food waste is a far-reaching problem with sizable financial, ethical, social, and environmental costs. Over 66 million tons of food waste is produced annually in the United States alone. This waste can be converted into valuable digestate by-products that promote a circular economy within agri-food systems. The present work investigated the use of two liquid digestates of microaerobic fermentation from mixed food waste and beer mash, respectively, as biostimulants for non-bearing citrus plants (nursery stock) grown in a zero-runoff greenhouse system with recirculating irrigation. The digestates' impact on the structure and diversity of the microbiota was determined on the irrigation water, soil, leaves, roots, and rhizosphere of citrus plants. A combination of culture-dependent (selective media) and culture-independent approaches (Next-Generation Sequencing) was used to assess the composition of the microbial communities and to single out the presence of foodborne pathogens. Our results suggest that the use of digestates is safe (i.e., no human or plant pathogens were present in the digestates or enriched in the plant production system following amendments). Digestates application to the irrigation water reduced the bacterial diversity within 24–48 h and selectively and significantly stimulated beneficial resident host-associated microorganisms (Pseudomonas putida) by two to three orders of magnitude. Carbon dynamics were analyzed in the nutrient solutions by measuring dissolved organic carbon and characterizing carbon species through gas chromatography-electron ionization-mass spectrometry. Our results indicate that dissolved organic carbon in the recirculating irrigation water spikes after each digestate amendment and it is quickly metabolized by bacteria, plateauing 24 h after application. Soil carbon, nitrogen, and nutrient dynamics were also analyzed, and results suggest that digestates increased the concentration of some plant nutrients in soils without causing a surge of potentially toxic elements. This study represents a proof-of-concept for the safe re-use of organic wastes, from farming and consumers, in agriculture. Implementing this type of integrated plant production system could reduce the environmental impact of food waste and benefit the public by improving soil health, reducing agricultural footprint, and increasing crop fitness by deploying a method based on a circular economy and sustainable food production approaches.


INTRODUCTION
The current use of natural resources, accelerated by anthropogenic pressures at many levels, threatens to exceed Earth's carrying capacity [World Health Organization Regional Office for Europe (WHO), 2018]. The way we produce, choose, consume, and waste food, poses a serious threat to our planet. As much as one-third of all food produced is either lost before it reaches a store or wasted before it reaches consumers' plates [Food and Agriculture Organization of the United Nations (FAO), 2013]. Food waste has deleterious effects on the environment, our society, and the economy. In the United States (US) alone, 40-50 percent (%) of all food is thrown away, with an estimated annual cost of about US $750 billion every year (FAO, 2019). Most of this waste is not recycled but instead ends up in landfills, taking up to 21% of America's landfill volume and contributing 9% to the total US greenhouse gas emissions [United States Environmental Protection Agency (EPA), 2017; Friedlingstein et al., 2019]. With agriculture utilizing 70% of the water throughout the world (FAO, 2017), food waste also represents a significant loss of freshwater resources. Nevertheless, issues associated with food waste expand beyond its environmental impact: with millions of low-income consumers struggling to gain access to nutrient-rich foods on a daily basis including fresh fruits, vegetables, and meat, while many smallholder farmers are living on the margins of food insecurity; hence better understanding of how to minimize food waste would have a significant societal impact as well. Improving the efficiency of the food supply chain can lower production costs to farmers, while also lowering the cost of food to the consumer and thus increase access to both (FAO, 2015;Horton et al., 2019;Wright et al., 2019).
In recent decades, the need to transition from a linear economy based on the "take-make-consume-dispose" model toward a circular economy has received increased attention worldwide (Ellen MacArthur Foundation, 2013;Wright et al., 2019;Gao et al., 2020). In the circular economy concept, the value of resources is maintained in the system as long as possible, and products are re-claimed to extract as much value as possible before safely returning them to the biosphere. Food waste is rich in energy, water, and nutrients. With the physical, chemical and biological degradation of agricultural soils intensifying, and the demand of mineral fertilizers and pesticides continuously rising, nutrients from organic wastes have potential to be recycled into high-value agricultural products and used as biostimulants (e.g., Plant Probiotic Microorganisms, bioprotectants, biocontrollers, and biofertilizers) (Xu and Geelen, 2018). This approach requires transforming organic wastes through treatments and their reuse in agriculture. While a variety of methods and food-waste management practices exist (e.g., aerated composting) the development of new processes for nutrient-reuse was identified as one of the main challenges in waste management within a circular economy (Mahanty et al., 2017;Ioannidou et al., 2020). Most of the current composting processes are time-, space-, and energy-consuming, and have sizable greenhouse gas (GHG) footprints, and poor carbon sequestration efficiency (Green and Popa, 2011). In recent years, anaerobic digestion for treatments of organic materials from agricultural and food waste have become more common [Monlau et al., 2015; USDA's Agricultural Marketing Services (AMS), 2017]. Controlled bioprocess producing valuable liquid or solid digestates bioproducts have also been developed (Cerda et al., 2019;Elsharkawy et al., 2019;Giwa et al., 2019). Within this concept, the Bokashi fermentation became increasingly popular worldwide, especially for smallholder farmers across Asia and Latin America (Xu et al., 2001;Gómez-Velasco et al., 2014). Unlike compost, Bokashi is mainly driven by microbial anaerobic to microaerophilic fermentation toward acidic conditions. Bokashi's positive impact on crop productivity and the environment has been shown (Möller and Müller, 2012) however, little is known on the effect of Bokashi liquid digestates on soil and plant microbiota and soil carbon dynamics (Christel, 2017). Efficient and cost-effective reuse of waste by-products depends on their relative agronomic efficiency, their stability, pollution risks and safety but more importantly on how they impact the soil and plant microbiome and their intricate network (Pandey et al., 2016;Bernal, 2017;Vanotti et al., 2019).
In recent years, sustainable methods have been developed for increasing crop production while simultaneously reducing chemical inputs such as synthetic chemical fertilizers and pesticides (Lampridi et al., 2019). Harnessing resident soil and plant-associated microbes will be instrumental to make agriculture production systems more sustainable in the longterm (Lugtenberg and Kamilova, 2009;Mendes et al., 2011;Berendsen et al., 2012;Turner et al., 2013;Schlaeppi and Bulgarelli, 2015;Ciancio et al., 2016;Lareen et al., 2016;Berg et al., 2017;Riera et al., 2017;Amaya-Gómez et al., 2020). Plant growth-promoting microbes (PGPM) (i.e., bacteria and fungi), have been widely reported for their beneficial roles in protecting plants from biotic and abiotic stresses and improving crop productivity (Ruzzi and Aroca, 2015;Abhilash et al., 2016). However, many issues need to be addressed before PGPMs can be used on large-scale to support sustainable food production. One of the main challenges is the inconsistent performance of bio-inoculants under field conditions (Timmusk et al., 2017) primarily due to their inability to maintain critical population thresholds (Ruzzi and Aroca, 2015;Berg et al., 2017) and to offset the antagonism from the native microbiome (Acea et al., 1988;van Veen et al., 1997). It has become evident that PGPMs will not be effective for an extended time unless the environment is modified to make it more conducive for their growth and survival (Pagliaccia et al., 2007(Pagliaccia et al., , 2008Zhou et al., 2020). On the other hand, it might not be necessary to introduce foreign microorganisms (i.e., inoculant) since soils already hold a vast reservoir of microbes, and plants can enrich their own root-associated microbiome which can function as biocontol agents and PGPMs (Pagliaccia et al., 2007;Lebeis, 2015;Reinhold-Hurek et al., 2015;Lareen et al., 2016;Timmusk et al., 2017;Zhang et al., 2017;Köhl et al., 2019;Raymaekers et al., 2020).
In California, the citrus industry is valued at $3.4 billion with a total economic impact of over $7 billion (Babcock, 2018). California's citriculture is facing historical drought and water management challenges [University of California Agriculture and Natural Resources (UCANR), 2015; Environmental Defense Fund (EDF), 2016] as well as pest (e.g., Asian citrus psyllid, Diaphorina citri Kuwayama) and disease (Huanglongbing, HLB) threats. As result, programs such as the "Citrus Pest and Disease Prevention Program" and the "Citrus Nursery Stock Pest Cleanliness Program" have transformed pest management efforts and citrus nursery production in the state [da Graça et al., 2016;Babcock, 2018;Dai et al., 2019;McRoberts et al., 2019; California Department of Food and Agriculture (CDFA), 2020]. This demanding and dynamic agricultural environment has created a unique opportunity to adopt alternative citrus fruit and tree production systems.
To develop a sustainable next-generation citriculture, it is imperative to capitalize on soil carbon dynamics, soil and plant microbiota, and microbial traits beneficial to citrus, the environment, and the economy. In this context, the development of novel practices that deploy targeted manipulation of microorganisms to stimulate soil life and harness internal ecosystem processes offers a promising approach to enhance agricultural sustainability (Philippot et al., 2013;Bender et al., 2016) and adopt circular economy-based solutions. As demonstrated in other cases, controlled management of liquid digestates through drip irrigation systems in open field or protected citrus orchards could be implemented to enhance native soil beneficial microbial communities. Furthermore, improvement of indoor citrus nursery production could come from adopting zero run-off irrigation systems where 100% of irrigation/fertigation water is recycled and reused however, these are understudied areas (Bailey and Lazarovits, 2003;Bonilla et al., 2012;Cao et al., 2016;Ciancio et al., 2016;Barzee et al., 2019).
Recent studies have focused on understanding the role of the citrus microbiome and its impact on plant health and productivity (Trivedi et al., 2010(Trivedi et al., , 2012. It has been suggested that members of the citrus microbiome influence the interaction between Candidatus Liberibacter asiaticus (the bacterium associated with HLB) and citrus (Blaustein et al., 2017;Wang et al., 2017;Blacutt et al., 2020). Several studies have investigated why some citrus trees identified in severely HLBdiseased citrus groves (with the same genotype as symptomatic trees), evaded the HLB disease. Preliminary published research suggests that the HLB-escaped trees differ from the diseased ones because of their associated native microbial community composition (Sagaram et al., 2009). The microbial community of escaped plants seems to be enriched in beneficial traits as compared to those of symptomatic trees (Trivedi et al., 2011(Trivedi et al., , 2012Riera et al., 2017).
Taken together these elements led to the theory that native soil and plant root-associated microorganisms, endowed with plant growth promoting capability and digestate-metabolizing abilities, can be used to recycle food waste while supporting sustainable agroecosystems with reduced synthetic chemical input. We formulated the following hypothesis for our study: (1) digestate products are suitable for delivery to citrus through drip irrigation lines; (2) digestates are safe to use (i.e., no plant or human pathogens); and (3) digestates will selectively increase certain bacterial groups (such as Pseudomonads) that are known to improve plant health. In this study, we tested two types of digestates, one derived from mixed post-consumer food waste (FW) and one from beer mash (BM).
To unravel the potential use of digestates in zero-runoff indoor citrus plant production, we investigated the effects of FW and BM liquid digestates products as biostimulants. Their impact on the structure and diversity of the microbiota in the recirculating nutrient solution and leaves, roots, and rhizosphere of citrus plants, was studied using a combination of culture-dependent (selective media substrates) and cultureindependent approaches (next-generation sequencing). The taxonomic composition of the microbial communities and the bacteria isolates enriched by FW and BM applications, were identified, and their beneficial traits investigated. Simultaneously, potting soil and nutrient solution chemical analyses were performed to determine whether FW and BM amendments significantly changed soil carbon and nitrogen content over the course of the treatments. Using gas chromatography-electron ionization-mass spectrometry (GC/EI-MS) we demonstrated the dependence of microbial growth on available dissolved organic carbon in the nutrient solution and identified overall variation in dissolved organic carbon speciation.

Materials, Analyses, and Processing
Two digestate products, one obtained from mixed food waste (FW) and one from beer mash (BM), were used as amendments (applied through drip irrigation) to the recirculating nutrient solution (Figure 1). The two microaerobic fermentation digestates, provided by River Road Research, were produced under a controlled anaerobic Bokashi-like fermentation technology in controlled conditions (Popa et al., 2017). Digestates were stored at 4 • C in anaerobic conditions in FIGURE 1 | On the left, citrus seedlings in drip-irrigated hydroponic units in a temperature-controlled greenhouse. On the right, simplified layout of a recirculating hydroponic system. FW and BM digestates, at 4 ml/L of nutrient solution, were amended to the reservoir at 7 days interval. collapsible plastic containers. We selected these two digestates to test because they are commonly available wastes in southern California urban surroundings; however, the spent grain in the beer mash represents a uniform feedstock while food waste composition is more variable and could contain undesirable elements such as sodium from processed food. Food Waste Digestate samples were sent via next day delivery to A&L Western Agricultural Labs, Modesto CA, an agricultural and environmental testing laboratory, for fertilizer and metal analyses. Fertilizer analyses were performed using several different methods (Please see Supplementary Table 3). Metal analyses were performed using an EPA Method SW846-6010 except for mercury for which the method used was EPA Method SW846-7471A. Fertilizer and metal analysis results are reported in Table 1 and Supplementary Table 3. Dissolved organic carbon species in the two digestates were analyzed by GC/EI-MS at UCR, and results are reported in Supplementary Table 1.
A&L Western Agricultural labs tested the two digestates for fecal coliforms and Salmonella presence. Further analysis using culture-dependent (King's B, Nutrient Agar, and Tryptic Soy Agar selective media substrates) and culture-independent approaches (next-generation sequencing) were carried out at UCR to assess the presence of pathogens and any active bacteria in the digestates. Upon arrival, each digestate was plated out on King's B, Nutrient Agar, and Tryptic Soy Agar selective media substrates to identify any possible alive culturable bacteria residual. Both stocks were also collected onto sterile nitrocellulose membranes of two different sizes, 0.22 and 8 µm (Millipore, nitrocellulose). The filter-concentrated bacterial cells were then processed as described below to extract DNA and carry out Illumina sequencing of 16S rRNA gene amplicons. Laboratory tests for salmonella (A&L Labs) were negative, while fecal coliform counts were less than the reporting limit of 3.0 MPN/g. Our lab culturing on the digestates were negative. NGS results are discussed and reported in Supplementary Figure 2.

Zero-Runoff Cultural Systems
Experiments were conducted in temperature-controlled greenhouses (20-32 • C) using 9 separate recirculating zerorunoff units with top drip irrigated systems. Pineapple sweet orange seedlings (Citrus sinensis L. Osbeck) grown in 1-gallon plastic pots containing a potting soil mix (made of 33% coarse silica sand, 33% redwood shavings and 33% fine coir) were used as a host plant to assess the change of citrus, soil and nutrient solution microbiota after waste by-products amendments through the drip-irrigated system. Each zero-runoff system (located on a bench ca. 40 centimeters (cm) above ground level) consisted of an upper plastic trough, for plant growth, and a lower reservoir where the nutrient solution is kept and recycled (Figure 1).
Each trough held 2 pots with 3 citrus seedlings each. The nutrient solution was pumped from the reservoir to the trough and distributed to each plant via drip tubing with 7.5 liters per hour (L/h) emitters. Standard oxygenation of the nutrient solution was accomplished by air pumps bubbling air into the nutrient solution 20 min every 2 h throughout the day. The excess nutrient solution from each pot drained by gravity back to the reservoir and was recirculated. Plants were irrigated for 30 minutes (min) twice a day between 8:00 A.M. until 8:00 P.M. The solution pH was adjusted manually to 5.8 ± 0.4 with H 3 PO 4 and KOH, every week. The electrical conductivity (EC 1.6 ± 0.2 mS/cm 2 ) and the nutrient solution volume in reservoirs (24 L) were monitored once a week and levels maintained as required. The final element concentrations in the nutrient solution were as follows: Total nitrogen, N 250 (Ammoniacal nitrogen, NH 4 86.9; Nitrate Nitrogen, NO 3 150; Urea Nitrogen 13); Phosphate, P 2 O 5 59.5; potassium oxide, K 2 0 238.1; boron, B 0.31; copper, Cu 0.31; iron, Fe 1.25; manganese, Mn 0.63; molybdenum, Mo 0.15; and zinc, Zn 0.63 µg/ml.
Citrus seedlings were obtained from the Citrus Clonal Protection Program (CCPP) at the University of California, Riverside. One-year-old seedlings were placed into the zerorunoff systems described above. Plants were kept for 4 days on the recirculating zero-runoff systems before the first treatment to allow the bacteria from the plant root and rhizosphere to be dispersed into the nutrient solution reservoirs. There were three experiments in which the reservoirs were amended with different treatments. Each experiment consisted of three treatments, 1) FW by-product, 2) BM by-product, and non-amended control nutrient solution, with 3 replications per treatment and 18 total plants (6 plants for each replication) per treatment. FW and BM digestates, at 4 ml/L of nutrient solution, were first applied 5 days after the start of the experiment and reapplied at 7 days interval for 63 days for a total of 9 treatments applications. In all greenhouse experiments, the by-products were amended to the nutrient solution 30 min before the start of the irrigation. Between experiments, each zero-runoff system was sterilized with 10% bleach and rinse with water.

Culture Dependent Characterization of Bacterial Populations in Nutrient Solutions
Ten-ml samples of the nutrient solution were collected at different intervals at 0, 24, and 48 hours (h) after the waste byproducts amendment from the reservoir of each treatment. After serial 10-fold dilutions, aliquots were plated in quadruplicates,

Bacterial Collection and Identification in Nutrient Solutions and Bacterial Beneficial Traits
To investigate the effect of the digestates amendments on the selection of specific culturable bacteria in the recirculating nutrient solution, 231 single bacterial colonies were isolated from the nutrient solutions from the greenhouse experiments. Fortytwo colonies of representative predominant bacteria were then selected for DNA extraction and sequencing to provide genetic identification. Thirty-four isolates were collected from FW and BM treated and nine from non-treated nutrient solutions. For the selection, Petri dishes containing 50-100 single bacterial colonies were chosen, and single colonies were randomly picked from a different sector of each plate. Each isolated colony was purified by repeated streaking onto King B agar. Purified isolates were stored at room temperature on KB agar plates or maintained for long-term storage at −70 • C in sterile water supplemented with 15% glycerol. Single colonies were picked, mixed, and vortexed in 100 µl sterile distilled deionized water, out of which 1.25 µl of the mixture was used as DNA template in PCR, amplifying subunits of ribosomal RNA genes or 16S rDNA (dos Santos et al., 2019). Forward 27F (5 ′ -AGAGTTTGATCMTGGCTCAG-3 ′ ) and reverse 1492R (5 ′ -TACGGYTACCTTGTTACGACTT-3 ′ ) primers were used for 16S rDNA PCR amplification (Suzuki and Giovannoni, 1996). The PCR mixture consisted of 1× PCR buffer, MgCl 2 (2.5 mM), 250 µM each dNTP, 0.4 µM each primer, and 2.5 U of DNA polymerase (Jumpstart TM Taq DNA polymerase, Sigma-Aldrich, USA). The thermal cycling conditions used were as follows: 94 • C for 5 min; 35 cycles of 94 • C for 20 s, 48 • C for 20s, and 72 • C for 1 min; then 72 • C for 10 min. The PCR products were purified using a QIAquick PCR purification kit (Qiagen, Chatsworth, CA, USA) and then separated in 2.0% agarose gel electrophoresis. Resulting DNA bands were ∼1,400 bp large. The purified PCR products were then sequenced using the chain termination method (Sanger sequencing). Nucleotide sequences of 43 PCR amplicons were edited and assembled utilizing the Vector NTI software package (Thermo Fisher Scientific, USA). The consensus sequences of individual contig (5 total; GenBank accession numbers MW144948-MW144952) with different lengths were subjected to a remote BLAST search (megablast) in the NCBI nr/nt database using an Entrez query with default setting. The top species IDs with their accession number, % coverage and identity are reported in Table 2.
To acquire preliminary information on potential beneficial traits possessed by the most abundant bacteria isolated from the treated nutrient solution (listed in Table 2) the Pathosystems Resource Integration Center (PATRIC) database was used to Bacteria were isolated from different treated nutrient solutions (Food Waste, Beer Mash, and Control with no amendment). Contig designation indicates the contig number in our studies along with its GenBank accession number. Contig composition includes information on the number of isolates and from which treatment the bacteria were isolated from, followed by Query Length in base pairs (bp), the top bacterial ID, the scores, and percentages of coverage and identity (%).
search within the 16S rDNA gene sequences. Briefly, consensus sequences obtained from the 43 16S rDNA-sequences described above were used to carry out de novo functional annotation analysis with PATRIC and obtain reference genomes that best match the 16S region. The reference whole genomes were then analyzed with the PATRIC Protein family sorter to obtain the heatmap of relevant genes. Results are reported in Supplementary Figure 1. The beneficial traits investigated here were the one typically reported to be involved in phosphate solubilization, siderophore production, and iron acquisition, VOC production, osmoprotection, and osmotic tolerance, phytohormone production, antagonism, and nutrient competition . Because members of the genus Pseudomonas dominated the bacterial community at 24 and 48 h, but the ITs analysis reported the most abundant Pseudomonads as "Unknown, " we further investigated these Amplicon sequence variants (ASV)s by performing a BLAST search against the nucleotide database using default parameters. The results are reported in Table 3.

Nutrient Solution Samples Collection for Culture-Independent Characterization of Microbial Communities
Following the preliminary results from the bacterial culturing on KB and TSA media, samples of the nutrient solution were collected at different times (0, 24, and 48 h after treatments) for 4 weeks (on 1st, 2nd, 8th, and 9th week) to extract bacterial DNA and provide next-generation sequencing (NGS) datasets of the bacteria associated with recirculating nutrient solution. One liter of nutrient solution from the reservoirs from each treatment was collected onto sterile nitrocellulose membranes of two different sizes, 0.22 and 8 µm (Millipore, nitrocellulose), using a sterile vacuum aspirator filtering apparatus. The filterconcentrated bacterial cells were resuspended in 10-30 ml buffer 50 mM Tris-HCl pH 9, 20 mM EDTA, 400 mM NaCl, 0.75 M sucrose, and stored until extraction procedure in −80 • C freezer. Before extraction, the filters were sonicated in their original resuspension buffers by placing the 50 ml falcon tubes in cold water for 40 min. Before sonication, ice water in the Branson CPX952518R sonicator (Branson Ultrasonics TM , Danbury, CT, USA) was degassed for 5 min. During sonication, the frequency was set on 40 kHz. Genomic DNA was extracted from 200 µl of bacterial suspension using the PowerSoil R DNA Extraction Kit (MO BIO, Carlsbad, CA, USA). The purity and the concentration of the DNA template were assessed with an Infinite M1000 Pro (Tecan, Männedorf, Switzerland). The final DNA was eluted in 100 µl of DNA Elution Buffer. The DNA was aliquoted into three 1.5 ml Eppendorf DNA LoBind Tubes with 32 µl in each tube and dried using a SpeedVac concentrator (Thermo Fisher Scientific). The dried DNA was stored at −20 • C until further use for bacterial Illumina library construction.

Citrus Leaves, Roots, and Soil Samples Collection for Culture-Independent Characterization of the Microbial Community
Two weeks after the last digestates amendment (when the total bacterial population in all nutrient solutions reservoirs had reached 10 4 -10 5 CFU/ml), plants were taken out from the zerorunoff system and leaves, roots, and soil samples taken. At least 10 mature and fully expanded leaves were randomly collected from each pot (containing 3 plants of Citrus sinensis each) and pooled into 50 ml falcon conical tubes. Feeder roots were sampled from two sides of the tree after the root system was taken out from the plastic pot and topsoil was removed. Feeder roots were pulled apart, shaken to remove soil, and placed into 50 ml conical tubes. Between each sampled tree, the clippers and forceps were cleaned by spraying a 20% household bleach solution on them, and then wiping them clean with paper towels. Researchers also changed their gloves between each potted tree that was sampled. Thirty grams (g) of potting soil samples taken ∼10 cm from the base of the trunk, was randomly collected from each pot, and pooled into 50 ml falcon tubes (Corning). All samples were immediately placed in a cooler with ice and then stored at 4 • C. All tissue samples were prepared and lyophilized within 24 h from sample collection. Samples were not surface-sterilized; thus, our microbial NGS datasets and analyses included both epiphytes and endophytes. Leaf samples were cut into 1 to 2 cm long pieces and ∼6 g of leaves were placed into 50 ml conical tubes and stored at −80 • C; root samples were rinsed with autoclaved purified water Bacteria were identified from different treated nutrient solutions (Food Waste, Beer Mash and Control with no amendment). Bacterial ID (Accession), relative abundance in percentage (%), E-value, Coverage (%) and % Identity values are reported. % identity values are from analyses using BLAST (NCBI) where coverage was at least 30%. Pseudomonas putida JB (CP026675.1) and Pseudomonas putida strain DLL-E4 (CP045553.1) were identified from the control treatment but their relative abundance was only 0.53 and 0.74%, respectively.
(Barnstead Mega-Pure System MP-6a, Thermo Fisher Scientific, Waltham, MA, USA) and ∼5 g of rinsed root tissue was placed into 50 ml conical tubes and stored at −80 • C. Tissue samples were then lyophilized with a bench-top freeze dryer (Labconco FreeZone 4.5L, Kansas City, MO, USA) for 16 to 20 h. The lyophilized tissues were stored at −80 • C. DNA from freeze-dried leaves and roots were obtained as described previously by Ginnan et al. (2018). The purity and the concentration of the DNA template were assessed as described above. The DNA was aliquoted into three 1.5 ml Eppendorf DNA LoBind Tubes with 32 µl in each tube and dried using a SpeedVac concentrator (Thermo Fisher Scientific). The dried DNA was stored at −20 • C until further use for bacterial Illumina library construction. DNA from 100 mg of freeze-dried soil was extracted using the MoBio PowerSoil R DNA kit (MoBio Laboratories Inc., CA, USA) and eluted as described above. . PCR primers ITS-1507F (GGTGAAGTCGTAACAAGGTA) and ITS-23SR (GGGTTBCCCCATTCRG) were used to target a small portion of the small subunit, the large subunit, and the hypervariable ITS region, with the reverse primers including a 12 bp barcode (Hunt et al., 2006;Frank et al., 2008). Thermal cycling parameters were 94 • C for 5 min; 35 cycles of 94 • C for 20 s, 56 • C for 20 s, and 72 • C for 40 s, and followed by 72 • C for 5 min. PCR products were purified using the Qiagen QIAquickPCR Purification Kit (Qiagen, Valencia, CA, USA).

Library Construction, Sequence Processing, and Bioinformatics Workflow
DNA sequencing was performed using an Illumina MiSeq (Illumina, Inc, San Diego, CA, USA). Sequences were demultiplexed using Qiime 1.9.1 (Caporaso et al., 2012). Further quality control of the reads was conducted through the DADA2 v 1.14.0 pipeline (Callahan et al., 2016) including removal of PhiX, quality control (maxEE=1, min length=135bp, maxN=0), dereplicating, and creation of the Amplicon sequence variants (ASV) table. In the dataset, there were 11,119 ASVs in the first experiment and 12,861 ASVs in the second experiment initially identified. After removing sequences which did not belong to the kingdom bacteria or had an unknown phylum identification and did not occur in the dataset at least 3 times (singletons and doubletons) 8,855 and 10,264 ASVs remained. Samples that did not have at least 1,000 sequences were removed leaving a total of 77 and 63 samples, respectively, in the two experiments. Taxonomy identification was assigned using custom scripts and BLAST. Phyloseq v 1.30.3 (McMurdie and Holmes, 2013) was used to assist in the graphical representation of the data. For taxonomy barplots, data were aggregated to the genus level and transformed to relative abundance. Taxa that did not occur > 1% relative abundance were removed. Shannon's index was used to assess alpha diversity. Statistical differences for the alpha diversity and the taxonomic data were assessed with Kruskal-Wallis and pairwise Wilcoxon tests. Beta diversity was assessed using Bray-Curtis dissimilarity, which was depicted using Principal Coordinates Analysis (PCoA), and statistically assessed by performing Adonis tests [999 permutations]. The bacterial sequences have been deposited in the National Center for Biotechnology Information (NCBI)'s Sequence Read Archive (SRA) under bioproject PRJNA670208.

Dissolved Organic Carbon (DOC) Analysis of the Nutrient Solution and Amendments
Samples of nutrient solution were taken from all 9 replicates-25 ml total per zero-runoff system, as well as samples of fresh nutrient solution before it went into the system, and samples of the original liquid amendment for each treatment -FW and BM. Samples were taken before the amendment, after amendment, and after 24, 48, and 72 h on weeks 1 and 4. Immediately after sampling, all samples were stored on ice and filtered using 0.22 µm syringe filters to filter out microorganisms and root tips, leaving only dissolved organic carbon compounds. 3 ml samples were analyzed by a DOC analyzer.

Gas Chromatography/Electron Ionization-Mass Spectrometry (GC/EI-MS) Analysis
Nutrient solution and amendments were analyzed by GC/EI-MS to characterize dissolved organic carbon species. A protocol similar to Chen et al. (2019) was followed with some modifications. 3 ml of filtered samples were dried using N 2 gas evaporation and then trimethylsilylated with 100 µL of N,O-bis(trimethylsilyl)trifluoroacetamide (BSTFA; 99 µL) with trimethylchlorosilane (TMCS; 1 µL) (≥99%, SUPELCO) and 50 µL anhydrous pyridine (≥99%, EMD Millipore) at 70 • C for 1 h. A blank (control) made of 50 µL of anhydrous pyridine and 100 µL of BSTFA was used. This derivatization procedure converts hydroxyl functional groups into volatile trimethylsilyl (TMS) derivatives that can be analyzed with GC/EI-MS (Agilent 6890N GC coupled with 5975 MSD). An Agilent 7683B series autosampler was used to inject 5 µL of each derivatized sample in splitless mode. A J&W Scientific DB-5 column (30 m × 0.25 mm i.d., 0.25 µm film) was used to separate the TMS derivatives. The carrier gas used was helium at a flow rate of 1 ml/min. The GC ran for an overall time of 65.17 min at the following temperature program: initial temperature held at 60 • C for 1 min, temperature ramp of 3 • C min−1 up to 200 • C, temperature held at 200 • C for 2 min, temperature ramp of 20 • C min−1 up to 310 • C, and finally temperature held at 310 • C for 10 min.; a solvent delay time of 10.5 min was programmed. The GC inlet temperature was set at 250 • C and detector at 280 • C. The MS scan was performed in the m/z range of 50-550. GC/EI-MS data, including peak ID, peak height and area, and likely compound or molecule match, was processed using Agilent GC/MSD ChemStation software. We then used the O/C vs. H/C ratios of elemental formulas to get an overview of the compound categories within each amendment and generated van Krevelen diagrams that allow for classifying compounds general compound classes such as lipids, proteins, amino sugars, carbohydrates, lignin, tannins, condensed aromatics, and unsaturated hydrocarbons (Minor et al., 2014;Rivas-Ubach et al., 2018). Dissolved organic carbon species characterized in Beer Mash and Food digestates GC/EI-MS analyses are reported in Supplementary Material.

Carbon, Nitrogen, and C:N Ratio Changes in the Potting Soil
Five grams of soil were sampled on Week 1 and Week 10, in composite from each of the 3 replicates for Control, Beer Mash, and Food Waste, then sieved (2 µm), dried, and ground for analysis. Three mg samples of each of the replicates from each of the treatments were prepared (9 samples per week, 18 samples total) for analysis of total C, total N, and C:N using the Costech ECS 4010-Delta V Plus Isotope Ratio Mass Spectrometer. These standards were used: Acet, Glycine, Peach, EM Soil, USGS64, and USGS66, and samples were combusted completely in EA at 88% content CO 2 and 1020C combustion T set.

Elemental Composition Changes in the Potting Soil
On week 1 and week 10, the first and last weeks of the experiment, composite samples were taken from each of the 3 replicates for each treatment (Control, Beer Mash and Food Waste), homogenized, sieved (2 µm) and dried for analysis by XRF (Xray fluorescence) to determine the elemental composition of the soils and whether it changed over time with treatment.

Statistical Analyses and Diversity Measurements
Data from the Dissolved Organic Carbon (DOC) changes in the nutrient solution at different time points (0, 0.1, 24, and 48 h) after digestates amendment and data from carbon, nitrogen, and C:N ratio changes in the potting soil were analyzed by one-way analysis of variance (ANOVA) followed by the all Pairwise Multiple Comparison Procedures (Holm-Sidak method), at P = 0.05 and at P = 0.01, respectively. Data from each elemental composition changes (µg/g) between week 1 and 10 in soil samples from different treatments were analyzed by ANOVA tests followed by multiple comparisons vs. control group procedures (Holm-Sidak method), at P=0.01. Cumulative data from NGS analysis on relative abundance of fluorescent pseudomonads in the nutrient solution samples taken from 4 weeks out of 9 weeks experiments were analyzed by means of Kruskal-Wallis test at p < 0.05. Data were analyzed using SigmaPlot version 14.0, from Systat Software, Inc., San Jose, CA, USA. Measurements included 3 technical replicates. All experiments were repeated at least 3 times, and results from a single representative experiment were chosen for presentation. The bioinformatics and statistical procedures for all Illumina MiSeq data are described within the "library construction, sequence processing, and bioinformatics workflow" section above.

Effects of Mixed Food Waste and Beer Mash Digestates on Culturable Bacterial Populations in the Nutrient Solution
The culturable bacterial populations were monitored in the nutrient solution before and after the amendment. Before the addition of FW and BM, the recirculating nutrient solution used to irrigate citrus plants harbored ca. 5 log CFU/ml of aerobic heterotrophic bacteria (AHB; Figure 2A) and ca. 4 log CFU/ml level of fluorescent pseudomonads ( Figure 2B). However, within 24 h after the addition of FW and BM, the average population of fluorescent pseudomonads in amended treatments (ca. 6-7 log CFU/ml) were two to three orders of magnitude greater than that in non-amended treatments (ca. 3 logs CFU/ml; Figure 2B). The population surge occurred after each consecutive addition (9 weeks) of FW and BM to the nutrient solution. Within 24 to 48 h after the addition of FW and BM, the average fluorescent Pseudomonads population accounted for ca. 97.51% (ranging from 13.74 to 88.13%) and 44.84% (ranging from 2.35 to 94.57%) of the total AHB population, respectively, in the FW and BM -amended solutions but only ca. 39.73% (ranging from 1.34 to 86.63%) in the non-amended nutrient solution.

Taxonomy-Based Bacterial Diversity in the Nutrient Solution and Bacterial Identification
The microbiome in the nutrient solution used in a zero-runoff production system was characterized from samples collected from FW and BM-treated and untreated nutrient solutions from 2 experiments. Our results from the culture-independent NGS method confirmed the results obtained by the culturedependent methods reported above. The bacterial composition of the nutrient solutions ( Figure 3A, FW and BM-treated and untreated) was distributed amongst 11 defined phyla with Proteobacteria being the most abundant. Barplots depict the top 15 defined genera. For data visualization, all other genera were grouped into "other" and colored gray ( Figure 3A). Overall, the relative abundance of the bacterial communities in the 24 h and 48 h untreated control samples was not significantly different compared with the 0 h samples. The largest number of bacterial sequences identified in the untreated nutrient solutions belonged to the genera Parachiamydia (19.4%), Massillia (17.3%), and Thiomicrospira (9.58%), with no significant differences in their relative abundance between 0, 24, and 48 h. Figure 3A shows that while the genera Bdellovibrio, Parachiamydia, Massillia, and Thiomicrospira were the most dominant in all 3 treatments (control, FW and BM) at time 0 h, their relative abundance decreased to undetectable levels through 24 and 48 h following the amendment of the nutrient solution with FW and BM, with the exception of Bdellovibrio, which increased again to 22.9 and 10.2%, respectively, in both FW and BM at 48 h. At the same time, the largest number of bacterial sequences identified in the FW and BM treated nutrient solutions at 24 and 48 h belonged to the genus Pseudomonas. In control samples and at time 0 h, Pseudomonas never exceeded more than 4% of the total relative abundance. In BM samples at 24 and 48 h, Pseudomonas represented 66.83 and 64.91%, respectively. While in FW, Pseudomonas represented 93.79 and 62.89% of the relative abundances for 24 and 48 h, respectively ( Figure 3A). Principal Coordinates Analysis (PCoA) of Bray-Curtis distances revealed that nutrient solution samples formed two distinct clusters according to sample time collection and amendment, with samples collected at time 24 and 48 h following FW and BM amendments forming one cluster and samples collected at time 0 h with all treatments in the other (P < 0.01 [Adonis]) ( Figure 3B).
Because members of the genus Pseudomonas dominated the bacterial community at 24 and 48 h, we analyzed their relative abundance within this genus and the taxa within the top 10 most abundant bacterial species (Figure 2C and  Figure 2C and Supplementary Figure 3). At the highest resolution, two strains closely related to Pseudomonas putida (strain JB, DLL-E4 and 15420352) were found to be the most abundant taxa, with a relative abundance of around 39% in both treatments, within the pseudomonas species identified among 103 ASVs present in both BM and FW Frontiers in Sustainable Food Systems | www.frontiersin.org FIGURE 2 | of 9 weeks experiments. Cumulative data from NGS were analyzed by Kruskall-Wallis tests to verify statistical differences among treatment groups (See Supplementary Figure 3). Head arrows indicate the timing of FW and BM applications (4 ml/L). Error bars represent the SE of the mean (n = 4) where variation was great enough to be presented. treated nutrient solutions. A Bdellovibrio sp. was the 2nd highest abundant taxa with a relative abundance of around 14% in both treatments ( Table 3). Bacterial species identified by 16S rRNA gene Sanger sequencing and a culturing approach confirmed that Pseudomonas sp. (strain HAL1, Contig 19) and Pseudomonas moorei (strain 28, Contig 14) as the most abundant isolates in both BM and FW treated nutrient solutions ( Table 2). When the taxonomic affiliation of the same 2 contigs (#14 and 19) was investigated, and the referenced whole genomes were analyzed against the PATRIC database, the closest strain to both contigs in both BM and FW treated nutrient solutions was again Pseudomonas putida (strain JB), which shared 99.0 and 100% gene sequence similarity, respectively (Supplementary Table 2).

Effect of Waste By-products on Citrus Leaves, Roots, Rhizosphere Soil Bacterial Populations
Samples of citrus leaves, roots, and soil around the roots were taken on week 11 from all 9 replicate zero-runoff production systems. Samples were taken 2 weeks after the last amendment of FW and BM into the system. The bacterial composition on citrus leaves, roots, and soil at the end of the greenhouse experiment was distributed amongst 18 defined phyla. For graphical representations, the top 15 genera from each tissue type are highlighted and the remaining genera are grouped in "other" and colored in gray (Figure 4). Overall, the relative abundance of the bacterial taxa in the untreated control samples were similar to the FW and BM-treated samples. The largest number of bacterial sequences identified in the leaf samples were Bacillus, comprising 17.2, 16.4, and 9.4% of the relative abundance of Control, BM, and FW samples, respectively. Gemmatimonas (16.6%), Burkholderia (0.09%), and Paraburkholderia (7.68%) were the most abundant genera in the untreated root samples.
Gemmatimonas abundance was enhanced in both BM (30.7%) and FW (26.0%) treated root samples. In root samples, the genus Thauera was found in much higher abundance in FW (14%) samples than BM or control (< 1%). Pseudomonas was found in control (5.34%), BM (5.23%), and FW (4.28%) root samples but not in an appreciable amount in soil (< 1%). Gemmatimonas abundance was also enhanced in both BM and FW treated soil samples, accounting for 44.9 and 29.6%, respectively, compared to the control (19.6%). In addition to Gemmatimonas, soil samples treated with BM (15.9%) also had a higher abundance of Brevundimonas compared to the FW (2.5%) and control (0.1%) samples. The samples' alpha diversity was analyzed using Shannon's diversity index (Figures 5A,B). ASV diversity was highly dependent on plant compartment (P < 0.01 [Kruskal-Wallis]) with leaf having a significantly lower Shannon index than root and soil (P < 0.05 [pairwise Wilcox]). High diversity values for soil (4.15 ± 0.198) and root (3.82 ± 0.0977) and consistently decreased diversity estimates in the leaf samples (3.04 ± 0.120) were observed. Shannon indices of the leaf control samples (2.94 ± 0.219) were comparable with those of the leaves of FWtreated and BM-treated samples (2.95 ± 0.277; 3.20 ± 0.105), with no significant difference among the treatments. While at the soil and root level there was a clear separation between FWtreated (4.44 ± 0.197; 3.59 ± 0.144) and BM-treated (3.38 ± 0.390; 3.71 ± 0.150) and untreated samples (4.77 ± 0.145; 4.26 ± 0.163) with the soil control having the highest diversity and being significantly higher than the BM soil (P < 0.05 [pairwise Wilcox]). Control roots had a significantly higher Shannon index than the FW-treated roots (P < 0.05 [pairwise Wilcox]).
Beta diversity analyses revealed that leaves, root and rhizosphere soil samples formed three distinct clusters that were significantly different (P < 0.01 [Adonis]) (Figures 5C,D). Both FW and BM-treated samples separated along the axis 1 with the soil samples separating on the axis 2 as well (P < 0.05 [Adonis]).

Carbon Dynamics in the Nutrient Solution
In the closed zero-runoff systems used in these experiments, there are only 4 main sources of organic carbon compounds in the recirculating nutrient solution: incoming water, plant roots, microbes, and organic growing medium (Waechter-Kristensen et al., 1999). Temporal trends in dissolved organic carbon concentrations are mirrored by changes in bacterial populations where DOC concentrations spike immediately after amendment followed by rapid consumption by microbes within the solution. After 24 h, DOC concentrations return to pre-treatment levels for both FW and BM treatments, demonstrating rapid and complete consumption of added organic carbon. DOC concentrations remained stable throughout the experiment in controls without amendments. These results demonstrate that dissolved organic carbon availability dictates bacterial population growth in the presence of both amendments (Figure 6).
In addition to total concentration of DOC having a strong influence on microbial population dynamics within the nutrient solution, the speciation of available organic carbon compounds provided by the two amendments may provide insight on selectivity of microbes for specific carbon sources. A subset of organic carbon compounds within the nutrient solution was characterized using GC/EI-MS after derivatization, where compounds can then be categorized into general classes (i.e., lipids, condensed hydrocarbons, proteins, lignin, amino sugars, and carbohydrates) based on their O/C and H/C composition (Figure 7). These compounds may have originated from the formation of microbial metabolites, from the potting soil, or the amendment.
The Van Krevelan diagram (Minor et al., 2014;Rivas-Ubach et al., 2018) of beer mash and food waste digestates prior to addition to nutrient solution shows that organic compounds are dominated by condensed hydrocarbons and   (Figure 7). Carbohydrates were only detected in controls, which is unexpected given that both amendments are expected to contain carbohydrate compounds (Popa et al., 2017). Lipids were detected in nutrient solution after BM amendment and were detected up to 48 h after amendment addition. Though lipids were also detected in nutrient solution after FW treatment, they were not detected after 24 h. Overall, organic compound composition generally became dominated by more oxidized compounds (increase in O:C ratio) and more saturated state (i.e., increase in H:C ratio over time).

Nutrient Dynamic in the Potting Soil
Nutrient dynamics were monitored in the potting soil overtime after amendment. First, the soil solid phase C and N content were measured by flash combustion (Figure 8). Both Food Waste and Beer Mash treatments significantly increased soil C and N concentrations as compared to the control (no treatment), with respective p-values of 0.026 (for N) and 0.005 (for C) through analysis of variance and pairwise comparisons between treatments. Beer Mash amendment increased average N by 166% and average C by 164% while the Food Waste amendment increased average N concentrations by 152% and C by 259% (Figures 8A,B). Overall, the treatments BM and FW did not vary significantly from one another with regards to C and N contributions. Average C:N ratios of all soils were near or within the range of 25-35, which is considered optimal for microbial use . Soil amendments with C:N ratio outside of this range require supplementation to prevent N loss (Ibid) ( Figure 8C). Figure 9 shows the changes in elemental composition of soil solid phase collected at Week 1 and 10. We categorized elements into two groups: plant nutrients (Mo, Mg, K, P, Se, Ca, Cl, Fe, Mn, S, and Zn) and potentially toxic elements (Na, Al, Co, Cu, Hg, Pb, V, Cr, As, Cd). While there were changes in the elemental composition with BM and FW treatment, the change in concentrations were not statistically significant for most of the elements analyzed. There were some outliers, notably magnesium (Mg), an essential plant nutrient, which increased in concentration significantly by 1,650% through treatment with FW and 3,242% through treatment with BM (p = 0.003). Phosphorus, Ca, and Co also increased significantly, but only with BM treatment (p = 0.008, 0.004, and 0.028, respectively), while K concentrations decreased significantly (p = 0.042). Overall, amending soils with BM lead to the highest increase in the majority of plant nutrient element concentrations.
Because some elements are toxic to plants and humans and thus undesirable in agricultural amendments, we then analyzed potentially toxic elements (PTEs) separately. Although there was a change in soil concentrations of the measured PTEs over time, there was no statistically significant increase or decrease from Beer Mash and Food Waste treatment as compared to untreated soils for most PTEs, except for Co concentration which decreased significantly (p = 0.028) with BM treatment (Figure 9).

DISCUSSION
The present study demonstrated the effect of digestates, derived from FW and BM, on the structure and diversity of bacterial populations in the microbiome of the recirculating nutrient solution, potting soil, roots, and leaves of non-bearing citrus plants. Data obtained using different approaches (culturedependent and -independent) were combined to gather insights into the distribution of specific taxa and on the ability of strains belonging to these taxa to promote PGPR beneficial traits. Additional analyses of carbon and nutrients dynamics provided preliminary understandings of how changes in the microbial composition, resulting from digestate application, influenced ecosystem functions and below-ground processes such as carbon and nutrient cycling. The implications of such effects and primarily the enhancement of native microorganisms in intensive indoor and potentially field production systems, as well as the safety of digestates, are also discussed here.
Our results from culture-dependent analyses indicated that, for nine consecutive weeks, both types of digestates selectively and significantly enhanced the native population of aerobic heterotrophic bacteria and fluorescent pseudomonads in the recirculating nutrient solution used to fertigate indoor-grown citrus plants (Figures 2A,B). These results agree with earlier findings on the effect of exotic substrates in selecting and stimulating a specific microbial taxon (i.e., Pseudomonas spp.) in the irrigation water that was used in non-woody plants grown in indoor closed production systems (Pagliaccia et al., 2007(Pagliaccia et al., , 2008. Other studies have also shown that nutritional amendments such as salicylate, ammonium sulfate, methionine, and olive oil, when added to the soil or recirculating nutrient solutions, enhanced the efficacy of several bacterial and fungal biological control agents (Colbert et al., 1993a,b;Yamada and Ogiso, 1997;Ji and Wilson, 2003;Luziatelli et al., 2019).
Soils and plant systems represent biologically diverse habitats. Their biota supports many ecosystem functions and they respond to natural perturbations and human interventions (Bender et al., 2016). Our results support the Bender et al. (2016) concept of soil ecological engineering for agricultural sustainability, where targeted manipulation of soil organisms stimulates soil life and harnesses internally regulated ecosystem processes (Philippot et al., 2013;Bender et al., 2016). Targeting and enhancing beneficial native rhizosphere microbial communities combined with currently available state-of-the-art farming technologies FIGURE 4 | Stacked column bar graph representing the relative abundance and distribution, at the end of the greenhouse experiment, of the 15 most abundant genera from each tissue type, leaf, root and potting soil, treated with food waste and beer mash digestate or not treated (control). Samples were taken 2 weeks after the last FW and BM digestate application into the system. FIGURE 5 | Alpha diversity estimates of the bacterial communities in leaves, root and potting soil samples, for experiment 1 (A) and 2 (B), were estimated using the Shannon diversity index. Box plots display the first (25%) and third (75%) quartiles, the median and the maximum and minimum observed values within each data set. Alpha diversity estimates represent biological replicates for the soil (N = 8; 17) root (N = 16; 17) and leaf samples (N = 8; 17). Data were analyzed by means of (Continued) FIGURE 5 | Kruskall-Wallis and pairwise Wilcoxon tests were run to verify statistical differences among groups. The overall plant compartment and wastes digestates treatment effects F (DFn, DFd) and P-value are displayed at the top of each graph. Significant differences (P < 0.05) across plant compartments are indicated with lowercase letters. Bray-Curtis beta diversity for experiment 1 (C) and experiment 2 (D). Points represent individual sample communities and are colored by treatment and shapes represent bio-compartment.
FIGURE 6 | Dissolved Organic Carbon (DOC) changes in the nutrient solution. Average DOC concentration vs. time detected by DOC analyzer in a greenhouse experiment with citrus as the host plant at week 4. We chose to report week 4 because it is representative and because it was a middle point in the experiment. Treatments consisted of nutrient solution receiving no digestate (Control) or amended with beer mash BM or food waste FW digestate at time zero. Time points (0, 0.1, 24, and 48 h) with different letters are significantly different at P = 0.01 according to the all pairwise multiple comparison procedures (Holm-Sidak method). Error bars represent the SE of the mean (n = 3) where variation was great enough to be presented.
(e.g., closed production systems and precision agriculture with cloud-based applications for optimized plant growth parameters) could help maintain yields while reducing the dependence on external inputs (e.g., chemical fertilizers and pesticides) in agrosystems (Delgado et al., 2019). Using a variety of experimental approaches, including cultureindependent analyses, we showed that the relative abundances of Pseudomonas spp. and to a lesser extent Bdellovibrio spp., increased in the irrigation water treated with both digestates, while the bacterial diversity decreased within 24-48 h compared to the non-treated systems. The low microbial diversity and high Pseudomonas spp. abundance in the nutrient solution persisted for 48 h, whereas no significant difference in the abundance of the Pseudomonads was observed 7 days after treatment (Figures 2C, 3A, and Supplementary Figure 3). These results and the ecological diversity measures ( Figure 3B) suggest that, once the nutrient sources are depleted, the bacterial community structures reset to a more diverse composition. In our study, Bdellovibrio was the only genus that survived overtime in the nutrient solution (e.g., longer than 72 h) and it took over when Pseudomonas spp. abundance decreased (Figure 3A). This is not surprising since Bdellovibrio has been shown to be an effective predator of Gram-negative bacteria; it also has biological control properties in environmental as well as medical microbiological settings (Harini et al., 2013) however, further studies are needed to evaluate its potential as PGPMs for agricultural use. Wu et al. (2018) and Wright et al. (2019) also observed a similar effect on the reduction of the bacterial biodiversity, studying the dynamics of the bacterial community using dissolved organic matter (DOM) extracted from sediments adjacent to groundwater and chitin, respectively. In both studies, bacteria able to use DOM and chitin as nutrient sources were able to proliferate and became FIGURE 7 | Gas-chromatography/electron ionization mass-spectrometry (GC/EI-MS) analysis of (A) beer mash (BM) and food waste (FW) digestates amendments and (B) nutrient solutions sampled from controls (Not treated), BM, and FW treatments. Samples were taken before and after application (-BA and -AA), 24 and 48 h after digestate application. O/C vs. H/C ratio of elemental formulas was plotted to provide an overview of the compound categories within each amendment. Data were assigned nominally into categories based on a van Krevelen empirically derived metabolite map depicting chemical compound classes and their relative boundaries as assigned in Figure 5 in Kim et al. (2003) study. FW-AA, FW-BA, and FW2 samples had values below detection-limit and are not reported here.
dominant. In agreement with these studies, we also identified bacteria belonging to Pseudomonas spp. that were enhanced by the organic components of the digestates.
The carbon dynamics assays supported the findings of the microbial analyses. DOC in the recirculating nutrient solution followed a similar trend as the bacterial populations, spiking FIGURE 8 | runoff plant production systems set up in a greenhouse where nutrient solution was amended and not amended (control) with beer mash BM and food waste FW digestate at 7 days interval for a total of 9 treatments applications. For (A,B) values are the means of three replicates. Bars with different letters are significantly different at P = 0.05 according to the all pairwise multiple comparison procedures (Holm-Sidak method). Error bars represent the SE of the mean (n = 3) where variation was great enough to be presented. For (C) data were analyzed by two-way analysis of variance. No statistically significant difference was found within treatments (P = 0.109), within weeks (P = 0.369) and no interaction between treatments and weeks (P = 0.123).
after each digestate amendment, then plateauing again 24 h post-treatment. These results suggest that the bacteria used organic carbon from FW and BM for metabolic processes and replication. In the recirculating nutrient solution of the closed zero-runoff systems, there are four primary sources of organic carbon: influent water, plant roots, microbes, and organic growing medium. The plateau in DOC concentration in the digestate treated nutrient solutions (Figure 6) likely indicates that depletion of labile carbon in the nutrient solution induced a change in microbial activity from growth to maintenance. Unlike the DOC data, changes in carbon speciation over time did not show a temporal trend, though there was a slight trend toward more saturated compounds dominating over time, and the amendments starting out highly reduced (electronrich) (Figure 7). Organic compounds exhibited increasing O:C over time likely due to microbial oxidation during aerobic respiration and some contribution from abiotic oxidation from the nutrient solution being open to the air. A previous study showed a gradual increase toward carbonyl/carboxyl structures attributed to metabolites such as nucleobases, lipids, and amino acids/peptides/protein products of microbial metabolisms (Wu et al., 2018). The authors compared van Krevelen diagrams of detected compounds with and without microbes and found a clear disappearance of protein-like compounds with elemental composition of carbon, hydrogen, oxygen, and nitrogen (CHON), indicating utilization or degradation of these compounds by microbes over the experimental period.
Many studies have demonstrated that soil microbiota drive a wide range of biogeochemical processes, with many species mediating organic matter transformation, which in turn shape microbial community composition over time (Jouraiphy et al., 2005;Kramer and Gleixner, 2008;Gougoulias et al., 2014;Boye et al., 2017). Shifts in microbial community composition have been significantly correlated to changes in the relative proportions of tannin-, protein-, condensed aromatic-, lipid-, and lignin-like compounds in other studies, where relatively labile tannin and protein compounds are correlated with the microbial community in early stages and microbial community composition changes over time to be more correlated with less labile lipid and lignin-like compounds (Wu et al., 2018). Gilmour et al. (2003) suggested that stable microbial decomposition of the organic matter by anaerobic digestion and composting produces a biosolid organic matter residue that is relatively recalcitrant. This supports the results in our study, where the relative FIGURE 9 | Elemental composition changes (µg/g) between week 1 and 10 in potting soil samples from different treatments as detected by X-ray fluorescence spectrometer. Treatment consisted of nutrient solution not amended (Control) in blue, or amended with Beer mash (BM) in red, food waste (FW) in yellow. Plant nutrient elements tested included: molybdenum (Mo), magnesium (Mg), potassium (K), phosphorus (P), selenium (Se), calcium (Ca), chlorine (Cl), iron (Fe), manganese (Mn), sulfur (S), zinc (Zn); Potentially toxic elements tested included: sodium (Na), aluminum (Al), cobalt (Co), copper (Cu), mercury (Hg), lead (Pb), vanadium (V), and chromium (Cr). Arsenic (As) and cadmium (Cd) concentrations were considered too low to include. Error bars represent standard deviations of the mean (n = 3) where variation was great enough to be presented. For each element, significant differences observed by ANOVA tests (at P = 0.05) followed by multiple comparisons vs. control group procedures (Holm-Sidak method) are indicated with asterisks *.
proportions of lignin-like and CHO compounds continuously increased during incubation (Figure 7), likely because of their relative resistance to biodegradation, contribution of microbial metabolites, and/or oxidation of certain DOM moieties (Fuchs et al., 2011).
Both our previous and current work suggest that, after amendments, microorganisms preferably utilize labile C in DOC (i.e., aromatic compounds), as indicated by the rapid increase in cell counts (Pagliaccia et al., 2007(Pagliaccia et al., , 2008. In Wu et al. (2018) study, microbial community structure analysis showed a rapid enrichment of Pseudomonas spp. in early stages. Pseudomonas spp. was often reported to be the dominant genus in microcosms amended with labile C, such as acetate (Cui et al., 2011) or glucose (Ghosh and Leff, 2013). Also, many Pseudomonas species are known to be efficient competitors for resources through secretion of antibiotics (Hibbing et al., 2010) and siderophorebased metal chelators (Stubbendieck and Straight, 2016). It is therefore not surprising that Pseudomonas became the most abundant species in the nutrient solution in the early incubation stage when labile C was still available. In our study, conventional identification (culturing plus Sanger sequencing) and Illumina MiSeq sequencing showed that both BM and FW applications changed the structure and composition of the Pseudomonas community in the nutrient solution.
Culture-dependent analysis identified Pseudomonas sp. strain HAL1 as the most abundant bacteria in the nutrient solution samples treated with either BM or FW after 24 and 48 h. P. moorei was also enriched after digestates amendment ( Table 2). Culture-independent analysis identified Pseudomonas putida (strain JB, DLL-E4, and 15420352) as the most abundant bacteria in both BM and FW treated nutrient solutions while no Pseudomonas spp. were found within the most abundant bacteria in the control nutrient solutions ( Table 3). These results are in agreement with current literature on phytoremediation principles, which indicates that bacteria with PGPR traits can facilitate the removal of inorganic contaminants stimulating either plant growth or phytoremediation activity (Glick, 2010;Mahanty et al., 2017;Luziatelli et al., 2019). This principle could be used to ensure the selection of a pseudomonas population from the rhizosphere and the nutrient solution, for beneficial purposes such as biocontrol, biofertilization, phytostimulation, and disease management in both indoor nursery and citrus field production.
Plant Growth-Promoting Microbes have a direct effect on phytopathogens through the production of metabolites with antibiotic activity, or by competition for space and nutrients; and an indirect beneficial effect on the plant by enhancing defense responses or promoting plant growth (Mendes et al., 2011;Ma et al., 2019). Our preliminary bioinformatics analysis (PATRIC) indicated that all the bacteria species enriched by BM and FW amendments possess genes involved in the production of siderophores and iron acquisition (Supplementary Figure 1A, e.g., Ferric Siderophore Transport system-ABC efflux system and TonB family proteins). Some of these bacteria also possess genes putatively involved in direct antagonism, phosphate solubilization, and nutrient facilitation (Supplementary Figure 1B, e.g., Acetolactate synthase, Transcriptional regulator PadR). Their enrichment in the recirculating nutrient solutions could be very valuable for the plant since as shown in this work, the fluorescent Pseudomonas spp. reached a consistent population level in the nutrient solution of 10 7 -10 8 CFU/ml after each digestate application. This level of bacterial cell concentration is wellknown to induce the activation of quorum sensing regulatory systems (Hibbing et al., 2010). The regulation of genes through quorum sensing controls essential biological functions in many bacteria, including several regulatory mechanisms for the production of beneficial secondary metabolites in Pseudomonas spp. (i.e., antibiotics and siderophores) (Sahu et al., 2018). This finding is particularly promising for the development of new approaches for integrated plant production and disease management practices, especially under zero runoff indoor production systems, where the beneficial bacterial metabolites produced in the nutrient solution would reach the plant roots at every irrigation cycle (i.e., once to several times a day). This approach could also be applied under field conditions where targeted management of soil community composition could be promoted around the plant root system with the application of digestates through drip fertigation systems. Barzee et al. (2019) investigated the feasibility of drip fertigation for the application of digestate biofertilizers (rich in N, P, K as well as Cu, Ni, Mg, and Zn) to tomato crops in farm-scale cultivation experiments. Their results demonstrated that digestate biofertilizers could also be effectively delivered to the tomato plants through a subsurface drip irrigation system, with increased crop yields.
Digestate materials could have potential agronomic value because they contain high amounts of nutrients and organic matter (Möller and Müller, 2012), as well as volatile organic acids and complex organics from undigested particles (Franke-Whittle et al., 2014). However, the literature is divergent about the effect of digestates on soil quality, and the fate of the digestates' C and N in soils is still an understudied topic (Möller, 2015). We carried out a preliminary analysis to elucidate the agronomic value of BM and FW and their potential impact on long-term sustainability in terms of soil fertility and environmental impact. From a review of the available literature on this subject, Möller (2015) stated that due to the lower portion of easily degradable compounds in digestate, the organic fractions in the digestate can contribute to soil organic matter (SOM) stabilization (i.e., increase stable soil organic carbon), which can lead to changes in biological, chemical and physical soil characteristics. In our study, both BM and FW significantly increased N and C in the potting soil around the roots (Figures 8A,B). Regarding the C/N ratio, a range of 25 to 30 is considered optimal for microbial activity, balanced composting, and reduced nitrogen loss . In our case, no statistically significant difference was found in the C/N ratio of samples analyzed between treatments and between the first and last time points; however, the C/N ratio did increase from 16 to 25 in soils treated with food waste digestates. These results suggest that in the short term the high N content in the BM and FW digestates can deliver immediate nutrients to citrus trees in the field, while in the long term, the high C content could provide a slower release of nutrients (i.e., through PGPR mediated processes) and help reduce nutrient leaching. Furthermore, the bokashi fermentation used to produce both BM and FW digestates significantly reduces gaseous emission with a GHG footprint of < 0.22 MtCO 2 equivalent (million metric tons) and 1% of CO 2 released during fermentation (Alattar et al., 2012) compared to standard aerobic composting which has a GHG footprint of ca. 1.2 MtCO 2 equivalent (million metric tons) and ca. 50% carbon released as CO 2 . Overall, our results are in agreement with Möller (2015) and support the hypothesis that although the digestate's potential as fertilizers is of minor relevance, their indirect effect on the microbial community, the stabilization of the organic matter, and indirect reduction of gaseous emissions have the greatest relevance. Appropriate application rate and time need to be investigated further to optimize the long-term effects of the digestates and to stabilize the positive microbial activities (i.e., keep the bacteria population-level above 10 5 CFU/ml). Further studies are also needed to investigate which types of feedstock would produce the most valuable digestates in terms of agronomic value.
Based on previous findings (Pagliaccia et al., 2007), and the recurring re-inoculation of the rhizosphere through the recirculating irrigation cycle, we would have expected the Pseudomonas populations enriched in the nutrient solution to colonize the citrus roots system. To our surprise, the bacterial analyses of citrus leaves, roots, and potting soil produced different patterns at the end of the greenhouse experiment. NGS sequencing showed that, compared with non-treated roots, roots treated with BM and FW had higher relative abundances of Gemmatimonas spp. or Brevundimonas, yet no significant difference in the abundance of Pseudomonas spp. was observed. Additionally, we found greater variation in beta-diversity for both root and potting soil microbiomes, while all leaf samples were very similar, but with the lowest diversity and evenness estimates within all treatments (Figure 5). This contrasting pattern was surprising but could reflect divergent microbial responses in different environments (nutrient solution vs. soil) but also within different plant compartments (root vs. leaves) and plant species. In their extensive review, Philippot et al. (2013) discussed how soil properties but also plant species influence microbial community composition and structure in the rhizosphere. Further work with other citrus species and varieties as well as other crops will be essential to decipher the mechanisms underlying host plant effects on the root microbiome, as well as any fitness benefits to the plant caused by the increased Pseudomonas abundance in the irrigation water. On the other hand, our results raise the intriguing hypothesis that the citrus genotype might have mediated changes in the host-microbiome under our abiotic induced conditions (i.e., digestates application). In our study, the citrus cultivar used ("Pineapple sweet orange, " Citrus sinensis (L.) Osbeck) seemed to have lost the capability to recruit the Pseudomonas onto the root system (i.e., Pseudomonads high population of 6-7 log CFU/ml from irrigation water did not colonize root). The "Pineapple sweet orange" has been continuously propagated and cultivated in commercial USA settings since the 1800's (Hodgson, 1967). It has been shown that, due to cultivation practices and intense artificial selection for yield and other desirable agronomic characteristics, quite a few plant cultivars have lost, to some extent, their ability to recruit and associate with beneficial soil biota (Sawers et al., 2008;Porter and Sachs, 2020). Mueller and Sachs (2015), discussed the native microbiome and anticipated that hostmediated microbiome engineering would be more efficient using wild hosts rather than hosts that have experienced domestication, or adaptation to microbially depauperated experimental environments.
In our study, we also tested changes in toxic vs. nontoxic micronutrients in the potting soil following BM and FW application, and found negligible signs of accumulation for both, except for Mg, which increased significantly with both digestates compared to the non-treated system. P and Ca also accumulated, but only with BM treatment while K decreased (Figure 9). Some of the results were unexpected since the FW digestate amendment contained a higher amount of K and Ca ( Table 1). Although the FW digestate contained a much higher amount of Na (Table 1), which is expected due to the higher content of Na typically present as a food preservative in food waste (Popa et al., 2017), our analysis showed a decrease in this element, which is a very desirable trait in agriculture and especially in zero-run off plant production systems. These results agree with Barzee et al. (2019), who showed that digestates do not cause a higher risk for salt accumulation. However, as in Barzee's study, longterm trials with drip irrigation systems need to be carried out to understand the impact of repeated digestates applications on the soil chemical properties in both greenhouse and in the field and, most importantly, their effect on the rhizospheric soil and citrus plant biota, especially the roots.
The validation of the digestates' agronomic values also demand that they do not contain or enhance human pathogens and therefore pose a threat to workers and consumers' health. Our results suggest that, although DNA of bacterial groups belonging to the genus Enterococcus, known to be human pathogens, were detected by NGS analysis in the digestates stocks (Supplementary Figure 2), such bacteria were never cultured or detected from the digestates stocks or ever identified by NGS analysis in any bio compartment tested (i.e., soil, root, leaves) in the plant production system after digestates application (Table 2, Figures 3A, 4). Testing for fecal coliform and Salmonella (A & L Western Agricultural Laboratories) confirmed the absence of known human pathogens (data not shown). Therefore, we concluded that under our experimental conditions, the use of Bokashi fermentation digestates is safe. However, additional physical or chemical processes might be necessary to treat digestates to ensure their safe use in diverse non-experimental agricultural systems. Additional investigations are also necessary to assess cause and effect, such as Koch's postulates experimentation, to further define Pseudomonas spp.'s roles in controlling plant diseases in planta. Regarding the potential use of these bioproducts to increase or accelerate citrus growth, the overall results indicate that their application could be beneficial for the citrus industry, inside and outside California, but could also be translated to other crops. However, because citrus plants are slow-growing woody plants, the current tests did not provide any evident biostimulant effect. Longer-term trials with additional citrus varieties will be required to answer this question. Alternatively, further tests with faster-growing plants could be used to evaluate the effect of the digestates on plant growth rate, fitness, etc. A fascinating question to answer is whether, within the same plant species, different crop varieties (new commercial variety vs. a wild relative) will be colonized differently by the digestate-enhanced Pseudomonas spp.
Based on our results, digestates derived from agriculture and food waste can be safely used at very low concentrations as biostimulants and offer a promising solution to enhance the native ecosystem while reducing environmental impact. Further research is needed to determine if such applications would have any impact on plant performance (e.g., yield increase or disease resistance) and to identify the underlying mechanisms of such effects.

CONCLUSIONS
More sustainable agri-food systems can be achieved by adopting biotechnological processes that allow nutrient recycling, and the implementation of climate-smart agricultural (CSA) approaches. This study addresses the pressing need to have easily deployable mitigation strategies that integrate the management of soil carbon dynamics, favor beneficial microbiota, and create alternative disposal for agricultural and food waste.
We showed that there is no need to introduce PGPMs into a cropping system (rhizosphere and nutrient solution). Developing novel agricultural practices that deploy targeted manipulation of soil organisms to stimulate soil life and harness internally regulated ecosystem processes offer a more sustainable approach to enhance crop production. Farmers can benefit from the resident microbiota because the plant rhizosphere already contains prospective members that can function as plant growth-promoting (PGP) and as plant disease controlling (PDC) agents when enhanced by digestate by-products obtained from agricultural and food wastes. Digestates can be used as biostimulants to generate a positive shift in the native microbial communities with a parallel increase in beneficial microbiota (i.e., fluorescent Pseudomonads) and their related bioprocesses. This could offer a very promising alternative to synthetic chemicals to stimulate plant growth traits and increase plant defense responses (priming) to pathogens.
Combined with their widespread sourcing from agricultural and industrial processes, waste digestates could help reduce dependency on synthetic chemicals and foster the establishment of eco-friendly, circular economies for agricultural production. More research is needed to elucidate the best feedstock to use for digestate production and the many indirect effects digestates have on plant pathogens, crop yield, N emissions in indoor and in the field production, and the possible long-term effects on the microbiome and soil fertility.

DATA AVAILABILITY STATEMENT
The original contributions generated for the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

AUTHOR CONTRIBUTIONS
DP wrote the manuscript, was primarily responsible for conceptualization and experimental design, contributed to data collection, analysis, and all manuscript drafts. SB and XC were responsible for greenhouse operation, carried out the data collection and analysis for all the microbial work, and assisted with manuscript revisions and tables creation. DS contributed to data collection and analysis for the carbon dynamic work, wrote with SY some manuscript sections for the carbon dynamic work, performed manuscript revisions, and assisted to create some figures and tables. ED performed the sequence processing, bioinformatics, created the figures on the NGS data, and contributed to the manuscript revisions. AD contributed extensively to the manuscript revisions, figures, and tables creation. SY was responsible for the conceptualization and experimental design and collection for the carbon dynamic work, revised, and corrected the manuscript. BP assisted SB and XC for the NGS Library construction. JB, PR, BP, NE, PR, RP, and GV assisted with the experimental design, revised and corrected the manuscript. JB, SY, and GV provided greenhouse and laboratory resources for this study. All authors contributed to the article and approved the submitted version.