The mathematics of sexual attraction
© BioMed Central Ltd 2010
Published: 29 March 2010
Skip to main content
© BioMed Central Ltd 2010
Published: 29 March 2010
Pollen tubes follow attractants secreted by the ovules. In a recent paper in BMC Plant Biology, Stewman and colleagues have quantified the parameters of this attraction and used them to calibrate a mathematical model that reproduces the process and enables predictions on the nature of the female attractant and the mechanisms of the male response.
See research article: http://www.biomedcentral.com/1471-2229/10/32
Darwin referred to the rapid success of flowering plants in evolution as an 'abominable mystery'. A great deal of this success relies on the peculiar ways of sexual reproduction that they have evolved. The facts that delight us and spark our curiosity about the sexual behavior of animals, in fact the essence of many novels - attraction, deception, chemistry, male competition, female selection, abortion, death and sacrifice - take place just below our noses every time we appreciate the fine scent of a flower.
Ever since Linnaeus, the process has intrigued biologists, and although the fundamental basis of the interactions were understood in the 19th century by researchers such as Robert Brown, Giovanni Battista Amici, Darwin and Sergei Nawaschin , the molecular nature of this fatal attraction is only now beginning to be understood . Mathematics is now being used to investigate the problem for the first time, in a ground-breaking article recently published in BMC Plant Biology . In this study, Stewman and colleagues use a semi-vivo system to better quantify the nature and range of the ovule's attractants. More importantly, they applied a stochastic mathematical model of whole-cell motility to pollen tube growth. This sort of model has previously been used with success to describe chemotaxis in various eukaryotic cells, such as leukocytes and Listeria. When calibrated with the new experimental data, the model enabled predictions to be made about the size and physical features of the attractant. The results suggest that the growth features observed greatly enhance the attraction efficiency of the ovules.
Most of the beauty that flowers may have for us through their combinations of colors, shapes and scents have, from a developmental point of view, only one purpose: attracting insects and other animals through deception to perform pollination. Reproduction occurs when a pollen grain lands on the receptive surface of pistil, the stigma (Figure 1), whether it arrives there through a sophisticated relationship between insects and orchids or through the mere dispersion of grass pollen by the wind. There, the word 'chemistry' has a literal meaning as well as a metaphorical one, because the bonding between the specialized outer layers of the cell walls of the pollen and the outer stigma cells has been calculated to be stronger than the strongest industrial superglue . Many receptors and ligands come into place to assure proper recognition, and if the match is compatible, the stigma cells nurture the highly dehydrated pollen grains by providing them with water and nutrients, allowing germination. The pollen tube then grows out of the pollen grain in an extreme example of polarized, apical cell growth that results from an unusual set of cell features . The pollen tube is one of the fastest growing cells in nature, and encodes a very specialized transcriptome for cell signaling and communication, which makes them efficient stimulus-perception machines . But how and why these cell-cell interactions take place and navigate the pollen tube so precisely to target the ovule's entrance - the micropyle - are still matters of debate.
These events all take place deep inside the pistil tissues, and it has only recently been possible to image them directly by means of two-photon microscopy . In the decades since the pioneering work of Rosen, Mascarenhas and others, evidence has accumulated that a combination of chemotropic molecules is probably needed for some of the steps. However, for most of the path (from adhesion to entry into the ovary; Figure 1), the anatomical arrangement of the tissues seems to be sufficient to mechanically provide a limited freedom for tube growth. The physical and chemical features of the tissues that make direct contact with pollen tubes (lipids, water, glycosylated proteins, and so on) provide the rest of the signals [1, 6]. With the advent of genetics, a number of screens in Arabidopsis isolated mutations for most steps of this so-called progamic phase of reproduction (a representative list of mutants for each step is shown in red in Figure 1). Some of these mutants have provided evidence for long-range targeting, with signaling occurring over perhaps as far as 500 μm. The list of genes involved is likely to get significantly larger as new screening methods reveal dozens of male- and female-specific mutations, especially in Arabidopsis , but so far none has brought us close to the specific chemotaxis molecule(s).
The development of a suitable semi-vivo system of fertilization for Arabidopsis (Figure 2b, c) has not only enabled the use of genetics, it has also been used to investigate the physical nature of the diffusing molecule in Arabidopsis, as it has in Torenia . Stewman et al.  have now carefully analyzed the growth curvature angles of Torenia pollen tubes (Figure 2d) in various experimental conditions, namely with different incubation times of the ovules (presumably corresponding to different levels of a standing gradient), and could thus determine important quantitative parameters of the attraction process. First, they found that the gradient action could extend a distance of 100 to 150 mm , a distance significantly longer than had been previously thought [6, 9] or than had been experimentally tested with artificial gradients of isolated molecules [2, 12]. This distance estimate probably means that various types of molecules with different ranges and actions come together to produce the biological reaction. There is no doubt that species specificity must be assigned by gene-encoded proteins, and that the isolated LUREs and the Zea mays egg apparatus 1 protein (ZmEA1) both produce positive tropic effects. However, the evidence seems to implicate many other non-specific, small, diffusible chemicals, such as ions or even the signaling gas nitric oxide (NO) . In fact, the involvement of NO is not a surprise, because various behaviors shown by pollen tubes when NO concentration is disturbed seem to indicate a slowing down of growth of the pollen tube as it approaches the diffusion source; these effects fit with those found in the model of Stewman and colleagues .
The apparent simplicity of pollen tube growth is leading to various approaches to mathematical modeling of its main features. Pollen tubes have been successfully described through mechanical or geometrical approaches . Looking from a different angle, the prominence of the signaling systems of pollen tube growth was recently modeled by assuming that soluble N-ethylmaleimide sensitive factor attachment protein receptors (SNAREs) and small signaling GTPases are the main effectors of growth .
Stewman et al.  took a step further and focused on the formal properties of the system as shown from the kinetic parameters precisely derived from the semi-vivo system. The problem of gradient sensing is not trivial, and the authors  followed the strategy of building on stochastic whole-cell motility models, which basically assume that pollen tubes can sense a difference in the fraction of the receptors bound to an attractant and change their growth depending on this fraction. One prediction of the model was that a slower growth rate inside a standing gradient of attractant (which is assumed to exist from the isotropic diffusion of the attractant from the ovule) would greatly enhance the ability of pollen tubes to successfully target ovules, a prediction that they could statistically validate with their experiments . Furthermore, the model describes the observed patterns of random and directed growth observed during growth of pollen tubes in vitro.
Although the assumption of a purely formal mathematical description of any biological phenomenon may lead to purely phenomenological descriptions of limited experimental value, the complexity of the system and the number of different cellular components that seem to be fundamental for growth  make these approaches highly informative, or even fundamental, to understanding the mechanistic basis of the macroscopic response of the system. For example, this model  may help us to understand how pollen tubes couple external guidance cues with intracellular ion gradients or other known cell steering mechanisms. In addition, the model assumes that there are at least two patches of receptors, which are separated in the pollen tube. Although this assumption says nothing about what the receptors are, it does suggest that a minimal model of a sensor with a strict apical point location would not work; instead, it favors a membrane or cytosolic spatially segregate or a receptor that stretches across (at least) the diameter of the tube. The model also assumes that the change in concentration across the tip of the tube is much less than the average concentration at the tip. Despite being affected by the limitation of the two-dimensional modeling performed , this is a strong quantitative prediction that can be validated when searching for novel molecules that might fit the profile. And how does slowing down of the tubes take place near the micropyle? NO, for example, is known for slowing down growth rates , and the formulation of the model does, in fact, allow for multi-factorial interactions affecting the growth process. Importantly, this mathematical formulation may allow discrimination of different effects during experimental procedures based both on the deviation of the growth angle and/or the relationship of the growth rate to successful targeting.
As with many other mathematical approaches to complex biological behavior, this new model from Stewman et al.  raises more questions than answers. But the fact that new approaches are contributing to a precise experimental description of the system [2, 11, 12] may make mathematical modeling an important tool for testing and selecting candidate molecules that may fit the in vivo biological profile of the final step of plant sexual attraction.
I thank Leonor Boavida for drawing and compiling Figure 1. Financial support for my laboratory is provided by FCT grants PTDC/BIA-BCM/108044/2008 and PTDC/QUI/64339/06.