Host genotype structures the microbiome of a globally dispersed marine phytoplankton

A mechanistic framework for complex microbe-host symbioses

Virtually all multicellular organisms on Earth live in symbiotic associations with complex microbial communities: the microbiome. This ancient relationship is of fundamental importance for both the host and the microbiome. Recently, the analyses of numerous microbiomes have revealed an incredible diversity and complexity of symbionts, with different mechanisms identified as potential drivers of this diversity. However, the interplay of ecological and evolutionary forces generating these complex associations is still poorly understood. Here we explore and summarise the suite of ecological and evolutionary mechanisms identified as relevant to different aspects of microbiome complexity and diversity. We argue that microbiome assembly is a dynamic product of ecology and evolution at various spatio-temporal scales. We propose a theoretical framework to classify mechanisms and build mechanistic host-microbiome models to link them to empirical patterns. We develop a cohesive foundation for the theoretical understanding of the combined effects of ecology and evolution on the assembly of complex symbioses.

Impact of host species on assembly, composition, and functional profiles of phycosphere microbiomes

Microalgal microbiomes play vital roles in the growth and health of their host, however, their composition and functions remain only partially characterized, especially across microalgal phyla. In this study, a natural seawater microbiome was introduced to three distinct, axenic species of microalgae, the haptophyte Isochrysis galbana, the chlorophyte Tetraselmis suecica, and the diatom Conticribra weissflogii (previously Thalassiosira), and its divergence and assembly under constant illumination was monitored over 49 days using 16S rRNA amplicon and metagenomic analyses. The microbiomes had a high degree of host specificity in terms of taxonomic composition and potential functions, including CAZymes profiles. Rhodobacteraceae and Flavobacteriaceae families were abundant across all microalgal hosts, but I. galbana microbiomes diverged further from T. suecica and C. weissflogii microbiomes. I. galbana microbiomes had a much higher relative abundance of Flavobacteriaceae, whereas the two other algal microbiomes had higher relative abundances of Rhodobacteraceae. This could be due to the bacterivorous mixotrophic nature of I. galbana affecting the carbohydrate composition available to the microbiomes, which was supported by the CAZymes profile of I. galbana microbiomes diverging further from those of T. suecica and C. weissflogii microbiomes. Finally, the presence of denitrification and other anaerobic pathways was found exclusively in the microbiomes of C. weissflogii, which we speculate could be a result of anoxic microenvironments forming in aggregates formed by this diatom during the experiment. These results underline the significant role of the microalgal host species on microbiome composition and functional profiles along with other factors, such as the trophic mode of the microalgal host.

IMPORTANCE

As the main primary producers of the oceans, microalgae serve as cornerstones of the ecosystems they are part of. Additionally, they are increasingly used for biotechnological purposes such as the production of nutraceuticals, pigments, and antioxidants. Since the bacterial microbiomes of microalgae can affect their hosts in beneficial and detrimental ways, understanding these microbiomes is crucial to both the ecological and applied roles of microalgae. The present study advances the understanding of microalgal microbiome assembly, composition, and functionality across microalgal phyla, which may inform the modeling and engineering of microalgal microbiomes for biotechnological purposes.

Rhodobacteraceae are key players in microbiome assembly of the diatom Asterionellopsis glacialis

The complex interactions between bacterioplankton and phytoplankton have prompted numerous studies that investigate phytoplankton microbiomes with the aim of characterizing beneficial or opportunistic taxa and elucidating core bacterial members. Oftentimes, this knowledge is garnered through 16S rRNA gene profiling of microbiomes from phytoplankton isolated across spatial and temporal scales, yet these studies do not offer insight into microbiome assembly and structuring. In this study, we aimed to identify taxa central to structuring and establishing the microbiome of the ubiquitous diatom Asterionellopsis glacialis. We introduced a diverse environmental bacterial community to A. glacialis in nutrient-rich or nutrient-poor media in a continuous dilution culture setup and profiled the bacterial community over 7 days. 16S rRNA amplicon sequencing showed that cyanobacteria (Coleofasciculaceae) and Rhodobacteraceae dominate the microbiome early on and maintain a persistent association throughout the experiment. Differential abundance, co-abundance networks, and differential association analyses revealed that specific members of the family Rhodobacteraceae, particularly Sulfitobacter amplicon sequence variants, become integral members in microbiome assembly. In the presence of the diatom, Sulfitobacter species and other Rhodobacteraceae developed positive associations with taxa that are typically in high abundance in marine ecosystems (Pelagibacter and Synechococcus), leading to restructuring of the microbiome compared to diatom-free controls. These positive associations developed predominantly under oligotrophic conditions, highlighting the importance of investigating phytoplankton microbiomes in as close to natural conditions as possible to avoid biases that develop under routine laboratory conditions. These findings offer further insight into phytoplankton-bacteria interactions and illustrate the importance of Rhodobacteraceae, not merely as phytoplankton symbionts but as key taxa involved in microbiome assembly.

