D R Stacy

A Krueger-Hadfield

Dr Myriam Valero

Guido Bonthond
email: gbonthond@geomar.de

Stacy A Krueger-Hadfield

Francisco Rafael Barboza

Masahiro Nakaoka

Myriam Valero

Gaoge Wang

Sven Künzel

Florian Weinberger

How do microbiota associated with an invasive seaweed vary across scales?

Keywords: microbiota, microbiome, haplodiplontic, holobiont, geographic scales, global diversity patterns, invasive species

Communities are shaped by scale dependent processes. To study the diversity and variation of microbial communities across scales, the invasive and widespread seaweed Agarophyton vermiculophyllum presents a unique opportunity. We characterized pro-and eukaryotic communities associated with this holobiont across its known distribution range, which stretches over the northern hemisphere. Our data reveal that community composition and diversity in the holobiont vary at local but also larger geographic scales. While processes acting at the local scale (i.e., within population) are the main structuring drivers of associated microbial communities, changes in community composition also depend on processes acting at larger geographic scales.

Interestingly, the largest analyzed scale (i.e., native and nonnative ranges) explained variation in the prevalence of predicted functional groups, which could suggest a functional shift in microbiota occurred over the course of the invasion process. While high variability in microbiota at the local scale supports A.

vermiculophyllum to be a generalist host, we also identified a number of core taxa. These geographically independent holobiont members imply that co-introduction of specific microbiota may have additionally promoted the invasion process.

Introduction

