Microscale ecology regulates particulate organic matter turnover in model marine microbial communities

The degradation of particulate organic matter in the ocean is a central process in the global carbon cycle, the mode and tempo of which is determined by the bacterial communities that assemble on particle surfaces. Here, we find that the capacity of communities to degrade particles is highly dependent on community composition using a collection of marine bacteria cultured from different stages of succession on chitin microparticles. Different particle degrading taxa display characteristic particle half-lives that differ by ~170 h, comparable to the residence time of particles in the ocean’s mixed layer. Particle half-lives are in general longer in multispecies communities, where the growth of obligate cross-feeders hinders the ability of degraders to colonize and consume particles in a dose dependent manner. Our results suggest that the microscale community ecology of bacteria on particle surfaces can impact the rates of carbon turnover in the ocean.

L earning how the composition of ecological communities impacts their function is one of the central challenges in ecology [1][2][3][4] . In the case of microbes, this problem is particularly complex, not only because of the extreme diversity of taxa and genes that make up microbial communities, but also because community function depends on microscale processes that are hard to measure such as aggregation, dispersal, and cell-cell interactions 5 . A prime example of the link between microscale community ecology and large-scale ecosystem function is found in the biological turnover of particulate organic matter. In the marine environment, biopolymer particles formed by aggregation of fragments of decaying organisms, fecal pellets, and extracellular polysaccharides are degraded and consumed by heterotrophic bacteria that attach to particle surfaces and form dense microbial communities of large taxonomic and metabolic diversity [6][7][8][9] . Because particulate matter tends to sink in the water column, its degradation in the upper layers of the ocean where oxygen abounds is crucial to sustain the marine food web and prevent the sequestration of carbon and nitrogen into the deep sea 9-11 . Therefore, particle-attached microbial communities play a fundamental role by closing the loop of the global carbon cycle and maintaining the balance of nutrients in marine ecosystems. Although many physical aspects of the bacteria-particle interaction such as attachment or the effects of flow 12,13 have been well characterized, the potential role that ecological interactions between microbes may play in controlling the dynamics of particle colonization and degradation-and thus the "mode and tempo" of the global carbon cycle-is much less clear.
Previous studies have suggested that ecological interactions between microbes can play a significant role in controlling the dynamics of community assembly on particles 14 . Competition for particle surface and thus primary resource access is likely to be strong among particle-attached bacteria and interference competition mediated by secondary metabolites can be a powerful strategy to deter competitors 15,16 . Moreover, over the time scales of particle turnover, trophic interactions mediated by byproducts of degradation and primary metabolism can strongly influence the overall dynamics of bacterial growth 17 . To release the carbon trapped in particulate matter, bacteria secrete hydrolytic enzymes -often regulated by quorum sensing signals 18 -that deconstruct complex biopolymers and release soluble sugars into the environment. The bioavailable sugars can in turn be taken up by nearby cells, thus unlocking a niche for cheaters that consume resources but do not contribute to degradation 17,19 . Likewise, byproducts of primary metabolism such as organic acids or amino acids that are released to the local environment can be consumed by cross-feeding bacteria that co-assemble on the particle. On chitin particles, these types of trophic interaction have been shown to lead to successional waves and invasion of secondary consumers, which eventually become the numerically dominant members of the community 17 . These findings led us to hypothesize that interactions across trophic levels at the microscale might alter the catabolism of chitin and consumption of byproducts, possibly affecting the rate of particle turnover and the conversion from particle to bacterial biomass.
To test this hypothesis, in this study we used an isolate collection obtained directly from particle-attached communities previously shown to colonize in microscale successions 17 . In brief, these communities were enriched on~50 µm paramagnetic chitin hydrogel particles incubated in seawater from the coastal ocean (Nahant, MA, USA). Bacteria were isolated directly from the particles, resulting in a collection that includes taxonomic orders such as Alteromonadales, Flavobacteriales, Rhodobacterales, Vibrionales, and Oceanospirillales. Notably, the composition of our collection coincides well with the taxonomic profiles of natural chitinous marine particles collected at 200-500 m depth in the North Pacific gyre 20 . This overlap between our isolate collection and the taxonomic composition of natural particleattached communities suggests that isolates obtained from model particles represent a relevant set of strains with which to study the effect of ecological interactions on particle turnover.
Bacterial isolates in our collection fall into two coarse-grained functional groups, defined on the basis of shared physiological characteristics and colonization dynamics 17 . The first group comprises primary degraders, which secrete chitinolytic enzymes, are motile, can grow rapidly on degradation byproducts, and belong to species that tend to appear early during particle colonization. The second group corresponds to secondary consumers, which in general do not secrete enzymes, cannot grow on chitin, grow poorly if at all on monomers, are not motile, and belong both to late and early successional species (Fig. 1a, Supplementary  Fig. 1). Although secondary consumers cannot grow on chitin particles alone, they can reach 100-1000 fold higher abundance in the presence of primary degraders 17 due to their ability to utilize carbon sources released by primary degraders during colonization.
Our goal in this study is to provide a quantitative description how particle degradation kinetics depend on the assembly of primary degraders and secondary consumers during particle colonization. To this end, we first studied how monocultures of primary degraders consumed particles by tracking changes in particle volume over time using high-throughput, high-resolution time-lapse microscopy ( Fig. 1b, Supplementary Movie 1), and guiding our analysis with simple mathematical models of colonization and resource consumption. Subsequently, we assembled two-strain communities of primary degraders and secondary consumers and developed a quantitative phenomenological characterization of the impact of secondary consumers on degradation. Our results reveal that early colonizing taxa can differ significantly in their hydrolytic power to breakdown chitin, that particle degradation is limited by the number of enzymesecreting bacteria that colonize the particle surface, and that secondary consumers effectively become parasites that increase in abundance at the cost of the primary degraders when cocolonizing on particle surfaces. Furthermore, the presence of parasitic secondary consumers can delay or even obstruct particle degradation. All these effects suggest that microscale community ecology on particle surfaces plays a major role in controlling community function by primarily slowing down resource turnover rates.

