Phytoplankton support complex bacterial microbiomes that rely on phytoplankton-derived extracellular compounds and perform functions necessary for algal growth. Recent work has revealed sophisticated interactions and exchanges of molecules between specific phytoplankton-bacteria pairs, but the role of host genotype in regulating those interactions is unknown. Here, we show how phytoplankton microbiomes are shaped by intraspecific genetic variation in the host using global environmental isolates of the model phytoplankton host Thalassiosira rotula and a laboratory common garden experiment. A set of 81 environmental T. rotula genotypes from three ocean basins and eight genetically distinct populations did not reveal a core microbiome. While no single bacterial phylotype was shared across all genotypes, we found strong genotypic influence of T. rotula, with microbiomes associating more strongly with host genetic population than with environmental factors. The microbiome association with host genetic population persisted across different ocean basins, suggesting that microbiomes may be associated with host populations for decades. To isolate the impact of host genotype on microbiomes, a common garden experiment using eight genotypes from three distinct host populations again found that host genotype influenced microbial community composition, suggesting that a process we describe as genotypic filtering, analogous to environmental filtering, shapes phytoplankton microbiomes. In both the environmental and laboratory studies, microbiome variation between genotypes suggests that other factors influenced microbiome composition but did not swamp the dominant signal of host genetic background. The long-term association of microbiomes with specific host genotypes reveals a possible mechanism explaining the evolution and maintenance of complex phytoplankton-bacteria chemical exchanges.

phytoplankton-bacteria interactions j microbiome j phycosphere j population genetics j host-microbe interactions I nteractions between marine phytoplankton and bacteria can exert a profound influence on ecosystem function and biogeochemical cycling, impacting rates of primary production, phytoplankton aggregation, organic carbon export, and nutrient cycling (1)(2)(3). As an aquatic analog of the plant rhizosphere, the most intimate relationships between phytoplankton and bacteria exist in the phycosphere, the region immediately surrounding a phytoplankton cell, where molecules can be exchanged despite the effects of turbulence and diffusion (4). Relationships between phytoplankton and bacteria in the phycosphere range from cooperative to competitive (5). For example, during exponential growth, phytoplankton actively secrete amino acids that are taken up by bacteria, despite potentially significant energy costs (6). In turn, some bacteria synthesize essential vitamins and growth hormones that stimulate phytoplankton productivity (7,8). During phytoplankton senescence, formerly "friendly" bacteria can become pathogenic, producing algicides that lyse phytoplankton cells and release organic carbon to the environment (9). These complex ecological interactions have been investigated in the laboratory for specific phytoplankton-bacteria pairs (3). However, the persistence and phylogenetic breadth of these relationships for both host and microbiome remain open questions (10)(11)(12).

In terrestrial habitats, clear linkages exist between the genetic background (i.e., host genotype and population genetic structure) of foundational plant species and the organisms that rely on them. For example, genetic variation within tree species can influence the structure of associated epiphytic, mycorrhizal, and invertebrate communities (13)(14)(15) through mutualism, parasitism, commensalism, facilitation, and competition (reviewed in ref. 16). These associations extend to plant-associated bacteria, whose abundance, composition, and diversity reflect intraspecific trait variation among host genotypes (17). Because microbes regulate processes such as decomposition, nutrient dynamics, and energy flow, the influence of intraspecific genetic variation in plants on their associated bacteria extends the effects of community genetics to ecosystem processes (18). In contrast to terrestrial habitats, seawater allows both bacteria and their phytoplankton hosts to drift with nearly unlimited dispersal across the global ocean. It is unknown whether planktonic communities mirror those in terrestrial habitats, where host genetics can shape bacterial community composition and influence ecosystem processes (16), or whether dispersal in the dilute marine environment overwhelms the formation of close and specific relationships between phytoplankton and their

