Skip Navigation

This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (1)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by HARRISON, L. G.
Right arrow Articles by VON ADERKAS, P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by HARRISON, L. G.
Right arrow Articles by VON ADERKAS, P.
Agricola
Right arrow Articles by HARRISON, L. G.
Right arrow Articles by VON ADERKAS, P.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

Annals of Botany 93: 423-434, 2004
© 2004 Annals of Botany Company

Spatially Quantitative Control of the Number of Cotyledons in a Clonal Population of Somatic Embryos of Hybrid Larch Larix x leptoeuropaea

LIONEL G. HARRISON*,1 and PATRICK VON ADERKAS2

1 Department of Chemistry, University of British Columbia, 2036 Main Mall, Vancouver, BC, Canada V6T 1Z1 and 2 Graduate Centre for Forest Biology, Department of Biology, University of Victoria, PO Box 1700 STN CSC, Victoria, BC, Canada V8W 3N5

* For correspondence. E-mail lionel{at}chem.ubc.ca

Received: 9 June 2003; Returned for revision: 30 September 2003; Accepted: 17 December 2003


   ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 LITERATURE CITED
 

Background and Aims Many conifer embryos, both in natural seeds and in clonal populations of somatic embryos, display variability in the number of cotyledons. In hybrid larch, Larix x leptoeuropaea (synonymous with L. x marschlinsii Coaz), such variability has previously been reported in somatic embryos, together with a decrease in the average cotyledon number when benzyladenine (BA) is applied exogenously. Described here is a spatially quantitative study with the aim of throwing some light on the way cotyledon number is determined, and hence the mechanism of cotyledon formation.

Methods Stock cultures of embryogenic tissue were maintained and later made embryogenically active by standard methods. Development through cotyledon formation was followed by optical microscopy with quantitative measurement of embryo diameter and number of cotyledons. SEMs of representative stages and cotyledon numbers were done for purposes of illustration in this account. Existing mathematics of waveforms on a disc were cast into a form suitable to compare with the quantitative data.

Key Results The number of cotyledons is linearly related to the diameter of the apical surface of the embryo (which approximates a circular disc) at the time of first appearance of the cotyledon primordia. This linearity is a constant-spacing phenomenon between adjacent primordia. Addition of BA to the medium restricts the range of apical diameters without changing inter-cotyledon spacing. Slope/intercept ratio of the linear plot matches expectation for initiation of cotyledon pattern as a harmonic waveform on a circular disc.

Conclusions The entire pattern of cotyledon primordia arises as a single entity coordinated by a mechanism with wave-forming properties. This is explicable by diverse mechanisms, especially either mechanical buckling (‘biophysical’) or reaction–diffusion kinetics (‘physicochemical’).

Key words: Larix x leptoeuropaea, larch, cotyledon formation, somatic embryos, pattern formation, spatial quantitation, disc waveforms.


   INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 LITERATURE CITED
 
Strategy of the project
We report here a study on quantitative spatial aspects of the initiation of whorled sets of cotyledons on the flat discoid apex of somatic embryos of a hybrid larch Larix x leptoeuropaea. Our strategy builds on work done over two decades in the group of one of the authors (Harrison et al., 1981, 1988, 1997; Harrison and Hillier, 1985) upon formation of whorled sets of hairs on the flattened growing tip of an unusually large unicellular alga, Acetabularia acetabulum (Dasycladales, Chlorophyta). Therefore, we are looking for features over-riding the distinction between unicellular and multicellular apices. We proceed expecting an ultimately affirmative response to the question: ‘When a definitive identification of the whorl-forming event is reached, will it throw light on the development of anything other than the dasyclads?’ And we adopt ‘the pattern formation viewpoint’, which envisages that ‘one must first establish such a morphogenetic sequence by identifying individual pattern-forming events leading to the final morphology.’ (quotations: Dumais and Harrison, 2000)

Our approach, at this stage, is somewhat independent of (but we hope ultimately complementary to) approaches that have led to current controversy (Barton and Poethig, 1993; Kaplan and Cooke, 1997; Bowman and Eshed, 2000) on whether cotyledons arise out of the shoot apical meristem. We look for physicochemical concepts that may form a unifying bridge between patterning mechanisms on multicellular apices with numerous spatially patterned gene activities, and those of a growing tip from which the single nucleus is some centimetres distant.

The conspicuous similarities between the phenomena studied in A. acetabulum and L. x leptoeuropaea are as follows.

1. In both cases, the morphogenetic apex of any one specimen can be observed from day to day or hour to hour, because the unicell is a free-living marine organism and the somatic embryos grow from the surface of a mass of embryogenic tissue and are not concealed within megagametophytes. The morphogenetic events can be observed through the first appearance of whorled organs. Measurements can be made very close in time to the pattern-forming event.

2. In both cases, the entire set of whorled structures arises strictly simultaneously, unlike many later-developing sets of whorled structures in plants. [For L. x leptoeuropaea, we see no sign of a sequential process, or relationship to a spiral, as mentioned by Buchholz (1919).] Thus we count the emergence of all the organ primordia in a whorl as one event.

3. In both cases, the apex flattens just before the new organs appear. This allows measurement of a well-defined apical diameter.

4. In both cases, a larger apex gives rise to a whorl with more organs. This is very far from being a trivial observation. Both quantitative studies have shown that the apex diameter at the time of the whorl-pattern-forming event is very sloppily controlled, but the spacing between adjacent organ primordia is rather precisely controlled at the time of their first morphological appearance. This focuses attention on the ability of a developing organism to set up quantitative scales of distance. How do the physics and chemistry of life achieve this?

