Recruitment, loss and coexistence in a guild of territorial coral reef fishes

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.

Poisson hurdle model-based method for clustering microbiome features

MOTIVATION

High-throughput sequencing technologies have greatly facilitated microbiome research and have generated a large volume of microbiome data with the potential to answer key questions regarding microbiome assembly, structure and function. Cluster analysis aims to group features that behave similarly across treatments, and such grouping helps to highlight the functional relationships among features and may provide biological insights into microbiome networks. However, clustering microbiome data are challenging due to the sparsity and high dimensionality.

RESULTS

We propose a model-based clustering method based on Poisson hurdle models for sparse microbiome count data. We describe an expectation-maximization algorithm and a modified version using simulated annealing to conduct the cluster analysis. Moreover, we provide algorithms for initialization and choosing the number of clusters. Simulation results demonstrate that our proposed methods provide better clustering results than alternative methods under a variety of settings. We also apply the proposed method to a sorghum rhizosphere microbiome dataset that results in interesting biological findings.

AVAILABILITY AND IMPLEMENTATION

R package is freely available for download at https://cran.r-project.org/package=PHclust.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The microbiome of the habitat-forming brown alga Fucus vesiculosus (Phaeophyceae) has similar cross-Atlantic structure that reflects past and present drivers1

Latitudinal diversity gradients have provided many insights into species differentiation and community processes. In the well-studied intertidal zone, however, little is known about latitudinal diversity in microbiomes associated with habitat-forming hosts. We investigated microbiomes of Fucus vesiculosus because of deep understanding of this model system and its latitudinally large, cross-Atlantic range. Given multiple effects of photoperiod, we predicted that cross-Atlantic microbiomes of the Fucus microbiome would be similar at similar latitudes and correlate with environmental factors. We found that community structure and individual amplicon sequencing variants (ASVs) showed distinctive latitudinal distributions, but alpha diversity did not. Latitudinal differentiation was mostly driven by ASVs that were more abundant in cold temperate to subarctic (e.g., Granulosicoccus_t3260, Burkholderia/Caballeronia/Paraburkholderia_t8371) or warm temperate (Pleurocapsa_t10392) latitudes. Their latitudinal distributions correlated with different humidity, tidal heights, and air/sea temperatures, but rarely with irradiance or photoperiod. Many ASVs in potentially symbiotic genera displayed novel phylogenetic biodiversity with differential distributions among tissues and regions, including closely related ASVs with differing north-south distributions that correlated with Fucus phylogeography. An apparent southern range contraction of F. vesiculosus in the NW Atlantic on the North Carolina coast mimics that recently observed in the NE Atlantic. We suggest cross-Atlantic microbial structure of F. vesiculosus is related to a combination of past (glacial-cycle) and contemporary environmental drivers.

Exploration of marine bacterioplankton community assembly mechanisms during chemical dispersant and surfactant-assisted oil biodegradation

Assembly processes in marine microbial communities amended with crude oil and chemical dispersant are poorly understood and even more so when biosurfactants are used. We set up a microcosm experiment in which microbiome structure was analyzed using 16S rRNA gene amplicon sequencing and six null models to better understand and quantify the mechanisms and patterns controlling the assembly of a marine crude oil degrading microbial community in the presence of chemical dispersant or rhamnolipid biosurfactant. Although each null model quantifies different aspects of the community assembly, there was a general agreement that neither purely stochastic nor purely deterministic processes dominated the microbial communities, and their influence was variable over time. Determinism was dominant in the early phase of incubation, while stochasticity was prevalent in the middle and late stages. There was faster recruitment of phylogenetically distant species in the dispersant-amended community compared to oil-only or rhamnolipid-amended communities. This analysis provides important insights of how chemical dispersants and rhamnolipid influence microbial communities' dynamics and identified which groups may be excluded-an important consideration for biodegradation process and oil spill response.

Robustness analysis of metabolic predictions in algal microbial communities based on different annotation pipelines

