Armand Cavé-Radet

Sara Correa-Garcia

Cécile Monard

Abdelhak El Amrani

Armel Salmon

Malika Ainouche

Étienne Yergeau
email: etienne.yergeau@iaf.inrs.ca

Phenanthrene contamination and ploidy level affect the rhizosphere bacterial communities of Spartina spp

Keywords: Phenanthrene, rhizosphere, rhizoremediation, Spartina, polyploidy, bacterial communities

Spartina spp. are widely distributed salt marsh plants that have a recent history of hybridization and polyploidization. These events have resulted in a heightened tolerance to hydrocarbon contaminants, but the effects of this phenomenon on the rhizosphere microbial communities is unknown. Here, we grew two parental Spartina species, their hybrid and the resulting allopolyploid in salt marsh sediments that were contaminated or not with phenanthrene. The DNA from the rhizosphere soil was extracted and the bacterial 16S rRNA gene was amplified and sequenced, whereas the abundances of the genes encoding for the PAH ring-hydroxylating dioxygenase (RHD) of gram-negative and gram-positive bacteria were quantified by real-time PCR. Both the contamination and the plant genotype significantly affected the bacterial communities. In particular, the allopolyploid S. anglica, harbored a more diverse bacterial community in its rhizosphere. The interspecific hybrid and the allopolyploid also harbored significantly more copies of the PAH-RHD gene of gramnegative bacteria in their rhizosphere than the parental species, irrespective of the contamination treatments. Overall, our results are showing that the recent polyploidization events in the Spartina affected its rhizosphere bacterial communities, both under normal and contaminated conditions, possibly increasing its phytoremediation potential.

Introduction

Polycyclic aromatic hydrocarbons (PAHs) are a group of ubiquitous organic pollutants, that are highly concerning in view of their potentially severe impact on natural ecosystems and public health. The remediation of sites contaminated by these compounds due to human activity, such as oil spill, is mainly carried out through excavation, soil leaching or various techniques based on microbial degradation [START_REF] Samanta | Polycyclic aromatic hydrocarbons: Environmental pollution and bioremediation[END_REF][START_REF] Wilson | Bioremediation of soil contaminated with polynuclear aromatic hydrocarbons (pahs): A review[END_REF]. Phytoremediation is an alternative low-cost and environmentally friendly technology that uses plants and their associated microorganisms to remove pollutants from the environment. There are various types of phytoremediation, depending on the contaminant targeted and its fate [START_REF] Pilon-Smits | Phytoremediation[END_REF]. For instance, phytoextraction mainly concerns the extraction of inorganic pollutant from the soil and its translocation and accumulation in plant tissues. Volatile compounds can also be extracted from the soil by plants, and after translocation be released into the atmosphere during phytovolatilization.

However, one of the most interesting form of phytoremediation is rhizodegradation, where plants stimulate a wide variety of root-associated bacteria and fungi [START_REF] Correa-García | Rhizoremediation of petroleum hydrocarbons: A model system for plant microbiome manipulation[END_REF][START_REF] Ghosal | Current state of knowledge in microbial degradation of polycyclic aromatic hydrocarbons (pahs): A review[END_REF] to degrade organic contaminants, often all the way to CO2 (mineralization). Close interactions between plants and soil microorganisms were reported to play a major role in detoxification and metabolization of xenobiotics [START_REF] El Amrani | Omics" insights into pah degradation toward improved green remediation biotechnologies[END_REF][START_REF] Oliveira | Hydrocarbon contamination and plant species determine the phylogenetic and functional diversity of endophytic degrading bacteria[END_REF], and recent advances in 'omics' approaches now allow to address the plant-microbe complex in more details (Bell et al. 2014b). Root exudates were demonstrated to have a central role, modifying the microbial PAH degraders community structure [START_REF] Cebron | Root exudates modify bacterial diversity of phenanthrene degraders in pah-polluted soil but not phenanthrene degradation rates[END_REF][START_REF] Rohrbacher | Root exudation: The ecological driver of hydrocarbon rhizoremediation[END_REF] and significantly increasing microbial biomass [START_REF] Esperschutz | Influence of chronic ozone stress on carbon translocation pattern into rhizosphere microbial communities of beech trees (fagus sylvatica l.) during a growing season[END_REF]. Rhizodegradation of PAH was suggested to be a model for understanding and manipulating plant-microbe interactions [START_REF] Correa-García | Rhizoremediation of petroleum hydrocarbons: A model system for plant microbiome manipulation[END_REF].