5. In both cases, the circle of centres of the organ primordia is inset from the edge of the discoid apex. In terms of the count of events, our inferences from the two cases diverge at this point. For A. acetabulum, two events were postulated to determine the inset and the spacing between adjacent organs. But the patterns of L. x leptoeuropaea fit very well, quantitatively, with still snapshots of the modes of vibration of a circular disc. This is strong circumstantial evidence that both the cotyledon primordia and the flat central region (where the shoot apical meristem will emerge) are at this time parts of a single patterning entity.

Variability in cotyledon number in the conifers
Polycotyly occurs naturally in angiosperms (e.g. species of Psitticanthus, Loranthaceae, have from three to nine cotyledons; Kuijt, 1967) but is far commoner among gymnosperms, in particular conifers; typically the initially circularly symmetric embryo is still rotationally symmetrical when capped by a whorl of cotyledons (Juguet, 2002).

The control of cotyledon number has been variously ascribed to embryo size, genetic heritability, phenotypic variability and gene expression. Although the number of cotyledons per whorl is generally characteristic of a particular taxon, it is not necessarily fixed (Butts and Buchholz, 1940). Embryo size has been observed to vary from year to year, and with it cotyledon number. In fir, Abies procera (Sorensen and Franklin, 1977), and in Cola acuminata (Oladokun, 1982), embryos were larger and had more cotyledons in good rather than poor seed years. In Pinus banksiana, cotyledon number has been shown to be heritable (Briand et al., 1998). Phenotypic variation can be brought about by altering hormone metabolism. In dicotyledonous plants, polycotyly is found in both pin1, a mutant affecting auxin transport (Aida et al., 2002), and amp1, a mutant that produces six times the normal endogenous cytokinin level (Chaudhury et al., 1993). Somatic embryo development is regulated by timed applications of exogenous plant growth regulating substances (PGRs). Somatic embryos show far greater variability in the number of cotyledons than their zygotic counterparts, and the variability can be controlled by PGRs, e.g. in vitro, by exogenous application of anti-auxin 2,3,5-triiodobenzoic zcid (Choi et al., 1997) and benzyladenine (von Aderkas, 2002). Genetic analysis of lateral organ formation during embryogenesis has been restricted to a few model dicotyledons. Polycotyly is found in some mutants, e.g. fass (Torres-Ruiz and Jürgens, 1994) and pin1 (Aida et al., 2002). Other mutations affect both the number of cotyledons initiated and the rate of leaf initiation, e.g. amp1, xtc (Conway and Poethig, 1997), häuptling and cop2 (Lehman et al., 1996). To date no cotyledon-specific genes have been isolated, reinforcing the view that events at the apex of the developing embryo involve progressive regulation and interplay of various sets of genes (Aida et al., 1997). Lateral organ formation, from cotyledons to leaf primordia, is thought to be coordinated by spatial and temporal gene expression of WUS, ANT, CUCL1, CUCL2, STM, PIN1 and CLAV (for review, see Bowman and Eshed, 2000). This model is not without controversy (Kaplan and Cooke, 1997). This centres on the timing of initiation of a shoot apical meristem capable of initiating leaves, and involves the question as to whether leaves and cotyledons are homologous. Our material and experimental design are not suited to addressing these questions. We focus on embryo size and cotyledon number to identify and characterize pattern-forming events.

In conifer embryos, the events preceding the appearance of cotyledons are as follows. The globular embryo grows so as to flatten its apex to a circular disc. This stage has been called ‘shoulders’ (Spurr, 1949) or ‘ring meristem’ (Juguet, 2002; von Aderkas, 2002). (We avoid ‘shoulders’ because it relates to a side view or longitudinal section. The important perspective for our study is the plan view of the disc.) Cotyledon initials arise from this flat disc (Spurr, 1949). There has been some discussion as to whether final numbers of cotyledons have been affected by fusion or splitting of cotyledon initials (Buchholz, 1919; Juguet, 1992). These studies were not of initial spacing and numbers of cotyledons. Typically, they used sampling from various stages, not direct observation of single embryos throughout development. Study of initiation in individual embryos is clearly required to characterize the pattern-forming event without later modifications.

We chose to study somatic embryos of hybrid larch (Larix x marschlinsii Coaz = L. decidua x L. leptolepis = L. x leptoeuropaea) because one of us had recently characterized the distributions of cotyledon number in maturation media without benzyladenine (BA) and with several concentration levels of BA, which changed the mean cotyledon number and distributions about the mean (von Aderkas, 2002).


   MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 LITERATURE CITED
 
Plant material
Larix x leptoeuropaea (clone 69-18) used in all experiments was provided by Marie-Anne Lelu. It had been initiated by secondary embryogenesis (Lelu et al., 1994). Larch somatic embryos were multiplied on maintenance medium and matured by transferring cultures to a medium that promoted maturation.

Maintenance.
Stock cultures were maintained on gelled, modified Murashige and Skoog (MS) medium known as MSG (Becwar et al., 1990), which was composed of MS salts supplemented with 60 mM sucrose, 9 µM 2,4-dichlorophenoxyacetic acid (2,4-D), 2·5 µM benzylamino (BA) and Gelrite (4 g L–1). Prior to autoclaving, the pH was adjusted to 5·8. Following autoclaving, 10 mM glutamine (filter-sterilized) was added. Every 14 d, the tissue was subcultured to fresh medium. Petri dishes containing the cultures were kept in the dark at 22·4 °C.

