How can microbial population genomics inform community ecology?

Recurring seasonality exposes dominant species and niche partitioning strategies of open ocean picoeukaryotic algae

Ocean spring phytoplankton blooms are dynamic periods important to global primary production. We document vertical patterns of a diverse suite of eukaryotic algae, the prasinophytes, in the North Atlantic Subtropical Gyre with monthly sampling over four years at the Bermuda Atlantic Time-series Study site. Water column structure was used to delineate seasonal stability periods more ecologically relevant than seasons defined by calendar dates. During winter mixing, tiny prasinophytes dominated by Class II comprise 46  ±  24% of eukaryotic algal (plastid-derived) 16S rRNA V1-V2 amplicons, specifically Ostreococcus Clade OII, Micromonas commoda, and Bathycoccus calidus. In contrast, Class VII are rare and Classes I and VI peak during warm stratified periods when surface eukaryotic phytoplankton abundances are low. Seasonality underpins a reservoir of genetic diversity from multiple prasinophyte classes during warm periods that harbor ephemeral taxa. Persistent Class II sub-species dominating the winter/spring bloom period retreat to the deep chlorophyll maximum in summer, poised to seed the mixed layer upon winter convection, exposing a mechanism for initiating high abundances at bloom onset. Comparisons to tropical oceans reveal broad distributions of the dominant sub-species herein. This unparalleled window into temporal and spatial niche partitioning of picoeukaryotic primary producers demonstrates how key prasinophytes prevail in warm oceans.

Adaptations to nitrogen availability drive ecological divergence of chemosynthetic symbionts

Bacterial symbionts, with their shorter generation times and capacity for horizontal gene transfer (HGT), play a critical role in allowing marine organisms to cope with environmental change. The closure of the Isthmus of Panama created distinct environmental conditions in the Tropical Eastern Pacific (TEP) and Caribbean, offering a "natural experiment" for studying how closely related animals evolve and adapt under environmental change. However, the role of bacterial symbionts in this process is often overlooked. We sequenced the genomes of endosymbiotic bacteria in two sets of sister species of chemosymbiotic bivalves from the genera Codakia and Ctena (family Lucinidae) collected on either side of the Isthmus, to investigate how differing environmental conditions have influenced the selection of symbionts and their metabolic capabilities. The lucinid sister species hosted different Candidatus Thiodiazotropha symbionts and only those from the Caribbean had the genetic potential for nitrogen fixation, while those from the TEP did not. Interestingly, this nitrogen-fixing ability did not correspond to symbiont phylogeny, suggesting convergent evolution of nitrogen fixation potential under nutrient-poor conditions. Reconstructing the evolutionary history of the nifHDKT operon by including other lucinid symbiont genomes from around the world further revealed that the last common ancestor (LCA) of Ca. Thiodiazotropha lacked nif genes, and populations in oligotrophic habitats later re-acquired the nif operon through HGT from the Sedimenticola symbiont lineage. Our study suggests that HGT of the nif operon has facilitated niche diversification of the globally distributed Ca. Thiodiazotropha endolucinida species clade. It highlights the importance of nitrogen availability in driving the ecological diversification of chemosynthetic symbiont species and the role that bacterial symbionts may play in the adaptation of marine organisms to changing environmental conditions.

Duplicated antibiotic resistance genes reveal ongoing selection and horizontal gene transfer in bacteria

Horizontal gene transfer (HGT) and gene duplication are often considered as separate mechanisms driving the evolution of new functions. However, the mobile genetic elements (MGEs) implicated in HGT can copy themselves, so positive selection on MGEs could drive gene duplications. Here, we use a combination of modeling and experimental evolution to examine this hypothesis and use long-read genome sequences of tens of thousands of bacterial isolates to examine its generality in nature. Modeling and experiments show that antibiotic selection can drive the evolution of duplicated antibiotic resistance genes (ARGs) through MGE transposition. A key implication is that duplicated ARGs should be enriched in environments associated with antibiotic use. To test this, we examined the distribution of duplicated ARGs in 18,938 complete bacterial genomes with ecological metadata. Duplicated ARGs are highly enriched in bacteria isolated from humans and livestock. Duplicated ARGs are further enriched in an independent set of 321 antibiotic-resistant clinical isolates. Our findings indicate that duplicated genes often encode functions undergoing positive selection and horizontal gene transfer in microbial communities.

