domain: sdv sde

text:
P Ory

S Palesse

D Delmas

H Montanié

In situ structuring of virioplankton through bacterial exoenzymatic activity: interaction with phytoplankton

The abundance and composition of microbial communities were investigated during 2007 in Marennes Oléron Bay (France) in order to characterize the biological relationships governing the planktonic food web. Compared to the results of previous years, there was a lower abundance and a lower magnitude of variations in autotrophic and heterotrophic microbial entities in 2007. Using a Spearman rank correlation analysis, a significant link for biological abundance was only noted between bacteria and ciliates.

However, an interesting relationship was highlighted between chlorophyll a and bacterial aminopeptidase and -glucosidase exoenzymatic activities. Two characteristic periods of virus-bacteria uncoupling were observed in spring and autumn 2007. Using multivariate analysis, the clustering of monthly viral community structures (defined by RAPDfingerprinting) was related to bacterial proteolysis activity and secondarily to flagellates.

Spring and autumn were characterized by phytoplanktonic blooms of large (>10 µm) and small cells (<10 µm), respectively, during the virus-bacteria uncoupling phases. High predation by flagellates in the autumn, during a period of microbial food web activity, involved the highest viral fingerprint richness, inferring a top-down control in driving the number of virus-host systems. We discuss the proteolysis level as a descriptor in structuring viral assemblages and argue that changes in the quantity and quality of DOM may bridge phytoplankton to bacteria and reflect changes in phytoplankton composition.

INTRODUCTION

Bacteria and viruses are the most abundant biological entities in the marine microbial world [START_REF] Borsheim | Enumeration and biomass estimation of planktonic bacteria and viruses by transmission electron microscopy[END_REF][START_REF] Hennes | Direct counts of viruses in natural waters and laboratory cultures by epifluorescence microscopy[END_REF]. Producers such as phytoplankton and bacterioplankton are involved in trade-offs between (1) bottom-up and (2) top-down mechanisms regulating both abundance and community structure (Winter et al. 2010).

(1) In terms of abundance and community structure, bacterioplankton are particularly dependent on the quantity and quality of inorganic nutrients and organic matter, which may impact their metabolism, growth and community structure [START_REF] Kirchman | Changes in bacterial activity and community structure in response to dissolved organic matter in the Hudson River, New York[END_REF][START_REF] Teira | Effects of inorganic and organic nutrients inputs on bacterioplankton community composition along a latitudinal transect in the Atlantic Ocean[END_REF]. Around 30% of organic matter in the ocean is composed of high molecular weight compounds (HMWC, >600 Da), and their uptake by bacteria is only possible after extracellular hydrolysis [START_REF] Karner | Extracellular enzymatic activity and secondary production in free-living and marine-snow-associated bacteria[END_REF][START_REF] Nagata | Organic matter-bacteria interactions in seawater[END_REF]. The major components of this organic matter being proteins and carbohydrates, the ability of a species to produce peptidase and glucosidase enzymatic activities influences its growth and acts as a structuring factor of the community [START_REF] Van Wambeke | Short scale variations in nutrients, ectoenzymatic activities and bottom-up effects on bacterial production and community structure during late summer-autumn transition in the open NW Mediterranean Sea[END_REF]).

(2) The development of bacterial populations is also controlled by predators, namely grazers such as unpigmented nanoflagellates and parasites such as viruses (Winter et al. 2010). With regards to their abundance and grazing efficiency, unpigmented nanoflagellates are able to sustain their growth with bacteria and can potentially control bacterial abundance [START_REF] Berglund | Relative strength of resource and predation limitation of heterotrophic nanoflagellates in a low-productive sea area[END_REF]. The virus-mediated mortality of bacteria is comparable to that of grazer-mediated mortality [START_REF] Fuhrman | Viruses and protists cause similar bacterial mortality in coastal seawater[END_REF][START_REF] Fuhrman | Marine viruses and their biogeochemical and ecological effects[END_REF]. However, while grazing leads to carbon transfer to higher trophic levels, except with sloppy feeding [START_REF] Cherrier | Utilization and turnover of labile dissolved organic matter by bacterial heterotrophs in eastern North Pacific surface waters[END_REF], viral lysis enhances the microbial loop through the viral shunt [START_REF] Wilhelm | Viruses and nutrient cycles in the sea[END_REF], where the release of dissolved organic matter (DOM) which is available for bacteria favors the recycling of organic and inorganic nutrients [START_REF] Suttle | Marine viruses-major players in the global ecosystem[END_REF]. Indeed, without the microbial loop activities of prokaryotes, which transform nutrients into biomass, DOM would not be transferred to higher trophic levels. Thus, the viral shunt may act as a lubricant for the microbial loop short term and could contribute actively to biochemical cycles long term [START_REF] Middelboe | Regeneration of dissolved organic matter by viral lysis in marine microbial communities[END_REF][START_REF] Winget | Tangential flow diafiltration: an improved technique for estimation of virioplankton production[END_REF], Weinbauer et al. 2009).