Embryogenesis.
Clumps of embryogenic tissue were removed from the maintenance medium and transferred to MSG medium supplemented with 200 mM sucrose and 4 g L–1 Gelrite, to which either 0 or 4·4 µM BAP had been added.

Microscopy, observation and measurements
Observation, photography and measurement of developing embryos in vivo were done with an Olympus SZ-Tr dissecting microscope furnished with PM-10-M camera and OSM 213120 eyepiece micrometer. As embryos appeared on the embryogenic tissue masses, they were selected individually for observation and measurements. Each embryo was followed regularly, usually day by day. A staging of the developmental sequence was devised, with particular attention to the stage at which circular symmetry was first seen morphologically to be broken. The main objective was to measure a clearly defined embryo diameter at a time as close as possible to the symmetry-breaking event.

For scanning electron microscopy, representative specimens were fixed for 1 d in 2·5 % glutaraldehyde in 0·05 M phosphate buffer. The embryos were washed in buffer three times (1 h per wash), then dehydrated in a series of ethanol solutions (30, 50, 70, 95, 95 %) for 30 min each, then left overnight in 95 % ethanol. The next morning they were transferred to 100 % ethanol before being critical-point dried (Bomar SPC-1500 Critical Point Dryer ), mounted on stubs, sputtercoated with gold (Edwards S150B Sputter Coater) and viewed on a Hitachi S-3500N scanning electron microscope in the Electron Microscope Unit of the University of Victoria.


   RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 LITERATURE CITED
 
Variability in the number of cotyledons
Figure 1 illustrates the variability in number (nc) of cotyledons, with examples from nc = 2 to nc = 8 shown at a fairly advanced stage of development of the cotyledons. For embryos grown in one author’s laboratory (L.G.H.) for quantitative measurements, Fig. 2 shows a histogram of the distribution of nc. This is almost identical to the histogram reported for embryos grown in the second author’s laboratory (PvonA) in von Aderkas (2002).



View larger version (161K):
[in this window]
[in a new window]
 
Fig. 1. SEMs of developing somatic embryos to show variability of cotyledon number nc. Examples with nc = 2, 3, 4, 5, 6 and 8 are shown in parts A–F. Scale bars = 250 µm (A–C); 500 µm (D); 200 µm (E and F). These samples are at a much later stage of development than that to which the d2 measurements in Fig. 7 refer, and are on average about 1·6 times larger in diameter than the Fig. 7 data. Part A shows that nc = 2 in fact occurs, though there happened to be no such samples in the 44 embryos measured for Fig. 7A.

 


View larger version (12K):
[in this window]
[in a new window]
 
Fig. 2. Distribution of number of cotyledons for embryos grown in control medium. The numbers here are from the data of Harrison (see text), as used also in Fig. 7A.

 
Staging adopted to locate symmetry-breaking events in the developmental sequence
Figure 3 illustrates the sequence of development of an embryo, with a notation for stages as explained in Fig. 4. For SEMs, a different embryo had to be used for each picture. In Fig. 3, examples have been shown to illustrate as closely as possible what the sequence would look like for one embryo. Hence from Fig. 3C onwards, all examples were chosen for nc = 3. Our in vivo observations under the light microscope led us to the following concept of a ‘normal’ sequence. The embryo starts as a shallow circular dome, advances to a shape close to hemispherical (sometimes topping a short cylinder) and then conspicuously flattens (Fig. 3B) without losing circular symmetry about its axis of growth. The cotyledons then appear as outgrowths from the flattened top surface (Fig. 3C). It seems most reasonable to suppose that the number, nc, of cotyledons has been established by a patterning event between the stages shown in parts B and C in Fig. 3. Usually, measurements of the spatial dimensions of the embryos have been made by following an embryo from day to day, making measurements each day. For determination of the spacing between cotyledons at the moment when their number is established, measurements were selected for stages like those shown in B and C in Fig. 3.



View larger version (151K):
[in this window]
[in a new window]
 
Fig. 3. SEMs of embryos to illustrate stages with terminology defined in Fig. 4. From (C) onwards, the examples shown are all for embryos forming three cotyledons. All scale bars = 250 µm. Stages: Do(H) (A); An(F) (B); Co(N) (C–F), from earliest at which cotyledons can be counted to advanced morphological development of cotyledons.

 


View larger version (29K):
[in this window]
[in a new window]
 
Fig. 4. Staging of somatic embryo development used in this work. The terms ‘annulus’ and ‘annular’ are used here for a shape which has lost hemispherical symmetry, while retaining circular symmetry, by flattening of the top. This has some implication that flattening arises from a shift of growth rate maximum from the pole of the dome to a ring near to its equator. Not every embryo goes through all the stages listed in sequence. In particular, the cratered morphology is often an aberrant form that fails to produce cotyledons. The normal sequence is direct from An(F) to Co(I). Change in cotyledon number [i.e. Co(N) to Co(Re)] happens in a minority of cases; but it indicates that cotyledons must be counted as early as possible for correlation with the symmetry-breaking event at Co(I).

 
Figure 4 lists the staging terminology adopted in this work. We have avoided some conventional terms such as ‘buttress’ or ‘shoulder’ in favour of terms that indicate the early circular symmetry of the embryo about its axis of growth and initiation of cotyledons out of a shape that otherwise retains this circularity.

