Long-Horizon Associative Learning Explains Human Sensitivity to Statistical and Network Structures in Auditory Sequences.

Networks are a useful mathematical tool for capturing the complexity of the world. In a previous behavioral study, we showed that human adults were sensitive to the high-level network structure underlying auditory sequences, even when presented with incomplete information. Their performance was best explained by a mathematical model compatible with associative learning principles, based on the integration of the transition probabilities between adjacent and nonadjacent elements with a memory decay. In the present study, we explored the neural correlates of this hypothesis via magnetoencephalography (MEG). Participants (N = 23, 16 females) passively listened to sequences of tones organized in a sparse community network structure comprising two communities. An early difference (∼150 ms) was observed in the brain responses to tone transitions with similar transition probability but occurring either within or between communities. This result implies a rapid and automatic encoding of the sequence structure. Using time-resolved decoding, we estimated the duration and overlap of the representation of each tone. The decoding performance exhibited exponential decay, resulting in a significant overlap between the representations of successive tones. Based on this extended decay profile, we estimated a long-horizon associative learning novelty index for each transition and found a correlation of this measure with the MEG signal. Overall, our study sheds light on the neural mechanisms underlying human sensitivity to network structures and highlights the potential role of Hebbian-like mechanisms in supporting learning at various temporal scales.

A social perspective on perceived distances reveals deep community structure.

Community structure, including relationships between and within groups, is foundational to our understanding of the world around us. For dissimilarity-based data, leveraging social concepts of conflict and alignment, we provide an approach for capturing meaningful structural information resulting from induced local comparisons. In particular, a measure of local (community) depth is introduced that leads directly to a probabilistic partitioning conveying locally interpreted closeness (or cohesion). A universal choice of threshold for distinguishing strongly and weakly cohesive pairs permits consideration of both local and global structure. Cases in which one might benefit from use of the approach include data with varying density such as that arising as snapshots of complex processes in which differing mechanisms drive evolution locally. The inherent recalibrating in response to density allows one to sidestep the need for localizing parameters, common to many existing methods. Mathematical results together with applications in linguistics, cultural psychology, and genetics, as well as to benchmark clustering data have been included. Together, these demonstrate how meaningful community structure can be identified without additional inputs (e.g., number of clusters or neighborhood size), optimization criteria, iterative procedures, or distributional assumptions.

The role of fish predators and their foraging traits in shaping zooplankton community structure.

Differentiation of foraging traits among predator populations may help explain observed variation in the structure of prey communities. However, few studies have investigated the phenotypic effects of predators on their prey in natural communities. Here, we use a comparative analysis of 78 Greenlandic lakes to examine how foraging trait variation among threespine stickleback populations can help explain variation in zooplankton community composition among lakes. We find that landscape-scale variation in zooplankton composition was jointly explained by lake properties, such as size and water chemistry, and the presence and absence of both stickleback and arctic char. Additional variation in zooplankton community structure can be explained by stickleback jaw protrusion, a trait with known utility for foraging on zooplankton, but only in lakes where stickleback co-occur with arctic char. Overall, our results illustrate how trait variation of predators, alongside other ecosystem properties, can influence the composition of prey communities in nature.

Phylogeny structures species' interactions in experimental ecological communities.

Species' traits and interactions are products of evolutionary history. Despite the long-standing hypothesis that closely related species possess similar traits, and thus experience stronger competition, measuring the effect of evolutionary history on the ecology of natural communities remains challenging. We propose a novel framework to test whether phylogeny influences patterns of coexistence and abundance of species assemblages. In our approach, phylogenetic trees are used to parameterize species' interactions, which in turn determine the abundance of species in a given assemblage. We use likelihoods to score models parameterized with a given phylogeny, and contrast them with models built using random trees, allowing us to test whether phylogenetic information helps to predict species' abundances. Our statistical framework reveals that interactions are indeed structured by phylogeny in a large set of experimental plant communities. Our results confirm that evolutionary history can help predict, and potentially manage or conserve, the structure and function of complex ecological communities.

Patterns of genetic connectedness between modern and medieval Estonian genomes reveal the origins of a major ancestry component of the Finnish population.