A c c e p t e d m a n u s c r i p t

One area of interest that is understudied is the use of phytoremediation for the cleanup of coastal salt marshes contaminated by hydrocarbons coming from oil spills at sea. In that context, Spartina is a particularly interesting genus, colonizing salt marshes all around the world and providing key ecosystem services, while being particularly tolerant to the effects of oil spills [START_REF] Alvarez | Transcriptome response to the deepwater horizon oil spill identifies novel candidate genes for oil tolerance in natural populations of the foundation plant spartina alterniflora[END_REF][START_REF] Lin | Impacts and recovery of the deepwater horizon oil spill on vegetation structure and function of coastal salt marshes in the northern gulf of mexico[END_REF][START_REF] Lin | Response of salt marshes to oiling from the deepwater horizon spill: Implications for plant growth, soil surface-erosion, and shoreline stability[END_REF][START_REF] Ramanarao | Differential expression of the transcripts of spartina alterniflora loisel (smooth cordgrass) induced in response to petroleum hydrocarbon[END_REF][START_REF] Silliman | Degradation and resilience in louisiana salt marshes after the bp-deepwater horizon oil spill[END_REF]. This genus is also characterized by numerous hybridization and genome doubling events (polyploidy) [START_REF] Aïnouche | Hybridization, polyploidy and invasion: Lessons from spartina (poaceae)[END_REF], which are major evolutionary mechanisms for eukaryotes, and most specially for plants [START_REF] Soltis | Polyploidy and angiosperm diversification[END_REF][START_REF] Wendel | Genome evolution in polyploids[END_REF].

