Hydrodynamics of sponge pumps and evolution of the sponge body plan
Abstract
Sponges are suspension feeders that filter vast amounts of water. Pumping is carried out by flagellated chambers that are connected to an inhalant and exhalant canal system. In ‘leucon’ sponges with relatively high-pressure resistance due to a complex and narrow canal system, pumping and filtering are only possible owing to the presence of a gasket-like structure (forming a canopy above the collar filters). Here, we combine numerical and experimental work and demonstrate how sponges that lack such sealing elements are able to efficiently pump and force the flagella-driven flow through their collar filter, thanks to the formation of a ‘hydrodynamic gasket’ above the collar. Our findings link the architecture of flagellated chambers to that of the canal system, and lend support to the current view that the sponge aquiferous system evolved from an open-type filtration system, and that the first metazoans were filter feeders.
Introduction
Many aquatic suspension feeders use cilia and flagella to generate feeding flows from which they capture prey particles. Unlike ciliary arrays that drive fluid tangential to their attachment surface (Brooks and Wallingford, 2014; Narematsu et al., 2015; Gilpin et al., 2020), flagellated collar-cells (choanocytes) in sponges must pump water perpendicular to and through a perforated surface, working against a pressure resistance. Several choanocytes are closely packed inside chambers that are connected to inhalant and exhalant canal systems. The canal architecture and its associated pressure resistance impacts the operating condition of choanocyte chambers, and dictates morphological adaptation in such pumping units. Most sponges possess an elaborate aquiferous system (leucon) (Reiswig, 1975a), but some Calcarea, a small class of the Porifera, have a less elaborate aquiferous system (ascon and syconoid) (Manuel, 2006) Figure 1.
Asconoid forms are organized as a single tube in which the choanocytes form a single layer on the inner surface of the tube wall. Syconoid forms, in turn, possess many ascon-shaped cylindrical chambers branching off a central cavity (atrium). Water enters through short canals or directly through openings (ostia) into the chambers and leaves through a single outlet (apopyle) to the atrium and from there out of the osculum. In leuconoid forms, the holes in the surface of the body are connected to a complex canal system with a multitude of small chamber pump units.
Although these are often depicted as a progressive gain in complexity in the evolution of leuconoid from ascon and sycon forms, phylogenetic analyses indicate that in Calcarea both leucon and sycon type sponges most likely arose independently from ascon forms (Manuel, 2006), and some evidence points to the opposite process in which ascon and sycon forms may have derived from a rigid leucon body architecture (Dohrmann et al., 2006). It is unknown what gave rise to the leuconoid forms of the other three sponge classes. One of the greatest modern puzzles is whether Porifera or Ctenophora branched first in the evolution of multicellular animals (Whelan et al., 2015; Telford et al., 2016; Feuda et al., 2017). Therefore, understanding how the sponge filtration system came about would go a long way to helping resolve whether sponges evolved from a colonial unicellular ancestor (Porifera first hypothesis) or from a tissue-grade metazoan (as in some interpretations of the Ctenophora first hypothesis).
Both choanocytes and choanoflagellates, a group of free-living unicellular and colonial flagellates, share ancestory to animals, and possess similar collar-flagellated cells, to generate adequate feeding currents. In both cells, presence of a flagellar vane, a sheet-like structure along the length of the flagellum, has been observed (Afzelius, 1961; Fjerdingstad, 1961; Mehl and Reiswig, 1991; Weissenfels, 1992; Leadbeater, 2006; Mah et al., 2014). Vanes differ from mastigonemes, which are extracellular nanometer-thick fibers that protrude perpendicular to the flagellum (Bouck, 1971; Fenchel, 1982), both morphologically and functionally. Mastigonemes appear to exist and remain in the beat plane and are thought to reverse the thrust generation by the flagellum (Holwill and Peters, 1974). Vanes, however, are horizontal fibers of glycocalyx that extend perpendicularly from the flagellar axis and appear to augment the thrust generation (Nielsen et al., 2017).
Choanocytes in sponges, in addition to driving the flow through the collar, must also overcome the pressure resistance arising from the canal system. The narrow canals in leuconoid body plans add significant pressure resistance to the pump, between 2 and 50 Pa for different sponge species (Larsen and Riisgåd, 1994; Leys et al., 2011; Ludeman et al., 2017). Our recent study (Asadzadeh et al., 2019a) demonstrated that directional flow and efficient filtration at such hydrodynamic conditions is only possible because the pumping unit is sealed against back flow, by three extracellular elements, namely, a tight mesh on the distal part of the collar, a vane on the flagellum, and finally a gasket-like structure called Sollas’ membrane (Dendy, 1888; Sollas, 1889), that connects the microvilli (collar elements) at their tips (Figure 1). These features divide the chamber into two regions of low and high pressure, and force the inflow through the collar slits at its base. In contrast, calcareous sponges have fewer canals and thus have less pressure resistance, and appear to lack the gasket and tight mesh that covers the collar (Dendy, 1888; Eerkes-Medrano and Leys, 2006).
Here, we ask how sponges with apparent absence of a gasket, and therefore an open aquiferous system, can nevertheless effectively pump and filter water for prey, and what the trade-offs and limitations in such pumping systems are. We combine modeling with experiments to study the pumping and feeding in calcareous sponges using Sycon coactum as a model organism. Because of the geometric complexity, rather than attempting to use singularities to model the low Reynolds flow, as done for example for flows near a hole in a plane wall (Davis et al., 1981), or flagella-driven flows (Roper et al., 2013), we use computational fluid dynamic (CFD) simulations to investigate the functionality of the pumping elements in the absence of a gasket. Specifically, we study the role of each pumping element in pumping sufficient water through the ostia and into the chamber and examine the mechanism that forces the inhalant flow to pass through the collar filter. The modeled pumping rates are compared to measured exhalant flow rates and the predicted retention efficiency of calcareous sponges is contrasted with actual retention efficiencies measured for different prey types. The limitations and trade-offs between the pumping rate and retention efficiency shed light on what may have been the first poriferan filtration system.
Morphology of the flagellated chamber
Each flagellated chamber in Sycon coactum contains thousands of choanocytes that pump water through hundreds of ostia (Leys and Eerkes-Medrano, 2006) Appendix 1—figure 1. To reduce the computational costs, we consider a cylindrical axisymmetric and axially periodic chamber. The computational domain is thereby reduced to a wedge section into the cylindrical chamber that includes a central ostium surrounded by 24 choanocytes (Figure 2). The boundary conditions (BCs) are periodic along the long axis of the chamber (due to the periodic occurrence of ostia) and laterally symmetric on the left and right sides of the section modeled (if one is facing the long axis), due to cylindrical symmetry of the tubular chamber. Once fluid leaves the computational domain vertically into the inner cylindrical core of the choanocyte chamber, it is directed axially toward the open end of the chamber (apopyle). Pressure variation inside the core of the chamber is negligible (Appendix 1—figure 2), hence a uniform pressure boundary is applied on the top surface of the computational domain.
The collars are modeled as a porous structure (Equation 6) with a porosity that corresponds to that of a network of parallel and equally spaced cylinders (Keller, 1964) representing the microvilli of 0.1 µm in diameter. Although recent observations (Appendix 1—figure 2) have shown some mesh structures appearing locally between collars, chambers in S. coactum appear to lack a solid gasket. Hence we do not include a physical gasket in the CFD model. However, to study the effect of the gasket, we also consider cases where a physical gasket is included in the model. (Appendix 2—figure 1). Each flagellum is modeled as a thin sheet of width W (resembling the vane, Appendix 1—figure 2) that beats in a plane and is subject to a no-slip and no-penetration boundary condition with a prescribed motion given in Equation 3.
We initially consider mean measured values for the dimension of the chamber elements: spacing between flagella of 5 µm, collars of diameter , and length with spacing between microvilli of , and an ostium of diameter and length . Flagella have a wavelength of , an amplitude of (with length scale , see Equation 3) and vane width of 0.7 µm, all beating in phase at a frequency of 30 Hz (base case). In subsequent experimental simulations, these dimensions are varied along with the dimensions of other elements involved in the pumping system (e.g. length of the collar, length and diameter of the ostium, porosity of the collar etc.).
Results and discussion
Pumping mechanism
To facilitate efficient feeding, particles carried in the inhalant flow should reach the collars where the prey particles are filtered (Reiswig, 1971; Reiswig, 1975b; Weissenfels, 1976; Imsiecke, 1993). To investigate how the flow is forced through the collar rather than bypassing it in the absence of a gasket, we examine the average velocity field in two perpendicular planes in the middle of the computational domain (Figure 3AB). The velocity fields reveal a large backflow in the middle of the domain above the ostium where the spacing between the flagella is large. The pressure drop after the ostium provides a suction drawing the water from the higher pressure region above the flagella to the spacing between the flagella. This flow meets the inhalant flow from the ostium nearly at the same height as that of the collar, resulting in a stagnation zone that acts as a ‘hydrodynamic gasket’ that forces the inhalant flow toward the collars. Videos 1 and 2 demonstrate how formation of the hydrodynamic gasket ensures an efficient encounter of particles with the collar and prevents particles from bypassing the collars in the absence of a physical gasket.
The flow passing through the collar (filtration rate) is several times higher than the flow through the ostium (pumping rate) (Figure 3C), suggesting that a large proportion of the water is re-filtered inside the chamber. Increasing the vane width increases both the pumping and filtration rate, but the vane width has a larger effect on the latter (Figure 3C), likely due to the confinement and hydrodynamic interaction between the individual flagellum and its associated collar (Nielsen et al., 2017; Asadzadeh et al., 2019b). For instance, increasing the vane width from 0.3 to 0.7 µm increases the pumping rate by 53%, while the filtration rate is more than doubled.
To elucidate the mechanism responsible for pumping water through the ostia and to differentiate it from the mechanism that drives the water through the collars, we consider two hypothetical cases: (1) with flagella that beat only above the collars (Figure 4A) and (2) with short flagella confined to the collar height (Figure 4B). Flagella that are only present outside and above the collars are as efficient as complete flagella in pumping water into the chamber but are not able to drive water through the collar Figure 4C. In contrast, short flagella that are confined inside the collars are almost as efficient as complete flagella in driving water through the collars but are not able to pump new water through the ostium into the chamber. These results demonstrate two separate contributions from the confined and unconfined part of the flagellum of calcareous sponges. This pumping mechanism is distinctly different from that of demosponges and glass sponges, both of which possess a tissue or mucus gasket, in which the confined part of the flagellum with vanes that span the full diameter of the fine-meshed collars (Mah et al., 2014) is the pumping unit (Asadzadeh et al., 2019a).
Pump characteristic curve
To analyze the functionality of the pump at different hydrodynamic conditions, we consider the pump characteristic curve. In so doing, one should differentiate the basic pumping unit (pressure-generating unit) from the pressure-resistive parts of the system (the canal system). The basic pumping units are the simplest subdivision of the choanocyte chamber exposed to the same pressure difference resulting from the canal system. Hence, these units all work in parallel. For sponges lacking a physical gasket, the unit is one hole (ostium of zero length) with several neighboring choanocytes (Figure 2B). Presence of a physical gasket, however, leads to two zones of high and low pressure inside the spherical choanocyte chambers (Weissenfels, 1992; Asadzadeh et al., 2019a), effectively exposing each individual choanocyte to the same pressure difference. Therefore, in sponges with a physical gasket, the basic pumping unit is only one choanocyte, irrespective of its proximity to the nearest ostium (Asadzadeh et al., 2019a).
Basic pumping units are characterized by a maximum pumping rate (at zero pressure load), and a maximum pressure (at zero net flow). Due to the linearity of the governing equations at low Reynolds numbers, the pump characteristic of the unit is linear, hence:
where and are normalized pressure and pumping rate, respectively. Assuming Poiseuille flow in the canal system, and a tubular ostium of length and diameter , the pressure resistance (of the system) is given by White, 2011:
where is the canal resistance, µ the dynamic water viscosity, and is a characteristic of the pumping unit. Figure 5 depicts the dimensionless pump (Equation 1) and system (Equation 2) characteristic curves for different sponge pumps. Intersection of the two characteristics defines the operating condition of the pumping unit. Employing typical dimensions of the ostium and of the syconoid pumping unit (e.g. for the pumping unit in Figure 2B) results in , indicating at the operating condition. The higher the resistance from the canal system, the higher the slope of the system curve, hence the lower the pumping rate. As shown in Figure 5, the functionality of the syconoid pumping unit would be impaired if connected to a complex and narrow canal system with two orders of magnitude higher pressure resistance (e.g. ). For such hydrodynamic conditions, a modification of the pumping unit to that of the leucon type employing a physical gasket is inevitable. Thanks to the sealing elements, in such pumps are much higher (typically 3 orders of magnitudes) than syconoid ones (Asadzadeh et al., 2019a), resulting in a much lower slope of the system curve () yielding at the operating condition. Note that the actual characteristic curve in the leucon pump is non-linear, which has been ascribed to the bending of the vane at the relatively higher working pressure (Asadzadeh et al., 2019a).
Although sponges with a complex canal system appear to have a 'gasket-present’ pumping unit, a ’gasket-absent’ unit can still function efficiently if connected to a complex, yet open and less resistive canal system, or inefficiently if connected to an open canal system, but with highly resistive ostia. Therefore, there could well be sponges characterized as leuconoid types yet missing a physical gasket, or as syconoid types but with a physical gasket. The appearance, structure, and nature of the gasket would also depend on the hydrodynamics dictated by the canal system. Therefore, further research, using better preservation and imaging techniques, is required to elucidate the exact nature of the sealant element, and to investigate the role that hydrodynamics is playing in forming the seal in different sponges.
Hydrodynamic trade offs in the morphology of sponges
The open architecture of the syconoid pump dictates a delicate trade-off between the pumping rate and the retention efficiency. Although a wider and shorter ostium provides less resistance to flow hence higher pumping rate, the reduced pressure loss comes with a cost of a weaker backflow in the core of the region and thus a weaker hydrodynamic gasket located further into the chamber and above the collar tips. Under such circumstances, a larger fraction of the flow bypasses the collars and the retention efficiency drops. Lower retention efficiency of the calcareous sponge Pericharax heteroraphis as compared to different species of demposponges has been previously ascribed to its wider ostia and large chambers (Wilkinson, 1978). In contrast, a narrow and long ostium with higher resistance to flow results in a reduced pumping rate but an increased retention efficiency as the hydrodynamic gasket is located below the tip of the collars and the stronger backflow is more efficient in preventing flow from bypassing the collar filter (Appendix 1—figure 4).
The dependency of the pumping rate into the chamber on the length of the collar suggests another trade-off: choanocytes with shorter collars and a longer unconfined part of the flagellum can pump more water into the chamber, but this comes at the cost of more flow bypassing the collar, and vice versa. As a result, regions of the chamber having shorter collars will be less efficient in filtering the inhalant flow. This effect may explain the presence of pseudopodial extensions extending from the cell surface in these regions (Leys and Eerkes-Medrano, 2006). The pseudopodia, which are on average twice as long as the collars, have been observed to reach beyond the collar and seem to be involved in particle capture (Leys and Eerkes-Medrano, 2006), potentially compensating for the decreased filtration efficiency of the collars.
On individual collars, the spacing between adjacent microvilli appears wider at the base but becomes smaller toward the tip of the collar where the microvilli tips are occasionally fused together (Leys and Eerkes-Medrano, 2006). To study the effect of non-uniform porosity along the collars on the performance of the sponge pump, we consider cases where the spacing between microvilli decreases from the base to the tip of the collar. We find that more widely spaced (higher porosity) microvilli on collars enhance the filtration rate with a minimal effect on the pumping rate into the chamber (compare the plots in Appendix 1—figure 5). This result further highlights the fact that the pumping rate into the chamber is nearly independent of the confined part of the flagellum, while the filtration rate is greatly affected by the local collar-flagella parameters.
Flagella kinematics
Thus far, we have considered a synchronous beat and constant amplitude. We now consider different scenarios, that is, phase shift among the flagella, flagella beating in different planes (Appendix 1—figure 6), different frequencies, and different wavelength (Appendix 1—figure 7 and 8). The results show that beat synchronization among choanocytes is not necessary for the functionality of the pump. Despite improved efficiency of pumping in ciliary arrays by synchronization (Niedermayer et al., 2008; Golestanian et al., 2011; Elgeti and Gompper, 2013), asynchronized flagella in sponges continue to pump efficiently into the chamber (discussed in details in the SI). Lack of synchronization has been reported both within colonies of the choanoflagellate Salpingoeca rosetta (Roper et al., 2013; Kirkegaard et al., 2016), and in choanocyte chambers of the freshwater sponge Spongilla lacustris (Mah et al., 2014).
The amplitude of the flagella waveform is limited inside the collars, but it can increase outside, a behavior observed in choanocytes of the S. lacustris (Mah et al., 2014). To study the effect of increased amplitude, we modify the beat form in Equation 3 by setting and , which results in the increased amplitude of the waveform to a maximum of 2 µm at the tip of the flagella. This modulation in the beat form improves pumping rate by 51% (to ) while keeping the stagnation area at the same height relative to the base (Figure 3A), but it is only 37% more energetically demanding, suggesting that such a modification is beneficial to the sponge.
To compare the CFD predictions with observations, we estimate the volume flow rate per ostium in Sycon coactum by dividing the experimentally measured exhalant flux rate by the estimated total number of ostia. This results in an estimated flow rate of per ostium, and per choanocyte (given 24 choanocytes). CFD results show that ascon and sycon type sponges with typical values and dimensions of the morphological elements, are able to pump more than through an ostium without any requirement for a physical gasket (e.g. for the case with increased amplitude and a vane width of 1.4 µm on the unconfined part of the flagellum, ). Furthermore, both experimental and CFD estimates of flow rate per choanocyte () are comparable to published estimates for leucon sponges that range from 17 to for different species of demosponges and glass sponges (Larsen and Riisgåd, 1994; Leys et al., 2011; Ludeman et al., 2017), suggesting similar pumping capacity despite different pumping mechanism.
Effect of a gasket
To study the effect of a physical gasket on the flow and pumping rate, we model this structure as an impermeable baffle and incorporate it into the computational domain (Appendix 2—figure 1). The inclusion of a physical gasket improves the pumping rate by 21% if the collar is long (4.8 µm), but it does not alter the pattern of the averaged flow much (Figure 6AB). However, with shorter collars, a physical gasket decreases the pumping rate, and the more so the shorter the collar (Figure 6C). This phenomenon may explain the shorter collars in calcareous sponges compared to those in demosponges (Leys and Eerkes-Medrano, 2006) .
Retention efficiency
Simulations of passive prey particles (Videos 1 and 2) show that the prey retention efficiency is potentially 100%. However, these simulations assume that all particles encountering the filter are retained. In reality, however, the retention efficiency can vary for different prey types, especially in the absence of a gasket sealing off the collar filter area from the rest of the chamber. We studied experimentally the retention efficiency in Sycon coactum of three different prey types: small micro-algae (Euk < 20 µm), and non-photosynthetic bacteria with high (HNA) and low (LNA) nucleic acid content. While Euk and HNA were retained with near 100% efficiency, the efficiency of retention for the LNA prey particles was considerably lower (60 ± 6%) (Figure 7). LNA is largely associated with SAR11, the smallest and most abundant bacterial clade in the ocean (Mary et al., 2006). These bacteria may slip through the mucus filter of tunicates (Dadon-Pilosof et al., 2017; Dadon‐Pilosof et al., 2019) and evade filtration in many sponges (Ribes et al., 2012).
Evolutionary implications
Our modeling and experimental analyses indicate that the morphology of different sponge body plans results in different hydrodynamics with associated trade-offs in the sponge pump in terms of volume filtered and retention efficiency. Although the presence of sealing elements is crucial for the functionality of the high-pressure leucon sponge pump, the open architecture of the ascon and sycon type pumps is as efficient at particle capture without these structures. Taken in this light, our findings have implications for interpreting the first poriferan body plans. Neoproterozoic oceans were food and oxygen poor, with smaller bodied animals (Sperling and Stockey, 2018). In a Neoproterozoic ocean, competition with colonial flagellates could have favored a filter-feeder that captured larger particles without clogging, as is the case in ascon and sycon forms. With the Cambrian explosion came additional food sources providing the energy to enable specialization of the pump and filter to capture food in a range of habitats. Our findings demonstrate that the presence of sealing elements around the sponge collar is directly related to the hydrodynamics of operating conditions associated with the body architecture, and comes with delicate trade-offs between the sponge pumping rate and retention of particles. These analyses support the view (Manuel et al., 2003; Cavalier-Smith, 2017; Nielsen, 2019) that the sponge aquiferous system evolved from an open-type ascon-like filtration system, and the idea that the first metazoans were filter feeders.
Materials and methods
Numerical simulations
Request a detailed protocolWe use computational fluid dynamics (CFD) to solve the governing Navier-Stokes equations of the flow inside the computational domain (Figure 2). A finite volume method is used to discretize and solve the equations on a discrete representation of the computational domain consisting of polyhedral cells by applying the commercial CFD code STAR-CCM+ (14.04.013-R8). We use mesh morphing along with the overset method to move the computational mesh (Appendix 2, Section 2). The morphing motion redistributes mesh vertices in response to the movement of the flagellum. Video 3 shows the moving mesh (viewed in the xy-plane), inside the computational domain. To ensure that the solutions are independent of the computational mesh sizes, the simulations (for the base case) have been repeated using four coarser and one finer meshes (Appendix 2—figure 5). To demonstrate the independence of the CFD results from the global boundary conditions (BC) of the flagellate chambers (i.e. the chamber being closed at one end and open at the other end (apopyle), one simulation was conducted in an extended slice of a cylindrical domain including five ostia, closed (no-slip BC) at one end and open (uniform pressure BC) at the other end (Appendix 2—figures 3 and 4)).
Measurement of the pumping rate and retention efficiency in live sponges
Request a detailed protocolTen specimens of Sycon coactum were retrieved by scuba divers and transported at controlled temperatures to the laboratory at the Bamfield Marine Sciences Center, Bamfield, BC, Canada. Specimens were cleaned of macro-epibionts and left in large tanks with a flow-through system until the experiments were performed.
The volume flow rate through the oscula was calculated as the product of the ex-current jet speed and osculum area assuming a plug flow profile across the oscula (Dye Speed, DS, methods see Morganti et al., 2017 and references therein). Briefly, the excurrent jet speed was measured by releasing small amounts of filtered (0.2 µm) seawater mixed with sodium fluorescein dye next to the sponge ostia and videotaping the dye front in the excurrent jet along with a known scale. This procedure was repeated 5–10 times per sponge, and the DS was measured by recording the time it took the dye front to travel a short distance (10–20 mm). The osculum area was measured from photographs of each osculum using ImageJ (Schindelin et al., 2012).
To measure the retention efficiency of the sponges, we used a direct comparison of prey cell concentration in the water inhaled and exhaled by the sponge as described by Yahel et al., 2006. For the In-Ex method experiments, samples were processed in pairs (inhaled and exhaled samples). All samples were spiked with 1.0 µm Polysciences Inc Fluoresbrite yellow-green beads (Cat # 17154) from a stock solution that had been pre-calibrated with Becton Dickinson Trucount Control beads (Cat # 340335). Instrument flow rate was determined using bead counts. For the current analysis, we excluded the picocyanobacteria as their numbers were low, and identification was uncertain. Retention efficiency (%) was calculated as: , where and are the cell concentration in the inhaled and exhaled water, respectively.
Appendix 1
The calcareous sponge Sycon coactum
The general anatomy of S. coactum is illustrated in Appendix 1—figure 1. The sponge forms a single tube, 8 cm long, with a wall thickness of ∼1 mm and an inner diameter of 8 mm. A 1 mm2 piece of the wall was estimated to contain 36 chambers. Each chamber, in turn, was estimated to posses ∼283 ostia, each surrounded by 20–25 choanocytes. An excurrent jet velocity of 7.5 mms−1 was measured by dye visualization of the jet from an osculum with a diameter of 2 mm (Appendix 1—video 1).
Flagella beat model
The waveform of each individual flagellum is modeled as:
where d is the local lateral displacement of the each flagellum, a is the amplitude modulation of the waveform, is the characteristic length scale of the amplitude, Y is the coordinate along the centerline axis of the flagellum, Yb is the coordinate of the flagellum at its base on the cell surface, is the wave speed, where f is the frequency, and the wavelength, and t is time. Note that in the model of Equation 3, the flagellum would be extensible and its length varies slightly during the beat cycle (ca. 2%), resulting in marginally wrong velocities on the flagellum. However, implementing Equation 3 in the numerical scheme is computationally cheaper than preserving the arclength. Re-running the simulation for the base case while the arclength is preserved shows that the flow field in the domain (Appendix 1—figure 3) is insignificantly affected with less than 4% difference in the pumping rate.
For the flagellate pumps, the Reynolds number (), the ratio of inertia to viscous forces, is small (2.9 × 10−3), employing water density and dynamic viscosity , flagella length of , flagella wave speed of . As a result, the unsteady inertia terms in the governing equations can be neglected and flow is quasi-steady and symmetrical in respect to the flagellum beat. Therefore, to obtain the average velocity and pressure field, the simulation was conducted only for half of a beat cycles. For the particle filtration, however, the simulation is conducted for a 145 beat cycles in order to track the particles inside the domain.
The mechanical power expenditure by the flagellum is given by:
where denotes the stress tensor, is the unit normal vector on the surface S pointing into the fluid, and is the surface area of the flagella.
Simulations of particle filtration
Videos 1 and 2 illustrate passive particles entering into the chamber through the ostium and carried by the flow. Particle color denotes its velocity according to the color scale at the bottom. Particles that arrive at the collar are removed from the simulation. Note that while most of the particles are carried almost directly toward the collars, some are carried upward, above the collars, and toward the tip of the flagella, but almost all the particles are eventually returned by the back eddies and are captured. Simulation duration is 145 beat cycles. Side view (Video 1). Top view (Video 2).
Effect of ostium length and diameter on the pumping
Assuming a fully developed flow in a tubular ostium of length and diameter , the pressure resistance is given by White, 2011: , where μ is the dynamic water viscosity (here , corresponding to seawater at 20°C with salinity of 30 PSU). Note that pressure resistance due to the apopyle is 2 orders of magnitude smaller than that of the ostium (hence negligible), given the averaged dimensions of the apopyle , and of the ostium , , and 283 ostia per apopyle. Due to the linearity of the governing equations at low Reynolds numbers, the pump characteristic of the unit is linear, hence: , where and are the maximum flow rate (at zero length ostia) and pressure delivered by the pump, respectively. Combining these two equations, the pumping rate of the unit as a function of the ostium dimension is:
Having and , Equation 5 accurately predicts the pumping rate for a given length of the ostium as tested against CFD results Appendix 1—figure 4C.
Porosity of the collar filter
We model the porosity of a collar filter composed of several microvilli of radius as
where is the velocity normal to the filter subject to pressure drop, , κ is the dimensionless porosity of a network of parallel and equidistantly spaced cylinders (Keller, 1964), and μ is the viscosity of the water.
For the case of non-uniform porosity along the length of the collars (Appendix 1—figure 5), the dimensionless porosity κ of the collar is modeled as , where κb is the dimensionless porosity at the base corresponding to the maximum spacing () between adjacent microvilli, the position along the collar with respect to the base, and the length of the collar.
Asynchronization among flagella
Choanocytes are closely packed in the chambers with spacing between neighboring flagella of ∼5–10 µm. At such relatively short distances, synchronization of flagella beat might be expected (Brumley et al., 2014), although lack of synchronization has been reported both within a colony of the choanoflagellate Salpingoeca rosetta (Kirkegaard et al., 2016; Roper et al., 2013), and in choanocyte chambers of the sponge Spongilla lacustris (Mah et al., 2014). To investigate the potential effect of synchronized and asynchronized flagella beat on the mechanical energy produced by the flagella, we consider three different scenarios where: (1) neighboring flagella beat slightly out of phase (phase shift angle between neighboring flagella of ), (2) neighboring flagella are completely out of phase (), and (3) neighboring flagella beat in different (perpendicular) planes.
The snapshots of flow fields (Appendix 1—figure 6AB) are highly affected by the position of the neighboring flagella with respect to each other, especially near and between the flagella. However, the pumping rate and flow through the collar filter - lower and upper data, respectively, in C - are insignificantly affected by asyncronization among flagella (Appendix 1—figure 6C).
Multi-cilia and ciliary arrays beat in phase (Wan, 2018) or slightly out of phase that is, metachronal waves (Niedermayer et al., 2008; Golestanian et al., 2011; Elgeti and Gompper, 2013) to produce strong flows tangential to the surfaces to which they are attached. When the cilia beating is completely out of phase, the tangential forces of neighboring cilia cancel each other and no flow is created. In contrast, the array of flagella that lines the inner walls of sponge chambers produces a flow that is moving perpendicular to the attachment surface (and through it). In addition, some flagella are located in proximity to the water sources (ostium) while others are more distanced. To investigate the relationships between flagella location, their synchronization, and the mechanical power produced we compare the base case in which all flagella are beating in phase to two cases where the beat of the flagella next to the ostium (eight proximal flagella) is not synchronized to that of the 16 distal flagella. In the first case, the proximal flagella are beating faster () than the distal flagella () and in the second case, the proximal flagella beat slower () and the distal flagella beat faster (Appendix 1—figure 7AB). Results show that pumping rate is not impaired by such asynchronous flagella beat. Moreover, pumping rate of the first case (where proximal flagella pump faster (Appendix 1—figure 7A)) is 22% higher than that of the second case where the distal flagella beat faster ( vs ) despite the smaller number of the proximal flagella. The improved pumping rate in the first case is ascribed to higher capability of the proximal flagella in pumping due to their proximity to the ostium (for the case of synchronized flagella beat they contribute 64% of the total pumping, despite being only a third of the flagella in the array). Additionally, faster proximal flagella leads to 33% less mechanical power expenditure ( vs , where is the averaged total power expenditure by all 24 flagella over the beat cycle for the -th case). Considering that pumping rate increases linearly with the beat frequency, whereas the power expenditure increases as the square of the beat frequency (, ), employing faster beating of the fewer proximal flagella can be a suitable compromise between the power expenditure and pumping rate.
The presence of both proximal and distal flagella, and their collective imparted upward momentum is essential for efficient pumping of the inflow through the collars. Without the presence of the distal flagella, the system continues to pump, but the operation of the hydrodynamic gasket is impaired and water flow through the collars is dramatically reduced. The vacant space created by the absence of distal flagella, shifts the location of the back flow from above ostium toward the vacant space. Consequently, a significant part of the inflow bypasses the collars and leaves the domain without being filtered (Appendix 1—figure 7C). Considering the minor effect of synchronous flagella beat in reducing the power expenditure, the required imparted momentum of the distal flagella can potentially also be achieved with flagella having a different waveform, for example different wavelength, as long as the averaged upward flow velocity due the presence of the flagella is able to maintain the hydrodynamic gasket in place (Appendix 1—figure 7D). We find that the volume flow rate is affected by the wavelength, with an optimum wavelength of 5 µm (Appendix 1—figure 8).
Appendix 2
Inclusion of gasket into the model
We model the gasket as an impermeable plate that forms a canopy at the tip of the collars as depicted in Appendix 2—figure 1.
Moving computational mesh
We use mesh morphing along with the overset method to move the computational mesh corresponding to the motion of the flagella. With the overset method, rather than moving the entire mesh, it deforms the mesh only around the flagella, so-called the overset region, which significantly reduces the computational cost. A background stationary mesh including the collars is also generated which is superimposed by the overlapping overset mesh. The two mesh regions are implicitly coupled, the field data are interpolated back and forth between the two meshes (i.e. the overset and background meshes) to generate a smooth solution at each iteration.
Visualization of the flow field in the inner core of the chamber
Once fluid leaves the computational domain vertically into the core of choanocyte chamber, it is accumulated and directed axially toward one end of the chamber (apopyle) where it is open. To illustrate this, we model the flow in the cylindrical core of the chamber Figure 2A in the main text. The flow enters the cylinder from the lateral surface and one end (velocity inlet boundary of 1 µms−1, in the range of the mean velocity of flow leaving the wedge geometry of Figure 2), and exit from the other end (uniform pressure boundary). Figure 2 shows the velocity and pressure fields in a cross-section oriented axially in the center of the domain. Lateral flow and shear close to sides (hence right above flagella tip) is relatively weak, and we assume that the flagella action is not strongly affected.
Independence of solution from the global boundary conditions
Appendix 2—figure 4A shows snapshots of the velocity field in three middle units (separated by the white dashed lines) of the domain. The velocity field is virtually periodic in all these units indicating the minor effect of the domain ends on each single units, independent of their location with respect to the ends. Appendix 2—figure 4B shows the snapshots of the pressure field in the middle unit. Despite high variation inside the ostium and between the flagella, the pressure is nearly uniform above the flagella.
Mesh convergence study
In order for the solution to be independent of the mesh sizes, the CFD simulations are conducted at different mesh size. Because the problem is quasi-steady, the mesh convergence study has been performed at one arbitrary time step. Appendix 2—figure 5 shows the independence of the solution from the mesh sizes for the base case.
Data availability
All data generated or analysed during this study are included in the manuscript and supporting files.
References
-
Hydrodynamics of the Leucon sponge pumpJournal of the Royal Society Interface 16:20180630.https://doi.org/10.1098/rsif.2018.0630
-
Hydrodynamic functionality of the Lorica in choanoflagellatesJournal of the Royal Society Interface 16:20180478.https://doi.org/10.1098/rsif.2018.0478
-
Origin of animal multicellularity: precursors, causes, consequences-the choanoflagellate/sponge transition, neurogenesis and the cambrian explosionPhilosophical Transactions of the Royal Society B: Biological Sciences 372:20150476.https://doi.org/10.1098/rstb.2015.0476
-
Surface properties of SAR11 Bacteria facilitate grazing avoidanceNature Microbiology 2:1608–1615.https://doi.org/10.1038/s41564-017-0030-5
-
Prey taxonomy rather than size determines salp dietsLimnology and Oceanography 64:1996–2010.https://doi.org/10.1002/lno.11165
-
Axisymmetric stokes flows due to a rotlet or stokeslet near a hole in a plane wall: filtration flowsJournal of Fluid Mechanics 103:183–205.https://doi.org/10.1017/S0022112081001298
-
Studies on the comparative anatomy of sponges. II. On the anatomy and histology of stelospongus flabelliformis, Carter, with notes on the developmentQuarterly Journal of Microscopical Science 29:325–358.
-
Ultrastructure and embryonic development of a syconoid calcareous spongeInvertebrate Biology 125:177–194.https://doi.org/10.1111/j.1744-7410.2006.00051.x
-
Ecology of heterotrophic microflagellates. I. some important forms and their functional morphologyMarine Ecology Progress Series 8:211–223.https://doi.org/10.3354/meps008211
-
The ultrastructure of choanocyte collars in Spongilla lacustris (L.)Zeitschrift für Zellforschung und Mikroskopische Anatomie 53:645–657.https://doi.org/10.1007/BF00339512
-
The multiscale physics of cilia and flagellaNature Reviews Physics 2:74–88.https://doi.org/10.1038/s42254-019-0129-0
-
Hydrodynamic synchronization at low Reynolds numberSoft Matter 7:3074–3082.https://doi.org/10.1039/c0sm01121e
-
Dynamics of the hispid flagellum of Ochromonas danica. The role of mastigonemesThe Journal of Cell Biology 62:322.https://doi.org/10.1083/jcb.62.2.322
-
Viscous flow through a grating or lattice of cylindersJournal of Fluid Mechanics 18:94–96.https://doi.org/10.1017/S0022112064000064
-
Motility of colonial choanoflagellates and the statistics of aggregate random walkersPhysical Review Letters 116:038102.https://doi.org/10.1103/PhysRevLett.116.038102
-
The mystery of the flagellar vane in choanoflagellatesNova Hedwigia. Beiheft 130:213.
-
Feeding in a calcareous sponge: particle uptake by pseudopodiaThe Biological Bulletin 211:157–171.https://doi.org/10.2307/4134590
-
The energetic cost of filtration by demosponges and their behavioural response to ambient currentsThe Journal of Experimental Biology 220:995–1007.https://doi.org/10.1242/jeb.146076
-
Choanoflagellate and choanocyte collar-flagellar systems and the assumption of homologyEvolution & Development 16:25–37.https://doi.org/10.1111/ede.12060
-
Phylogeny and evolution of calcareous spongesCanadian Journal of Zoology 84:225–241.https://doi.org/10.1139/z06-005
-
The presence of flagellar vanes in choanomeres of porifera and their possible phylogenetic implicationsJournal of Zoological Systematics and Evolutionary Research 29:312–319.https://doi.org/10.1111/j.1439-0469.1991.tb00676.x
-
Synchronization, phase locking, and metachronal wave formation in ciliary chainsChaos: An Interdisciplinary Journal of Nonlinear Science 18:037128.https://doi.org/10.1063/1.2956984
-
Early animal evolution: a morphologist's viewRoyal Society Open Science 6:190638.https://doi.org/10.1098/rsos.190638
-
Particle feeding in natural populations of three marine demospongesThe Biological Bulletin 141:568–591.https://doi.org/10.2307/1540270
-
The aquiferous systems of three marine demospongiaeJournal of Morphology 145:493–502.https://doi.org/10.1002/jmor.1051450407
-
Bacteria as food for temperate-water marine spongesCanadian Journal of Zoology 53:582–589.https://doi.org/10.1139/z75-072
-
Functional convergence of microbes associated with temperate marine spongesEnvironmental Microbiology 14:1224–1239.https://doi.org/10.1111/j.1462-2920.2012.02701.x
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
The temporal and environmental context of early animal evolution: considering all the ingredients of an "Explosion"Integrative and Comparative Biology 58:605–622.https://doi.org/10.1093/icb/icy088
-
Coordination of eukaryotic cilia and flagellaEssays in Biochemistry 62:829–838.https://doi.org/10.1042/EBC20180029
-
Size independent selective filtration of ultraplankton by hexactinellid glass spongesAquatic Microbial Ecology 45:181–194.https://doi.org/10.3354/ame045181
Article and author information
Author details
Funding
Villum Fonden (9278)
- Seyed Saeed Asadzadeh
- Poul Scheel Larsen
- Jens H Walther
NSERC (2016-05446)
- Sally P Leys
Villum Fonden
- Seyed Saeed Asadzadeh
Danish Council for Independent Research Natural Sciences (7014-00033B)
- Thomas Kiørboe
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We acknowledge support from the Danish Council for independent Research (7014-00033B) to TK, from the Villum Foundation for JHW. through research grant no. 9278 and from NSERC Discovery Grant 2016–05446 to SPL. The Centre for Ocean Life is supported by the Villum Foundation.
Copyright
© 2020, Asadzadeh et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 4,038
- views
-
- 396
- downloads
-
- 17
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
The rise of angiosperms to ecological dominance and the breakup of Gondwana during the Mesozoic marked major transitions in the evolutionary history of insect-plant interactions. To elucidate how contemporary trophic interactions were influenced by host plant shifts and palaeogeographical events, we integrated molecular data with information from the fossil record to construct a time tree for ancient phytophagous weevils of the beetle family Belidae. Our analyses indicate that crown-group Belidae originated approximately 138 Ma ago in Gondwana, associated with Pinopsida (conifer) host plants, with larvae likely developing in dead/decaying branches. Belids tracked their host plants as major plate movements occurred during Gondwana’s breakup, surviving on distant, disjunct landmasses. Some belids shifted to Angiospermae and Cycadopsida when and where conifers declined, evolving new trophic interactions, including brood-pollination mutualisms with cycads and associations with achlorophyllous parasitic angiosperms. Extant radiations of belids in the genera Rhinotia (Australian region) and Proterhinus (Hawaiian Islands) have relatively recent origins.
-
- Evolutionary Biology
Studying the fecal microbiota of wild baboons helps provide new insight into the factors that influence biological aging.