The Finnish population is a unique example of a genetic isolate affected by a recent founder event. Previous studies have suggested that the ancestors of Finnic-speaking Finns and Estonians reached the circum-Baltic region by the 1st millennium BC. However, high linguistic similarity points to a more recent split of their languages. To study genetic connectedness between Finns and Estonians directly, we first assessed the efficacy of imputation of low-coverage ancient genomes by sequencing a medieval Estonian genome to high depth (23×) and evaluated the performance of its down-sampled replicas. We find that ancient genomes imputed from >0.1× coverage can be reliably used in principal-component analyses without projection. By searching for long shared allele intervals (LSAIs; similar to identity-by-descent segments) in unphased data for >143,000 present-day Estonians, 99 Finns, and 14 imputed ancient genomes from Estonia, we find unexpectedly high levels of individual connectedness between Estonians and Finns for the last eight centuries in contrast to their clear differentiation by allele frequencies. High levels of sharing of these segments between Estonians and Finns predate the demographic expansion and late settlement process of Finland. One plausible source of this extensive sharing is the 8th-10th centuries AD migration event from North Estonia to Finland that has been proposed to explain uniquely shared linguistic features between the Finnish language and the northern dialect of Estonian and shared Christianity-related loanwords from Slavic. These results suggest that LSAI detection provides a computationally tractable way to detect fine-scale structure in large cohorts.

Cortical timescales and the modular organization of structural and functional brain networks.

Recent years have seen growing interest in characterizing the properties of regional brain dynamics and their relationship to other features of brain structure and function. In particular, multiple studies have observed regional differences in the "timescale" over which activity fluctuates during periods of quiet rest. In the cerebral cortex, these timescales have been associated with both local circuit properties as well as patterns of inter-regional connectivity, including the extent to which each region exhibits widespread connectivity to other brain areas. In the current study, we build on prior observations of an association between connectivity and dynamics in the cerebral cortex by investigating the relationship between BOLD fMRI timescales and the modular organization of structural and functional brain networks. We characterize network community structure across multiple scales and find that longer timescales are associated with greater within-community functional connectivity and diverse structural connectivity. We also replicate prior observations of a positive correlation between timescales and structural connectivity degree. Finally, we find evidence for preferential functional connectivity between cortical areas with similar timescales. We replicate these findings in an independent dataset. These results contribute to our understanding of functional brain organization and structure-function relationships in the human brain, and support the notion that regional differences in cortical dynamics may in part reflect the topological role of each region within macroscale brain networks.

Differentiating spillover: an examination of cross-habitat movement in ecology spillover in ecology.

Organisms that immigrate into a recipient habitat generate a movement pattern that affects local population dynamics and the environment. Spillover is the pattern of unidirectional movement from a donor habitat to a different, adjacent recipient habitat. However, ecological definitions are often generalized to include any cross-habitat movement, which limits within- and cross-discipline collaboration. To assess spillover nomenclature, we reviewed 337 studies within the agriculture, disease, fisheries and habitat fragmentation disciplines. Each study's definition of spillover and the methods used were analysed. We identified four descriptors (movement, habitat type and arrangement, and effect) used that differentiate spillover from other cross-habitat movement patterns (dispersal, foray loops and edge movement). Studies often define spillover as movement (45%) but rarely measure it as such (4%), particularly in disease and habitat fragmentation disciplines. Consequently, 98% of studies could not distinguish linear from returning movement out of a donor habitat, which can overestimate movement distance. Overall, few studies (12%) included methods that matched their own definition, revealing a distinct mismatch. Because theory shows that long-term impacts of the different movement patterns can vary, differentiating spillover from other movement patterns is necessary for effective long-term and inter-disciplinary management of organisms that use heterogeneous landscapes.

Consequences of kelp forest ecosystem shifts and predictors of persistence through multiple stressors.

Ecological communities can be stable over multiple generations, or rapidly shift into structurally and functionally different configurations. In kelp forest ecosystems, overgrazing by sea urchins can abruptly shift forests into alternative states that are void of macroalgae and primarily dominated by actively grazing sea urchins. Beginning in 2014, a sea urchin outbreak along the central coast of California resulted in a patchy mosaic of remnant forests interspersed with sea urchin barrens. In this study, we used a 14-year subtidal monitoring dataset of invertebrates, algae, and fishes to explore changes in community structure associated with the loss of forests. We found that the spatial mosaic of barrens and forests resulted in a region-wide shift in community structure. However, the magnitude of kelp forest loss and taxonomic-level consequences were spatially heterogeneous. Taxonomic diversity declined across the region, but there were no declines in richness for any group, suggesting compositional redistribution. Baseline ecological and environmental conditions, and sea urchin behaviour, explained the persistence of forests through multiple stressors. These results indicate that spatial heterogeneity in preexisting ecological and environmental conditions can explain patterns of community change.