Genomic basis of environmental adaptation in the widespread poly-extremophilic Exiguobacterium group

Delineating cohesive ecological units and determining the genetic basis for their environmental adaptation are among the most important objectives in microbiology. In the last decade, many studies have been devoted to characterizing the genetic diversity in microbial populations to address these issues. However, the impact of extreme environmental conditions, such as temperature and salinity, on microbial ecology and evolution remains unclear so far. In order to better understand the mechanisms of adaptation, we studied the (pan)genome of Exiguobacterium, a poly-extremophile bacterium able to grow in a wide range of environments, from permafrost to hot springs. To have the genome for all known Exiguobacterium type strains, we first sequenced those that were not yet available. Using a reverse-ecology approach, we showed how the integration of phylogenomic information, genomic features, gene and pathway enrichment data, regulatory element analyses, protein amino acid composition, and protein structure analyses of the entire Exiguobacterium pangenome allows to sharply delineate ecological units consisting of mesophilic, psychrophilic, halophilic-mesophilic, and halophilic-thermophilic ecotypes. This in-depth study clarified the genetic basis of the defined ecotypes and identified some key mechanisms driving the environmental adaptation to extreme environments. Our study points the way to organizing the vast microbial diversity into meaningful ecologically units, which, in turn, provides insight into how microbial communities adapt and respond to different environmental conditions in a changing world.

The active free-living bathypelagic microbiome is largely dominated by rare surface taxa

A persistent microbial seed bank is postulated to sustain the marine biosphere, and recent findings show that prokaryotic taxa present in the ocean's surface dominate prokaryotic communities throughout the water column. Yet, environmental conditions exert a tight control on the activity of prokaryotes, and drastic changes in these conditions are known to occur from the surface to deep waters. The simultaneous characterization of the total (DNA) and active (i.e. with potential for protein synthesis, RNA) free-living communities in 13 stations distributed across the tropical and subtropical global ocean allowed us to assess their change in structure and diversity along the water column. We observed that active communities were surprisingly more similar along the vertical gradient than total communities. Looking at the vertical connectivity of the active vs. the total communities, we found that taxa detected in the surface sometimes accounted for more than 75% of the active microbiome of bathypelagic waters (50% on average). These active taxa were generally rare in the surface, representing a small fraction of all the surface taxa. Our findings show that the drastic vertical change in environmental conditions leads to the inactivation and disappearance of a large proportion of surface taxa, but some surface-rare taxa remain active (or with potential for protein synthesis) and dominate the bathypelagic active microbiome.

Russian Biodiversity Collections: A Professional Opinion Survey

Biodiversity collections are important vehicles for protecting endangered wildlife in situations of adverse anthropogenic influence. In Russia, there are currently a number of institution- and museum-based biological collections, but there are no nation-wide centres of biodiversity collections. In this paper, we report on the results of our survey of 324 bioconservation, big-data, and ecology specialists from different regions of Russia in regard to the necessity to create several large national biodiversity centres of wildlife protection. The survey revealed specific goals that have to be fulfilled during the development of these centres for the protection and restoration of endangered wildlife species. The top three problems/tasks (topics) are the following: (1) the necessity to create large national centres for different types of specimens; (2) the full sequencing and creation of different "omic" (genomic, proteomic, transcriptomic, etc.) databases; (3) full digitisation of a biodiversity collection/centre. These goals may constitute a guideline for the future of biodiversity collections in Russia that would be targeted at protecting and restoring endangered species. With the due network service level, the translation of the website into English, and permission from the regulator (Ministry of Science and Higher Education of Russian Federation), it can also become an international project.

Vibrio cyclitrophicus population-specific biofilm formation and epibiotic growth on marine copepods

