J A Ezzedine

• A Janicot

• S Rasconi

• I Domaizon

• S Jacquet

Short-Term Dynamics of Bdellovibrio and Like Organisms in Lake Geneva in Response to a Simulated Climatic Extreme Event

Keywords: Bdellovibrio and like organisms, Dynamics, Prokaryotes, Lake Geneva, Ecological stress

The short time-scale dynamics of three families of Bdellovibrio and like organisms (i.e. Bdellovibrionaceae, Peredibacteraceae, and Bacteriovoracaceae) were studied on the surface waters of Lake Geneva in summer. Using mesocosms deployed nearshore in July 2019, we simulated an extreme climatic event (an input of carbon from the watershed in response to runoff from the catchment, light reduction, and mixing in response to stormy conditions) and aimed to study the impact of both abiotic and biotic factors on their dynamics. The three families of Bdellovibrio and like organisms (BALOs) showed different dynamics during the experiment. Peredibacteraceae was the most abundant group, whereas Bacteriovoracaceae was the least abundant. Compared with the other two families, the abundance of Bdellovibrionaceae did not fluctuate, remaining relatively stable over time. Environmental variables only partially explained the dynamics of these families; in particular, temperature, pH, and chloride concentrations were positively correlated with Bacteriovoracaceae, Bdellovibrionaceae, and Peredibacteraceae abundance, respectively. Prokaryote-like particles (PLPs), such as those with high DNA content (HDNA), were strongly and positively correlated with Peredibacteraceae and Bacteriovoracaceae. In contrast, no relationships were found between Bdellovibrionaceae and PLP abundance, nor between the virus-like particles (VLPs) and the different BALOs. Overall, the experiment revealed that predation was stable in the face of the simulated climatic events. In addition, we observed that Peredibacteraceae and Bacteriovoracaceae share common traits, while Bdellovibrionaceae seems to constitute a distinct category.

Introduction

Among microbial predators, Bdellovibrio and like organisms (BALOs) are obligate predatory bacteria with a broad spectrum of prey that are mainly Gram-negative bacteria. BALOs have been reported to act as "population-balancers" [START_REF] Iebba | Higher prevalence and abundance of Bdellovibrio bacteriovorus in the human gut of healthy subjects[END_REF] and "microbial alpha diversity drivers" [START_REF] Johnke | Bdellovibrio and like organisms are predictors of microbiome diversity in distinct host groups[END_REF], highlighting their potential important roles in microbial ecosystems. They are classified into two distinct polyphyletic clades, Oligoflexia and Alphaproteobacteria, and are characterised by two well-described modes of reproduction [START_REF] Hahn | Silvanigrella aquatica gen. nov., sp. nov., isolated from a freshwater lake, description of Silvanigrellaceae fam. nov. and Silvanigrellales ord. nov., reclassification of the order Bdellovibrionales in the class Oligoflexia, reclassification of the families Bacteriovoracaceae and Halobacteriovoraceae in the new order Bacteriovoracales ord. nov., and reclassification of the family Pseudobacteriovoracaceae in the order Oligoflexales[END_REF][START_REF] Koval | Reclassification of Bacteriovorax marinus as Halobacteriovorax marinus gen. nov., comb. nov. and Bacteriovorax litoralis as Halobacteriovorax litoralis comb. nov.; description of Halobacteriovoraceae fam. nov. in the class Deltaproteobacteria[END_REF][START_REF] Jurkevitch | Phylogenetic diversity and evolution of predatory prokaryotes[END_REF]. BALOs are found in many habitats such as soil and aquatic environments [START_REF] Jurkevitch | Phylogenetic diversity and evolution of predatory prokaryotes[END_REF][START_REF] Williams | Data report: exploring the presence of Bdellovibrio and like organisms in deep-sea sediment by culture-independent and culture-dependent methods[END_REF][START_REF] Sutton | Ecology and characteristics of Bdellovibrios from three tropical marine habitats[END_REF], and they have recently been discovered and studied in peri-alpine lakes [START_REF] Paix | Diversity, dynamics and distribution of Bdellovibrio and like organisms in peri-alpine lakes[END_REF][START_REF] Ezzedine | New 16S rRNA primers to uncover Bdellovibrio and like organisms diversity and abundance[END_REF][START_REF] Ezzedine | Bdellovibrio and like organisms in Lake Geneva: an unseen elephant in the room?[END_REF][START_REF] Ezzedine | Bdellovibrio sp.: an important bacterial predator in Lake Geneva?[END_REF]. In general, the water column is not a favourable environment for the multiplication of BALOs, which thrive better in biofilms, sediments, or closed environments such as in aquaculture [START_REF] Williams | Ecology of the predatory Bdellovibrio and like organisms[END_REF][START_REF] Li | Bdellovibrio and like organisms enhanced growth and survival of Penaeus monodon and altered bacterial community structures in its rearing water[END_REF][START_REF] Kandel | Abundance, diversity and seasonal dynamics of predatory bacteria in aquaculture zero discharge systems[END_REF]. Given their predation obligation and distribution in nature and diversity, BALOs may have a strong impact on the dynamics and structure of bacterial communities. However, studies related to the ecological impact of BALOs on microbial communities are scarce.

BALOs have been seldom studied in freshwater lakes. Recently, these bacteria have been shown to favour eutrophic conditions, both in terms of abundance and diversity, due to a higher number and diversity of potential prey [START_REF] Chauhan | Increased diversity of predacious Bdellovibrio-like organisms (BLOs) as a function of eutrophication in Kumaon Lakes of India[END_REF]. It is noteworthy that BALOs have also been detected in less nutrient-replete lakes, such as Lake Geneva, Bourget, and Annecy, characterised by meso-to oligotrophic status [START_REF] Paix | Diversity, dynamics and distribution of Bdellovibrio and like organisms in peri-alpine lakes[END_REF][START_REF] Ezzedine | New 16S rRNA primers to uncover Bdellovibrio and like organisms diversity and abundance[END_REF][START_REF] Ezzedine | Bdellovibrio and like organisms in Lake Geneva: an unseen elephant in the room?[END_REF]. However, in these lakes, BALOs have been studied away from shore, at reference stations 1 3 of the lakes, on a monthly or seasonal time scale, i.e. at a temporal scale unsuited to studying microbial predator dynamics. In Lake Geneva, probably one of the most studied lakes regarding BALOs, Bdellovibrionaceae seems to be the most diverse family in terms of OTUs, while Peredibacteraceae are the most abundant, as measured by quantitative PCR (qPCR). It was also found that potential prey mainly belonged to the genus Pseudomonas [START_REF] Ezzedine | Bdellovibrio sp.: an important bacterial predator in Lake Geneva?[END_REF][START_REF] Chauhan | Increased diversity of predacious Bdellovibrio-like organisms (BLOs) as a function of eutrophication in Kumaon Lakes of India[END_REF] are predated by Bdellovibrionaceae.