Continental scale comparison of mycobiomes in Parmelia and Peltigera lichens from Turkey and South Korea.

BACKGROUND: Lichens, traditionally considered as a simple partnership primarily between mycobiont and photobiont, are, in reality, complex holobionts comprised of a multitude of microorganisms. Lichen mycobiome represents fungal community residing within lichen thalli. While it is acknowledged that factors like the host lichen species and environmental conditions influence the structure of the lichen mycobiome, the existing research remains insufficient. To investigate which factor, host genus or location, has a greater impact on the lichen mycobiome, we conducted a comparative analysis of mycobiomes within Parmelia and Peltigera collected from both Turkey and South Korea, using high-throughput sequencing based on internal transcribed spacer region amplification. RESULTS: Overall, the lichen mycobiome was dominated by Capnodiales (Dothideomycetes), regardless of host or location. At the order level, the taxonomic composition was not significantly different according to lichen genus host or geographical distance. Hierarchical clustering of the top 100 abundant ASVs did not clearly indicate whether the lichen mycobiome was more influenced by host genus or location. Analyses of community similarity and partitioning variables revealed that the structure of the lichen mycobiome is more significantly influenced by location than by host genus. When analyzing the core mycobiome by host genus, the Peltigera mycobiome contained more ASV members than the Parmelia mycobiome. These two core mycobiomes also share common fungal strains, including basidiomycete yeast. Additionally, we used chi-squared tests to identify host genus-specialists and location-specialists. CONCLUSIONS: By comparing lichen mycobiomes of the same genera across different countries, our study advances our comprehension of these microbial communities. Our study elucidates that, although host species play a contributory role, geographic distance exerts a more pronounced impact on the structure of lichen mycobiome. We have made foundational contributions to understanding the lichen mycobiome occupying ecologically crucial niches. We anticipate that broader global-scale investigations into the fungal community structures will provide more detailed insights into fungal residents within lichens.

Neutral processes and taxonomic scale drive beta species-genetic diversity correlations in a submesophotic tropical reef fish.

If similar evolutionary forces maintain intra- and interspecific diversity, patterns of diversity at both levels of biological organization can be expected to covary across space. Although this prediction of a positive species-genetic diversity correlation (SGDC) has been tested for several taxa in natural landscapes, no study has yet evaluated the influence of the community delineation on these SGDCs. In this study, we focused on tropical fishes of the Indo-Pacific Ocean, using range-wide single nucleotide polymorphism data for a deep-sea fish (Etelis coruscans) and species presence data of 4878 Teleostei species. We investigated whether a diversity continuum occurred, for different community delineations (subfamily, family, order and class) and spatial extents, and which processes explained these diversity patterns. We found no association between genetic diversity and species richness (α-SGDC), regardless of the community and spatial extent. In contrast, we evidenced a positive relationship between genetic and species dissimilarities (ß-SGDC) when the community was defined at the subfamily or family level of the species of interest, and when the Western Indian Ocean was excluded. This relationship was related to the imprint of dispersal processes across levels of biological organization in Lutjanidae. However, this positive ß-SGDC was lost when considering higher taxonomic communities and at the scale of the entire Indo-Pacific, suggesting different responses of populations and communities to evolutionary processes at these scales. This study provides evidence that the taxonomic scale at which communities are defined and the spatial extent are pivotal to better understand the processes shaping diversity across levels of biological organization.

Shifts in soil ammonia-oxidizing community maintain the nitrogen stimulation of nitrification across climatic conditions.

