Experimental setup and sampling
The study was conducted in a decoupled aquaponic system where four hydroponic units with plants were combined with one RAS unit with rainbow trout (Oncorhynchus mykiss) and operated for 30 days. The whole experiment was repeated three times to serve as replicates in time. The hydroponic system consisted of four deep water culture (DWC) units (1 m × 1 m × 0.35 m) with a polystyrene floating cover with 30-mm holes for holding plant baskets. Before the aquaponics experiment, lettuce and companion plant seeds were germinated in a greenhouse (for details see Supplementary information). After the germination period, 15 seedlings of lettuce with 15 seedlings of companion plants were transplanted in each DWC. When lettuce was the companion plant as a control treatment, 15 + 15 lettuce seedlings were transplanted. Thus, in one DWC unit, we had one of the following treatments: (i) lettuce + companion plant spearmint (Mentha spicata) Lett-MINT, (ii) lettuce + companion plant wormwood (Artemisia absinthium) Lett-WORM, (iii) lettuce + companion plant rucola (Diplotaxis tenuifolia) Lett-RUCO, (iv) lettuce + companion plant lettuce (Lactuca sativa) Lett-LETT (control treatment).
Sixteen hours of light was provided to the plants during the whole experiment by using LED lights (King wua Bright®, China). Humidity varied according to the temperature (Fig. S1) but was over 50% during the whole experiment. The four DWC units were completely separate and did not interact with each other. At the start of the experiment, 135 L water from the RAS unit was given to each DWC and later 100 L water was given daily. Water was aerated in the fish tank, biofilter tank, sump tank, and each DWC unit. Each DWC was supplied with additional trace elements (1.25 ml micronutrient solution (Ingestad 2006) and macronutrients (2 g of K2SO4 and 1.6 g of MgSO4⋅ 7 H2O) daily.
The RAS unit consisted of a fish tank (1 m3), swirl pool solid filter (500 L), fixed bed biofilter (500 L), and sump tank (500 L). The total volume of water in the RAS was approximately 2500 L. In the RAS unit, 30 rainbow trouts (Oncorhynchus mykiss) each of approximately 200–250 g in weight were stocked at the start of each time replicate. Fish were purchased from a commercial supplier in Finland. The fish were not used in any procedures and they only provided water to grow plants in the hydroponic systems. In the RAS water, the concentrations of total ammonia nitrogen (TAN), nitrite, nitrate, and pH were recorded twice a week while temperature was recorded daily (Fig. S2). Tap water was added to the RAS unit to replace the water given to the DWC units. The DWC and related tubing were disinfected with 90% ethanol at the start of each time replicate. The position of the treatments in the DWC units was randomized between the time replicates. The details of the operation of the aquaponic system are provided in Supplementary information.
Plants were harvested at the start and on the 30th day of each time replicate. Start point samples contained the whole seedling, i.e., root and shoot together due to the small size of the plants. For the endpoint samples, roots and shoots were sampled separately. Each plant sampling consisted of three replicates of five pooled plants.
Lettuce and companion plant dry weight (60 °C, 72 h) was recorded at the end of each time replicate. For bacterial community analysis, three samples of five pooled plants from each plant species were surface-washed with sterile water, and samples of approximately 100 mg of fresh weight were collected. RAS water was sampled by taking three 50-ml water samples on the day of the harvest from the sump tank of RAS and filtered through a membrane filter (Millipore, 0.22 µm pore size, Ø 47 mm). All samples were stored at − 80°C.
DNA extraction
DNA for bacterial community analysis was extracted from 100 mg of plant material with the NucleoSpin Plant II kit (Macherey-Nagel, Düren, Germany) by using PL1 lysis buffer. Samples were homogenized using two 3.2 mm chrome steel beads and approximately 0.6 g of 0.1-mm zirkonia-silica beads (Biospec) in an Omni Bead Ruptor Elite connected to an Omni BR-Cryo Cooling Unit with speed 4 m/s for 2 × 45 s under cooling. DNA from the water sample filters was extracted with the Nucleospin Soil kit (Macherey-Nagel, Düren, Germany). DNA was stored at − 80°C.
PCR amplification
Bacterial 16S ribosomal RNA (rRNA) gene PCR was carried out as a nested PCR to limit the co-amplification of plant chloroplasts and mitochondria. The primers in the first PCR step were 799F (5´-AAC MGG ATT AGA TAC CCK G-3´) and 1492R (5 ´-GGY TAC CTT GTT ACG ACT T-3´) (Chelius and Triplett 2001). The 1st PCR step was run with the reaction components 1 × GoTaq buffer, 0.2 mM dNTPs, 0.24 µM primers, 0.75 U GoTaq polymerase (Promega) and 2–10 ng DNA template with 25 µl of total volume per reaction. The cycling program was 95°C for 3 min, followed by 25 cycles of 95°C for 30 s, 53°C for 40 s, 72°C for 60 s and finally 72°C for 5 min.
The primers in the 2nd PCR step were M13-1062F (M13 linker for attaching barcodes and sequencing adapters 5’-TGT AAA ACG ACG GCC AGT-3’ (Mäki et al. 2016) followed by 1062F 5’-GTC AGC TCG TGY YGT GAG-3’) (Ghyselinck et al. 2013)) and 1390R (5´- ACG GGC GGT GTG TRC AA- 3´) (Zheng et al. 1996). In the 2nd PCR step, 2 µl of the 1st PCR product was used as a template with the same reaction components and cycling program as in the 1st PCR with 20 cycles.
Barcodes and sequencing adapters were added in a separate PCR step. The forward primer contained sequencing adapter IonA, barcode and M13 linker. The reverse primer was P1-1390R, where P1 is a sequencing adapter. The other reaction components and the cycling program were the same as in the 1st PCR but with 8 cycles and 2 µl of the 2nd PCR product as template. PCR products were purified with sparQ PureMag beads (Quantabio) and the DNA concentration of purified products was measured with Quant-iT PicoGreen dsDNA Assay (Invitrogen). The proportion of bacterial (ca. 350 bp) and plant mitochondrial (ca. 700 bp) amplicons in the purified PCR products were estimated based on image analysis of agarose gels with ImageJ software, and equimolar amounts of the products were pooled based on the bacterial product. The mitochondrial product was removed by gel extraction of the bacterial amplicon pool with Monarch DNA Gel Extraction Kit (New England BioLabs).
Sequencing and sequence processing
Pooled bacterial amplicons were sequenced using Ion Torrent PGM with Ion PGM Hi-Q View OT2 Kit, PGM Hi-Q View Sequencing Kit, and Ion 316 Chip v2. The sequencing data was analyzed with mothur software package v.1.43.0 (Schloss et al. 2009 following the relevant parts of MiSeq SOP protocol (https://mothur.org/wiki/MiSeq_SOP; accessed in April 2020; Kozich et al. 2013). Sequences were quality filtered using average quality of 25 and a window size of 50 bases, a minimum sequence length of 200 bp, a maximum length of 400 bp, maximum homopolymer length maxhomop = 8, maximum number of ambiguous bases maxambig = 0, maximum number of differences to primer sequence pdiffs = 1, and maximum number of differences to barcode sequence bdiffs = 0. The sequences were aligned against the Silva database v.1.38 (Quast et al. 2013). Chimeras were detected with UCHIME in de novo mode (Edgar et al. 2011) in mothur. After quality filtering, alignment, and removal of chimeras and nontarget sequences, 1485253 bacterial sequence reads were recorded. Pairwise distances of the sequence reads were calculated for the clustering step and operational taxonomic units (OTUs) were defined with a 97% similarity cutoff by the opticlust method. The OTUs were classified against Silva v. 1.38 to get their taxonomic affiliations. Bacterial richness (number of OTUs) and Shannon diversity were calculated in mothur with the command summary.single. One start sample, four companion plant samples, and 8 lettuce samples were removed due to few reads. The sequence data were submitted to NCBI under BioProject accession no. PRJNA814115.
Statistical analyses
Total plant production and lettuce biomass were calculated as the average weight of three composite samples consisting of 5 plants per species by subtracting the start weight from the end weight in each time replicate (n = 9). The total production of lettuce grown with lettuce was calculated the same way using three composite samples consisting of 5 plants per category (5 target lettuce and 5 companion lettuce) in each time replicate.
To assess the differences between treatments in total dry weight of lettuce (end lettuce weight – start lettuce weight) and total plant production in the systems, a two-factor analysis of variance (ANOVA) followed by Tukey’s post hoc test was performed where time replicate was considered as a random factor and treatment (companion plant species) as a fixed factor. In the case the time replicate and the treatment effect interacted, a one-factor ANOVA followed by Tukey’s post hoc test was performed, where data were split into time replicates and treatments were considered as a fixed factor. Root and shoot OTU numbers and Shannon diversity between treatments were compared with one-factor ANOVA followed by Tukey’s post hoc test as above. The normality of variables was checked by analyzing Q-Q plots and histograms. The homogeneity of variance was tested using Levene’s test. All statistical analyses related to biomass, OTU number, and Shannon diversity index were conducted using the software SPSS (version 26).
The community composition of bacteria was assessed at the OTU level. All the microbial community analyses were carried out in R (v. 3.6.1, R Core Team 2019) and RStudio. Plots were prepared with the package ggplot2 (Wickham 2016). The R package phyloseq (McMurdie and Holmes 2013) was used to prepare taxonomy plots. Bacterial OTU data were rarefied with the function rrarefy in the vegan package (v. 2.4.2, Oksanen et al. 2019) to the median of read counts (7366 reads per sample). For the samples with less than median reads, all the reads were included (minimum 559 reads). Only OTUs with more than 5 reads were included in the multivariate analyses. The differences in bacterial communities between treatments were tested with permutational analysis of variance (PERMANOVA, McArdle and Anderson 2001) with the function adonis2 in vegan with Bray-Curtis dissimilarity and 999 permutations for significance testing. Time replicate was used as a blocking factor in permutations. To check if the data meet the assumptions of PERMANOVA, multivariate dispersion was tested by the function betadisper. The microbial community differences between the treatments were visualized by non-metric multidimensional scaling (NMDS) plots with Bray-Curtis dissimilarities with the function metaMDS in vegan. To identify OTUs that were differentially abundant in lettuce grown with companion plants compared to lettuce grown alone, we used DESeq2 (v. 1.28.1, Love et al. 2014), ANCOM (v. 2.1, Mandal et al. 2015; Kaul et al. 2017) and ALDEx2 (v. 1.20.0, Fernandes et al. 2014). An unrarefied OTU table was used, and only OTUs that occurred in 25% of samples in each comparison (lettuce vs. lettuce with companion) were included. An OTU was considered differentially abundant when at least two of the three tools detected it. In DESeq2, the design was ‘time replicate + treatment’, and cutoffs of p < 0.05 and log2 fold change > 2 were used for differentially abundant OTUs. In ANCOM, the analysis included the time replicate as a covariate, and treatment was used as a group indicator for detecting structural zeros. In ALDEx2, we used 256 Monte Carlo instances and effect size > 1 to detect differentially abundant OTUs. Redundancy analysis (RDA) was performed on Hellinger transformed OTU data (rarefied, only OTUs with > 5 reads) in vegan to test for correlation between bacterial community composition and lettuce biomass with the time replicate as a conditional variable. The significance of the RDA model was tested as permutations with the function anova.cca. Biomass-affected OTUs were defined as OTUs with RDA axis 1 scores ± 2 standard deviations of the median score. To estimate the potential sources of microbes in lettuce, we used SourceTracker (Knights et al. 2011). Only OTUs with more than 15 reads in the rarefied OTU table were included. Companion plants, lettuce itself, companion plant and lettuce seedlings, RAS water, and DNA extraction and PCR controls were considered as potential sources of microbes, and lettuce grown together with companion plants as the sink.