IMPORTANCE

Most, if not all, microeukaryotic organisms harbor an associated microbial community, termed the microbiome. The microscale interactions that occur between these partners have global-scale consequences, influencing marine primary productivity, carbon cycling, and harmful algal blooms to name but a few. Over the last decade, there has been a growing interest in the study of phytoplankton microbiomes, particularly within the context of bloom dynamics. However, long-standing questions remain regarding the process of phytoplankton microbiome assembly. The significance of our research is to tease apart the mechanism of microbiome assembly with a particular focus on identifying bacterial taxa, which may not merely be symbionts but architects of the phytoplankton microbiome. Our results strengthen the understanding of the ecological mechanisms that underpin phytoplankton-bacteria interactions in order to accurately predict marine ecosystem responses to environmental perturbations.

Mitigating salinity stress through interactions between microalgae and different forms (free-living & alginate gel-encapsulated) of bacteria isolated from estuarine environments

Salinity stress in estuarine environments poses a significant challenge for microalgal survival and proliferation. The interaction between microalgae and bacteria shows promise in alleviating the detrimental impacts of salinity stress on microalgae. Our study investigates this interaction by co-cultivating Chlorella sorokiniana, a freshwater microalga, with a marine growth-promoting bacterium Pseudomonas gessardii, both of which were isolated from estuary. In this study, bacteria were encapsulated using sodium alginate microspheres to establish an isolated co-culture system, preventing direct exposure between microalgae and bacteria. We evaluated microalgal responses to different salinities (5 PSU, 15 PSU) and interaction modes (free-living, gel-encapsulated), focusing on growth, photosynthesis, cellular metabolism, and extracellular polymeric substances (EPS) properties. High salinity inhibited microalgal proliferation, while gel-fixed interaction boosted Chlorella growth rate by 50.7 %. Both attached and free-living bacteria restored Chlorella's NPQ to normal levels under salt stress. Microalgae in the free-living interaction group exhibited a significantly lower respiratory rate compared to the pure algae group (-17.2 %). Increased salinity led to enhanced EPS polysaccharide secretion by microalgae, particularly in interaction groups (19.7 %). Both salt stress and interaction increased the proportion of aromatic proteins in microalgae's EPS, enhancing its stability by modulating EPS glycosidic bond C-O-C and protein vibrations. This alteration caused microalgal cells to aggregate, free-living bacteria co-culture group, and fixed co-culture group increasing by 427.5 %, 567.1 %, and 704.1 %, respectively. In gel-fixed bacteria groups, reduced neutral lipids don't accumulate starch instead, carbon redirects to cellular growth, aiding salt stress mitigation. These synergistic activities between salinity and bacterial interactions are vital in mitigating salinity stress, improving the resilience and growth of microalgae in saline conditions. Our research sheds light on the mechanisms of microalgal-bacterial interactions in coping with salt stress, offering insights into the response of estuarine microorganisms to global environmental changes and their ecological stability.

Antibiotics alter development and gene expression in the model cnidarian Nematostella vectensis

Background

Antibiotics are commonly used for controlling microbial growth in diseased organisms. However, antibiotic treatments during early developmental stages can have negative impacts on development and physiology that could offset the positive effects of reducing or eliminating pathogens. Similarly, antibiotics can shift the microbial community due to differential effectiveness on resistant and susceptible bacteria. Though antibiotic application does not typically result in mortality of marine invertebrates, little is known about the developmental and transcriptional effects. These sublethal effects could reduce the fitness of the host organism and lead to negative changes after removal of the antibiotics. Here, we quantify the impact of antibiotic treatment on development, gene expression, and the culturable bacterial community of a model cnidarian, Nematostella vectensis.

Methods

Ampicillin, streptomycin, rifampicin, and neomycin were compared individually at two concentrations, 50 and 200 µg mL-1, and in combination at 50 µg mL-1 each, to assess their impact on N. vectensis. First, we determined the impact antibiotics have on larval development. Next Amplicon 16S rDNA gene sequencing was used to compare the culturable bacteria that persist after antibiotic treatment to determine how these treatments may differentially select against the native microbiome. Lastly, we determined how acute (3-day) and chronic (8-day) antibiotic treatments impact gene expression of adult anemones.