Assessing the dynamics of BALOs and the impact of environmental forcing on this functional group can be significantly facilitated by analysing them on short time scales and under relatively controlled experimental conditions. This is what in situ deployed mesocosms can typically provide, as they constitute an ideal experimental tool, intermediate between microcosms commonly used in laboratories and in situ ecological surveys [START_REF] Mcallister | Measurements of primary production in coastal sea water using a large volume plastic sphere[END_REF][START_REF] Antia | Further measurements of primary production using a large-volume plastic sphere[END_REF][START_REF] Egge | Blooms of phytoplankton including Emiliania huxleyi (Haptophyta). Effects of nutrient supply in different N:P ratios[END_REF]. In the past, in situ deployed mesocosms have proved to be a relevant type of instrumentation for observing changes in complex planktonic communities while controlling environmental conditions such as brightness, concentrations of nutrients or dissolved gases, water mixing, and predatory effects, among others. [START_REF] Mostajir | Experimental test of the effect of ultraviolet-B radiation in a planktonic community[END_REF][START_REF] Thingstad | Conceptual models for the biogeochemical role of the photic zone microbial food web, with particular reference to the Mediterranean Sea[END_REF][START_REF] Jacquet | Short-timescale variability of picophytoplankton abundance and cellular parameters in surface waters of the Alboran Sea (western Mediterranean)[END_REF][START_REF] Havskum | Routine quantification of phytoplankton groups-microscopy or pigment analyses?[END_REF]. We conducted a mesocosm experiment in Lake Geneva, designated project "MESOLAC," which aimed at analysing the response of surface planktonic populations to extreme events caused by global climatic change. The latter has already resulted in an increase in the number and strength of the so-called extreme events, including storms, heavy rainfall, floods, or drought [START_REF] Lavell | Climate change: new dimensions in disaster risk, exposure, vulnerability, and resilience[END_REF]. Ecosystems are likely to be massively impacted by such events, and lakes are no exception [START_REF] Jenny | Scientists' warning to humanity: rapid degradation of the world's large lakes[END_REF]. For example, Woolway and Merchant [START_REF] Woolway | Worldwide alteration of lake mixing regimes in response to climate change[END_REF] reported alterations in the thermal regime of lakes, which could result in planktonic diversity reduction and the inability of lake ecosystems to return to their original state due to the loss of their functional redundancy. Therefore, addressing issues regarding lacustrine functions and their resilience is of major importance to better understand the impact of climate change, including global warming, on their state and evolution in the near future.

The focus of this work was to study (i) the dynamics and diversity of three main BALO families belonging to the Oligoflexia group (e.g. Bdellovibrionaceae, Bacteriovoracaceae and Peredibacteraceae) on a short time scale and (ii) their response to simulated extreme events and abiotic and biotic factors. We hypothesised that the addition of dissolved organic carbon (simulating an input of carbon from the watershed in response to runoff from the catchment), light reduction, and mixing (simulating stormy conditions) would stimulate bacterial growth and thus the predation by bacterial predators such as BALOs. Using flow cytometry, allowing us to discriminate between two groups of bacteria, namely low DNA content [START_REF] Gasol | Significance of size and nucleic acid content heterogeneity as measured by flow cytometry in natural planktonic bacteria[END_REF] or low nucleic acid content [START_REF] Lebaron | Does the high nucleic acid content of individual bacterial cells allow us to discriminate between active cells and inactive cells in aquatic systems?[END_REF] (LDNA) and high DNA content [START_REF] Gasol | Significance of size and nucleic acid content heterogeneity as measured by flow cytometry in natural planktonic bacteria[END_REF] or high nucleic acid content [START_REF] Lebaron | Does the high nucleic acid content of individual bacterial cells allow us to discriminate between active cells and inactive cells in aquatic systems?[END_REF] (HDNA), we also hypothesised that the supply of nutrients may favour the HDNA, likely to constitute preferential prey for BALOs.

Materials and Methods

Experimental Setup

The MESOLAC experiment was conducted on Lake Geneva over 4 weeks in July 2019. Nine experimental mesocosms (Supplementary Fig. 1) were placed approximately 100 m away from the coast and at a depth of 7 m near the Alpine Centre for Research on Lake Ecosystems and Food Webs (CARRTEL, Thonon-les-Bains; 46° 22′ 09.1″ N, 6° 27′ 15.3″ E). An ecological anchoring system, set up by professional scuba divers, including the corresponding author, ensured the fixation of each mesocosm and good preservation of the bottom (i.e. sediments and plants) of the lake. Mesocosms consisted of reinforced polythene bags of 4.5 m height and 1.4 m diameter at their upper level with a finishing conical part towards the bottom (Insinööritoimisto Haikonen Oy, Norra Paipis, Finland). The mesocosms were supported at every metre by a plastic frame and a double system of buoys on top. All bags were filled with water simultaneously on the same day (1 July) within a few hours and were untouched for 3 days to acclimate. The total volume of water in each mesocosm was approximately 4 m 3 . The experiment started on 4 July, designated T0. The experimental design included three treatments, each replicated thrice. Each bag was covered by filters (Lee filters) applied on the surface of the mesocosm, to protect embedded water from bird droppings and other external elements, so that they were all treated similarly. The control treatment (C) consisted of minimal light reduction (i.e. 5%). Treatment M consisted of a medium stress situation, i.e. 30% light reduction, dissolved organic carbon (i.e. total DOC ~ 2 mg L -1 ; using autoclaved peat soil extract) at a concentration × 1.5 times that of the control and regular manual mixing for 5 min daily for 2 weeks. Treatment H consisted of a shorter but stronger intensity stress simulation. For 5 days, light was reduced by ~ 85%, DOC concentration was increased fivefold (i.e. total DOC ~ 6 mg L -1 ), and daily mixing was performed for 15 min. After 5 days, treatment H was exposed to the control conditions. More details can be found in [START_REF] Tran-Khac | In situ pelagic dataset from continuous monitoring: a mesocosm experiment in Lake Geneva (MESOLAC)[END_REF].

Sampling Strategy and Analysed Parameters

Water samples were collected at a depth of 2 m in each bag, with a 2-4-day interval from 4 July to 22, July 2019. The sampling dates are 4th, 8th, 11th, 16th, 18th, and 22nd of July and are, respectively, relabelled as 0 day, 4 days, 7 days, 12 days, 14 days, and 18 days. Physicochemical parameters such as pH, dissolved oxygen concentration, turbidity, phosphorus, and nitrogen concentrations were measured following standard protocols [START_REF] Rimet | The Observatory on LAkes (OLA) database: sixty years of environmental data accessible to the public[END_REF]. Chlorophyll a concentration was obtained using spectrophotometry after pigment extraction in 90% acetone [START_REF] Tran-Khac | In situ pelagic dataset from continuous monitoring: a mesocosm experiment in Lake Geneva (MESOLAC)[END_REF]. Meteorological data were collected daily via the 1 3