From a qualitative point of view, prey-predator or parasite-host trophic interactions between microbial entities (prey-predator or parasite-host; [START_REF] Miki | Complex interactions in the microbial world: underexplored key links between viruses, bacteria and protozoan grazers in aquatic environments[END_REF] induce changes in the structure of each microbial assemblage: the size-selective grazing of nanoflagellates enhances the development of either small or large bacteria [START_REF] Jürgens | Predation as a shaping force for the phenotypic and genotypic composition of planktonic bacteria[END_REF] and impacts the morphometric (occurrence of resistant morphotypes) and taxonomic structures of bacterial communities [START_REF] Suzuki | Effect of protistan bacterivory on coastal bacterioplankton diversity[END_REF][START_REF] Langenheder | Regulation of bacterial biomass and community structure by metazoan and protozoan predation[END_REF]. Viruses, which have a higher specificity for their host, induce modifications in the composition of the bacterial community as described in the 'killing the winner' theory [START_REF] Thingstad | Elements of a theory for the mechanisms controlling abundance, diversity, and biogeochemical role of lytic bacterial viruses in aquatic systems[END_REF], which has been confirmed with experimental investigations and field data [START_REF] Fuhrman | Viral influence on aquatic bacterial communities[END_REF][START_REF] Hewson | Bacterial diversity in shallow oligotrophic marine benthos and overlying waters: effects of virus infection, containment, and nutrient enrichment[END_REF], Weinbauer & Rassoulzadegan 2003[START_REF] Schwalbach | Viral effects on bacterial community composition in marine plankton microcosms[END_REF], Winter et al. 2004b[START_REF] Bouvier | Key role of selective viral-induced mortality in determining marine bacterial community composition[END_REF], as reviewed and discussed in Winter et al. (2010). Hence, viral lysis of the abundant host group would favor rare groups and help maintain bacterial specific richness.

Viral lysis can also indirectly affect bacterial community composition: lysis, releasing new substrates, modifies DOM composition and enhances, via feedback, the growth of competitive bacterial species [START_REF] Middelboe | Virus-induced transfer of organic carbon between marine bacteria in a model community[END_REF]. Bacterial diversity and viral abundance are either positively or negatively linked (Winter et al. 2004a, Auguet et al. 2008). The presence of grazers can also stimulate viral multiplication in the bacterioplankton through a cascading effect from grazing-mediated resource enrichment (Šimek et al. 1997[START_REF] Jürgens | Morphological and compositional changes in a planktonic bacterial community in response to enhanced protozoan grazing[END_REF][START_REF] Ram | Functional responses of prokaryotes and viruses to grazer effects and nutrient additions in freshwater microcosms[END_REF].

Bottom-up and top-down controls act together in shaping virus-host system (VHS) diversity. As a result of trade-offs between competition and defense specialists (Winter et al. 2010), the number of simultaneously co-existing VHS are selected by top-down mechanisms (grazers, viruses), being more or less based on host-specificity, while bottom-up mechanisms (resources) drive the community structure of the VHS [START_REF] Sandaa | Viral control of bacterial biodiversity -evidence from a nutrient-enriched marine mesocosm experiment[END_REF]. Indeed, changes in nutritive conditions may affect the selection for dominant members of the bacterial community, and thus viral and bacterial community structures could change in parallel. First, bacterial phylotypes that are good competitors for nutrients would be lysed and then the virusresistant phylotypes would compete with each other. However, top-down and bottom-up controls are linked because viruses and grazers could also have a strong indirect bottom-up effect by modifying the quality and quantity of nutritive resources [START_REF] Sandaa | Viral control of bacterial biodiversity -evidence from a nutrient-enriched marine mesocosm experiment[END_REF][START_REF] Miki | Indirect interactions in the microbial world: specificities and similarities to plant-insect systems[END_REF].

However, aside from bacterioplankton, virioplankton and grazers, phytoplankton act not only as a provider of DOM but also as a competitor for inorganic nutrients [START_REF] Kirchman | The uptake of inorganic nutrients by heterotrophic bacteria[END_REF][START_REF] Rivkin | Inorganic nutrient limitation of oceanic bacterioplankton[END_REF][START_REF] Joint | Competition for inorganic nutrients between phytoplankton and bacterioplankton in nutrient manipulated mesocosms[END_REF]. Consequently, phytoplankton drive the development and structuring of other microbial entities and their interactions. Their link with viruses through viral infection has been clearly identified, particularly during bloom periods [START_REF] Bratbak | Viral mortality of the marine alga Emiliania huxleyi (Haptophyceae) and termination of algal blooms[END_REF][START_REF] Brussaard | Effects of grazing, sedimentation and phytoplankton cell lysis on the structure of a coastal pelagic food web[END_REF][START_REF] Agusti | Dissolved esterase activity as a tracer of phytoplankton lysis: evidence of high phytoplankton lysis rates in the northwestern Mediterranean[END_REF], with from 10 to 50% of algal cells visibly infected at the end of a specific bloom [START_REF] Gastrich | Viruses as potential regulators of regional brown tide blooms caused by the alga, Aureococcus anophagefferens[END_REF]). The release of phytoplankton cell content by viral lysis would provide enough DOM to meet the carbon demand of bacteria in the case of a phytoplankton decline (59 to 70%; [START_REF] Brussaard | A mesocosm study of phaeocystis globosa (Prymnesiophyceae) population dynamics -II. Significance for the microbial community[END_REF].

The relationship between phytoplankton and bacterial biomass or production has not been clearly established in field studies, probably due to the uncoupling of production periods relying on: (1) a consumption of allochthonous dissolved organic matter by bacteria, (2) a time lag between DOM exudation by phytoplankton and its subsequent consumption by bacteria, and (3) a nutrient limitation of consumption [START_REF] Alonso-Saez | Factors controlling the year-round variablity in the carbon flux through bacteria in a coastal marine system[END_REF]).

Exoenzymatic activities have been hypothesized to be indicators of the bacterial consumption of HMWC and thus could also characterize peaks of primary production. This is supported by the following arguments: (1) the proteolysis:glycolysis ratio is higher during productive periods, with low specific -glucosidase activity; (2) when phytoplanktonic cells collapse, the large amount of polysaccharides released induces the synthesis of -glucosidase;

(3) during periods of nutrient limitation, more HMWC and DOM may be released [START_REF] Alonso-Saez | Factors controlling the year-round variablity in the carbon flux through bacteria in a coastal marine system[END_REF]. However, exoenzymatic activities may also depend upon an influx of allochthonous C compounds (from river discharge or the resuspension of benthic matter).

These could be induced by the release of viral lysis products from bacteria, and activity changes could reflect shifts in bacterial community composition [START_REF] Martinez | Variability in ectohydrolytic enzyme activities of pelagic marine bacteria and its significance for substrate processing in the sea[END_REF].

In Marennes Oléron Bay (France), previous studies on pelagic microbial organisms and the control of their dynamics highlighted a link between virus abundance and chlorophyll a (chl a) concentration [START_REF] Ory | Pelagic food web patterns: do they modulate virus and nanoflagellate effects on picoplankton during the phytoplankton spring bloom ?[END_REF]) and described chl a biomass as a significant factor driving the variability of bacterial activities (secondary production and exoenzymatic activities; [START_REF] Auguet | Dynamic of virioplankton abundance and its environmental control in the Charente estuary (France)[END_REF] Sampling strategy and field measurements. Samplings were performed monthly from January to December 2007 at high tide, during periods of medium tidal coefficient (60) at 1 northern station located in the external part of the Charente River estuary (45°59'16.7''N, 1°10'0.7''W, 'Stn E' previously monitored by [START_REF] Auguet | Dynamic of virioplankton abundance and its environmental control in the Charente estuary (France)[END_REF] and by [START_REF] Ory | Pelagic food web patterns: do they modulate virus and nanoflagellate effects on picoplankton during the phytoplankton spring bloom ?[END_REF]); in February, sampling was performed closer to the mouth of the estuary (45°59'5.06''N, 1°06'32''W), because bad weather prohibited sampling at Stn E. Samples were collected near the surface, at <1 m depth. Subsamples for microorganism counts (except ciliates) were fixed on board according to the specificity of each organism. The remainder was kept in 8 l polypropylene dark bottles for <2 h at the in situ temperature until laboratory analyses were performed. Physical parameters (temperature, salinity and pH) were recorded on board (YSI 6600EDS-M multiparameter probe).

MATERIALS AND METHODS

Study

Zooplankton were collected using a 200 m WP-2 net equipped with a digital flowmeter (Model 2030, General Oceanics). The catch was fixed in formaldehyde (4% final concentration).

Physical and chemical parameters. Mineral and organic seston were assessed by filtering 200 to 500 ml of water, depending on the turbidity, under gentle vacuum pressure (GF/F Whatman filters; [START_REF] Auguet | Dynamic of virioplankton abundance and its environmental control in the Charente estuary (France)[END_REF].

Dissolved inorganic nutrients (phosphate [PO

4 ], urea, ammonium [NH 4 ], nitrite [NO 2 ]
and nitrate [NO 3 ]) were filtered (20 ml, GF/F Whatman) into a glass flask and stored for 2 to 3 wk at -20°C. Analyses were performed using an autoanalyser (Skalar; [START_REF] Strickland | A practical handbook of seawater analysis[END_REF].

Biological parameters. Chlorophyll a and pheopigment: Water samples of 200 ml were vacuum filtered onto 25 mm Whatman GF/F filters. The filters were stored for <1 mo at -20°C until assessment by extraction with 90% acetone (overnight at 4°C, under agitation) and fluorimetric analysis (Turner, TD 700) [START_REF] Lorenzen | A method for the continuous measurement of in vivo chlorophyll concentration Deep[END_REF].

Microbiological abundance:

Triplicate subsamples for viral and bacterial abundance assessments were fixed on board with 0.02 m filtered formaldehyde (2% final concentration) and stored for <24 h at 4°C. Samples were enumerated by epifluorescence microscopy [START_REF] Noble | Use of SYBR Green I for rapid epifluorescence counts of marine viruses and bacteria[END_REF] as follows: 2 ml of MilliQ water-diluted samples (50-fold on average for accurate microscopic observations) were retained on 0.02 m silicate membranes (Anodisc) followed by SYBR Green I staining (Invitrogen-Molecular Probes) for 30 min, then mounted on a glass slide and stored at -20°C for from 1 to 2 d until bacteria and viruses could be counted under blue excitation at 1000× magnification on a Zeiss Axioskop 2 Mot Plus microscope (100× Plan APO oil objective).

Small phytoplanktonic cells were quantified by flow cytometry using a FACSCalibur flow cytometer (size < 10 m) from samples flash-freezed and stored at -80°C (for <1 yr).

Synechococcus were distinguished from photosynthetic eukaryotes by their orange fluorescence; picoeukaryotes and nanoeukaryotes produced a red autofluorescence [START_REF] Joux | Distribution of picophytoplankton and bacterioplankton along a river plume gradient in the Mediterranean Sea[END_REF].

For nanoflagellates, duplicate samples of seawater were fixed on board (1% paraformaldehyde final concentration) and stored at 4°C for <1 wk. Samples were filtered through 0.8 µm black polycarbonate membranes (Nuclepore); the cells were stained with DAPI [START_REF] Porter | The use of DAPI for identifying and counting aquatic microflora[END_REF]) and frozen at -20°C for from 1 to 2 wk until counting under ultraviolet excitation. Counting was performed at 1000× magnification on a Zeiss Axioskop 2 Mot Plus microscope (100× Plan APO oil objective).

For ciliates, 1 l seawater was sampled and stored for 6 h in the dark at the in situ temperature in order to decant the suspended matter. Thereafter, the top 850 ml were siphoned off, gently mixed, and a subsample was fixed in alkaline lugol (1% final concentration) and stored at 4°C for <1 mo. Ciliates were counted in Utermöhl settling chambers, after decantation of from 20 to 50 ml, under white light at 630× magnification on a Leica DMIRB microscope (630× Fluotar objective). A sample of microphytoplankton was directly fixed with a combination of formaldehyde (1% final concentration) and alkalin lugol (1% final concentration), and from 10 to 20 ml subsamples were directly decanted in Utermöhl settling chambers before counting at 400× magnification on a Nikon diaphot microscope (200× APO Plan objective). Mesozooplankton samples were analyzed with a binocular loupe [START_REF] Sautour | Distribution of zooplankton populations in Marennes-Oléron Bay (France), structure and grazing impact of copepod communities[END_REF].

Bacterial activity: Leucine aminopeptidase (AMPase) and -glucosidase (-GLCase) activities were measured using L-leucine-7-amino-4-methylcoumarin hydrochloride (Leu-MCA, Sigma) and 4-methylumbelliferyl -D-glucopyranoside (MUF-Glc, Sigma) as a protein and carbohydrate model substrate, according to [START_REF] Hoppe | Use of fluorogenic model substrates for extracellular enzyme activity (EEA) of bacteria[END_REF]. Different substrate concentrations were added to aliquot fractions of sampled water to determine the enzymatic kinetics as follows: 10 final concentrations, ranging between 1.9 and 1000 µM for AMPase and 7 final concentrations ranging between 0.39 and 100 µM for -GLCase.

Incubations were performed in the dark and at the in situ temperature for 1 to 15 and 3 to 15 h for the Leu-MCA and MUF-Glc assays, respectively, depending on the sea temperature and consequently on the level of enzymatic activity. Sodiumdodecyl sulfate (SDS 1% final concentration) and NH 4 -glycin buffer (pH 10.5) were added to stop the reactions of the AMPase and -GLCase activity assays (Chróst 1992[START_REF] Delmas | SDS-preservation for deferred measurement of exoproteolytic kinetics in marine samples[END_REF]. Negative controls were performed for each substrate concentration by addition of an inhibitor (SDS or NH 4 -glycin, respectively) at the beginning of the incubation. Samples were immediately frozen at -20°C and stored for 1 to 2 wk before the fluorescence reading was performed (Kontron SFM 25 spectrofluorometer: excitation, 380 and 364 nm, emission 440 and 450 nm for Leu-MCA and MUF-Glc, respectively). Solutions of 7-amino-4-methylcoumarin (AMC, Sigma) (20 to 2000 nM) and 4-methylumbelliferon (MUF, Sigma) (2 to 2000 nM) were used as calibration standards. The enzymatic activity, expressed as the maximum of potential activity (V max , nmol l -1 h -1 ), was estimated through non-linear regression (Hyper32 software, JS Easterby). The activity of AMPase is considered to be an indicator of bacterial production in MOB based on the significant relationship shown between AMPase (V max ) and bacterial production (BP), which was assessed using tritiated thymidine incorporation during a previous survey (2002 to 2003, same sampling station): log BP = 0.72 logAMPase -1.14, r² = 0.67, n = 10, p = 0.003 (D. Delmas, unpubl. data). Therefore, the ratio of AMPase:-GLCase V max could be calculated and considered to be the expression of the nutritional quality indicator for each substrate [START_REF] Polymenakou | Benthic microbial abundance and activities in an intensively trawled ecosystem (Thermaikos Gulf, Aegean Sea)[END_REF].

Bacterial community composition: Bacterial community composition was determined by the in situ fluorescent hybridization method [START_REF] Delong | Phylogenetic stains: Ribosomal RNA-based probes for the identification of single cells[END_REF]) using 8 Cy-3 labeled probes (MWG-Biotech): Eub 338 I, II and III, specific for Eubacteria), Pla5a for Planctomycetes, CF319a for the Cytophaga-Flavobacterium cluster (CFB), Alf1b for Alphaproteobacteria, Bet43a for Betaproteobacteria, Gam42a for Gammaproteobacteria, and a non-specific control probe (NON) [START_REF] Amann | Phylogenetic identification and in situ detection of individual microbial cells without cultivation[END_REF][START_REF] Bouvier | Key role of selective viral-induced mortality in determining marine bacterial community composition[END_REF].

Archaea was detected using Arch915 [START_REF] Amann | Phylogenetic identification and in situ detection of individual microbial cells without cultivation[END_REF]. Samples fixed in paraformaldehyde (2% final concentration) were kept overnight at 4°C before being stored at -20°C. Thawed subsamples were filtered onto 0.2 m polycarbonate membranes (Whatman), rinsed with 0.2 m filtered Milli-Q water, air-dried and stored at -20°C. The hybridization of probes (2.5 ng l -1 ) to bacteria was carried out for 4.5 h at 46°C before the membrane was washed at 48°C in a pre-warmed (48°C) wash solution for 15 min (102 mM NaCl [except for CFB: 80 mM NaCl]; 20 mM Tris, pH 8; 5 mM EDTA, pH 8; 0.01% SDS). The membrane was then rinsed, first with 0.2 m filtered Milli-Q water, then with 80% ethanol, and then airdried. Bacterial cells were counter-stained with DAPI (0.25 g l -1 ) in a mounting solution of Citifluor (Biovalley) and Vectashield (Ablys) mixed to a ratio of 4:1. A minimum of 250 cells were enumerated using an epifluorescence microscope (Axioskop, Zeiss) equipped with Cy3 (blue) and DAPI (UV) filters.

Viral community structure (RAPD-PCR fingerprinting).

Approximately 12 l water samples were filtered using a peristaltic pump through a succession of nylon fiber sieves (from 250 down to 25 µm) and sequentially though 2 membranes barriers (3 and 0.8 µm pore size cellulose acetate filters, 142 mm diameter; Sartorius) to remove sediment particles and large microorganisms; a final 0.2 m filtration (cellulose acetate membrane, 142 mm; Sartorius) resulted in a bacteria-free sample (Auguet et al. 2008). Based on [START_REF] Suttle | Use of ultrafiltration to isolate viruses from seawater which are pathogens of marine phytoplankton[END_REF] and using the protocol of Auguet et al. (2008), viruses were then concentrated by a tangential ultrafiltration through a 30 kDa polysulfon cartridge (Ultraslice support, Sartorius) into 200 to 300 ml. After storage at -20°C for 1 yr for the oldest samples, viral concentration was achieved by ultracentrifugation at 150000g at r max for 3. RNase-free DNase buffer or 1× of HotStart Taq polymerase reaction buffer, 1 µl Lambda DNA (0.46 µg) or T2 virus suspension or the monthly samples (viral titer: 1 × 10 7 viruses ml - 1 ) and nuclease-free water for a final volume 20 µl. The DNase reaction was performed over 30 min at 37°C and stopped by heating at 75°C for 15 min. One microliter of each DNase reaction was used as template DNA for RAPD-PCR amplification.

Prior to RAPD-PCR, virions need to be decapsidated by heating at 60°C for 10 min and rapidly cooling to 4°C until amplification [START_REF] Fuhrman | Prokaryotic and viral diversity patterns in marine plankton[END_REF]. The efficiency of decapsidation was also tested with a DNase treatment performed on both encapsidated and decapsidated T2 virions. Decapsidation efficiency was confirmed with the hydrolysis of DNA of decapsidated T2 as opposed to encapsidated T2 (positive control). RAPD-PCR was then performed with, as template DNA, results from DNase and decapsidation tests, monthly samples DNase treated and the 12 monthly samples.

Adapted from [START_REF] Winget | Randomly amplified polymorphic DNA PCR as a tool for assessment of marine viral richness[END_REF], amplification was performed with 2 primers (OPA13 [5'-CAG CAC CCA-3'] and CRA22 [5'-CCG CAG CCA-3'], MWG-Biotech Inc.) and using 10 M dNTP, 4 M primer (acting forwardly and reversely), 1× Taq polymerase reaction buffer, 2 U Taq polymerase (HotStart, Qiagen) and 1 l viral concentrate compared to a constant viral titer (2 × 10 6 viruses) with CRA22. The PCR was carried out as follows: (1) 94°C for 15 min, (2) 35°C for 3 min, (3) 72°C for 1 min, (4) 94°C for 30 s, (5) repeat Steps 2 to 4 for 30 cycles, (6) 72°C for 10 min, and (7) hold at 10°C (PTC-100, MJ Research). The PCR products were resolved by 2% MetaPhor agarose gel (Lonza) electrophoresis in 0.5× TBE buffer (Tris-Borate-EDTA) at 4 V cm -1 , stained with Gel Star (Lonza) and visualized with Gel Doc 2000 (Bio-Rad).

Banding patterns were analyzed with the Diversity Database program (Version 2.2.0, Bio-Rad) which provided banding pattern alignment analyses, qualitative analyses of the number of bands for each sample (presence versus absence) and quantitative analyses using densitometry. The total number of observed bands, the RAPD-PCR amplicons, was defined as the fingerprint richness. For each month, triplicates of RAPD-PCR amplifications and gel migrations were used to test the reproducibility of the RAPD-PCR amplifications.

Statistical analyses. Regression and correlation analyses (Spearman rank correlation,

) were performed with Xlstat (the correlation significance was tested using a 2-tailed t-test).

RAPD-PCR banding patterns were first compared with the 'unweighted pair group method with arithmetic mean' in the Diversity Database program to assess the reproducibility of RAPD-PCR banding patterns and the effect of DNase treatment. Second, viral diversity data were standardized for 'species' in order to reduce the dominant contribution of abundant species in a Bray-Curtis similarity matrix and then normalized by fourth-root transformation.

Environmental data were log-transformed under draftsman plot information before being represented in a similarity matrix (Euclidean distance). Non-metric multi-dimensional scaling (MDS) ordination was performed from both similarity matrices (viral diversity and environmental data). Non-parametric analysis of similarity (ANOSIM) tested the robustness of MDS clustering (good if R > 0.6). Third, attempts to link viral assemblages to environmental data were performed to find the best explanatory variables of the observed RAPD-PCR banding pattern. The link was carried out using the BEST procedure in the Primer 6 software [START_REF] Clarke | A method of linking multivariate community structure to environmental variables[END_REF][START_REF] Clarke | Evidence that viral abundance across oceans and lakes is driven by different biological factors[END_REF]. This procedure combined the BIO-ENV and Stepwise procedures and consisted in performing a Spearman rank correlation analysis based on the computation between the elements of 2 triangular matrices: a resemblance matrix of the viral assemblage (Bray-Curtis distance) and an explanatory data matrix of environmental data (Euclidean distance). The BIO-ENV procedure selects a subset of environmental data to maximize the matching Spearman rank correlation ( coefficient) with the 2 assemblage similarities. It identifies numerically an environmental data subset sufficient to explain all the biotic structure of the viral assemblage. A global BEST match permutation test assesses the significance of rho. A linkage tree was then applied with the LINKTREE algorithm, corresponding to binary dividing clustering of the monthly viral assemblage structures. Each cluster was chosen to maximize the ANOSIM R statistic and was interpreted in terms of inequalities in explicative environmental variables. An absolute measure of cluster differences is given by the B% value, where low values correspond to closer samples on the MDS plot. In fine, this LINKTREE procedure was validated by the similarity profile test (SIMPROF test), a permutation test used to demonstrate the significant evidence of an internal group structure. In our study case, this complete procedure selected the optimum combination in an 'environmental' variable subset, maximizing rank-correlation () between the viral assemblage structure and environmental matrices. The 'environmental matrix' was composed of 19 biotic and abiotic variables, and we chose to test the best combination of 1 to 9 out of the 19 available.

RESULTS

Physical and chemical parameters

Water temperature and salinity data are presented in Table 1 together with measurements of inorganic nutrients (for details see Fig. S1 in the supplement at www.interres.com/articles/suppl/ a064p223_ supp.pdf): NO 3 remained above 20 mol l -1 most of the year, except in the autumn, and peaked at 115 mol l -1 in February (riverine sampling) and at 67 mol l -1 in March. PO 4 fell below 0.5 mol l -1 from February until the end of June, while it peaked towards the end of August and reached its highest value (2.0 mol l -1 ) in November.

Seston was composed on average of 20% organic and 80% inorganic matter. With an annual mean of 18.7 ± 13.4 mg l -1 , it peaked twice, in February (53.8 mg l -1 ) and July (27.5 mg l -1 ).

The annual mean chl a concentration was 2.2 ± 1.2 µg l -1 . Three major peaks occurred, in February, April and July (2.1, 3.5 and 4.5 µg l -1 , respectively).

Microorganism abundance

Large phytoplankton (>10 m, 2.1 [±1.9] × 10 4 cells l -1 annual mean) and small phytoplankton (<10 µm, 1.1 [±0.9] × 10 7 cells l -1 annual mean) showed a simultaneous increase only in October. Three other major peaks occurred earlier in the year, in April (5.2 × 10 4 cells l -1 , mainly diatoms), June (4.8 × 10 4 cells l -1 , diatoms and Dinophyceae) and August (4.0 × 10 4 cells l -1 , almost exclusively Dinophyceae) (Fig. 1). Small phytoplankton peaked in March (mainly Synechococcus), July (mainly picophytoplankton) and a major peak occurred in October (due to both).

The maximum variation in virus abundance was the 2.9-fold difference observed between 1.3 × 10 7 viruses ml -1 in April and 3.75 × 10 7 viruses ml -1 in February. There were 3 periods when an increase in abundance occurred: in May, August and November. With an annual mean of 2.67 (±0.95) × 10 6 cells ml -1 , bacterial abundance varied from 1.53 × 10 6 cells ml -1 in December to 4.54 × 10 6 cells ml -1 in August. The dynamics of abundance featured 3 peaks, in February, April and August. The February and August peaks coincided with high values of viral abundance. In April it was decoupled from the viral dynamics and instead coincided with increases in diatom abundance (Fig. 1). The relationship between bacteria and viruses was weak ( = 0.15, p = 0.63 and r² = 0.051, p = 0.48 in linear regression where viruses were considered as a dependant variable), resulting notably from the clear uncoupling periods observed in April-May-June and October-November (grey area in Fig. 1, Table 2).

The annual mean heterotrophic nanoflagellate (HNF) concentration was 3.31 (±4.95) × 10 6 cells l -1 . Composed of 82% small size organisms (<3 m), HNF had 2 periods of growth, one in April (2.54 × 10 6 cells l -1 ) and a major one at the end of October, reaching 1.77 × 10 7 cells l -1 (Fig. 1). This high autumnal peak, 5 times higher than the annual mean, was characterized by the exclusive development of 2 m HNF. Pigmented nanoflagellates represented <10% of the total nanoflagellates, except in April and May when they reached 18 and 17%, respectively.

The dynamics of ciliates (3.78 [±3.48] × 10 3 cells l -1 annual mean) was characterized by 3 growth periods, in February (4.45 × 10 3 cells l -1 ), July (1.23 × 10 4 cells l -1 ) and November (4.4 × 10 3 cells l -1 ) (Fig. 1). No synchronization with HNF dynamics occurred, except during the small April peak. Ciliates were significantly related to total bacterial abundance ( = 0.82, p = 0.001; Table 2). Mesozooplankton (0.43 [±1.2] × 10 4 ind. m -3 annual mean) peaked mainly in April (3.9 × 10 4 ind. m -3 ) and at the end of October (2.1 × 10 3 ind. m -3 ) (Fig. 1).

Bacterial enzymatic activity

Total and specific exoenzymatic activities evolved similarly ( = 0.89 and 0.95 for AMPase and -GLCase, respectively, p < 0.001) so only the total activity results were taken into account. Both potential AMPase and -GLCase activities were significantly related ( = 0.79, p = 0.002) and had a concomitant major peak in April (V max = 630 and 10 nmol l -1 h -1 , respectively) (Fig. 2). During the rest of the year, V max AMPase peaked once in July (394 nmol l -1 h -1 ), while a -GLCase peak occurred in August. With an annual mean of 189 ± 187 nmol l -1 h -1 , AMPase activity was significantly higher than -GLCase activity (4.1 ± 2.9 nmol l -1 h -1 ). The AMPase:-GLCase ratio varied between 13 in January and 78 in July; it increased during the spring period and at the end of October. AMPase and -GLCase activities were closely and positively related to the chl a level ( = 0.90, p <0.001 and  = 0.82, p = 0.001, respectively) and also to an abundance of large phytoplanktonic cells (diatoms:  = 0.73, p = 0.007 for AMPase and  = 0.73, p = 0.007 for -GLCase).

Nevertheless, the major peak of Dinophyceae in August coincided with the second peak of -GLCase activity. The AMPase:-GLCase ratio was positively correlated with the level of chl a ( = 0.71, p = 0.015) (Table 2).

Bacterial community diversity

Based on fluorescence in situ hybridization (FISH) results, Bacteria (Eub probe) represented between 40% in May and 81% in April of the total prokaryote community (DAPI counts) (Fig. S2 in the supplement at www.int-res.com/articels/suppl/a0664p233_supp.pdf). The higher percentages were observed during winter and the first part of spring. Cells of the Archaea domain were scarce (4.6 [±6.9] × 10 4 cells ml -1 annual mean), representing <1% of the total prokaryote abundance for most of the year, except in January and at the end of October (9.7 and 9.2% of the prokaryotic community, respectively). Among all prokaryotes, from 42 to 81.5% of DAPI-stained cells were consistently detected. The frequency of FISHtargeted cells, belonging to the domains Eubacteria and Archaea, was in agreement with previous observations from many marine and estuarine ecosystems (reviewed in [START_REF] Bouvier | Factors influencing the detection of bacterial cells using fluorescence in situ hybridization (FISH): a quantitative review of published reports[END_REF]. Within the Eubacteria domain, Betaproteobacteria and Gammaproteobacteria were the dominant bacterial groups, representing up to 12 and 13% of total bacteria, respectively, while Alphaproteobacteria cells generally remained below 1%, except in July with 1.41% (Fig. S2). Gammaproteobacteria had an annual mean of 2.2 (±1.9) × 10 5 cells ml -1 and dominated from June to September, reaching up to 5.7 × 10 5 cells ml -1 , and in November (2.5 × 10 5 cells ml -1 ). Betaproteobacteria (7.9 [±7.6] × 10 4 cells ml -1 annual mean) were dominant in winter (2.4 × 10 5 cells ml -1 maximum reached in March) and at the beginning of October (up to 9.1 × 10 4 cells ml -1 ). Betaproteobacteria was the only bacterial group that showed a significant correlation with the Charente River runoff ( = 0.88, p < 0.001). Finally, the CFB group and Planctomycetes always represented <0.5% of the total community.

Viral community composition

The repeated profiles of both amplification and gel triplicates showed an average similarity of 84.3%, ranging from 73.3 to 100% homology (n = 12, data not shown).

The viral origin of the RAPD-PCR products was confirmed by comparing the banding patterns of amplicons from viral concentrates with or without DNase treatment. The efficiency of both DNase treatment and decapsidation was standardized with a lambda virus genome and T2 particles, respectively (Fig. S3A,B in the supplement at www.intres.com/articels/suppl/a0664p233_supp.pdf). Notably, since no difference was observed between the use of both buffers in the DNase-treatment, we chose the Taq-polymerase buffer to avoid any possible inhibition in the PCR. Moreover, the higher the DNase quantities in contact with decapsidated T2 DNA, the weaker the subsequent amplification of free DNA (Fig. S3B). Our DNase test procedure was efficient at hydrolyzing T2 DNA only after decapsidation of the genome, even if some DNase-resistant templates tended to remain, inferring the successful liberation of the viral genome during decapsidation. Concerning environmental viral samples, DNase treatment on viral particles did not affect the banding profiles of the templates (Fig S3C). The use of 1 and 5 U of DNase induced a drop in banding pattern intensity, but without significant band losses compared to the control (80 and 90.9% of similarity for April and November samples, respectively; Dice Index).

RAPD fingerprint profiles showed up to 46 bands with the CRA22 primer, ranging from 130 to 1600 bp, and up to 18 bands with the OPA13 primer, ranging from 200 to 1000 bp (Fig. 3). The number of distinct amplicons produced in each sample was higher with CRA22 than with OPA13 (14 ± 2.8 vs. 5 ± 2.4 respective annual means). The comparison of profiles, generated by both primers over time, revealed a similar trend in the changes of fingerprint richness. However, no significant correlation was observed for the number of observed bands between the primers ( = 0.43, p = 0.13).

The fingerprints of the last 4 mo of the year were clearly different from those of the other 8 mo of the year with both CRA22 and OPA13. The mean size of fingerprint amplicons was significantly higher during these months than during the rest of the year (630 ± 74 and 655 ± 24 bp vs. 503 ± 56 and 498 ± 153 bp with CRA22 and OPA13 primers, respectively; ttest, p < 0.05). The viral community composition obtained with CRA22 was characterized by the occurrence of specific bands in April (170 and 220 bp amplicons), May (150,275,300 and 1400 bp amplicons) and June (160 bp amplicon) (Fig. 3). Moreover, amplicons of around 300 bp were of higher intensity for these 3 mo. Concerning the OPA13 profiles, a 200 bp amplicon characterized the month of April, while a 260 bp band was present throughout the year.

Changes in viral community composition

The MDS ordination performed either on viral community composition defined with CRA22 or on environmental variables described a good ordination with no real prospect of misleading interpretation (plot stress value of 0.1). The MDS ordination of viral community composition (Fig. 4A1) showed 3 different clusters separated at the 35% similarity threshold (2) May, June and July; and (3) October, 22nd October, November and December (Fig. 4A3).

Using the extracted explanatory data, a clustering dendrogram (LINKTREE procedure) of the monthly viral community structure showed 4 significant clusters (Fig. 4B): (1) the spring months (4, 5, 6 and 7) with an 84% difference (B%) compared with the rest of the year (R = 0.53, p = 0.001); (2) the winter months (1, 2, 3) with a 69% difference with the rest of the samples (R = 0.56, p = 0.001); and (3) August ( 8) distinguished (86%) from the last cluster (4) composed of the autumn months (9 to 12) (R = 1, p = 0.001).

The clustering, based on the distribution of viral amplicon band sizes, can therefore be explained by examining the explanatory variables. The winter cluster (1, 2 and 3, CRA22 primer) was characterized by a low level of AMPase activity (<127 nmol l -1 h -1 ) combined with few HNFs (<7.9 × 10 5 cells l -1 ). A lower AMPase level (<127 nmol l -1 h -1 ) occurred in late summer and autumn (8 to 12), combined with a higher HNF abundance (>2.2 × 10 6 cells l -1 ). The August sample was distinct from the autumnal months (9 to 12) because of its high bacterial abundance (>4.5 × 10 6 cells ml -1 ) and its high AMPase concentration (>126 nmol l -1 h -1 ). Spring and early summer months were distinguished differently according to the primer used. With CRA22, 4, 5, 6 and 7 were grouped without distinction, and their viral community composition was related exclusively to high AMPase (>297 nmol l -1 h -1 ). With OPA13, nanophytoplankton abundance appeared to be an explanatory parameter of viral diversity, distinguishing June from the rest of the samples (>3.2 × 10 6 cells l -1 ).

The OPA13 primer LINKTREE analysis was performed, but the results should be treated carefully because of the very low number of bands detected in each sample. However, AMPase and HNF were confirmed as explanatory variables, although, interestingly, salinity and river runoffs interceded in the changes in viral community composition, setting apart February and March (data not shown).

DISCUSSION

Methodological considerations

On a temporal scale, the entire microbial community structure is likely to change within a day or week, with different generation times and temporal variability in biological responses to biological and environmental factors (Hewson & Fuhrman 2006[START_REF] Fuhrman | Microbial community structure and its functional implications[END_REF]). However, considering the number of parameters examined here, a long-term survey using a weekly or even a shorter time scale was inconceivable.

Fixation and conservation of viral samples are known to be quite uncertain processes that may induce variable losses of total viral counts compared to in situ concentrations [START_REF] Brussaard | Optimization of procedures for counting viruses by flow cytometry[END_REF]). Drops of from 20 to 35% in viral abundance have been reported after a few hours in the presence of fixatives (glutaraldehyde or formaldehyde) and storage at 4°C with pelagic samples [START_REF] Wen | Accurate estimation of viral abundance by epifluorescence microscopy[END_REF]), as well as with viruses from estuarine sediments [START_REF] Helton | Assessment of factors influencing direct enumeration of viruses within estuarine sediments[END_REF]. However, in the present study, all samples were treated in the same way, allowing valuable cross-comparisons.