Virtually all multicellular organisms, including seaweeds, are holobionts [START_REF] Margulis | Words as battle cries: symbiogenesis and the new field of endocytobiology[END_REF], but redefined in Bordenstein & Theis 2015) and interact with complex and variable microbiota. Unlike terrestrial plants, seaweeds lack a sophisticated root system and utilize their resources from or through the surrounding medium (i.e., seawater). The community on the algal surface, mainly composed of microbes, but including multicellular eukaryotes as well, occupies a delicate niche. It colonizes the interface between organism and environment and forms a barrier which both threats (e.g., pathogens) and resources (e.g., light) need to cross to reach the host [START_REF] Wahl | The second skin: ecological role of epibiotic biofilms on marine organisms[END_REF]Egan et al., 2013). Red macroalgae are known in particular to manipulate the composition of associated epibiota by the production of metabolites (e.g., Weinberger & Friedlander 2000;[START_REF] Harder | Chemical mediation of ternary interactions between marine holobionts and their environment as exemplified by the red alga Delisea pulchra[END_REF][START_REF] Saha | Microbial "gardening" by a seaweed holobiont: Surface metabolites attract protective and deter pathogenic epibacterial settlement[END_REF]. The associated community structure is also dependent on local microbial availability, which in turn is shaped by interactions within the microbial community and environmental variables. In other words, microbial communities associated with seaweeds and other holobionts are shaped by processes that act on different ecological scales, such as the chemistry of the host and environmental conditions (e.g., salinity, temperature and spatial scales with which such physical parameters vary). While species abundance and distribution patterns are controlled by scale-dependent processes [START_REF] Mcgill | Ecology. Matters of scale[END_REF], microbial diversity patterns across scales have only

Accepted Article

This article is protected by copyright. All rights reserved recently received interest (e.g., [START_REF] Martiny | Drivers of bacterial beta-diversity depend on spatial scale[END_REF]Tedersoo et al., 2014;[START_REF] Locey | Scaling laws predict global microbial diversity[END_REF][START_REF] Lindström | Local and regional factors influencing bacterial community assembly[END_REF].

The invasive seaweed Agarophyton vermiculophyllum (Ohmi) Gurgel et al. (commonly known under the synonym Gracilaria vermiculophylla) represents a unique opportunity to study variation in microbiota across different ecological scales. Like most Rhodophyta, this alga has a haplodiplontic lifecycle in which haploid male gametophytes fertilize haploid female gametophytes, giving rise to diploid tetrasporophytes that in turn produce male and female gametophytes through meiosis. Both gametophytes and the tetrasporophytes develop into independent adults that are morphologically only distinguishable by reproductive structures [START_REF] Krueger-Hadfield | Invasion of novel habitats uncouples haplo-diplontic life cycles[END_REF]. Agarophyton vermiculophyllum originates from the Northwest Pacific, where its distribution range stretches from near the tropics (21.4° N; [START_REF] Hu | A unique genetic lineage at the southern coast of China in the agar-producing Gracilaria vermiculophylla (Gracilariales, Florideophyceae)[END_REF] to high latitudes up to ~52° N in the cold temperate zone [START_REF] Skriptsova | Red algae of the genus Gracilaria in the south of the Russian Far East[END_REF]. Over the course of the invasion process the distribution range expanded to the Eastern Pacific (Bellorin, Oliveira & Oliveira 2004), Eastern Atlantic [START_REF] Rueness | Life history and molecular sequences of Gracilaria vermiculophylla (Gracilariales, Rhodophyta), a new introduction to European waters[END_REF] and Western Atlantic (Thomsen, Gurgel, Fredericq & McGlathery 2006), where it presently also covers considerable latitudinal ranges (Krueger-Hadfield et al., 2017 and references therein). This alga thrives in a diversity of environments that differ in temperature, salinity, nutrient availability, sedimentation, grazing and epiphytic overgrowth (reviewed in [START_REF] Hu | Adaptation mechanisms and ecological consequences of seaweed invasions: a review case of agarophyte Gracilaria vermiculophylla[END_REF]. As an ecosystem engineer it has profound consequences on the local environment and can for example increase habitat availability and faunal biomass (Byers, Gribben, Yeager & Sotka 2012), change ecosystem functions [START_REF] Ramus | An invasive foundation species enhances multifunctionality in a coastal ecosystem[END_REF] compete with native ecosystem engineers (Hammann, Buchholz, Karez & Weinberger 2013), alter nutrient cycling (Gonzalez, Smyth, Piehler & McGlathery 2013) and affect commercial fishing (Freshwater et al., 2006). The distribution range covers local to global geographic scales, across which different processes influence the associated microbial community. For instance, different micro-environments within one single algal individual may harbor different communities (e.g., [START_REF] Morrissey | Disentangling the influence of environment, host specificity and thallus differentiation on bacterial communities in siphonous green seaweeds[END_REF] which at the same time depend on environmental parameters acting at local (e.g., Campbell, Marzinelli, Gelber & Steinberg 2015) and regional (e.g., [START_REF] Lindström | Local and regional factors influencing bacterial community assembly[END_REF] scales. Microbial availability and connectivity may still be important at the coastal or continental scale (e.g., [START_REF] Martiny | Drivers of bacterial beta-diversity depend on spatial scale[END_REF][START_REF] Sunagawa | Ocean plankton. Structure and function of the global ocean microbiome[END_REF] and effects related to the invasion process, such as genetic diversity of the host or potential adaptations, are likely more relevant at continental and global scales (e.g., Arnaud-Haond et al., 2017). Since the nonnative range is spatially far larger and the invasion process has been facilitated by multiple introductions [START_REF] Krueger-Hadfield | ) Genetic identification of source and likely vector of a widespread marine invader[END_REF] it is probable that one or multiple selection processes have occurred, both in the native and nonnative ranges. Although such adaptions have not been identified genetically, several lines of evidence corroborate that populations from the nonnative range tend to be more tolerant to several stressors, including grazing (Hammann et al., 2016a), temperature and salinity extremes (Hammann, Wang, Boo, Aguilar-Rosas & Weinberger 2016b;[START_REF] Sotka | Combining niche shift and population genetic analyses predicts rapid phenotypic evolution during invasion[END_REF] and epiphytic overgrowth (Wang et al., 2017a;Wang et al., 2017b).

Accepted Article

This article is protected by copyright. All rights reserved Experimental work has indicated that A. vermiculophyllum populations from native and nonnative ranges are chemically better defended against epiphytic bacterial settlers from their own range [START_REF] Saha | Rapid adaptation to controlling new microbial epibionts in the invaded range promotes invasiveness of an exotic seaweed[END_REF]. This does not only suggest that A. vermiculophyllum is (chemically) well-equipped to manipulate epibiota, but also implies that the defense is fine-tuned to the environment at some spatial scale. The synthesis of molecules that interact with the microbial community is plastic, as the production or increase in production may also depend on external conditions. These traits might importantly aid a seaweed host in maintaining an epiphytic community that is overall beneficial. Whether such plasticity, increased plasticity through adaptation, or more specific adaptations have contributed to its success as an invader remains unclear. However, the abovementioned studies strongly suggest that interactions with the associated microbial community have been important in the invasion process. Whether changes in the interaction between host and microbes (such as shown in Wang et al., 2017a;Wang et al., 2017b;[START_REF] Saha | Rapid adaptation to controlling new microbial epibionts in the invaded range promotes invasiveness of an exotic seaweed[END_REF], reflect a single adaptation of the host that occurred during the invasion, or rather adaptation or acclimation acting on more local scales is an important question to which the answer may importantly contribute to our understanding of the invasion success of A. vermiculophyllum and other invasive species.

While specific host-microbe interactions between A. vermiculophyllum and epibiota have been studied [START_REF] Saha | Microbial "gardening" by a seaweed holobiont: Surface metabolites attract protective and deter pathogenic epibacterial settlement[END_REF]Wang et al., 2017a;Wang et al., 2017b;[START_REF] Saha | Rapid adaptation to controlling new microbial epibionts in the invaded range promotes invasiveness of an exotic seaweed[END_REF], the composition of the epiphytic and endophytic microbial communities and variation therein across the native and nonnative ranges and smaller scales has not been characterized. It is unknown what commonly the most abundant microbial endoand epiphytic taxa are and whether A. vermiculophyllum hosts core microbiota (i.e., a subset of taxa that is persistently associated, independent of the geographic origin; [START_REF] Shade | Beyond the Venn diagramm: the hunt for a core microbiome. Beyond the Venn diagram: the hunt for a core microbiomeemi_2585[END_REF]. Characterizing the putative core microbiota and how associated microbial communities vary across scales is crucial to identifying the ecological processes and symbionts that are relevant to the invasion process.

To this end, we collected individuals from fourteen populations and characterized microbiota using highthroughput amplicon sequencing of the prokaryotic 16S-V4 and the eukaryotic 18S-V7 ribosomal DNA regions aiming to answer the following two primary questions: (i) Do microbial communities differ between the algal surface (epibiota) and the internal tissue (endobiota) micro-environments? (ii) How do epiphytic communities vary across geographic scales? Furthermore, we expanded subsets of the sampling with a second year and the collection of haploid individuals, to address two additional questions: (iii) How do microbial communities vary between years within populations? (iv) Do microbial communities differ among host ploidies and sexes (sensu [START_REF] Hughes | Ecology and the Evolution of Biphasic Life Cycles[END_REF]?

Methods

Sampling design

Accepted Article

This article is protected by copyright. All rights reserved Since A. vermiculophyllum spores require hard substratum to settle and germinate and adults lack the ability to produce a new holdfast [START_REF] Krueger-Hadfield | Invasion of novel habitats uncouples haplo-diplontic life cycles[END_REF], we specifically collected individuals that were fixed to hard substratum indicating spore recruitment to avoid sampling fragments recently separated from the same individual or clonally reproduced drifting individuals. To avoid confounding possible effects of life cycle stage (i.e., ploidy and/or sex) with other effects related to the geographic scales, we identified life cycle stages and compared micro-environments and geographic scales using only diploid individuals, processing six diploids per collection site. To sample groups that can be considered populations, only sites where individuals that were presumably sexually reproductive were visited (i.e., sites known to include diploids and haploid females and males fixed to hard substrata; [START_REF] Krueger-Hadfield | Invasion of novel habitats uncouples haplo-diplontic life cycles[END_REF]. We selected the collection sites such that the overall design was organized into five nested scales, which in ascending hierarchy were: (i) individuals, (ii) populations, (iii) ecoregions (as defined in [START_REF] Spalding | Marine ecoregions of the world: a bioregionalization of coastal and shelf areas[END_REF], (iv) continental coasts, and (v) ranges (native and nonnative; Figure 1). To distinguish local from regional effects, all populations were sampled in pairs, separated by approximately 100 km. Three ecoregions from the Asian coast (Yellow Sea, Oyashio Current and Northeastern Honshu), two from the European coast (North Sea and Celtic Seas) and one from the Atlantic and Pacific coasts of the North American continent (Virginian and Northern California, respectively) were visited. The highest hierarchical scale (i.e., range) included one native and three nonnative continental coasts (Figure 1).

Collection and sample processing

Agarophyton vermiculophyllum specimens were collected from 14 September to 21 September 2016 and 24 August to 23 September 2017. At each site at least ten individuals and three water samples of 50 mL were taken. Algae were sampled with gloves, directly stored in new and separate freezer storage bags in a cooling box until arrival at the laboratory where they were immediately stored at 4 ⁰C until processing. All processing occurred within a maximum of 12 hours after collection.

One gram of each alga was taken from the terminal part of the thallus (including branches with apices) and transferred to a 50 mL tube containing 15 sterile glass beads of 4 mm in 7.5 mL of artificial seawater. The artificial seawater was prepared from sterile distilled water and 24 gL -1 NaCl. To separate the epiphytic community, we followed the same method as [START_REF] Saha | Rapid adaptation to controlling new microbial epibionts in the invaded range promotes invasiveness of an exotic seaweed[END_REF] with slight modifications. Tubes were placed on a vortexer for 3 min after which the water was filtered through a 0.2 µL PCTE filter with a vacuum pump.

Subsequently, a new volume of 7.5 mL sterile distilled water was added, vortexed for 3 min and filtered through the same filter. Water samples were filtered directly. Filters were stored in absolute ethanol until DNA extraction. From the remaining algal tissue, a small fragment (2-4 cm) was preserved in DESS solution (25% DMSO, 2.5M EDTA and NaCl saturated) for DNA extraction of endobiota. One additional fragment was preserved in silica gel for ploidy identification. Diploids were identified morphologically with a dissecting

Accepted Article

This article is protected by copyright. All rights reserved microscope before sample processing or by post-hoc microsatellite genotyping, following the amplification and genotyping methods described in [START_REF] Krueger-Hadfield | Invasion of novel habitats uncouples haplo-diplontic life cycles[END_REF]. In brief, DNA was extracted from the thalli preserved in silica and ploidy was determined using ten polymorphic microsatellite loci [START_REF] Kollars | ) Development and characterization of microsatellite loci for the haploid-diploid red seaweed Gracilaria vermiculophylla[END_REF].

We considered an individual diploid if at least one locus was heterozygous [START_REF] Krueger-Hadfield | Invasion of novel habitats uncouples haplo-diplontic life cycles[END_REF].

Based on either morphological or microsatellite identification of life cycle stage, six diploids were selected for amplicon sequencing. In order to test differences in associated microbiota among life cycle stages, we collected at least thirty individuals from the four North American populations and identified diploids, haploid females and haploid males based on the presence of tetrasporangial sori (diploid tetrasporophytes), cystocarps (female gametophytes) or spermatangial sori (male gametophytes) using a dissecting microscope (see [START_REF] Krueger-Hadfield | Everywhere you look, everywhere you go, there's an estuary invaded by the red seaweed Gracilaria vermiculophylla (Ohmi) Papenfuss[END_REF]. From these sites, six diploids and five of each haploid stage (i.e., sexes) were processed for characterizing microbial communities. Finally, from the four populations collected in 2016, ploidies and sexes were not identified (morphologically nor molecularly) and the microsatellite method was unsuccessful for one of the Japanese populations from 2017 (Futatsuiwa). Therefore, instead of six diploids all ten individuals were included in the sequencing for the populations from 2016 and the Futatsuiwa population from 2017 (Table S1).

DNA extraction and sequencing to characterize microbial communities

In the laboratory, ethanol was evaporated from preserved filters with a vacuum centrifuge at 30 °C for one to two hours. The tissue fragments were taken from the DESS solution and rinsed two times with DNA free water. We note that it is possible that certain firmly attached epibiota could have remained associated with the surface throughout this treatment and that differences between endo-and epiphytic communities may be underestimated. DNA was extracted with the ZYMO Fecal/soil microbe kit (D6102; ZYMO Research, Irvine, CA, USA) following the manufacturer's protocol. The 16S-V4 region was amplified with the primers 515F (S-*-Univ-0515-a-S-19) and 806R (S-D-Arch-0786-a-A-20; see [START_REF] Klindworth | Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies[END_REF] and 18S-V7 region with F-1183mod and R-1443mod (Ray et al., 2016). To prepare amplicon libraries, we applied the two-step PCR strategy from Gohl et al. (2016), using the KAPA HIFI HotStart polymerase (Roche, Basel, Switzerland) and the same indexing primers. The procedure for the 18S-V7 fragment was identical, except for the target primer sequences. Amplicon concentrations from the second PCR were estimated from gel pictures and amplicons were pooled accordingly into 16S and 18S libraries which were purified by a gel extraction step, combined in a 5:1 ratio and sequenced at the Max-Plank-Institute for Evolutionary Biology (Plön, Germany) on the Illumina 2x300 MiSeq platform. Libraries included three negative controls from the DNA extractions, three negative PCR controls and two mock communities as positive control (HM-782D, Bei Resources, Manassas, VA, USA).

The fastq files were demultiplexed (0 mismatches) after which sequence assemblage and quality filtering was performed using Mothur software (v1.40.5; [START_REF] Schloss | Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities[END_REF] and the SILVA alignment [START_REF] Quast | The SILVA ribosomal RNA gene database project: improved data processing and web-based tools[END_REF] release 132) to denoise and classify the full set of reads. Unique sequences were clustered into OTUs with the

Accepted Article

This article is protected by copyright. All rights reserved opticlust algorithm based on the traditional 3% dissimilarity criterion. Mitochondrial, chloroplast, eukaryotic and unclassified sequences were removed from the 16S dataset. We also removed sequences from the 18S dataset that were unclassified, classified as Bacteria, Archaea, higher plants and the genus of the host (Gracilaria in the SILVA 132 release). Since we were primarily interested in epiphytes (that is, sessile taxa growing on the surface of the alga) and sessile animals are not known as prevalent epifauna of A. vermiculophyllum (Nyberg, Thomson & Wallentinus, 2009), sequences classified to Animalia were also discarded from the 18S dataset. While in the present work referred to as the eukaryotic microbiota, this constituent of the community also contains many multicellular epiphytes. Finally, OTUs that were singletons in the full dataset and samples with less than a 1,000 16S sequence counts and samples with less than a 100 18S sequence counts were also excluded from downstream analyses. Raw de-multiplexed amplicon reads and metadata were deposited in the SRA database (accession: PRJNA564581).

Functional profiling

Functional profiles were obtained with tax4fun software (Aßhauer, Wemheuer, Daniel & Meinicke 2015) which uses the SILVA reference alignment and Kyoto encyclopedia of Genes and Genomes (KEGG; [START_REF] Kanehisa | Data, information, knowledge and principle: back to metabolism in KEGG[END_REF] to predict metagenomics profiles. Since the package was not yet compatible with the Silva 132 release at the time of analysis, OTUs were reclassified in Mothur using the Silva 123 release. Aiming to characterize general functional groups within the communities we selected genes essential to autotrophy (ribulose-1,5biphosphate carboxylase oxygenase; K01601), aerobic heterotrophy (cytochrome c oxidase subunit III; K02276) and nitrogen fixation (nitrogenase iron protein; K02586). We also included anaerobic heterotrophs as a functional group by combining the mean gene count values for sulfate reduction (adenylylsulfate reductase; K00394), denitrification (methane/ammonia monooxygenase; K10944), acetogenesis (formate-tetrahydrofolate ligase; K01938) and fumarate respiration (fumarate reductase; K00244), which we presumed constitute the majority of anaerobic respiration metabolisms.

Data analysis

For a total of 273 samples the 16S-V4 and 18S-V7 regions were sequenced and the models explained below were applied to both datasets. Importantly, to adjust for the effect of differences in sequencing depth, inherent to amplicon sequencing based community data, we included the natural logarithm of the sequencing depth (LSD) as a continuous variable in all multi-and univariate models. To reduce complexity and computation time, datasets were trimmed for all multivariate analyses to the 95% most abundant OTUs. Sequence counts were expressed in proportion to the total sample count (relative counts) and sorted in descending order after which cumulative percentages were computed. To compare community abundance-based compositions, multivariate

Accepted Article

This article is protected by copyright. All rights reserved generalized linear models (mGLMs) were fitted with the manyglm() function from the R package mvabund [START_REF] Wang | mvabund-an R package for model-based analysis of multivariate abundance data[END_REF]Warton, Wright & Wang 2012).

The mGLMs assumed a negative binomial distribution with a natural logarithm as link function. In addition to providing a deviance and p-value at the multivariate level (the 'community response'), the mGLM computes the same statistics at the univariate level (each OTU) which we used to count OTUs that were differentially abundant across factors of interest (p unadjusted < 0.05). Furthermore, by obtaining both uni-and multivariate outputs simultaneously, we could use the mGLMs to compare communities and characterize core microbiota coherently using a single model. To confirm whether mGLMs satisfied the distribution and mean-variance assumptions, QQ-plots and residuals vs. fitted-plots were visually inspected.

To address the first question -whether microbial communities are different among micro-environments -we tested the alternative hypothesis that microbial abundance-based community composition varies among microenvironments. Since mixed models are currently not available in the mvabund package for mGLMs, we used a two-step modelling approach. To first adjust the data for geographic effects we included the factor 'population identity' in the first model along with the LSD. On the residuals of this first model we ran a second model with micro-environment (including water, surface and tissue) as response variable and a Gaussian distribution (the residuals from this first model were in the link function and therefore normally distributed). The p-value of the overall effect of micro-environment was obtained using the anova.manylm() function, using 500 bootstrap iterations. We obtained the between-group comparisons and corresponding p-values of the levels within microenvironment by running summaries with different reference levels using the summary.manylm() function and also 500 bootstrap iterations (see Table S2a for the full output). Multivariate dissimilarities were visualized with non-metric dimensional scaling (nMDS, with the R package vegan; [START_REF] Oksanen | Package 'vegan'. Community ecology package, version 2[END_REF], using the raw data and rescaled residuals from the mGLMs.

OTU richness, evenness and predicted functionality were compared with mixed linear models using the lme4 R package (Bates, Mächler, Bolker & Walker 2015). Since richness and evenness are the variables of which diversity is essentially a function, we used both parameters as indicators of structural differences in communities. We calculated asymptotic richness (S chao , based on the chao1 estimator, Chao & Bunge 2002) and rarified richness (S n ). For S n all samples were rarified to the read count of the sample with the lowest sequencing depth (1,008 for prokaryotes and 110 for eukaryotes). From a strictly theoretical point of view, both S chao and S n are expected to be independent of the sequencing depth (but the accuracy of the asymptotic estimate should increase with more reads) and are not necessarily correlated. We used the probability of interspecific encounter (PIE; Hurlbert 1971) as a measure of evenness [START_REF] Mcglinn | Measurement of Biodiversity (MoB): A method to separate the scale-dependent effects of species abundance distribution, density, and aggregation on diversity change[END_REF]. To compare functional group abundances between micro-environments and variation in group abundances across geographic scales the same model structures were used. LSD and micro-environment were included as fixed variables and population and individual identity as random variables. To meet the assumption of normality PIE

Accepted Article

This article is protected by copyright. All rights reserved was logit transformed and the gene counts representing the predicted functional group abundances were log transformed. Assumptions on the distribution were verified and p-values for differences among water, epi-and endobiota were obtained by post-hoc TukeyHSD tests (see Table S2b for the full statistical output). Marginal and conditional R 2 values (variation explained by fixed effects and fixed plus random effects, respectively) were computed with the r.squaredGLMM function from the lme4 R package [START_REF] Nakagawa | A general and simple method for obtaining R2 from generalized linear mixed-effects models[END_REF].

To characterize how variation in epibiotic communities is partitioned across geographic scales, we used the subset of data representing diploid algal surface samples, comprising all fourteen populations. We fitted mGLMs including the LSD and the nested structure of spatial scales (range/coast/region/population) as explanatory variables. OTU richness, PIE and predicted functional group abundances were analyzed using mixed linear models (Table S2c). The LSD was included as fixed effect and the nested structure of spatial scales as random effect, resulting in the following model structure: ~ LSD + (1|range/coast/region/population), based on the approach used by [START_REF] Messier | How do traits vary across ecological scales? A case for trait-based ecology[END_REF]. We extracted variance estimates and calculated the corresponding 95% confidence intervals by bootstrapping the models with a 1,000 iterations.

Multivariate GLMs with the same model structure were also fitted on subsets of the data containing only the water samples to verify the assumption that environmental microbiota vary across geographic scales (see Table S2a).

The four populations collected in both 2016 and 2017 were first analyzed with an mGLM. The lowest spatial scale (population identity) was included to represent spatial effects, under which the year was nested.

Multivariate distances were visualized in dendrograms using a hierarchical clustering approach for which we used the residuals from an mGLM with LSD as single response variable to adjust for the effect of sequencing depth. The clustering was conducted with the R package pvclust, resampling with multiscale bootstrapping and a 1,000 iterations (Suzuki & Shimodaira 2006). The nested structure was specified as random in univariate models to analyze how much variation was explained by the year after explaining the geographic variation (see Table S2d).

To answer the question whether communities vary among host ploidies and sexes we tested the alternative hypothesis that different lifecycle stages host different microbial communities, using the four North American populations, for which all lifecycle stages had been collected. The mGLMs included the variables LSD, microenvironment (from which the water level was removed), life cycle stage and population identity. Based on the lowest AICsum we did not consider the interaction between micro-environment and lifecycle stage in the mGLMs for both pro-and eukaryotic datasets. In the univariate models 'population' and 'individual identity' were specified as random variables while 'LSD', 'micro-environment' and 'lifecycle stage' were fixed. As the interaction between micro-environment and lifecycle stage yielded a lower AICc for richness estimates, we retained the interaction in all univariate models (see Table S2e for statistical output and Table S2f for model comparisons).

Accepted Article

This article is protected by copyright. All rights reserved Since the concept of 'cores' is somewhat flexible [START_REF] Shade | Beyond the Venn diagramm: the hunt for a core microbiome. Beyond the Venn diagram: the hunt for a core microbiomeemi_2585[END_REF], we characterized core microbiota based on two approaches, applying an occurrence criterion (occurrence core) and a model-based composition criterion (composition core). Both approaches were performed on the subset of our data containing the six populations for which tissue and surface samples were available. We considered OTUs detected in strictly 100% of the water, algal surface or tissue samples as members of the water, epi-and endophytic occurrence cores. Those present in 100% of both tissue and surface samples were considered as part of the algal core. The composition core was determined based on the univariate results of the mGLMs applied on the same dataset to compare microbial communities between micro-environments. As these mGLMs were run on data that had been adjusted for spatial effects (by using the residuals from mGLMs including population identity as explanatory variable) we considered OTUs that were significantly more abundant (p unadjusted < 0.05) to be micro-environment specific and geographically independent. OTUs fulfilling this criterion were identified using the coefficients and p-values from the mGLM summaries and classified accordingly to either water, surface, tissue or alga (no difference between surface and tissue but in both more abundant than in the water).

We explored correlations between OTUs within the holobiont following the joint modeling approach (Warton et al., 2015). Under the assumption that the residuals have been adjusted for the effects of micro-environments (water, surface and tissue) and population identity (and all physical-chemical variables confounded in them), these correlations would represent the variation driven by biological variables and could thus hint at patterns of co-occurrence or mutual exclusion among microbes. The untransformed residuals of pro-and eukaryotic mGLM models (normally distributed in the scale of the link function) were concatenated and all Pearson correlation coefficients and corresponding p-values were computed between the 25 most abundant pro-and eukaryotic OTUs. Correlations with p-values > 0.05 were discarded.

Results

Sequencing summary

A total of 4,537,292 prokaryotic reads and 24,793 OTUs remained from 40 water, 42 endo-and 147 epiphytic samples after removal of low quality sequences, singletons and 16S samples with less than a 1,000 reads.

Samples ranged from 1008 to 112,921 in read counts with a median of 13,534. Of the eukaryotic dataset 576,509 sequences and 1,919 OTUs from 37 water, 40 endo-and 138 epiphytic samples remained after quality filtering and removal of 18S samples with less than a 100 reads. Samples ranged from 110 to 15,883 in read counts with a median of 1,868. The three most abundant prokaryotic OTUs were classified to the genera Pleurocapsa (Cyanobacteria) and Granulosicoccus (γ-Proteobacteria) and to the family Rhodobacteraceae (α-Proteobacteria). The three most abundant eukaryotic OTUs were classified to the family Bacillariophyceae

Accepted Article

This article is protected by copyright. All rights reserved (Diatomea), the order Ectocarpales (Phaeophyceae) and the genus Navicula (Diatomea). The most abundant phyla were Proteobacteria, Bacteriodetes and Cyanobacteria. While in epiphytic communities, Bacteriodetes was the second most abundant and Cyanobacteria were third, in endophytic communities this was reversed.

The most abundant phylum in both endo-and epiphytic eukaryotic communities was the Ochrophyta (including Diatomea) and was followed by Florideophycidae. While in endophytic communities the third most abundant phylum was Chlorophyta, Peronosporomycetes (Oomycetes) was the third phylum in epiphytic communities (Figure S1). The removal of low read count samples reduced the number of individuals (initially 6) in a few populations (see Table S3).

(i) Are microbial communities different among micro-environments?

For both pro-and eukaryotic communities, the null-hypothesis that there are no differences in abundancebased community composition among micro-environments was rejected (p pro and p euk < 0.01). Post-hoc comparisons, using summaries with different reference levels resolved significant p-values for all pairwise combinations (p pro and p euk < 0.01). These differences are also apparent in the nMDS plots (Figure 2a andb). Overall, 28.9% and 29.6% of 3,054 and 226 prokaryotic and eukaryotic OTUs were differentially abundant, respectively (Figure 3a). Also for the diversity metrics the null-hypothesis was rejected (p pro and p euk < 0.01).

Asymptotic and rarefied OTU richness of epibiota were higher than of endobiota in both pro-and eukaryotic communities. Post-hoc pairwise comparisons showed that in prokaryotic communities, richness of the water and epibiotic communities was similar, while richness of eukaryotic communities was higher in the water than both endo-and epiphytic communities (Figure 2c, e, only asymptotic richness shown, see Table S2b for the statistical output). The sequencing depth had a significant effect on S chao in both pro-and eukaryotic communities. Despite S chao being an estimate of the asymptotic richness (species count of the overall theoretical community), this is not unexpected for sequencing based community data (e.g., [START_REF] Smith | Sequence depth, not PCR replication, improves ecological inference from next generation DNA sequencing[END_REF]. The sequencing depth had no effect on the rarified richness (p pro = 0.30, p euk = 0.73) or the evenness parameter PIE (which is mainly depending on the shape of the rarefaction curve; [START_REF] Mcglinn | Measurement of Biodiversity (MoB): A method to separate the scale-dependent effects of species abundance distribution, density, and aggregation on diversity change[END_REF]; p pro = 0.34, p euk = 0.24). Differences in predicted functional group abundances were found for autotroph, anaerobic heterotroph and diazotroph groups (p < 0.01). Whereas in auto-and diazotroph groups, counts increased from the water to epiphytic communities and from epi-to endophytic communities, predicted anaerobic heterotrophy was highest in the water and decreased in epi-and endophytic communities. The null-hypothesis was, however, not rejected for aerobic heterotrophy (p = 0.48; Figure 2g-j).

(ii) How do epiphyte communities vary across geographic scales?

Differences in epiphyte community composition were found across all geographic scales (p range = 0.03, p coast = 0.03, p region < 0.01, p population < 0.01 for the prokaryotic communities and p range < 0.01, p coast < 0.01, p region < 0.01, p population < 0.01 for eukaryotic communities). Differentially abundant OTUs ranged from 35.9% to 69.5% of 3251

Accepted Article

This article is protected by copyright. All rights reserved prokaryotic OTUs and 36.4% to 64.9% of 231 eukaryotic OTUs, with for both datasets more differentially abundant OTUs at lower spatial scales (Figure 3b). A similar pattern can be observed from nMDS plots (Figure S2) where between-level differences are more pronounced at lower scales. Also across scales the sequencing depth was only significant for S chao (p LSD < 0.01 in both pro-and eukaryotic models) and explained considerable variation (47.0% in prokaryotic and 27.9% in eukaryotic communities). In line with the multivariate results, the smallest geographic scale (population) was the most important, explaining variation ranging from 27.5% prokaryotic S chao to 61.7% in prokaryotic S n (Figure 4). The highest scale (range) was estimated to explain a minor fraction of variation in prokaryotic S chao and S n (4.59% and 4.16% percent, respectively). For eukaryotic S n , coast was estimated to explain 18.1% and range to explain 2.94% in eukaryotic PIE. Most variation in predicted functional group abundance was explained at the population level with percentages ranging from 25.2% (phototrophs) to 60.7% (anaerobic heterotrophs; without variation explained by LSD). Predicted autotroph and diazotroph abundances also varied at the highest hierarchical level (range) and explained 42.1% of the variation in autotrophs and 33.0% in diazotrophs (Figure 4). Within the water, pro-and eukaryotic microbiota varied compositionally at all scales (p < 0.01) and this was most pronounced at the population level where a number of OTUs (106 pro-and 11 eukaryotic) differed significantly (Table S2c, Figure S2i).

(iii) How do microbial communities vary between years?

The mGLMs resolved significant p-values for both population (representing the geographic effects) and year (in all cases p < 0.01). However, in prokaryotic communities geographic effects were relatively more pronounced with 64.9% of 2,849 OTUs differentially abundant among populations compared to 14.0% between years. In the eukaryotic communities 65.4% of 179 OTUs were differentially abundant among populations and within populations 26.8% between years (Figure 3c). In the hierarchical clustering of the prokaryotic OTU table all individuals clustered within their respective populations with high bootstrap values (≥ 99, Figure S3). Within the population clusters, with some exceptions, samples from different years clustered adjacently. For the eukaryotic data, the overall pattern was similar but with less bootstrap support for the nodes combining both years of the same populations and a few more individuals clustering within different populations. The effect of the higher geographic scales detected with the subset of the data including all 14 populations was not visible in these dendrograms. Variation in prokaryotic OTU richness and evenness was explained majorly by geographic effects (S chao = 42.9%, S n = 63.6% and PIE = 65.6%), while variation explained by year within population was little to none (S chao = 0%, S n = 4.2% and PIE = 5.7%). In eukaryotic communities, however, more variation was explained by the inter-annual changes (Figure S4, Table S2e).

(iv) Are microbial communities different among host ploidies and sexes?

Accepted Article

This article is protected by copyright. All rights reserved

In the mGLMs for prokaryotes and eukaryotes life cycle stage explained the least deviance of all included variables and was only significant for eukaryotic communities (p eukaryote = 0.02). However, while in this case the alternative hypothesis was not rejected, post-hoc comparisons did not resolve any significant pairwise differences (see Table S2a). The hypothesis that S chao , S n and PIE differ among life cycle stages was in all cases rejected. The interaction between life cycle stage and micro-environment was not significant either in any of the models (Figure S5).

Core microbiota and correlations

For each of the micro-environments (i.e., water, surface and tissue) a number of prokaryotic OTUs was present in 100% of the samples from the six populations from which both endo-and epiphytic samples were sequenced (Table S4). Based on the occurrence approach we counted two endophytic, five epiphytic, seven algal (present in 100% endo-and epiphytic samples) and two water core OTUs, which were all prokaryotes. The composition core counted 141 endophytic, 26 epiphytic, 123 algal and 166 water prokaryotic OTUs and one endophytic, 6 epiphytic, 6 algal and 32 eukaryotic OTUs (Figure 3a, Table S4). After removing the effects of geography and micro-environment, the correlations among the 25 most abundant pro-and eukaryotic OTUs with p-values < 0.05 were mostly positive for within domain comparisons (within prokaryotes 61 positive and 7 negative correlations and within eukaryotes 19 positive and 10 negative correlations). Between domain correlations were majorly negative (28 positive and 34 negative correlations; Figure S6).

Discussion

Endo-and epiphytic microbial communities

The tissue and surface of A. vermiculophyllum harbor specific communities that differ in terms of composition, diversity and predicted functionality from each other and from the surrounding water. This is perhaps not surprising as distinct endo-and epibiota were also recognized in green algae [START_REF] Morrissey | Disentangling the influence of environment, host specificity and thallus differentiation on bacterial communities in siphonous green seaweeds[END_REF]Aires, Moalic, Serrão & Arnaud-Haond 2015). Overall, epiphytic communities are characterized by higher richness (two-fold or more) and evenness compared to endophytic communities. The decreased richness and evenness in the tissue compared to the surface may reflect that the surface is more exposed to the environment and more prone to community coalescence [START_REF] Rillig | Interchange of entire communities: microbial community coalescence[END_REF]. Further, chemical control by the host might be higher within the tissue and the biofilm on the algal surface could act as a biotic filter that prevents potential symbionts from reaching the endophytic microbial niche. Although richness and evenness are merely diversity parameters, both have been associated with functional community properties, including productivity (Bell, Newman, Silverman, Turner & Lilley 2005) and community stability (Wittebolle et al., 2009;Coyte, Schluter & Foster 2015). Differences in the predicted functional groups between tissue and surface suggest that endophytic communities are more autotrophic and diazotrophic, while epiphytic communities are more

Accepted Article

This article is protected by copyright. All rights reserved heterotophic (anaerobically). It is noteworthy that one of the most abundant OTUs, which was also a member of the algal core in both the occurrence and compositional approaches, was classified to the Pleurocapsa (Cyanobacteria; Figure S1a). Speculatively, such cyanobacteria might be the main group of potentially autotrophic endophytes and may even be able to fix inorganic nitrogen. Many terrestrial plants (e.g., legumes) are known to form mutualistic relationships with diazotrophic bacteria and therefore have an advantage in nitrogen limiting environments. For seaweeds, such associations have only been found in few occasions [START_REF] Rosenberg | Nitrogen fixation by blue-green algae associated with the siphonous green seaweed Codium decorticatum: effects on ammonium uptake[END_REF]Gerard, Dunham & Rosenberg 1990;[START_REF] Raut | Diazotrophic macroalgal associations with living and decomposing Sargassum[END_REF]. Some However, specialized structures such as nodules like in M. faya have not been identified in A. vermiculophyllum and future experimental work is required to determine the functional activity of these prevalent Cyanobacteria and other members of endo-and epiphytic cores.

Core microbiota

Within the high variability of microbial communities, core microbiota (or core microbiomes) represent a signature of stability. Different concepts are used to define microbial community cores based on occurrence, composition, persistence or connectivity (see [START_REF] Shade | Beyond the Venn diagramm: the hunt for a core microbiome. Beyond the Venn diagram: the hunt for a core microbiomeemi_2585[END_REF]. For occurrence cores, a range of percentage thresholds has been applied (e.g., 12% to 100%, Astudillo-García et al., 2017), and OTUs may be collapsed into higher taxonomic levels. We applied a highly conservative 100% occurrence threshold without collapsing OTUs. Since the caveat of an occurrence core is that the probability of detecting core taxa is directly correlated with abundance we also defined a composition core with a model based approach, using the univariate results of the mGLMs comparing composition between micro-environments. With this approach we detected substantially more core taxa (303 compared to 14 associated to tissue, surface or alga). While this may be less strict compared to the occurrence criterion, the model controls for effects from sampling location and the increased abundances of these OTUs endophytically, epiphytically or in the alga compared to the water, thus representing a geographically stable signature within the high amount of variation.

Despite the considerable geographic distribution, with both occurrence and composition based criteria a number of core taxa was detected. This contrasts somewhat with recent studies on seaweed associated microbiota in the green algal genera Caulerpa and Ulva (Burke, Steinberg, Rusch, Kjelleberg & Thomas 2011;[START_REF] Roth-Schulze | Functional biogeography and host specificity of bacterial communities associated with the Marine Green Alga Ulva spp[END_REF]but see Arnaud-Haond et al., 2017). In these algae functional -but not taxonomiccores were detected and it was suggested that the assemblage of microbial communities in seaweeds might be shaped by stochasticity and abiotic environmental filtering, selecting for functions rather than specific taxa

Accepted Article

This article is protected by copyright. All rights reserved (Burke et al., 2011;[START_REF] Roth-Schulze | Functional biogeography and host specificity of bacterial communities associated with the Marine Green Alga Ulva spp[END_REF]. In contrast, a recent study comparing microbiota associated with the kelp Ecklonia radiata showed that the microbial community composition was more dependent on the host-condition rather than the local and regional environments [START_REF] Marzinelli | Continental-scale variation in seaweed host-associated bacterial communities is a function of host condition, not geography[END_REF]. For red algae, taxonomic and functional cores have, to our knowledge, not been characterized, nor have associated microbiota been compared across multiple geographic scales. Our results may reflect that the environment (i.e., a composite of the host and the local environment) filters taxa more strongly in A. vermiculophyllum than in green algae that have been studied previously. However, function and taxonomy are correlated and the filtering of specific functions might indirectly also filter for certain taxa and vice versa.

The detected core microbes are either globally available taxa or represent associations that predate the invasion process and originate from the native range. The co-introduction of microbes with specific functions in the holobiont is thought to be an important facilitator of successful invasion (Arnaud-Haond et al., 2017;[START_REF] Rodríguez-Echeverría | Jack-of-all-trades and master of many? How does associated rhizobial diversity influence the colonization success of Australian Acacia species?[END_REF]. Only few of the detected core microbiota were found in the water column and A. vermiculophyllum's core may thus largely constitute co-introduced microbes.

The ability of an invasive seaweed to maintain specific taxa with beneficial functions might provide an advantage over competitors and perhaps even protect against disadvantageous microbiota from the environment that are unable to settle in the already populated niche. As shown in [START_REF] Saha | Microbial "gardening" by a seaweed holobiont: Surface metabolites attract protective and deter pathogenic epibacterial settlement[END_REF], A. vermiculophyllum is associated with cultivatable microorganisms that provide it with protection from microbial pathogens. At the same time the host produces metabolites that selectively promote the settlement of these protective bacteria. Such mechanisms may be general and directed at taxonomic or functional groups or be specific and targeted at microbes that are part of the core microbiota. However, to test such hypotheses it would be necessary to disentangle the filtering by different environmental components from one another (i.e., that from the local abiotic and biotic environments and that of the host itself).

Variation in epibiota across scales

Given that previous studies indicated that interactions between host and epibiota might have been involved in the invasion process of A. vermiculophyllum [START_REF] Saha | Rapid adaptation to controlling new microbial epibionts in the invaded range promotes invasiveness of an exotic seaweed[END_REF]Wang et al., 2017a;Wang et al., 2017b), we were specifically interested in epiphytic communities from a geographic perspective. While microbial communities have been studied across local and regional scales in other seaweeds (e.g., [START_REF] Lindström | Local and regional factors influencing bacterial community assembly[END_REF][START_REF] Marzinelli | Continental-scale variation in seaweed host-associated bacterial communities is a function of host condition, not geography[END_REF][START_REF] Roth-Schulze | Functional biogeography and host specificity of bacterial communities associated with the Marine Green Alga Ulva spp[END_REF], in this study we also compared variation in microbiota across more global geographic scales, reaching to that of a hemisphere. To a certain extent our results are in line with previous studies [START_REF] Marzinelli | Continental-scale variation in seaweed host-associated bacterial communities is a function of host condition, not geography[END_REF][START_REF] Roth-Schulze | Functional biogeography and host specificity of bacterial communities associated with the Marine Green Alga Ulva spp[END_REF], which suggest that seaweed microbiota strongly depend on the local scale as we found for both pro-and eukaryotic datasets that most variation in OTU richness and evenness is explained among populations. In addition, abundances of all considered predicted functional groups varied mostly at the local scale. However, some functional groups

Accepted Article

This article is protected by copyright. All rights reserved (autotrophs and diazotrophs) varied also substantially between ranges and the overall community composition differed among levels of all geographic scales. This indicates that epibiotic communities, in terms of composition, diversity and function, are mostly defined at the local scale, but also shows that processes acting at larger scales contribute to shaping the microbial community. What these processes are and from where they originate (the environment, the host or an interaction between environment and host) is important to identify as they may represent general mechanisms of relevance to other microbiota or holobionts. While compositionally all scales are important, in terms of predicted function we found that only the most local (population) and the most global scale (the native and nonnative ranges) mattered. On this scale, the invasion process itself is most relevant and may thus hint at a functional shift in the epibiotic community that occurred during the process.

Given the strong genetic signature of the host itself between the ranges [START_REF] Krueger-Hadfield | ) Genetic identification of source and likely vector of a widespread marine invader[END_REF], it is tempting to speculate that genetic differentiation may be an important process acting at this scale, underlying functional differences within epibiota. However, to test this hypothesis, the confounding genetic and geographic patterns need to be disentangled. To achieve this, an experimental approach is required where different populations are translocated or compared in the context of a common garden.

In addition to comparing geographic scales, we collected four populations in two consecutive summers and were able to characterize how these populations varied between years in terms of pro-and eukaryotic epibiota (Figure S3 and S4). Prokaryotic communities were relatively similar between years and these results suggest a temporally stable signature within the A. vermiculophyllum holobiont. Our result is in line with Lachnit, Meske, [START_REF] Lachnit | Epibacterial community patterns on marine macroalgae are host-specific but temporally variable[END_REF], where denaturing gradient gel electrophoresis was used to characterize A. vermiculophyllum associated epibiota. Although based on only twelve samples from a single location [START_REF] Lachnit | Epibacterial community patterns on marine macroalgae are host-specific but temporally variable[END_REF] observed that A. vermiculophyllum communities differed between summer and winter, but returned to similar compositions during the summer and winter of the subsequent year. While this pattern was especially clear within prokaryotic communities in our data, richness and evenness were less stable between years within eukaryotic communities. Therefore, seasonal stability in the structure of the associated community might be considerably lower for eukaryotic epibiota and the presence of these taxa may be less specific and more sensitive to stochasticity and local environmental conditions. The results from [START_REF] Lachnit | Epibacterial community patterns on marine macroalgae are host-specific but temporally variable[END_REF] and this study emphasize the need to also consider temporal scales. Important environmental variables such as temperature that can strongly impact microbial communities [START_REF] Sunagawa | Ocean plankton. Structure and function of the global ocean microbiome[END_REF]) vary both at geographic (e.g., with latitude) and temporal scales (e.g., with season). However, the effect of temperature on microbial communities is scale-dependent and may be unambiguous [START_REF] Nottingham | Microbes follow Humboldt: temperature drives plant and soil microbial diversity patterns from the Amazon to the Andes[END_REF] or not [START_REF] Hendershot | Consistently inconsistent drivers of microbial diversity and abundance at macroecological scales[END_REF]. In the case of holobionts it is also important to identify processes originating from the host that affect diversity and composition of associated microbiota. To identify, disentangle and classify the importance of these abiotic and biotic processes, more scales must be considered and microbial communities should be studied along natural gradients in for example temperature, salinity and latitudinal space.

Accepted Article

This article is protected by copyright. All rights reserved Agarophyton vermiculophyllum life cycle stages are characterized by morphological and ecological differences [START_REF] Krueger-Hadfield | Invasion of novel habitats uncouples haplo-diplontic life cycles[END_REF]. For the haplodiplontic life cycle to be stabilized over evolutionary timescales, ecological differentiation is essential [START_REF] Hughes | Ecology and the Evolution of Biphasic Life Cycles[END_REF]. However, we found no differences in associated microbiota among reproductive diploids, haploid females and haploid males. This could indicate that traits related to chemical control of endo-and epiphytes are conserved across life cycle stages and the dominance of diploids in the nonnative range [START_REF] Krueger-Hadfield | Invasion of novel habitats uncouples haplo-diplontic life cycles[END_REF] may thus be related to physiological traits (e.g., thallus strength, [START_REF] Lees | Nonnative Gracilaria vermiculophylla tetrasporophytes are more difficult to debranch and are less nutritious than gametophytes[END_REF] rather than traits associated to manipulation of microbiota. However, we sampled during a time of the year when populations tend to be healthy and with substantial biomass (Weinberger, Buchholz, Karez & Wahl 2008;[START_REF] Muangmai | Seasonal fluctuation in a marine red alga, Gracilaria vermiculophylla (Gracilariales, Rhodophyta), from Nokonoshima Island, Southern Japan[END_REF] and we cannot predict whether microbiota remain similar among lifecycle stages under stressful conditions.

Holobiont and invasion biology

Communities are structured by processes that are scale dependent (McGill 2010 and references therein). In the case of communities associated with a holobiont these processes may originate from the host, the local environment in which the host resides or the combined environment of host and local conditions. Evidently, a microbial community that is benign, or even beneficial, contributes to the health and success of the host. An ability to exert a certain degree of control over the associated community could therefore be a valuable strategy. Adaptations to associate specific beneficial symbionts could have occurred over time and if such (core) microbes accompany the host to nonnative ranges, they may provide an advantage and promote invasiveness (i.e., accompanying mutualist hypothesis; [START_REF] Rodríguez-Echeverría | Rhizobial hitchhikers from Down Under: invasional meltdown in a plantbacteria mutualism[END_REF]. Given the tremendous environmental variability of microbes it might be advantageous to adapt simultaneously a more flexible strategy in associating beneficial microbial taxa and invest energy in mechanisms targeting more general groups (generalist host hypothesis; [START_REF] Rodríguez-Echeverría | Jack-of-all-trades and master of many? How does associated rhizobial diversity influence the colonization success of Australian Acacia species?[END_REF]. While maintaining essential associations by hosting certain microbes as core members, a high degree of flexibility towards microbial associations may thus be of additional benefit for an invasive species. Arnaud-Haond et al. (2017) identified core microbiota from the invasive green alga Caulerpa taxifolia which they showed were associated with (or had accompanied) the host across native and nonnative ranges. Our data show high variability at the local scale and suggest that A. vermiculophyllum is a generalist host. However, the substantial number of core microbiota we found implies that this alga was also accompanied over the course of the invasion. Agarophyton vermiculophyllum may therefore qualify for both a generalist and accompanied host. Although originally based on legumes and rhizobial symbionts, [START_REF] Rodríguez-Echeverría | Jack-of-all-trades and master of many? How does associated rhizobial diversity influence the colonization success of Australian Acacia species?[END_REF] classified such accompanied generalist hosts as species with a high probability to invade, and with a typically fast invasion process. Given its current wide latitudinal distribution ranges along the nonnative coasts and in some cases rapid regional expansion (e.g., [START_REF] Nyberg | Flora and fauna associated with the introduced red alga Gracilaria vermiculophylla[END_REF]Byers et al., 2012), the invasion process of A. vermiculophyllum may well

Accepted Article

This article is protected by copyright. All rights reserved have been accelerated by accompanying symbionts. However, to characterize the functions of these core microbiota and to assess whether they have promoted the invasion, the A. vermiculophyllum holobiont needs to be studied in an experimental context.

Conclusions

This is the first study addressing seaweed microbiota with a sampling design at the scale of the known distribution range. The results presented here show that A. vermiculophyllum associated microbial communities differ within the individual between surface and tissue in composition, diversity and predicted functionality. Our data also reveal that A. vermiculophyllum hosts microbial core taxa, some of which are specific to the tissue or the surface, while others are associated to both niches. These taxa are conserved across the distribution range and may represent functionally important symbionts. Variation in epibiota is compositionally partitioned across all local and global scales, but mostly at the local scale (i.e., among populations), at which OTU richness and evenness varied most substantially. At the scale of the invasion (i.e., native and nonnative ranges), however, differences in community composition and predicted functional groups could suggest that a shift in epibiota occurred over the course of the invasion process. Finally, with microbiota that are locally highly variable and a geographically conserved core A. vermiculophyllum matches criteria of both the generalist host and the accompanying mutualist hypotheses, which may in part explain its remarkable success as an invasive species.

Accepted Article

This article is protected by copyright. All rights reserved numbers corresponding to the biogeographic framework from [START_REF] Spalding | Marine ecoregions of the world: a bioregionalization of coastal and shelf areas[END_REF]. The insets are ordered by coast and show the sampling sites (labelled as in Table S1) inside the corresponding ecoregions. Sites from where only epibiotic samples were taken are labeled in green and those from which also endobiotic samples were acquired are green-red. Populations collected in both 2016 and 2017 are indicated by squares. Symbols in the lower right corner indicate which ploidies and sexes (⊕, ♀, ♂) were collected. Variance associated to geographic scales, sequencing depth (blue) and residuals (grey) were extracted from univariate mixed models and the respective 95% confidence intervals were computed by bootstrapping with a 1,000 iterations. (a-c) Prokaryotic diversity parameters. (d-f) Eukaryotic diversity parameters. (g-j) Abundance of predicted functional groups (variation explained by the LSD is not shown).

Accepted Article

  terrestrial plants hosting diazotrophs are highly successful invaders (e.g., Myrica faya in Hawai; Vitousek, Walker, Whiteaker, Mueller-Dombois & Matson 1987). Hypothetically, an intimate relationship with Cyanobacteria capable of fixing dinitrogen -coastal marine waters are commonly nitrogen limited (Vitousek & Howarth 1991) -could provide A. vermiculophyllum a substantial advantage to survive in nonnative habitats.

FiguresFigure

  Figures

  Figure 1. Geographic scales included in the sampling of this study: ranges (native and nonnative), coasts, ecoregions and populations. Sampled ecoregions are in blue (native) and orange (nonnative) and labelled with

Figure 2 .

 2 Figure 2. Differences among micro-environments. (a-b) Non-metric dimensional scaling (nMDS) plots and effects for water, epi-and endophytic microbial communities, adjusted for the effects of sequencing depth and population identity. Upper insets on the right display nMDS plots based on the raw data and lower right insets on the mGLM residuals also adjusted for effects of micro-environment (ME). (a) Prokaryotic communities. (b) Eukaryotic communities. Effects from univariate models are displayed with 95% confidence intervals for prokaryotic S chao (c), prokaryotic logit PIE (d), eukaryotic asymptotic S chao (e), eukaryotic logit PIE (f) and predicted abundance of functional groups (g-j). Effects are labeled with letters to indicate significant differences (p < 0.05). The marginal and conditional R 2 values (R 2 m and R 2 c ) are displayed in the upper left or right corners.

Figure 3 .

 3 Figure 3. Differentially abundant OTUs from mGLMs (a) among micro-environments (tissue, surface and water), (b) within different geographic scales (range, coast, region and population) and (c) populations and years (population and year). Bars represent percentages of differentially abundant OTUs and numbers next to the bars to the corresponding OTU counts. Prokaryotic and eukaryotic results are shown in red and yellow, respectively. Grey bars in (a) indicate the amount of OTUs significantly associated with water, epiphytic, algal (epi and endo) and endophytic samples (composition cores, see main text for details). The results are based on a significance level of p < 0.05.

Figure 4 .

 4 Figure 4. Variation in microbial communities across scales, visualized as percentages of explained variance.

Acknowledgements

We thank Caitlin Cox and Mike Crowley at the Heflin Center for Genomic Sciences (University of Alabama at Birmingham) for use of the capillary sequencer for fragment analysis. We are also grateful to Nadja Stärck and Anna Neu for support during the fieldwork. This study was funded by grants from the Deutsche Forschungsgemeinschaft awarded to FW (DFG grant number WE2700/5-1) and TB (DFG grant number BA5508/2-1). Start-up funds from the University of Alabama at Birmingham to SAKH supported North American field work and fragment analysis. FRB acknowledges the financial support of the German Academic Exchange Service (DAAD) through the project Doctoral Programs in Germany 2015/16 (57129429).

Accepted Article

This article is protected by copyright. All rights reserved Aires T., Moalic Y, Serrão EA, Arnaud-Haond S (2015) Hologenome theory supported by cooccurrence networks of species-specific bacterial communities in siphonous algae (Caulerpa). FEMS microbiology ecology, 91, fiv067. https://doi.10.1093/femsec/fiv067. Arnaud-Haond S., Aires T, Candeias R, Teixeira S, Duarte CM, Valero M, ... Serrão E (2017) Entangled fates of holobiont genomes during invasion: nested bacterial and host diversities in Caulerpa taxifolia. Molecular ecology, 26, 2379-2391. Aßhauer K. P., Wemheuer B, Daniel R, Meinicke P (2015) Tax4Fun: predicting functional profiles from metagenomic 16S rRNA data. Bioinformatics, 31, 2882-2884. https://doi.10.1093/bioinformatics/btv287. Astudillo-García C., Bell JJ, Webster NS, Glasl B, Jompa J, Montoya JM, ... Taylor MW (2017) Evaluating the core microbiota in complex communities: a systematic investigation. Environmental microbiology, 19, 1450-1462. https://doi.10.1111/1462-2920.13647. Bates D., Mächler M, Bolker B, Walker S (2015) Fitting linear mixed-effects models using lme4. arXiv preprint arXiv:1406.5823, 67, 1-48. https://doi.10.18637/jss.v067.i01. Bell T., Newman JA, Silverman BW, Turner SL, Lilley AK (2005) The contribution of species richness and composition to bacterial services. Nature, 436, 1157. https://doi.10.1038/nature03891. Bellorin A. M., Oliveira MC, Oliveira EC (2004) Gracilaria vermiculophylla: a western Pacific species of Gracilariaceae (Rhodophyta) first recorded from the eastern Pacific. Phycological Research, 52, 69-79. https://doi.10.1111/j.1440-183.2004.00330.x. Bordenstein S. R., Theis KR (2015) Host biology in light of the microbiome: ten principles of holobionts and hologenomes. PLoS biology, 13, e1002226. https://doi.10.1371/journal.pbio.1002226. Burke C., Steinberg P, Rusch D, Kjelleberg S, Thomas T (2011) Bacterial community assembly based on functional genes rather than species. Proceedings of the National Academy of Sciences of the United States of America, 108, 14288-14293. https://doi.10.1073/pnas.1101591108. Byers J. E., Gribben PE, Yeager C, Sotka EE (2012) Impacts of an abundant introduced ecosystem engineer within mudflats of the southeastern US coast. Biological Invasions, 14, 2587-2600. https://doi.10.1007/s10530-012-0254-5.

Accepted Article

This article is protected by copyright. All rights reserved Campbell A. H., Marzinelli EM, Gelber J, Steinberg PD (2015) Spatial variability of microbial assemblages associated with a dominant habitat-forming seaweed. Frontiers in microbiology, 6, 230. https://doi.10.3389/fmicb.2015.00230. Chao A., Bunge J (2002) Estimating the number of species in a stochastic abundance model. Biometrics, 58, 531-539. https://doi.10.1111/j.0006-341X.2002.00531.x. Coyte K. Z., Schluter J, Foster KR (2015) The ecology of the microbiome: Networks, competition, and stability. Science, 350, 663-666. https://doi.10.1126/science.aad2602. Egan S., Harder T., Burke C., Steinberg P., Kjelleberg S., Thomas T. (2013) The seaweed holobiont: understanding seaweed-bacteria interactions. FEMS microbiology reviews, 37, 462-476. https://doi:10.1111/1574-6976.12011. Freshwater D. W., Montgomery F, Greene JK, Hamner RM, Williams M, Whitfield PE (2006) Distribution and identification of an invasive Gracilaria species that is hampering commercial fishing operations in southeastern North Carolina, USA. Biological Invasions, 8, 631-637. https://doi.10.1007/s10530-005-1809-5. Gerard V., Dunham S, Rosenberg G (1990) Nitrogen-fixation by cyanobacteria associated with Codium fragile (Chlorophyta): environmental effects and transfer of fixed nitrogen. Marine Biology, 105, 1-8. Gohl D. M., Vangay P, Garbe J, MacLean A, Hauge A, Becker A, ... Hunter R (2016) Systematic improvement of amplicon marker gene methods for increased accuracy in microbiome studies. Nature biotechnology, 34, 942. https://doi.10.1038/nbt.3601. Gonzalez D. J., Smyth AR, Piehler MF, McGlathery KJ (2013) Mats of the nonnative macroalga, Gracilaria vermiculophylla, alter net denitrification rates and nutrient fluxes on intertidal mudflats. Limnology and Oceanography, 58, 2101-2108. https://doi.10.4319/lo.2013.58.6.2101. Gurgel C. F. D., Norris JN, Schmidt WE, Le HN, Fredericq S (2018) Systematics of the Gracilariales (Rhodophyta) including new subfamilies, tribes, subgenera, and two new genera, Agarophyton gen. nov. and Crassa gen. nov. Phytotaxa, 374, 1-23. https://doi.10.11646/phytotaxa.374.1.1. Hammann M., Buchholz B, Karez R, Weinberger F (2013) Direct and indirect effects of Gracilaria vermiculophylla on native Fucus vesiculosus. Aquatic Invasions, 8, 121-132. https://doi.10.3391/ai.2013.8.2.01.

Accepted Article

This article is protected by copyright. All rights reserved Hammann M., Rempt M, Pohnert G, Wang G, Boo SM, Weinberger F (2016a) Increased potential for wound Accepted Article This article is protected by copyright. All rights reserved Suzuki R., Shimodaira H (2006) Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics, 22, 1540-1542. https://doi.10.1093/bioinformatics/btl117. Tedersoo L., Bahram M, Polme S, Koljalg U, Yorou NS, Wijesundera R, ... Abarenkov K (2014) Fungal biogeography. Global diversity and geography of soil fungi. Science, 346, 1256688. https://doi.10.1126/science.1256688. Thomsen M. S., Gurgel CFD, Fredericq S, McGlathery KJ (2006) Gracilaria vermiculophylla (Rhodophyta, Gracilariales) in Hog Island Bay, Virginia: A cryptic alien and invasive macroalga and taxonomic correction. Journal of Phycology, 42, 139-141. https://doi.10.1111/j.1529-8817.2006.00160.x. Vitousek P. M., Howarth RW (1991) Nitrogen limitation on land and in the sea: how can it occur? Biogeochemistry, 13, 87-115. Vitousek P. M., Walker LR, Whiteaker LD, Mueller-Dombois D, Matson PA (1987) Biological Invasion by Myrica faya alters ecosystem development in Hawaii. Science, 238, 802-804. https://doi.10.1126/science.238.4828.802.

Accepted Article

This article is protected by copyright. All rights reserved Warton D. I., Blanchet FG, O'Hara RB, Ovaskainen O, Taskinen S, Walker SC, ... Hui FK (2015) So many variables: joint modeling in community ecology. Trends in Ecology & Evolution, 30, 766-779. https://doi.10.1016/j.tree.2015.09.007. Weinberger F., Buchholz B, Karez R, Wahl M (2008) The invasive red alga Gracilaria vermiculophylla in the Baltic Sea: adaptation to brackish water may compensate for light limitation. Aquatic Biology, 3, 251-264. https://doi.10.3354/ab00083. Weinberger F., Friedlander M (2000) Response of Gracilaria conferta (Rhodophyta) to oligoagars results in defense against agar-degrading epiphytes. Journal of Phycology, 36, 1079-1086. https://doi.10.1046/j.1529-8817.2000.00003.x. Wittebolle L., Marzorati M, Clement L, Balloi A, Daffonchio D, Heylen K, ... Boon N (2009) Initial community evenness favours functionality under selective stress. Nature, 458, 623. https://doi.10.1038/nature07840. [dataset] Bonthond G., Bayer T., Krueger-Hadfield S. A., Barboza F. R., Nakaoka M., … Weinberger F.; 2019; PRJNA564581; Short Read Archive; https://www.ncbi.nlm.nih.gov/bioproject/PRJNA564581/ Data Accessibility De-multiplexed V4-16S and V7-18S rRNA gene amplicon reads and associated metadata are deposited in the SRA database under the Bioproject accession number PRJNA564581.

Author Contributions

GB, TB and FW conceptualized the study. Field collections were conducted by GB, TB, SAKH, MN, MV, GW and FW. GB, SAKH and SK conducted laboratory work. GB processed and analyzed the data. FRB helped analyze and interpret the data. GB drafted the manuscript. TB, SAKH, FRB, MN, MV, GW, SK and FW helped writing the manuscript.