Asteraceae, one of the largest flowering plant families, are adapted to a vast range of ecological niches. Their adaptability is partially based on their strong ability to reproduce. The initial, yet challenging, step for the reproduction of animal-pollinated plants is to transport pollen to flower-visiting pollinators. We adopted Hypochaeris radicata as a model species to investigate the functional morphology of the typical floral feature of Asteraceae, a pollen-bearing style. Using quantitative experiments and numerical simulations, here we show that the pollen-bearing style can serve as a ballistic lever for catapulting pollen grains to pollinators. This can potentially be a pollen dispersal strategy to propel pollen to safe sites on pollinators' bodies, which are beyond the physical reach of the styles. Our results suggest that the specific morphology of the floret and the pollen adhesion avoid pollen waste by catapulting pollen within a specific range equal to the size of a flowerhead. The insights into the functional floral oscillation may shed light on the superficially unremarkable, but ubiquitous functional floral design of Asteraceae.

Animal-mediated pollination is a complex interaction between plants and pollinators with conflicting interests. Plants strive for efficient and reliable pollen dispersal, whereas pollinators seek floral rewards that can be harvested as efficiently as possible (van der Kooi et al., 2021; Westerkamp, 1996). Plants employ various strategies to attract pollinators and ensure efficient pollen dispersal. One common approach is to offer floral rewards, such as nectar and pollen, to entice pollinators and facilitate direct physical contact between the exposed pollen grains and the floral visitors. To enhance this process, many plants have diversified their morphologies to include features such as a corolla tube (Gurung et al., 2021; Katzer et al., 2019), staminal lever (Claßen-Bockhoff et al., 2004), trapping device (Bröderbauer et al., 2013) or even sexual deception through the use of petals (Vereecken et al., 2012). In contrast, there are other dynamic pollination strategies that do not involve direct contact between the presented pollen and pollinators, such as explosive pollen release (Edwards et al., 2005; Switzer et al., 2018; Taylor et al., 2006) and buzz pollination (Michener, 1962; King et al., 1996; King and Buchmann, 2011; Vallejo-Marín, 2019, 2022; Jankauski et al., 2022). In this paper, we focus on a ubiquitous floral feature of Asteraceae, with the aim of shedding light on previously unexplored dynamic pollen dispersal mechanics in the plant family.

The Asteraceae family, with over 23,000 species making up roughly 10% of all flowering plants, is widely distributed on every continent except Antarctica (Barreda et al., 2015). These plants are characterized by their tightly packed inflorescences, or flower heads, composed of numerous individual flowers, known as florets (Fig. 1). Pollination in the Asteraceae family is initiated by the release of pollen grains into the anther tube, which pushes or brushes the pollen out and exposes it at the stylar surface (Erbar and Leins, 2015; Leins and Erbar, 2006). The flowers of Asteraceae are visited by a diverse range of pollinators, and thus they are considered to be generalists (Alarcón et al., 2008). However, some visitors utilize their elongated mouthparts to obtain nectar without frequently contacting the exposed pollen grains (Fig. 2F–H), potentially leading to insufficient pollen transfer for successful pollination. Corbiculate bees have refined pollen-collection methods using adhesive saliva and a pollen-collecting apparatus (Parker et al., 2015; Matherne et al., 2021; Roberts and Vallespir, 2015; Snodgrass, 1910), but without strategies to prevent excess pollen collection, the costs of attracting such pollinators can outweigh the benefits. While dynamic pollination strategies, such as explosive pollen release and buzz pollination, have been documented in other plant families, they have not been documented in Asteraceae. In this study, we examined the feasibility of transferring pollen to pollinators without direct physical contact through the oscillation of a pollen-bearing style. Our objective was to uncover a previously undocumented mode of pollen dispersal in a typical flower of the Asteraceae family through quantitative experiments and numerical simulations.

Fig. 1.

Floral structures of Hypochaeris radicata. A whole flowerhead (left) composed of many florets (middle), at the base of which are the so-called filaments. Adapted from Ito and Gorb (2019a).

Fig. 1.

Floral structures of Hypochaeris radicata. A whole flowerhead (left) composed of many florets (middle), at the base of which are the so-called filaments. Adapted from Ito and Gorb (2019a).

Fig. 2.

Conceptualization of the current study. (A,B) Stylar deflection toward a pollinator results in direct pollen delivery to the pollinator (A), whereas deflection in the opposite direction enables ballistic pollen delivery to the pollinator (as in this study; B). (C–H) Interaction between insects and styles. A style (position indicated by blue triangles) is deflected by a limb of a bee, Apis mellifera (C–E), or that of a fly from the family Bombyliidae (F–H). Upon release, the style snaps back towards the insect. (I–K) Snapshots from a high-speed video of an oscillating style showing pollen dispersal after a larger deflection. DF, deflection; OS, oscillation; DU, pollen dispersal units.

Fig. 2.

Conceptualization of the current study. (A,B) Stylar deflection toward a pollinator results in direct pollen delivery to the pollinator (A), whereas deflection in the opposite direction enables ballistic pollen delivery to the pollinator (as in this study; B). (C–H) Interaction between insects and styles. A style (position indicated by blue triangles) is deflected by a limb of a bee, Apis mellifera (C–E), or that of a fly from the family Bombyliidae (F–H). Upon release, the style snaps back towards the insect. (I–K) Snapshots from a high-speed video of an oscillating style showing pollen dispersal after a larger deflection. DF, deflection; OS, oscillation; DU, pollen dispersal units.

Plant species