Microscale interactions between marine phytoplankton and their bacterial microbiomes can influence ecosystem functioning and global biogeochemical cycling through complex exchanges of metabolites and sophisticated ecological processes. Previous investigation of the phytoplankton microbiome has not focused on the role of a host's underlying genetic background. Through examination of a single phytoplankton species' microbiome across the global ocean, we found that host genotype strongly influenced microbiome community composition, with associations that potentially persist across generations and ocean basins but assemble rapidly (within days). The long-term association of microbiomes with host genetic background could explain the evolution and maintenance of intricate phytoplankton-bacteria interactions. associated microbiota. Although phytoplankton drift freely across the global ocean with nearly unlimited dispersal and divide primarily asexually, they still possess clear genetic structure. Phytoplankton species are organized into genetically distinct populations that possess phenotypic trait variation, are associated with specific environmental conditions, and undergo large increases in abundance, known as blooms (19)(20)(21)(22)(23). Furthermore, genetically distinct phytoplankton populations persist on time scales of decades to centuries (23,24), providing ample opportunity for populations to develop specific relationships with other microbes. Few studies have evaluated the microbiomes of multiple strains within a phytoplankton species. Some studies found that individual phytoplankton species may possess a core microbiome, with a consistent set of bacterial phylotypes and metabolic potentials (25,26), while others found that phytoplankton strains supported distinct microbiomes (10) and differed in their growth responses to the same bacterial strain (7,27,28). These intriguing findings have not been rigorously examined in light of the genetic background of the host phytoplankton species. Given that plant genotype often (29)(30)(31) but not always (32) drives host microbiomes and given the differences between terrestrial and planktonic habitats, understanding to what extent the genetic background of phytoplankton species influences microbiome composition is critical to understanding the nature of their interactions and parsing the roles of both partners in global biogeochemical cycles.

Here, we assessed host-microbe interactions using the model marine phytoplankton Thalassiosira rotula, a cosmopolitan species characterized by high genotypic and phenotypic diversity, which is subdivided into genetically distinct populations (24,33). We examined whether the T. rotula microbiome was influenced by host genetic background, either at the genotype or population level. We combined a study of the microbiomes of 81 environmental genotypes (representing eight genetically distinct populations) sampled from around the world (Fig. 1 and Table 1) with a common garden experiment explicitly designed to evaluate the influence of host genetic background on microbiome assembly in the absence of other environmental changes. The union of these approaches allowed us to examine the role of host genetics in setting bacterial assemblage composition and to determine whether a core microbiome might exist in a globally distributed, pelagic phytoplankton.

Characterizing the T. rotula Microbiome. Marine phytoplankton and bacteria have well-established biogeochemical linkages through bulk pools of nutrients and organic carbon, although the specificity of exchanges and phylogenetic breadth of the partnerships between these organisms remain areas of active investigation (3). Here, we use a rich dataset of 81 single-cell environmental isolates of the phytoplankton species T. rotula and their co-isolated bacterial microbiomes paired with controlled laboratory experiments to examine the relationship between phytoplankton host genetic background and bacterial microbiome. Because environmental T. rotula microbiomes were obtained from washed single cells or chains, the microbiomes we recovered were largely from bacterial cells physically associated with the phytoplankton host. In order to increase the likelihood that bacterial sequences examined were derived from the in situ algal microbiome, we instituted a number of controls, which likely eliminated some real microbial diversity but reduced false relationships and removed potential contaminants. These included washing phytoplankton cells prior to culturing to reduce free-living bacteria to an estimated <1 cell per freshly-isolated isolated phytoplankton, using a low-organic matter medium to support algal-derived bacterial growth and employing strict quality controls for sequence analyses (SI Appendix, SI Results and Discussion).

Each environmental T. rotula isolate represented a unique multilocus genotype (Dataset S1) and possessed a diverse microbiome; after strict quality controls, we identified 278 bacterial amplicon sequence variants (ASVs) that clustered at the 99.9% identity level, with a range of eight to 51 bacterial ASVs and an average of 24.4 ± 9.6 ASVs per T. rotula environmental genotype (Dataset S2). Our results are consistent with previous reports of bacterial richness ranging from 1 to 125 operational taxonomic units (OTUs; roughly equivalent to ASVs) per phytoplankton isolate (11,12,34), including eight to 80 OTUs per T. rotula isolate (10,35). Furthermore, the T. rotula-associated bacterial communities were significantly different from the in situ bacterioplankton communities (SI Appendix, Fig. S1 and Dataset S2; Analysis of Similarities (ANOSIM) R = 0.80, P = 0.001), highlighting that microbiome composition does not simply reflect water column bacterioplankton composition (3).