Anthropogenic nitrogen (N) loading alters soil ammonia-oxidizing archaea (AOA) and bacteria (AOB) abundances, likely leading to substantial changes in soil nitrification. However, the factors and mechanisms determining the responses of soil AOA:AOB and nitrification to N loading are still unclear, making it difficult to predict future changes in soil nitrification. Herein, we synthesize 68 field studies around the world to evaluate the impacts of N loading on soil ammonia oxidizers and nitrification. Across a wide range of biotic and abiotic factors, climate is the most important driver of the responses of AOA:AOB to N loading. Climate does not directly affect the N-stimulation of nitrification, but does so via climate-related shifts in AOA:AOB. Specifically, climate modulates the responses of AOA:AOB to N loading by affecting soil pH, N-availability and moisture. AOB play a dominant role in affecting nitrification in dry climates, while the impacts from AOA can exceed AOB in humid climates. Together, these results suggest that climate-related shifts in soil ammonia-oxidizing community maintain the N-stimulation of nitrification, highlighting the importance of microbial community composition in mediating the responses of the soil N cycle to N loading.

Insight into physico-chemical properties and microbial community structure of biogas slurry from household biogas plants of sub-Himalaya for its implications in improved biogas production.

Numerous metagenomics studies, conducted in both full-scale anaerobic digesters and household biogas plants, have shed light on the composition and activity of microbial flora essential for optimizing the performance of biogas reactors, underscoring the significance of microbial community composition in biogas plant efficiency. Although the efficiency of household biogas plants in the sub-Himalayan region has been reported, there is no literature evidence on the microbial community structure of such household biogas plants in the sub-Himalayan region. The current study evaluated the physico-chemical properties and bacterial community structure from the slurry samples of household biogas plants prevalent in the sub-Himalayan region. The slurry samples were observed to be rich in nutrients; however, their carbon and nitrogen contents were higher than the recommended standard values of liquid-fermented organic manure. The species richness and diversity indices (Chao1, Shannon, and Simpson) of household biogas plants were quite similar to the advanced biogas reactors operating at mesophilic conditions. 16S rRNA gene amplicon sequencing reveals microbial diversity, showing a higher abundance of Firmicutes (70.9%) and Euryarchaeota (9.52%) in advanced biogas reactors compared to household biogas plants. Microbial analysis shows a lack of beneficial microbes for anaerobic digestion, which might be the reason for inefficient biogas production in household biogas plants of the sub-Himalayan region. The lack of efficient bacterial biomass may also be attributed to the digester design, feedstock, and ambient temperatures. This study emphasized the establishment of efficient microbial consortia for enhanced degradation rates that may increase the methane yield in biogas plants.

Serial cultures in invert emulsion and monophase systems for microbial community shaping and propagation.

BACKGROUND: Microbial communities harbor important biotechnological potential in diverse domains, however, the engineering and propagation of such communities still face both knowledge and know-how gaps. More specifically, culturing tools are needed to propagate and shape microbial communities, to obtain desired properties, and to exploit them. Previous work suggested that micro-confinement and segregation of microorganisms using invert (water-in-oil, w/o) emulsion broth can shape communities during propagation, by alleviating biotic interactions and inducing physiological changes in cultured bacteria. The present work aimed at evaluating invert emulsion and simple broth monophasic cultures for the propagation and shaping of bacterial communities derived from raw milk in a serial propagation design. RESULTS: The monophasic setup resulted in stable community structures during serial propagation, whereas the invert emulsion system resulted in only transiently stable structures. In addition, different communities with different taxonomic compositions could be obtained from a single inoculum. Furthermore, the implementation of invert emulsion systems has allowed for the enrichment of less abundant microorganisms and consequently facilitated their isolation on culture agar plates. CONCLUSIONS: The monophasic system enables communities to be propagated in a stable manner, whereas the invert emulsion system allowed for the isolation of less abundant microorganisms and the generation of diverse taxonomic compositions from a single inoculum.

Disentangling the role of environmental filtering and biotic resistance on alien invasions in a reservoir area.

Large-scale water conservancy projects benefit human life but have modified the landscape and provided opportunities for alien plant invasions. Understanding the environmental (e.g., climate), human-related (e.g., population density, proximity to human activities), and biotic (e.g., native plant, community structure) factors driving invasions is essential in the management of alien plants and biodiversity conservation in areas with intense human pressure. To this end, we investigated the spatial patterns of alien plant species distribution in the Three Gorges Reservoir Area (TGRA) of China and distinguished the role of the external environment and community characteristics in determining the occurrence of alien plants with differing levels of known invasion impacts in China using random forest analyses and structural equation models. A total of 102 alien plant species belonging to 30 families and 67 genera were recorded, the majority being annual and biennial herbs (65.7%). The results showed a negative diversity-invasibility relationship and supported the biotic resistance hypothesis. Moreover, percentage coverage of native plants was found to interact with native species richness and had a predominant role in resisting alien plant species. We found alien dominance was mainly the result of disturbance (e.g., changes in hydrological regime), which drove native plant loss. Our results also demonstrated that disturbance and temperature were more important for the occurrence of malignant invaders than all alien plants. Overall, our study highlights the importance of restoring diverse and productive native communities in resistance to invasion.