Figure 5 shows some abnormalities displayed in the development of a minority of embryos. Figure 5A and B shows the number, nc, of cotyledons changing either by fusion of two cotyledons or splitting of one into two (usually the latter, from the time-sequence observations of LGH). These pictures illustrate that nc should be measured at the earliest possible stage of appearance of the cotyledons, if nc is to be correlated with the dimensions of the embryo at the moment of pattern formation. Figure 5C shows that the cotyledon primordia are not necessarily in contact with each other. They can be spaced out with unchanged regions of the flat top circular surface between them. Figure 5D shows the phenomena we call polygonization and facetting, in which, before appearance of cotyledons, the embryo loses circular symmetry to become a frustum of a pyramid, therefore with a polygonal top surface (in this case, five facets and pentagonal top). Figure 5E and F shows cratering of the top surface before cotyledons appeared. The crater sometimes had a rounded rim and sometimes a sharp one. In the latter case, the embryo most commonly stopped developing before cotyledons formed.



View larger version (157K):
[in this window]
[in a new window]
 
Fig. 5. SEMs (except D which is from light microscopy) illustrating a few abnormalities that were observed in the developmental sequences. (A and B) Examples of cotyledons only partly divided, either one splitting into two or two fusing into one. (C) Cotyledons growing not close to contact with each other. (D) Breaking of circular symmetry by polygonization of the top and facetting of the sides, in this case turning the embryo into a frustum of a pentagonal pyramid. This symmetry-breaking did not in each case correspond exactly to the number of cotyledons afterwards initiated, but did on average (Table 2). (E and F) Examples of crater formation on the top surface, with a rounded rim in (E) and a sharp rim in (F). The latter shape usually does not go on to form cotyledons. Scale bars = 250 µm (A); 500 µm (B, C and F); 200 µm (D and E).

 
Spatially quantitative measurements and the intercotyledon spacing {lambda}
Measurements were made to determine whether or not the spacing between adjacent cotyledons was constant, as was found for the hairs of Acetabularia vegetative whorls (Harrison et al., 1981). It is not clear, a priori, whether one should expect a constant-spacing phenomenon, whatever its mechanism, to be established round the edge of the flattened circular top of the embryo (measured diameter d2) or around a smaller circle on that top surface. Figure 6 illustrates the terminology used in presenting the results. In some of the experiments, a diameter, d1, was also measured at the base of the embryo; it did not correlate significantly with the number of cotyledons. The degree of tapering of the embryo between base and flattened top is quite variable. Some embryos are almost precisely cylindrical.



View larger version (12K):
[in this window]
[in a new window]
 
Fig. 6. Definition of the measured quantities d2, {lambda} and c/2; see also legend to Fig. 7.

 
To look for a constant-spacing effect in a way that does not prejudge its precise location or its mechanism, we tested for a linear relationship between the diameter d2 of the flattened top and the number nc of cotyledons:

d2 = ({lambda}/{pi})nc + c(1)

If such a relationship is found, the line may or may not go through the origin. If it does, the constant-spacing event is around the outside edge of the top surface. If there is a positive intercept at nc = 0, then the constant spacing {lambda} is established around a smaller circle, at a distance c/2 from the outside edge. This phenomenological description should not be read as assigning a special mechanistic role to events located at that inner circle; see the section ‘How many developmental events?’ in the Discussion.

In the present work, such a linear relationship was found (Fig. 7), both for embryos grown in control medium and for those grown in medium with 4·4 µm BA. Table 1 gives the values of the parameters in eqn (1) from the plots in Fig. 7. The spacing {lambda} is not significantly affected by the concentration of BA used, although the addition of BA certainly reduces the average number of cotyledons, as observed in this work and in von Aderkas (2002). The reason for this change is a change in the distribution of embryo diameters d2, as shown in Fig. 8. Essentially, the addition of BA greatly reduces the fraction of the embryos with d2 above about 300 µm.



View larger version (8K):
[in this window]
[in a new window]
 
Fig. 7. Plots of mean diameter, d2, vs. number, nc, of cotyledons for embryos grown at 22·4 °C with: control medium (A); medium with 4·4 µm BAP (B). For (A) the numbers of embryos for which measurements have been averaged in each data point are given by the heights of the bars in Fig. 2. In each graph, the solid line is a visually judged best fit to the points. In (B) the open circle is a value from only one embryo and has been ignored as an outlier in drawing the line. The dashed lines have a ratio slope/intercept = 0·226, corresponding to a fit of harmonic waveforms on a disc to its diameter with a particular boundary condition, as described in the Appendix. Those lines could be moved arbitrarily up or down, provided that slope and intercept are changed by the same factor.

 

View this table:
[in this window]
[in a new window]
 
Table 1. Spatial quantitation of the cotyledon-initiating event from the plots in Fig. 7
 


View larger version (17K):
[in this window]
[in a new window]
 
Fig. 8. Distributions of diameter d2 for embryos grown in: control medium (open bars) or medium with 4·4 µm BA (black bars). These show that the mean number of cotyledons found for the BAP medium is smaller than that for the control medium because the range of embryo diameters is smaller for the BA medium, while the spacing {lambda} (Fig. 7) is very little affected by BA.

 
For the minority of embryos that showed breaking of circular symmetry by polygonization, sectoring or facetting before appearance of cotyledons, the number of cotyledons eventually formed was not always exactly the same as the number of sides of the polygon, etc., but there was a fairly good correlation of average numbers (Table 2).


View this table:
[in this window]
[in a new window]
 
Table 2. For embryos that broke symmetry by forming polygonal or sectored top surface, relation of number of sides to number of cotyledons formed later
 


   DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 LITERATURE CITED
 