Animals, plants, and algae rely on symbiotic microorganisms for their development and functioning. Genome sequencing and genomic analyses of these microorganisms provide opportunities to construct metabolic networks and to analyze the metabolism of the symbiotic communities they constitute. Genome-scale metabolic network reconstructions rest on information gained from genome annotation. As there are multiple annotation pipelines available, the question arises to what extent differences in annotation pipelines impact outcomes of these analyses. Here, we compare five commonly used pipelines (Prokka, MaGe, IMG, DFAST, RAST) from predicted annotation features (coding sequences, Enzyme Commission numbers, hypothetical proteins) to the metabolic network-based analysis of symbiotic communities (biochemical reactions, producible compounds, and selection of minimal complementary bacterial communities). While Prokka and IMG produced the most extensive networks, RAST and DFAST networks produced the fewest false positives and the most connected networks with the fewest dead-end metabolites. Our results underline differences between the outputs of the tested pipelines at all examined levels, with small differences in the draft metabolic networks resulting in the selection of different microbial consortia to expand the metabolic capabilities of the algal host. However, the consortia generated yielded similar predicted producible compounds and could therefore be considered functionally interchangeable. This contrast between selected communities and community functions depending on the annotation pipeline needs to be taken into consideration when interpreting the results of metabolic complementarity analyses. In the future, experimental validation of bioinformatic predictions will likely be crucial to both evaluate and refine the pipelines and needs to be coupled with increased efforts to expand and improve annotations in reference databases.

Clonal structure through space and time: High stability in the holothurian Stichopus chloronotus (Echinodermata)

Sea cucumbers are increasingly exploited for human consumption and for their curative properties, and many wild populations are now depleted or in danger of extinction. While aquaculture is seen as an alternative to fisheries and as a mean to restore wild populations, more knowledge is needed on their reproductive strategies to render this practice efficient, notably for fissiparous holothurians, which are some of the mobile animals able of asexual reproduction by transverse fission. Little information is available on their population genetic diversity and structure. Here, the clonal structure of populations of the fissiparous sea cucumber Stichopus chloronotus has been investigated using nine microsatellite loci and a random sampling, at different spatial (intra‐reef and inter‐reef) and temporal (inter‐season and inter‐year) scales. Our findings highlight the importance of asexual reproduction in maintaining these populations, and the prevalence of the “initial seedling recruitment” strategy (ISR), leading to a high stability of clonal composition over seasons and years. It also seemed that clonal propagation was limited to the reef scale (<10 km) while reefs were connected by sexual dispersal. This is the first time that clonal structure in sea cucumbers has been studied at such a fine scale, with a specific sampling strategy. It provides key findings on the genetic diversity and structure of fissiparous sea cucumbers, which will be useful for the management of wild populations and aquaculture.

Microbial community function in the bleaching disease of the marine macroalgae Delisea pulchra

Disease is increasingly viewed as a major factor impacting the health of both natural and cultured populations of marine organisms, including macroalgae. The red macroalga Delisea pulchra suffers from a bleaching disease resulting from host stress and infection by opportunistic bacterial pathogens. However, how pathogens cause the disease and how the entire macro algal-associated community is involved in the process is unclear. Here, we perform a metagenomic analysis of microbial communities associated with diseased and healthy D. pulchra across multiple bleaching events. Analysis of reconstructed 16S rRNA gene sequences showed that bacteria belonging to the families Rhodobacteraceae, Saprospiraceae and Flavobacteriaceae, including bacteria previously implicated in algal bleaching, to be enriched in diseased D. pulchra. Genes with predicted functions related to chemotaxis, motility, oxidative stress response, vitamin biosynthesis and nutrient acquisition were also prevalent in microbiomes of bleached algae, which may have a role in pathogenicity. Reconstruction of genomes that were abundant on bleached samples revealed that no single organism contains all bleaching-enriched functional genes. This observation indicates that potential virulence traits are distributed across multiple bacteria and that the disease in D. pulchra may result from a consortium of opportunistic pathogens, analogous to dysbiotic or polymicrobial diseases.

Immigration Rates during Population Density Reduction in a Coral Reef Fish