Results

Under most exposures, the time of larval settlement extended as the concentration of antibiotics increased and had the longest delay of 3 days in the combination treatment. Culturable bacteria persisted through a majority of exposures where we identified 359 amplicon sequence variants (ASVs). The largest proportion of bacteria belonged to Gammaproteobacteria, and the most common ASVs were identified as Microbacterium and Vibrio. The acute antibiotic exposure resulted in differential expression of genes related to epigenetic mechanisms and neural processes, while constant application resulted in upregulation of chaperones and downregulation of mitochondrial genes when compared to controls. Gene Ontology analyses identified overall depletion of terms related to development and metabolism in both antibiotic treatments.

Discussion

Antibiotics resulted in a significant increase to settlement time of N. vectensis larvae. Culturable bacterial species after antibiotic treatments were taxonomically diverse. Additionally, the transcriptional effects of antibiotics, and after their removal result in significant differences in gene expression that may impact the physiology of the anemone, which may include removal of bacterial signaling on anemone gene expression. Our research suggests that impacts of antibiotics beyond the reduction of bacteria may be important to consider when they are applied to aquatic invertebrates including reef building corals.

Marine particle microbiomes during a spring diatom bloom contain active sulfate-reducing bacteria

Phytoplankton blooms fuel marine food webs with labile dissolved carbon and also lead to the formation of particulate organic matter composed of living and dead algal cells. These particles contribute to carbon sequestration and are sites of intense algal-bacterial interactions, providing diverse niches for microbes to thrive. We analyzed 16S and 18S ribosomal RNA gene amplicon sequences obtained from 51 time points and metaproteomes from 3 time points during a spring phytoplankton bloom in a shallow location (6-10 m depth) in the North Sea. Particulate fractions larger than 10 µm diameter were collected at near daily intervals between early March and late May in 2018. Network analysis identified two major modules representing bacteria co-occurring with diatoms and with dinoflagellates, respectively. The diatom network module included known sulfate-reducing Desulfobacterota as well as potentially sulfur-oxidizing Ectothiorhodospiraceae. Metaproteome analyses confirmed presence of key enzymes involved in dissimilatory sulfate reduction, a process known to occur in sinking particles at greater depths and in sediments. Our results indicate the presence of sufficiently anoxic niches in the particle fraction of an active phytoplankton bloom to sustain sulfate reduction, and an important role of benthic-pelagic coupling for microbiomes in shallow environments. Our findings may have implications for the understanding of algal-bacterial interactions and carbon export during blooms in shallow-water coastal areas.

Light-dark fluctuated metabolic features of diazotrophic and non-diazotrophic cyanobacteria and their coexisting bacteria

Cyanobacteria, the most abundant photosynthetic organisms in oceans, are tightly associated with diverse microbiota. However, the relationships between heterotrophic bacteria and cyanobacteria, particularly the diazotrophic group, are not fully understood. Here, we compared diel gene expressions of N2 fixing cyanobacteria Crocosphaera watsonii WH0003 and non-diazotrophic Synechococcus sp. RS9902 and their associated bacteria using metatranscriptomics approach. WH0003 showed significant up-regulation of O2 restriction and oxidative phosphorylation related genes at nighttime due to large carbon and energy investments for active N2 fixation. In contrast, RS9902 had higher expression for those genes at daytime. The two cyanobacteria hosted distinct bacterial communities with clear separate substrate utilization niches to reduce competition. Light-dark partitioning of nutrient acquisition among the dominant bacterial groups likely contributed to the dynamic balance for community coexistence. Moreover, particle-attached (PA) bacteria in RS9902 largely expressed glycoside hydrolases to hydrolyze complex carbohydrate compounds, while free-living (FL) bacteria priorly assimilated soluble, diffusible molecules. Spatial partitioning of nutrient acquisition between PA and FL bacteria implied that location initially influenced metabolic features of host associated bacteria. Our results advance knowledge on light-dark regulated metabolic activities of diazotrophic and non-diazotrophic cyanobacteria, and provide new insights into the coexisting strategies of different bacterial groups.

The diatom Thalassiosira rotula induces distinct growth responses and colonization patterns of Roseobacteraceae, Flavobacteria and Gammaproteobacteria