Seasonal variations in bulk viral abundance (2.75 [±0.8] × 10 7 viruses ml -1 ) and fluctuations due to the concentration processes (coefficient of variation of 76%) resulted in variable viral titers in the target samples (from 2 × 10 6 to 4 × 10 7 viruses ml -1 ) and may have introduced a possible bias in the amplification of rare species. In order to evaluate this bias and given that PCR inhibitors prohibited the addition of >1 µl of viral suspension, we also performed RAPD-PCR by adding a constant abundance of viruses (2 × 10 6 viruses) based on the least concentrated sample (Fig. 5). Both amplification reactions using the CRA22 primer provided the same results for the specific richness and relative quantity of amplicons between the 2 patterns (data not shown). Considering that the use of a specific viral titer was no more representative of the bulk diversity than the use of a unique volume and that the use of a constant percentage of bulk abundance was not possible, we regard the choice of identical volumes (1 µl) as a valid way of studying the structure of viral communities. RAPD-PCR fingerprinting was considered as a useful analysis to assess viral community composition, assuming the uniqueness of viral species origin for each amplicon. However, recent in silico investigations confirmed current thinking by reporting that 26 to 33% of targeted phage genome (of a total of 85 complete genomes) may be amplified in 2 or more fragments (Winter & Weinbauer 2010). Nevertheless, this fingerprinting method is sufficiently accurate to investigate changes in viral community structure, while the choice of primer is preponderant in discriminating between different subsets of environmental viruses. We confirmed that CRA22 is more discriminative than OPA13 (in terms of number of bands), although the latter may be useful for the characterization of the potentially different responses of targeted population subsets.