Mowing increased community stability in semiarid grasslands more than either fencing or grazing.

A substantial body of empirical evidence suggests that anthropogenic disturbance can affect the structure and function of grassland ecosystems. Despite this, few studies have elucidated the mechanisms through which grazing and mowing, the two most widespread land management practices, affect the stability of natural grassland communities. In this study, we draw upon 9 years of field data from natural grasslands in northern China to investigate the effects of gazing and mowing on community stability, specifically focusing on community aboveground net primary productivity (ANPP) and dominance, which are two major biodiversity mechanisms known to characterize community fluctuations. We found that both grazing and mowing reduced ANPP in comparison to areas enclosed by fencing. Grazing reduced community stability by increasing the likelihood of single-species dominance and decreasing the relative proportion of nondominant species. In contrast, mowing reduced the productivity of the dominant species but increased the productivity of nondominant species. As a consequence, mowing improved the overall community stability by increasing the stability of nondominant species. Our study provides novel insight into understanding of the relationship between community species fluctuation-stability, with implications for ecological research and ecosystem management in natural grasslands.

Responses of the mesozooplankton community to marine heatwaves: Challenges and solutions based on a long-term time series.

Marine heatwaves (MHWs) are extreme weather events that have major impacts on the structure and functioning of marine ecosystems worldwide. Due to anthropogenic climate change, the occurrence of MHWs is predicted to increase in future. There is already evidence linking MHWs with reductions in biodiversity and incidence of mass mortality events in coastal ecosystems. However, because MHWs are unpredictable, the quantification of their effects on communities is challenging. Here, we use the Helgoland Roads long-term time series (German Bight, North Sea), one of the richest marine time series in the world, and implement a modified before-after control-impact (BACI) design to evaluate MHW effect on mesozooplankton communities. Mesozooplankton play an essential role in connecting primary producers to higher trophic levels, and any changes in their community structure could have far-reaching impacts on the entire ecosystem. The responses of mesozooplankton community to MHWs in terms of community structure and densities occurred mainly in spring and autumn. Abundances of seven taxa, including some of the most abundant groups (e.g. copepods), were affected either positively or negatively in response to MHWs. In contrast, we observed no clear evidence of an impact of summer and winter MHWs; instead, the density of the most common taxa remained unchanged. Our results highlight the seasonally dependent impacts of MHWs on mesozooplankton communities and the challenges in evaluating those impacts. Long-term monitoring is an important contributor to the quantification of effects of MHWs on natural populations.

Non-genetic phenotypic variability affects populations and communities in protist microcosms.

Intraspecific trait variation (ITV), potentially driven by genetic and non-genetic mechanisms, can underlie variability in resource acquisition, individual fitness and ecological interactions. Impacts of ITV at higher levels of biological organizations are hence likely, but up-scaling our knowledge about ITV importance to communities and comparing its relative effects at population and community levels has rarely been investigated. Here, we tested the effects of genetic and non-genetic ITV on morphological traits in microcosms of protist communities by contrasting the effects of strains showing different ITV levels (i.e. trait averages and variance) on population growth, community composition and biomass production. We found that genetic and non-genetic ITV can lead to different effects on populations and communities across several generations. Furthermore, the effects of ITV declined across levels of biological organization: ITV directly altered population performance, with cascading but indirect consequences for community composition and biomass productivity. Overall, these results show that the drivers of ITV can have distinct effects on populations and communities, with cascading impacts on higher levels of biological organization that might mediate biodiversity-ecosystem functioning relationships.

Nutritional Quality of Basal Resource in Stream Food Webs Increased with Light Reduction-Implications for Riparian Revegetation.