Diatoms as important phytoplankton components interact with and are colonized by heterotrophic bacteria. This colonization has been studied extensively in the past but a distinction between the bacterial colonization directly on diatom cells or on the aggregated organic material, exopolymeric substances (EPS), was little addressed. Here we show that the diatom Thalassiosira rotula and EPS were differently colonized by strains of Roseobacteraceae and Flavobacteriaceae in two and tree partner treatments and an enriched natural bacterial community as inoculum. In two partner treatments, the algae and EPS were generally less colonized than in the three partner treatments. Two strains benefitted greatly from the presence of another partner as the proportions of their subpopulations colonizing the diatom cell and the EPS were much enhanced relative to their two partner treatments. Highest proportions of bacteria colonizing the diatom and EPS occurred in the treatment inoculated with the enriched natural bacterial community. Dissolved organic carbon, amino acids and carbohydrates produced by T. rotula were differently used by the bacteria in the two and three partner treatments and most efficiently by the enriched natural bacterial community. Our approach is a valid model system to study physico-chemical bacteria-diatom interactions with increasing complexity.

Does ecological drift explain variation in microbiome composition among groups in a social host species?

Within a given species, considerable inter-individual, spatial, and temporal variation in the composition of the host microbiome exists. In group-living animals, social interactions homogenize microbiome composition among group members, nevertheless divergence in microbiome composition among related groups arise. Such variation can result from deterministic and stochastic processes. Stochastic changes, or ecological drift, can occur among symbionts with potential for colonizing a host and within individual hosts, and drive divergence in microbiome composition among hosts or host groups. We tested whether ecological drift associated with dispersal and foundation of new groups cause divergence in microbiome composition between natal and newly formed groups in the social spider Stegodyphus dumicola. We simulated the initiation of new groups by splitting field-collected nests into groups of 1, 3, and 10 individuals respectively, and compared variation in microbiome composition among and within groups after 6 weeks using 16S rRNA gene sequencing. Theory predicts that ecological drift increases with decreasing group size. We found that microbiome composition among single founders was more dissimilar than among individuals kept in groups, supporting this prediction. Divergence in microbiome composition from the natal nest was mainly driven by a higher number of non-core symbionts. This suggests that stochastic divergence in host microbiomes can arise during the process of group formation by individual founders, which could explain the existence of among-group variation in microbiome composition in the wild. Individual founders appear to harbour higher relative abundances of non-core symbionts compared with founders in small groups, some of which are possible pathogens. These symbionts vary in occurrence with group size, indicating that group dynamics influence various core and non-core symbionts differently.

A Gulf Stream frontal eddy harbors a distinct microbiome compared to adjacent waters

Mesoscale oceanographic features, including eddies, have the potential to alter productivity and other biogeochemical rates in the ocean. Here, we examine the microbiome of a cyclonic, Gulf Stream frontal eddy, with a distinct origin and environmental parameters compared to surrounding waters, in order to better understand the processes dominating microbial community assembly in the dynamic coastal ocean. Our microbiome-based approach identified the eddy as distinct from the surround Gulf Stream waters. The eddy-associated microbial community occupied a larger area than identified by temperature and salinity alone, increasing the predicted extent of eddy-associated biogeochemical processes. While the eddy formed on the continental shelf, after two weeks both environmental parameters and microbiome composition of the eddy were most similar to the Gulf Stream, suggesting the effect of environmental filtering on community assembly or physical mixing with adjacent Gulf Stream waters. In spite of the potential for eddy-driven upwelling to introduce nutrients and stimulate primary production, eddy surface waters exhibit lower chlorophyll a along with a distinct and less even microbial community, compared to the Gulf Stream. At the population level, the eddy microbiome exhibited differences among the cyanobacteria (e.g. lower Trichodesmium and higher Prochlorococcus) and in the heterotrophic alpha Proteobacteria (e.g. lower relative abundances of specific SAR11 phylotypes) versus the Gulf Stream. However, better delineation of the relative roles of processes driving eddy community assembly will likely require following the eddy and surrounding waters since inception. Additionally, sampling throughout the water column could better clarify the contribution of these mesoscale features to primary production and carbon export in the oceans.

Microbiome and photoperiod interactively determine thermal sensitivity of polar and temperate diatoms