How many developmental events?
Three features need explanation in terms of pattern-forming theory. (1) The hemispherical embryo, stage Do(H), flattens to give a top surface in the form of a flat disc, stage An(F), with a well-defined diameter, d2. (2) Cotyledons form with a constant spacing between adjacent ones. (3) That spacing is not around any visible morphological feature of the embryo, but around a circle of smaller radius, inset a fixed (i.e. independent of nc) distance from the edge of the disc. How many developmental events need to be postulated to account for all these features: one, two or three? Chronologically, (1) happens first while (2) and (3) are simultaneous. In the previously studied case of Acetabularia, one of us has maintained (e.g. Harrison et al., 1981, 1988, 2001) that a two-stage mechanism is necessary for a system of many repeated structures (most commonly 7–22) to arise in a whorled arrangement. The first stage defines an annulus around which the second stage generates a repeating pattern. This kind of theory ascribes the origin of the inset c/2 and the spacing {lambda} to two separate events.

But for small numbers of repeats in a pattern, two stages may not be necessary. For dichotomous branching, a single-stage mechanism generating wave patterns on a hemispherical tip (surface spherical harmonics) gives a good account (Harrison et al., 2001). The present example, however, with branchings ranging from one to eight parts, spans the intermediate range between small and large numbers of parts in a pattern and directly addresses the question of how complex a pattern must be for a one-stage mechanism to be inadequate.

We take the flattening of the tip as having occurred and consider features (2) and (3), the patternings on a flat disc, together. If these are produced by any kind of wave-forming mechanism, the appropriate waveforms to consider are the patterns of vibration of a circular drumskin attached to a rigid outer border. These combine circular ripples with pie-slice patterns. A priori, it is a promising possibility that the mathematics of such shapes will explain simultaneously the circle at inset c/2 and the repeated maxima around that circle at spacing {lambda}. A brief description of how these waveforms will fit patterns of increasing numbers of whorled structures to increasingly large discs is given in the Appendix. This shows that a plot of d2 vs. nc should indeed be a straight line with a positive intercept. The ratio of slope to intercept is a significant quantity. The fit of one complete up-and-down ripple to the radius of the disc (second zero of the Bessel function at the edge) gives the ratio slope/intercept 0·226 (broken lines, Fig. 7A and B), compared with 0·19 and 0·21 for the experimental data (points and solid lines).

From this analysis, the answer to ‘how many developmental events?’ is that inset c/2 and spacing {lambda} can be determined by any one mechanism that is able to form two-dimensional wave patterns on a circular disc. This could be, for example, chemical reaction–diffusion (Lacalli, 1981), mechanical buckling (Green, 1999) or a more complex mechano-chemical interaction. The fit of our data to waveforms on a disc is concordant with any of these types of mechanism, and strengthens the concept that the pattern of the emerging cotyledons is a single macroscopic entity governed by interactions coordinated on the spatial scale of the whole disc.

This view is supported by the absence of obvious structural features linking or separating the cotyledons in the earliest stages of their development. Gutmann et al. (1996) show longitudinal sections of somatic embryos of L. x leptoeuropaea; the embryonal mass consists of one cell type with somewhat irregular cell size and distribution. Indeed, the cotyledons are more regular in size and shape than the detailed cellular structure of the tissue from which they arise, further supporting the view that the organization of a plant is organismal and transcends the order of the cells (Kaplan and Hagemann, 1991; Cooke and Lu, 1992).

As to feature (1), the transformation of a hemisphere into a flat-topped shape: this requires slow growth at the centre and fastest growth in the outermost regions, without loss of circular symmetry. We think that this is probably initiated by a separate patterning process from the later cotyledon-forming one. In that case, the three listed features need two patterning mechanisms to explain them. But, unlike the Acetabularia patterning, the first process leading to tip flattening is not being called upon to generate the inset circle around which the cotyledons are spaced. It has only to turn a hemisphere into a disc, and leave the rest to the second process.

Significance of a constant-spacing phenomenon
A constant spacing between adjacent repeats of an organ when it first becomes identifiable is a very significant phenomenon, indicating that the organism is able to set up a quantitative scale of distance. This is a common feature of biological development, but far from easy to explain definitively, although theories having this property have been known for at least half a century (Turing, 1952). Up to the present, constant-spacing effects have been explained by reaction–diffusion theory but not yet by mechanical or mechano-chemical theory. Morphogen distributions on a disc (Fig. 9), computed for Brusselator-type reaction–diffusion dynamics (Lacalli, 1981), may be compared with expressed cotyledon patterns in Figs 1E and 5A.



View larger version (48K):
[in this window]
[in a new window]
 
Fig. 9. Two computations of spatial distributions of a reaction-diffusion morphogen on a disc (from Lacalli, 1981, with permission). (A) An n = 6 pattern. Compare it with Fig. 1E. (B) A pattern with secondary branching. Compare it with Fig. 5A.

 
That rates of reaction and diffusion can set up a measure of distance (e.g. the wavelength of a pattern) may be shown by dimensional analysis. Reaction rates provide a measure of time; a first-order rate constant (k) is in (time)–1, usual units s–1. Diffusion gives a proportionality between distance and (time)1/2; diffusivity (D) has dimensions (distance)2 (time)–1, e.g. cm2 s–1. Diffusivity and rate constant combined in the form (D/k)1/2 give (cm2 s–1/s–1)1/2 = cm. Reaction–diffusion theories can be mathematically complex; but they lead to wavelength expressions that are always variations on the theme (D/k)1/2.

Thus in the Brusselator (Prigogine and Lefever, 1968), a promising mechanism for branching processes in plants (e.g. Holloway and Harrison, 1999), reactants A and B give rise to diffusible intermediates X and Y. The kinetics (Harrison, 1993) cause X and Y to be patterned with wavelength {lambda} approximately:

{lambda} = 2{pi}(DX DY/a2 bcd–2 A4 B)1/4(2)

a, b, c and d are rate constants (of first, second and third orders); A and B are the concentrations of A and B. The denominator has units s–2, and the dimensional analysis is like that of (D/k)1/2, giving {lambda} in cm.

What kinds of substances may be involved in patterning?
Diverse classes of biochemicals might be involved in reaction–diffusion mechanisms. X and Y could be microRNAs, transcriptional regulators or other proteins, or the smallest molecules and ions. In Acetabularia, an effect of calcium concentration on {lambda} led Goodwin to propose calcium as a morphogen, while Harrison indicated it as an activator of a protein to become morphogen precursor A [see Chapter 10 of Harrison (1993) for this and other examples]. The diversity increases with the possibility of mechanical properties playing a role in patterning. We expect ultimate linkage to genes, but the link may be a very long bridge upon which our current work is at the other end. We have studied a clonal population, so that aberrant development does not disclose mutations, e.g. the ‘crater’ morphology is not related to the ‘cup-shaped cotyledon’ genes CUC1 or 2 (Aida et al., 1997). More probably, it is part of a statistical distribution of different patterns made by the same dynamic mechanism (Harrison et al., 2001; Nagata et al., 2003). In a very general sense, however, genes first expressed as chemical patterns, such as annuli, in any branching process in plants, are of potential interest. These include a few flower development genes intermediate between meristem identity and organ identity genes: FIMBRIATA (FIM) in Antirrhinum (Simon et al., 1994); UNUSUAL FLORAL ORGANS (UFO) in Arabidopsis (Ingram et al., 1995); and BLH2, a homeodomain gene in Arabidopsis, expressed in various tissues at different times in development, including expression in the floral apex in an annular pattern where the sepal whorl is about to form (M. Pidcowich, D. Godt, K. Kushalappa and G. Haughn, pers. comm.). More plan-view studies of gene expressions on apices would help in bridging the gap between genetic work and our macroscopic work. For embryos, current controversy on the relationship of lateral organ formation to the initiation of a shoot apical meristem connects to our present work only in that we adduce evidence that the cotyledon-forming region and the central region where the apical meristem will arise participate in the same patterning mechanism at the time of cotyledon initiation.

In biochemistry remote from genetics but likely to be close to developmental dynamics, significant classes of substances include: (1) for mechanical buckling, cell wall constituents that govern mechanical properties; (2) for mechano-chemistry, the cytoskeleton, e.g. actin and actin-associated proteins; (3) for reaction–diffusion, integral membrane proteins such as those heading signal-transduction cascades (for possible relation of these to Brusselator A and X, see Harrison et al., 1997), and also enzymes affecting rate constants, e.g. Brusselator a, b, c and d; (4) various types of plant growth regulators (PGRs). We had hoped that BA would turn out to be a {lambda}-controller; but Figs 7 and 8 show that it reduces the average number of cotyledons by restricting the size range of the embryos, without changing {lambda}. Other PGRs should, however, be tried.


   ACKNOWLEDGEMENTS
 
We thank the following: for culture maintenance, the UBC Chemistry Department’s Biological Service, Elena Polishchuk, Director, and Susanna Ip, and at UVic, Andrea Coulter, Stephen Burgess, Lindsay White and Heather Milligan; for assistance with the SEM, Dr C. Singla; for comments on the manuscript of this paper, David M. Holloway, and for research grants to both of us, NSERC Canada.


   APPENDIX
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 LITERATURE CITED
 
Matching the data for d2 vs. nc to Bessel function waveforms on a disc
[For a clear introductory text on classical wave theory in one, two and three dimensions, see Chapter 3 of Kauzmann (1957); for application of Bessel functions to morphogenetic patterns in unicellular algae, see Lacalli (1981).] These are essentially products (in the mathematical sense) of undulating radial functions Jn (r), having circular maxima, minima and zeros (nodal circles) with simple sine wave functions of angle {theta} around the disc:

U(r,{theta}) = A Jn (xnl r /a)sin(n{theta})(A1)

where U is the displacement that has wave character, e.g. actual up-or-down distance for a vibrating drumskin, or increase or decrease of a chemical concentration from a spatially uniform value for reaction–diffusion dynamics; r,{theta} are circular polar coordinates with origin at the centre of the disc; A is a time-dependent amplitude, oscillatory for a drumskin vibration, growing monotonically in time to a final steady value for a reaction–diffusion pattern; a is the outer radius of the disc, corresponding to d2/2; xnl r/a = x is the independent variable for the Bessel function, Jn, representing radius scaled so that the radius a corresponds to the lth zero of the function of order n (Fig. 10A); and the angular function has a number of maxima, n, around the circle corresponding to nc for the cotyledon pattern.



View larger version (13K):
[in this window]
[in a new window]
 