Monthly dynamics

This 2007 survey gave unexpected results concerning micro-organism dynamics and their relationships (Fig. 1, Table 2). First, the annual means of primary and secondary producers were consistent with previous data only for Dinophyceae and Synechococcus, but were otherwise lower [START_REF] Auguet | Dynamic of virioplankton abundance and its environmental control in the Charente estuary (France)[END_REF][START_REF] Ory | Pelagic food web patterns: do they modulate virus and nanoflagellate effects on picoplankton during the phytoplankton spring bloom ?[END_REF]. Indeed, among small phytoplankton cells, pico-and nano-phytoeukaryotes showed an annual mean that was about 3 times lower than that of the previous year. The most important difference concerned diatoms, which had a 6-fold lower mean abundance in 2007. Mean bacterial and viral abundances (1.5-fold lower than in 2003 and 2006) remained quite constant over the year, without the usual large summer peaks [START_REF] Auguet | Dynamic of virioplankton abundance and its environmental control in the Charente estuary (France)[END_REF][START_REF] Ory | Pelagic food web patterns: do they modulate virus and nanoflagellate effects on picoplankton during the phytoplankton spring bloom ?[END_REF]. However, the periods of viral increase were close to those obtained in the 2006 survey, although this was not the case for bacterial dynamics.