Vibrio spp. form a part of the microbiome of copepods-an abundant component of marine mesozooplankton. The biological mechanisms of the Vibrio-copepod association are largely unknown. In this study we compared biofilm formation of V. cyclitrophicus isolated from copepods (L-strains related to other particle-associated strains) and closely related strains originating from seawater (S-strains), and visualized and quantified their attachment and growth on copepods. The S- and L-strains formed similar biofilms in the presence of complete sea salts, suggesting previously unknown biofilm mechanisms in the S-strains. No biofilms formed if sodium chloride was present as the only salt but added calcium significantly enhanced biofilms in the L-strains. GFP-L-strain cells attached to live copepods at higher numbers than the S-strains, suggesting distinct mechanisms, potentially including calcium, support their colonization of copepods. The cells grew on live copepods after attachment, demonstrating that copepods sustain epibiotic V. cyclitrophicus growth in situ. The results demonstrate that in spite of their 99.1% average nucleotide identity, these V. cyclitrophicus strains have a differential capacity to colonize marine copepods. The introduced V. cyclitrophicus-A. tonsa model could be informative in future studies on Vibrio-copepod association.

Discordant population structure among rhizobium divided genomes and their legume hosts

Symbiosis often occurs between partners with distinct life history characteristics and dispersal mechanisms. Many bacterial symbionts have genomes comprising multiple replicons with distinct rates of evolution and horizontal transmission. Such differences might drive differences in population structure between hosts and symbionts and among the elements of the divided genomes of bacterial symbionts. These differences might, in turn, shape the evolution of symbiotic interactions and bacterial evolution. Here we use whole genome resequencing of a hierarchically structured sample of 191 strains of Sinorhizobium meliloti collected from 21 locations in southern Europe to characterize population structures of this bacterial symbiont, which forms a root nodule symbiosis with the host plant Medicago truncatula. S. meliloti genomes showed high local (within-site) variation and little isolation by distance. This was particularly true for the two symbiosis elements, pSymA and pSymB, which have population structures that are similar to each other, but distinct from both the bacterial chromosome and the host plant. Given limited recombination on the chromosome, compared to the symbiosis elements, distinct population structures may result from differences in effective gene flow. Alternatively, positive or purifying selection, with little recombination, may explain distinct geographical patterns at the chromosome. Discordant population structure between hosts and symbionts indicates that geographically and genetically distinct host populations in different parts of the range might interact with genetically similar symbionts, potentially minimizing local specialization.

Genome sequencing and multifaceted taxonomic analysis of novel strains of violacein-producing bacteria and non-violacein-producing close relatives

Violacein is a water-insoluble violet pigment produced by various Gram-negative bacteria. The compound and the bacteria that produce it have been gaining attention due to the antimicrobial and proposed antitumour properties of violacein and the possibility that strains producing it may have broad industrial uses. Bacteria that produce violacein have been isolated from diverse environments including fresh and ocean waters, glaciers, tropical soils, trees, fish and the skin of amphibians. We report here the isolation and characterization of six violacein-producing bacterial strains and three non-violacein-producing close relatives, each isolated from either an aquatic environment or moist food materials in northern California, USA. For each isolate, we characterized traditional phenotypes, generated and analysed draft genome sequences, and carried out multiple types of taxonomic, phylogenetic and phylogenomic analyses. Based on these analyses we assign putative identifications to the nine isolates, which include representatives of the genera Chromobacterium, Aquitalea, Iodobacter, Duganella, Massilia and Janthinobacterium. In addition, we discuss the utility of various metrics for taxonomic assignment in these groups including average nucleotide identity, whole genome phylogenetic analysis and extent of recent homologous recombination using the software program PopCOGenT.

Differential carbon utilization enables co-existence of recently speciated Campylobacteraceae in the cow rumen epithelial microbiome

The activities of different microbes in the cow rumen have been shown to modulate the host's ability to utilize plant biomass, while the host-rumen interface has received little attention. As datasets collected worldwide have pointed to Campylobacteraceae as particularly abundant members of the rumen epithelial microbiome, we targeted this group in a subset of seven cows with meta- and isolate genome analysis. We show that the dominant Campylobacteraceae lineage has recently speciated into two populations that were structured by genome-wide selective sweeps followed by population-specific gene import and recombination. These processes led to differences in gene expression and enzyme domain composition that correspond to the ability to utilize acetate, the main carbon source for the host, at the cost of inhibition by propionate. This trade-off in competitive ability further manifests itself in differential dynamics of the two populations in vivo. By exploring population-level adaptations that otherwise remain cryptic in culture-independent analyses, our results highlight how recent evolutionary dynamics can shape key functional roles in the rumen microbiome.

Functional responses of key marine bacteria to environmental change – toward genetic counselling for coastal waters