Recent in vitro work from our group has shown that genome doubling increased tolerance to phenanthrene in S. anglica as compared to its single-genome parents [START_REF] Cavé-Radet | Increased tolerance to organic xenobiotics following recent allopolyploidy in spartina (poaceae)[END_REF]. Indeed, the genomic shock induced by the merging of divergent genomes results in massive shifts in genetic and epigenetic pathways, that may lead to the emergence of new adaptive phenotypes to environmental constraints [START_REF] Comai | The advantages and disadvantages of being polyploid[END_REF]. Although previous studies have shown a link between plant genotype and the rhizosphere microbial community structure and gene expression during phytoremediation (Bell et al. 2014a;[START_REF] Yergeau | Soil contamination alters the willow rhizosphere and root metatranscriptomes and the root-rhizosphere interactome[END_REF], nothing is known about the effects of genome doubling events on the root-associated microorganisms and their response to PAH contamination.

Here, we hypothesized that the heightened tolerance of polyploid Spartina species to PAH contaminants observed in vitro also results in an increased abundance of PAH degraders within the root-associated microorganisms. We selected four Spartina species: the parental species Spartina alterniflora and S. maritima, their interspecific sterile hybrid S. x townsendii and the allopolyploid derivative S. anglica. These different Spartina species were grown in salt marsh sediments contaminated or not with phenanthrene, after which the bacterial 16S rRNA gene of the rhizosphere soil was amplified and sequenced, and the PAH ring-

Material and methods

Plant and soil material

Plants were collected in 2016 on their natural habitats along coastlines of France and England.

We selected four Spartina species: the parental species Spartina alterniflora (2n=6x=62) and S. maritima (2n=6x=60), their interspecific sterile hybrid S. x townsendii (2n=6x=62) that appeared and emerged at the end of the 19 th century and the allopolyploid derivative S. anglica that subsequently appeared following genome doubling (2n=12x=120, 122, 124). For the parental species, S. alterniflora was sampled at Le Faou (Roadstead of Brest, France), whereas S. maritima was sampled in Brillac-Sarzeau and Le Hezo (Gulf of Morbihan, France). The homoploid hybrid S. x townsendii, which distribution range is not extended to France was collected in Hythe (Romney Marsh, England), and the allopolyploid S. anglica was sampled at La Guimorais (Saint-Coulomb, France), where bulk sediments (not associated with plants) were also collected. The plant roots and rhizomes were washed abundantly with tap water before their introduction to the new sediment, in order to eliminate foreign plant residues, stones, and most of the rhizosphere sediments. After that, all the plants were acclimated in the La Guimorais sediments for three weeks before the start of the experiment.

During plant acclimation, the same sediments were used to prepare phenanthrene treated and untreated sediments for the experiment. For this, wet sediments were rinsed several times with tap water and sieved through a 5 mm sieve to reduce salt content and to remove plant residues and stones. We based our phenanthrene sediment contamination level (150 mg kg -1 ) on the work of [START_REF] Hong | A comprehensive study of the impact of polycyclic aromatic hydrocarbons (pahs) contamination on salt marsh plants spartina alterniflora: Implication for plant-microbe interactions in phytoremediation[END_REF], where it was reported that 100 mg kg -1 of phenanthrene was not toxic to S. alterniflora but affected plant-microorganism interactions. Previous studies have shown that the PAH concentration in the sediments of polluted rivers, estuaries or salt marshes around the world ranged from 2 mg kg -1 all the way to 1943 mg kg -1 [START_REF] Budzinski | Evaluation of sediment contamination by polycyclic aromatic hydrocarbons in the gironde estuary[END_REF][START_REF] Jiang | Characterization and distribution of polycyclic aromatic hydrocarbon in sediments of haihe river, tianjin, china[END_REF][START_REF] Shi | Contamination of rivers in tianjin, china by polycyclic aromatic hydrocarbons[END_REF][START_REF] Watts | Uptake of polycyclic aromatic hydrocarbons (pahs) in salt marsh plants spartina alterniflora grown in contaminated sediments[END_REF]. We selected a A c c e p t e d m a n u s c r i p t concentration of 150 mg kg -1 of phenanthrene, which represents a non-phytotoxic, moderate exposure in the range of concentrations reported in polluted environments. To contaminate the soil, 300 g of wet sediments was air-dried and then split in two 150 g subsamples. One 150 g subsample was spiked with 100 mM of phenanthrene (1.95 g of phenanthrene diluted in 109.4 mL of absolute ethanol), whereas the other 150 g subsample (control) was spiked with 109.4 mL of absolute ethanol. After total evaporation of ethanol for one day, these sediment samples were vigorously mixed by hand with an additional 12.85 kg of wet and rinsed sediments to reach a final concentration of 150 mg phenanthrene kg -1 substrate for the contaminated treatment. Some contaminated and control soils were sampled at this step and kept at -20°C to serve as T0 samples.

Experimental design

Both polluted and control sediments were supplemented with sterile vermiculite before the start of the experiment (volume 1/3) for a better substrate breathability. Individual Spartina plants (four species) were then transplanted in pot containing 1 kg of contaminated or control substrates, and rhizomes were carefully placed into rhizobags of 3 cm of diameter made from plastic mesh (1-cm mesh size) inspired from the experimental design of [START_REF] Hong | A comprehensive study of the impact of polycyclic aromatic hydrocarbons (pahs) contamination on salt marsh plants spartina alterniflora: Implication for plant-microbe interactions in phytoremediation[END_REF].

The rhizobags were filled with the exact same substrate as the remainder of the pot.

Rhizobags were used to create a physical separation between the bulk soil and the rhizosphere because of the difficulty to retrieve an intact root system when uprooting Spartina plants from very wet sediments. The rhizosphere was thus defined as the soil inside the rhizobags, where roots and rhizomes were also found. The experiment was replicated three times in a completely randomized design, resulting in a total of 24 pots that were placed in a phytotronic chamber with a light/dark regime of 16/8h, in an average ambient temperature of 20°C. Pots were watered with 200 mL of half-strength Hoagland's nutrient solution (Hoagland and HotSart ReadyMix, 0.4 mg.ml -1 bovine serum albumin, 0.6 µM forward and reverse primers, and 1 µL of DNA. PCR conditions were as follows: 5 min of initial denaturation (95°C), 25 cycles of 30s denaturation (95°C), 30s for primer annealing (55°C), and 45s elongation (72°C), followed by final extension for 10 min at 72°C. PCR amplicons were then purified using AMPure XP beads. Amplicon indexing was carried out by a second PCR in 25 µL final reaction volumes containing 12.5 µL 2X KAPA HiFi HotSart ReadyMix, 2.5 µL of specific index primers 1 and 2 from the Nextera Index kit, and 5 µL of purified amplicon. The second step PCR conditions were as follows: 3 min of initial denaturation (95°C), 8 cycles of 30s denaturation (95°C), 30s for primer annealing (55°C), and 45s elongation (72°C), followed by final extension for 5 min at 72°C. Index PCR cleanup was performed using AMPure XP beads. Finally, PCR amplicons were quantified using PicoGreen, normalized at 1 ng µl

Real-time quantitative PCR

The qPCR experiments were conducted on a Stratagene Mx3005P qPCR system (Agilent Technologies), associated with the corresponding software MxPro Mx3005P (v4.10; Agilent).

The qPCR reactions were performed using the primers designed by [START_REF] Cébron | Real-time pcr quantification of pah-ring hydroxylating dioxygenase (pah-rhdα) genes from gram positive and gram negative bacteria in soil and sediment samples[END_REF]. The qPCR reactions were performed in 20µl total volume containing 1x iTaq universal SYBR® Green reaction mix supplemented with 300µM of each primer (Integrated DNA Technologies™) and 5µl of DNA template at a concentration around 5-10ng µL -1 or water (negative control). The amplifications were carried out following the protocol provided in [START_REF] Cébron | Real-time pcr quantification of pah-ring hydroxylating dioxygenase (pah-rhdα) genes from gram positive and gram negative bacteria in soil and sediment samples[END_REF] with some modifications. Briefly, the first step consisted on heating to

Statistical analyses

Statistical analyses were performed in R (v 3.5.1). Inversed Simpson index was calculated after rarefaction through random subsampling based on the size of the smallest library (14,671 sequences). Comparisons between α-diversity of bacterial communities, the abundance of the PAH-RHD genes, the relative abundance of phyla and the relative abundance of selected putative hydrocarbon degraders according to the Spartina species and the treatments (polluted or control substrate) were conducted using ANOVA with post-hoc Tukey HSD tests (aov and tukeyHSD functions of the stats package). Prior to ANOVA the normality, and the homoscedasticity of the data to be tested were confirmed using the shapiro.test and bartlett.test of the stats package, respectively. Spearman correlations were calculated using the cor.test of the stats package. Principal coordinate analysis (PCoA) to describe β-diversity was performed on normalized OTU tables using Bray-Curtis dissimilarity index in the vegan A c c e p t e d m a n u s c r i p t

Results

Bacterial alpha diversity and community structure

The presence of phenanthrene led to a general decrease of bacterial alpha diversity in the rhizosphere of all the plant species (two-way ANOVA: F = 12.41, P = 0.002) (Fig. 1a). The two-way ANOVA also highlighted a significant effect of plant species on the bacterial alpha diversity (F = 6.19 P = 0.005), and Tukey HSD post-hoc tests showed that the rhizosphere of S. anglica was significantly more diverse than the rhizosphere of S. alterniflora and S. maritima (Fig. 1a). The interaction effect contamination × species was not significant for the bacterial alpha diversity in two-way ANOVA tests. PCoA ordinations based on bacterial OTUs showed some level of clustering based on plant species and contamination (Fig. 1b).

Permanova tests confirmed that Spartina species (F = 2.09, R 2 = 0.22, P < 0.001) and phenanthrene (F = 2.48, R 2 = 0.09, P < 0.001) had highly significant influences on the bacterial community structure.

Bacterial community composition

The relative abundance of various dominant bacterial phyla in the rhizosphere of Spartina varied between the different plant species and between the contaminated and noncontaminated soils (Fig. 2). These shifts were not significant for some phyla, such as the Acidobacteria, the Actinobacteria, and the Bacteroidetes. In contrast, plant species significantly influenced the rhizosphere relative abundance of Chloroflexi (F = 6.21, P = 0.05, significantly higher in the rhizosphere of S. anglica as compared to S. alterniflora and S. x townsendii in Tukey HSD tests), Planctomycetes (F = 8.86, P = 0.00108, significantly higher in the rhizosphere of S. anglica as compared to S. alterniflora and S. x townsendii, and in the rhizosphere of S. maritima vs. S. alterniflora) and total Proteobacteria (F = 4.54, P = 0.02, significantly lower in the rhizosphere of S. anglica as compared to S. alterniflora and S. x The summed relative abundance of all these genera varied between around 2% to over 3% of all reads and was significantly higher in the rhizosphere of plants growing in contaminated soil (F = 8.43, P = 0.01, Fig. 3), but showed no significant differences between plant species.

For the individual genera, the relative abundance of Paenibacillus, Bacillus, Burkholderia, Acinetobacter, Alcanivorax, Rhodococcus and Pseudomonas did not vary significantly between the plant species and treatments. However, the relative abundances of Sphingomonas (F = 9.85, P < 0.001, significantly higher in the rhizosphere of S. x townsendii as compared to all other species), Sphingobacterium (F = 5.41, P = 0.009, significantly higher in the rhizosphere of S. maritima as compared to all other species), Nocardia (F = 7.56, P = 0.002; significantly higher in the rhizosphere of S. anglica as compared to all other species) showed significant differences between the rhizosphere of the different plant species. The relative abundances of Massilia (F = 17.54, P < 0.001), Cycloclasticus (F = 4.96, P = 0.04) and

Mycobacterium (F = 7.49, P = 0.01) were significantly affected by the phenanthrene treatment, with significantly higher relative abundance in the rhizosphere of plants growing in the contaminated soil (Fig. 3). In addition, the interaction between contamination and plant species was significant for Sphingobacterium (F = 4.73, P = 0.001). Indicator species analyses confirmed some of the trends observed in ANOVA tests, as it identified bacterial OTUs related to Massilia and Cycloclasticus as the OTUs with the highest indicator power for contaminated rhizospheres among all OTUs (Table 1 and Supplementary Table S1). The relative abundance

Abundance of PAH ring-hydroxylating dioxygenases

The abundance of the gene encoding for the alpha subunit of the PAH-RHD of gram-positive bacteria was not significantly different between the rhizosphere of the four Spartina species and was not affected by sediment contamination, nor by their interaction (Fig. 4a). In contrast, the abundance of the gram-negative version of the same gene showed significant differences between the four Spartina species (F = 63.2, P = 4.31 × 10 -9 ), with the allopolyploid S.

anglica and the hybrid S. x townsendii harboring on average 354 % more copies of the gene in their rhizospheres than the parental species (Fig. 4b). Contamination and the contamination × species interaction did not have any significant effect on the abundance of the gram-negative version of the PAH-RHD.

A c c e p t e d m a n u s c r i p t

Discussion

Recently, we have shown enhanced in vitro tolerance to phenanthrene in the allopolyploid in S. anglica as compared to the parental species [START_REF] Cavé-Radet | Increased tolerance to organic xenobiotics following recent allopolyploidy in spartina (poaceae)[END_REF]. The hypothesis behind the present pot study was that this increased tolerance of the allopolyploid S. anglica to hydrocarbon contamination would also result in an increased abundance of PAH degraders within its root-associated microorganisms. The results presented here are coherent with this hypothesis, as we found significant differences between the different pot-grown Spartina species in terms of their bacterial community composition, structure and diversity, and, more importantly, in the abundance of PAH ring-hydroxylating dioxygenase (RHD) genes in their rhizosphere. The allopolyploid S. anglica harbored a more diverse bacterial community, composed of relatively more Nocardia, Chloroflexi and Planctomycetes and less Proteobacteria in its rhizosphere as compared to its diploid parents and their hybrid. In addition, the allopolyploid (S. anglica) and the hybrid (S. x townsendii) had significantly more PAH-RHD genes related to gram-negative bacteria in their rhizosphere than the parental species, suggesting a heightened capacity for PAH degradation by the rhizosphere microbial community. Previous studies from our group using willows have shown that the phylogeny of the host plant significantly influenced the fungal community composition under high levels of contaminant (Bell et al. 2014a). We also showed that the metatranscriptomic response of the rhizosphere microbial communities to contamination varied between different willow genotypes, and that this response was mirrored in the growth of the genotypes in contaminated soil [START_REF] Yergeau | Soil contamination alters the willow rhizosphere and root metatranscriptomes and the root-rhizosphere interactome[END_REF]. However, this is the first time, to our knowledge, that functional and taxonomical differences between the rhizosphere bacterial communities of recently naturally speciated plants are reported in the context of soil/sediment contamination.

The rhizosphere bacterial communities of the four pot-grown Spartina species tested showed significant responses to the presence of the contaminant. In the presence of PAH, salt

A c c e p t e d m a n u s c r i p t marsh plants were reported to harbor different microbial communities, favoring hydrocarbondegrading microorganisms [START_REF] Ribeiro | Hydrocarbon degradation potential of salt marsh plant-microorganisms associations[END_REF]. Similarly, we found here an increase in the relative abundance of putative PAH degraders, and more specifically for the genera Mycobacterium, Cycloclasticus and Massilia. These genera represent large bacterial groups that are able to metabolize and degrade PAH [START_REF] Johnsen | Principles of microbial pah-degradation in soil[END_REF][START_REF] Kanaly | Biodegradation of high-molecular-weight polycyclic aromatic hydrocarbons by bacteria[END_REF] and are consistent with previous results about the PAH-degrading bacteria associated to Spartina in salt marshes [START_REF] Daane | Isolation and characterization of polycyclic aromatic hydrocarbon-degrading bacteria associated with the rhizosphere of salt marsh plants[END_REF][START_REF] Launen | Characterization of the indigenous pah-degrading bacteria of spartina dominated salt marshes in the new york/new jersey harbor[END_REF]. Interestingly, Cycloclasticus spp. were often reported as one of the predominant PAH degraders in seawater [START_REF] Tremblay | Chemical dispersants enhance the activity of oil-and gas condensate-degrading marine bacteria[END_REF]Yergeau et al. 2015b), whereas Mycobacterium are typical soil PAH degraders [START_REF] Hennessee | Effects of polycyclic aromatic hydrocarbon mixtures on degradation, gene expression, and metabolite production in four mycobacterium species[END_REF] and Massilia were found to degrade PAHs in soils or associated to roots [START_REF] Liu | Colonization on root surface by a phenanthrene-degrading endophytic bacterium and its application for reducing plant phenanthrene contamination[END_REF][START_REF] Wang | Efficient biodegradation of phenanthrene by a novel strain massilia sp. Wf1 isolated from a pah-contaminated soil[END_REF]. This diversity of putative PAH degraders associated with the rhizosphere of Spartina might be linked to the nature of its habitat, at the interface of terrestrial, plant and marine ecosystems. These bacterial genera probably contain candidates of interest for salt marsh remediation, but their potential would have to be confirmed.

However, these taxonomical shifts in the relative abundance of these putative microbial hydrocarbon degraders were not reflected in the abundance of the PAH-RHD genes.

Indeed, contamination with phenanthrene did not shift the abundance of PAH-RHD genes in the rhizosphere of the pot-grown Spartina, both for the gram-positive and -negative variants of the gene. Previous studies had alluded to a similar phenomenon, where genes and microorganisms related to the degradation of hydrocarbon contaminants were already enriched in the rhizosphere of plants growing in non-contaminated soil [START_REF] Daane | Isolation and characterization of polycyclic aromatic hydrocarbon-degrading bacteria associated with the rhizosphere of salt marsh plants[END_REF].

These hydrocarbon degrading genes were also expressed by microbial isolates in response to root exudates [START_REF] Mark | Transcriptome profiling of bacterial responses to root exudates identifies genes involved in microbe-plant interactions[END_REF][START_REF] Matilla | Genomic analysis reveals the major driving forces of bacterial life in the rhizosphere[END_REF][START_REF] Ramachandran | Adaptation of rhizobium leguminosarum to pea, alfalfa and sugar beet rhizospheres investigated by comparative transcriptomics[END_REF], and found in the rhizosphere metatranscriptome of willows growing in not-contaminated soils (Yergeau et [START_REF] Yergeau | Soil contamination alters the willow rhizosphere and root metatranscriptomes and the root-rhizosphere interactome[END_REF]. One of the explanations put forward was the presence of secondary metabolites in the rhizosphere environment that were structurally analog to hydrocarbon contaminants [START_REF] Singer | Secondary plant metabolites in phytoremediation and biotransformation[END_REF], resulting in an enrichment and activation of PAH degradation genes even in the absence of exogenous PAH.

Alternatively, the significant enrichment of PAH-RHD GN genes in the rhizospheres of the hybrid and the allopolyploid could be due to a larger total bacterial community as recently reported for Solidago polyploids [START_REF] Wu | Polyploidy in invasive solidago canadensis increased plant nitrogen uptake, and abundance and activity of microbes and nematodes in soil[END_REF]. Regardless of the underlying mechanism, the higher abundance of PAH-RHD genes in the rhizosphere of the hybrid and the allopolyploid suggest a heightened capacity for contaminant tolerance and degradation. Of course, this higher abundance of PAH-RHD genes does not mean that they were actively expressed, and a study based on the quantification of transcripts is needed to complete this observation.

The discrepancy between the qPCR of functional genes and the relative abundance of putative hydrocarbon degraders could be explained by two reasons: 1) unknown microbial groups could harbour these genes, since part of soil microbial diversity is still unexplored; and

2) it is possible that not all members of the genera singled out carried out PAH-RHD degradation genes, as intra-genus and intra-species variability in that trait is expected. Recent studies showed that bacterial isolates undistinguishable based on their 16S rRNA gene harbored different gene pools, had different ecological behaviors and functional capacity [START_REF] Irshad | Bacterial subspecies variation and nematode grazing change p dynamics in the wheat rhizosphere[END_REF][START_REF] Lopes | Genome variations between rhizosphere and bulk soil ecotypes of a pseudomonas koreensis population[END_REF]. These results are suggesting that one should exert caution when trying to extrapolate functional information from amplicon sequencing data.