The effect of temperature on ectothermic organisms in the context of climate change has long been considered in isolation (i.e. as a single driver). This is challenged by observations demonstrating that temperature-dependent growth is correlated to further factors. However, little is known how the chronobiological history of an organism reflected in its adaptation to re-occurring cyclic patterns in its environment (e.g. annual range of photoperiods in its habitat) and biotic interactions with its microbiome, contribute to shaping its realized niche. To address this, we conducted a full-factorial microcosm multi-stressor experiment with the marine diatoms Thalassiosira gravida (polar) and Thalassiosira rotula (temperate) across multiple levels of temperature (4°C; 9°C; 13.5°C) and photoperiod (4 h; 16 h; 24 h), both in the presence or absence of their microbiomes. While temperature-dependent growth of the temperate diatom was constrained by short and long photoperiods, the polar diatom coped with a 24 h photoperiod up to its thermal optimum (9°C). The algal microbiomes particularly supported host growth at the margins of their respective fundamental niches except for the combination of the warmest temperature tested at 24 h photoperiod. Overall, this study demonstrates that temperature tolerances may have evolved interactively and that the mutualistic effect of the microbiome can only be determined once the multifactorial abiotic niche is defined.

Interaction between genetics and inulin affects host metabolism in rainbow trout fed a sustainable all plant-based diet

Inulin affects nutrition and metabolism in many animals. Although inulin is widely used in the diet of teleosts, its mechanism of action is unknown. Here, we investigated the effect of inulin (2 %) on the intestinal microbiome and metabolism in rainbow trout (Oncorhynchus mykiss) selected for growth and survival when fed a 100 % plant-based diet (suave) and a control line (temoin). Metabolic responses to the two factors (line and inulin) in liver, intestine, muscle and adipose were tissue-specific, with line and interaction between the two factors influencing overall expression in liver. In the intestine, inulin and line and in muscle, line influenced the expression of metabolic genes. Microbiota between the mucus and digestive contents was significantly different, with genera from Proteobacteria being more abundant in the mucus, whereas genera from the Firmicutes and Planctomycetes being more abundant in contents. Effect of inulin and interaction between factors on the microbiome was evident in contents. The significant taxa of control and inulin-fed groups differed greatly with Streptococcus and Weissella being significantly abundant in the inulin-fed group. There was a general trend showing higher levels of all SCFA in temoin group with propionic acid levels being significantly higher. An operational taxonomic unit (OTU) belonging to the Ruminococcaceae was significantly abundant in suave. The tissue-specific correlations between OTU and gene expression may indicate the link between microbiome and metabolism. Together, these results suggest that line and inulin impact the gene expression in a tissue-specific manner, possibly driven by specific OTUs enriched in inulin-fed groups and suave.

Antibiotic-Induced Recruitment of Specific Algae-Associated Microbiome Enhances the Adaptability of Chlorella vulgaris to Antibiotic Stress and Incidence of Antibiotic Resistance

Insights into the symbiotic relation between eukaryotic hosts and their microbiome lift the curtain on the crucial roles of microbes in host fitness, behavior, and ecology. However, it remains unclear whether and how abiotic stress shapes the microbiome and further affects host adaptability. This study first investigated the effect of antibiotic exposure on behavior across varying algae taxa at the community level. Chlorophyta, in particular Chlorella vulgaris, exhibited remarkable adaptability to antibiotic stress, leading to their dominance in phytoplankton communities. Accordingly, we isolated C. vulgaris strains and compared the growth of axenic and nonaxenic ones under antibiotic conditions. The positive roles of antibiotics in algal growth were apparent only in the presence of bacteria. Results of 16S rRNA sequencing further revealed that antibiotic challenges resulted in the recruitment of specific bacterial consortia in the phycosphere, whose functions were tightly linked to the host growth promotion and adaptability enhancement. In addition, the algal phycosphere was characterized with 47-fold higher enrichment capability of antibiotic resistance genes (ARGs) than the surrounding water. Under antibiotic stress, specific ARG profiles were recruited in C. vulgaris phycosphere, presumably driven by the specific assembly of bacterial consortia and mobile genetic elements induced by antibiotics. Moreover, the antibiotics even enhanced the dissemination potential of the bacteria carrying ARGs from the algal phycosphere to broader environmental niches. Overall, this study provides an in-depth understanding into the potential functional significance of antibiotic-mediated recruitment of specific algae-associated bacteria for algae adaptability and ARG proliferation in antibiotic-polluted waters.

Emergent community architecture despite distinct diversity in the global whale shark (Rhincodon typus) epidermal microbiome