Coastal ecosystems deteriorate globally due to human-induced stress factors, like nutrient loading and pollution. Bacteria are critical to marine ecosystems, e.g., by regulating nutrient cycles, synthesizing vitamins, or degrading pollutants, thereby providing essential ecosystem services ultimately affecting economic activities. Yet, until now bacteria are overlooked both as mediators and indicators of ecosystem health, mainly due to methodological limitations in assessing bacterial ecosystem functions. However, these limitations are largely overcome by the advances in molecular biology and bioinformatics methods for characterizing the genetics that underlie functional traits of key bacterial populations – “key” in providing important ecosystem services, being abundant, or by possessing high metabolic rates. It is therefore timely to analyze and define the functional responses of bacteria to human-induced effects on coastal ecosystem health. We posit that categorizing the responses of key marine bacterial populations to changes in environmental conditions through modern microbial oceanography methods will allow establishing the nascent field of genetic counselling for our coastal waters. This requires systematic field studies of linkages between functional traits of key bacterial populations and their ecosystem functions in coastal seas, complemented with systematic experimental analyses of the responses to different stressors. Research and training in environmental management along with dissemination of results and dialogue with societal actors are equally important to ensure the role of bacteria is understood as fundamentally important for coastal ecosystems. Using the responses of microorganisms as a tool to develop genetic counselling for coastal ecosystems can ultimately allow for integrating bacteria as indicators of environmental change.

Gene flow and introgression are pervasive forces shaping the evolution of bacterial species

BACKGROUND

Although originally thought to evolve clonally, studies have revealed that most bacteria exchange DNA. However, it remains unclear to what extent gene flow shapes the evolution of bacterial genomes and maintains the cohesion of species.

RESULTS

Here, we analyze the patterns of gene flow within and between >2600 bacterial species. Our results show that fewer than 10% of bacterial species are truly clonal, indicating that purely asexual species are rare in nature. We further demonstrate that the taxonomic criterion of ~95% genome sequence identity routinely used to define bacterial species does not accurately represent a level of divergence that imposes an effective barrier to gene flow across bacterial species. Interruption of gene flow can occur at various sequence identities across lineages, generally from 90 to 98% genome identity. This likely explains why a ~95% genome sequence identity threshold has empirically been judged as a good approximation to define bacterial species. Our results support a universal mechanism where the availability of identical genomic DNA segments required to initiate homologous recombination is the primary determinant of gene flow and species boundaries in bacteria. We show that these barriers of gene flow remain porous since many distinct species maintain some level of gene flow, similar to introgression in sexual organisms.

CONCLUSIONS

Overall, bacterial evolution and speciation are likely shaped by similar forces driving the evolution of sexual organisms. Our findings support a model where the interruption of gene flow-although not necessarily the initial cause of speciation-leads to the establishment of permanent and irreversible species borders.

Prochlorococcus have low global mutation rate and small effective population size

Prochlorococcus are the most abundant free-living photosynthetic carbon-fixing organisms in the ocean. Prochlorococcus show small genome sizes, low genomic G+C content, reduced DNA repair gene pool and fast evolutionary rates, which are typical features of endosymbiotic bacteria. Nevertheless, their evolutionary mechanisms are believed to be different. Evolution of endosymbiotic bacteria is dominated by genetic drift owing to repeated population bottlenecks, whereas Prochlorococcus are postulated to have extremely large effective population sizes (N_e) and thus drift has rarely been considered. However, accurately extrapolating N_e requires measuring an unbiased global mutation rate through mutation accumulation, which is challenging for Prochlorococcus. Here, we managed this experiment over 1,065 days using Prochlorococcus marinus AS9601, sequenced genomes of 141 mutant lines and determined its mutation rate to be 3.50 × 10^-10 per site per generation. Extrapolating N_e additionally requires identifying population boundaries, which we defined using PopCOGenT and over 400 genomes related to AS9601. Accordingly, we calculated its N_e to be 1.68 × 10⁷, which is only reasonably greater than that of endosymbiotic bacteria but surprisingly smaller than that of many free-living bacteria extrapolated using the same approach. Our results therefore suggest that genetic drift is a key driver of Prochlorococcus evolution.

Seasonal niche differentiation among closely related marine bacteria