As bacteria are thought to be the major players in organic contaminant degradation during rhizoremediation (Bell et al. 2014a;Bell et al. 2014b;[START_REF] Correa-García | Rhizoremediation of petroleum hydrocarbons: A model system for plant microbiome manipulation[END_REF][START_REF] El Amrani | Omics" insights into pah degradation toward improved green remediation biotechnologies[END_REF], and that recent plant-microbe metatranscriptomic studies confirmed that A c c e p t e d m a n u s c r i p t the hydrocarbon degradation genes expressed in the root-rhizosphere environment were mostly linked to bacteria [START_REF] Gonzalez | Trees, fungi and bacteria: Tripartite metatranscriptomics of a root microbiome responding to soil contamination[END_REF][START_REF] Yergeau | Soil contamination alters the willow rhizosphere and root metatranscriptomes and the root-rhizosphere interactome[END_REF], it is interesting to note the higher bacterial diversity in the rhizosphere of the allopolyploid S. anglica as compared to the parental Spartina species. Upon a contamination event, the plant-associated microbial diversity could have a central importance, as higher diversity generally results in functional redundancy, which would enable the microbial communities to cope with a wider variety of environmental conditions while still providing essential services to the plant. Previous studies from our group have shown that the initial soil diversity explained better the difference in willow growth under highly contaminated conditions than diversity at the end of the experiment (Yergeau et al. 2015a). It was further suggested in that study that restoring microbial diversity of degraded environments could be the key for successful phytoremediation. This potential link between the higher rhizosphere diversity and the increased resilience of S. anglica to contamination is intriguing and warrants further research.

Even though we acclimated the plants in the sediments for three weeks, we cannot exclude the possibility that the differences in the microbial communities in the rhizosphere of S. anglica were because the bulk soil used for the experiment was from a site where S. anglica predominantly grew. This might also have affected indirectly the microbial communities through shifts in the physiology of the parental plants and of the hybrid that were adapted to the conditions found in other sediments. However, some commonalities in the responses of the allopolyploid and of the hybrid suggest that the effect of the soil origin on some of the results presented here was not necessarily that strong.

We described here for the first time the root-associated bacterial community of four pot-grown Spartina species in a context of allopolyploidization and PAH contamination.

Significant differences were observed between the plant species and between contaminated and control rhizospheres. The most salient differences were the higher abundance of gram-A c c e p t e d m a n u s c r i p t negative PAH-RHD genes and the higher bacterial diversity in the rhizosphere of the allopolyploid as compared to the parental species, which could be a major factor contributing to the increased tolerance of S. anglica to contaminant stress. The underlying cause for these shifts could be related to changes in the root exudation patterns following hybridization, but more work would be necessary to test this hypothesis.

  s c r i p t hydroxylating dioxygenase (PAH-RHD) genes of gram-positive (GP) and gram-negative 82 (GN) bacteria were quantified by real-time quantitative PCR.

  Arnon 1950) every five days. After 60 days of growth, plants were uprooted and the rhizosphere soil (inside the rhizobags) was collected and stored at -20°C until DNA extraction.Bacterial 16S rRNA gene amplificationGenomic DNA from soil samples was extracted using a MoBio PowerSoil DNA extraction kit. Amplicons were prepared from total DNA by PCR targeting the bacterial 16S rRNA gene[START_REF] Caporaso | Ultra-high-throughput microbial community analysis on the illumina hiseq and miseq platforms[END_REF]. PCR amplifications were performed in 25 µL final reaction volumes containing final concentrations of 1X KAPA HiFi

  ) followed by 40 cycles of denaturation at 95°C for 30s, annealing at either 57°C (PAH-RHD GN) or 54°C (PAH-RHD GP) for 35s and elongation at 72°C for 75s, after which SYBR Green fluorescence was measured. At the end of the run, a melting curve analysis was performed where fluorescence was measured at 0.5°C temperature increment every 5s from 51 to 95°C. Standard curves for the two genes were made from 10-fold dilutions of linearized plasmid containing the gene fragment of interest, cloned from amplified soil DNA[START_REF] Yergeau | Microarray and real-time pcr analyses of the responses of high arctic soil bacteria to hydrocarbon pollution and bioremediation treatments[END_REF]).