The evaluation of bacterial cell activity based on the measurement of exoenzymatic activity was consistent with results obtained previously: in 2003, at the same sampling station, potential AMPase activity ranged between 0.02 and 1.3 µmol l -1 h -1 (D. Delmas unpubl. data). Although bacterial activities in 2007 did not reach the level attained in 2003, cellspecific proteolysis activity was similar during these 2 yr (mean of 66 ± 63 amol cell -1 h -1 vs.

115 ± 121 amol cell -1 h -1 ). In addition, the values of the AMPase and -GLCase activities described here fell within the range of the published data from coastal or estuarine waters (Chróst & Velimirov 1991[START_REF] Karner | Extracellular enzymatic activity and secondary production in free-living and marine-snow-associated bacteria[END_REF][START_REF] Cunha | Inorganic nutrient regulation of bacterioplankton heterotrophic activity in an estuarine system (Ria de Aveiro, Portugal)[END_REF]. Moreover, -GLCase activity was between 10-and 100-fold lower than AMPase activity, while the proteolysis:glycolysis ratio tended to be lower in the inner part of an estuary (8 to 160; [START_REF] Polymenakou | Benthic microbial abundance and activities in an intensively trawled ecosystem (Thermaikos Gulf, Aegean Sea)[END_REF]) compared to coastal areas (about 1000; [START_REF] Christian | Bacterial ectoenzymes in marine waters: Activity ratios and temperature responses in three oceanographic provinces[END_REF].

However, the glycolysis potential was not as strongly linked to microbial communities as that of AMPase activity (Table 2, Fig. 4). Indeed, the low values of -GLCase V max could be below the sensitivity threshold for a study of the co-variations between community structure and activity. Moreover, a negligible importance of carbohydrates in bacterial diet has been shown [START_REF] Zoppini | Extracellular enzyme activity and dynamics of bacterial community in mucilaginous aggregates of the northern Adriatic Sea[END_REF]. Nevertheless, higher uptakes of carbohydrates by bacteria have been described for members of the Alphaproteobacteria, at low ambient glucose concentrations [START_REF] Alonso | Roseobacter and SAR11 dominate microbial glucose uptake in coastal North Sea waters[END_REF]. Indeed our sampling site, located at the river mouth, is characterized by few alpha-Proteobacteria (Fig. S2).

Interestingly, ciliates were the only predators that correlated with bacterial abundance in the present study (Table 2). There is still a question mark over the prey selectivity of ciliates on bacteria, but it may be ciliate species specific and depend on the type and size of bacterial prey [START_REF] Bernard | Bacteria or microflagellates as a major source for marine ciliates: possible implications for the microzooplankton[END_REF][START_REF] Ayo | Grazing rates of diverse morphotypes of bacterivorous ciliates feeding on four allochtonous bacteria[END_REF]. Abundance and dynamics of grazers (ciliates and HNF) generally followed the same pattern of increase and decline as in the 2006 MOB study [START_REF] Ory | Pelagic food web patterns: do they modulate virus and nanoflagellate effects on picoplankton during the phytoplankton spring bloom ?[END_REF], except that a much greater peak of HNF was observed in the autumn of 2007 (Fig. 1). Indeed, particularly strong increases were observed during this period for components of the microbial food web, a trophic pathway comprising small cellproducers and their protozoan predators [START_REF] Legendre | Plankton and nutrient dynamics in marine waters[END_REF]. Primary production in the autumn appeared to be based on small autotrophic cells such as picophytoeukaryotes and Synechococcus. Considering that HNF are major picoplankton grazers [START_REF] Monger | Feeding selection of heterotrophic marine nanoflagellates based on the surface hydrophobicity of their picoplankton prey[END_REF], the high peak of picophytoeukaryotes, as well as the small increase in bacteria, could sustain a strong increase in HNF, ciliates and lastly mesozooplankton (Fig. 1). During this productive period, Synechococcus were less abundant than picophytoeukaryotes, but their growth phase lasted longer. While picophytoeukaryotes may represent the main source of carbon for microzooplankton (ciliates), not only did Synechococcus appear to be weakly impacted by grazing, but also low viral lysis rates have been reported [START_REF] Suttle | Dynamics and distribution of Cyanophages and their effect on marine Synechococcus spp[END_REF][START_REF] Worden | Assessing the dynamics and ecology of marine picophytoplankton: the importance of the eukaryotic component[END_REF][START_REF] Baudoux | Estimates of virus-vs. grazing induced mortality of picophytoplankton in the North Sea during summer[END_REF].

Virus-Bacteria relationship

No significant relationship existed between the dynamics of virus and bacteria in contrast to what is usually observed in most marine (r² = 0.94 to 0.67; [START_REF] Boehme | Viruses, bacterioplankton, and phytoplankton in the southeastern Gulf of mexico: distribution and contribution to oceanic DNA pools[END_REF][START_REF] Cochlan | Spatial distribution of viruses, bacteria and chlorophyll a in neritic, oceanic and estuarine environments[END_REF], freshwater and estuarine areas (r² = 0.85; [START_REF] Auguet | Dynamic of virioplankton abundance and its environmental control in the Charente estuary (France)[END_REF]. However, 3 reviews of studies performed in freshwater sites and marine sites described the large variation in the strength of the relationship (r²) between viral and bacterial abundance and also in terms of the power-slope of the regression [START_REF] Maranger | Viral abundance in aquatic systems: a comparison between marine and fresh waters[END_REF][START_REF] Gasol | Comparative analyses in aquatic microbial ecoogy: how far do they go?[END_REF][START_REF] Clarke | Evidence that viral abundance across oceans and lakes is driven by different biological factors[END_REF]). In the present study, relatively weaker abundance oscillations and uncoupled viral and bacterial dynamics for 5 samples out of 12 explained this poor relationship (Fig. 1). This uncoupling may result from a relatively greater lytic induction of lysogenic bacteria and/or variations in the viral burst size according to temperature or trophic status, leading to a release of more virions [START_REF] Parada | Viral burst size of heterotrophic prokaryotes in aquatic systems[END_REF]. Recently, [START_REF] Motegi | Viral control of bacterial growth efficiency in marine pelagic environments[END_REF] suggested that nutrient-rich conditions entailed increasing the expression of the antiphage traits of bacteria and/or the destruction of viral capsids by bacterial extracellular proteases.

Furthermore, the complexity of the MOB ecosystem is under the temporal variations of inputs from the Charente River (allochthonous viruses and bacteria; [START_REF] Auguet | Dynamic of virioplankton abundance and its environmental control in the Charente estuary (France)[END_REF].

Changes in morphological (transmission electron microscopy) and virus-like genome structures (pulsed-field gel electrophoresis, PFGE) were reported along the salinity gradient of the Charente River [START_REF] Auguet | Structure of virioplankton in the Charente estuary (France): Transmission electron microscopy versus pulsed field gel electrophoresis[END_REF]). In addition, riverine viruses transferred to MOB could express their infectious potential through estuarine bacterial communities (Auguet et al. 2008). Here, the use of the 2 primers, CRA22 and OPA13, may describe 2 temporal changes of viral community structures since they target 2 different subsets of the viral community (Winter & Weinbauer 2010). The link between the viral community subset amplified with OPA13 and the salinity could infer the presence of some virus/host systems of riverine origin (data not shown). Bacterial community composition was also under the influence of Charente River inputs. During winter, when nutrients were not limiting, the dominant group tended to be Betaproteobacteria (Fig. S2). The Betaproteobacteria group is known to belong mainly to the freshwater prokaryote community, but its presence has been reported in coastal waters [START_REF] Giovannoni | Evolution, diversity, and molecular ecology of marine prokaryotes[END_REF]. Currently, it is unclear if this community represents an active marine bacterial group or if it is passively carried with freshwater outputs to the coastal zone.

The strong correlation between the presence of Betaproteobacteria and the Charente River outputs corroborates the second hypothesis, but their constant presence and their reactivity to the manipulation of flagellates or viruses [START_REF] Ory | Pelagic food web patterns: do they modulate virus and nanoflagellate effects on picoplankton during the phytoplankton spring bloom ?[END_REF] support their involvement in the coastal microbial loop.

Temporal dynamics of the viral community structure

The potential link between viral and bacterial abundance or bacterial community composition (BCC) was not clearly demonstrated here. The viral phylogenetic level analyzed may have been too detailed to detect significant correlations with the taxonomic groups of prokaryotes (Bacteria and Archaea). Interestingly, a strong statistical link was found between the variation in viral community structure and the potential of bacterial exoenzymatic activity, particularly AMPase, despite the difference in analysis level (species level for viral community structure and overall community level for bacterial exoenzymatic activities; Fig. 4). Thus, bacterial exoenzymatic activity may be considered as an integrator of different processes acting at different trophic levels. The link between the viral fingerprints and AMPase could reflect direct and indirect interactions between viruses and bacteria, taking into account the influences of phytoplankton and grazer activities.

We therefore hypothesized that viral community structure could partly depend on phytoplankton, either directly through phytoplankton viruses, or indirectly through the uptake or transformation by bacteria of DOM or organic compounds excreted from autotrophic cells.

However, this assumption had to be moderated since DOM may also originate from viral lysis, grazers and from allochthonous inputs. Thus, the weaker relationship between viral and bacterial abundance could be the result of viruses infecting other organisms. Indeed, links between viral abundance and autotrophic parameters have already been described [START_REF] Boehme | Viruses, bacterioplankton, and phytoplankton in the southeastern Gulf of mexico: distribution and contribution to oceanic DNA pools[END_REF][START_REF] Maranger | Viral and bacterial dynamics in Arctic sea ice during the spring algal bloom near Resolute, N.W.T., Canada[END_REF][START_REF] Maranger | Viral abundance in aquatic systems: a comparison between marine and fresh waters[END_REF], Guixa-Boixereu et al. 1999[START_REF] Clarke | Evidence that viral abundance across oceans and lakes is driven by different biological factors[END_REF][START_REF] Jackson | Viruses in wetland ecosystems[END_REF][START_REF] Ory | Pelagic food web patterns: do they modulate virus and nanoflagellate effects on picoplankton during the phytoplankton spring bloom ?[END_REF]. Moreover, phytoplankton community composition has already been proposed as a proxy for evaluating the organic carbon resources available to aquatic bacteria and driving shifts in bacterial community dynamics and structure, although bacteria may be synchronizing phytoplankton dynamics through nutrients [START_REF] Kent | Synchrony in aquatic microbial community dynamics[END_REF]. Proteolysis activity may, therefore, reflect phytoplankton-bacteria interactions and their relationships with viruses, supported by the strong statistical correlation between chl a and AMPase activity (Table 2).

Among the factors explaining the variation in viral community structure (Fig. 4), the potential AMPase activity differentiated spring and autumnal phytoplankton blooms, concomitant with the uncoupling in virus-bacteria abundance. Both periods were characterized by high viral richness (as shown by the number of bands) and the emergence of specific bands, but differed in proteolysis level (high in spring and low in autumn). Would changes in bacterial community structure be synchronized or not with changes in viral community structure? Our data were not sufficient to clarify this point.

Virus-bacteria interactions

Assuming that viruses in coastal environments are mainly bacteriophages [START_REF] Wommack | Distribution of viruses in the Chesapeake Bay[END_REF][START_REF] Maranger | Viral abundance in aquatic systems: a comparison between marine and fresh waters[END_REF], their multiplication depends on the density of sensitive and permissive bacterial cells, whose abundance varies according to different factors such as bacterial community structure or physiological status. Thus, the virioplankton lytic activity is known to impact bacterial community composition [START_REF] Hewson | Covariation of viral parameters with bacterial assemblage richness and diversity in the water column and sediments[END_REF]. Many hypotheses suggest that the species diversity of the hosts (specific richness) is not only sustained by adaptation to ecological niches [START_REF] Fuhrman | Annually reoccurring bacterial communities are predictable from ocean conditions The[END_REF]) but also by viral infection.

(1) According to the 'kill the winner' theory, since viral infection depends on bacterial density, viruses may modulate bacterial community composition by regulating the dominant bacterial groups and thus may favor the presence of rare and less competitive populations [START_REF] Thingstad | Elements of a theory for the mechanisms controlling abundance, diversity, and biogeochemical role of lytic bacterial viruses in aquatic systems[END_REF]. (2) Moreover, [START_REF] Bouvier | Key role of selective viral-induced mortality in determining marine bacterial community composition[END_REF], on the basis of in vitro experiments, concluded that virus-sensitive groups may be the rare bacterial groups that are more metabolically competitive for resources. Their potential bloom may be regulated by viral lysis. Nevertheless, the specificity of viral infection and the vulnerability of prokaryotic hosts to >1 viral species remains in question (Winter et al. 2010). Intrinsic changes of bacterial cell-specific activity could induce an increase in bacterial sensitivity and permissivity to viruses by inducing changes in receptor regions of the bacterial cell membranes involved in nutrient uptake and also related to viral adsorption [START_REF] Middelboe | Bacterial growth rate and marine virus-host dynamics[END_REF], in fine resulting in an increase in viral yield [START_REF] Bouvier | Key role of selective viral-induced mortality in determining marine bacterial community composition[END_REF], Auguet et al. 2008). This assumption is relevant to [START_REF] Sandaa | Viral control of bacterial biodiversity -evidence from a nutrient-enriched marine mesocosm experiment[END_REF], who hypothesized that the number of different virus-host pairs is controlled by viral lysis, and their identity, by substrate availability. Thus, viral community structure could be directly related to the most productive bacterial clones and, consequently, to bacterial proteolysis (Fig. 4B).

Virus-phytoplankton interactions

During the spring phytoplanktonic period, when the strongest virus-bacteria uncoupling occurs, phytoplankton viruses could partly regulate the large-phytoplankton bloom, which involves mainly diatoms. From our current knowledge, genomes of eukaryotic microalgal viruses range in size from 200 to 500 kb, those of cyanophages from 30 to 380 kb and those of bacteriophages from 39 to 60 kb [START_REF] Larsen | Population dynamics and diversity of phytoplankton, bacteria and viruses in a seawater enclosure[END_REF][START_REF] Baudoux | Estimates of virus-vs. grazing induced mortality of picophytoplankton in the North Sea during summer[END_REF][START_REF] Sandaa | Burden or benefit? Virus-host interactions in the marine environment[END_REF]). For our study site using a PFGE analysis, a clear relationship was demonstrated in data from a 7 mo survey in 2002/2003: (1) between large virus-like genomes (VLG > 200 kb) and large capsid size (>65 nm) and (2) between small VLGs (<100 kb) and small capsid size (<65 nm) [START_REF] Auguet | Structure of virioplankton in the Charente estuary (France): Transmission electron microscopy versus pulsed field gel electrophoresis[END_REF]). Thus, the occurrence of large VLGs (261 to 300 kb) and the 50% contribution of large tailed phages (105 nm) to total tailed phages, reported in April 2003 [START_REF] Auguet | Structure of virioplankton in the Charente estuary (France): Transmission electron microscopy versus pulsed field gel electrophoresis[END_REF]) during an uncoupling period with a low number of viruses, as in April 2007, may come from phycoviruses or cyanophages.