Bacteria display dynamic abundance fluctuations over time in marine environments, where they play key biogeochemical roles. Here, we characterized the seasonal dynamics of marine bacteria in a coastal oligotrophic time series station, tested how similar the temporal niche of closely related taxa is, and what are the environmental parameters modulating their seasonal abundance patterns. We further explored how conserved the niche is at higher taxonomic levels. The community presented recurrent patterns of seasonality for 297 out of 6825 amplicon sequence variants (ASVs), which constituted almost half of the total relative abundance (47%). For certain genera, niche similarity decreased as nucleotide divergence in the 16S rRNA gene increased, a pattern compatible with the selection of similar taxa through environmental filtering. Additionally, we observed evidence of seasonal differentiation within various genera as seen by the distinct seasonal patterns of closely related taxa. At broader taxonomic levels, coherent seasonal trends did not exist at the class level, while the order and family ranks depended on the patterns that existed at the genus level. This study identifies the coexistence of closely related taxa for some bacterial groups and seasonal differentiation for others in a coastal marine environment subjected to a strong seasonality.

Microcoleus (Cyanobacteria) form watershed-wide populations without strong gradients in population structure

The relative importance of separation by distance and by environment to population genetic diversity can be conveniently tested in river networks, where these two drivers are often independently distributed over space. To evaluate the importance of dispersal and environmental conditions in shaping microbial population structures, we performed genome‐resolved metagenomic analyses of benthic Microcoleus‐dominated cyanobacterial mats collected in the Eel and Russian River networks (California, USA). The 64 Microcoleus genomes were clustered into three species that shared >96.5% average nucleotide identity (ANI). Most mats were dominated by one strain, but minor alleles within mats were often shared, even over large spatial distances (>300 km). Within the most common Microcoleus species, the ANI between the dominant strains within mats decreased with increasing spatial separation. However, over shorter spatial distances (tens of kilometres), mats from different subwatersheds had lower ANI than mats from the same subwatershed, suggesting that at shorter spatial distances environmental differences between subwatersheds in factors like canopy cover, conductivity, and mean annual temperature decreases ANI. Since mats in smaller creeks had similar levels of nucleotide diversity (π) as mats in larger downstream subwatersheds, within‐mat genetic diversity does not appear to depend on the downstream accumulation of upstream‐derived strains. The four‐gamete test and sequence length bias suggest recombination occurs between almost all strains within each species, even between populations separated by large distances or living in different habitats. Overall, our results show that, despite some isolation by distance and environmental conditions, sufficient gene‐flow occurs among cyanobacterial strains to prevent either driver from producing distinctive population structures across the watershed.

The Rhizobial Microbiome from the Tropical Savannah Zones in Northern Côte d'Ivoire

Over the past decade, many projects have been initiated worldwide to decipher the composition and function of the soil microbiome, including the African Soil Microbiome (AfSM) project that aims at providing new insights into the presence and distribution of key groups of soil bacteria from across the African continent. In this national study, carried out under the auspices of the AfSM project, we assessed the taxonomy, diversity and distribution of rhizobial genera in soils from the tropical savannah zones in Northern Côte d’Ivoire. Genomic DNA extracted from seven sampled soils was analyzed by sequencing the V4-V5 variable region of the 16S rDNA using Illumina’s MiSeq platform. Subsequent bioinformatic and phylogenetic analyses showed that these soils harbored 12 out of 18 genera of Proteobacteria harboring rhizobia species validly published to date and revealed for the first time that the Bradyrhizobium genus dominates in tropical savannah soils, together with Microvirga and Paraburkholderia. In silico comparisons of different 16S rRNA gene variable regions suggested that the V5-V7 region could be suitable for differentiating rhizobia at the genus level, possibly replacing the use of the V4-V5 region. These data could serve as indicators for future rhizobial microbiome explorations and for land-use decision-making.

Population Genomics of Plant Viruses: The Ecology and Evolution of Virus Emergence

The genomics era has revolutionized studies of adaptive evolution by monitoring large numbers of loci throughout the genomes of many individuals. Ideally, the investigation of emergence in plant viruses requires examining the population dynamics of both virus and host, their interactions with each other, with other organisms and the abiotic environment. Genetic mechanisms that affect demographic processes are now being studied with high-throughput technologies, traditional genetics methods, and new computational tools for big-data. In this review, we discuss the utility of these approaches to monitor and detect changes in virus populations within cells and individuals, and over wider areas across species and communities of ecosystems. The advent of genomics in virology has fostered a multidisciplinary approach to tackling disease risk. The ability to make sense of the information now generated in this integrated setting is by far the most substantial obstacle to the ultimate goal of plant virology to minimize the threats to food security posed by disease. To achieve this goal, it is imperative to understand and forecast how populations respond to future changes in complex natural systems.