Results
Primary degraders display differences in hydrolytic power. We tracked the dynamics of particle consumption by measuring changes in particle volume over time, V(t), using high-throughput time-lapse microscopy of individual chitin microbeads. We chose an initial concentration of degrader cells of 5 × 10 5 cells/ml-an upper-bound estimate of the concentration of degrading bacteria in coastal waters 21 -and quantified V(t) over a period of 240 h, for four primary degraders and four secondary consumers incubated in media with no carbon source other than the particle. As expected, secondary consumers, which are unable to use chitin as a carbon source, did not grow on particles in monoculture and therefore did not affect V(t) over the course of the 10-day timelapse. For primary degraders, instead, V(t) was characterized by a long period of no detectable change, followed by an abrupt collapse (Fig. 1c, Supplementary Movie 2). Measurements of bacterial growth during degradation showed that bacteria grew steadily from the beginning of the incubation and that, despite no apparent change in particle volume, depolymerization was a continuous process that culminated in particle collapse (Fig. 1c, Supplementary Fig. 2). In some cases particles seemed to increase in volume immediately before collapse, possibly indicating that the degradation of cross-linked chitin in the hydrogel allowed water molecules to expand the matrix 22 . We interpreted the swelling-collapse transition as the point where a certain critical amount of polymer was degraded and used the time to reach this transition to quantify the ability of bacteria to consume particles. More precisely, we measured the hydrolytic power of degraders in terms of the particle half-life, τ 1/2 , i.e., the time it took for the particle to decrease to half its volume (see Methods).
We found a remarkable variation in τ 1/2 among the four different primary degraders, despite the fact that all of these isolates appeared early on in the ecological succession on chitin particles (Fig. 1a). At an initial cell concentration of 5 × 10 5 cells/ ml for all primary degraders, particle half-lives varied from~30 h for the fastest degrader (a strain of the genus Psychromonas, named psych6C06) to~200 h for the slow degraders (a strain of Vibrio nigripulchritudo named vnigr6D03) (Fig. 1d). The large number of chitinase copies in psych6C06 (19 copies) suggested that gene dosage played a role in controlling the hydrolytic power of the strains. However, overall the differences between τ 1/2 among primary degraders could not be clearly correlated to variation in gene content, suggesting instead that expression levels and the properties of extracellular enzymes played a more significant role. Gene content did, however, distinguish primary and secondary consumers: GlcNAc-specific phosphotransferase (PTS) transporters, chitin binding proteins, and GlcNAc-specific chemotaxis genes where all enriched among degraders compared to secondary consumers (Fig. 1e, Supplementary Table 1).
Chitin degradation is intrinsically linked to the production of public goods such as chitinases and as such can be subject to cooperative growth dynamics 23 , i.e., a positive dependency between cell densities and growth or depolymerization rates. If cooperativity does play a role, half degradation times would be highly sensitive to cell numbers, increasing disproportionally in cases where cell load is low. To test the relevance of this phenomenon and, in general, to study how τ 1/2 depended on initial conditions, we measured degradation kinetics as a function of the initial concentration of primary degrader, [B p ] 0 , which until now was arbitrarily set to 5 × 10 5 cells/ml. In addition, we guided our analysis with simple models of particle degradation and bacterial growth. To construct these models, we assumed that particle depolymerization was proportional to the density of bacteria, B. We studied two possibilities: that bacteria displayed cooperative growth kinetics, i.e., with growth rate proportional to B n with n > 1 and that cooperativity played no significant role and growth and occurred at fixed, density independent percapita rates. Assuming that τ 1/2 depends linearly on the speed of depolymerization, under a scenario of no cooperativity we predict that τ 1/2 should scale linearly with −ln[B p ] 0 , whereas cooperative growth would make τ 1/2 scale nonlinearly with ln [B p ] 0 (see Fig. 2f and Methods).
In agreement with the simplest model with no cooperativity, we find a linear relation between τ 1/2 and ln[B p ] 0 (Fig. 2, Supplementary Table 2). This behavior implies that the particle half-life is controlled by simple mass action kinetics 24 that-at least in the conditions of our experiment-are not influenced by cooperativity. More precisely, we find that τ 1/2 is well described by the following expression: where t 0 is the intercept of the lines in Fig. 2e and represents a timescale to degradation that is intrinsic to each strain, β is the slope and represents the per-capita contribution to the degradation process and ln[B p ] 0 captures the effect of the primary  Table 1 degrader concentration in the local environment, akin to a chemical potential for the cell-particle reaction. The relationship found in Eq. (1) shows that the turnover of particulate organic matter depends on the load of primary degraders in the milieu in a simple, predictable manner. The lack of an observed cooperativity suggests that the possible benefits that bacteria may derive from sharing hydrolysis products are effectively offset by local competition for resources between neighbors. Overall, our results indicate that variation in the composition and abundance of primary degraders can have a significant impact on the rate of particulate organic matter turnover.
Secondary consumers behave as parasites. To understand how ecological interactions between primary degraders and secondary consumers influence particle degradation, we focused our analysis on two primary degraders and one secondary consumer. We chose the relatively "slow" degrader, vsple1A01 (Fig. 1c, d) a member of the Vibrio splendidus clade, the most abundant group of marine vibrios in coastal seawaters 25 , and the relatively "fast" degrader, Pseudoalteromonas sp. palte3D05, a common member of heterotrophic bacterioplankton communities 26,27 . Secondary consumers, or strains unable to degrade chitin, have previously been found to invade particle-attached communities and to become numerically dominant during community assembly 17 (Fig. 1a). We focused our efforts on a secondary consumer cultivated from seawater-incubated chitin particles, a strain of the genus Maribacter (a type of marine Flavobacteria), that we here call marib6B07. As with other secondary consumers, marib6B07 is able to cross-feed when grown in coculture with degraders 17 . Interestingly, despite their inability to degrade chitin under laboratory conditions, the genomes of marib6B07 and the other secondary consumers all contain chitinases (e.g., marib6B07 has two). Yet, in contrast to degraders, none of these genomes contain genes for PTS transport of chitin monomers and only one out of the four consumers have genes for GlcNAc-specific chemotaxis and chitin-binding proteins (Supplementary Table 1). These differences in the genomes of primary degraders and secondary consumers suggest that their functional roles in the community may be determined by the interaction between multiple traits, such as the ability to swim toward breakdown products of chitin and to transport them into the periplasm.
Co-incubation of mari6B07 with vsple1A01 and palte3D05 showed that mari6B07 increased τ 1/2 relative to primary degrader monocultures (Fig. 3a), implying that the cross-feeder impaired the ability of degrader populations to depolymerize the particle. To study this phenomenon in a quantitative manner, we measured how τ 1/2 responded to changes in the initial concentration of secondary consumer, [B s ] 0 , with the number of cells of the primary degrader fixed at a given concentration ([B p ] 0 ≈ 1.25 × 10 5 cells/ml) (Fig. 3a). We found that for small values of [B s ] 0 , τ 1/2 increased roughly linearly, such that a onefold increase in the secondary consumer [B s ] 0 had approximately the same effect as a tenfold reduction of the primary degrader [B p ] 0 in monoculture.
Surprisingly, at a threshold [B s ] 0 we observed an abrupt increase in τ 1/2 , to the extent that particle degradation did not occur within the 240 h imaging period, suggesting that the population of primary producers might have been inhibited from colonization and/or growth. To investigate how this phenomenon depended on the composition of the two-strain community, we varied the abundance of the primary degraders, [B p ] 0 and secondary consumer [B s ] 0 , in order to obtain degradation phase  (Fig. 3b). The degradation phase planes show that complete inhibition did not depend on the total concentration of the secondary consumer, [B s ] 0 , but on the ratio of secondary consumer to primary degrader γ = [B s ] 0 /[B p ] 0 , (Fig. 3b, c). For the slow degrader, vsple1A01, degradation was blocked at γ >~1, whereas for the fast degrader, palte3D05, degradation was blocked above a ratio of γ >~16, showing that the slow degrader was more sensitive to the inhibitory effects of secondary consumer marib6B07 than the fast degrader. This analysis indicates that the balance between the relative abundances of secondary consumers to primary degraders in the environment, in addition to the degradation kinetics of the primary consumer, may be an important parameter that dictates the turnover rates of carbon over short time scales (see Discussion). Quantification of the abundance of each strain in coculture before and after particle degradation showed that the interaction between primary degrader and secondary consumer is parasitic, i.e., positive for the secondary consumer, negative for the degrader. CFU counts during the time course of degradation in cocultures of vsple1A01 and marib6B07 showed that primary degrader growth rate and yield were lower than in monoculture, and that the loss of degrader cells was compensated by the growth of secondary consumers (Fig. 3c). Secondary consumers doubled approximately five times by the time of particle collapse, in contrast to their zero doublings in monoculture. Notably, the total yield of the coculture was always equal or lower to the yield of the monoculture, highlighting the parasitic nature of the interaction. Similar results were observed with cocultures of palte3D05 and marib6B07 ( Supplementary Fig. 3). Thus, secondary consumers, whose growth is facilitated by primary degraders, exert a negative feedback on degraders, limiting their ability to hydrolyze resources, and thereby potentially limiting their own growth.
Given the higher ratio of secondary consumer to degrader (γ) required to inhibit palte3D05 compared to vsple1A01, we hypothesized that "slow degraders" might be more susceptible to the detrimental effect of secondary consumers. To test this hypothesis as well as whether the observed parasitic interactions can be generalized to other primary degrader-secondary consumer pairs, we measured the effect of coculture at γ = 1 on particle degradation for all primary degraders (Fig. 1d) with four different secondary consumers (including marib6B07) of diverse taxonomic origins, all of which were co-isolated from the same chitin-attached communities. The results showed that while the fast degraders psych6C06 and palte3D05 were only mildly affected by coculturing with secondary consumers at γ = 1, the slow degraders vsple1A01 and vnigr6D03 were highly susceptible (Fig. 4). The slowest degrader, vnigr6D03, was in fact inhibited by all four secondary consumers, three of which caused total blockage of particle consumption. These data further indicate that parasitic interactions between degraders and consumers may not depend on specific taxa, but rather on the hydrolytic power of the degrader.
Consistent with the observation that interactions are not specific to strains or species but to functional roles (i.e., secondary consumer, primary degrader), we did not find evidence of Agar plate assays designed to detect secreted inhibitory factors showed no interaction between the secondary consumer and primary degraders. Moreover, cocultures of vsple1A01 and palte3D05 with marib6B07 in liquid media supplemented with GlcNAc (the monomer of chitin) as sole carbon source showed no decrease in growth rates (Supplementary Fig. 4). This suggests that either an antagonistic factor is only secreted in the particle environment or, what is more likely, that the observed inhibition of primary degrader growth is based on interference with physical processes that only take place when resources are concentrated on particles (e.g., colonization and attachment).