Although the importance of density-dependent dispersal has been recognized in theory, few empirical studies have examined how immigration changes over a wide range of densities. In a replicated experiment using a novel approach allowing within-site comparison, we examined changes in immigration rate following the gradual removal of territorial damselfish from a limited area within a much larger patch of continuous habitat. In all sites, immigration occurred at intermediate densities but did not occur before the start of removals and only rarely as density approached zero. In the combined data and in 5 of 7 sites, the number of immigrants was a hump-shaped function of density. This is the first experimental evidence for hump-shaped, density-dependent immigration. This pattern may be more widespread than previously recognized because studies over more limited density ranges have identified positive density dependence at low densities and negative density dependence at high densities. Positive density dependence at low density can arise from limits to the number of potential immigrants and from behavioral preferences for settling near conspecifics. Negative density dependence at high density can arise from competition for resources, especially high quality territories. The potential for non-linear effects of local density on immigration needs to be recognized for robust predictions of conservation reserve function, harvest impacts, pest control, and the dynamics of fragmented populations.

Co-occurrence Analysis of Microbial Taxa in the Atlantic Ocean Reveals High Connectivity in the Free-Living Bacterioplankton

We determined the taxonomic composition of the bacterioplankton of the epipelagic zone of the Atlantic Ocean along a latitudinal transect (51°S–47°N) using Illumina sequencing of the V5-V6 region of the 16S rRNA gene and inferred co-occurrence networks. Bacterioplankon community composition was distinct for Longhurstian provinces and water depth. Free-living microbial communities (between 0.22 and 3 μm) were dominated by highly abundant and ubiquitous taxa with streamlined genomes (e.g., SAR11, SAR86, OM1, Prochlorococcus) and could clearly be separated from particle-associated communities which were dominated by Bacteroidetes, Planktomycetes, Verrucomicrobia, and Roseobacters. From a total of 369 different communities we then inferred co-occurrence networks for each size fraction and depth layer of the plankton between bacteria and between bacteria and phototrophic micro-eukaryotes. The inferred networks showed a reduction of edges in the deepest layer of the photic zone. Networks comprised of free-living bacteria had a larger amount of connections per OTU when compared to the particle associated communities throughout the water column. Negative correlations accounted for roughly one third of the total edges in the free-living communities at all depths, while they decreased with depth in the particle associated communities where they amounted for roughly 10% of the total in the last part of the epipelagic zone. Co-occurrence networks of bacteria with phototrophic micro-eukaryotes were not taxon-specific, and dominated by mutual exclusion (~60%). The data show a high degree of specialization to micro-environments in the water column and highlight the importance of interdependencies particularly between free-living bacteria in the upper layers of the epipelagic zone.

Host-microbe interactions as a driver of acclimation to salinity gradients in brown algal cultures

Like most eukaryotes, brown algae live in association with bacterial communities that frequently have beneficial effects on their development. Ectocarpus is a genus of small filamentous brown algae, which comprises a strain that has recently colonized freshwater, a rare transition in this lineage. We generated an inventory of bacteria in Ectocarpus cultures and examined the effect they have on acclimation to an environmental change, that is, the transition from seawater to freshwater medium. Our results demonstrate that Ectocarpus depends on bacteria for this transition: cultures that have been deprived of their associated microbiome do not survive a transfer to freshwater, but restoring their microflora also restores the capacity to acclimate to this change. Furthermore, the transition between the two culture media strongly affects the bacterial community composition. Examining a range of other closely related algal strains, we observed that the presence of two bacterial operational taxonomic units correlated significantly with an increase in low salinity tolerance of the algal culture. Despite differences in the community composition, no indications were found for functional differences in the bacterial metagenomes predicted to be associated with algae in the salinities tested, suggesting functional redundancy in the associated bacterial community. Our study provides an example of how microbial communities may impact the acclimation and physiological response of algae to different environments, and thus possibly act as facilitators of speciation. It paves the way for functional examinations of the underlying host-microbe interactions, both in controlled laboratory and natural conditions.

Conyza canadensis suppresses plant diversity in its nonnative ranges but not at home: a transcontinental comparison