Fig. 10. (A) The radial part, Bessel function Jn (x), of a waveform on a disc, for n = 2. (B– E) For n = 1, 2, 3 and 10, Jn functions truncated at the second zero (l = 2), with broken lines to stress the approximate constancy in width of the patterned region as the disc becomes larger. (F) Plot of the second-zero radius versus n (or nc for the postulated correlation with our data). The broken lines in Fig. 7 are drawn with the same ratio of slope to intercept as this plot. Rotate (F) 90° anticlockwise for comparison with Fig. 7. As shown here, comparison of (F) with (B–E) shows that the intercept c/2 corresponds roughly to the distance from peak to boundary of the Bessel function. Lettering B, C, D, E on the points plotted in (F) corresponds to the plots thus lettered above.

 
To clarify the nature of the function Jn, the first cycle of this function is plotted in Fig. 10 for n = 1, 2, 3 and 10. If one chooses the boundary condition that the edge of the disc is at the second zero (l = 2) of Jn, each plot ends at the disc outer radius. Evidently, as n increases, patterning is moving towards the outer regions of the disc, with an increasing flat unpatterned space in the centre. (Compare the pattern for n = 10 with the shapes shown for nc = 6 and 8 in Fig. 1E and F.) A conspicuous feature of Fig. 10B–E is that, as the radius of the disc grows, the width of the patterned region stays almost constant. In particular, the distance from the first maximum to the second zero remains close to 6 on the scale of the radial variable x (distance between the two dotted lines). This distance corresponds to c/2. (Fig. 7 has diameters as ordinates, and intercept c; Fig. 10F, rotated 90°, has radii as ordinates, and intercept c/2.)