Discussion
Despite the significant efforts put into understanding the factors that drive the turnover of organic matter in the ocean 28,29 , the potential role that microbial interactions may play in this process has remained relatively unexplored. Our study leveraged a simplified model based on wild isolates that naturally colonize chitin particles to dissect this question. We provided evidence that both differences in primary degrader type and the ratio of primary degrader to secondary consumer can significantly alter particle degradation kinetics. Remarkably, we show that even in the ideal conditions of our experiments (no N limitation and high number of cells pregrown in rich media) particle turnover times can be as high as 200 h or more, that is, in the same range as the residence time of particles in the ocean's mixed layer. Moreover, we showed interactions between primary degraders and secondary consumers lead to a significant increase in particle turnover times. We observed this phenomenon with different secondary consumers, suggesting that strain-independent processes such as blocking of the particle surface or alteration of resource gradients could mediate inhibition. The slowdown of the particle degradation kinetics imposed by the growth of secondary consumers is consistent with culture-independent data on colonization dynamics in natural seawater, which showed that secondary consumers invade particle-attached communities in a successional manner, displacing primary degraders and becoming numerically dominant in the community 17 . Taken together these results suggest that the microscale community ecology of particle-attached bacteria can play an important role in controlling rates of carbon turnover in the ocean. Further work will be needed to assess the precise mechanisms of interaction between primary degraders and secondary consumers on particle surfaces. Models with increased realism, including more complex particle compositions and predation by grazers or phages 14,30,31 will be valuable to build a more complete picture of the community ecology processes that control the carbon cycle in the ocean.