Microbiomes confer beneficial physiological traits to their host, but microbial diversity is inherently variable, challenging the relationship between microbes and their contribution to host health. Here, we compare the diversity and architectural complexity of the epidermal microbiome from 74 individual whale sharks (Rhincodon typus) across five aggregations globally to determine if network properties may be more indicative of the microbiome-host relationship. On the premise that microbes are expected to exhibit biogeographic patterns globally and that distantly related microbial groups can perform similar functions, we hypothesized that microbiome co-occurrence patterns would occur independently of diversity trends and that keystone microbes would vary across locations. We found that whale shark aggregation was the most important factor in discriminating taxonomic diversity patterns. Further, microbiome network architecture was similar across all aggregations, with degree distributions matching Erdos-Renyi-type networks. The microbiome-derived networks, however, display modularity indicating a definitive microbiome structure on the epidermis of whale sharks. In addition, whale sharks hosted 35 high-quality metagenome assembled genomes (MAGs) of which 25 were present from all sample locations, termed the abundant 'core'. Two main MAG groups formed, defined here as Ecogroup 1 and 2, based on the number of genes present in metabolic pathways, suggesting there are at least two important metabolic niches within the whale shark microbiome. Therefore, while variability in microbiome diversity is high, network structure and core taxa are inherent characteristics of the epidermal microbiome in whale sharks. We suggest the host-microbiome and microbe-microbe interactions that drive the self-assembly of the microbiome help support a functionally redundant abundant core and that network characteristics should be considered when linking microbiomes with host health.

The microbiome of the dinoflagellate Prorocentrum cordatum in laboratory culture and its changes at higher temperatures

In the ocean, phytoplankton are dependent on communities of bacteria living in the phycosphere, a hot spot of metabolic and genetic exchange. Many types of interactions between phytoplankton and phycosphere bacteria have been shown, but it is unclear if the microbial communities associated with microalgae strains in culture collections are beneficial or harmful to the host strain. Here, we studied the microbial communities associated with four strains of the dinoflagellate Prorocentrum cordatum that had been isolated from distant geographical locations and maintained in culture collection for hundreds of generations. Community composition was determined by 16S rRNA gene amplicon sequencing. The dinoflagellate host strain was the strongest parameter separating communities, while growth phase, lifestyle (particle-attached versus free-living) and temperature had only a modulating effect. Although the strains had been isolated from distant locations in the Atlantic and Pacific Ocean, 14 ASVs were shared among all strains, the most abundant ones being Gilvibacter, Marivita, uncultivated Rhodobacteraceae, Marinobacter, Hyphomonadaceae, Cupriavidus, Variovorax, and Paucibacter. Adaptation to higher temperatures resulted in specific changes in each phycosphere microbiome, including increased abundance of rare community members. We then compared the growth of the four xenic cultures to that of the axenic P. cordatum CCMP1329. At 20°C, growth of the xenic cultures was similar or slower than that of CCMP1329. At 26°C, all four xenic cultures experienced a death phase, while the axenic culture stably remained in the stationary phase. At 30°C, only two of the xenic cultures were able to grow. A shift of dinoflagellate metabolism from autotrophy to mixotrophy and competition between dinoflagellate and bacteria for limiting nutrients, including essential vitamins, may contribute to these differences in growth patterns.

Selection constrains lottery assembly in the microbiomes of closely related diatom species

It is generally recognised that interactions between microalgae and bacteria play an important role in the functioning of marine ecosystems. In this context, increasing attention is paid to the processes that shape microalga-associated microbiomes. In recent years, conflicting evidence has been reported with respect to the relative importance of selective vs neutral processes in the assembly process. Whereas some studies report strong selection imposed by the host, others propose a more neutral, lottery-like assembly model according to which the chance of bacteria becoming part of the microbiome is proportional to their abundance in the environment and not driven by the selectional pressure created by the host. In the present study, we investigated to what degree selective vs neutral assembly processes constrain taxonomic, phylogenetic and functional variation within and between microbiomes associated with 69 isolates belonging to the Cylindrotheca closterium benthic marine diatom complex. The diatom cultures were initiated from non-axenic clonal isolates from different marine environments and geographic locations, and were then reared in a common garden (lab) environment. An important environmental imprint, likely due to in situ lottery dynamics, was apparent in the diatom microbiomes. However, microbiome assembly was also phylogenetically and functionally constrained through selective filtering related to the host microhabitat. Randomised microbiome assembly simulations revealed evidence for phylogenetic overdispersion in the observed microbiomes, reflecting an important role in the assembly process for competition between bacteria on the one hand and predominantly genetically driven differences between the hosts on the other hand. Our study thus shows that even between closely related diatom strains, host selection affects microbiome assembly, superimposing the predominantly stochastically driven recruitment process.