(

  [START_REF] Dixon | Vegan, a package of r functions for community ecology[END_REF] and ape[START_REF] Paradis | Ape: Analyses of phylogenetics and evolution in r language[END_REF] packages. Permanova testing the impact of Spartina species, treatments, and their interactions were conducted through 999 permutations. Indicative species among OTUs were investigated using default parameters provided by the indicspecies package[START_REF] Cáceres | Associations between species and groups of sites: Indices and statistical inference[END_REF].

  A significant effect of contamination was also observed for the Chloroflexi (F = 20.93, P < 0.001), total Proteobacteria (F = 14.23, P = 0.002), Verrucomicrobia (F = 4.85, P = 0.43) and Gammaproteobacteria (F = 11.39, P = 0.004). The genera Sphingobacterium, Acinetobacter, Nocardia, Pseudomonas, Mycobacterium, Burkholderia, Bacillus, Sphingomonas, Rhodococcus, Paenibacillus, Massilia, Alcanivorax, Cycloclasticus were singled out as putative hydrocarbon degraders based on a survey of the available literature.

of

  Massilia was significantly and negatively correlated to bacterial alpha diversity (inverse Simpson index) (rs = -0.48, P = 0.017). Apart from a positive correlation for Rhodococcus (rs = 0.41, P = 0.047), the other genera were not significantly correlated to bacterial alpha diversity.

Table 1 .

 1 Top 10 strongest indicator bacterial OTUs for the phenanthrene-contaminated Spartina rhizospheres.

	OTU	stat	P-value Taxonomy
	Otu00188	0.998 0.001	Massilia
	Otu00140	0.996 0.001	Cycloclasticus
	Otu00370	0.991 0.001	Alcaligenaceae unclassified
	Otu00552	0.989 0.049	Flavobacteriaceae unclassified
	Otu02955	0.960 0.002	Salinirepens
	Otu01642	0.954 0.006	Sphingobium
	Otu02528	0.954 0.005	Alcaligenaceae unclassified
	Otu01089	0.941 0.005	Peredibacter
	Otu04988	0.921 0.008	unassigned
	Otu03519	0.896 0.006	Gammaproteobacteria unclassified

Acknowledgements A.C.-R. was supported by an Erasmus+ mobility scholarship from the European Commission. This work was supported by the Ministère de l'Enseignement Supérieur et de la Recherche, by the CNRS, the Observatoire des Sciences et de l'Univers de Rennes (OSUR) and by a NSERC Discovery Grant (2014-05274 to EY). This research was enabled in part by support provided by Calcul Québec (www.calculquebec.ca) and Compute Canada (www.computecanada.ca).