CLIMATIK platform, which is reserved for INRAE research. For DNA analysis, 200-350 mL of water was filtered through 0.2-µm PC filters. All filters were stored at -20 °C until DNA extraction. Water samples were also taken for delayed flow cytometry analyses and consisted of 2 mL of water fixed with glutaraldehyde (15 min, 0.5%), flash-frozen in liquid nitrogen, and stored at -80 °C.

DNA Extraction

The 0.2-µm filters were subjected to DNA extraction using a homemade protocol. In a 2-mL Eppendorf tube, 300 µL of TE buffer (TRIS, 1 M [pH 8] and EDTA, 0.5 M [pH 8]) was added to each filter. Next, a lysis step was performed by adding 200 µL of lysis solution (TRIS, 1 M [pH 8]; EDTA, 0.5 M [pH 8], and sucrose: 0.7 M). After a thermic shock at -80 °C for 15 min and at 55 °C for 2 min, 50 µL of 10% sodium dodecyl sulphate (SDS), and 10 µL of proteinase K (20 mg/mL) were added. The samples were then incubated at 37 °C for 1 h with gentle stirring and placed on a heating block at 55 °C for 20 min. After a quick centrifugation step (13,000 rpm at 4 °C for 3 min), the supernatant was collected. Afterwards, 50 µL of sodium acetate (3 M [pH 5.2]) and 1.5 µL of 25 µg/µL GenEluteTM-LPA (Sigma-Aldrich, St Louis, MO, USA) were added. Subsequently, one volume of isopropanol was added, and the tubes were centrifuged for 10 min at 12,000×g at 4 °C. Following this step, two rounds of ethanol (80%) washing were carried out to clean the DNA pellet. The remaining ethanol was evaporated using a SpeedVac for 20 min. Finally, 30 µL of TE was added, and samples were incubated at 37 °C for 1 h to allow the pellet to gently dissolve into the TE buffer. The DNA concentration was measured using a NanoDrop 1000 spectrophotometer. For DNA concentrations greater than 25 ng µL -1 , a dilution was performed. All DNA preparations were stored at -20 °C until further analysis.

Quantification of BALOs by qPCR

Qiagen's Rotor-Gene Q machine and QuantiTect SYBR Green PCR kit (Qiagen, Hilden, Germany) were used to amplify the BALOs. Standard curves were prepared using identified clones of Bdellovibrionaceae, Peredibacteraceae, and Bacteriovoracaceae from our previous study [START_REF] Ezzedine | New 16S rRNA primers to uncover Bdellovibrio and like organisms diversity and abundance[END_REF]. Briefly, plasmids were extracted and purified using a NucleoSpin Plasmid kit (Macherey-Nagel GmbH & Co. KG, Düren, Germany) according to the manufacturer's instructions. Plasmids were then digested using a BamHI restriction enzyme (Takara Bio Inc., Shiga, Japan) following the manufacturer's instructions. The digested plasmid concentrations were measured using a Quant-iT PicoGreen ds DNA Reagent kit (Invitrogen, Carlsbad, CA, USA), and fluorescence was read using a Fluoroskan Ascent FL plate reader. The number of copies for each BALO clone was calculated using the following formula: Serial dilutions were then conducted from 10 9 to 10 0 copies. Diluted DNA from 10 7 to 10 0 was duplicated and amplified by the BALOs' qPCR set of primers to construct the standard curve. The list of primers used to amplify the 16S rDNA gene is shown in Supplementary Table 1 [START_REF] Ezzedine | New 16S rRNA primers to uncover Bdellovibrio and like organisms diversity and abundance[END_REF][START_REF] Van Essche | Development and performance of a quantitative PCR for the enumeration of Bdellovibrionaceae[END_REF]. All primer sets for BALO amplification were optimised in our previous study [START_REF] Ezzedine | New 16S rRNA primers to uncover Bdellovibrio and like organisms diversity and abundance[END_REF]. Two controls were added at each time point. The Bdellovibrionaceae standard curve had an efficiency of 0.92, R 2 of 0.996, and the threshold was set to 0.02. For Peredibacteraceae, the efficiency was 0.87, R 2 was 0.87, and the threshold was fixed at 0.015. For Bacteriovoracaceae, the efficiency was 0.96, R 2 was 0.993, and the threshold was set to 0.02. The qPCR mixture volume was 25 µL and consisted of final concentration: 1 X Master Mix (QuantiTect SYBR Green PCR kit, Qiagen), 0.3 mg mL -1 of BSA, 0.2 µM of forward and reverse primers, and 1 µL of template DNA (25 ng µL -1 ). The program was as follows: 95 °C for 15 min followed by 40 cycles of 95 °C for 45 s, 60 °C for 45 s, and 72 °C for 45 s, with + 1 °C every 5 s from 60 to 95 °C. For all environmental samples, those that failed to amplify or were outside the standard curve limits were not considered in the analysis. BALO abundances were obtained in copy per reaction and were transformed to copy per millilitre of filtered water using the following formula:

Flow Cytometry Analysis

We used a FACSCalibur flow cytometer (BD BioSciences, Franklin Lakes, NJ, USA) to determine the total prokaryote abundance. After each water sample was thawed at ambient temperature, 2.5 µL was added to 245 µL filtered (< 0.02 µm) TE buffer and 2.5 µL of SYBR Green I (diluted 10,000 times). The sample was then heated for 10 min at 75 °C before the FCM analysis [START_REF] Jacquet | A few tests prior to flow cytometry and epifluorescence analyses of freshwater bacterioand virioplankton communities[END_REF]. "Listmode" files were exported and analysed using CYTOWIN [START_REF] Vaulot | A simple method to preserve oceanic phytoplankton for flow cytometric analyses[END_REF]. The analysis provided information on prokaryotelike particles (PLPs), which encompassed two subgroups designated "high" (HDNA) and "low" DNA contents (LDNA), indicating possibly active and less productive/ active bacterial populations, respectively [START_REF] Lebaron | Does the high nucleic acid content of individual bacterial cells allow us to discriminate between active cells and inactive cells in aquatic systems?[END_REF]. Virus-like particles (VLPs) and the virus-to-prokaryote ratio (VPR) were also analysed [START_REF] Marie | Enumeration of marine viruses in culture and natural samples by flow cytometry[END_REF].

Statistics