Differential recruitment of opportunistic taxa leads to contrasting abilities in carbon processing by bathypelagic and surface microbial communities

Different factors affect the way dissolved organic matter (DOM) is processed in the ocean water column, including environmental conditions and the functional capabilities of the communities. Recent studies have shown that bathypelagic prokaryotes are metabolically flexible, but whether this versatility translates into a higher ability to process DOM has been barely explored. Here we performed a multifactorial transplant experiment to compare the growth, activity and changes in DOM quality in surface and bathypelagic waters inoculated with either surface or bathypelagic prokaryotic communities. The effect of nutrient additions to surface waters was also explored. Despite no differences in the cell abundance of surface and deep ocean prokaryotes were observed in any of the treatments, in surface waters with nutrients the heterotrophic production of surface prokaryotes rapidly decreased. Conversely, bathypelagic communities displayed a sustained production throughout the experiment. Incubations with surface prokaryotes always led to a significant accumulation of recalcitrant compounds, which did not occur with bathypelagic prokaryotes, suggesting they have a higher ability to process DOM. These contrasting abilities could be explained by the recruitment of a comparatively larger number of opportunistic taxa within the bathypelagic assemblages, which likely resulted in a broader community capability of substrate utilization.

A multi-omics approach to solving problems in plant disease ecology

The swift rise of omics-approaches allows for investigating microbial diversity and plant-microbe interactions across diverse ecological communities and spatio-temporal scales. The environment, however, is rapidly changing. The introduction of invasive species and the effects of climate change have particular impact on emerging plant diseases and managing current epidemics. It is critical, therefore, to take a holistic approach to understand how and why pathogenesis occurs in order to effectively manage for diseases given the synergies of changing environmental conditions. A multi-omics approach allows for a detailed picture of plant-microbial interactions and can ultimately allow us to build predictive models for how microbes and plants will respond to stress under environmental change. This article is designed as a primer for those interested in integrating -omic approaches into their plant disease research. We review -omics technologies salient to pathology including metabolomics, genomics, metagenomics, volatilomics, and spectranomics, and present cases where multi-omics have been successfully used for plant disease ecology. We then discuss additional limitations and pitfalls to be wary of prior to conducting an integrated research project as well as provide information about promising future directions.

A core phyllosphere microbiome exists across distant populations of a tree species indigenous to New Zealand

The phyllosphere microbiome is increasingly recognised as an influential component of plant physiology, yet it remains unclear whether stable host-microbe associations generally exist in the phyllosphere. Leptospermum scoparium (mānuka) is a tea tree indigenous to New Zealand, and honey derived from mānuka is widely known to possess unique antimicrobial properties. However, the host physiological traits associated with these antimicrobial properties vary widely, and the specific cause of such variation has eluded scientists despite decades of research. Notably, the mānuka phyllosphere microbiome remains uncharacterised, and its potential role in mediating host physiology has not been considered. Working within the prevailing core microbiome conceptual framework, we hypothesise that the phyllosphere microbiome of mānuka exhibits specific host association patterns congruent with those of a microbial community under host selective pressure (null hypothesis: the mānuka phyllosphere microbiome is recruited stochastically from the surrounding environment). To examine our hypothesis, we characterised the phyllosphere and associated soil microbiomes of five distinct and geographically distant mānuka populations across the North Island of New Zealand. We identified a habitat-specific and relatively abundant core microbiome in the mānuka phyllosphere, which was persistent across all samples. In contrast, non-core phyllosphere microorganisms exhibited significant variation across individual host trees and populations that was strongly driven by environmental and spatial factors. Our results demonstrate the existence of a dominant and ubiquitous core microbiome in the phyllosphere of mānuka, supporting our hypothesis that phyllosphere microorganisms of mānuka exhibit specific host association and potentially mediate physiological traits of this nationally and culturally treasured indigenous plant. In addition, our results illustrate biogeographical patterns in mānuka phyllosphere microbiomes and offer insight into factors contributing to phyllosphere microbiome assembly.