The impact of invasive species across their native and nonnative ranges is poorly quantified and this impedes a complete understanding of biological invasions. We compared the impact of the native North American plant, Conyza canadensis, which is invasive to Eurasia, on species richness at home and in a number of introduced regions through well replicated transcontinental field studies, glasshouse experiments and individual-based models. Our results demonstrated mostly negative relationships between C. canadensis abundance and native species richness in nonnative ranges, but either positive or no relationships in its native North American range. In glasshouse experiments, the total biomass of Conyza was suppressed more by species from its native range than by species from regions where it is nonnative, but the effects of Conyza on other species did not show a consistent biogeographical pattern. Finally, individual-based models led to the exclusion of Conyza from North American scenarios but to high abundances in scenarios with species from the nonnative ranges of Conyza. We illustrate biogeographical differences in the impact of an invader across regional scales and suggest that inherent differences in one specific aspect of competitive ability, tolerance to the effects of other species, may play some role in these differences.

Models of experimentally derived competitive effects predict biogeographical differences in the abundance of invasive and native plant species

Mono-dominance by invasive species provides opportunities to explore determinants of plant distributions and abundance; however, linking mechanistic results from small scale experiments to patterns in nature is difficult. We used experimentally derived competitive effects of an invader in North America, Acroptilon repens, on species with which it co-occurs in its native range of Uzbekistan and on species with which it occurs in its non-native ranges in North America, in individual-based models. We found that competitive effects yielded relative abundances of Acroptilon and other species in models that were qualitatively similar to those observed in the field in the two ranges. In its non-native range, Acroptilon can occur in nearly pure monocultures at local scales, whereas such nearly pure stands of Acroptilon appear to be much less common in its native range. Experimentally derived competitive effects of Acroptilon on other species predicted Acroptilon to be 4-9 times more proportionally abundant than natives in the North American models, but proportionally equal to or less than the abundance of natives in the Eurasian models. Our results suggest a novel way to integrate complex combinations of interactions simultaneously, and that biogeographical differences in the competitive effects of an invader correspond well with biogeographical differences in abundance and impact.

Interspecific territoriality in gibbons (Hylobates lar and H. pileatus) and its effects on the dynamics of interspecies contact zones

We investigated the ecology and interspecific interactions of the two gibbon species (Hylobates lar and H. pileatus) that overlap in distribution within a narrow zone of contact in the headwaters of the Takhong River at Khao Yai National Park in central Thailand. The zone is about 10-km wide, with phenotypic hybrids comprising 6.5% of the adult population. We compared species with respect to diet, territory size, intra- and interspecific group encounters, and territory quality. The two gibbon species exploited the same types of resources within their territories despite variation in the relative abundance of food-plant species between territories. The gibbons were interspecifically territorial, and males of both species displayed aggressive behaviors at both intraspecific and interspecific territorial boundaries. There was no marked difference in the amount of overlap between territories of conspecific and heterospecific pairs of groups. Although the habitat was not homogeneous, territory quality did not vary significantly between species. The species have not diverged in habitat preference or in diet. Neither species dominated in interspecific encounters, and both were reproducing well in the contact zone. We analyzed the potential consequences of several types of interspecific interactions on individual dispersal options and on the structure of the contact zone. Interference competition through interspecific territoriality affects the dispersal of individuals into the range of the other species. In general, territorial competition coupled with limited hybridization leads to predictions of a narrow contact zone or parapatry between species; thus, behavioral and ecological interactions between species need to be considered as potential factors in explaining range borders of primate species.

Bacterial community assembly based on functional genes rather than species

The principles underlying the assembly and structure of complex microbial communities are an issue of long-standing concern to the field of microbial ecology. We previously analyzed the community membership of bacterial communities associated with the green macroalga Ulva australis, and proposed a competitive lottery model for colonization of the algal surface in an attempt to explain the surprising lack of similarity in species composition across different algal samples. Here we extend the previous study by investigating the link between community structure and function in these communities, using metagenomic sequence analysis. Despite the high phylogenetic variability in microbial species composition on different U. australis (only 15% similarity between samples), similarity in functional composition was high (70%), and a core of functional genes present across all algal-associated communities was identified that were consistent with the ecology of surface- and host-associated bacteria. These functions were distributed widely across a variety of taxa or phylogenetic groups. This observation of similarity in habitat (niche) use with respect to functional genes, but not species, together with the relative ease with which bacteria share genetic material, suggests that the key level at which to address the assembly and structure of bacterial communities may not be "species" (by means of rRNA taxonomy), but rather the more functional level of genes.