Statistical analyses were performed using R software version 3.6.3 [START_REF] Core | R: A language and environment for statistical computing[END_REF], and figures were produced using the ggplot2 package [START_REF] Wickham | ggplot2: Create Elegant Data Visualisations Using the Grammar of Graphics[END_REF] or ggpubr [START_REF] Kassambara | ggplot2" Based Publication Ready Plots[END_REF]. One-way ANOVA (p < 0.05) or Kruskal-Wallis test followed by Dunn's test [START_REF] Dinno | Dunn's Test of Multiple Comparisons Using Rank Sums[END_REF] or Tukey's HSD post hoc tests were performed to compare the mean values of each treatment throughout the experiment. NMDS and Adonis tests were also performed to test the differences between groups, i.e. replicates, treatments, and dates. Spearman's rank correlation test was conducted to examine the correlation between variables, and a correlogram was drawn using the corrplot package [START_REF] Wei | Visualization of a correlation matrix[END_REF] following the Antoine Soetewey code (https:// stats andr. com/ blog/ corre lation-coeffi cient-and-corre lation-test-in-r/.). The correlogram shows the correlation coefficient for all pairs of variables (with more intense colours for more correlations), and correlations not statistically significant are represented by a white box. Moreover, to visualise the variation in microbial community structure and the relationship between environmental variables and sample clusters, we performed a redundancy analysis (RDA). All RDAs performed were significant according to the ANOVA test and chosen according to the detrended correspondence analysis (DCA) [START_REF] Ramette | Multivariate analyses in microbial ecology[END_REF]. RDA was performed using the Vegan package [START_REF] Oksanen | Package "vegan" title community ecology package[END_REF]. Environmental parameters were selected using the envfit function (p < 0.05). Multicollinear variables were removed using variance inflation factors (VIF) [START_REF] Oksanen | Package "vegan" title community ecology package[END_REF].

Results

BALOs' Abundance and Dynamics

qPCR measurements indicated that Peredibacteraceae abundance was the highest, followed by Bdellovibrionaceae and Bacteriovoracaceae. The abundance of Bdellovibrionaceae (Fig. 1A) in the H and C treatments seemed to follow the same decreasing pattern until day 14, with Bdellovibrionaceae being less abundant in treatment H. The abundance of Bdellovibrionaceae slightly increased in treatment M until day 4 before decreasing as the others. Bdellovibrionaceae abundance varied from 4.27 × 10 2 to 4.30 × 10 2 copies mL -1 in C, from 6.24 × 10 1 to 1.91 × 10 2 copies mL -1 in H, and from 1.68 × 10 2 to 2.06 × 10 2 copies mL -1 in M. While statistical analyses revealed that treatments H and M were different (Dunn's test p value = 0.02), no significant difference was observed between treatments C and H, or between C and M (Fig. 2A). Peredibacteraceae (Fig. 1B) and Bacteriovoracaceae (Fig. 1C) dynamics showed a similar pattern to those of Bdellovibrionaceae, with a decrease in abundance that was more marked in Bacteriovoracaceae, Peredibacteraceae varied from 4.74 × 10 3 to 2.59 × 10 2 copies mL -1 in C, from 9.61 × 10 2 to 8.27 × 10 1 copies mL -1 in H, and from 2.46 × 10 3 to 1.88 × 10 2 copies mL -1 in M. For Bacteriovoracaceae abundance varied from 2.88 × 10 3 to 5.51 × 10 0 copies mL -1 in C, from 8.38 × 10 2 to 5.41 × 10 0 copies mL -1 in H, and from 2.12 × 10 3 to 8.39 × 10 0 copies mL -1 in M. Kruskal-Wallis test showed no significant difference between treatments for both Peredibacteraceae (p value = 0.19) and Bacteriovoracaceae (p value = 0.34). The Adonis test confirmed that the treatments were not responsible for the disparity (p value = 0.196). Moreover, the same test confirmed that the replicate mesocosms did not differ (p value = 0.588). However, differences existed between dates (p value < 0.01), and three different clusters could be defined: days 0-4, 7, and 12-18 (Fig. 3).

Relationships with Abiotic Factors

The dynamics of the main environmental parameters during the experiment are described in Supplementary Text Box 1 and presented in Supplementary Fig. 2. On the one hand, the correlogram (white boxes for non-significant correlation; Fig. 4) shows that Bdellovibrionaceae, unlike the other two BALOs, displayed very few relationships with measured environmental variables. Correlations were observed between pH (ρ = 0.30) and water hardness (ρ = 0.27). On the other hand, Peredibacteraceae and Bacteriovoracaceae were positively and significantly correlated with water hardness, nitrite (NO

), total nitrogen (N tot ), phosphate (PO

), total phosphorus (P tot ), chloride (Cl -), temperature, and particulate phosphorus (P part ), and negatively correlated with ammonium (NH 4 + ), silica (SiO 2 ), and chlorophyll a (Chl a). The RDA, which explained 65% of the variance (Fig. 5A), showed that Bdellovibrionaceae correlated positively with water hardness, pH, and day 18; Peredibacteraceae with nitrite (NO 2 -) and chloride (Cl -), and Bacteriovoracaceae with temperature and day 0. The correlogram confirmed these observations. Overall, environmental parameters explained, to some extent, the variation in abundance

Relationships Between Biotic Variables

The correlogram (Fig. 4) suggested that Bdellovibrionaceae was not significantly correlated with Prokaryote-like particles (PLP), HDNA, or LDNA. Contrastingly, both Peredibacteraceae and Bacteriovoracaceae had a significant positive relationship with PLP and HDNA and a significant negative relationship with LDNA. More specifically, Peredibacteraceae was more correlated with HDNA (ρ = 0.81) than Bacteriovoracaceae (ρ = 0.60), whereas the opposite was true for LDNA (ρ = -0.51 vs -0.66). Furthermore, Peredibacteraceae had a superior relationship with PLP (ρ = 0.53) than Bacteriovoracaceae (ρ = 0.30). Considering the evolution of PLP, HDNA, and LDNA during the experiment (Fig. 1D-F), we observed that the trend was generally similar to that observed for BALOs. The curve of the M treatment was slightly higher than that of the C and H treatments. Moreover, except for LDNA, abundance tended to decrease during the experiment-LDNA abundance increased during the experiment until days 13-14. PLPs (Fig. 1D) decreased during the experiment in all mesocosms, from 5.44 × 10 6 to 4.10 × 10 6 cells mL -1 for C, from 5.14 × 10 6 to 2.63 × 10 6 cells mL -1 for H, and from 5.13 × 10 6 to 4.26 × 10 6 cells mL -1 for M. The same pattern was also observed for the HDNAs (Fig. 1E) with a decrease in abundance from 3.80 × 10 6 to 1.65 × 10 6 cells mL -1 for C, from 3.67 × 10 6 to 1.10 × 10 6 cells mL -1 for H, and from 3.65 × 10 6 to 1.18 × 10 6 cells mL -1 for M. The LDNA dynamics were completely different (Fig. 1F) since the proportion of this group increased significantly during the experiment (from 1.63 × 10 6 to 2.5 × 10 6 cells mL -1 for C, from 1.47 × 10 6 to 1.58 × 10 6 cells mL -1 for H, and from 1.59 × 10 6 to 2.00 × 10 6 cells mL -1 for M). It is also noteworthy that the HDNAs were more abundant than the LDNAs at the beginning of the experiment for all treatments, which was reversed at the end. The ANOVA test showed that HDNA and LDNA mean abundances were not significantly different between the treatments (p value = 0.37 and 0.11, respectively). However, treatments impacted PLP, and differences were significant between C and H (p value = 0.01) and between M and H (p value = 0.007). As shown in Figs. 1D and2B, the H treatment had less PLP than C and M. Moreover, the correlogram in Fig. 4 shows that HDNA and PLP were positively and significantly correlated (ρ = 0.77). In contrast, LDNA was not significantly correlated with PLP and was significantly and negatively correlated with HDNA. PLP, LDNA, and HDNA had 8, 12, and 14 positive or negative significant relationships with environmental variables, respectively. LDNA had mostly negative relationships with environmental parameters, whereas HDNA had mostly positive relations. The RDA (Fig. 5B), which explained 70% of the variance in PLP, HDNA, and LDNA, revealed that PLP correlated positively with water hardness and day, HDNA with nitrite and chloride, and LDNA with chlorophyll a and day 12.

The present study also allowed us to examine VLP dynamics that were found to be opposite to those of BALOs and prokaryotes (Supplementary Fig. 3). VLP abundance increased more in H than in the C and M treatments; however, the ANOVA test revealed that the differences were not significant (p value = 0.5). Moreover, the correlogram ), total carbon (C tot ), total phosphorus (P tot ), total nitrogen (N tot ), temperature (temp), water hardness (Wt hardness), prokaryote-like particles (PLP), high DNA content (HDNA), low DNA content (LDNA), virus-like particles (VLP), and virus-to-prokaryote ratio (VPR) 1 3