Biofilms are considered a basal resource with high nutritional quality in stream food webs, as periphytic algae are abundant of polyunsaturated fatty acids (PUFAs). PUFAs are essential for growth and reproduction of consumers who cannot or have very limited capacity to biosynthesize. Yet, how the nutritional quality based on PUFA of basal food sources changes with light intensity remains unclear. We conducted a manipulative experiment in mesocosms to explore the response and mechanisms of nutritional quality to shading, simulating riparian restoration. We found a significant increase in PUFA% (including arachidonic acid, ARA) under shading conditions. The increased PUFA is caused by the algal community succession from Cyanobacteria and Chlorophyta to Bacillariophyta which is abundant of PUFA (especially eicosapentaenoic acid, EPA; docosahexaenoic acid, DHA). On the other hand, shading increased PUFA via upregulating enzymes such as Δ12 desaturase (FAD2, EC:1.14.19.6) and 3-ketoacyl-CoA synthase (KCS, EC:2.3.1.199) in the biosynthesis of unsaturated fatty acid elongation pathways. Our findings imply that riparian reforestation by decreasing light intensity increases the nutritional quality of basal resources in streams, which may enhance transfer of good quality carbon to consumers in higher trophic levels through bottom-up effects.

Multi-scale beta-diversity patterns in testate amoeba communities: species turnover and nestedness along a latitudinal gradient.

The relationship between species diversity and spatial scale is a central topic in spatial community ecology. Latitudinal gradient is among the core mechanisms driving biodiversity distribution on most scales. Patterns of ß-diversity along latitudinal gradient have been well studied for aboveground terrestrial and marine communities, whereas soil organisms remain poorly investigated in this regard. The West Siberian Plain is a good model to address diversity scale-dependence since the latitudinal gradient does not overlap with other possible factors such as elevational or maritime. Here, we collected 111 samples following hierarchical sampling (sub-zones, ecosystem types, microhabitat and replicate samples) and performed multi-scale partitioning of ß-diversity of testate amoeba assemblages as a model of study. We found that among-ecosystem ß-diversity is a leading scale in testate amoeba assemblages variation. Rare species determine ß-diversity at all scale levels especially in the northern regions, where rare taxa almost exclusively accounted for the diversity at the ecosystem level. ß-Diversity is generally dominated by the turnover component at all scales in lower latitudes, whereas nestedness prevailed at among-ecosystem scale in higher latitudes. These findings indicate that microbial assemblages in northern latitudes are spatially homogeneous and constrained by historical drivers at larger scales, whereas in southern regions, it is dominated by the turnover component both at the microhabitat and ecosystem scales and therefore determined by recent vegetation and environmental heterogeneity. Overall, we have provided the evidence for the existence of negative latitudinal gradient for among-ecosystem ß-diversity but not for among-microhabitat and among-sample ß-diversity for terrestrial testate amoeba communities.

Methylparaben changes the community composition, structure, and assembly processes of free-living bacteria, phytoplankton, and zooplankton.

Parabens are common contaminants in river and lake environments. However, few studies have been conducted to determine the effects of parabens on bacteria, phytoplankton, and zooplankton communities in aquatic environments. In this study, the effect of methylparaben (MP) on the diversity and community structure of the aquatic plankton microbiome was investigated by incubating a microcosm with MP at 0.1, 1, 10, and 100 µg/L for 7 days. The results of the Simpson index showed that MP treatment altered the α-diversity of free-living bacteria (FL), phytoplankton, and zooplankton but had no significant effect on the α-diversity of particle-attached bacteria (PA). Further, the relative abundances of the sensitive bacteria Chitinophaga and Vibrionimonas declined after MP addition. Moreover, the relative abundances of Desmodesmus sp. HSJ717 and Scenedesmus armatus, of the phylum Chlorophyta, were significantly lower in the MP treatment group than in the control group. In addition, the relative abundance of Stoeckeria sp. SSMS0806, of the Dinophyta phylum, was higher than that in the control group. MP addition also increased the relative abundance of Arthropoda but decreased the relative abundance of Rotifera and Ciliophora. The ß-diversity analysis showed that FL and phytoplankton communities were clustered separately after treatment with different MP concentrations. MP addition changed community assembly mechanisms in the microcosm, including increasing the stochastic processes for FL and the deterministic processes for PA and phytoplankton. Structural equation modeling analysis showed a significant negative relationship between bacteria richness and phytoplankton richness, and a significant positive relationship between phytoplankton (richness and community composition) and zooplankton. Overall, this study emphasizes that MP, at environmental concentrations, can change the diversity and structure of plankton microbial communities, which might have a negative effect on ecological systems.