<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2023.1166322</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Coordination and divergence in community assembly processes across co-occurring microbial groups separated by cell size</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Li</surname> <given-names>Xinghao</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/2304357/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Stegen</surname> <given-names>James C.</given-names></name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/182319/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Yu</surname> <given-names>Yuhe</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/210697/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Huang</surname> <given-names>Jie</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/1197685/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Hubei Key Laboratory of Regional Development and Environmental Response, Hubei Engineering Research Center for Rural Drinking Water Safety, Hubei University</institution>, <addr-line>Wuhan</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>Donghu Experimental Station of Lake Ecosystems, Key Laboratory of Aquatic Biodiversity and Conservation of Chinese Academy of Sciences, Institute of Hydrobiology, Chinese Academy of Sciences</institution>, <addr-line>Wuhan</addr-line>, <country>China</country></aff>
<aff id="aff3"><sup>3</sup><institution>Fundamental and Computational Sciences Directorate, Biological Sciences Division, Pacific Northwest National Laboratory</institution>, <addr-line>Richland, WA</addr-line>, <country>United States</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Xiaobo Liu, Nanjing University of Science and Technology, China</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Shouliang Huo, Chinese Research Academy of Environmental Sciences, China; Yuyi Yang, Chinese Academy of Sciences (CAS), China</p></fn>
<corresp id="c001">&#x002A;Correspondence: Jie Huang, <email>jhuang@ihb.ac.cn</email></corresp>
</author-notes>
<pub-date pub-type="epub">
<day>02</day>
<month>06</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>14</volume>
<elocation-id>1166322</elocation-id>
<history>
<date date-type="received">
<day>20</day>
<month>02</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>02</day>
<month>05</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2023 Li, Stegen, Yu and Huang.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Li, Stegen, Yu and Huang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>Setting the pace of life and constraining the role of members in food webs, body size can affect the structure and dynamics of communities across multiple scales of biological organization (e.g., from the individual to the ecosystem). However, its effects on shaping microbial communities, as well as underlying assembly processes, remain poorly known. Here, we analyzed microbial diversity in the largest urban lake in China and disentangled the ecological processes governing microbial eukaryotes and prokaryotes using 16S and 18S amplicon sequencing. We found that pico/nano-eukaryotes (0.22&#x2212;20 &#x03BC;m) and micro-eukaryotes (20&#x2212;200 &#x03BC;m) showed significant differences in terms of both community composition and assembly processes even though they were characterized by similar phylotype diversity. We also found scale dependencies whereby micro-eukaryotes were strongly governed by environmental selection at the local scale and dispersal limitation at the regional scale. Interestingly, it was the micro-eukaryotes, rather than the pico/nano-eukaryotes, that shared similar distribution and community assembly patterns with the prokaryotes. This indicated that assembly processes of eukaryotes may be coupled or decoupled from prokaryotes&#x2019; assembly processes based on eukaryote cell size. While the results support the important influence of cell size, there may be other factors leading to different levels of assembly process coupling across size classes. Additional studies are needed to quantitatively parse the influence of cell size versus other factors as drivers of coordinated and divergent community assembly processes across microbial groups. Regardless of the governing mechanisms, our results show that there are clear patterns in how assembly processes are coupled across sub-communities defined by cell size. These size-structured patterns could be used to help predict shifts in microbial food webs in response to future disturbance.</p>
</abstract>
<kwd-group>
<kwd>ecological process</kwd>
<kwd>cell size</kwd>
<kwd>high-throughput sequencing</kwd>
<kwd>prokaryotes</kwd>
<kwd>microbial eukaryotes</kwd>
</kwd-group>
<counts>
<fig-count count="5"/>
<table-count count="0"/>
<equation-count count="4"/>
<ref-count count="48"/>
<page-count count="10"/>
<word-count count="6570"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Systems Microbiology</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="S1" sec-type="intro">
<title>Introduction</title>
<p>Understanding the processes that structure microbial communities in ecosystems is a fundamental goal in ecology (<xref ref-type="bibr" rid="B19">Gleason, 1927</xref>; <xref ref-type="bibr" rid="B39">Stegen et al., 2012</xref>). Analytical frameworks have been developed to identify the ecological process (<xref ref-type="bibr" rid="B38">Stegen et al., 2013</xref>) and distribution patterns of microbial communities (<xref ref-type="bibr" rid="B10">Chen et al., 2020</xref>). Two types of processes, namely deterministic and stochastic, can influence microbial assembly (<xref ref-type="bibr" rid="B12">Dini-Andreote et al., 2015</xref>). Deterministic processes are generally associated with effects of abiotic conditions (environmental filtering) and biotic interactions (<xref ref-type="bibr" rid="B41">Vellend, 2010</xref>; <xref ref-type="bibr" rid="B48">Zhou and Ning, 2017</xref>). Stochastic processes include probabilistic dispersal and random changes in species abundances that are not the consequence of environmentally determined fitness (<xref ref-type="bibr" rid="B34">Rosindell et al., 2011</xref>).</p>
<p>Microbial community assembly was originally studied from a deterministic perspective (<xref ref-type="bibr" rid="B26">Lozupone and Knight, 2007</xref>; <xref ref-type="bibr" rid="B8">Chase and Myers, 2011</xref>), where empirical evidence shows that a variety of environmental factors such as temperature, pH, salinity, and organic carbon influence community establishment at different scales (<xref ref-type="bibr" rid="B3">Bell et al., 2009</xref>; <xref ref-type="bibr" rid="B18">Gilbert et al., 2012</xref>; <xref ref-type="bibr" rid="B43">Wagner et al., 2017</xref>). Recent studies, however, provide increasing support for a significant role of stochastic process in microbial systems (<xref ref-type="bibr" rid="B38">Stegen et al., 2013</xref>). These studies have resulted in the need to quantify the relative influence of deterministic and stochastic processes on community assembly (<xref ref-type="bibr" rid="B30">Ofiteru et al., 2010</xref>). However, most studies that examine assembly processes look at the whole community. The variation in assembly processes within subsets of communities has been rarely reported. To the best of our knowledge, sub-community analyses have primarily been performed on the basis of abundant and rare taxa, which has greatly improved our understanding of the ecological process governing microbial communities (<xref ref-type="bibr" rid="B29">Mo et al., 2018</xref>; <xref ref-type="bibr" rid="B46">Xue et al., 2018</xref>). Relative abundances are emergent properties of ecological assembly, however, and there is a need to further understand how assembly processes differ across intrinsic organismal attributes such as cell size.</p>
<p>Cell or body size is one of the most fundamental attributes of an organism and Elton proposed that the body sizes of organisms may strongly influence the ecological processes governing community structure (<xref ref-type="bibr" rid="B16">Elton, 1927</xref>). Biodiversity (e.g., the number of species) and species interactions (e.g., predator-prey feeding links within food webs) (<xref ref-type="bibr" rid="B16">Elton, 1927</xref>) are key attributes of ecological communities that can be influenced by organismal body sizes (<xref ref-type="bibr" rid="B4">Blanchard et al., 2017</xref>). More broadly, body size is considered a &#x201C;master trait&#x201D; that sets the pace of life, constrains the roles of organisms within food webs, and influences the structure and dynamics of various ecological networks (<xref ref-type="bibr" rid="B40">Trebilco et al., 2013</xref>; <xref ref-type="bibr" rid="B4">Blanchard et al., 2017</xref>). Measuring the relationships between cell size and microbial community assembly processes can, therefore, provide new insight into the fundamental principles governing microbial communities.</p>
<p>Here, we investigate community assembly across size fractionated samples from Lake Donghu, one of the largest freshwater lake complexes (approximately 33.9 km<sup>2</sup>) in China (<xref ref-type="bibr" rid="B47">Yan et al., 2017</xref>). Recently, a series of ecological improvements (e.g., the &#x201C;Great Lake Donghu Ecological Water Network&#x201D;) have been implemented to prevent the deterioration of water quality. This lake complex has experienced large inputs of nutrients and other pollutants from urban sewage and intensive aquaculture for nearly 30 years (<xref ref-type="bibr" rid="B21">Gong and Xie, 2001</xref>). Given the large spatial extent and large range of environmental conditions in Lake Donghu, we expected to observe significant influences of both dispersal limitation and environmental selection over microbial communities. Given that microbes interact with each other across size classes and taxonomic groups, we also expected assembly processes to be coordinated across the classes/groups studied here. Specific aims of this study were (i) to reveal the diversity and spatial distribution patterns of differently sized microbes across environmental gradients in a typical freshwater lake, and (ii) to analyze whether co-occurring groups of microbes with distinct cell sizes are assembled by distinct or coordinated ecological processes.</p>
</sec>
<sec id="S2" sec-type="materials|methods">
<title>Materials and methods</title>
<p>Microorganisms can be classified into different size ranges by filtration of water samples through membranes with different pore sizes (<xref ref-type="bibr" rid="B22">Hasegawa et al., 2003</xref>; <xref ref-type="bibr" rid="B27">Ma et al., 2017</xref>). The use of membrane filter techniques and DNA sequencing methods has improved the ability to study microbial diversity (<xref ref-type="bibr" rid="B17">Filker et al., 2016</xref>), especially in identifying rare, novel and atypical taxa (<xref ref-type="bibr" rid="B25">Liu et al., 2015</xref>).</p>
<sec id="S2.SS1">
<title>Size-based separation of microorganisms from water</title>
<p>Ten sites spanning an environmental gradient in the Lake Donghu complex were sampled in November 2016 according to the same protocol (<xref ref-type="fig" rid="F1">Figure 1</xref>). Briefly, 5 L of water at three depths (i.e., surface, 1 m, 1.5 m or deeper depending on the depths of the sampling sites) were collected and filtered through a nylon net with 200 &#x03BC;m pore size. Filtrate from these three depths was then pooled to represent the whole water column (<xref ref-type="bibr" rid="B17">Filker et al., 2016</xref>). The pooled water samples were processed immediately in the laboratory, whereby 1 L of water was serially filtered through 20 and 0.22 &#x03BC;m pore size Durapore membranes (Millipore) using a vacuum pump. One 20 &#x03BC;m sized filter and 3&#x2212;5 0.22 &#x03BC;m filters were obtained for each sample. Thus, two groups of microbes with different cell size ranges (i.e., micro-sized: 20&#x223C;200 &#x03BC;m and pico/nano-sized: 0.22&#x223C;20 &#x03BC;m, <xref ref-type="bibr" rid="B11">de Vargas et al., 2015</xref>) were obtained. For each site, three replicates were obtained and analyzed. After filtration, filters were immediately placed into cryovials and stored at &#x2212;80&#x00B0;C until DNA extraction.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption><p>The location of sampling sites in Lake Donghu complex. Color bar indicates the environmental gradient (eutrophication) from the most heavily polluted region to the less polluted region.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1166322-g001.tif"/>
</fig>
</sec>
<sec id="S2.SS2">
<title>Measurement of geographic distances and environmental variables</title>
<p>Geographic distance between each pair of sampling sites was estimated as the shortest distance over the Earth&#x2019;s surface and was determined based on the longitude/latitude coordinates of each site (<xref ref-type="supplementary-material" rid="TS1">Supplementary Table 1</xref>) using the &#x201C;distm&#x201D; function in &#x201C;geosphere&#x201D; R package (v3.6.1) (<xref ref-type="bibr" rid="B32">R Core Team, 2019</xref>). Water temperature (Temp), dissolved oxygen concentration (DO), electrical conductivity (CON), and pH of the pooled water samples were measured using a multi-parameter probe (WTW 3430). Chemical parameters, including dissolved organic carbon (DOC), nitrate nitrogen (NO<sub>3</sub>-N), nitrite nitrogen (NO<sub>2</sub>-N), ammonium nitrogen (NH<sub>4</sub>-N), and orthophosphate (PO<sub>4</sub>-P) were determined as previously described (<xref ref-type="bibr" rid="B36">Simon et al., 2015</xref>). Chlorophyll a (Chl <italic>a</italic>) was also quantified by spectrophotometry after filtration of the pooled samples through a GF/C membrane (Whatman, NJ, USA).</p>
</sec>
<sec id="S2.SS3">
<title>DNA extraction, PCR, and sequencing</title>
<p>Microbial genomic DNA was extracted using PowerWater DNA isolation kit (MoBio Laboratories, CA, USA) according to the manufacturer&#x2019;s instructions. For detecting micro-eukaryotes, the hypervariable V4 region of the SSU rDNA was amplified using the primer pair EK-565F (5&#x2032;-GCAGTTAAAAAGCTCGTAGT-3&#x2032;) and primer 18 s-EUK-1134-R-UNonMet (5&#x2032;-TTTAAGTTTCAGCCTTGCG-3&#x2032;) (<xref ref-type="bibr" rid="B6">Bower et al., 2004</xref>; <xref ref-type="bibr" rid="B36">Simon et al., 2015</xref>). The amplification conditions were as follows: 35 cycles (95&#x00B0;C for 30 s, 55&#x00B0;C for 30 s and 72&#x00B0;C for 45 s), preceded by 3 min of denaturation at 95&#x00B0;C, and ending with a 10 min final extension step at 72&#x00B0;C. For detecting prokaryotes, the V4 region of the 16S rRNA gene was amplified with the primer set 515f (5&#x2032;-GTGCCAGCMGCCGCGGTAA-3&#x2032;) and 806r (5&#x2032;-GGACTACHVGGGTWTCTAAT-3&#x2032;) (<xref ref-type="bibr" rid="B7">Caporaso et al., 2011</xref>). Triplicate purified PCR products were pooled for library preparation. Sequencing was performed on the Illumina MiSeq PE300 platform (Illumina, San Diego, USA) at Majorbio Bio-Pharm Technology Co. Ltd., (Shanghai, China).</p>
<p>Raw reads were assembled using FLASH (<xref ref-type="bibr" rid="B28">Magoc and Salzberg, 2011</xref>) and quality-trimmed with Trimmomatic (<xref ref-type="bibr" rid="B5">Bolger et al., 2014</xref>). UCHIME was used to identify and remove chimeric sequences (<xref ref-type="bibr" rid="B15">Edgar et al., 2011</xref>). Quality-checked sequences were then clustered into operational taxonomic units (OTUs; 97% cutoff) using UPARSE v7.1 (<xref ref-type="bibr" rid="B14">Edgar, 2013</xref>).<sup><xref ref-type="fn" rid="footnote1">1</xref></sup> Taxonomic classification of the OTUs was conducted by Blast alignments (<xref ref-type="bibr" rid="B1">Altschul et al., 1990</xref>) against the Silva database version 138<sup><xref ref-type="fn" rid="footnote2">2</xref></sup> and PR2 database for microbial prokaryotes and eukaryotes, respectively. Fungal and metazoan OTUs were excluded prior to statistical analyses. Each sample was rarefied to the same sequence depth (<italic>n</italic> = 32,029 sequences for prokaryotes; <italic>n</italic> = 16,313 sequences for micro-sized eukaryotes; <italic>n</italic> = 13,975 sequences for pico/nano-sized eukaryotes).</p>
</sec>
<sec id="S2.SS4">
<title>Description of the co-structure among subcommunities</title>
<p>Redundancy analysis (RDA) and canonical correspondence analysis (CCA) have been widely used to reveal variation in microbial community composition (<xref ref-type="bibr" rid="B13">Dray et al., 2003</xref>). These methods do not, however, directly couple the composition of co-occuring microbial groups (<xref ref-type="bibr" rid="B31">Palmer, 1993</xref>). To enable this coupling we used co-inertia analysis (COIA), which was specifically used to study relationships between each pair of subcommunity types (<xref ref-type="bibr" rid="B13">Dray et al., 2003</xref>).</p>
</sec>
<sec id="S2.SS5">
<title>Analytical framework for studying community assembly processes</title>
<sec id="S2.SS5.SSS1">
<title>Testing for phylogenetic signals</title>
<p>To quantify the relative influences of different community assembly processes, we applied the analytical framework developed by <xref ref-type="bibr" rid="B38">Stegen et al. (2013)</xref> with slight modifications. This framework estimates the influences of variable selection, homogeneous selection, homogenizing dispersal, dispersal limitation and an &#x201C;undominated&#x201D; fraction. This framework uses phylogenetic and compositional turnover across communities. This requires an &#x201C;OTU table&#x201D; based on community composition and a &#x201C;phylogenetic tree&#x201D; based on OTU phylogenetic relationships.</p>
<p>Using phylogenetic turnover to infer ecological processes requires the presence of phylogenetic signal, whereby more closely related taxa are more similar ecologically. To test for phylogenetic signal, we modified the procedures in <xref ref-type="bibr" rid="B38">Stegen et al. (2013)</xref> by using the between-OTU version of Bray-Curtis distances instead of the between-OTU niche differences to identify the microbial presence similarities. In this case, phylogenetic signal is indicated when the distribution patterns of closely related microbes are more similar to each other than to the distribution patterns of distant relatives. Nevertheless, it should be noted that lack of relationships using our method does not necessarily mean there is a lack of phylogenetic signal. This is because the method is more conservative than the standard approach described in <xref ref-type="bibr" rid="B38">Stegen et al. (2013)</xref>. Therefore, if no significant relationships can be found using our method, further testing using niche values is recommended. To test for phylogenetic signal, a linear regression with permutation to control for non-independence of the data (i.e., Mantel test) was performed using the between-OTU version of Bray-Curtis distances and the phylogenetic distances. This showed that phylogenetic signal existed across relatively short phylogenetic distances, consistent with previous work (<xref ref-type="fig" rid="F2">Figure 2</xref>). It is therefore most appropriate to quantify phylogenetic turnover among closest relatives (<xref ref-type="bibr" rid="B39">Stegen et al., 2012</xref>).</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption><p>Phylogenetic Mantel correlogram showing significant phylogenetic signals across short phylogenetic distances. Solid and open symbols denote significant and non-significant correlations, respectively, relating the between-OTU version of Bray-Curtis distances of <bold>(A)</bold> prokaryotes, <bold>(B)</bold> pico/nano-eukaryotes (0.22&#x2013;20 &#x03BC;m), and <bold>(C)</bold> micro-eukaryotes (20&#x2013;200 &#x03BC;m) to between-OTU phylogenetic distances across a given phylogenetic distance. Significantly positive correlations indicate that Bray-Curtis distance between OTUs increases with their phylogenetic distance, but only across the phylogenetic distance class being evaluated (that is, there is phylogenetic signal in OTU distributions).</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1166322-g002.tif"/>
</fig>
</sec>
<sec id="S2.SS5.SSS2">
<title>Quantifying community assembly processes</title>
<p>To characterize phylogenetic community composition at the local scale (for each individual community), we quantified the mean-nearest-taxon-distance (MNTD) and the nearest-taxon-index (NTI) using &#x201C;mntd&#x201D; and &#x201C;ses.mntd&#x201D; in the R package &#x201C;picante.&#x201D; NTI quantifies the number of standard deviations that the observed MNTD is from the mean of the null distribution (999 randomizations) (<xref ref-type="bibr" rid="B39">Stegen et al., 2012</xref>, <xref ref-type="bibr" rid="B38">2013</xref>). Note that NTI is the negative of &#x201C;mntd.obs.z&#x201D; contained in the output of &#x201C;ses.mntd&#x201D; (<xref ref-type="bibr" rid="B23">Kembel et al., 2010</xref>). MNTD and NTI are calculated as follow.</p>
<disp-formula id="S2.Ex1">
<mml:math id="M1">
<mml:mrow>
<mml:mpadded width="+3.3pt">
<mml:mtext>MNTD</mml:mtext>
</mml:mpadded>
<mml:mo rspace="5.8pt">=</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mo largeop="true" movablelimits="false" symmetric="true">&#x2211;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mtext>f</mml:mtext>
<mml:msub>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
</mml:msub>
<mml:mo>&#x2062;</mml:mo>
<mml:mtext>min</mml:mtext>
<mml:mo>&#x2062;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mtext>&#x0394;</mml:mtext>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
<mml:mo>&#x2062;</mml:mo>
<mml:msub>
<mml:mi mathvariant="normal">j</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where f<sub>i<sub>k</sub></sub> is the relative abundance of OTU <italic>i</italic> in community k, <italic>n</italic><sub><italic>k</italic></sub> is the number of OTUs in <italic>k</italic>, and min (&#x0394;<sub>i<sub>k</sub>j<sub>k</sub></sub>) is the minimum phylogenetic distance between OTU <italic>i</italic> and all other OTUs <italic>j</italic> that are also in <italic>k</italic>.</p>
<disp-formula id="S2.Ex2">
<mml:math id="M2">
<mml:mrow>
<mml:mpadded width="+3.3pt">
<mml:mtext>NTI</mml:mtext>
</mml:mpadded>
<mml:mo rspace="5.8pt">=</mml:mo>
<mml:mrow>
<mml:mo>-</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2062;</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi>MNTD</mml:mi>
<mml:mrow>
<mml:mtext>obs</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>-</mml:mo>
<mml:mover accent="true">
<mml:msub>
<mml:mtext>MNTD</mml:mtext>
<mml:mrow>
<mml:mtext>null</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>&#x00AF;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mtext>sd</mml:mtext>
</mml:mrow>
<mml:mo>&#x2062;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mtext>MNTD</mml:mtext>
<mml:mrow>
<mml:mtext>null</mml:mtext>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where MNTD<sub>obs</sub> denotes the observed mean of phylogenetic distance and MNTD<sub>null</sub> are null values of MNTD, and sd denotes the standard deviation. For a single community, NTI greater than +2 indicates coexisting taxa are more closely related than expected by chance (phylogenetic clustering). NTI less than &#x2212;2 indicates coexisting taxa are more distantly related than expected by chance (phylogenetic overdispersion). A mean NTI taken across all communities that is significantly different from zero indicates clustering (NTI&#x003E;0) or overdispersion (NTI&#x003C;0) on average.</p>
<p>To further quantify phylogenetic turnover among closely related OTUs, but between community pairs, we used the beta-nearest taxon index (&#x03B2;NTI), which quantifies the degree that the beta-mean-nearest taxon distance (&#x03B2;MNTD) deviates from a null model expectation, as described in <xref ref-type="bibr" rid="B38">Stegen et al.</xref> (<xref ref-type="bibr" rid="B38">2013</xref>, <xref ref-type="bibr" rid="B37">2015</xref>).</p>
<disp-formula id="S2.Ex3">
<mml:math id="M3">
<mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">&#x03B2;</mml:mi>
<mml:mpadded width="+3.3pt">
<mml:mtext>MNTD</mml:mtext>
</mml:mpadded>
</mml:mrow>
<mml:mo rspace="5.8pt">=</mml:mo>
<mml:mrow>
<mml:mn>0.5</mml:mn>
<mml:mo>&#x2062;</mml:mo>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mo largeop="true" movablelimits="false" symmetric="true">&#x2211;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mtext>f</mml:mtext>
<mml:msub>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
</mml:msub>
<mml:mo>&#x2062;</mml:mo>
<mml:mtext>min</mml:mtext>
<mml:mo>&#x2062;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mtext>&#x0394;</mml:mtext>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
<mml:mo>&#x2062;</mml:mo>
<mml:msub>
<mml:mi mathvariant="normal">j</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo rspace="5.8pt">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo rspace="5.8pt">+</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mo largeop="true" movablelimits="false" symmetric="true">&#x2211;</mml:mo>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:msub>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mtext>f</mml:mtext>
<mml:msub>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:msub>
</mml:msub>
<mml:mo>&#x2062;</mml:mo>
<mml:mtext>min</mml:mtext>
<mml:mo>&#x2062;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:msub>
<mml:mtext>&#x0394;</mml:mtext>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:msub>
<mml:mo>&#x2062;</mml:mo>
<mml:msub>
<mml:mi mathvariant="normal">j</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:msub>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where f<sub>i<sub>k</sub></sub> is the relative abundance of OTU <italic>i</italic> in community <italic>k</italic>, <italic>n</italic><sub><italic>k</italic></sub> is the number of OTUs in <italic>k</italic> and min (&#x0394;<sub>i<sub>m</sub>j<sub>k</sub></sub>) is the minimum phylogenetic distance between OTU <italic>i</italic> in community <italic>k</italic> and all OTUs <italic>j</italic> in community <italic>m</italic>. &#x03B2;MNTD was calculated using R function &#x201C;comdistnt&#x201D; (abundance. weighted = TRUE) (<xref ref-type="bibr" rid="B23">Kembel et al., 2010</xref>).</p>
<disp-formula id="S2.Ex4">
<mml:math id="M4">
<mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">&#x03B2;</mml:mi>
<mml:mtext>NTI</mml:mtext>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">&#x03B2;</mml:mi>
<mml:msub>
<mml:mtext>MNTD</mml:mtext>
<mml:mrow>
<mml:mtext>obs</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>-</mml:mo>
<mml:mover accent="true">
<mml:mrow>
<mml:mi mathvariant="normal">&#x03B2;</mml:mi>
<mml:msub>
<mml:mtext>MNTD</mml:mtext>
<mml:mrow>
<mml:mtext>null</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>&#x00AF;</mml:mo>
</mml:mover>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mo>/</mml:mo>
<mml:mtext>sd</mml:mtext>
</mml:mrow>
<mml:mo>&#x2062;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mi mathvariant="normal">&#x03B2;</mml:mi>
<mml:msub>
<mml:mtext>MNTD</mml:mtext>
<mml:mrow>
<mml:mtext>null</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where &#x03B2;MNTD<sub>obs</sub> is observed &#x03B2;MNTD, &#x03B2;MNTD<sub>null</sub> are null values of &#x03B2;MNTD, and sd indicates the standard deviation of the &#x03B2;MNTD<sub>null</sub> distribution. We quantified &#x03B2;NTI for all pairwise comparisons, using a separate null model for each comparison.</p>
<p>Raup-Crick was used to identify whether the observed degree of turnover deviates from that expected for pairwise comparisons of phylogenetic community turnover that were not significant (see below) based on &#x03B2;NTI. These non-significant &#x03B2;NTI values indicate the potential dominance of stochastic processes (<xref ref-type="bibr" rid="B9">Chase et al., 2011</xref>; <xref ref-type="bibr" rid="B37">Stegen et al., 2015</xref>). Some applications of Raup-Crick do not consider OTU relative abundances, which carry information useful for understanding ecological processes. Here, we performed Raup-Crick as <xref ref-type="bibr" rid="B37">Stegen et al. (2015)</xref> described using Bray-Curtis dissimilarity to quantify compositional turnover. In short, a null distribution of Bray-Curtis values was first generated. Then, deviations between empirically observed Bray-Curtis and the null distribution were standardized to range from &#x2212;1 to +1, which indicate whether local communities are more dissimilar (approaching 1), as dissimilar (approaching 0), or more similar (approaching &#x2212;1), than expected by random chance. The resulting metric is referred as RC<sub><italic>bray</italic></sub>.</p>
<p>Overall, for a given community we estimated the relative influence of variable selection or homogeneous selection as the fraction of its comparisons to other communities with &#x03B2;NTI&#x003E;+2 or &#x03B2;NTI&#x003C;-2, respectively. Selection is excluded as the dominant process when | &#x03B2;NTI| &#x003C;2; in these cases, RC<sub><italic>bray</italic></sub>&#x003E;+0.95 or&#x003C;-0.95 were taken as evidence that dispersal limitation or homogenizing dispersal, respectively, was the dominant process. For a given local community the relative influence of dispersal limitation was estimated as the fraction of its between-community comparisons with | &#x03B2;NTI| &#x003C;2 and RC<sub><italic>bray</italic></sub>&#x003E;+0.95. Similarly, the relative influence of homogenizing dispersal was estimated as the fraction of comparisons with | &#x03B2;NTI| &#x003C;2 and RC<sub><italic>bray</italic></sub>&#x003C;-0.95. The &#x201C;undominated&#x201D; scenario where | &#x03B2;NTI| &#x003C;2 and | RC<sub><italic>bray</italic></sub>| &#x003C;0.95 indicates that neither selection nor dispersal strongly drive compositional turnover and its relative contribution was estimated as the fraction of comparisons characterized by | &#x03B2;NTI| &#x003C;2 and | RC<sub><italic>bray</italic></sub> | &#x003C;0.95.</p>
</sec>
</sec>
</sec>
<sec id="S3" sec-type="results|discussion">
<title>Results and discussion</title>
<p>To reveal variations in ecological processes between prokaryotes and microbial eukaryotes, as well as between eukaryotes with different cell sizes, we compared quantitative estimates of community assembly process across these three groups of microorganisms. Our results show that the assembly processes are clearly distinct across the three co-occurring microbial groups that are operationally defined by cell size. This is consistent with an influence of organismal cell size on community structure (<xref ref-type="bibr" rid="B16">Elton, 1927</xref>). Given the dominant influence of cell size over nearly all aspects of organismal ecology, cell size is likely to have direct or indirect causal influences over the assembly process patterns observed here.</p>
<sec id="S3.SS1">
<title>Phylotype diversity varied significantly across sub-communities</title>
<p>Using high-throughput sequencing analysis, we offer a comprehensive evaluation of microbial diversity across the sampled lake complex. We obtained between 67,978 and 119,131 high-quality reads for each sampling site after Illumina sequencing and data cleaning (<xref ref-type="supplementary-material" rid="TS2">Supplementary Table 2</xref>). We recovered more than 85% of the total phylotype diversity (measured as OTU richness based on full data sets) with an average of 1,036 &#x00B1; 85 OTUs per sample (<xref ref-type="supplementary-material" rid="TS3">Supplementary Table 3</xref>). As expected, the prokaryotes (16 S) were generally more diverse (834 &#x00B1; 71 OTUs) than the eukaryotes (101 &#x00B1; 23 OTUs across all eukaryote samples/communities; One-way ANOVA, <italic>P</italic> &#x003C; 0.05) (<xref ref-type="fig" rid="F3">Figures 3A-C</xref>). Rarefaction analyses indicated that all samples were sampled to near saturation (<xref ref-type="fig" rid="F3">Figures 3D-F</xref>). Consistent with previous studies (<xref ref-type="bibr" rid="B47">Yan et al., 2017</xref>; <xref ref-type="bibr" rid="B45">Xiao et al., 2019</xref>), our results suggest that Proteobacteria (33.38%; 685 out of 2052 OTUs), Bacteroidetes (14.57%; 299/2052 OTUs) and Cyanobacteria (10.48%; 215/2052 OTUs) were the three most diverse prokaryotic phyla. Other phyla accounting for more than 5% of the total phylotype diversity include Actinobacteria (9.06%) and Chloroflexi (5.12%).</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption><p>Broad view of the microbial diversity in different microbial groups. The phylogenetic trees of <bold>(A)</bold> prokaryotes, <bold>(B)</bold> pico/nano-eukaryotes (0.22&#x2013;20 &#x03BC;m), and <bold>(C)</bold> micro-eukaryotes (20&#x2013;200 &#x03BC;m) were constructed on OTU-level and were colored according to phyla. Only the top 10 abundant phyla are shown in color <bold>(D&#x2013;F)</bold>. Rarefaction analyses predicted that all samples were sampled to near saturation. Prokaryotic diversity (estimated as number of OTUs) was much higher than eukaryotic diversity.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1166322-g003.tif"/>
</fig>
<p>Micro-eukaryotes (20&#x2212;200 &#x03BC;m) and pico/nano-eukaryotes (0.22&#x2212;20 &#x03BC;m) were characterized by similar phylotype diversity (i.e., 102 &#x00B1; 27 OTUs and 99 &#x00B1; 21 OTUs; One-way ANOVA, <italic>P</italic> &#x003E; &#x003E; 0.05), even though their community structures can be clearly distinguished from each other in the NMDS plot (phylogenetic composition, <xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 1</xref>). Nearly 50% of the phylotype diversity showed unique patterns in both groups of micro-eukaryotes (<xref ref-type="supplementary-material" rid="DS1">Supplementary Figures 2A, B</xref>), which was also revealed by morphological observations (<xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 3</xref>). Furthermore, the micro-eukaryotes were characterized by a high diversity of Ciliophora, Chlorophyta and Dinophyta, which accounted for up to 25.54%, 19.91%, and 9.52% of the total phylotype diversity (231 OTUs), respectively. The pico/nano-eukaryotes included Ciliophora, Chlorophyta, and Cercozoa, which accounted for up to 21.61%,16.53%, and 12.29% of the total phylotype diversity (236 OTUs), respectively.</p>
</sec>
<sec id="S3.SS2">
<title>The sub-communities co-varied with each other, potentially due to trophic interactions</title>
<p>Differences in eukaryotic community structures were visualized by NMDS ordination (Stress: 0.035), in which micro-eukaryotes and pico/nano-eukaryotes formed two distinct clusters (PERMANOVA, <italic>P</italic> &#x003C; 0.05, for details see <xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 1</xref> and <xref ref-type="supplementary-material" rid="TS4">Supplementary Table 4</xref>). As expected, we also successfully captured the known co-variation between eukaryotic and prokaryotic communities (<xref ref-type="fig" rid="F4">Figure 4</xref>; Monte-Carlo Test, <italic>P</italic> &#x003C; 0.05). This is potentially explained by their close relationships within aquatic microbial food webs (<xref ref-type="bibr" rid="B35">Sherr and Sherr, 2002</xref>), whereby prokaryotes generally act as the main food resource of micro-eukaryotes and thus strongly shape the microbial eukaryotic distribution via bottom-up effects (<xref ref-type="bibr" rid="B24">Kuparinen and Bjornsen, 1992</xref>). More importantly, we found that the co-variation differed across cell size classes as the micro-eukaryotes exhibited a closer relationship with the prokaryotes (<italic>R</italic><sup>2</sup> = 0.75) than the pico/nano-eukaryotes did (<italic>R</italic><sup>2</sup> = 0.58) (<xref ref-type="fig" rid="F4">Figures 4A, B</xref>). This result, along with studies on cladocera (<xref ref-type="bibr" rid="B20">Gliwicz, 1990</xref>; <xref ref-type="bibr" rid="B33">Rizo et al., 2019</xref>) and copepoda (<xref ref-type="bibr" rid="B42">Viitasalo et al., 1995</xref>), highlights the importance of cell size effects within aquatic microbial food webs. This pattern may be due to direct and indirect influences of cell size because cell size influences many aspects of organismal ecology (e.g., population abundance, growth rate, spatial niche) (<xref ref-type="bibr" rid="B44">Woodward et al., 2005</xref>). In terms of a direct influence of cell size, the most likely is via size-structured predator-prey feeding links (<xref ref-type="bibr" rid="B40">Trebilco et al., 2013</xref>). That is, close trophic interactions between microbial eukaryotes and prokaryotes likely only took place when the eukaryotes were big enough to be predators.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption><p>Description of the co-structure between the microbial groups. Co-inertia analysis (COIA) was applied to estimate the relationship <bold>(A)</bold> between prokaryotes and pico/nano-eukaryotes (0.22&#x2013;20 &#x03BC;m), <bold>(B)</bold> between prokaryotes and micro-eukaryotes (20&#x2013;200 &#x03BC;m), and <bold>(C)</bold> between pico/nano-eukaryotes and micro-eukarytoes. Each arrow of the COIA plots corresponds to a single sampling site, denoting the difference on the community structures between the compared microbial groups. Eigenvalues corresponding to the two co-inertia axes are shown in the histograms. Significance of co-structure was assessed by a Monte-Carlo test (RV coefficient, <italic>P</italic>).</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1166322-g004.tif"/>
</fig>
</sec>
<sec id="S3.SS3">
<title>The influence of deterministic assembly on local composition varied across sub-communities</title>
<p>NTI across all local communities was significantly different from zero (<italic>t</italic><sub>30</sub> = 7.0, <italic>P</italic> &#x003C; 0.0001), with a mean of +4.46 (<xref ref-type="fig" rid="F5">Figure 5</xref>). This result, coupled with testing for phylogenetic signal, revealed that, in the Lake Donghu complex, closely related and ecologically similar taxa coexist to a greater degree than expected by chance (<xref ref-type="fig" rid="F5">Figure 5A</xref>). Furthermore, the mean values of NTI varied significantly across different groups of microbes (One-way ANOVA, between the two eukaryotic groups: <italic>F</italic> = 4.6, <italic>P</italic> = 0.045; Overall: <italic>F</italic> = 162.5, <italic>P</italic> &#x003C; 0.0001). One of the key inferences is that the relative influences of deterministic and stochastic processes over community structure at the local scale was different among the three microbial sub-communities.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption><p>The relative influence of different ecological processes on microbial communities. <bold>(A)</bold> Summary of the analytical framework developed by <xref ref-type="bibr" rid="B39">Stegen et al. (2012</xref>, <xref ref-type="bibr" rid="B38">2013</xref>). Local-scale phylogenetic community composition was characterized using the nearest-taxon-index (NTI). Turnover in phylogenetic community composition and OTU composition were quantified using the beta nearest-taxon-index (&#x03B2;NTI) and Bray-Curtis based Raup Crick (RC<sub>bray</sub>), respectively. <bold>(B)</bold> Kernel density estimates for distributions of abundance weighted NTI. As indicated by the dashed lines, individual values below &#x2013;2 (negative) or above +2 (positive) are statistically significant, revealing phylogenetic overdispersion or clustering, respectively. Values of &#x03B2;NTI and RC<sub>bray</sub> are shown using heatmaps for <bold>(C)</bold> prokaryotes, <bold>(D)</bold> pico/nano-eukaryotes (0.22&#x2013;20 &#x03BC;m), and <bold>(E)</bold> micro-eukaryotes (20&#x2013;200 &#x03BC;m). Values are indicated by color bars. Insignificant values are indicated as gray. Significant values are indicated in red or blue. For RC<sub>bray</sub>, rectangles surrounded by a dashed line indicate that the corresponding &#x03B2;NTI value was significant.</p></caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fmicb-14-1166322-g005.tif"/>
</fig>
<p>Deterministic environmental filtering had a particularly strong influence over prokaryotic community composition, as indicated by an NTI distribution stretching from +4 to +14 with a median of +9.38 (<xref ref-type="fig" rid="F5">Figure 5B</xref>). These NTI values are high in an absolute sense and are also high relative to the other two sub-communities. This provides good evidence that membership within the prokaryotic communities in Lake Donghu was strongly limited by abiotic and/or biotic factors. We suggest that top-down forces imposed by eukaryotes are one such biotic factor that likely played an important role. The microbial eukaryotic communities were also largely shaped by environmental filtering, but NTI values were smaller relative to the prokaryotes and varied across the two cell-sized classes. The micro-sized eukaryotes experienced stronger influences of environmental filtering with 70% of local communities having significant NTI values, while only 30% of pico/nano-sized eukaryotic communities had significant NTI values. These results combined with the co-variation analyses indicate a highly deterministic sub-system based on interactions between prokaryotes and micro-eukaryotes, and with pico/nano-eukaryotes assembling more stochastically and independently of this deterministic sub-system.</p>
</sec>
<sec id="S3.SS4">
<title>Pico/nano-eukaryotes assembled independently from the other two sub-communities</title>
<p>To uncover the assembly processes influencing microbial community we coupled &#x03B2;NTI with RC<sub><italic>bray</italic></sub> (<xref ref-type="bibr" rid="B38">Stegen et al., 2013</xref>). We examined 270 pairwise comparisons (135 for &#x03B2;NTI and 135 for RC<sub><italic>bray</italic></sub>) to quantify microbial community turnover in the Lake Donghu complex (<xref ref-type="fig" rid="F5">Figures 5C-E</xref>). As previously shown, we found that environmental filtering has detectable influences on microbial assembly (<xref ref-type="bibr" rid="B2">Andersson et al., 2010</xref>; <xref ref-type="bibr" rid="B30">Ofiteru et al., 2010</xref>; <xref ref-type="bibr" rid="B12">Dini-Andreote et al., 2015</xref>). However, most of the &#x03B2;NTI values (122/135) were consistent with stochastic turnover (&#x2212;2&#x003C;&#x03B2;NTI&#x003C;+2). We then applied RC<sub><italic>bray</italic></sub> to pairwise comparisons with phylogenetic turnover that did not deviate from the null expectation (i.e., | &#x03B2;NTI| &#x003C;2). About 50% (i.e., 63/122) of these pairwise comparisons did not deviate from the null expectation (i.e., | RC<sub><italic>bray</italic></sub>| &#x003C;0.95). As such, the magnitude of community turnover in 47% of all pairwise comparisons (63/135) was not dominated by stochasticity or determinism.</p>
<p>These turnover results contrast with the within-community analyses that indicated strong influences of deterministic assembly. For example, only 7% (i.e., 3/45) of phylogenetic community turnover in the micro-sized eukaryotes was primarily due to environmental filtering (| &#x03B2;NTI| &#x003E;2). This fraction, however, is much smaller than the fraction of significant NTI values (70%). Similar trends were also observed for the prokaryotes whereby 22% of turnover was deterministic vs. 100% determinism at the within-community scale. This suggests that deterministic processes strongly influenced local community composition, while stochastic processes have a stronger influence on community turnover (<xref ref-type="bibr" rid="B29">Mo et al., 2018</xref>).</p>
<p>While nearly half of all pairwise comparisons were not dominated by a single type of process, we found that dispersal limitation was the primary driver for 51% of community turnover in micro-sized eukaryotes, and the &#x201C;undominated&#x201D; fraction had an important but secondary role (42%) (<xref ref-type="fig" rid="F5">Figure 5E</xref>). Similar results were observed for prokaryotes (<xref ref-type="fig" rid="F5">Figure 5C</xref>), which were also strongly governed by dispersal limitation (&#x223C;69%). These results, combined with the significant relationships between community dissimilarity and geographic distance (<xref ref-type="supplementary-material" rid="DS1">Supplementary Figure 4</xref>), suggest that larger spatial distance increase the influence of dispersal limitation in the sampled lake complex.</p>
<p>These findings contrast with the classic paradigm that &#x201C;all microbes are everywhere.&#x201D; Moreover, the dispersal limitation effect on structuring microbial communities may also be related to the features of the studied lake system, which has been divided into several lake regions by artificial dykes. Interestingly, we noted that over 90% of spatial turnover in the community composition of the pico/nano-eukaryotes was undominated. That is, neither deterministic selection or spatial processes had a strong or consistent influence in the assembly of pico/nano-eukaryotes. This can be due to many reasons such as weak selection, or selective processes that change frequently through space or time, and/or spatially variable rates of dispersal. The much larger undominated fraction for the pico/nano-eukaryotes, relative to the other two sub-communities, provides further evidence that community assembly of the pico/nano-eukaryotes operated independently from the other two sub-communities.</p>
</sec>
</sec>
<sec id="S4" sec-type="conclusion">
<title>Conclusion</title>
<p>This study showed that co-occurring groups of microbes that are operationally defined by cell size are governed by divergent assembly processes. However, we also observed multiple patterns indicating that the community assembly processes influencing micro-eukaryotes are linked to community assembly processes influencing prokaryotes. We propose that this link is due, in part, to size-structured trophic interactions that lead to within-community deterministic assembly. We also provide multiple patterns indicating that the community assembly of pico/nano-eukaryotes is largely independent of the sub-system generated by the coupling between micro-eukaryotes and prokaryotes. This points to a potential influence of cell size whereby larger cell size enables trophic interactions with prokaryotes, and these interactions link community assembly of micro-eukaryotes and prokaryotes. While the results are consistent with an influence of cell size, there are likely other factors at work that may be indirectly associated with cell size or linked to other size-independent attributes of organismal ecology. Thus, additional studies should be conducted to quantitatively parse the direct and indirect influences of cell size over assembly processes, as well as the influences of size-independent mechanisms. Moreover, the results in this work add to a growing literature showing an important influence of dispersal limitation, which is now broadly recognized as working alongside environmental filtering to shape microbial communities.</p>
</sec>
<sec id="S5" sec-type="data-availability">
<title>Data availability statement</title>
<p>Quality-checked assembled sequence data are deposited in NCBI Sequence Reads Archive under Bio-Project <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="PRJNA848972">PRJNA848972</ext-link>.</p>
</sec>
<sec id="S6" sec-type="author-contributions">
<title>Author contributions</title>
<p>JH and YY conceived and designed the study. XL analyzed the data and wrote the first version of the manuscript. JS and JH wrote and revised the manuscript. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec id="S7" sec-type="funding-information">
<title>Funding</title>
<p>This work was supported by the National Key Research and Development Program of China (No. 2018YFD0900701), the National Natural Science Foundation of China (Nos. 32030015 and 31801984), and the Youth Innovation Promotion Association of the Chinese Academy of Sciences (No. 2019333). JS was supported by an Early Career Research Program award (grant no. 74193) from the United States Department of Energy (US DOE) at the Pacific Northwest National Laboratory, a multiprogram national laboratory operated by Battelle for US DOE under contract DE-AC05-76RL01830.</p>
</sec>
<ack><p>We are grateful to Xueping Xiong and Wengen Zhu for sample collection, Changyu Zhu, Chengzhi Mao, and Yunfeng Wang for data analysis assistance, and Dr. Xiaotian Luo for microscopic observation.</p>
</ack>
<sec id="S8" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="S9" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="S10" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fmicb.2023.1166322/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fmicb.2023.1166322/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Table_1.XLSX" id="TS1" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_2.XLSX" id="TS2" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_3.XLSX" id="TS3" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_4.XLSX" id="TS4" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Data_Sheet_1.docx" id="DS1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<fn-group>
<fn id="footnote1">
<label>1</label>
<p><ext-link ext-link-type="uri" xlink:href="http://drive5.com/uparse/">http://drive5.com/uparse/</ext-link></p></fn>
<fn id="footnote2">
<label>2</label>
<p><ext-link ext-link-type="uri" xlink:href="http://www.arb-silva.de">http://www.arb-silva.de</ext-link></p></fn>
</fn-group>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Altschul</surname> <given-names>S. F.</given-names></name> <name><surname>Gish</surname> <given-names>W.</given-names></name> <name><surname>Miller</surname> <given-names>W.</given-names></name> <name><surname>Myers</surname> <given-names>E. W.</given-names></name> <name><surname>Lipman</surname> <given-names>D. J.</given-names></name></person-group> (<year>1990</year>). <article-title>Basic local alignment search tool.</article-title> <source><italic>J. Mol. Biol.</italic></source> <volume>215</volume> <fpage>403</fpage>&#x2013;<lpage>410</lpage>.</citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Andersson</surname> <given-names>A. F.</given-names></name> <name><surname>Riemann</surname> <given-names>L.</given-names></name> <name><surname>Bertilsson</surname> <given-names>S.</given-names></name></person-group> (<year>2010</year>). <article-title>Pyrosequencing reveals contrasting seasonal dynamics of taxa within baltic sea bacterioplankton communities.</article-title> <source><italic>ISME J.</italic></source> <volume>4</volume> <fpage>171</fpage>&#x2013;<lpage>181</lpage>. <pub-id pub-id-type="doi">10.1038/ismej.2009.108</pub-id> <pub-id pub-id-type="pmid">19829318</pub-id></citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bell</surname> <given-names>C. W.</given-names></name> <name><surname>Acosta-Martinez</surname> <given-names>V.</given-names></name> <name><surname>McIntyre</surname> <given-names>N. E.</given-names></name> <name><surname>Cox</surname> <given-names>S.</given-names></name> <name><surname>Tissue</surname> <given-names>D. T.</given-names></name> <name><surname>Zak</surname> <given-names>J. C.</given-names></name></person-group> (<year>2009</year>). <article-title>Linking microbial community structure and function to seasonal differences in soil moisture and temperature in a chihuahuan desert grassland.</article-title> <source><italic>Microb. Ecol.</italic></source> <volume>58</volume> <fpage>827</fpage>&#x2013;<lpage>842</lpage>. <pub-id pub-id-type="doi">10.1007/s00248-009-9529-5</pub-id> <pub-id pub-id-type="pmid">19466479</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Blanchard</surname> <given-names>J. L.</given-names></name> <name><surname>Heneghan</surname> <given-names>R. F.</given-names></name> <name><surname>Everett</surname> <given-names>J. D.</given-names></name> <name><surname>Trebilco</surname> <given-names>R.</given-names></name> <name><surname>Richardson</surname> <given-names>A. J.</given-names></name></person-group> (<year>2017</year>). <article-title>From bacteria to whales: using functional size spectra to model marine ecosystems.</article-title> <source><italic>Trends Ecol. Evol.</italic></source> <volume>32</volume> <fpage>174</fpage>&#x2013;<lpage>186</lpage>. <pub-id pub-id-type="doi">10.1016/j.tree.2016.12.003</pub-id> <pub-id pub-id-type="pmid">28109686</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bolger</surname> <given-names>A. M.</given-names></name> <name><surname>Lohse</surname> <given-names>M.</given-names></name> <name><surname>Usadel</surname> <given-names>B.</given-names></name></person-group> (<year>2014</year>). <article-title>Trimmomatic: a flexible trimmer for illumina sequence data.</article-title> <source><italic>Bioinformatics</italic></source> <volume>30</volume> <fpage>2114</fpage>&#x2013;<lpage>2120</lpage>.</citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bower</surname> <given-names>S. M.</given-names></name> <name><surname>Carnegie</surname> <given-names>R. B.</given-names></name> <name><surname>Goh</surname> <given-names>B.</given-names></name> <name><surname>Jones</surname> <given-names>S. R. M.</given-names></name> <name><surname>Lowe</surname> <given-names>G. J.</given-names></name> <name><surname>Mak</surname> <given-names>M. W. S.</given-names></name></person-group> (<year>2004</year>). <article-title>Preferential PCR amplification of parasitic protistan small subunit rDNA from metazoan tissues.</article-title> <source><italic>J. Eukaryot Microbiol.</italic></source> <volume>51</volume> <fpage>325</fpage>&#x2013;<lpage>332</lpage>. <pub-id pub-id-type="doi">10.1111/j.1550-7408.2004.tb00574.x</pub-id> <pub-id pub-id-type="pmid">15218702</pub-id></citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Caporaso</surname> <given-names>J. G.</given-names></name> <name><surname>Lauber</surname> <given-names>C. L.</given-names></name> <name><surname>Walters</surname> <given-names>W. A.</given-names></name> <name><surname>Berg-Lyons</surname> <given-names>D.</given-names></name> <name><surname>Lozupone</surname> <given-names>C. A.</given-names></name> <name><surname>Turnbaugh</surname> <given-names>P. J.</given-names></name><etal/></person-group> (<year>2011</year>). <article-title>Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>108</volume> <fpage>4516</fpage>&#x2013;<lpage>4522</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1000080107</pub-id> <pub-id pub-id-type="pmid">20534432</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chase</surname> <given-names>J. M.</given-names></name> <name><surname>Myers</surname> <given-names>J. A.</given-names></name></person-group> (<year>2011</year>). <article-title>Disentangling the importance of ecological niches from stochastic processes across scales.</article-title> <source><italic>Proc. R. Soc. B.</italic></source> <volume>366</volume> <fpage>2351</fpage>&#x2013;<lpage>2363</lpage>.</citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chase</surname> <given-names>J. M.</given-names></name> <name><surname>Kraft</surname> <given-names>N. J. B.</given-names></name> <name><surname>Smith</surname> <given-names>K. G.</given-names></name> <name><surname>Vellend</surname> <given-names>M.</given-names></name> <name><surname>Inouye</surname> <given-names>B. D.</given-names></name></person-group> (<year>2011</year>). <article-title>Using null models to disentangle variation in community dissimilarity from variation in alpha-diversity.</article-title> <source><italic>Ecosphere</italic></source> <volume>2</volume> <fpage>1</fpage>&#x2013;<lpage>11</lpage>.</citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>J.</given-names></name> <name><surname>Wang</surname> <given-names>P.</given-names></name> <name><surname>Wang</surname> <given-names>C.</given-names></name> <name><surname>Wang</surname> <given-names>X.</given-names></name> <name><surname>Miao</surname> <given-names>L.</given-names></name> <name><surname>Liu</surname> <given-names>S.</given-names></name><etal/></person-group> (<year>2020</year>). <article-title>Distinct assembly mechanisms underlie similar biogeographic patterns of rare and abundant bacterioplankton in cascade reservoirs of a large river.</article-title> <source><italic>Front. Microbiol.</italic></source> <volume>11</volume>:<issue>158</issue>. <pub-id pub-id-type="doi">10.3389/fmicb.2020.00158</pub-id> <pub-id pub-id-type="pmid">32117173</pub-id></citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>de Vargas</surname> <given-names>C.</given-names></name> <name><surname>Audic</surname> <given-names>S.</given-names></name> <name><surname>Henry</surname> <given-names>N.</given-names></name> <name><surname>Decelle</surname> <given-names>J.</given-names></name> <name><surname>Mahe</surname> <given-names>F.</given-names></name> <name><surname>Logares</surname> <given-names>R.</given-names></name><etal/></person-group> (<year>2015</year>). <article-title>Eukaryotic plankton diversity in the sunlit ocean.</article-title> <source><italic>Science</italic></source> <volume>348</volume>:<issue>1261605</issue>.</citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dini-Andreote</surname> <given-names>F.</given-names></name> <name><surname>Stegen</surname> <given-names>J. C.</given-names></name> <name><surname>van Elsas</surname> <given-names>J. D.</given-names></name> <name><surname>Salles</surname> <given-names>J. F.</given-names></name></person-group> (<year>2015</year>). <article-title>Disentangling mechanisms that mediate the balance between stochastic and deterministic processes in microbial succession.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>112</volume> <fpage>E1326</fpage>&#x2013;<lpage>E1332</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1414261112</pub-id> <pub-id pub-id-type="pmid">25733885</pub-id></citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dray</surname> <given-names>S.</given-names></name> <name><surname>Chessel</surname> <given-names>D.</given-names></name> <name><surname>Thioulouse</surname> <given-names>J.</given-names></name></person-group> (<year>2003</year>). <article-title>Co-inertia analysis and the linking of ecological data tables.</article-title> <source><italic>Ecology</italic></source> <volume>84</volume> <fpage>3078</fpage>&#x2013;<lpage>3089</lpage>.</citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Edgar</surname> <given-names>R. C.</given-names></name></person-group> (<year>2013</year>). <article-title>UPARSE: highly accurate OTU sequences from microbial amplicon reads.</article-title> <source><italic>Nat. Methods</italic></source> <volume>10</volume> <fpage>996</fpage>&#x2013;<lpage>998</lpage>. <pub-id pub-id-type="doi">10.1038/nmeth.2604</pub-id> <pub-id pub-id-type="pmid">23955772</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Edgar</surname> <given-names>R. C.</given-names></name> <name><surname>Haas</surname> <given-names>B. J.</given-names></name> <name><surname>Clemente</surname> <given-names>J. C.</given-names></name> <name><surname>Quince</surname> <given-names>C.</given-names></name> <name><surname>Knight</surname> <given-names>R.</given-names></name></person-group> (<year>2011</year>). <article-title>UCHIME improves sensitivity and speed of chimera detection.</article-title> <source><italic>Bioinformatics</italic></source> <volume>27</volume> <fpage>2194</fpage>&#x2013;<lpage>2200</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btr381</pub-id> <pub-id pub-id-type="pmid">21700674</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Elton</surname> <given-names>C. S.</given-names></name></person-group> (<year>1927</year>). <source><italic>Animal ecology.</italic></source> <publisher-loc>London</publisher-loc>: <publisher-name>Sidgwick and Jackson</publisher-name>.</citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Filker</surname> <given-names>S.</given-names></name> <name><surname>Sommaruga</surname> <given-names>R.</given-names></name> <name><surname>Vila</surname> <given-names>I.</given-names></name> <name><surname>Stoeck</surname> <given-names>T.</given-names></name></person-group> (<year>2016</year>). <article-title>Microbial eukaryote plankton communities of high-mountain lakes from three continents exhibit strong biogeographic patterns.</article-title> <source><italic>Mol. Ecol.</italic></source> <volume>25</volume> <fpage>2286</fpage>&#x2013;<lpage>2301</lpage>. <pub-id pub-id-type="doi">10.1111/mec.13633</pub-id> <pub-id pub-id-type="pmid">27029537</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gilbert</surname> <given-names>J. A.</given-names></name> <name><surname>Steele</surname> <given-names>J. A.</given-names></name> <name><surname>Caporaso</surname> <given-names>J. G.</given-names></name> <name><surname>Steinbruck</surname> <given-names>L.</given-names></name> <name><surname>Reeder</surname> <given-names>J.</given-names></name> <name><surname>Temperton</surname> <given-names>B.</given-names></name><etal/></person-group> (<year>2012</year>). <article-title>Defining seasonal marine microbial community dynamics.</article-title> <source><italic>ISME J.</italic></source> <volume>6</volume> <fpage>298</fpage>&#x2013;<lpage>308</lpage>.</citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gleason</surname> <given-names>H. A.</given-names></name></person-group> (<year>1927</year>). <article-title>Further views on the succession-concept.</article-title> <source><italic>Ecology</italic></source> <volume>8</volume> <fpage>299</fpage>&#x2013;<lpage>326</lpage>. <pub-id pub-id-type="doi">10.2307/1929332</pub-id></citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gliwicz</surname> <given-names>Z. M.</given-names></name></person-group> (<year>1990</year>). <article-title>Food thresholds and body size in cladocerans.</article-title> <source><italic>Nature</italic></source> <volume>343</volume> <fpage>638</fpage>&#x2013;<lpage>640</lpage>. <pub-id pub-id-type="doi">10.1038/343638a0</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gong</surname> <given-names>Z.</given-names></name> <name><surname>Xie</surname> <given-names>P.</given-names></name></person-group> (<year>2001</year>). <article-title>Impact of eutrophication on biodiversity of the macrozoobenthos community in a Chinese shallow lake.</article-title> <source><italic>J. Fresh. Ecol.</italic></source> <volume>16</volume> <fpage>171</fpage>&#x2013;<lpage>178</lpage>. <pub-id pub-id-type="doi">10.1080/02705060.2001.9663802</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hasegawa</surname> <given-names>H.</given-names></name> <name><surname>Naganuma</surname> <given-names>K.</given-names></name> <name><surname>Nakagawa</surname> <given-names>Y.</given-names></name> <name><surname>Matsuyama</surname> <given-names>T.</given-names></name></person-group> (<year>2003</year>). <article-title>Membrane filter (pore size, 0.22-0.45 &#x03BC;m; thickness, 150 &#x03BC;m) passing-through activity of <italic>Pseudomonas aeruginosa</italic> and other bacterial species with indigenous infiltration ability.</article-title> <source><italic>FEMS Microbiol. Lett.</italic></source> <volume>223</volume> <fpage>41</fpage>&#x2013;<lpage>46</lpage>. <pub-id pub-id-type="doi">10.1016/S0378-1097(03)00327-6</pub-id> <pub-id pub-id-type="pmid">12798998</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kembel</surname> <given-names>S. W.</given-names></name> <name><surname>Cowan</surname> <given-names>P. D.</given-names></name> <name><surname>Helmus</surname> <given-names>M. R.</given-names></name> <name><surname>Cornwell</surname> <given-names>W. K.</given-names></name> <name><surname>Morlon</surname> <given-names>H.</given-names></name> <name><surname>Ackerly</surname> <given-names>D. D.</given-names></name><etal/></person-group> (<year>2010</year>). <article-title>Picante: r tools for integrating phylogenies and ecology.</article-title> <source><italic>Bioinformatics</italic></source> <volume>26</volume> <fpage>1463</fpage>&#x2013;<lpage>1464</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btq166</pub-id> <pub-id pub-id-type="pmid">20395285</pub-id></citation></ref>
<ref id="B24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kuparinen</surname> <given-names>J.</given-names></name> <name><surname>Bjornsen</surname> <given-names>P. K.</given-names></name></person-group> (<year>1992</year>). <source><italic>Weddell sea ecology.</italic></source> <publisher-loc>Berlin</publisher-loc>: <publisher-name>Springer</publisher-name>.</citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>L.</given-names></name> <name><surname>Yang</surname> <given-names>J.</given-names></name> <name><surname>Yu</surname> <given-names>Z.</given-names></name> <name><surname>Wilkinson</surname> <given-names>D. M.</given-names></name></person-group> (<year>2015</year>). <article-title>The biogeography of abundant and rare bacterioplankton in the lakes and reservoirs of China.</article-title> <source><italic>ISME J.</italic></source> <volume>9</volume> <fpage>2068</fpage>&#x2013;<lpage>2077</lpage>. <pub-id pub-id-type="doi">10.1038/ismej.2015.29</pub-id> <pub-id pub-id-type="pmid">25748371</pub-id></citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lozupone</surname> <given-names>C. A.</given-names></name> <name><surname>Knight</surname> <given-names>R.</given-names></name></person-group> (<year>2007</year>). <article-title>Global patterns in bacterial diversity.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>104</volume> <fpage>11436</fpage>&#x2013;<lpage>11440</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.0611525104</pub-id> <pub-id pub-id-type="pmid">17592124</pub-id></citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ma</surname> <given-names>J.</given-names></name> <name><surname>Ping</surname> <given-names>D.</given-names></name> <name><surname>Dong</surname> <given-names>X.</given-names></name></person-group> (<year>2017</year>). <article-title>Recent developments of graphene oxide-based membranes: a review.</article-title> <source><italic>Membranes</italic></source> <volume>7</volume>:<issue>52</issue>. <pub-id pub-id-type="doi">10.3390/membranes7030052</pub-id> <pub-id pub-id-type="pmid">28895877</pub-id></citation></ref>
<ref id="B28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Magoc</surname> <given-names>T.</given-names></name> <name><surname>Salzberg</surname> <given-names>S. L.</given-names></name></person-group> (<year>2011</year>). <article-title>FLASH: fast length adjustment of short reads to improve genome assemblies.</article-title> <source><italic>Bioinformatics</italic></source> <volume>27</volume> <fpage>2957</fpage>&#x2013;<lpage>2963</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/btr507</pub-id> <pub-id pub-id-type="pmid">21903629</pub-id></citation></ref>
<ref id="B29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mo</surname> <given-names>Y.</given-names></name> <name><surname>Zhang</surname> <given-names>W.</given-names></name> <name><surname>Yang</surname> <given-names>J.</given-names></name> <name><surname>Lin</surname> <given-names>Y.</given-names></name> <name><surname>Yu</surname> <given-names>Z.</given-names></name> <name><surname>Lin</surname> <given-names>S.</given-names></name></person-group> (<year>2018</year>). <article-title>Biogeographic patterns of abundant and rare bacterioplankton in three subtropical bays resulting from selective and neutral processes.</article-title> <source><italic>ISME J.</italic></source> <volume>12</volume> <fpage>2198</fpage>&#x2013;<lpage>2210</lpage>. <pub-id pub-id-type="doi">10.1038/s41396-018-0153-6</pub-id> <pub-id pub-id-type="pmid">29880912</pub-id></citation></ref>
<ref id="B30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ofiteru</surname> <given-names>I. D.</given-names></name> <name><surname>Lunn</surname> <given-names>M.</given-names></name> <name><surname>Curtis</surname> <given-names>T. P.</given-names></name> <name><surname>Wells</surname> <given-names>G. F.</given-names></name> <name><surname>Criddle</surname> <given-names>C. S.</given-names></name> <name><surname>Francis</surname> <given-names>C. A.</given-names></name><etal/></person-group> (<year>2010</year>). <article-title>Combined niche and neutral effects in a microbial wastewater treatment community.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>107</volume> <fpage>15345</fpage>&#x2013;<lpage>15350</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1000604107</pub-id> <pub-id pub-id-type="pmid">20705897</pub-id></citation></ref>
<ref id="B31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Palmer</surname> <given-names>M. W.</given-names></name></person-group> (<year>1993</year>). <article-title>Putting things in even better order &#x2013; the advantages of canonical correspondence-analysis.</article-title> <source><italic>Ecology</italic></source> <volume>74</volume> <fpage>2215</fpage>&#x2013;<lpage>2230</lpage>.</citation></ref>
<ref id="B32"><citation citation-type="journal"><collab>R Core Team.</collab> (<year>2019</year>). <source><italic>R: a language and environment for statistical computing.</italic></source> <publisher-loc>Vienna</publisher-loc>: <publisher-name>R Foundation for Statistical Computing</publisher-name>.</citation></ref>
<ref id="B33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rizo</surname> <given-names>E. Z.</given-names></name> <name><surname>Xu</surname> <given-names>S. L.</given-names></name> <name><surname>Tang</surname> <given-names>Q. H.</given-names></name> <name><surname>Papa</surname> <given-names>R. D. S.</given-names></name> <name><surname>Dumont</surname> <given-names>H. J.</given-names></name> <name><surname>Qian</surname> <given-names>S. S.</given-names></name><etal/></person-group> (<year>2019</year>). <article-title>A global analysis of cladoceran body size and its variation linking to habitat, distribution and taxonomy.</article-title> <source><italic>Zool J. Linn Soc.</italic></source> <volume>187</volume> <fpage>1119</fpage>&#x2013;<lpage>1130</lpage>.</citation></ref>
<ref id="B34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rosindell</surname> <given-names>J.</given-names></name> <name><surname>Hubbell</surname> <given-names>S. P.</given-names></name> <name><surname>Etienne</surname> <given-names>R. S.</given-names></name></person-group> (<year>2011</year>). <article-title>The unified neutral theory of biodiversity and biogeography at age ten.</article-title> <source><italic>Trends Ecol. Evol.</italic></source> <volume>26</volume> <fpage>340</fpage>&#x2013;<lpage>348</lpage>. <pub-id pub-id-type="doi">10.1016/j.tree.2011.03.024</pub-id> <pub-id pub-id-type="pmid">21561679</pub-id></citation></ref>
<ref id="B35"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sherr</surname> <given-names>E. B.</given-names></name> <name><surname>Sherr</surname> <given-names>B. F.</given-names></name></person-group> (<year>2002</year>). <article-title>Significance of predation by protists in aquatic microbial food webs.</article-title> <source><italic>Antonie Van Leeuwenhoek</italic></source> <volume>81</volume> <fpage>293</fpage>&#x2013;<lpage>308</lpage>.</citation></ref>
<ref id="B36"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Simon</surname> <given-names>M.</given-names></name> <name><surname>Jardillier</surname> <given-names>L.</given-names></name> <name><surname>Deschamps</surname> <given-names>P.</given-names></name> <name><surname>Moreira</surname> <given-names>D.</given-names></name> <name><surname>Restoux</surname> <given-names>G.</given-names></name> <name><surname>Bertolino</surname> <given-names>P.</given-names></name><etal/></person-group> (<year>2015</year>). <article-title>Complex communities of small protists and unexpected occurrence of typical marine lineages in shallow freshwater systems.</article-title> <source><italic>Environ. Microbiol.</italic></source> <volume>17</volume> <fpage>3610</fpage>&#x2013;<lpage>3627</lpage>. <pub-id pub-id-type="doi">10.1111/1462-2920.12591</pub-id> <pub-id pub-id-type="pmid">25115943</pub-id></citation></ref>
<ref id="B37"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stegen</surname> <given-names>J. C.</given-names></name> <name><surname>Lin</surname> <given-names>X. J.</given-names></name> <name><surname>Fredrickson</surname> <given-names>J. K.</given-names></name> <name><surname>Konopka</surname> <given-names>A. E.</given-names></name></person-group> (<year>2015</year>). <article-title>Estimating and mapping ecological processes influencing microbial community assembly.</article-title> <source><italic>Front. Microbiol.</italic></source> <volume>6</volume>:<issue>370</issue>. <pub-id pub-id-type="doi">10.3389/fmicb.2015.00370</pub-id> <pub-id pub-id-type="pmid">25983725</pub-id></citation></ref>
<ref id="B38"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stegen</surname> <given-names>J. C.</given-names></name> <name><surname>Lin</surname> <given-names>X. J.</given-names></name> <name><surname>Fredrickson</surname> <given-names>J. K.</given-names></name> <name><surname>Chen</surname> <given-names>X. Y.</given-names></name> <name><surname>Kennedy</surname> <given-names>D. W.</given-names></name> <name><surname>Murray</surname> <given-names>C. J.</given-names></name><etal/></person-group> (<year>2013</year>). <article-title>Quantifying community assembly processes and identifying features that impose them.</article-title> <source><italic>ISME J.</italic></source> <volume>7</volume> <fpage>2069</fpage>&#x2013;<lpage>2079</lpage>.</citation></ref>
<ref id="B39"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Stegen</surname> <given-names>J. C.</given-names></name> <name><surname>Lin</surname> <given-names>X. J.</given-names></name> <name><surname>Konopka</surname> <given-names>A. E.</given-names></name> <name><surname>Fredrickson</surname> <given-names>J. K.</given-names></name></person-group> (<year>2012</year>). <article-title>Stochastic and deterministic assembly processes in subsurface microbial communities.</article-title> <source><italic>ISME J.</italic></source> <volume>6</volume> <fpage>1653</fpage>&#x2013;<lpage>1664</lpage>. <pub-id pub-id-type="doi">10.1038/ismej.2012.22</pub-id> <pub-id pub-id-type="pmid">22456445</pub-id></citation></ref>
<ref id="B40"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Trebilco</surname> <given-names>R.</given-names></name> <name><surname>Baum</surname> <given-names>J. K.</given-names></name> <name><surname>Salomon</surname> <given-names>A. K.</given-names></name> <name><surname>Dulvy</surname> <given-names>N. K.</given-names></name></person-group> (<year>2013</year>). <article-title>Ecosystem ecology: size-based constraints on the pyramids of life.</article-title> <source><italic>Trends Ecol. Evol.</italic></source> <volume>28</volume> <fpage>423</fpage>&#x2013;<lpage>431</lpage>. <pub-id pub-id-type="doi">10.1016/j.tree.2013.03.008</pub-id> <pub-id pub-id-type="pmid">23623003</pub-id></citation></ref>
<ref id="B41"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Vellend</surname> <given-names>M.</given-names></name></person-group> (<year>2010</year>). <article-title>Conceptual synthesis in community ecology.</article-title> <source><italic>Q. Rev. Biol.</italic></source> <volume>85</volume> <fpage>183</fpage>&#x2013;<lpage>206</lpage>. <pub-id pub-id-type="doi">10.1086/652373</pub-id> <pub-id pub-id-type="pmid">20565040</pub-id></citation></ref>
<ref id="B42"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Viitasalo</surname> <given-names>M.</given-names></name> <name><surname>Koski</surname> <given-names>M.</given-names></name> <name><surname>Pellikka</surname> <given-names>K.</given-names></name> <name><surname>Johansson</surname> <given-names>S.</given-names></name></person-group> (<year>1995</year>). <article-title>Seasonal and long-term variations in the body size of planktonic copepods in the northern baltic sea.</article-title> <source><italic>Mar. Biol.</italic></source> <volume>123</volume> <fpage>241</fpage>&#x2013;<lpage>250</lpage>. <pub-id pub-id-type="doi">10.1007/BF00353615</pub-id></citation></ref>
<ref id="B43"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wagner</surname> <given-names>R.</given-names></name> <name><surname>Zona</surname> <given-names>D.</given-names></name> <name><surname>Oechel</surname> <given-names>W.</given-names></name> <name><surname>Lipson</surname> <given-names>D.</given-names></name></person-group> (<year>2017</year>). <article-title>Microbial community structure and soil pH correspond to methane production in arctic alaska soils.</article-title> <source><italic>Environ. Microbiol.</italic></source> <volume>19</volume> <fpage>3398</fpage>&#x2013;<lpage>3410</lpage>. <pub-id pub-id-type="doi">10.1111/1462-2920.13854</pub-id> <pub-id pub-id-type="pmid">28677239</pub-id></citation></ref>
<ref id="B44"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Woodward</surname> <given-names>G.</given-names></name> <name><surname>Ebenman</surname> <given-names>B.</given-names></name> <name><surname>Emmerson</surname> <given-names>M.</given-names></name> <name><surname>Montoya</surname> <given-names>J. M.</given-names></name> <name><surname>Olesen</surname> <given-names>J. M.</given-names></name> <name><surname>Valido</surname> <given-names>A.</given-names></name><etal/></person-group> (<year>2005</year>). <article-title>Body size in ecological networks.</article-title> <source><italic>Trends Ecol. Evol.</italic></source> <volume>20</volume> <fpage>402</fpage>&#x2013;<lpage>409</lpage>. <pub-id pub-id-type="doi">10.1016/j.tree.2005.04.005</pub-id> <pub-id pub-id-type="pmid">16701403</pub-id></citation></ref>
<ref id="B45"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xiao</surname> <given-names>F.</given-names></name> <name><surname>Bi</surname> <given-names>Y.</given-names></name> <name><surname>Li</surname> <given-names>X.</given-names></name> <name><surname>Huang</surname> <given-names>J.</given-names></name> <name><surname>Yu</surname> <given-names>Y.</given-names></name> <name><surname>Xie</surname> <given-names>Z.</given-names></name><etal/></person-group> (<year>2019</year>). <article-title>The impact of anthropogenic disturbance on bacterioplankton communities during the construction of donghu tunnel (Wuhan, China).</article-title> <source><italic>Microb. Ecol.</italic></source> <volume>77</volume> <fpage>277</fpage>&#x2013;<lpage>287</lpage>. <pub-id pub-id-type="doi">10.1007/s00248-018-1222-0</pub-id> <pub-id pub-id-type="pmid">29951743</pub-id></citation></ref>
<ref id="B46"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Xue</surname> <given-names>Y.</given-names></name> <name><surname>Chen</surname> <given-names>H.</given-names></name> <name><surname>Yang</surname> <given-names>J. R.</given-names></name> <name><surname>Liu</surname> <given-names>M.</given-names></name> <name><surname>Huang</surname> <given-names>B.</given-names></name> <name><surname>Yang</surname> <given-names>J.</given-names></name></person-group> (<year>2018</year>). <article-title>Distinct patterns and processes of abundant and rare eukaryotic plankton communities following a reservoir cyanobacterial bloom.</article-title> <source><italic>ISME J.</italic></source> <volume>12</volume> <fpage>2263</fpage>&#x2013;<lpage>2277</lpage>. <pub-id pub-id-type="doi">10.1038/s41396-018-0159-0</pub-id> <pub-id pub-id-type="pmid">29899512</pub-id></citation></ref>
<ref id="B47"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yan</surname> <given-names>Q.</given-names></name> <name><surname>Stegen</surname> <given-names>J. C.</given-names></name> <name><surname>Yu</surname> <given-names>Y.</given-names></name> <name><surname>Deng</surname> <given-names>Y.</given-names></name> <name><surname>Li</surname> <given-names>X.</given-names></name> <name><surname>Wu</surname> <given-names>S.</given-names></name><etal/></person-group> (<year>2017</year>). <article-title>Nearly a decade-long repeatable seasonal diversity patterns of bacterioplankton communities in the eutrophic lake donghu (Wuhan, China).</article-title> <source><italic>Mol. Ecol.</italic></source> <volume>26</volume> <fpage>3839</fpage>&#x2013;<lpage>3850</lpage>. <pub-id pub-id-type="pmid">28437572</pub-id></citation></ref>
<ref id="B48"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhou</surname> <given-names>J.</given-names></name> <name><surname>Ning</surname> <given-names>D.</given-names></name></person-group> (<year>2017</year>). <article-title>Stochastic community assembly: does it matter in microbial ecology?</article-title> <source><italic>Microbiol. Mol. Biol. R</italic></source> <volume>81</volume>:<issue>e17</issue>. <pub-id pub-id-type="doi">10.1128/MMBR.00002-17</pub-id> <pub-id pub-id-type="pmid">29021219</pub-id></citation></ref>
</ref-list>
</back>
</article>