(Fig. 4) suggested the absence of a significant correlation between VLP and BALOs. No significant correlations were observed between VLP and PLP, HDNA, or LDNA. When considering the virus-to-prokaryote ratio (VPR), H favoured an increase in VPR at the beginning of the experiment (until day 12; Supplementary Fig. 3B), and statistical analyses indicated a significant difference between H and M, but not between C and H or between C and M. The correlogram (Fig. 4) indicated significant and negative correlations between VPR and Peredibacteraceae (ρ = -0.49) and Bacteriovoracaceae (ρ = -0.32), but not with Bdellovibrionaceae. VPR also displayed a strong negative relationship with PLP (ρ = -0.86) and HDNA (ρ = -0.67), but not with LDNA.

Discussion

With global climate change, episodic extreme events (such as floods and violent storms) are expected to become more frequent, causing changes in lake subsidies and exposing them to nutrient pulses and transfer imbalances. The effects of allochthonous dissolved organic matter inputs on aquatic organisms, biotic interactions and metabolic processes feedback, and the impact on the littoral versus pelagic communities remain largely underexplored. Our experimental approach made it possible to simulate such a scenario with the main aim of investigating the response of a specific functional group of bacteria.

To study microbial interactions, especially predator-prey relationships, the time scale is important. For BALOs, predation is known to be rapid under controlled and favourable conditions, but the impacts of bacterial predators on bacterial prey dynamics, distribution, and composition under natural conditions have been poorly investigated until now. Using a 2-4-day sampling strategy and an outdoor mesocosm approach, we aimed to assess, for the first time, the dynamics of the main BALO families belonging to Oligoflexia, highlighting their response to abiotic and biotic factors.

After observing the BALO dynamics and performing statistical analyses, it appears that the simulated effects did not have a strong impact on BALOs' abundance. Said differently, BALOs were not affected by the stresses or modifications caused by extreme events. As predation is an important contributor to community functional dynamics and diversity, observing that the predators' dynamics is stable in the face of such disturbance is rather reassuring. We are aware, however, that further testing would be useful to confirm such results. Our result may be an extra observation that lake responses to extremes disturbance can be dependent on lake condition antecedents [START_REF] Havens | The importance of considering biological processes when setting total maximum daily loads (TMDL) for phosphorus in shallow lakes and reservoirs[END_REF][START_REF] Havens | Composition, size, and biomass of zooplankton in large productive Florida lakes[END_REF][START_REF] Havens | Effects of climate variability on cladoceran zooplankton and cyanobacteria in a shallow subtropical lake[END_REF][START_REF] Perga | Storm impacts on alpine lakes: antecedent weather conditions matter more than the event intensity[END_REF][START_REF] Stockwell | Storm impacts on phytoplankton community dynamics in lakes[END_REF] and, in some situations, although exposed to severe climatic conditions (e.g. storms), lakes may not be strongly modified because lake condition antecedents shape resistance and resilience of the ecosystem [START_REF] Thayne | Antecedent lake conditions shape resistance and resilience of a shallow lake ecosystem following extreme wind storms[END_REF].

When the treatment effects are not considered, it is clear that each family of BALOs displayed a very different dynamic, suggesting a variety of interactions and responses to forcing (abiotic and biotic interactions). Many factors and/ or processes could be responsible for these differences (e.g. the predation spectrum, prey presence/absence, resistance genes or other attributes, etc.). The Peredibacteraceae family was the most abundant among the studied BALOs. Despite their abundance decreasing slightly during the experiment, this was consistent with our previous studies on BALOs in Lake Geneva [START_REF] Paix | Diversity, dynamics and distribution of Bdellovibrio and like organisms in peri-alpine lakes[END_REF][START_REF] Ezzedine | New 16S rRNA primers to uncover Bdellovibrio and like organisms diversity and abundance[END_REF][START_REF] Ezzedine | Bdellovibrio and like organisms in Lake Geneva: an unseen elephant in the room?[END_REF][START_REF] Ezzedine | Bdellovibrio sp.: an important bacterial predator in Lake Geneva?[END_REF], suggesting that this group is the dominant one among the targeted BALOs, especially, in lake surface waters. In contrast, Bacteriovoracaceae were less abundant and decreased strongly during the experiment. This is also in accordance with another study of our team [START_REF] Ezzedine | Exploring the diversity, abundance, structure of Bdellovibrio and like organisms in four contrasted ecosystems over a year[END_REF]. Bdellovibrionaceae abundance remained relatively stable.

