Large Diversity and Original Structures of Acyl-Homoserine Lactones in Strain MOLA 401, a Marine Rhodobacteraceae Bacterium

Quorum sensing (QS) is a density-dependent mechanism allowing bacteria to synchronize their physiological activities, mediated by a wide range of signaling molecules including N-acyl-homoserine lactones (AHLs). Production of AHL has been identified in various marine strains of Proteobacteria. However, the chemical diversity of these molecules still needs to be further explored. In this study, we examined the diversity of AHLs produced by strain MOLA 401, a marine Alphaproteobacterium that belongs to the ubiquitous Rhodobacteraceae family. We combined an original biosensors-based guided screening of extract microfractions with liquid chromatography coupled to mass spectrometry (MS), High Resolution MS/MS and Nuclear Magnetic Resonance. This approach revealed the unsuspected capacity of a single Rhodobacteraceae strain to synthesize 20 different compounds, which are most likely AHLs. Also, some of these AHLs possessed original features that have never been previously observed, including long (up to 19 carbons) and poly-hydroxylated acyl side chains, revealing new molecular adaptations of QS to planktonic life and a larger molecular diversity than expected of molecules involved in cell–cell signaling within a single strain.


INTRODUCTION
Quorum sensing (QS) allows bacteria to sense their population density (Nealson, 1977) and coordinate their gene expression (Bassler, 1999;Fuqua and Greenberg, 2002) and physiology (Miller and Bassler, 2001). QS communication is based on the secretion and detection of small molecules by bacteria in their nearby environment (Atkinson and Williams, 2009). A large number of studies have demonstrated that QS regulates many different bacterial features including biofilm production (Parsek and Greenberg, 2005;Dickschat, 2010), nodulation (Cha et al., 1998;Loh et al., 2002), bioluminescence (Waters and Bassler, 2005), virulence factor production (Smith and Iglewski, 2003) among others (Diggle et al., 2007). The coordination of bacterial community activities is supposed to confer an ecological advantage to the population (Case et al., 2008).
Among the various molecular signals used in QS systems, AHLs (acyl-homoserine lactone or autoinducer type-1, AI-1) constitute the major class of semiochemicals (Williams et al., 2007;Papenfort and Bassler, 2016) which has been widely studied (Fuqua et al., 1994). AHLs are homoserine lactone (HSL) linked to fatty acyl chains through an amide bond. The acyl chain length can vary from 4 to 18 carbons and sometimes includes a 3-oxo or a 3-hydroxy functional group (Fuqua and Greenberg, 2002). AHLs are usually saturated, but some unsaturated bonds in the fatty acyl chains are known. AHLs are synthetized by AHLsynthases, which catalyze the amide bond formation between the acyl chain carried by the ACP (acyl carrier protein) and the amine moiety precursor SAM (S-adenosyl-methionine) (Fuqua and Greenberg, 2002;Parsek and Greenberg, 2005). Three AHL synthase genes are currently known: ainS-like (Gilson et al., 1995), luxI-like (Engebrecht and Silverman, 1984), and hdtSlike (Loh et al., 2002). Of these three, ainS-like genes are found only in Vibrio, luxI-like genes are the most well-studied and are present in many Proteobacteria genomes (Gelencser et al., 2012) and little is known about hdtS-like genes (Doberva et al., 2015).
Many marine bacteria regulate some of their physiological traits using QS systems, among them the Rhodobacteraceae, a key bacterial family in marine environments that drives important biogeochemical reactions (Gram et al., 2002;Schaefer et al., 2002;Wagner-Döbler et al., 2005;Wagner-Döbler and Biebl, 2006). Rhodobacteraceae are abundant in the ocean and it has been demonstrated that 87% of completely sequenced genomes in this group encode luxI-like genes (Cude and Buchan, 2013;Zan et al., 2014). Known AHLproducing Rhodobacteraceae are diverse (Wagner-Döbler et al., 2005). Among them, Ruegeria spp, is found associated with sponges (Mohamed et al., 2008;Zan et al., 2012Zan et al., , 2015, Dinoroseobacter spp, lives in association with dinoflagellates , Sulfitobacter sp. is a diatom-associated bacteria. (Amin et al., 2015;Limardo and Worden, 2015), and Phaeobacter gallaeciensis proliferates in coastal waters (Wagner-Döbler et al., 2005;Berger et al., 2011). Interestingly, QS has been mainly identified in strains isolated from niches where bacteria can reach high concentrations (phycosphere, sponges tissues) (Rolland et al., 2016). More generally, it is commonly thought that QS is uncommon in marine oligotrophic strains as bacterial concentrations in oligotrophic environments (approximately 10 5 cells per mL) were presumably too low to trigger QS behaviors. However, a few recent publications report the occurrence of QS in marine bacterial strains isolated in oligotrophic environments. In their pioneering work, Moran et al. (2004) sequenced the full genome of Silicibacter pomeroyi, an oligotrophic Roseobacter and detected two QS systems. Similarly, in a previous study, we reported many translated gene sequences affiliated to Roseobacters luxI and hdtS from oligotrophic environments in the predicted proteome of Global Ocean Sampling metagenomic dataset (Doberva et al., 2015).
Collectively, these preliminary observations suggest that QS could constitute an important physiological trait of Rhodobacteraceae in all types of aquatic environments. New AHLs are regularly described in this group (Schaefer et al., 2008;Thiel et al., 2009;Ziesche et al., 2015). This suggests that the real extent of AHL chemical diversity in marine bacteria is still unknown. However, to our knowledge, no studies had yet used bioguided microfractionation combined with thorough mass spectrometry-based and Nuclear Magnetic resonance spectroscopy (NMR) for an in-depth description of AHLs diversity emitted marine Rhodobacteraceae. In this study, we investigated the potential of strain MOLA 401, isolated in an oligotrophic lagoon, to produce different types of AHLs. The strain MOLA 401 was isolated in a tropical oligotrophic lagoon located in New Caledonia. A closely related strain (Maribius pelagius B5-6 T ; 96% 16S rRNA sequence identity) has been isolated in the oligotrophic Sargasso Sea (Atlantic Ocean) (Choi et al., 2007). We had sequenced the full genome of the strain MOLA 401, and previously reported the presence of luxI, luxR, and hdtS genes, revealing the potential of this strain to communicate by QS (Doberva et al., 2014). We report here on the chemical diversity of strain MOLA 401 AHLs.

Culture of Strain MOLA 401
The strain MOLA 401 is from the MOLA culture collection (WDCM911 1 ) and available (strain code BBCC401) upon request 2 . This strain was isolated on December 3, 2004 at 4 m depth from marine oligotrophic waters in the southwest lagoon of New Caledonia (France; 22 • 21.23 S/166 • 23.43 E) (Conan et al., 2008). Sampled waters harbored a Chl a concentration of 1.07 µg ml −1 (F. Joux, pers. communication). All culturing steps were performed using Marine Broth (MB) 2216 (BD Difco, Sparks, MD, United States of America). The draft genome sequence has been published under accession number JQEY00000000 and revealed a quorum-sensing dependent physiology (Doberva et al., 2014).

Phylogenetic Analyses of the Strain MOLA 401 16S rRNA and Quorum Sensing Genes
Phylogenetic analyses were conducted using MEGA 6 software (Tamura et al., 2013) to assess the position of strain MOLA 401 QS genes with respect to already published phylogenies and to infer their putative functional role. Four different genes were considered: 16S rRNA, luxR (an AHL receptor), luxI (an AHL synthase), and hdtS (another AHL synthase). All alignments were performed using ClustalW (Larkin et al., 2007) and trimmed manually. Phylogenetic trees were constructed using the Neighbor-Joining (NJ) method with p-distance correction and 1000 bootstrap replicates.

AHL Standards Used in This Study
N-acyl-homoserine lactones were obtained from Cayman Chemical (Ann Arbor, MI, United States). Stock solutions (10 mmol L −1 ) of the analytes were prepared in dimethylsulfoxide (DMSO). Oxo-C14:1-HSL was dissolved in acetonitrile (CH 3 CN). The list of standards AHLs is provided in Supplementary Table S1.

Extraction of AHL from Strain MOLA 401 Supernatants
The strain MOLA401 was pre-cultured in 30 mL (96 h, 25 • C, 100 rpm) and then cultivated under aerobic conditions in 3 L of MB in 6 Erlenmeyer flasks under continuous shaking (200 rpm, 25 • C, 72 h, 5 mL from the preculture). At late exponential cell growth phase (72 h), tert-butyl methyl ether (600 ml) was added in each flask. This mixture was shaken overnight at room temperature (150 rpm). The two phases were then separated and the organic phase was dried with MgSO 4 . Solvent was filtered and removed with a rotary evaporator. The crude extract (170.6 mg) was dissolved in DMSO (57 mg mL −1 ). A Phenomenex Strata C18, 55 µm, 5 g column was equilibrated with CH 3 CN (100 mL), H 2 O:CH 3 CN 75:25 (35 mL), then H 2 O (100 mL). The crude extract dissolved in DMSO was deposited on top of the column. Elution was carried out with H 2 O (50 mL), CH 3 CN (50 mL, Fraction M). Fraction M was then evaporated in a rotary evaporator yielding 40.6 mg of material.

HPLC Fractionation
Fraction M was then dissolved in DMSO (40 mg mL −1 ) and fractionated on a preparative HPLC system with 2 Varian Prep Star pumps, a manual injector, a Dionex Ultimate 3000 RS variable wavelength detector and a Dionex Ultimate 3000 fraction collector. The column was a Phenomenex Luna C18, 5 µm, 21.2 × 250 mm, and the flow rate was set to 20 mL min −1 . The solvent was gradient grade H 2 O and CH 3 CN (70:30 for 3 min, followed by a 12 min linear gradient from 70:30 to 0:100, followed by 100% CH 3 CN for 10 min). The eluents were monitored at 214, 254, 274, and 280 nm, and were collected between 3 and 25 min (1 fraction min −1 , 22 fractions total referenced as M1-M22). The solvent was removed from each fraction with a genevac HT-4X system. Each fraction was dissolved in 100 µL DMSO to perform biosensor tests.

Culture of Biosensor Strains
The 22 fractions (M1-M22) previously prepared were tested in the AHL biosensor assay following previously described protocols using Pseudomonas putida and Escherichia coli based biosensors (Andersen et al., 2001;Riedel et al., 2001;Steindler and Venturi, 2007). Briefly, P. putida F117 (pRK-C12; Kmr; ppuI::npt) was used for the detection of long-chain AHLs (Andersen et al., 2001) and E. coli MT102 (pJBA132) for the detection of short chain AHLs (Riedel et al., 2001). E. coli MT102 and P. putida F117 were cultivated in Luria-Bertani (LB) Broth (Sigma L3022) overnight with continuous shaking (200 rpm), at 37 • C supplemented with tetracycline (25 µg mL −1 ) and at 30 • C supplemented with gentamicin (20 µg mL −1 ), respectively. An overnight culture of each biosensor strain (200 µL) was inoculated in 9.8 mL of fresh LB medium with the adapted antibiotics. This fresh biosensor culture was dispensed into 96-well microplates (180 µL per well). Then, the microfractions in DMSO (20 µL) were added in each well in triplicate. Microplates were incubated at 30 • C and 37 • C depending of growth optimum of the selected biosensor strain, without shaking. After 0, 5 and 24 h of incubation, fluorescence was determined with a Victor1420 Multilabel Counter (Perkin-Elmer) at an excitation wavelength of 485 nm and a detection wavelength of 535 nm. OD620 was also measured to control for biosensor cell growth. Negative controls were biosensor cultures without extract, and sterile LB medium. Biosensor cultures with addition of commercial AHLs (C6-HSL for E. coli MT102 and oxo-C10-HSL for P. putida F117) were used as a positive control.

LC-MS Analyses
UHPLC-MS analyses were performed with a Waters (Milford, CT, United States) Acquity UPLC-TQD (Triple Quadrupole Detector) system controlled by the MassLynx 4.1 software. Column was an Acquity HSS C18 (2.1 × 50 mm) with 1.8 µm particle size (Waters). The column oven was set to 40 • C. The flow rate was maintained at 0.6 mL min −1 and the injection volume was 2 µL. The mobile phase was composed of 0.1% formic acid in water (eluent A) and 0.1% formic acid in acetonitrile (B). A gradient profile was used, starting with 95% of A, keeping this composition constant for 0.5 min. Proportion of B was linearly increased to 100% in 6.5 min, and was left at 100% for 3 min.
The T.Q. Detector operated in ElectroSpray Ionization (ESI) in the positive and negative modes. First, the third quadrupole (Q3) has been used in scanning mode on the m/z 50-800 mass range in order to confirm the molecular weight and the purity of our 26 standard AHLs (2 mg mL −1 in DMSO, 2 µL injected), but also to determine their retention time (RT) under our chromatographic conditions. Two cone voltages (30 and 60 volts) were applied both in ESI + and ESI − modes. The other ion source parameters were as follows: capillary voltage 3.2 kV for positive mode (3 kV in negative mode), the source temperature was set at 150 • C and the desolvation temperature was 450 • C. Nitrogen was used as desolvation gas at a flow rate of 800 L h −1 and as cone gas at a flow rate of 50 L h −1 . The analytical approach first involved the study of mass spectra obtained for our standard molecules. These compounds ionized significantly better in the positive mode and the signal of the protonated molecule ([M+H] + ) appeared more abundant when the cone voltage involved was lower. A peak corresponding to the cationized AHL with ubiquitous sodium ([M+Na] + ) often had a significant intensity too. Applying a higher cone voltage led to fragmentation in the ion source. In particular, a fragment ion at m/z 102 was specific for the HSL moiety. This signal was chosen as the specific ion indicating the presence of HSL-type compounds. In a second step, the first quadrupole (Q1) of the TQD instrument was used in scanning mode from m/z 50-500 as mass range and several cone voltages (10, 15, 20, and 25 volts) were applied in order to determine the best value to observe the more intense [M+H] + signal for each standard AHL. The [M+H] + ions were later used as the precursor ion for MS/MS experiments. Each ion of interest was selected by the first quadrupole (Q1) and then focused in the collision cell (Q2) where fragmentation reactions occurred. The resulting fragment ions were finally analyzed by the third quadrupole (Q3). The collision gas (argon) was introduced into the collision cell to maintain a pressure near to 4.5 × 10 −3 mbar. The collision energy was optimized to lead to an attenuation of the precursor ion beam of almost 85%. The fragmentation pattern of each [M+H] + standard ion (MS/MS spectrum) has been recorded with the most suitable parameters for a later comparison with those obtained for the signals of interest observed in samples.

Molecular Formula Determination and High Resolution MS/MS
High-resolution MS/MS analyses were conducted with a Thermo UHPLC-HRMS system. Analyses of microfractions and standards (1.0 µL injected) were performed in electrospray positive ionization mode in the 133.4-2000 Da range in centroid mode. The mass detector was an Orbitrap MS/MS FT Q-Exactive focus mass spectrometer. The analysis was conducted in FullMS data dependent MS2 mode. In FullMS, resolution was set to 70,000 and AGC target was 3.10 6 . In MS2, resolution was 17,500, AGC target 10 5 , isolation window 0.4 Da, normalized collision energy 30, with 15 s dynamic exclusion. UHPLC column was a Phenomenex Luna Omega polar C-18 150 × 2.1 mm, 1.6 µm. The column temperature was set to 42 • C, and the flow rate was 0.5 mL min −1 . The solvent system was a mixture of water (A) with increasing proportions of acetonitrile (B), both solvents modified with 0.1% formic acid. The gradient was as follows: 2% B 3 min before injection, then from 1 to 13 min, a shark fin gradient increase of B up to 100% (curve 2), followed by 100% B for 5 min. The flow was discarded (not injected into the mass spectrometer) before injection and up to 1 min after injection. The exact masses and corresponding molecular formulas are reported in Table 2. A full list of standards along with RTs and exact masses is provided in Supplementary Table S1.

Molecular Networking
A molecular network was constructed based on UHPLC-HRMS/MS analyses using the GNPS platform 3 . Nodes from MOLA 401 microfractions are in yellow, those from standards appear in blue, and those detected in both are in green. The number of compared ions was set to 8, and the minimum cosine for linking two parent ions was set to 0.7. With these parameters, only short side chains AHLs were not clustered. Since the strain microfractions contained so many AHLs, the detection limit was set to 1000 in order to simplify the cluster.

NMR Analyses
Nuclear magnetic resonance spectra were recorded in DMSO-d 6 on a Bruker 600 MHz NMR spectrometer equipped with a 1 mm

RESULTS AND DISCUSSION
Presence and Chemical Features of AHLs in Strain MOLA 401 MOLA 401 culture supernatant was extracted with tert-butyl methyl ether and fractionated in 22 fractions, M1-M22. These fractions were tested for AHL production using the biosensor strains E. coli MT102 and P. putida F117, which are GFP-based biosensors emitting light in presence of AHLs. Interestingly, 10 fractions were positive in these assays with both biosensors (M9, M10, M11, M12, M13, M14, M15, M16, M17, M18) (Supplementary Figures S1, S2). For detection of AHLs, we initially focused on LC-MS profiling with single ion recording at m/z 102, which corresponds to the mass of the protonated homoserine moiety. MS ionization conditions were optimized in order to favor the formation of this fragment. Then the full MS scan at the RTs pointed out in SIR102 allowed us to propose a list of pseudomolecular ion masses of putative AHLs. UHPLC coupled to high resolution MS/MS analyses were then conducted in the discovery mode. This allowed us to calculate the molecular formulas of the putative AHLs based high resolution masses, to obtain high resolution fragmentation analyses, and to obtain a molecular network including all the microfractions along with the 26 AHLs standards (Figure 1; Wang et al., 2016). All MS spectrum, SIR102 chromatograms, TIC chromatograms, MS spectrum, High Resolution MS spectrum are provided in Supplementary  Figures S4-S61.
Overall, it was demonstrated that MOLA 401 produced at least 20 different AHLs out of the 21 putative ones detected by SIR102 ( Table 1). The confirmation of the presence of the HSL subunit was obtained by the method described by Patel et al (Patel et al., 2016). In our case, all the AHLs had all 4 diagnostic fragments at m/z 102.055, 84.045, 74.061, and 56.050 in MS/MS. Also, compound Q was identical to standard 23 [N-hexadec-11(Z)enoyl-L-homoserine lactone]. The molecular network shown in Figure 1 further demonstrated that all MOLA 401 compounds identified in this study as potential AHLs clustered with the network defined by the standards. Much to our surprise, the molecular networking analysis uncovered AHLs in the strain although the detection limit was set to 1000, making the cluster much simpler (note that only B was not detected with this value). For the present article, restrained to the ones found in the SIR102 analysis, but many AHLs detected in the cluster did in fact present HSL diagnostic ions.
Eventually, it turned out that microfraction M17 essentially contained AHLs (P, Q, R). This fraction was analyzed by 1D and 2D NMR ( Table 2). The 1 H NMR spectrum showed the presence of a methylene at δ H 2.36 (m, 1H, 4a) and δ H 2.11 (m, 1H, 4b), an oxomethylene at δ H 4.33 (td, J = 8.9, 1.8, 1H, 5b) et δ H 4.20 (m, 1H, 5a), and a methyne at δ H 4.56 (m, 1H, 3). Long-range 1 H-13 C correlations between H-3/H-5a and carbonyle C-2 at δ C 175.3, as well as the sequence of FIGURE 1 | Molecular network containing the AHLs detected in strain MOLA401 microfractions (yellow), and in analytical standards (blue). Identical matches are in green. In the graph, the parent mass is reported on the node, and for each node the corresponding COZY correlations between H-3, H-4, and H-5 confirmed the presence of a lactone ring (Figure 2). Then, the 1 H-13 C HMBC correlations of amide proton at δ H 8.29 (td, J = 8.9, 1.8, 1H) with carbons C-3 (δ C 47.4) and C-2 (δ C 170.8) allowed us to position an acylamino group in C-3, therefore confirming that the 3 major compounds of fraction M17 were AHLs. For compounds with a hydroxyl group on the side chain (P, R), it was possible to ascertain the CH 2 CH(OH)CH 2 partial sequence based on 1 H-1 H correlations of the oxomethine at δ H 3.78 (m, 1H, 4 ) with the methylenes at δ H 2.19 (m, 2H, 3 ) and at δ H 1.37 (m, 1H, 5 a) / δ H 1.30 (m, 1H, 5 b), and based on the long-range 1 H-13 C correlation of H-3 with C-2 . The rest of the side chain cannot be attributed due to extensive overlapping of the signals. Nevertheless, vinyl protons give key information on the double bonds in P and Q. The attribution of protons and carbons a-d in fragment B based on HSQC and COZY experiments was straightforward (Figure 2). Fragment B was constituted of two methylenes at δ H 1.98 (m, 4H, a and d) and two vinyl protons at δ H 5.32 (m, 2H, b and d). The shape of the vinyl protons signal confirmed the Z configuration of the double bond in P and Q (Frost and Gunstone, 1975). All NMR spectrum are provided in Supplementary Figures  S62-S66.
The length of acyl chains in the detected AHLs ranged between 15 and 19 carbons. To our knowledge, this also the first report of acyl chains longer than 18 carbons. Also, 5 AHLs presented an odd number of carbons in their acyl side chain ( Table 1). This observation also constitutes an interesting feature, as very few AHLs with acyl side chain presenting an odd number of carbons have been previously identified. More frequently, such AHLs were present as trace elements (C13:0-HSL, C15:0-HSL, C15:1-HSL, C15:2-HSL) (Wagner-Döbler et al., 2005), except in Sulfitobacter sp. D13 where the 9-C17:1-HSL is an AHL which appears synthesized in large quantities (Ziesche et al., 2015).
We also detected at least 6 AHLs with two or three hydroxyl groups along the acyl side chain. An examination of previously characterized AHLs revealed only single hydroxylation per acyl chain (Churchill and Chen, 2011) located at C-3. This is the case for the AHL detected in the marine Roseobacter strains Phaeobacter sp. BS107 or Loktanella sp. F14 who produces 3-OH-C12:1-HSL (Ziesche et al., 2015). Thus, we report here another interesting new feature of marine AHLs, which is the existence of poly-hydroxylation of the acyl chain ( Table 1). The position of the hydroxyl groups along the acyl chain could not be determined as these groups did not induce fragmentation of the side chain in MS/MS. NMR of the microfractions were very difficult to interpret due to the relatively low proportion of each AHL in these fractions. However, despite these limitations, our data unambiguously indicate that the strain MOLA 401 is able to synthesize a wide diversity of AHLs. Also, we detected at least 2 AHLs presenting one double bond in their acyl side chain ( Table 1). The position and configuration of the double bound chain was confirmed in compound Q by the analytical standard 23. When there was oxygen and double bonds detected in the side chain, it was not possible to distinguish a carbonyl group or a hydroxyl and a carbon-carbon double bond, as the two would lead to the same molecular formula. Short acyl chain molecules are more polar and soluble in seawater than those presenting long aliphatic chains, which are thus less hydrophilic. However, it appears that marine bacteria produce AHLs with long chains (Wagner-Döbler et al., 2005;Zan et al., 2012). Thus, our data confirm these previous observations. Also, our technical approach revealed that many AHL acyl chains were oxidized. Such observation indicates that these AHLs are adapted for signal release and diffusion in marine environments as acyl side chain modifications would increase water solubility and compatibility with active efflux pumps (Pearson et al., 1999). Most of Rhodobacteraceae bacteria produce long chain AHLs with additional modifications (Cude and Buchan, 2013). For example, the marine free-living strain Rhodobacter sphaeroides produces C14:1-HSL (Puskas et al., 1997), the marine dinoflagellate associated bacterium Dinoroseobacter shibae synthetizes mainly C18:2-HSL and C18:1-HSL, but also traces of C16-HSL, C15-HSL and C14-HSL (Wagner-Döbler et al., 2005;Neumann et al., 2013;Patzelt et al., 2013), the sponge symbiont Ruegeria sp. emits OH-C14-HSL, OH-C14:1-HSL and OH-C12-HSL (Zan et al., 2012). S. pomeroyi produces the p-coumaroyl-HSL, a non-conventional AHL in which the acyl side chain is replaced by a coumaroyl moiety (Schaefer et al., 2008). Nevertheless, the poly-hydroxylation of acyl chain observed in strain MOLA 401 combined with the presence of unsaturation appears to be an original feature. We hypothesize that the AHL synthase produces a molecule with acyl chain containing 15 to 19 carbons, and that additional modifications of the acyl chain are mediated by cytochrome P450 (Chowdhary et al., 2007) (WP_036181863.1), which oxidizes aliphatic chains, and by desaturases which produce double bonds (Aguilar and de Mendoza, 2006). Interestingly, we detected a cytochrome P450 homolog in the genome of the strain MOLA 401 (Doberva et al., 2014).

Linking Genetic and Chemical Features
Phylogenetic analyses based on 16S rRNA and putative LuxI sequences confirm the position of strain MOLA 401 in the Rhodobacteraceae family and the Proteobacteria phylum ( Figure 3A). The position of this strain, close to two Maribius isolates was well supported (BP NJ = 100) ( Figure 3A) and confirmed affiliation to the Rhodobacteraceae family. The strain MOLA 401 putative LuxI protein sequence is closely related to other LuxI sequences of Rhodobacteraceae strains within Alphaproteobacteria (Figure 3B). Clustering of the Rhodobacteraceae LuxI sequences (group 1 includes Ruegeria pomeroyi, Roseobacter denitrificans, P. inhibens; group 2 includes D. shibae, Maribius sp., Jannaschia sp.) were well supported ( Figure 3B). Similarly, the phylogenetic tree based on the AHL receptor LuxR placed the strain MOLA 401 putative LuxR within the Rhodobacteraceae (Figure 3C). These data clearly indicate that the strain MOLA 401 strain belongs to the Roseobacter group with respect to its 16S rRNA or the genes encoding for AHL production and reception. This makes strain MOLA 401 an ideal model strain for future studies of QS in marine environments. Also these data confirmed previous observation based only on 16Sr RNA genes (Choi et al., 2007).
Another protein potentially involved in AHL production is HdtS, of which two homologs have been detected in the full genome sequence of the strain MOLA 401 (Doberva et al., 2014). HdtS is a member of the lysophosphatidic acid acyltransferase family (Laue et al., 2000) and has a dual functionality, acylation of lysophosphatidic acid (Cullinane et al., 2005) and AHL synthesis (Laue et al., 2000). The HdtS-mediated production of AHL has been demonstrated experimentally in P. fluorescens (Laue et al., 2000) and Acidithiobacillus ferrooxidans (Rivas et al., 2007). P. fluorescens produces 3-OH-C14:1-HSL, C10-HSL and C6-HSL, while A. ferrooxidans produces a C14-HSL. The strain MOLA 401 putative HdtS sequences clustered into two groups, both of which were related to putative HdtS found in other Rhodobacteraceae, with strong bootstrap supports ( Figure 3D). Interestingly, one homolog was clustered with HdtS from the Gammaproteobacteria A. ferrooxidans and P. fluorescens, the only HdtS enzymes with confirmed AHL synthesis activity (acyltransferase2 sequences, Figure 3D). This suggests that the MOLA 401 putative HdtS is similarly contributing to the AHL pool produced by MOLA 401, in cooperation with putative LuxI. However, an experimental confirmation of such HdtS based AHL production in strain MOLA 401 is required in future studies, also because the MOLA 401 strain does not produce similar AHL as those found in A. ferrooxidans and P. fluorescens.
The specificity of LuxI synthases varies, especially in regards to the type of acyl side chain recognized as substrate (Gould et al., 2004). For example, the LasI (a LuxI homolog) in P. aeruginosa produces different AHLs depending on the growth conditions and the host. By contrast, YspI and EsaI, respectively, found in Yersinia pestis and Erwinia stewartii, are specific to one type of acyl-ACP (Gould et al., 2006) producing defined AHLs. The strain MOLA401 putative LuxI synthase has the conserved the arginine and phenylalanine in positions 25 and 29 (two key aminoacids residues in this protein), respectively, similar to the P. aeruginosa LuxI (Supplementary Figure S3). Thus MOLA 401 putative LuxI sequence is consistent with a capacity to produce a large number of AHLs. One possible hypothesis is that the same LuxI-synthase may produce several AHLs with low side chain length specificity, as demonstrated by Neumann et al. (2013).

Culture of Strain MOLA 401 and QS Abilities
The strain MOLA 401 is a bacterium from the Rhodobacteraceae family isolated in an oligotrophic lagoon. Phylogenetically close Maribius strains have also been isolated in such oligotrophic waters, like in the Sargasso Sea (Choi et al., 2007). The ability of bacteria isolated from oligotrophic waters to communicate could appear paradoxical (Moran et al., 2004) as cell densities in such environments are below the expected threshold that enables QS. However, our study demonstrates the ability of the strain MOLA 401 to synthesize diverse types of AHLs. We experimented on a MOLA 401 strain cultured under rich nutrient conditions (Marine Broth). Thus, we could suggest that the large spectrum of AHL produced by MOLA 401 might give this strain the ability to exploit organic matter by a complex coordination of the bacterial population (Rolland et al., 2016). This observation is in line with previous hypothesis suggesting that such coordination allows particleattached bacteria to exploit marine organic matter (Moran et al., 2016). Future studies need to be conducted to evaluate the capacity of Rhodobacteraceae to produce AHLs when cultured in oligotrophic media.
Collectively, our technical approach based on a bioguided search of AHL in bacterial extracts and the obtained data reveal that the Rhodobacteraceae strain MOLA 401 isolated in an oligotrophic lagoon is able to produce a very large number of different AHLs. The AHLs characterized in this study possessed interesting and original features including variable acyl chain length and multiple-hydroxylation sites. The strain MOLA 401 strain provides new insights into the breadth of possible AHL diversity, suggesting the existence of original adaptations of bacterial dialogs to marine environments.

AUTHOR CONTRIBUTIONS
MD, DS, JS, NH, SS-F, VE, YF, and RL conducted the experimental work. MD, DS, PL, SS-F, and RL designed the experiments. All authors wrote the manuscript.