In autumn 2007, high peaks of small phytoplankton, mainly composed of picophytoeukaryotes and Synechococcus, may have allowed the development of viruses specific for these autotrophic groups (Fig. 1). During this period, the uniqueness of the viral fingerprint pattern also suggested the occurrence of phytoplankton viruses (Figs. 3 & 4).

Previously, in October 2003, 9.8% of the VLP were 105 nm in size with 1% tailed phages (H. Montanié per. obs.), compensating for the presence of phytoplankton or cyanobacterial viruses during this autumnal period when virus numbers were low. As in their response to grazing, small phytoplankton could react differently to viral lysis depending on the cell group:

Synechococcus would be less affected by viral lysis than picophytoeukaryotes [START_REF] Baudoux | Estimates of virus-vs. grazing induced mortality of picophytoplankton in the North Sea during summer[END_REF]. The 2 mo peak of Synechococcus lasted longer than the picophytoeukaryote peak and may corroborate the difference in the response to viruses (sensitivity, latency time and burst size) and in fine in the decline rate of microalgae. Indeed, it is established that the Synechococcus population resisted well to co-occurring phages due to variations in clonal composition according to acquisition of resistance by successive species [START_REF] Waterbury | Resistance to co-occurring phages enables marine Synechococcus communities to coexist with cyanophages abundant in seawater[END_REF][START_REF] Marston | Genetic diversity and temporal variation in the cyanophage community infecting marine Synechococcus species in Rhode Island's coastal waters[END_REF]. In fine, their abundance regulation may therefore rely mainly on microzooplankton [START_REF] Christaki | Growth and grazing of Prochlorococcus and Synechococcus by two marine ciliates[END_REF].