Bdellovibrionaceae displayed only relationships with pH and water hardness. It has been established in the literature that some bacterial species are vastly dependent on pH [48,[START_REF] Stockner | Autotrophic picoplankton: community composition, abundance and distribution across a gradient of oligotrophic British Columbia and Yukon Territory Lakes[END_REF]. In laboratory conditions, Bdellovibrio sp. motility has been shown to be inhibited at pH < 5 and pH > 9. The optimum attachment and growth for Bdellovibrio sp. ranged from neutrality to slightly alkaline pH [START_REF] Dashiff | Predation of human pathogens by the predatory bacteria Micavibrio aeruginosavorus and Bdellovibrio bacteriovorus[END_REF][START_REF] Varon | Interaction of Bdellovibrio bacteriovorus and Host Bacteria[END_REF][START_REF] Seidler | Isolation and characterization of hostindependent Bdellovibrios[END_REF]. In Lake Geneva, pH varied slightly, i.e. between 7.3 and 8.7, corresponding to the optimal range for Bdellovibrio growth. Further laboratory experiments will be welcomed to provide more pieces of information on pH optima for Bdellovibrionaceae and other BALOs of Lake Geneva. Peredibacteraceae had marked relationships with nitrite and chloride ions, what no other studies have shown before. These variables also varied very little during the experiment. For the chloride ion, [START_REF] Roeßler | Chloride dependence of growth in bacteria[END_REF] [START_REF] Roeßler | Chloride dependence of growth in bacteria[END_REF] demonstrated that chloride ions are not necessary for growth in 44 different bacterial strains, except in very salty environments. Naturally, this is not the case for Lake Geneva. The only ions reported to be beneficial to the growth of BALOs are calcium and magnesium ions because they are important for BALO penetration into prey [START_REF] Huang | Effects of calcium and magnesium ions and host viability on growth of bdellovibrios[END_REF]. For nitrite, Peredibacteraceae may be associated with this variable because they could consume nitrite-oxidising bacteria or use nitrite to survive under anoxic conditions. In fact, other BALOs such as Bdellovibrio bacteriovorus UP and Micavibrio spp. are capable of preying on Nitrospira in sludge flocs [START_REF] Feng | Predation by Bdellovibrio bacteriovorus significantly reduces viability and alters the microbial community composition of activated sludge flocs and granules[END_REF][START_REF] Dolinšek | Interactions of nitrifying bacteria and heterotrophs: identification of a Micavibrio-like putative predator of Nitrospira spp[END_REF]. Also, [START_REF] Sockett | Bdellovibrio as therapeutic agents: a predatory renaissance?[END_REF] [START_REF] Sockett | Bdellovibrio as therapeutic agents: a predatory renaissance?[END_REF] found that within the B. bacteriovorus HD100 genome there are genes that encode for nitrite reductase and nitric oxide reductase. This result suggests that B. bacteriovorus HD100 strain could use nitrite as a terminal electron acceptor instead of oxygen. In fact, Monnappa et al. (2013) [START_REF] Monnappa | Application of bacterial predation to mitigate recombinant bacterial populations and their DNA[END_REF] confirmed that B. bacteriovorus HD100 predated upon Escherichia coli under anoxic conditions when at least 1 mM of nitrate was available in the culture medium. Looking into the genome of Peredibacteraceae of Lake Geneva could help to understand better such a relationship with nitrogen. For Bacteriovoracaceae temperature seemed to be important, and literature has indicated temperature's impact on the growth and abundance of BALOs [START_REF] Williams | Ecology of the predatory Bdellovibrio and like organisms[END_REF][START_REF] Kandel | Abundance, diversity and seasonal dynamics of predatory bacteria in aquaculture zero discharge systems[END_REF][START_REF] Schoeffield | A comparison of the survival of intraperiplasmic and attack phase bdellovibrios with reduced oxygen[END_REF]. Bacteriovoracaceae abundance decreased with increasing temperature; however, the temperature difference in the present study was 3 °C, which is small considering that BALOs generally tolerate a wide temperature range [START_REF] Jurkevitch | Phylogenetic diversity and evolution of predatory prokaryotes[END_REF][START_REF] Davidov | Diversity and evolution of Bdellovibrio-and-like organisms (BALOs), reclassification of Bacteriovorax starrii as Peredibacter starrii gen. nov., comb. nov., and description of the Bacteriovorax-Peredibacter clade as Bacteriovoracaceae fam. nov[END_REF]. A study by [START_REF] Kandel | Abundance, diversity and seasonal dynamics of predatory bacteria in aquaculture zero discharge systems[END_REF] [START_REF] Kandel | Abundance, diversity and seasonal dynamics of predatory bacteria in aquaculture zero discharge systems[END_REF] showed that a variation of 4 °C in a closed system (aquaculture) does not affect BALO abundance. In general, it takes large temperature fluctuations, as in natural systems, to observe effects on the BALO population.

We know very little about BALOs' prey preference and interactions with other microorganisms in general. This study did not cover this issue, among the Prokaryotelike particles (PLP), we discriminated bacteria with high (HDNA) vs low (LDNA) DNA or nucleic acid content, some of which might be considered as potential BALO prey. Indeed, BALOs have a prey minimum size [START_REF] Summers | Predation strategies of the bacterium Bdellovibrio bacteriovorus result in overexploitation and bottlenecks[END_REF] and benefit more from prey with high amount of nutrient (more progeny) such as nucleic acids. Hence, we hypothesised that BALOs may preferably prey on HDNA, which would also benefit from the initial load of nutrients in treatments H and M. It is also noteworthy that the obligate oligotrophic LDNA cells could get richer in nucleic acids and shift to HDNA category [START_REF] Santos | Flow cytometry analysis of low/high DNA content (LNA/HNA) bacteria as bioindicator of water quality evaluation[END_REF][START_REF] Wang | Isolation and characterization of low nucleic acid (LNA)-content bacteria[END_REF]. After eliminating the treatment effect, the statistical tests' results did not show significant differences in the abundances of HDNA and LDNA. However, during the experiment, LDNA tended to increase while the reverse trend was observed for HDNA. Clearly, different dynamics occurred on different days as shown in Fig. 3. On the other hand, there was a significant difference in PLP abundances between C and H and M and H. The events that the H mesocosms were exposed to induced a significant decrease in PLPs compared with the M treatment and the control. This result was the opposite of what we expected. Our initial hypothesis was that our experimental conditions, especially the increase in dissolved organic carbon, would enhance bacterial growth and reproduction if temperature and inorganic nutrients were not limiting at that time, and by extension, promote the growth of BALOs and other bacterial predators. In fact, BALOs have been reported to be more abundant and diverse in eutrophic lakes owing to the abundance and diversity of the population of prey bacteria [START_REF] Chauhan | Increased diversity of predacious Bdellovibrio-like organisms (BLOs) as a function of eutrophication in Kumaon Lakes of India[END_REF]. To explain the overall absence of significant differences between treatments, we first surmised that one cannot rule out the possibility that a weak increase in organic matter input will have little or no effect on prokaryotic communities, as already shown elsewhere [START_REF] Kankaala | Impacts of added dissolved organic carbon on boreal freshwater pelagic metabolism and food webs in mesocosm experiments[END_REF]. Secondly, in another study with smaller mesocosms than the ones deployed here [START_REF] Kim | The influence of dissolved organic carbon on the microbial community associated with Tetraselmis striata for bio-diesel production[END_REF], a difference in microorganisms population could be observed but with DOC concentration considerably higher, i.e. 428.7 mg L -1 . Thirdly, it is possible that other elements were limiting for prokaryotic growth [START_REF] Jacquet | Effects of inorganic and organic nutrient addition on a coastal microbial community (Isefjord, Denmark)[END_REF]. However, it was not the case for inorganic nutrients because ammonium, nitrite, nitrate, nitrogen, and phosphorus concentrations were at sufficient concentrations to sustain microorganism growth [START_REF] Tran-Khac | In situ pelagic dataset from continuous monitoring: a mesocosm experiment in Lake Geneva (MESOLAC)[END_REF]. We do not think that autoclaving the peat extract soil did any damage on the nutrient composition as reported by Williams-Linera and Ewel (1984) [START_REF] Williams-Linera | Effect of autoclave sterilization of a tropical andept on seed germination and seedling growth[END_REF]. However, heterotrophic bacteria might not have efficiently assimilated the organic carbon from the soil peat extract. Indeed, [START_REF] Guillemette | Selective consumption and metabolic allocation of terrestrial and algal carbon determine allochthony in lake bacteria[END_REF] [START_REF] Guillemette | Selective consumption and metabolic allocation of terrestrial and algal carbon determine allochthony in lake bacteria[END_REF] showed that the terrestrial organic carbon which contains more humic substances and polysaccharides is more recalcitrant to bacterial degradation, while on the contrary algal-derived organic carbon is more labile and easier to incorporate. Additionally, Zhou et al. (2020) [START_REF] Zhou | Resource aromaticity affects bacterial community successions in response to different sources of dissolved organic matter[END_REF] reported that soilderived dissolved organic matter has a higher environmental filtering on bacteria where only those who have the adequate arsenal to degrade such compounds can benefit. Finally, we assumed that the frequency of mixing the mesocosms might not have been enough to let microorganisms benefit from DOC, as pointed out by [START_REF] Striebel | Different mixing techniques in experimental mesocosms-does mixing affect plankton biomass and community composition?[END_REF] [70] who advised to mix three times the mesocosms when using a disc and to preferably instal an automatic system to mix the water more frequently. Regardless of the small decrease in PLPs, their abundance remained relatively high during the experiment. In general, the PLPs and HDNAs decreased slightly over time. Conversely, LDNA slightly increased. This was unexpected as stated above. Regarding the decrease in PLPs and HDNAs, we surmised that the disruption in the number of prokaryotes, and by extension of potential prey, especially the most active one, HDNA, caused a decrease in BALO abundance. BALOs also obviously did not benefit from the increase in LDNAs, which contain less nucleic acid. This was supported by the correlogram results (Fig. 4). Bacteriovoracaceae and Peredibacteraceae showed a strong positive correlation with PLPs, especially HDNAs, and conversely, they both had a negative relationship with LDNAs, which is consistent with our hypothesis that HDNAs could be more preyed upon. The observed decrease in PLPs, particularly HDNAs, could thus have disturbed BALO dynamics in the mesocosms. Bdellovibrionaceae was not significantly correlated with PLP, HDNA, or LDNA. Other factors seemed to explain the dynamics over time. Furthermore, BALOs are known to mainly hunt for Gram-negative bacteria; however, not all Gram-negative bacteria are susceptible to BALO predation. Moreover, BALOs' species and strains have different prey spectrum and can be either generalist, specialist, or versatilist [START_REF] Shemesh | Samll eats big: ecology and diversity of Bdellovibrio and like organisms, and their dynamics in predator-prey interactions[END_REF][START_REF] Chen | Prey bacteria shape the community structure of their predators[END_REF]. Unfortunately, the prey spectrum of BALOs in Lake Geneva is unknown. However, we could suggest that the strategy of each studied BALOs resulted probably in a different outcome with an important decline of Bacteriovoracaceae, a relative stability of Bdellovibrionaceae, and an entre-deux for Peredibacteraceae. In addition to the predation spectrum issue, another explanation can be proposed for the striking absence of BALOs' concentration increase i.e. decoy microorganisms. It is known that non-prey and decoy cells can impact negatively BALOs' attachment and growth [73][START_REF] Shemesh | Plastic phenotypic resistance to predation by Bdellovibrio and like organisms in bacterial prey[END_REF][START_REF] Hobley | Bdellovibrio predation in the presence of decoys: three-way bacterial interactions revealed by mathematical and experimental analyses[END_REF]. The peat soil extract could have introduced decoy cells that hindered BALOs' predation, typically heated cells with intact cell that could be infected by BALOs, likely to result in various BALOs' species growth efficiency [START_REF] Varon | Interaction of Bdellovibrio bacteriovorus and Host Bacteria[END_REF].

On the other hand, the decrease in HDNA and especially PLP was partly explained by the treatments, but it seemed that other factors were also at play, particularly for the increase in LDNA. The increase in LDNA was strongly and positively correlated with chlorophyll a and negatively correlated with temperature. These two factors seemed to play a role in LDNA dynamics, despite their variations being minimal. Furthermore, LDNA was negatively correlated with HDNA, implying that the dynamics of one probably influenced that of the other. For the decrease of PLP, especially HDNA, the environmental variables of the correlogram and RDA did not appear to show a clear pattern. Both were positively correlated and showed several other correlations with other factors. In contrast, LDNA had more negative relationships with the environmental variables.

Other interactions may explain the dynamics of Bdellovibrionaceae, PLP, and HDNA. Among the other parameters, we were able to study the dynamics of viruses (PLP). According to literature, some bacteriophages are capable of infecting BALOs [START_REF] Hashimoto | Isolation of a Bacteriophage for Bdellovibrio bacteriovorus[END_REF], and in Lake Geneva PLPs and HDNAs [START_REF] Personnic | Seasonal variations of microbial abundances and virus-versus flagellate-induced mortality of picoplankton in three peri-alpine lakes[END_REF][START_REF] Jacquet | Are small grazers and/or viruses a structuring factor of the free-living bacterial community in Lake Geneva?[END_REF]. PLP dynamics were observed to increase over time. Statistical tests did not show any treatment effect on PLP abundance. Despite the relative increase in VLPs, there was no correlation with BALOs or PLPs, HDNA, and LDNA. The presence and abundance of VLPs had no visible impact on the prokaryotic community and BALOs. In contrast, the virus-to-prokaryote ratio (VPR), which can be used as a proxy for host-parasite interactions [START_REF] Parikka | Deciphering the virus-to-prokaryote ratio (VPR): insights into virushost relationships in a variety of ecosystems[END_REF], showed the opposite trend. This ratio increased strongly in the H treatment, with statistical analysis showing a significant difference between the H and M treatments. This suggests that phages affected prokaryotic growth and dynamics in the H treatment through cell lysis. In contrast, the M and C treatments did not seem to favour a "boom" of phages.

We also noted a negative correlation between VPRs and Peredibacteraceae and Bacteriovoracaceae but no significant correlation with Bdellovibrionaceae. Finally, we observed that Peredibacteraceae and Bacteriovoracaceae appeared to share common trends, and these two families were positively correlated with each other, and Bdellovibrionaceae had a weaker positive correlation. The current study could not explain Bdellovibrionaceae dynamics.

We are aware that mesocosms can have a significant detrimental effect on planktonic populations but we do not think that Bacteriovoracaceae could be more susceptible to such effect, compared to the other predators. Firstly, we did not measure inhibitors that could have affected Bacteriovoracaceae DNA detection but also all the other BALOs' DNA. Secondly, from previous studies, we know that Bacteriovoracaceae are less abundant than Peredibacteraceae and Bdellovibrionaceae in the water column [START_REF] Paix | Diversity, dynamics and distribution of Bdellovibrio and like organisms in peri-alpine lakes[END_REF][START_REF] Ezzedine | New 16S rRNA primers to uncover Bdellovibrio and like organisms diversity and abundance[END_REF][START_REF] Ezzedine | Bdellovibrio and like organisms in Lake Geneva: an unseen elephant in the room?[END_REF]. It is possible that variations are likely to be more important and visible for such low abundances. We can wonder, however, why Bacteriovoracaceae are less numerous than Peredibacteraceae and Bdellovibrionaceae? As pointed out above, Bacteriovoracaceae may prefer other types of habitats such as biofilm, sediment, wastewater treatment plant effluent, etc.… Another nonexclusive possibility is that these bacteria are less competitive than the other two BALOs and/or have a limited prey range. Also, they could be preferably more predated or lysed than other BALOs.

We hypothesise that factors other than those measured here, such as the importance and role of other predators, e.g. heterotrophic nanoflagellates, ciliates, and metazooplankton, could be determinants of the observed effect. These groups could likely exert considerable control over bacteria, explaining the observed dynamics of both PLPs and BALOs. It would be interesting to specifically target BALO prey using an adequate set of qPCR primers instead of looking at the general trend of prokaryotes. However, we only identified one BALO prey species in Lake Geneva, i.e. Pseudomonas spp. [START_REF] Ezzedine | Bdellovibrio sp.: an important bacterial predator in Lake Geneva?[END_REF]. Further efforts to pinpoint prey strain and type are needed to fully understand BALO dynamics in Lake Geneva and peri-alpine lakes.

Conclusions

Bacteria play key roles in ecosystem functioning, as they are involved in nutrient cycling and are prey for a variety of predators, among other functions. Therefore, it is essential to accurately predict their relationships with other species and general patterns governing their adaptation to the environment. Here, we shed some light on a group of bacteria that are still poorly understood in terms of their ecology, particularly, in natural freshwater systems such as lakes. Although BALOs were not significantly affected by our treatment conditions, which were supposed to mirror the advent of extreme events on the surface waters of Lake Geneva, they were shown to be a dynamic community, and that constitutive families of this functional group exhibited different patterns. As the main environmental 1 3 physicochemical parameters measured in the present study did not seem responsible for the observed dynamics, the importance of biotic interactions is proposed. These interactions were likely predation and parasitism through the action of phages and protozoans, but also possibly the scarcity of specific prey for BALOs. Next steps should include more (biotic) factors, focus on the quantity and quality of prey, and examine different types of predation on both total and specific groups of bacteria.

Fig. 1

 1 Fig. 1 Dynamics of BALOs as measured by qPCR using adequate set of primers that amplify the 16S rRNA gene of each family (A-C) and heterotrophic prokaryotes as measured using Flow cytometry (D-F)

  Number of copies =DNA concentration (Insert size + Plasmid size) × 660 × 6.02 × 10 23 .Copy per milliliter =Copy per reaction Dilution factor for 25 ng μL -1 × DNA elution volume Filtered volume on 0.2 μm .

Fig. 2 Fig. 3 3

 233 Fig. 2 Box plot of Bdellovibrionaceae abundances measured by qPCR (A) and prokaryote-like particles (PLP) measured by Flow cytometry (B) for each treatment over time. The solid line in each box plot corresponds to the median value. The letters C, M, and H denote the treatment types, i.e. control, moderate, and high, respectively. The statistical tests used are displayed at the top of the figures. The significance of the tests is represented by stars and the non-significance by the letters "ns"

Fig. 4 2 -

 42 Fig. 4 This correlogram is based on 21 observations that represent environmental variables, PLPs, VLPs, and BALOs. Positive and negative correlations are displayed in blue and red, respectively. The intensity of the colour is proportional to the correlation coefficient. The white boxes indicate that the correlation is not significantly different from 0 at the 5% level tested with a correlation test. Abbreviation are as follows: ammonium (NH 4 + ), chloride (Cl -), chlorophyll a

Fig. 5

 5 Fig. 5 RDA representation indicating the relationship between environmental variables (arrows) and BALOs (A) or PLP (B) sampled for 6 days. 0 day, 4 days, 7 days, 12 days, 14 days, and 18 days correspond, respectively, to 4th, 8th, 11th, 16th, 18th, and 22nd of July. Treatments are represented with shapes i.e. circle for C, triangle for M, and square for H. The days are highlighted with colours. Abbreviations are as follows: ammonium (NH 4 + ), chloride (Cl -), chloro-

Acknowledgements

The authors would like to thank all those who participated in the MESOLAC project: Laurent Espinat, Clementine Gallot, Pascal Perney, Jean-Christophe Hustache, Philippe Quetin, Viet Tran-Khac, Laura Crepin, and Mathilde Chevallay. We also thank Cécilia Barouillet and Teofana Chonova for their statistical assistance. Finally, professional scuba divers (Jean-Marc Bel, Stéphanie Kocca, Stéphan Jacquet and Eric Mocelin) are gratefully acknowledged for their help in anchoring the eco-friendly (e.g. ecological mooring) mesocosms. We would like to thank Editage (www. edita ge. com) for English language editing.

Data Availability All data can be made available. It is noteworthy that two data papers have been published on the MESOLAC data (https:// doi. org/ 10. 1016/j. dib. 2020. 106255 and https:// doi. org/ 10. 1016/j. dib. 2021. 107150).

Supplementary Information

The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s00248-021-01875-9.

Author Contributions JE participated to data analysis and writing. AJ helped to analyse the data. SR and ID co-organised the MESOLAC project with SJ. SJ co-organised the MESOLAC project with SR and ID, planned this study, participated to writing, and submitted the manuscript.

Funding

Declarations

Conflict of interest The authors declare that they have no conflict of interest.

Ethical Approval Not applicable.