The T. rotula-associated bacterial communities from the environmental genotypes were dominated by Proteobacteria (83.5% of all sequences), including γ-proteobacteria (61.3%) and α-proteobacteria (21.3%) (Dataset S2). However, no single ASV was present in all microbiomes. The most frequent ASVs in the T. rotula microbiomes were two α-proteobacteria Loktanella (ASVs 743 and 1,938) and Actinobacteria Rhodococcus (ASV 24), present in 59%, 44%, and 56% of all environmental genotypes, respectively (SI Appendix, Fig. S2). The relative abundance of individual ASVs per host genotype ranged widely, from 0 to >90% in a microbiome (SI Appendix, Fig. S3). For example, when present, the ASV most likely to occur in a microbiome, Loktanella 743, had a relative median abundance of 4.41% but a range of 0 to >80% in any given microbiome. The absence of a common ASV across microbiomes and high degree of variability in microbiome composition indicates that the cosmopolitan phytoplankton T. rotula sampled from across the global ocean lacks a "core" microbiome. Instead, T. rotula cells may associate with closely related lineages of bacteria that serve similar ecological roles but are adapted to distinct environmental conditions. With few exceptions (25), the lack of a core microbiome stands in contrast to most surveys of microbiome communities from multicellular eukaryotes (36) and marine algae (12).