Thus the mathematics of harmonic waveforms on a disc accounts simultaneously for the inset circle and the spacing of cotyledons around it as illustrated in Fig. 6, provided that whatever wave-making mechanism is beginning to operate at developmental stage An(F) has a boundary condition U = 0 at the edge of the disc, i.e. that edge is a nodal circle; and the fit is quantitatively good (in terms of slope/intercept in Fig. 7) if the mechanism selects the second nodal circle. For any postulated mechanism, mechanical or physicochemical, it would be a substantial theoretical project to determine which nodal circle would be selected.


   LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 APPENDIX
 LITERATURE CITED
 

    Aida M, Ishida T, Fukaki H, Fujisawa H, Tasaka M. 1997. Genes involved in organ separation in Arabidopsis: an analysis of cup-shaped cotyledon mutant. Plant Cell 9: 841–857.[Abstract/Free Full Text]

    Aida M, Vernoux T, Furutani M, Traas J, Tasaka M. 2002. Roles of PIN-FORMED1 and MONOPTEROS in pattern formation of the apical region of the Arabidopsis embryo. Development 129: 3965–3974.[Abstract/Free Full Text]

    Barton MK, Poethig RS. 1993. Formation of the shoot apical meristem in Arabidopsis thaliana: an analysis in the wild type and in the shoot meristemless mutant. Development 119: 823–831.[Abstract/Free Full Text]

    Becwar MR, Nagmani R, Wann SR. 1990. Initiation of embryogenic cultures and somatic embryo development in loblolly pine (Pinus taeda). Canadian Journal of Forest Research 20: 810–817.[CrossRef]

    Bowman JL, Eshed Y. 2000. Formation and maintenance of the shoot apical meristem. Trends in Plant Science 5: 110–115.[CrossRef][Web of Science][Medline]

    Briand CH, Gauthier S, Lemay JF, Lalibert S, Tremblay M-F. 1998. Relationship between field performance, family, embryo morphology, and isozyme heterozygosity, and in vitro reactivity in jack pine. Canadian Journal of Forest Research 28: 98–105.[CrossRef]

    Buchholz JT. 1919. Studies concerning the evolutionary status of polycotyledony. American Journal of Botany 6: 106–119.[CrossRef][Web of Science]

    Butts D, Buchholz JT. 1940. Cotyledon numbers in conifers. Transactions of the Illinois Academy of Science 33: 58–62.

    Chaudhury AM, Letham S, Craig S, Dennis ES. 1993. amp1 – a mutant with high cytokinin levels and altered embryonic pattern, faster vegetative growth, constitutive photomorphogenesis and precocious flowering. Plant Journal 4: 907–916.[CrossRef][Web of Science]

    Choi YE, Kim HS, Soh WY, Yang DC. 1997. Developmental and structural aspects of somatic embryos formed on medium containing 2,3,5-triiodobenzoic acid. Plant Cell Reports 16: 738–744.[CrossRef][Web of Science]

    Conway LJ, Poethig RS. 1997. Mutations of Arabidopsis thaliana that transform leaves into cotyledons. Proceedings of the National Academy of Sciences, USA 94: 10209–10214.[Abstract/Free Full Text]

    Cooke TJ, Lu B. 1992. The independence of cell shape and overall form in multicellular algae and land plants: cells do not act as building blocks for constructing plant organs. International Journal of Plant Sciences 153: S7–S27.[CrossRef][Web of Science]

    Dumais J, Harrison LG. 2000. Whorl morphogenesis in the dasycladalean algae: the pattern formation viewpoint. Philosophical Transactions of the Royal Society, London B355: 281–305.[Abstract/Free Full Text]

    Green PB. 1999. Expression of pattern in plants: combining molecular and calculus-based biophysical paradigms. American Journal of Botany 86: 1059–1076.[Abstract/Free Full Text]

    Gutmann M, von Aderkas P, Label P, Lelu M-A. 1996. Effects of abscisic acid on somatic embryo maturation of hybrid larch. Journal of Experimental Botany 47: 1905–1917.[Abstract/Free Full Text]

    Harrison LG. 1993. Kinetic theory of living pattern. Cambridge and New York: Cambridge University Press.

    Harrison LG, Hillier NA. 1985. Quantitative control of Acetabularia morphogenesis by extracellular calcium: a test of kinetic theory. Journal of Theoretical Biology 114: 177–192.[CrossRef]

    Harrison LG, Donaldson G, Lau W, Lee M, Lin BP, Lohachitranont S, Setyawati I, Yue J. 1997. CaEGTA uncompetitively inhibits calcium activation of whorl morphogenesis in Acetabularia. Protoplasma 196: 190–196.[CrossRef][Web of Science]

    Harrison LG, Graham KT, Lakowski BC. 1988. Calcium localization during Acetabularia whorl formation: evidence supporting a two-stage hierarchical mechanism. Development 104: 255–262.[Abstract]

    Harrison LG, Snell J, Verdi R, Vogt DE, Zeiss GD, Green BR. 1981. Hair morphogenesis in Acetabularia mediterranea: temperature-dependent spacing and models of morphogen waves. Protoplasma 106: 211–221.[CrossRef][Web of Science]

    Harrison LG, Wehner S, Holloway DM. 2001. Complex morphogenesis of surfaces: theory and experiment on coupling of reaction-diffusion patterning to growth. Faraday Discussions, Royal Society of Chemistry, London 120: 277–294.

    Holloway DM, Harrison LG. 1999. Algal morphogenesis: modelling interspecific variation in Micrasterias with reaction-diffusion patterned catalysis of cell surface growth. Philosophical Transactions of the Royal Society London B354: 417–433.[CrossRef]

    Ingram GC, Goodrich J, Wikinson MD, Simon R, Haughn GW, Coen ES. 1995. Parallels between UNUSUAL FLORAL ORGANS and FIMBRIATA, genes controlling flower development in Arabidopsis and Antirrhinum. Plant Cell 7: 1501–1510.[Abstract]

    Juguet M. 1992. A scanning electron microscope study of cotyledon development in Pinus caribaea Morelet, with some remarks on the relationships between polycotyly and dicotyly. Bulletin de la Société Botanique de France 139: 235–245.

    Juguet M. 2002. Origine et étapes de la monocotylie chez les Monocotylédones. IV – La zone cotylédogène annulaire des Spermaphytes et l’origine de la monocotylie. Acta Botanica. Gallica 149: 3–33.[Web of Science]

    Kaplan DR, Cooke TJ. 1997. Fundamental concepts in the embryogenesis of dicotyledons: a morphological interpretation of embryo mutants. Plant Cell 9: 1903–1919.[CrossRef][Web of Science][Medline]

    Kaplan DR, Hagemann W. 1991. The relationship of cell and organism in vascular plants. BioScience 41: 693–703.[CrossRef][Web of Science]

    Kauzmann W. 1957. Quantum chemistry. New York: Academic Press, Ch. 3.

    Kuijt J. 1967. On the structure and origin of the seedling of Psittacanthus schiedeanus (Loranthaceae). Canadian Journal of Botany 45: 1497–1506.[CrossRef]

    Lacalli TC. 1981. Dissipative structures and morphogenetic pattern in unicellular algae. Philosophical Transactions of the Royal Society, London B294: 547–588.[CrossRef]

    Lehman A, Black R, Ecker JR. 1996. HOOKLESS1, an ethylene response gene, is required for differential cell elongation in the Arabidopsis hypocotyl. Cell 85: 183–194.[CrossRef][Web of Science][Medline]

    Lelu MA, Bastien C, Klimaszewska K, Ward C, Charest PJ. 1994. An improved method for somatic plantlet production in hybrid larch (Larix x leptoeuropaea). Part I. Somatic embryo maturation. Plant Cell, Tissue and Organ Culture 36, 107–115.

    Nagata W, Harrison LG, Wehner S. 2003. Reaction-diffusion models of growing plant tips: bifurcations on hemispheres. Bulletin of Mathematical Biology 65: 571–607. [CrossRef][Web of Science][Medline]

    Oladokun M. 1982. Fruit characteristics in Cola acuminata. 2. Nut size, cotyledon number and color. Turrialba 32: 275–281.

    Prigogine I, Lefever R. 1968. Symmetry-breaking instabilities in chemical systems. II. Journal of Chemical Physics 48: 1695–1700.[CrossRef]

    Simon R, Carpenter R, Doyle S, Coen E. 1994. FIMBRIATA controls flower development by mediating between meristem and organ identity genes. Cell 78: 99–107. [CrossRef][Web of Science][Medline]

    Sorensen FC, Franklin JF. 1977. Influence of year of cone collection on seed weight and cotyledon number in Abies procera. Silvae Genetica 26: 41–43.[Web of Science]

    Spurr AR. 1949. Histogenesis and organizatin of the embryo in Pinus atrobus L. American Journal of Botany 36: 629–641.[CrossRef][Web of Science]

    Torres-Ruiz RA, Jürgens G.1994. Mutations in the FASS gene uncouple pattern formation and morphogenesis in Arabidopsis development. Development 120: 2967–2978.

    Turing AM. 1952. The chemical basis of morphogenesis. Philosophical Transactions of the Royal Society, London B237: 37–72.[CrossRef]

    von Aderkas P. 2002. In vitro phenotypic variation in larch cotyledon number. International Journal of Plant Scences 163: 301–307.[CrossRef]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
J Exp BotHome page
J. W. Chandler
Cotyledon organogenesis
J. Exp. Bot., August 1, 2008; 59(11): 2917 - 2931.
[Abstract] [Full Text] [PDF]


Home page
ANN BOT (LOND)Home page
D. M. Holloway and L. G. Harrison
Pattern Selection in Plants: Coupling Chemical Dynamics to Surface Growth in Three Dimensions
Ann. Bot., February 1, 2008; 101(3): 361 - 374.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (1)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by HARRISON, L. G.
Right arrow Articles by VON ADERKAS, P.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by HARRISON, L. G.
Right arrow Articles by VON ADERKAS, P.
Agricola
Right arrow Articles by HARRISON, L. G.
Right arrow Articles by VON ADERKAS, P.
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?