Top-down control by grazers

BCC and bacterial activity varied experimentally with grazing pressure (Šimek et al. 1999(Šimek et al. , Jardillier et al. 2005)). [START_REF] Zöllner | Effect of zooplankton-mediated trophic cascades on marine food web components (bacteria, nanoflagellates, ciliates)[END_REF], discussing the importance of specificity of predation, suggested that mesozooplankton, indirectly via the trophic cascade, sloppy-feeding and phytoplankton top-down control [START_REF] Kragh | Production and bioavailability of autochthonous dissolved organic carbon: effects of mesozooplankton[END_REF], and HNF, directly via predation, could affect the parameters of bacterial community activities more strongly than bacterial abundance (Van Wambeke 1994). The size-specific selectivity of HNF modified bacterial community composition and structure (Šimek et al. 2001). Similarly, in the present study, the apparently size-specific predation of ciliates on bacteria quantitatively related ciliates to the bacterial community, while the presumed host-specificity of viruses may induce qualitative links within the community (Table 2, Fig. 4). Phytoplanktonic blooms were also regulated by top-down control. In spring, this control was driven mainly by mesozooplankton, as happens during a typical phase of the herbivorous food web trophic model [START_REF] Legendre | Plankton and nutrient dynamics in marine waters[END_REF]. In autumn, the small-phytoplankters may be controlled sequentially by HNF, ciliates and mesozooplankton, as in a trophic model of the microbial food web. In 2007, the unusually high peak of nanoflagellates could have affected the number of VHS [START_REF] Sandaa | Viral control of bacterial biodiversity -evidence from a nutrient-enriched marine mesocosm experiment[END_REF]), as suggested by the higher viral richness observed in autumn.

Resource-mediated control

Bacterial activities were also under bottom-up pressure via the quantity and quality of available labile matter, mainly released by phytoplankton (Winter et al. 2010). Indeed, bacterial exoenzymatic activity depends on: (1) the quantity and quality of the labile substrate, mainly excreted by phytoplankters, by sloppy-feeding of grazers or by the bursting of virusinfected phytoplankton or bacterial cells [START_REF] Cherrier | Utilization and turnover of labile dissolved organic matter by bacterial heterotrophs in eastern North Pacific surface waters[END_REF][START_REF] Middelboe | Regeneration of dissolved organic matter by viral lysis in marine microbial communities[END_REF], Winter et al. 2004a[START_REF] Brussaard | A mesocosm study of phaeocystis globosa (Prymnesiophyceae) population dynamics -II. Significance for the microbial community[END_REF]) and (2) the BCC-comprising species, which may have different proteolysis or glycolysis capacities [START_REF] Riemann | Dynamics of bacterial community composition and activity during mesocosm diatom bloom[END_REF]. Thus, the supplies of organic matter may affect the magnitude of coupling between viruses and bacteria in the marine environment [START_REF] Motegi | Viral control of bacterial growth efficiency in marine pelagic environments[END_REF].

The phytoplankton bloom enhanced the production of fresh organic matter (largely composed of polysaccharides and proteins; Ittekkot 1982) both directly through excretion or indirectly through grazing and sloppy-feeding [START_REF] Jumars | Closing the microbial loop: dissolved carbon pathway to heterotrophic bacteria from incomplete ingestion, digestion and absorption in animals[END_REF]), which in turn induced higher potential bacterial enzymatic activity [START_REF] Biddanda | Carbon, nitrogen, and carbohydrate fluxes during the production of particulate and dissolved organic matter by marine phytoplankton[END_REF][START_REF] Misic | Ectoenzymatic activity and its relationship to chlorophyll-a and bacteria in the Gulf of Genoa (Ligurian Sea, NW Mediterranean)[END_REF]. In MOB, dissolved organic carbon (DOC) was mainly a large fraction of labile matter (around 20% assuming that both dissolved carbohydrates and proteins were labile; D. Delmas unpubl. data). The conservative aspect of DOC was proved all along the salinity gradient in the Charente River (i.e. DOC dilution is proportional to water mixing) and inferred significant riverine inputs into the basin [START_REF] Auguet | Dynamic of virioplankton abundance and its environmental control in the Charente estuary (France)[END_REF]. Moreover, the DOC quality was characterized in MOB by lower carbohydrate levels (monomeric as well as polymeric) and a higher dissolved free amino acid concentration, as compared to the Charente River [START_REF] Auguet | Dynamic of virioplankton abundance and its environmental control in the Charente estuary (France)[END_REF]). Thus, in a complex ecosystem like MOB, DOM composition may change according to: (1) high freshwater inputs that increase the carbohydrate concentration in the outer estuary (our sampling site) and (2) tidal resuspension of benthic sediment [START_REF] Wainright | Sediment-to-water fluxes of particulate material and microbes by resuspension and their contribution to planktonic food web[END_REF][START_REF] Hopkinson | Influence of the benthos on growth of planktonic estuarine bacteria[END_REF]. Annual variations in bacterial exoenzymatic activities in MOB reflected the DOC quality and the shifts in consumed DOM, as previously shown by [START_REF] Alonso-Saez | Factors controlling the year-round variablity in the carbon flux through bacteria in a coastal marine system[END_REF]. The proteolysis:glycolysis ratio, which was weak in August and November, corroborated previous time series results in Spain [START_REF] Alonso-Saez | Factors controlling the year-round variablity in the carbon flux through bacteria in a coastal marine system[END_REF]) and the preferential degradation of polysaccharides during these 2 mo under a temperate climate.