Hypochaeris radicata (Asteraceae) was previously adopted as a model species to investigate pollen adhesion (Ito and Gorb, 2019a,b). In this study, we used this species to examine the motion of pollen-bearing styles and the resulting pollen dispersal. Hypochaeris radicata is a perennial plant, native to Europe, and is a cosmopolitan invasive species occurring in a wide range of temperate zones, including America, Japan and Australia (Doi et al., 2006). Hypochaeris radicata is known to be self-incompatible. This means that the successful transport of pollen grains to different individual plants of the same species is necessary to enable its healthy reproduction (Picó et al., 2004). Flowering stems of H. radicata were collected in Kiel, Germany, for experimentation. They were placed in water until the youngest florets exposed fresh pollen.

Field observations

To observe the pollen-collecting behaviors of pollinators on flowerheads of H. radicata, we filmed videos in slow motion (120 frames s−1) using an iPhone 7 (Apple Inc.) together with a ×30 magnifying glass (Fig. 2C–E).

Mechanical characterization of florets

We collected the newly opened florets from the flower heads and analyzed their morphology under a light microscope (Leica Microsystems, Wetzlar, Germany) (Fig. 1). In order to characterize the mechanical properties of the florets, we measured the spring constant of the following segments: (1) styles, (2) anther tubes and (3) filaments. Fig. 3A illustrates the experimental setup used for this purpose. Each freshly opened floret was first horizontally positioned between two wooden blocks. We fixed the floret specimens at different positions to measure the spring constant of their different segments (Fig. 3B, right): (1) the entire anther tube was fixed to test the style (fixation at f1), (2) the basal part of the anther tube including the filaments was fixed to test the anther tube (fixation at f2), and (3) the basal part of the petal was fixed to test the filaments (fixation at f3). Floret specimens were deflected using a thin piece of metal mounted on a force transducer (10 g capacity; World Precision Instruments Inc., Sarasota, FL, USA). The deflections were always applied 1 mm distal to the fixation positions at a controlled displacement speed of 0.01 mm s−1. The force required to deflect the specimens was continuously recorded using AcqKnowledge 3.7.0 software (Biopac Systems Ltd, Goleta, CA, USA). To measure the spring constant of the floret specimens, we always used the data from the initial part of the obtained force–displacement curves (limited to a displacement equal to 50 μm). In total, we tested 34 florets including 9 tests on styles, 10 on anther tubes and 15 on filaments. Each floret was subjected to a single test and was not used again.

Fig. 3.

Mechanical characterization of florets. (A) The experimental setup. (B) Left: box plot of the spring constant of three different segments of florets (median, upper and lower quartiles and 1.5× interquartile range). Right: the fixation position (red) and the location of the applied deflection (arrows). Specimens were fixed at f1, f2 and f3 to measure the spring constants of the style, anther tube and filaments, respectively.

Fig. 3.

Mechanical characterization of florets. (A) The experimental setup. (B) Left: box plot of the spring constant of three different segments of florets (median, upper and lower quartiles and 1.5× interquartile range). Right: the fixation position (red) and the location of the applied deflection (arrows). Specimens were fixed at f1, f2 and f3 to measure the spring constants of the style, anther tube and filaments, respectively.

Stylar oscillation experiments