Although no single bacterial ASV was present in all microbiomes, similarities existed at the taxonomic order level, with (24), and black indicates populations identified in this study] and their associated microbiomes (Table 1 and SI Appendix, Table S1). Three sites (Wa, Nb, and Fr) were resampled (symbols with two colors), and whole seawater was collected twice from Narragansett Bay (Nb) for microbiome comparisons with the in situ whole seawater bacterial community (asterisks). The base map is a composite of log annual average chlorophyll a concentrations (milligrams/meter À3 , 2010) (https://oceandata.sci.gsfc.nasa.gov/directaccess/ MODIS-Aqua/Mapped/Annual/9km/chlor_a/). the order Alteromonadales observed in 80 of the 81 environmental isolates. The Alteromonadales are known to associate with a number of phytoplankton lineages (12,37). Here, we supplied all required substrates for algal growth in the media, including vitamins and iron; thus, T. rotula may have benefited from the bacterial reduction of both oxidative stress and waste compounds (5). Alternately, under these replete conditions, T. rotula may have derived no benefit from the bacteria, with the relationship being driven by bacterial metabolism of phytoplankton-derived organic material. Overall, these data suggest a strong relationship of the Alteromonadales with T. rotula.

The Role of Host Genetic Background in Structuring T. rotula Microbiomes. In order to better understand the factors that shape T. rotula-associated microbiomes, we examined the roles of stochasticity, environmental selection, and host genetic background. To test how stochastic effects, such as dispersal, ecological drift, and random death and birth events (38) influenced the microbiomes of environmental isolates, we used the normalized stochastic ratio of the Jaccard matrix (NST jac ), which quantifies the data structure ranging from deterministic (0%) to completely stochastic (100%). The NST jac averaged 36.2 ± 18.4% across all diatom microbiomes, indicating while some stochasticity was evident in the data structure, deterministic effects strongly impacted the microbiomes of environmental isolates.

To determine whether host genotype or environmental factors explained the composition of host microbiomes, we first examined microbiome composition in light of the T. rotula genetic background. Even at the bacterial class level, there appeared to be stark differences in the microbiomes of different host genotypes (SI Appendix, Fig. S4). Variation among microbiomes was not evenly spread across the global dataset. Instead, we found that phytoplankton isolates from the same water sample consistently coisolated with a more similar bacterial community than isolates from different water samples (ANOSIM R statistic = 0.83, P = 0.001), with strong separation between cocultured and in situ bacterial communities (SI Appendix, Fig. S1). Of the microbiomes we analyzed, the host was subdivided into eight genetically distinct populations from three ocean basins (Fig. 2A and SI Appendix, Tables S1 andS2). By grouping T. rotula genotypes into populations, we observed the same T. rotula population in different locations (e.g., North Sea [NS] and Olympic Peninsula B [WaB]) (24) as well as single locations from which different populations were sampled (e.g., Narragansett Bay Plankton Time Series A and K [NbA and Nbk]). We then performed multiple tests to determine how host population and environmental factors shaped the composition of the T. rotula microbiome.

First, microbiome composition was most strongly and equally correlated with host population membership and sea surface temperature (SST) (bioenv correlations both 0.51). Interestingly, the genetic composition of the host population was also strongly influenced by temperature (24), suggesting possible collinearity in the dataset. To further examine the relationship between host microbiome, host population, and SST, variance partitioning was conducted using two tests with different assumptions. We found that regardless of the variance-partitioning test used, host population explained by far the largest portion of variance in microbiome composition (adonis R 2 = 0.47 P = 0.001; varpart R 2 = 0.35), followed distantly by SST (adonis R 2 = 0.057 P = 0.001; varpart R 2 ¼0.04), and the interaction between SST and phytoplankton population (adonis R 2 = 0.038 P = 0.001; varpart R 2 = 0.01). We found no correlation between microbiome composition and chlorophyll a concentration, salinity, and T. rotula concentration. We cannot rule out other unmeasured factors (e.g., nutrients) that could affect both the microbiome and algal population structure (20,26). Although temperature is known to shape microbiome composition and physiology (39), genetic background of the algal host has not previously been shown to play a role in microbiome community composition.

We observed a significant relationship between host population and the taxonomic composition of its microbiome (Fig. 2B, ANOSIM R Statistic = 0.66, P = 0.001). Environmental genotypes from the same host population had microbiomes with compositional similarities even when those genotypes originated from different ocean basins (e.g., populations 1 and 4, Fig. 2B). For example, the microbiomes of T. rotula population 1 sampled from the north Pacific and north Atlantic (SOMLIT-Astan Time Series, France [FrA] and Olympic Peninsula A [WaA]) shared 20 ASVs (43 and 30% of the ASVs in FrA and WaA microbiomes, respectively). Two ASVs, α-proteobacteria Loktanella (ASV 743) and Actinobacteria Rhodococcus (ASV 24), were present in all population 1 genotypes from both ocean basins. The microbiome similarities among environmental genotypes from different ocean basins are notable, because rates of surface water circulation in the global ocean can translate into decades-long separation of individuals from the same host population (40).

Overall, compositional divergence between microbiomes was significantly correlated with genetic divergence between host phytoplankton populations (F ST ) (Fig. 3A) (mantel r = 0.45, P = 0.01; linear regression R 2 = 0.17, F = 6.63, P = 0.02). In contrast, geographic distance had no significant impact on the phytoplankton microbiome (Fig. 3B) (mantel r = 0.18, P = 0.12; linear regression R 2 = 0.02, F = 1.82, P = 0.18). Positive Additional sample metadata are in SI Appendix, Table S1. * Whole seawater was collected from Narragansett Bay on two occasions to determine bacterioplankton community composition. †

Host populations genotyped in this study.

geographic distance-decay relationships are signatures of neutral processes (41), evidenced by the biogeographic patterns of the global bacterioplankton community (42), and have been shown to occur in some plankton microbiomes (43). In contrast, our findings suggest that phytoplankton microbiomes sampled at different times or places (10) may reflect the genetic background of the host phytoplankton and are not simply the result of colonization by the bacteria responding most rapidly to phytoplankton exudates, as models have predicted (44).

To explicitly test whether distinct bacterial communities could assemble in response to T. rotula with different genetic backgrounds in the absence of potentially confounding environmental differences, we performed a common garden experiment under controlled conditions using host cells freshly isolated from the field (Table 2 and Dataset S3) (35,45,46). Eight host genotypes were made axenic, inoculated with a <1-μm natural seawater bacterial community, incubated for 5 d, and then the entire contents of the flasks were harvested by filtering. The final concentrations of host T. rotula were high (22,575 ± 10,129 cells/mL À1 ), similar to levels observed during phytoplankton blooms (47). Subsequent sequencing of treatment flasks revealed a total of 419 ASVs, with an average of 330. By evaluating replicates for each genotype as well as multiple genotypes within T. rotula populations, the common garden experiment allowed us to examine the relative importance of stochastic effects as well as the influence of two levels of host genetic background on bacterial composition: the underlying genotype of the host and host membership in a particular genetic population. Deterministic processes appeared to play a larger role than stochasticity in shaping host microbiomes in the common garden experiment (NST jac 23.0% [PopA], 39.0% [PopB], and 37.8% [PopC]). Consistent with a key role for selection, we found that replicate flasks from the same T. rotula common garden genotype yielded highly reproducible bacterial communities (Fig. 4B; ANOSIM, R = 0.93, significance = 0.001; adonis R 2 = 0.76, P = 0.001). Next, we examined the role of population. The T. rotula used in the common garden experiment originated from three genetically distinct populations (SI Appendix, Tables S3 andS4). We found that genotypes from the same population had more similar bacterial communities than those from other populations (Fig. 4B and Dataset S4; ANO-SIM, R = 0.52, significance = 0.001; adonis R 2 = 0.38, P = 0.001). Remarkably, the population-specific effects were evident rapidly, within the 5-d incubation period, suggesting that even stronger differences might have emerged over longer time frames but with concurrent increases in "bottle effects" (48,

Fig. 4. Bacterial community composition associated with T. rotula genotypes in a common garden experiment. (A) Percentage of relative abundance of bacterial taxonomic orders, obtained using 16S rRNA gene sequencing, after a 5-d common garden experiment. Eight T. rotula genotypes were grown in triplicate (lower x-axis labels), selected from three genetically distinct populations (upper x-axis), and inoculated (Inoc) with a <1 μm natural seawater bacterial community (Bact_T0) taken from a whole seawater sample (WSW_T0). Bacterial controls (Bacteria, upper x-axis) were analyzed on day 5 and were run in quadruplicate. (B) NMDS of a normalized Bray-Curtis dissimilarity matrix of 99.9% identity 16S rRNA gene ASV composition of bacterial communities from different T. rotula genotypes and populations (shapes, colors) and from bacterial controls (Bacteria) at the time of harvest. Ellipses represent 80% CI for the microbiomes from each genetically distinct phytoplankton population (stress = 0.08).

ECOLOGY 49). While the inoculum communities were dominated by SAR11, the microbiome cocultured with T. rotula was enriched in fastgrowing organisms (e.g., Vibrio, Marinomonas) similar to those in Fu (46). Marine bacterioplankton are known to respond rapidly to primary producer-derived microenvironments including nutrient plumes (50), lysed phytoplankton (44), and both real and synthetic phytoplankton (35,51). Here, we showed that marine bacteria rapidly differentiate their composition in the presence of phytoplankton with different genetic backgrounds. Despite similarities among the microbiomes of phytoplankton isolates from the same genetic population, T. rotula-associated bacterial communities exhibited variability both among genotypes from a single population and among replicates of a single genotype, suggesting that other mechanisms, such as stochastic and priority effects (52)(53)(54), also likely played a role in community assembly. As we showed, stochastic effects, while apparent, were not the dominant drivers shaping T. rotula microbiome composition. Priority effects could influence the microbiome through historical contingencies shaping the pool of available microbes as well as through colonization of phytoplankton surfaces (55). The common garden experiment could have been especially susceptible to priority effects, since we examined community assembly using axenic genotypes, a scenario in which abundant and fastgrowing phylotypes may dominate over phycosphere-adapted taxa.

Microbiome variation among genotypes within the same host population may well be explained by a phenomenon that we describe here as genotypic filtering, similar to the process of environmental filtering (56) or host filtering (57). Genotypic variation is associated with high levels of phenotypic variation in phytoplankton, including T. rotula (33,58). Prior observations of extensive phenotypic variation of T. rotula together with recent evidence that exometabolites from different phytoplankton species predictably influence microbiome composition (51), suggest that different T. rotula genotypes may release distinct organic compounds that select for or attract different bacterial associates. Phytoplankton like T. rotula have been shown to promote the growth and attachment of beneficial bacteria and inhibit colonization of opportunistic bacteria through secondary metabolites (59). Phenotypic variation among genotypes may thus lead to genotypic filtering, a mechanism not previously attributed to phytoplankton-associated bacterial communities. The taxonomic extent of genotypic filtering in the phytoplankton is unknown. Intriguingly, there is evidence of phytoplanktonintraspecific variation in exudate profiles, which are well-known to influence bacterial composition and activities (60) and evidence of genotype-specific responses to bacteria (7), both of which suggest that genotypic filtering may be a widespread phenomenon in the phytoplankton.

Microbiome stability among genotypes separated by globalscale geographic distances may be maintained in part by the phytoplankton life cycle. Prior work has posited that phytoplankton microbiomes are short lived with a duration of just 1 to 2 d, equivalent to the lifespan of phytoplankton cells (51). While it is possible that microbiomes reassemble on daily timescales, our results suggest that additional mechanisms may act to promote the long-term similarity of the phytoplankton microbiome, such as vertical transfer of bacteria during reproduction (61). The host T. rotula divides primarily asexually, yielding two daughter cells that each receive from the mother cell one-half of a hard, silicified cell-wall covering known as a frustule (62), which presumably contains a near-complete microbiome able to colonize the newly-formed half frustule. Sexual reproduction in phytoplankton like T. rotula occurs on time scales of once per year to once every 40 y, although the details are not well known (62). As sexually generated offspring cast off both halves of the parent frustule and a new frustule is formed, microbiome colonization would occur from bacteria in the surrounding water column. It is interesting to note that in colony-forming phytoplankton like T. rotula, it has been observed that only a few cells in the monoclonal colony tend to undergo sexual reproduction (62), providing a route for microbiome colonization from neighboring cells and long-term stability of microbiome community composition.

The process of T. rotula microbiome community assembly can be explained through an ecological framework of selection, drift, dispersal, and diversification (53,63). For selection and drift, the NST jac (38) of both the environmental isolates (36.2 ± 18.4%) and common-garden experiment (33.3 ± 8.9%) indicates that the genotype selection effect is dominated by deterministic rather than stochastic factors, which provides a foundation for exploring assembly mechanisms in future studies. With regards to dispersal, we find that the same genetic populations of T. rotula sampled from different locations carry common ASVs, suggesting reassociation of microbiomes and maintenance of certain microbiome partners despite geographic separation. Finally, we found evidence of genotypic filtering in T. rotula microbiomes, which could lead to the diversification and maintenance of different populations of T. rotula.

Here, across both global environmental samples and a common-garden experiment, we observed that host genotype shapes the microbiomes of a single-celled phytoplankton. The association of bacteria with host genetic background observed here has the potential to impact host ecology and biogeochemical cycling. For example, microbiomes are known to alter growth and survival rates of the host (64,65) and expand host metabolic potential (26). The evolutionary fingerprint of the phytoplankton host on bacterial microbiomes opens routes of inquiry into factors that initiate and maintain microbiome composition and function. In particular, the long-term association (or frequent reassociation) of bacteria to a geneticallyrestricted subset of a species sets up the potential for significant eco-evolutionary dynamics to occur in both partners and provides a mechanism for the evolution and maintenance of recently identified sophisticated molecular exchanges between phytoplankton and bacteria (3).

Global Environmental Sampling. Surface water containing T. rotula cells was collected 12 times from nine locations in the Atlantic, Pacific, and Indian Oceans (Fig. 1, Table 1, and SI Appendix, Table S1) (24) and single cells or chains of T. rotula were isolated and genotyped (SI Appendix, SI Materials and Methods). We amplified the V3-V4 region of the 16S ribosomal ribonucleic acid (rRNA) gene (SI Appendix, SI Materials and Methods) and paired-end sequencing (2 × 250 bp) of the amplicons was performed on a MiSeq (Illumina) using version 2 chemistry at the Duke Genome Sequencing and Analysis Core Facility. Quality control, ASV clustering at 99.9% identity, and statistical analysis are described in SI Appendix, SI Materials and Methods. 2), genotyped using microsatellite markers, and maintained in axenic culture (SI Appendix, SI Materials and Methods). Direct comparisons were not made between the common garden and global environmental T. rotula genotypes, because we modified PCR and fragment analysis protocols compared to Whittaker and Rynearson (24) and used different 16S rRNA gene primers between the two datasets (SI Appendix, SI Materials and Methods). Eight axenic genotypes were chosen to undergo the common garden experiment described in SI Appendix, SI Materials and Methods. After the 5-d experiment, we amplified the V4-V5 region of the 16S rRNA gene of the triplicate common garden genotypes, quadruplicate bacterial controls, and starting inoculum, sequenced the amplicons on a MiSeq. 2 × 250 bp run (Illumina), and performed bioinformatic and statistical analysis (SI Appendix, SI Materials and Methods).

Data Availability. All 16S rRNA gene sequences from global environmental genotype microbiomes are available in the National Center for Biotechnology Information (NCBI) Sequence Read Archive (BioProject PRJNA494816) under accession nos. SRR8046252 through SRR8046343. All 16S rRNA gene sequences from the common garden genotype microbiomes are available in the NCBI Sequence Read Archive (BioProject PRJNA705708) under accession nos. SRX10227248 through SRX10227277. All environmental data are deposited in the Biological and Chemical Oceanography Data Management Office Database for the environmental isolates (66) (https://www.bco-dmo.org/dataset/ 860347) and the common garden experiment (67) (https://www.bco-dmo.org/ dataset/860381). All remaining study data are included in the article and/or supporting information.