Methods
Bacterial culturing conditions. Bacterial strains used in this study were previously isolated from model chitin particles 17  Particle degradation experiments. Particle degradation experiments were performed in clear 96-well plates (VWR 10062-900). Each well contained 180 µl Tibbles-Rawling minimal media, bacterial cells at defined concentrations prepared as described above, and approximately 100 chitin magnetic beads (New England Biolabs #E8036L). Before being used in the experiments, the chitin magnetic beads storage buffer was removed using a neodymium magnet (McMaster-Carr #5862K38) to retain the beads. Beads were washed twice in Tibbles-Rawling minimal media and size selected using 100 µm and 40 µm strainers (VWR, #10199-658 and #10199-654, respectively).
For Fig. 1b, the colonized particle was stained in the well after 24 by adding Syto9 (Thermo Fisher, S34854), 500 nM final concentration for 1 h at room temperature in the dark. Microscopy was performed on an EVOS FL Auto Imaging System (Fisher #AMAFD1000) using a GFP lightcube (Thermo Fisher AMEP4651) and a 20× fluorite, long working distance objective (Fisher #AMEP4682, NA 0.40, WD 3.1 mm) and the software's (revision 31201) Z-stack function. 3Dreconstruction was done using the ImageJ distribution Fiji (ImageJ 1.51N).
Time-lapse imaging. Phase-contrast time-lapse images were acquired with an EVOS FL Auto Imaging System (Fisher #AMAFD1000) using the EVOS software (revision 31201) and a 20× fluorite, long working distance, phase-contrast objective (Fisher #AMEP4682, NA 0.40, WD 3.1 mm). Images were manually focused for each particle to capture the maximum cross-section area (see Fig. 1c, upper panel). Time lapses ran a maximum of 240 h, with images acquired every 2 h. To minimize evaporation effects, culturing plates were wrapped in para film during the timelapse experiments and outer wells filled with 200 µl water.
Image processing and volume quantification. Phase-contrast images were analyzed using the ImageJ distribution Fiji (ImageJ 1.51N). A polygonal shape was manually drawn around the particle to determine the area of the particles' crosssection. To convert from cross-section area in square pixel (1 pixel = 0.4545 µm) to volume (in µm 3 ), we assumed a spherical shape of the particles. Volumes were  Table 3 for raw data normalized to initial volume at t = 0 h to account for variation in particle sizes. In order to estimate the particle half-life, we fitted a sigmoidal function 1 1þe k xÀτ 1=2 ð Þ ð Þ using MATLAB (Version R2016b) and the "fit" function with initial values for k (0.5) and τ 1/2 (initial estimates vary for each strain), constraining both variables to positive values (see also Supplementary Fig. 5 for examples of sigmoidal fits to the data).
Coculture experiments. Cell counts were obtained by sampling 100 µl from 96well culture plates (inoculated with 170 µl minimal media, 2 × 10 µl of the normalized bacterial culture, and 10 µl particles as described above). Imaging was performed as described above. For CFU counts, samples were vortexed thoroughly to detach cells from particles and 10 µl were plated in 10 −2 and 10 −3 dilutions in replicates on MB2216 agar plates using rattler beads (Zymo S1001). After 72 h, colonies were counted to obtain CFUs.
DNA quantification. To quantify DNA as a proxy for biomass from monocultures in 96-well plates, wells were mixed thoroughly by pipetting and 100 µl of each well (including the particles) were sampled and frozen at −20°C for subsequent analysis. Cells were lysed by thawing and boiling (95°C, 10 min) 10 µl of each sample. Lysed samples were diluted 1:10 in TE buffer and quantified using Quant-it pico green (Fisher # P7589) standard protocols.
Modeling of particle half-life for noncooperative degraders. The degradation of a chitin particle by bacteria can be modeled in a simple way by taking into account two processes: free-living bacteria attach to the particle surface at a rate, a, proportional to their planktonic concentration, [B] 0 , such that a = a 0 [B] 0 , where a 0 is the attachment rate per bacterial cell, and attached cells degrade the particle at a rate p, and chitin monomers are converted to bacterial biomass at a conversion factor r. Note, that the conversion factor r may take into account the loss of monomers to the environment. This results in a set of differential equations for the amount of bacterial biomass, B(t), and the total amount of particle, R(t), In the above parametrization, the degradation of a chitin particle is described by four independent parameters: the total size of the particle, R 0 ; the attachment rate of bacteria, a; the biomass conversion rate, r; and (iv) the degradation rate; p. By rescaling the total number of bacteria by the growth rate, we obtain a more canonical parametrization: let b(t) = B(t)/r, α = a/r, and β = rp, then, which, for initial conditions b(0) = B(0) = 0 and R(0) = R 0 , are solved by the equations b t ð Þ ¼ α β e βt À 1 À Á ð6Þ From Eq. (6) the time required to fully degrade the particle can be found by numerically solving the transcendental equation for R(T) = R 0 . Additional analytical insight can be gained, however, by assuming that attachment is slow compared to growth, and hence R(T) ≈ R 0 − b(T). Furthermore, we additionally assume that rpT ⪢ 0 and, hence ln(e rpT −1) ≈ rpT. This leads to a simple expression for the total degradation time, T, required to fully degrade a particle, Hence, the degradation time, T, depends linearly on the logarithm of the attachment rate a, and hence the planktonic concentration of bacteria (see Supplementary Fig. 6).
There is an intuitive interpretation of the degradation time from Eq. (7): the increase in bacteria on the particle (Eq. 3) can be subdivided into two regimes. In the first regime, b ≈ 0, such that db dt % α and attachment dominates. In this regime, the number of bacteria on the particle increases linearly with time, b(t) ≈ αt. In the second regime, once a certain density of bacteria has been reached, i.e., b(t) ⪢ α/β, growth takes over as the dominant contributor, such that db dt % βb t ð Þ, and b(t)~e βt . These two regimes correspond to the terms ln 1 α and lnβ in Eq. (7).
Note, that in this formalism the number of bacteria could increase indefinitely. In reality, however, it is plausible that there exists a maximum carrying capacity for bacteria on the particle. Any effect of total carrying capacity will only kick in in the second regime, and will thus not affect the contribution of attachment to the degradation time.
Expected half-life for cooperative degraders. To calculate half-lives in the presence of cooperative growth we used a simple extension to the model, We solve these equations numerically to find the time it takes to consume half of the resources, R(τ 1/2 ) = R(0)/2, as a function of α (R(0) = 10 5 , β = 0.005). Simulations were performed in R using the deSolve package.
Genome sequencing and annotation. Genomes were sequenced on an Illumina MiSeq instrument, assembled with the CLC Assembly Cell and annotated using the RAST 32 annotation pipeline with default parameters.
Data availability. Sequence data that supports the findings of this study have been deposited in the NCBI database under Bioproject # PRJNA414740 with genome accession numbers "PIZM00000000", "PDUR00000000", "PDUS00000000", "PIZL00000000", "PDUT00000000", "PDUV00000000", "PDUU00000000", "PIZK00000000". All other data supporting the findings of this study are available from the corresponding author upon request.