The pool of organic matter available to bacteria was continuously fed from viral lysis and phytoplankton, and its quality varied according to the phytoplankton composition, which is also under the influence of viral lysis and grazing. In the present study, the difference in phytoplankton dynamics and composition could thus be responsible for the temporal changes in bacterial AMPase activity, as suggested by the concordance of changes in diatom abundance and proteolysis level (Fig. 2). Bacteria can respond rapidly to dissolved or particulate material associated with the phytoplankton bloom [START_REF] Lancelot | Activity of heterotrophic bacteria and its coupling to primary production during the spring phytoplankton bloom in the southern bight of the North Sea[END_REF]. Indeed, chl a concentrations and potential exoenzymatic activity have already been clearly related [START_REF] Sala | Functional diversity of bacterioplankton assemblages in western Antarctic seawaters during late spring[END_REF]. Synchronous changes in diatom abundance and in proteolysis level in MOB may indicate, as suggested by [START_REF] Kent | Synchrony in aquatic microbial community dynamics[END_REF], that the resource-mediated (bottom-up) control is the mechanism by which an environmental signal is passed along from the phytoplankton to the bacterial communities.

In spring, the diatom bloom could induce the release of HMW compounds, including polypeptides, which stimulate bacterial cell activity via the production of the AMPase enzyme. In contrast, the small phytoplankton autumnal bloom: (1) did not seem to provide sufficient matter to strongly stimulate the potential exoenzymatic activity of bacteria, and their abundance remained low; (2) may repress this activity because of the low quality of released matter (more refractory HMW-DOM); or (3) the excreted LMW-DOM, derived from viral or grazing activity and/or of allochthonous origin, would be preferentially assimilated by bacteria [START_REF] Azua | Influence of organic matter quality in the cleavage of polymers by marine bacterial communities[END_REF][START_REF] Berggren | Efficient aquatic bacterial metabolism of dissolved low-molecular-weight compounds from terrestrial sources[END_REF]. Our results argue in favor of an indirect role for phytoplankton in the resource-mediated control of the type of VHS that is present.

CONCLUSIONS

The in situ survey of 2007 in MOB featured unexpected dynamics of pelagic microorganisms, notably an absence of large viral and bacterial peaks in the summer. Two periods of high viral specific richness with a concomitant uncoupling between viral and bacterial dynamics were clearly distinguished: (1) during the herbivorous spring period, largephytoplankton blooms may be mainly regulated by mesozooplankton yet could sustain bacterial activity; (2) the trophic structure of the autumn months was characterized by the presence of a strong microbial food web promoting the transfer of organic matter from picophytoplankton to HNF and ciliates, but only weakly enhancing the bacterial standing stock. In both cases, the virus-bacteria uncoupling may be explained first by the simultaneous presence of microalgal viruses and the stimulation of bacterial growth by DOM released by phytoplankton and then by standard prey-predator interactions (virus and HNF as predators).

Interestingly, spring and autumnal shifts in viral community structure were related to the level of AMPase activity. We thus provide conclusive evidence on the importance of AMPase activity as a direct indicator of the phytoplankton influence over bacterioplankton, which indirectly attests to the ecosystem trophic pathway (herbivorous vs. microbial food web, concordant with the prevalence of large vs. small autotrophic cells). Overall, we demonstrated that proteolysis activity may reflect the phytoplankton-bacteria interactions and their relationship with viruses. For these reasons, we suggest using the AMPase activity level as a proxy in structuring virioplankton. However, it would be useful in a future work to approach bacterial community structure at the species level using fingerprinting techniques and to specify the origin of DOM using carbon stable isotope analysis. A3) subset of environmental variables 'best matching' the viral pattern, as selected by the Bio-Env procedure. (B) A linkage tree of the viral community composition with a similarity scale bar defined as decreasing B%. Each division was described by inequalities in the explanatory data best matching the viral structure selected by the the Bio-Env procedure. The dotted line describes non-significant sub-groups: the separation of 5 from 4, 6 and 7 was based on bacterial abundance, and 9 from 10, 11 and 12 was based on nanophytoeukaryotic abundance. AMPase: aminopeptidase activity; HNF: heterotrophic nanoflagellate abundance; Nano: nanophytoeukaryote abundance; Bact: bacterial abundance.

Legends

  . By analyzing a new dataset collected during a 1 yr survey in the mesotrophic coastal ecosystem of the Marennes Oléron Bay (MOB), the present study addresses viral assemblage structuring. The objectives were: (1) to assess the structure of the virioplankton assemblage by RAPD-PCR (randomly amplified polymorphism DNApolymerase chain reaction) analysis and (2) to relate changes in viral assemblage with changes in bacterial community composition, together with the dynamics of other microbial entities (grazers and phytoplankton) and environmental parameters, notably the potential consumption of high molecular weight-DOM, by measuring bacterial exoenzymatic activity.

  area. MOB is located in the north of the Bay of Biscay (46°00'N, 1°10'W), along the French Atlantic coast. The bay (surface area = 175 km²; mean depth = 8.6 m) is composed of 60% intertidal mudflats and exposed to a macrotidal system (amplitude = 2 to 6 m). MOB is open to the Atlantic Ocean through a broad northern pass (Pertuis d'Antioche) and a narrowed southern pass (Pertuis de Maumusson). The importance of tidal exchanges through the northern pass results in a northward to southward circulation of oceanic waters from the Bay of Biscay (water residence time = 11 d). The bay is polyhaline to euryhaline (mean salinity = 30.9 ± 4.2‰; Struski 2005), and freshwater arrives mainly from the Charente River (90%), with annual mean runoffs of 100 m 3 s -1 (data: DDE-Charente Maritime/HYDRO-MEDD/DE; http://hydro.eaufrance.fr).

  5 h at 4°C (LE 70 Beckmann ultracentrifuge, SW 28.1 rotor). The pellet of viral concentrate was resuspended in 100 l of SM buffer (0.1 M NaCl, 8 mM MgSO 4 H 2 O, 50 mM Tris-HCl, 0.005% [w/v] glycerol) and stored at 4°C. To examine the potential presence of non-viral DNA in viral concentrate samples, DNase treatment tests were performed with the April 2007 and November 2007 samples and analyzed with RAPD-PCR amplification. Preliminary, optimization of the DNAse test was performed with a Lambda genome (as a positive genome control) and a T2 virus suspension (as a purified virion target), with the aim of optimizing buffer choice (hot-start PCR buffer [Qiagen] versus RQ1 RNase-free DNase buffer) and the number of DNase units used (1 or 5 U; Promega). DNase reaction mixtures contained 1 or 5 U of RQ1 RNase-free DNase, 1× RQ1

(

  based on clustering analysis from the Bray-Curtis similarity matrix): (1) January, February, March, April and August; (2) May, June and July; and (3) 4 October, 22 October, November and December. Non-parametric ANOSIM confirmed the robustness of the 3 clusters (R = 0.863, p = 0.001). The BEST procedure revealed 4 environmental variables that explained the changes in viral community composition: AMPase activity, bacterial abundance, HNF abundance and nanoeukaryote abundance ( = 0.425, p = 0.07, number of permutations = 99). MDS ordination based on the 4 selected experimental variables clustered 5 different groups, including 2 isolated monthly samples (April and August): (1) January, February and March;

Fig. 1 .Fig. 2 .Fig. 3 .Fig. 4 .

 1234 Fig. 1. Temporal dynamics of abiotic and biotic parameters. Survey performed during 2007 in Marennes Oléron Bay: (A) salinity; (B) seawater (SW) temperature; and abundance of (C) Synechococcus, (D) viruses, (E) picophytoeukaryotes, (F) bacteria, (G) nanophytoeukaryotes, (H) heterotrophic nanoflagellates (HNF), (I) Dinophyceae, (J) ciliates, (K) diatoms and (L) mesozooplankton. Only HNF, viral and bacterial abundance have a standard deviation scale. Organisms are classified as monitored with autotrophs on the left part and heterotrophs on the right part of the figure and according to their increasing size. Grey background in Panels D & F marks the 2 periods of uncoupling. VP: viral particles

Fig. 5 .

 5 Fig. 5.Validation of the RAPD-PCR method. Electrophoregrams resulting from amplifications performed with the CRA22 primer of either a constant abundance of viruses (A: 2 × 10 6 viruses) or a constant volume of viral suspension (B: 1 µl). M: 100 bp ladder (Fermentas); numbers from 1 to 12: monthly samples; -: negative control

Acknowledgements. The present study was partly supported by the French PNEC-Littoral Atlantique program. This paper was also a contribution to the AQUAPHAGE project funded by the French ANR Biodiversity program (ANR 07 BDIV 015-05 awarded to H.M.). P. Ory's work was supported by a joint PhD fellowship from the CNRS and the Région Poitou-Charente. The authors thank P. Lebaron and P. Catala for cytometry analyses, M. Bréret and F. Mornet for microplankton counts, B. Sautour for mesozooplankton counts, P. Pineau, N. Lachaussée and D. Vilday for their sampling assistance and H. J. Hartmann for improving the manuscript. Two anonymous reviewers are gratefully acknowledged for their helpful comments.