Black circular polyethylene plates with a hole at the center were prepared for visualizing the distribution of pollen grains catapulted by stylar oscillations. Each newly opened floret was vertically fixed through a hole in the plate onto the vise to stand upright (Fig. 4). To examine whether the pollen dispersal caused by the stylar oscillation depends on the magnitude of an initial deflection, an insect pin attached to a micro-manipulator (World Precision Instruments) was brought into contact with either the style or the anther tube of the standing floret (Fig. 4). Contact with the anther tube resulted in a large deflection of the floret, whereas contact with the style caused a smaller deflection. The insect pin was kept moving horizontally until the deflected floret was released (Figs 2I–K and 4). After release, the floret started to oscillate and this led to the release of clumps of pollen grains, here called ‘dispersal units’ (DU), from the style (Fig. 2J,K; Movie 1). The oscillations were filmed by a high-speed camera (Olympus, Tokyo, Japan) at 5000 frames s−1 and tracked using open-source tracking software (Tracker by Douglas Brown; https://physlets.org/tracker/) (Brown and Cox, 2009). In total, we analyzed the stylar oscillations of 55 specimens, among which 24 specimens were used to analyze the distribution and the number of pollen grains catapulted by the oscillations.

Fig. 4.

Experimental setup for the oscillation experiments. The initial deflections, which triggered the stylar oscillation, were applied to two different positions (style and anther tube) to cause two distinct deflection magnitudes (a smaller and larger deflection, respectively). The stylar oscillations were recorded by using a high-speed camera. The dispersed pollen grains, catapulted by the stylar oscillations, landed on the circular substrates and were photographed from above.

Fig. 4.

Experimental setup for the oscillation experiments. The initial deflections, which triggered the stylar oscillation, were applied to two different positions (style and anther tube) to cause two distinct deflection magnitudes (a smaller and larger deflection, respectively). The stylar oscillations were recorded by using a high-speed camera. The dispersed pollen grains, catapulted by the stylar oscillations, landed on the circular substrates and were photographed from above.

Assuming linearity of the oscillatory system, the damping ratio, ζ, of a floret was obtained based on the logarithmic decrement, δ:
(1)
where x1 and xi are the amplitudes of the first peak and the peak that is i−1 periods away, respectively. By using Eqn 1, we obtained the damping ratio, ζ:
(2)
The dispersal units that landed on the plate were photographed under a microscope (Keyence, Osaka, Japan). The images were analyzed by a custom-written script in MATLAB (Mathworks, Natick, MA, USA) to obtain dispersal distances. To count the number of dispersed pollen grains, the grains on the plate were collected, placed in a droplet of mineral oil on a glass slide, and sandwiched with a glass coverslip to spread out as a single layer. They were then photographed under a microscope (Leica Microsystems). Using the images, the grains were counted by a custom-written script in MATLAB.

Numerical simulation of trajectories of dispersal units

Because of the technical difficulties of tracking single pollen grains based on the high-speed videos, we computed trajectories of dispersal units and their dispersal distances. During stylar acceleration, dispersal units on the style experience inertia. They could leave the style if the inertial force exceeds the attachment force of the dispersal unit on the stylar surface. The inertial force, Fi, required to detach a dispersal unit from the style can be obtained from:
(3)
where ac is the critical acceleration, Fa is the adhesion force of the pollen on the style, c is the number of contact points between a dispersal unit and the style, m is the mass of single pollen grains, and n is the number of pollen grains forming a single dispersal unit.

The adhesion of pollen grains has previously been measured for H. radicata on the stylar surface: median adhesion force, Fa, was 98 nN (N=50) (Ito and Gorb, 2019a). The mean (±s.d.) mass of individual pollen grains, m, was 15.0±0.4 ng (Ito and Gorb, 2019b). The number of pollen grains forming a single dispersal unit, n, was determined by dividing the number of dispersed grains by the number of dispersal units (n=12.9±5.3, mean±s.d., N=24). The number of contact points, c, was set as 2 to restrict major pollen dispersal to the first cycle of oscillations, as observed in experiments. Using these data, the critical acceleration that can detach a dispersal unit from the style was calculated, and is shown as a boxplot in Fig. 5C. In the numerical simulation, we used the median value of the critical acceleration (ac=1005 m s−2).

Fig. 5.

Mechanics of pollen dispersal. (A) Correlation analysis between the number of dispersed grains and each variable extracted from the motion data: acceleration a1 (left), frequency f (middle) and damping ratio ζ (right). (B) Time series of the stylar displacement x, linear acceleration a and velocity v. The first, second and third peak of linear acceleration is denoted as a1, a2 and a3, respectively. The corresponding velocity at the moment of each acceleration peak is denoted as v1, v2 and v3, respectively. The shaded half-cycles show the conditions where pollen should be catapulted towards a pollinator if pollen detachment from the style occurs. (C) Bar plots of the acceleration peaks in stylar oscillations after small or large deflections (left) and the corresponding velocity at each acceleration peak (right). The shaded region and horizontal dashed line on the left show the first and third quantile, and the median value of the critical acceleration, ac. The critical acceleration is the acceleration required to detach dispersal units from the style, calculated based on the measured pollen adhesion on the style and the average number of pollen grains in a dispersal unit (n=13). *P<0.05, **P<0.01, ***P<0.001.

Fig. 5.

Mechanics of pollen dispersal. (A) Correlation analysis between the number of dispersed grains and each variable extracted from the motion data: acceleration a1 (left), frequency f (middle) and damping ratio ζ (right). (B) Time series of the stylar displacement x, linear acceleration a and velocity v. The first, second and third peak of linear acceleration is denoted as a1, a2 and a3, respectively. The corresponding velocity at the moment of each acceleration peak is denoted as v1, v2 and v3, respectively. The shaded half-cycles show the conditions where pollen should be catapulted towards a pollinator if pollen detachment from the style occurs. (C) Bar plots of the acceleration peaks in stylar oscillations after small or large deflections (left) and the corresponding velocity at each acceleration peak (right). The shaded region and horizontal dashed line on the left show the first and third quantile, and the median value of the critical acceleration, ac. The critical acceleration is the acceleration required to detach dispersal units from the style, calculated based on the measured pollen adhesion on the style and the average number of pollen grains in a dispersal unit (n=13). *P<0.05, **P<0.01, ***P<0.001.

We measured the stylar velocity, v, and acceleration, a, at different time frames using the recorded high-speed videos (N=55) (Fig. 6A, red asterisks). When the stylar acceleration was greater than the critical acceleration, the trajectory of the dispersal unit was simulated (Fig. 6A, solid curves). The initial position and velocity of each dispersal unit were set to be equal to those of the styles. The initial position of a dispersal unit (x, y) is referred to here as the detachment point of that unit (Fig. 6B).

Fig. 6.

Trajectories of dispersal units and detachment points for stylar oscillations. (A) Calculated trajectories of dispersal units (solid lines, colored differently for clarity) based on the measured trajectories of stylar oscillations (red asterisks) initiated by different deflection magnitudes (small and large). (B) Estimated detachment points for stylar oscillations initiated by deflections of different magnitude.

Fig. 6.

Trajectories of dispersal units and detachment points for stylar oscillations. (A) Calculated trajectories of dispersal units (solid lines, colored differently for clarity) based on the measured trajectories of stylar oscillations (red asterisks) initiated by different deflection magnitudes (small and large). (B) Estimated detachment points for stylar oscillations initiated by deflections of different magnitude.

Here, we assumed the dispersal units to consist of 13 pollen grains (n=13) based on the measurements mentioned above. The diameter of the dispersal units can be calculated by (Huang and Yu, 2012):
(4)
where the radius of a single pollen grain, r, is equal to 15 μm (Ito and Gorb, 2019b).
Reynolds number at a given time, Re, is given as:
(5)
where v is the velocity of a dispersal unit at a given time, ρ is the density of air (1.204 kg m−3) and ϕ is the dynamic viscosity of air at 20°C (1.825×10−7 kg m−1 s−1). Considering the low range of Reynolds number (0.3<Re<8.5), the drag coefficient CD is not constant, and instead is calculated by the following empirical relationship (Jones and Knudsen, 1961):
(6)
Once CD is obtained, the drag force, D(t), can be calculated using the following equation:
(7)
where A is the frontal area of the dispersal unit and equal to πr2. Then, we obtain a system of first-order differential equations:
(8)
where the gravitational acceleration, g, is equal to 9.81 m s−2. We solved the system of the differential equations by using the Runge–Kutta fifth-order accurate (RK45) solver in the Python package SciPy. The dispersal distance of a dispersal unit was defined as the landing position on the x-axis with reference to the origin, i.e. the position of the style.

Functional segments of floret

Florets of H. radicata can be subdivided into three functional segments from top to bottom: an exposed distal style, an anther tube in the middle and five basal filaments (Figs 1 and 3B). Mechanical testing of the florets revealed distinct stiffness of the three segments. The style, which passes through the anther tube and is exposed at its distal segment, has the lowest spring constant (K=0.2±0.1 N m−1, N=9) among other segments. It is embraced by the stiff anther tube, which features the highest spring constant (K=6.5±2.7 N m−1, N=10). The anther tube is connected to the petal at its proximal part by five filaments, which altogether have an intermediate spring constant (K=2.31.9 N m−1, N=15). The spring constants of the three segments are significantly different from each other (Tukey multiple comparisons of means, filaments versus style: P=0.01; anther tube versus style: P=2.0×10−7; anther tube versus filaments: P=1.3×10−4).

Mechanics of pollen dispersal

The high-speed video analysis enabled us to investigate stylar oscillations. Fig. 2I–K shows three snapshots of the stylar oscillation from release to return. When the floret was deflected, the elastic energy was mainly stored in the proximal part of the flexible style and the filaments. However, no obvious deformation was observed in the anther tube (Fig. 2I). Upon the release, the style snapped back in the opposite direction of the applied displacement, initiating the first half-cycle of oscillation. The oscillation decayed quickly, and the style returned to the resting position (N=55, damping ratio ζ=0.182±0.05, frequency f=119±23 Hz) (Fig. 5B).

Based on the high-speed videos, we found that the first half-cycle caused the major pollen dispersal (Fig. 2K). The maximum acceleration in each half cycle decreased over time (Fig. 5B,C), so that the values were significantly different from each other (Tukey multiple comparisons of means, large deflections, a1 versus a2: P<0.001, a2 versus a3: P=0.0013, a3 versus a1: P<0.001; small deflections, a1 versus a2: P<2.2×10−6, a2 versus a3: P=0.0098, a3 versus a1: P<0.001).

As shown in Fig. 5C (left), for both large and small deflections, the maximum stylar acceleration in the first half-cycle, a1, exceeded the critical acceleration, ac. However, the maximum acceleration in the subsequent half-cycles (i.e. a2 and a3) largely overlapped with the critical acceleration, and the third peak, a3, mostly became lower than the median critical acceleration.

The velocity corresponding to the maximum acceleration in the first half cycle, v1, was by far the highest and significantly different from velocities corresponding to acceleration in the second and third peaks, v2 and v3 (Fig. 5C, right) (Tukey multiple comparisons of means, large deflections, v1 versus v2: P<0.001, v2 versus v3: P=0.97, v3 versus v1: P<0.001; small deflections, v1 versus v2: P<0.001, v2 versus v3: P=0.77, v3 versus v1: P<0.001).

We examined the relationship between the number of dispersed grains, nd, and three potentially influential parameters: frequency, damping ratio and maximum acceleration of the styles (Fig. 5A). While no correlation was found between the frequency or damping ratio and the number of dispersed grains (Pearson's product-moment correlation, f versus nd: t=−0.22, d.f.=22, P=0.82; ζ versus nd: t=−0.74, d.f.=22, P=0.47), we found a significant correlation between maximum acceleration and the number of dispersed grains (Pearson's product-moment correlation, maximum a versus nd: t=3.2, d.f.=22, P=0.004). The relatively low coefficient of determination (R2=0.32) can be attributed to variation in the number of pollen grains on the examined styles.

Effects of deflection magnitude on pollen dispersal

Large deflections resulted in significantly higher acceleration in the first two peaks, but not in the subsequent peaks (t-test, a1: t=3.90, d.f.=52, P=0.00028; a2: t=2.03, d.f.=52, P=0.047; a3: t=0.96, d.f.=52, P=0.34) (Fig. 5C, left). They also led to a significantly higher number of dispersed grains, in comparison with small deflections (t-test, t=−2.34, d.f.=16, P=0.034). This is reflected by the positive correlation between the maximum acceleration in the first half-cycle, a1, and the number of dispersed grains, nd, demonstrated in Fig. 5A (see above). In contrast, large deflection caused higher corresponding velocity at only the first acceleration peak (t-test, v1: t=6.10, d.f.=52, P<0.001; v2: t=1.13, d.f.=52, P=0.26; v3: t=1.66, d.f.=52, P=0.10) (Fig. 5C, right).

However, as shown in Fig. 7A, the different deflection magnitudes had no significant effect on pollen distribution (t-test, t=−0.96, d.f.=397, P=0.34). The stylar oscillation tended to catapult the dispersal units against the deflection direction with a median dispersal distance of 5.5 mm, for both large and small deflections. Surprisingly, 93.5% of the dispersal units landed within 15.8 mm of the origin, which is equal to the mean radius of flowerheads of H. radicata.

Fig. 7.

Distribution of dispersal units catapulted by stylar oscillations. Stacked histograms show distribution data for large and small deflections from quantitative experiments (N=55) (A) and numerical simulations (B).

Fig. 7.

Distribution of dispersal units catapulted by stylar oscillations. Stacked histograms show distribution data for large and small deflections from quantitative experiments (N=55) (A) and numerical simulations (B).

Although our simulations showed a broader peak of the dispersal distance (Fig. 7B), the simulated pollen distribution exhibited similar trends to the measurements: a tendency for pollen dispersal to be directed toward a pollinator and dispersal distances that were restrained within the size of the flowerhead.

Fig. 6A shows the representative simulated trajectories of dispersal units based on the tracked stylar oscillations (red asterisks), initiated by small and large deflection of the styles. Fig. 6B shows estimated detachment points, where Eqn 3 was satisfied. Large deflections widened the range of detachment points, in comparison to small deflections. They also shifted the detachment points towards lower values in both the x- and y-axes, compared with the small deflections (Kruskal test, x-axis: χ2=29.9, P<0.001; y-axis: χ2=45.2, P<0.001). This means that the larger deflections caused pollen detachment closer to the ground and further from a pollinator than the smaller deflections.

This study focused on the oscillation of the style triggered by an initial deflection and its potential as a ballistic lever that catapults pollen towards a flower-visiting pollinator. As the pollen catapult is a passive spring-driven system, adaption to a variety of deflection magnitudes is likely to be key to the successful delivery of pollen to a pollinator and the minimization of waste for dispersed pollen. Below, we discuss the potential role of morphological and mechanical features of the reproductive part of a floret as well as the adhesive properties of pollen in the pollen catapult system.

Floret: a functional composite

Similar to other species in the Asteraceae family (Meiri and Dulberger, 1986), H. radicata florets can be regarded as composites made from three morphologically distinctive segments (Figs 1 and 3B). It is known that each of the segments plays an essential role in pollen presentation (Erbar and Leins, 2015, 1995; Leins and Erbar, 2006). Do they also contribute to catapulting the pollen towards visiting pollinators, after pollen presentation? To present pollen and successfully deliver it to visiting pollinators, the floret with the pollen-bearing style should meet these requirements: (1) high mechanical compliance during contact to store elastic energy and prevent damage and (2) catapult pollen towards a pollinator by adapting to varying magnitudes of initial deflections.

At the distal end of a floret, the flexibility of the pollen-bearing style enables high compliance under contact as well as energy storage for the pollen catapult after release (Fig. 2I–K). The middle segment of the floret is reinforced by the anther tube, which functions as a support for the slender style during contact and stabilizes the shape of the floret during non-contact periods. Previous studies have shown that the deformation of a flexible cantilever with a high aspect ratio occurs proximally to the deflection points (Rajabi et al., 2018). Therefore, in the case of physical contact at the style, the elastic energy would be stored in the broad proximal area of the floret. However, the reinforcement of the proximal segment of the floret by the anther tube prevents the deformation of this segment, limiting the amount of energy stored and released at the distal end of the style. This is likely to restrict excessive stylar acceleration, and thus reduce the number of catapulted pollen grains.

The other energy-storing segment, consisting of filaments, is located at the base of the floret. This segment, with its intermediate stiffness, functions as a flexible joint. It enhances the mechanical compliance of the whole structure in contact, and therefore enables adaptation to the different deflection magnitudes. When a large initial deflection is applied, the flexible joint helps to recline the whole floret (Figs 2I and 6A). As mentioned earlier, larger deflections lower the detachment points and shift them further from a pollinator (Fig. 6B). This is likely to counterbalance the increased initial velocity of catapulted dispersal units (Fig. 5C, v1) so that pollen is still distributed within the same range of distances, independent of the initial deflection magnitude (Fig. 7A).

Pollen adhesion

Some animals, such as insects and elks, clean themselves by vibrating their hairs to remove accumulated particles (Amador and Hu, 2015; Amador et al., 2017). In such systems, the natural frequency governs particle detachment:
(9)
where a is the acceleration, S is the spacing between the hairs, which dictates the maximum deflection of the hairs, and f is the natural frequency (Amador and Hu, 2015). However, here, we did not find correlations between the natural frequency and the number of dispersed grains (Fig. 5A) or between the natural frequency and maximum linear acceleration (Fig. S1).

It is important to mention that the major pollen detachment occurs only in the first half-oscillation cycle. Hence, the general oscillatory frequency has very little influence on pollen detachment. To restrict the pollen detachment only to the first half-cycle, the critical acceleration, ac, should be less than the first maximum acceleration, a1, but greater than maximum acceleration in any other oscillation cycles (Fig. 5C, left).

As shown in Eqn 3, the critical acceleration depends not only on the pollen-style adhesion but also on the pollen–pollen adhesion, which governs the number of pollen grains in a dispersal unit. The weaker pollen–pollen adhesion, therefore, decreases the number of grains in a dispersal unit, in comparison to the n=13 measured here. As the inertial force is proportional to the mass of the dispersal unit, a smaller number of pollen grains in a dispersal unit can challenge its detachment from the style. If the number of grains is as small as n=3, the median critical acceleration would be greater than most of the measured maximum stylar accelerations, and this would completely incapacitate the pollen catapult. In contrast, higher pollen–pollen adhesion would result in a higher number of grains in the dispersal units. This would, in turn, lead to a prolonged and probably isotropic pollen dispersal (both towards and against the pollinator) and, therefore, increase the number of wasted pollen grains. Hence, for guided pollen dispersal towards a pollinator, pollen adhesion and the size of dispersal units should be kept within a certain range.

Unlike wind-pollinated species, insect-pollinated ones are covered by a thick layer of a viscous oily substance called pollenkitt (Pacini and Hesse, 2005), which helps to form large dispersal units by bonding pollen grains together (Pacini and Hesse, 2005; Lisci et al., 1996). In addition to the general function of pollenkitt as a ‘pollen adhesive’, we recently showed that it inhibits pollen adhesion by weakening water capillary attraction (Ito and Gorb, 2019b). Therefore, we can assume that pollenkitt may play a role in maintaining pollen adhesion within a specific range. Further studies are needed to quantitatively investigate the role of pollenkitt in pollen–pollen adhesion.

Phase shift in the first half-oscillation cycle

Once pollen detachment occurs, pollen dispersal is mainly governed by the initial velocity of the dispersal units, rather than their acceleration. Hence, the velocity corresponding to each maximum stylar acceleration (Fig. 5C, right) is also key to guided pollen dispersal towards a pollinator.

Unlike the stylar acceleration peaks, the corresponding stylar velocities abruptly decayed, so that the velocities corresponding to the second and third maximum accelerations (i.e. v2 and v3, respectively) were not significantly different from each other (Fig. 5C, right). Only the first corresponding velocity, v1, was significantly higher than the others, because of the phase shift between the stylar acceleration and velocity in the first half-cycle.

In a free oscillation, the displacement, x(t), is given as:
(10)
where A is an arbitrary constant and ω is the angular frequency. Differentiating Eqn 10 with respect to t gives the velocity v(t) as:
(11)
In the same manner, the acceleration, a(t), is given as:
(12)
The phase difference between the acceleration and velocity is equal to π/2, or a quarter of an oscillation cycle. Therefore, in the free oscillation conditions, the corresponding velocities at the moment of acceleration peaks, ideally, are the local minima. As shown in Fig. 5B, this typical phase difference appeared in the stylar oscillations, resulting in the extremely small values of v2 and v3. However, the first corresponding velocity, v1, was shifted from the local minimum; therefore, v1 was by far the highest of the corresponding velocities. The cause of the phase shift is seemingly the style sliding over the insect pin during the release, and therefore a detailed study on this topic will be required in the future.

Biological significance

Asteraceae are known for attracting various insect pollinators (Alarcón et al., 2008). Their robust adaptability makes a number of them notorious invasive species (Thebaud and Abbott, 1995; Hao et al., 2011). To enable their versatile reproduction, Asteraceae may use multiple strategies to secure pollen transfer to diverse pollinators. Each pollinator species may have its own difficulties in terms of being a functional pollen vector for the benefit of plants. Bees, for example, are both pollen-transporting vectors and pollen consumers. Actively collected pollen grains are soaked with their saliva and packed in a specialized pollen-carrying apparatus, and are therefore no longer available for pollination (Parker et al., 2015; Thorp, 2000). Some pollinators, such as flies from the family Bombyliidae, with their specialized long mouthparts and slender limbs, less frequently make physical contact with styles on flowers, yet are able to obtain nectary rewards (Kastinger and Weber, 2001) (Fig. 2F–H).

In this paper, we discussed that the superficially unremarkable floret of H. radicata, by optimizing its (1) morphological and mechanical properties, (2) adhesive properties of pollen and (3) phase difference between the stylar acceleration and velocity, functions as a projectile tool to catapult pollen towards pollinators. The airborne pollen dispersal system is advantageous for delivering pollen to various body parts of pollinators, including their safe sites, as well as engaging the long-limbed insects beyond the physical reach of the styles (Macior, 1964; Koch et al., 2017). The attached pollen on the pollinators is potentially transferred to the stigma of other flowers by electrostatic attraction (Vaknin et al., 2001; Eisikowitch et al., 2005; Bowker and Crenshaw, 2007a,b) or grooming (Amador et al., 2017). It is also remarkable that more than 90% of the dispersal units landed within a region of the same size as that of the flowerhead. Therefore, even if the dispersal units fail to reach the target pollinator, they remain within the flowerhead, waiting to be transferred to forthcoming pollinators.

It is well known that florets of Asteraceae consecutively mature from the outer row towards the central row. For many species of this family, including H. radicata, this means that young and most pollen-bearing styles are situated in the center of the flowerhead (Harris, 1995). The tuned travel distance of the dispersal units, therefore, may explain the order of floret maturity in Asteraceae. The working principles of the pollen catapult in Asteraceae are rooted in a common floral feature of the family. If this mechanism is found to be present in other species within the family, it could potentially contribute to the reproductive success of one of the largest families of flowering plants. However, future research is necessary to determine the prevalence of similar pollen dispersal mechanics in other species within the family.

Conclusions

In this study, we investigated the pollen dispersal strategy resulting from the stylar oscillation of H. radicata. Based on the combination of high-speed motion analysis, mechanical tests and numerical simulations, we found that the morphologically and mechanically distinctive segments comprising the floret, as well as optimal pollen adhesion, potentially contribute to (1) catapulting pollen towards a visiting pollinator while (2) minimizing wasteful pollen dispersal by (3) adapting to different degrees of oscillation-inducing external deflection. As ballistic pollen dispersal arises from a standard floral morphology of Asteraceae, it can potentially play a role in the pollination of the most ubiquitous flowering family on the planet.

The authors thank Alexander Köhnsen for the help in motion tracking, and the YouTube channel ‘JustNature’ for permission to use video footage of a Bombyliidae. Yuka Ito is acknowledged for permission to use her iPhone for fieldwork, and Aria Ito for inspiration for this study.

Author contributions

Conceptualization: S.I., H.R.; Methodology: S.I., H.R.; Validation: S.I., H.R.; Investigation: S.I., H.R.; Data curation: S.I., H.R.; Writing - original draft: S.I.; Writing - review & editing: S.I., H.R., S.N.G.; Visualization: S.I., H.R.; Supervision: H.R., S.N.G.; Project administration: H.R., S.N.G.; Funding acquisition: S.I.

Funding

This work was supported by a Deutscher Akademischer Austauschdienst (DAAD) research grant to S.I. “DAAD Research grant – doctoral programs in Germany” to S.I.

Data availability

The code used in this study is available from: https://biomimetica.notion.site/pollen-catapult-c75f8abd93e44ab7b03f6226b47d4446

Alarcón
,
R.
,
Waser
,
N. M.
and
Ollerton
,
J.
(
2008
).
Year-to-year variation in the topology of a plant-pollinator interaction network
.
Oikos
117
,
1796
-
1807
.
Amador
,
G. J.
and
Hu
,
D. L.
(
2015
).
Cleanliness is next to godliness: mechanisms for staying clean
.
J. Exp. Biol.
218
,
3164
-
3174
.
Amador
,
G. J.
,
Matherne
,
M.
,
Waller
,
D.
,
Mathews
,
M.
,
Gorb
,
S. N.
and
Hu
,
D. L.
(
2017
).
Honey bee hairs and pollenkitt are essential for pollen capture and removal
.
Bioinspir. Biomim.
12
,
026015
.
Barreda
,
V. D.
,
Palazzesi
,
L.
,
Tellería
,
M. C.
,
Olivero
,
E. B.
,
Raine
,
J. I.
and
Forest
,
F.
(
2015
).
Early evolution of the angiosperm clade Asteraceae in the Cretaceous of Antarctica
.
Proc. Natl. Acad. Sci. USA
112
,
10989
-
10994
.
Bowker
,
G. E.
and
Crenshaw
,
H. C.
(
2007a
).
Electrostatic forces in wind-pollination-Part 1: measurement of the electrostatic charge on pollen
.
Atmos. Environ.
41
,
1587
-
1595
.
Bowker
,
G. E.
and
Crenshaw
,
H. C.
(
2007b
).
Electrostatic forces in wind-pollination-Part 2: simulations of pollen capture
.
Atmos. Environ.
41
,
1596
-
1603
.
Bröderbauer
,
D.
,
Weber
,
A.
and
Diaz
,
A.
(
2013
).
The design of trapping devices in pollination traps of the genus Arum (Araceae) is related to insect type
.
Bot. J. Linn. Soc
172
,
385
-
397
.
Brown
,
D.
and
Cox
,
A. J.
(
2009
).
Innovative uses of video analysis
.
Phys. Teach
47
,
145
-
150
.
Claßen-Bockhoff
,
R.
,
Speck
,
T.
,
Tweraser
,
E.
,
Wester
,
P.
,
Thimm
,
S.
and
Reith
,
M.
(
2004
).
The staminal lever mechanism in Salvia L. (Lamiaceae): a key innovation for adaptive radiation?
Org. Divers. Evol
4
,
189
-
205
.
Doi
,
M.
,
Ito
,
M.
and
Auld
,
B. A.
(
2006
).
Growth and reproduction of Hypochoeris radicata L
.
Weed Biol. Manag.
6
,
18
-
24
.
Edwards
,
J.
,
Whitaker
,
D.
,
Klionsky
,
S.
and
Laskowski
,
M.
(
2005
).
A record-breaking pollen catapult
.
Nature
435
,
2005
.
Eisikowitch
,
D.
,
Vaknin
,
Y.
,
Bechar
,
A.
,
Ronen
,
B.
and
Gan-Mor
,
S.
(
2005
).
The role of electrostatic forces in pollination
.
Plant Syst. Evol.
222
,
133
-
142
.
Erbar
,
C.
and
Leins
,
P.
(
1995
).
Portioned pollen release and the syndromes of secondary pollen presentation in the Campanulales-Asterales-complex
.
Flora
190
,
323
-
338
.
Erbar
,
C.
and
Leins
,
P.
(
2015
).
Diversity of styles and mechanisms of secondary pollen presentation in basal Asteraceae-New insights in phylogeny and function
.
Flora
217
,
109
-
130
.
Gurung
,
V.
,
Yuan
,
Y. W.
and
Diggle
,
P. K.
(
2021
).
Comparative analysis of corolla tube development across three closely related Mimulus species with different pollination syndromes
.
Evol. Dev
23
,
244
-
255
.
Hao
,
J. H.
,
Qiang
,
S.
,
Chrobock
,
T.
,
Van Kleunen
,
M.
and
Liu
,
Q. Q.
(
2011
).
A test of baker's law: Breeding systems of invasive species of Asteraceae in China
.
Biol. Invasions
13
,
571
-
580
.
Harris
,
E. M.
(
1995
).
Inflorescence and floral ontogeny in Asteraceae: A synthesis of historical and current concepts
.
Bot. Rev.
61
,
93
.
Huang
,
W.
and
Yu
,
L.
(
2012
).
Serial symmetrical relocation algorithm for the equal sphere packing problem
.
arXiv
1202.4149
. .
Ito
,
S.
and
Gorb
,
S. N.
(
2019a
).
Attachment-based mechanisms underlying capture and release of pollen grains
.
J. R. Soc. Interface
16
,
20190269
.
Ito
,
S.
and
Gorb
,
S. N.
(
2019b
).
Fresh ‘pollen adhesive’ weakens humidity-dependent pollen adhesion
.
ACS Appl. Mater. Interfaces
11
,
24691
-
24698
.
Jankauski
,
M.
,
Ferguson
,
R.
,
Russell
,
A.
and
Buchmann
,
S.
(
2022
).
Structural dynamics of real and modelled Solanum stamens: Implications for pollen ejection by buzzing bees
.
J. R. Soc. Interface
19
,
20220040
.
Jones
,
A. M.
and
Knudsen
,
J. G.
(
1961
).
Drag coefficients at low Reynolds numbers for flow past immersed bodies
.
AIChE J.
7
,
20
-
25
.
Kastinger
,
C.
and
Weber
,
A.
(
2001
).
Bee-flies (Bombylius spp., Bombyliidae, Diptera) and the pollination of flowers
.
Flora
196
,
3
-
25
.
Katzer
,
A. M.
,
Wessinger
,
C. A.
and
Hileman
,
L. C.
(
2019
).
Nectary size is a pollination syndrome trait in Penstemon
.
New Phytol.
223
,
377
-
384
.
King
,
M. J.
and
Buchmann
,
S. L.
(
2011
).
Bumble bee-initiated vibration release mechanism of rhododendron pollen
.
America (NY)
82
,
1407
-
1411
.
King
,
M. J.
,
Buchmann
,
S. L.
and
Spangler
,
H.
(
1996
).
Activity of asynchronous flight muscle from two bee families during sonication (buzzing)
.
J. Exp. Biol
199
,
2317
-
2321
.
Koch
,
L.
,
Lunau
,
K.
and
Wester
,
P.
(
2017
).
To be on the safe site – Ungroomed spots on the bee's body and their importance for pollination
.
PLoS One
12
,
1
-
16
.
Leins
,
P.
and
Erbar
,
C.
(
2006
).
Secondary pollen presentation syndromes of the Asterales — a phylogenetic perspective
.
Bot. Jahrbücher
127
,
83
-
103
.
Lisci
,
M.
,
Cardinali
,
G.
and
Pacini
,
E.
(
1996
).
Pollen dispersal and role of pollenkitt in Mercurialis annua L
.
Euphorbiaceae). Flora
191
,
385
-
391
.
Macior
,
L. W.
(
1964
).
An experimental study of the floral ecology of Dodecatheon meadia
.
Am. J. Bot.
51
,
96
-
108
.
Matherne
,
M.
,
Dowell-Esquivel
,
C.
,
Howington
,
O.
,
Lenaghan
,
O.
,
Steinbach
,
G.
,
Yunker
,
P. J.
and
Hu
,
D. L.
(
2021
).
Biomechanics of pollen pellet removal by the honey bee
.
J. R. Soc. Interface
18
,
20210549
.
Meiri
,
L.
and
Dulberger
,
R.
(
1986
).
Stamen filament structure in the Asteraceae: the anther collar
.
New Phytol.
104
,
693
-
701
.
Michener
,
C. D.
(
1962
).
An interesting method of pollen collecting by bees from flowers with tubular anthers
.
Biol. Trop.
10
,
167
-
175
.
Pacini
,
E.
and
Hesse
,
M.
(
2005
).
Pollenkitt - Its composition, forms and functions
.
Flora
200
,
399
-
415
.
Parker
,
A. J.
,
Tran
,
J. L.
,
Ison
,
J. L.
,
Bai
,
J. D. K.
,
Weis
,
A. E.
and
Thomson
,
J. D.
(
2015
).
Pollen packing affects the function of pollen on corbiculate bees but not non-corbiculate bees
.
Arthropod. Plant. Interact
9
,
197
-
203
.
Picó
,
F. X.
,
Ouborg
,
N. J.
and
Van Groenendael
,
J.
(
2004
).
Influence of selfing and maternal effects on life-cycle traits and dispersal ability in the herb Hypochaeris radicata (Asteraceae)
.
Bot. J. Linn. Soc.
146
,
163
-
170
.
Rajabi
,
H.
,
Shafiei
,
A.
,
Darvizeh
,
A.
,
Gorb
,
S. N.
,
Dürr
,
V.
and
Dirks
,
J.-H.
(
2018
).
Both stiff and compliant: Morphological and biomechanical adaptations of stick insect antennae for tactile exploration
.
J. R. Soc. Interface
15
,
20180246
.
Roberts
,
R. B.
and
Vallespir
,
S. R.
(
2015
).
Specialization of hairs bearing pollen and oil on the legs of bees (Apoidea: Hymenoptera)1
.
Ann. Entomol. Soc. Am.
71
,
619
-
627
.
Snodgrass
,
R.
(
1910
).
The Anatomy of the Honey Bee
.
Cornell University Press
.
Switzer
,
C. M.
,
Combes
,
S. A.
and
Hopkins
,
R.
(
2018
).
Dispensing pollen via catapult: explosive pollen release in mountain laurel (Kalmia latifolia)
.
Am. Nat
191
,
767
-
776
.
Taylor
,
P. E.
,
Card
,
G.
,
House
,
J.
,
Dickinson
,
M. H.
and
Flagan
,
R. C.
(
2006
).
High-speed pollen release in the white mulberry tree, Morus alba L
.
Sex. Plant Reprod.
19
,
19
-
24
.
Thebaud
,
C.
and
Abbott
,
R. J.
(
1995
).
Characterization of invasive Conyza species (Asteraceae) in Europe: quantitative trait and isozyme analysis
.
Am. J. Bot
82
,
360
-
368
.
Thorp
,
A. R. W.
(
2000
).
The collection of pollen by bees
.
Plant Syst. Evol.
222
,
211
-
223
.
Vaknin
,
Y.
,
Gan-Mor
,
S.
,
Bechar
,
A.
,
Ronen
,
B.
and
Eisikowitch
,
D.
(
2001
).
Are flowers morphologically adapted to take advantage of electrostatic forces in pollination?
New Phytol.
152
,
301
-
306
.
Vallejo-Marín
,
M.
(
2019
).
Buzz pollination: studying bee vibrations on flowers
.
New Phytol.
224
,
1068
-
1074
.
Vallejo-Marín
,
M.
,
Pereira Nunes
,
C. E.
and
Russell
,
A. L.
(
2022
).
Anther cones increase pollen release in buzz–pollinated Solanum flowers
.
Evolution
76
,
931
-
945
.
Van Der Kooi
,
C. J.
,
Vallejo-Marín
,
M.
and
Leonhardt
,
S. D.
(
2021
).
Mutualisms and (A)symmetry in plant–pollinator interactions
.
Curr. Biol.
31
,
R91
-
R99
.
Vereecken
,
N. J.
,
Wilson
,
C. A.
,
Hötling
,
S.
,
Schulz
,
S.
,
Banketov
,
S. A.
and
Mardulyn
,
P.
(
2012
).
Pre-adaptations and the evolution of pollination by sexual deception: Cope's rule of specialization revisited
.
Proc. R. Soc. B Biol. Sci.
279
,
4786
-
4794
.
Westerkamp
,
C.
(
1996
).
Pollen in bee-flower relations
.
Bot. Acta
109
,
325
-
332
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information