In this viewpoint paper, we ask the following questions. Are we modelling the correct processes that drive plant growth, and is growth driven mostly by sink or source activity? In current models, is the importance of soil resources (nutrients, water, temperature and their interaction with meristematic activity) considered adequately? Do classic models account for architectural adjustment as well as integrating the fundamental principles of development? Whilst answering these questions with the available data in the literature, we put forward the opinion that plant architecture and sink activity must be pushed to the centre of plant growth models. In natural conditions, sinks will more often drive growth than source activity, because sink activity is often controlled by finite soil resources or developmental constraints.

This viewpoint paper also serves as an introduction to this Special Issue devoted to plant growth modelling, which includes new research covering areas stretching from cell growth to biomechanics. All papers were presented at the Second International Symposium on Plant Growth Modeling, Simulation, Visualization and Applications (PMA06), held in Beijing, China, from 13–17 November, 2006. Although a large number of papers are devoted to FS models of agricultural and forest crop species, physiological and genetic processes have recently been included and point the way to a new direction in plant modelling research.

Here we develop a concept for the role of models in forest ecosystem management based on historical analyses. Five paradigms of forest management are identified: (1) multiple uses, (2) dominant use, (3) environmentally sensitive multiple uses, (4) full ecosystem approach and (5) eco-regional perspective. An overview of model approaches is given that is dedicated to this purpose and to developments of different kinds of approaches. It is discussed how these models can contribute to goal setting, decision support and development of guidelines for forestry operations. Furthermore, it is shown how scenario analysis, including stand and landscape visualization, can be used to depict alternatives, make long-term consequences of different options transparent, and ease participation of different stakeholder groups and education.

In our opinion, the current challenge of forest ecosystem management in Europe is to integrate system knowledge from different temporal and spatial scales and from various disciplines. For this purpose, using a set of models with different focus that can be selected from a kind of toolbox according to particular needs is more promising than developing one overarching model, covering ecological, production and landscape issues equally well.

The modelling approach for plant internal resource allocation is based on a priority scheme assuming that growth processes have priority over allocation to secondary compounds and that growth-related metabolism is more strongly affected by nitrogen deficiency than defence-related secondary metabolism.

It is shown that the model can reproduce the effect of nitrogen fertilization on allocation patterns in apple trees and the effects of elevated CO₂ and competition in juvenile beech and spruce trees. The analysis of model behaviour reveals that large fluctuations in plant internal availability of carbon and nitrogen are possible within a single vegetation period. Furthermore, the model displays a non-linear allocation behaviour to carbon-based secondary compounds.

The simulation results corroborate the underlying assumptions of the presented modelling approach for resource partitioning between growth-related primary metabolism and defence-related secondary metabolism. Thus, the dynamical modelling approach, which considers variable source and sink strengths of plant internal resources within different phenological growth stages, presents a successful translation of existing concepts into a dynamic mathematical model.

An experiment was conducted on spring wheat (Triticum aestivum, ‘Minaret’), in a growth chamber from, 2004 to, 2005. Four harvests were made of six plants each to determine the size and mass of individual organs, including the root system, leaf blades, sheaths, internodes and ears of the main stem and different tillers. Leaf status (appearance, expansion, maturity and death) of these 24 plants was recorded. With the structures and mass of organs of four individual sample plants, the GREENLAB model was calibrated using a non-linear least-square-root fitting method, the aim of which was to minimize the difference in mass of the organs between measured data and model output, and to provide the parameter values of the model (the sink strengths of organs of each type, age and tiller order, and two empirical parameters linked to biomass production).

The masses of all measured organs from one plant from each harvest were fitted simultaneously. With estimated parameters for sink and source functions, the model predicted the mass and size of individual organs at each position of the wheat structure in a mechanistic way. In addition, there was close agreement between experimentally observed and simulated values of leaf area index.

In order to implement and test RGG, a new modelling language, the eXtended L-system language (XL) was created. Models using XL are interpreted by the interactive, Java-based modelling platform GroIMP. Three models, a semi-quantitative gibberellic acid (GA) signal transduction model, and a phytochrome-based shade detection and object avoidance model, both coupled to an existing morphogenetic structural model of barley (Hordeum vulgare L.), serve as examples to demonstrate the versatility and suitability of RGG and XL to represent the interaction of diverse biological processes across hierarchical scales.

The dynamics of the concentrations in the signal transduction network could be modelled qualitatively and the phenotypes of GA-response mutants faithfully reproduced. The light model used here was simple to use yet effective enough to carry out local measurement of red:far-red ratios. Suppression of tillering at low red:far-red ratios could be simulated.

The RGG formalism is suitable for implementation of multi-scaled FSPM of plants interacting with their environment via hormonal control. However, their ensuing complexity requires careful design. On the positive side, such an FSPM displays knowledge gaps better thereby guiding future experimental design.

The simulation of plant topology is based on the growth of a set of virtual buds whose activity is modelled using stochastic processes. The geometry of the resulting axes is modelled by simple descriptive functions. The potential growth of each bud is represented by means of a numerical value called physiological age, which controls the value for each parameter in the model. The set of possible values for physiological ages is called the reference axis. In order to mimic morphological and architectural metamorphosis, the value allocated for the physiological age of buds evolves along this reference axis according to an oriented finite state automaton whose occupation and transition law follows a semi-Markovian function.

Simulations were performed on tomato plants to demostrate how the AmapSim simulator can interface external modules, e.g. a GREENLAB growth model and a radiosity model.

The algorithmic ability provided by AmapSim, e.g. the reference axis, enables unified control to be exercised over plant development parameter values, depending on the biological process target: how to affect the local pertinent process, i.e. the pertinent parameter(s), while keeping the rest unchanged. This opening up to external functions also offers a broadened field of applications and thus allows feedback between plant growth and the physical environment.

The amount of photosynthetically active radiation absorbed by a plant is estimated on a daily or hourly basis through precise characterization of the light environment and three-dimensional virtual plants built using AMAP software. Several treatments are performed over four experiments and on two genotypes to test the model, quantify the contribution of different organs to light interception and evaluate the impact of heliotropism.

This approach is used to simulate the amount of light absorbed at organ and plant scales from crop emergence to maturity. Blades and capitula were the major contributors to light interception, whereas that by petioles and stem was negligible. Light regimen simulations showed that heliotropism decreased the cumulated light intercepted at the plant scale by close to 2·2 % over one day.

The approach is useful in characterizing the light environment of organs and the whole plant, especially for studies on heterogeneous canopies or for quantifying genotypic or environmental impacts on plant architecture, where conventional approaches are ineffective. This model paves the way to analyses of genotype–environment interactions and could help establish new selection criteria based on architectural improvement, enhancing plant light interception.

3-D virtual plants were reconstructed from field observations and 3-D digitization and were used to simulate the light regime in cotton stands of different densities.

All densities showed the same linear relationship between LAI and the sum of light intercepted by the canopy, from seedling emergence up to flowering. Simulated R : FR ratio profiles can very likely explain (1) the longer first internodes on main stem and branches and (2) the azimuthal re-orientation of branches toward the inter-row.

Simulation tools were used to analyse plant plasticity in terms of light quantity and quality. The methodology applied here at the stand scale will now be continued at the plant scale to further strengthen the above hypotheses.

This model is original in combining a turbid-medium-like envelope enclosing the volume occupied by vine shoots with the use of discrete geometric polygons representing leaves randomly located within this volume to represent plant structure. Reconstruction rules were adapted to capture the main determinants of grapevine shoot architecture and their variability. Using a simplified set of parameters, it was possible to describe (1) the 3D path of the main shoot, (2) the volume occupied by the foliage around this path and (3) the orientation of individual leaf surfaces. Model parameterization (estimation of the probability distribution for each parameter) was carried out for eight contrasting C x T pairs.

The parameter values obtained in each situation were consistent with our knowledge of grapevine architecture. Quantitative assessments for the generated virtual scenes were carried out at the canopy and plant scales. Light interception efficiency and local variations of light transmittance within and between experimental plots were correctly simulated for all canopies studied. The approach predicted these key ecophysiological variables significantly more accurately than the classical complete digitization method with a limited number of plants. In addition, this model accurately reproduced the characteristics of a wide range of individual digitized plants. Simulated leaf area density and the distribution of light interception among leaves were consistent with measurements. However, at the level of individual organs, the model tended to underestimate light interception.

Two field experiments were conducted on irrigated fields in the North China Plain with a block design of four replications. Each experiment included three PPDs: 2·8, 5·6 and 11·1 plants m^–2. Detailed observations were made on the dimensions and fresh biomass of above-ground plant organs for each phytomer throughout the seasons. Growth stage-specific target files (a description of plant organ weight and dimension according to plant topological structure) were established from the measured data required for GREENLAB parameterization. Parameter optimization was conducted using a generalized least square method for the entire growth cycles for all PPDs and years. Data from in situ plant digitization were used to establish geometrical symbol files for organs that were then applied to translate model output directly into 3-D representation for each time step of the model execution.

The analysis indicated that the parameter values of organ sink variation function, and the values of most of the relative sink strength parameters varied little among years and PPDs, but the biomass production parameter, computed plant projection surface and internode relative sink strength varied with PPD. Simulations of maize plant growth based on the fitted parameters were reasonably good as indicated by the linearity and slopes similar to unity for the comparison of simulated and observed values. Based on the parameter values fitted from different PPDs, shoot (including vegetative and reproductive parts of the plant) and cob fresh biomass for other PPDs were simulated. Three-dimensional representation of individual plant and plant stand from the model output with two contrasting PPDs were presented with which the PPD effect on plant growth can be easily recognized.

This study showed that GREENLAB model has the ability to capture plant plasticity induced by PPD. The relatively stable parameter values strengthened the hypothesis that one set of equations can govern dynamic organ growth. With further validation, this model can be used for agronomic applications such as yield optimization.

A greenhouse experiment was carried out with three homogeneous planting densities (plant spacing = 0·3, 0·6 and 1 m). Detailed records of plant development, plant architecture and organ growth were made throughout the growing period. Model calibration was performed for each situation using a statistical optimization procedure (multi-fitting).

Obvious limitations of the present version of the model appeared to account fully for the plant plasticity induced by inter-plant competition for light. A lack of stability was identified for some model parameters at very high planting density. In particular, those parameters characterizing organ sink strengths and governing light interception proved to be environment-dependent. Remarkably, however, responses of the parameter values concerned were consistent with actual growth measurements and with previously reported results. Furthermore, modifications of total biomass production and of allocation patterns induced by the planting-density treatments were accurately simulated using the sets of optimized parameters. These results demonstrate that the overall model structure is potentially able to reproduce the observed plant plasticity and suggest that sound biologically based adaptations could overcome the present model limitations. Potential options for model improvement are proposed, and the possibility of using the kernel algorithm currently available as a fitting tool to build up more sophisticated model versions is advocated.

The empirical production equation of GREENLAB is extrapolated to stands by computing the exposed photosynthetic foliage area of each plant. The computation is based on the combination of Poisson models of leaf distribution for all the neighbouring plants whose crown projection surfaces overlap. To study the effects of density on architectural development, we link the proposed competition model to the model of interaction between functional growth and structural development introduced by Mathieu (2006, PhD Thesis, Ecole Centrale de Paris, France).

The model is applied to mono-specific field crops and forest stands. For high-density crops at full cover, the model is shown to be equivalent to the classical equation of field crop production ( Howell and Musick, 1985, in Les besoins en eau des cultures; Paris: INRA Editions). However, our method is more accurate at the early stages of growth (before cover) or in the case of intermediate densities. It may potentially account for local effects, such as uneven spacing, variation in the time of plant emergence or variation in seed biomass. The application of the model to trees illustrates the expression of plant plasticity in response to competition for light. Density strongly impacts on tree architectural development through interactions with the source–sink balances during growth. The effects of density on tree height and radial growth that are commonly observed in real stands appear as emerging properties of the model.

A three-dimensional forest simulator is used to explore the impact of crown shape plasticity on tree growth. Crown deformation is mediated through the local response to light and overall allometric constraints governing tree dimensions. By altering shape response parameters of Hevea brasiliensis the impact of increased or decreased plasticity is explored in a variety of competitive environments defined by various combinations of tree density and relative frequency of different strategies. The possible interactions between plasticity and growth rate and plasticity and below-ground competition are also explored.

Crown plasticity confers competitive superiority in all cases studied. Interactions with other processes may downplay or enhance this competitive advantage.

Simulation results strongly suggest that crown plasticity does have a significant impact on tree performance in nature and that commonly observed crown shape deformation response of trees is of adaptive value.

In the model, a single variable, the source–sink ratio controls different events in plant architecture. In particular, the number of fruits and branch formation are determined as increasing functions of this ratio. For some sets of well-chosen parameters of the model, the dynamical evolution of the ratio during plant growth generates rhythms.

Cyclic patterns in branch formation or fruit appearance emerge without being forced by the model. The model is based on the theory of discrete dynamical systems. The mathematical formalism helps us to explain rhythm generation and to control the behaviour of the system. Rhythms can appear during both the exponential and stabilized phases of growth, but the causes are different as shown by an analytical study of the system. Simulated plant behaviours are very close to those observed on real plants. With a small number of parameters, the model gives very interesting results from a qualitative point of view. It will soon be subjected to experimental data to estimate the model parameters.

The GREENLAB model was selected as a reasonable choice to link growth model parameters to QTL. Virtual genes and virtual chromosomes were defined to build a simple genetic model that drove the settings of the species-specific parameters of the model. The QTL Cartographer software was used to study QTL detection of simulated plant traits. A genetic algorithm was implemented to define the ideotype for yield maximization based on the model parameters and the associated allelic combination.

By keeping the environmental factors constant and using a virtual population with a large number of individuals generated by a Mendelian genetic model, results for an ideal case could be simulated. Virtual QTL detection was compared in the case of phenotypic traits – such as cob weight – and when traits were model parameters, and was found to be more accurate in the latter case. The practical interest of this approach is illustrated by calculating the parameters (and the corresponding genotype) associated with yield optimization of a GREENLAB maize model. The paper discusses the potentials of GREENLAB to represent environment x genotype interactions, in particular through its main state variable, the ratio of biomass supply over demand.

A new generic model of plant cellular morphogenesis is described that expresses interactions amongst cellular entities explicitly: the plant is described as a multi-scale structure, and interactions between distinct entities is established through a topological neighbourhood. Tissues are represented as 2D biphasic systems where the cell wall responds to turgor pressure through a viscous yielding of the cell wall.

This principle was used in the development of the CellModeller software, a generic tool dedicated to the analysis and modelling of plant morphogenesis. The system was applied to three contrasting study cases illustrating genetic, hormonal and mechanical factors involved in plant morphogenesis.

Plant morphogenesis is fundamentally a cellular process and the CellModeller software, through its underlying generic model, provides an advanced research tool to analyse coupled physical and biological morphogenetic mechanisms.

The FEM was used to carry out 2-D simulations of tree uprooting in saturated soft clay and loamy sand-like soil. The anchorage model consisted of a root system embedded in a soil block. Two root patterns were used and individual roots removed to determine their contribution to anchorage.

In clay-like soil the size of the root–soil plate formed during overturning was defined by the longest roots. Consequently, all other roots localized within this plate had no influence on anchorage strength. In sand-like soil, removing individual root elements altered anchorage resistance. This result was due to a modification of the shape and size of the root–soil plate, as well as the location of the rotation axis. The tap root and deeper roots had more influence on overturning resistance in sand-like soil compared with clay-like soil. Mechanical stresses were higher in the most superficial roots and also in leeward roots in sand-like soil. The relative difference in stresses between the upper and lower sides of lateral roots was sensitive to root insertion angle. Assuming that root eccentricity is a response to mechanical stresses, these results explain why eccentricity differs depending on root architecture.

A simple 2-D Finite Element model was developed to better understand the mechanisms involved during tree overturning. It has been shown how root system morphology and soil mechanical properties can modify the shape of the root plate slip surface as well as the position of the rotation axis, which are major components of tree anchorage.

Three-dimensional digitizing in situ was used to obtain accurate root system architecture data for mature Quercus alba in two forest stands. These data were used as input to tools developed, which analyse the spatial position of roots, topology and geometry. The contribution of roots to soil reinforcement was determined by calculating additional soil cohesion using the limit equilibrium model, and the factor of safety (FOS) using an existing slope stability model, Slip4Ex.

Existing models may incorrectly estimate the additional soil cohesion provided by roots, as the spatial position of roots crossing the potential slip surface is usually not taken into account. However, most soil reinforcement by roots occurs close to the tree stem and is negligible at a distance >1·0 m from the tree, and therefore global values of FOS for a slope do not take into account local slippage along the slope.

Within a forest stand on a landslide-prone slope, soil fixation by roots can be minimal between uniform rows of trees, leading to local soil slippage. Therefore, staggered rows of trees would improve overall slope stability, as trees would arrest the downward movement of soil. The chain of tools consisting of both software (free for non-commercial use) and functions available from the first author will enable a more accurate description and use of root architectural parameters in standard slope stability analyses.

After a general introduction to the problems of water availability, this review focuses on a critical evaluation of recent progress in unravelling mechanisms for modifying plant growth responses to drought.

Investigations of key regulatory mechanisms integrating plant growth responses to water deficits at the whole-organism, cellular and genomic levels continue to provide novel and exiting research findings. For example, recent reports contradict the widespread conception that root-derived abscisic acid is necessarily involved in signalling for stomatal and shoot-growth responses to soil water deficits. The findings bring into question the theoretical basis for alternate-side root-irrigation techniques. Similarly, recent reports indicate that increased ABA production or increased aquaporin expression did not lead to improved drought resistance. Other reports have concerned key genes and proteins involved in regulation of flowering (FT), vegetative growth (DELLA), leaf senescence (IPT) and desiccation tolerance (LEA). Introgression of such genes, with suitable promoters, can greatly impact on whole-plant responses to drought. Further developments could facilitate the introduction by breeders of new crop varieties with growth physiologies tailored to improved field performance under drought. Parallel efforts to encourage the introduction of supplementary irrigation with water made available by improved conservation measures and by sea- or brackish-water desalination, will probably provide comprehensive solutions to coping with drought-prone environments.

Eleven Hypochaeris species from 18 populations were studied using fluorescence in situ hybridization (FISH) with 35S and 5S rDNA probes. A phylogenetic framework was established from neighbour-net analysis of amplified fragment length polymorphism (AFLP) fingerprint data.

A single 5S rDNA locus is invariably found on the short arm of chromosome 2. Using 35S rDNA loci, based on number (one or two) and localization (interstitial on the long arm of chromosome 2, but sometimes lacking, and terminal or interstitial on the short arm of chromosome 3, only very rarely lacking), seven karyotype groups can be distinguished; five of these include polyploids. Karyotype groups with more than one species do not form monophyletic groups.

Early evolution of Hypochaeris in South America was characterized by considerable karyotype differentiation resulting from independent derivations from an ancestral karyotype. There was marked diversification with respect to the position and evolution of the 35S rDNA locus on chromosome 3, probably involving inversions and/or transpositions, and on chromosome 2 (rarely 3) concerning inactivation and loss. Among these different karyotype assemblages, the apargioides group and its derivatives constitute by far the majority of species.

Genomic in situ hybridization (GISH) and fluorescent in situ hybridization (FISH) with ribosomal DNA (rDNA) probes, together with sequencing of the internal transcribed spacer (ITS) region of the rDNA, were used to identify the origin of P. egaliksensis and to explore its genomic organization, particularly at rDNA loci.

GISH showed that P. egaliksensis inherited all chromosomes from P. mistassinica and P. nutans and did not reveal major intergenomic rearrangements between the parental genomes (e.g. interchromosomal translocations). However, karyological comparisons and FISH experiments suggested small-scale rearrangements, particularly at rDNA sites. Primula egaliksensis lacked the ITS-bearing heterochromatic knobs characteristic of the maternal parent P. mistassinica and maintained only the rDNA loci of P. nutans. These results corroborated sequence data indicating that most ITS sequences of P. egaliksensis were of the paternal repeat type.

The lack of major rearrangements may be a consequence of the considerable genetic divergence between the putative parents, while the rapid elimination of the ITS repeats from the maternal progenitor may be explained by the subterminal location of ITS loci or a potential role of nucleolar dominance in chromosome stabilization. These small-scale rearrangements may be indicative of genome diploidization, but further investigations are needed to confirm this assumption.

Auxin supply, light conditions and sucrose concentration were modified with the aim of establishing efficient root proliferation conditions. The morphological variability among 62 established hairy root clones was phenotyped by scanning the roots and analysing the images using ‘whinRHIZO’ software procedures. PCR analysis of integration in transformed root cells of rol and aux oncogenes from the T-DNA of the Ri plasmid was used to study the molecular variability among clones.

Auxin supply was necessary to obtain and stimulate growth and branching, and IBA applied at 0·5 µm was the most efficient auxin. Significant differences were shown among the 62 clones for total root length and for the percentage of fine roots. These variables were stable across subcultures and could hence be used for efficient characterization of hairy root clones. The majority of hairy root clones (86 %) exhibited non-significant phenotype differences with non-transformed roots. Eight clones were significantly different from the non-transformed controls in that they possessed a low proportion of fine roots. Two other hairy root clones grew significantly faster than the other clones. The PCR analysis revealed a low variability in the integration of rol and aux oncogenes in transformed root cells. The T_R-DNA was never integrated as aux1 and aux2 genes were not found, although rolB and rolC genes from the T_L-DNA were always present.

The discovery of low morphological variability among coffee hairy roots together with the identification of morphological variables allowing easy identification of phenotypically altered clones represent two important results. They make hairy roots a possible, and efficient, tool for functional-genomic studies of coffee root genes.

Using light microscopy, ovule and embryo sac development are described in five (of 12) species of Trithuria, the sole genus of Hydatellaceae, and compared with those of Cabombaceae and Nymphaeaceae.

The ovule of Trithuria is bitegmic and tenuinucellate, rather than bitegmic and crassinucellate as in most other Nymphaeales. The seed is operculate and possesses a perisperm that develops precociously, which are both key features of Nymphaeales. However, in the Indian species T. konkanensis, perisperm is relatively poorly developed by the time of fertilization. Perisperm cells in Trithuria become multinucleate during development, a feature observed also in other Nymphaeales. The outer integument is semi-annular (‘hood-shaped’), as in Cabombaceae and some Nymphaeaceae, in contrast to the annular (‘cap-shaped’) outer integument of some other Nymphaeaceae (e.g. Barclaya) and Amborella. The megagametophyte in Trithuria is monosporic and four-nucleate; at the two-nucleate stage both nuclei occur in the micropylar domain. Double megagametophytes were frequently observed, probably developed from different megaspores of the same tetrad. Indirect, but strong evidence is presented for apomictic embryo development in T. filamentosa.

Most features of the ovule and embryo sac of Trithuria are consistent with a close relationship with other Nymphaeales, especially Cabombaceae. The frequent occurrence of double megagametophytes in the same ovule indicates a high degree of developmental flexibility, and could provide a clue to the evolutionary origin of the Polygonum-type of angiosperm embryo sac.

-control treatments were monitored for 10 years using six replicated experiments and analysed using meta-analysis. Meta-regressions were used to explore heterogeneity between sites.

The effectiveness of treatments varied between sites depending on the measure used to assess P. aquilinum performance. In general, cutting twice per year was the most successful treatment but on some sites other, less expensive treatments were as good. The effectiveness of treatments at different sites was not related to rhizome mass, but the effectiveness of most applied treatments were inversely related to post-control litter. Effective treatment was also associated with high species richness.

It is concluded that successful development of national weed control programs requires multi-site experimental approaches. Here, meta-analyses demonstrate that variation in effectiveness between sites could be explained in part by pre-specified variables. Reliance on data from a single site for policy formulation is therefore clearly dangerous.

Corn (Zea mays ‘Hycorn 82’) and mungbean (Vigna radiata ‘Berken’) were grown in continuous flowing nutrient solutions containing 0, 0·2, 1·0 and 5·0 µm La or Ce. At harvest plants were divided into roots and shoots, dried, weighed and analysed for macro- and micronutrients, as well as for La and Ce.

La and Ce did not increase the growth of corn or mungbean. The dry weight of corn shoots was decreased by 32 % in the presence of 5·0 µm Ce; the other La and Ce concentrations had no effect. La and Ce concentrations of 0·9 and 5·0 µm decreased the shoot dry weight of mungbean by 75 or 95 %, the two elements having closely similar effects. Decreases in the uptake of Ca, Na, Zn and Mn by corn were observed with increases in solution La and Ce. For mungbean, the uptake rates of all measured elements decreased with increases in solution La and Ce. The concentrations of La and Ce in the roots of both species were higher than in the shoots and increased strongly with increasing concentrations of La or Ce in solution. The La and Ce concentrations in mungbean shoots were always higher than in corn shoots.

La and Ce did not enhance the growth of corn or mungbean, but decreased the growth, root function and consequently the nutritional status of mungbean at concentrations >0·2 µm in solution. It is concluded that if La or Ce have positive effects on corn and mungbean growth, they can only occur at solution concentrations below 0·2 µm.

Transmission electron microscopy (TEM), scanning electron microscopy (SEM) and both bright-field and polarized-light microscopy were used to study a range of developmental stages.

Raphide crystals are initiated very early in plant development. They are consistently present in most species and have a fairly uniform distribution within mature tissues. Individual raphides may be formed by calcium oxalate deposition within individual crystal chambers in the vacuole of an idioblast. Druse crystals form later in the true leaves, and are absent from some species. Distribution of druses within leaves is more variable. Druses initially develop at leaf tips and then increase basipetally as the leaf ages. Druse development may also be initiated in crystal chambers.

The unusual grooved raphides in Amorphophallus species probably result from an unusual crystal chamber morphology. There are multiple systems of transport and biomineralization of calcium into the vacuole of the idioblast. Differences between raphide and druse idioblasts indicate different levels of cellular regulation. The relatively early development of raphides provides a defensive function in soft, growing tissues, and restricts build-up of dangerously high levels of calcium in tissues that lack the ability to adequately regulate calcium. The later development of druses could be primarily for calcium sequestration.

Two F₂ populations were generated from crosses between ‘Angel’ and varieties of intolerant M. truncatula, the male-sterile mutant tap and the cultivar ‘Caliph’. Genetic mapping with SSR (single sequence repeat) and gene-based markers allowed identification of the trait-defining gene. Quantitative gene expression studies showed the activity of the respective alleles.

Segregation ratios indicated the control of SU-herbicide tolerance by a single dominant gene. SU herbicides inhibit the biosynthesis of the branched-chain amino acids by targeting the acetolactate synthase enzyme, allowing the choice of a mapping approach using acetolactate synthase (ALS) gene homologues as candidates. SSR-marker analysis suggested the ALS-gene homologue on chromosome 3 in M. truncatula. The ALS-gene sequences from ‘Angel’ and intolerant genotypes were sequenced. In ‘Angel’, a single point mutation from C to T translating into an amino acid change from proline to leucine was identified. The polymorphism was used to develop a diagnostic marker for the tolerance trait. Expression of the mutant ALS allele was confirmed by quantitative RT-PCR and showed no differences at various seedling stages and treatments to the corresponding wild-type allele.

The identification of the trait-defining gene and the development of a diagnostic marker enable efficient introgression of this economically important trait in annual medic improvement programs.

Four experiments were performed on two different species (Arabidopsis thaliana and Helianthus annuus) in which the area of leaves that had apparently reached their final size was measured upon reversal of water stresses of different intensities and durations.

Re-growth of leaves that had apparently reached their final size occurred in both species, and its magnitude depended only on the time elapsed from growth cessation to rewatering. Leaf area increased up to 186% in A. thaliana and up to 88% in H. annuus after rewatering, with respect to the leaves of plants that remained under water deficit. Re-growth was accounted for by cell expansion. Increase in leaf area represented actual growth and not only a reversible change due to increased turgor.

After the leaf has ceased to grow, leaf cells retain their ability to expand for several days before leaf size becomes fixed. A response window was identified in both species, during which the extent of leaf area recovery decreased with time after the ‘initial’ leaf growth cessation. These results suggest that re-growth after rewatering of leaves having apparently attained their final size could be a generalized phenomenon, at least in dicotyledonous plants.

Several components of male fertility were analysed. Pollination rates in the male phase were recorded across one flowering period. Pollen viability was compared among morphs through hand pollinations both with pollen from a single sexual morph and also simulating a situation of pollen competition; in the latter experiment, paternity was assessed with microsatellite markers. It was also determined whether effects of genetic relatedness between pollen donors and recipients could influence the siring success. Finally, paternal effects occurring beyond the fertilization process were tested for by measuring the height reached by seedlings with different sires over three consecutive growing seasons.

The males and protandrous morphs had higher pollination rates than the protogynous morph, and the seedlings they sired grew taller. No differences in male fertility were found between males and protandrous individuals. Departures from random mating due to effects of genetic relatedness among sires and pollen recipients were also ruled out.

Males and protandrous individuals are probably better sires than protogynous individuals, as shown by the higher pollination rates and the differential growth of the seedlings sired by these morphs. In contrast, the fertility of males was not higher than the male fertility of the protandrous morph. While the appearance of males in sexually specialized heterodichogamous populations is possible, even in the absence of a fitness advantage, it is not clear that males can be maintained at an evolutionary equilibrium with two classes of heterodichogamous hermaphrodites.

Inflorescences of P. billietiae at different stages of development were observed using SEM. The number of appendices in male, female and sterile flowers were counted on 11 young inflorescences (5–6 flowers per inflorescence). The number of staminodes and carpels on ABFs were counted on 19 inflorescences (n = 143). These data were used for regression and ANOVA analyses.

There was an average of 4·1 stamens per male flower, 9·8 carpels per female flower and 6·8 staminodes per sterile male flower. There was an average of 7·3 floral appendices per atypical flower. Staminodes and carpels are inserted on the same whorl in ABFs. A negative exponential relationship was found between the average number of staminodes and the number of carpels in ABFs. If only the ABFs consisting of less than six carpels are considered, there is a linear relationship between the number of carpels and the average number of staminodes. The value of the slope of the regression equation indicates that on average, in P. billietiae, 1·36 carpels are replaced by one staminode.

In P. billietiae, the number of appendages in female flowers imposes a constraint on the maximum total number of appendages (carpels and staminodes) that can develop on ABFs. The quantitative analyses indicate that the average number of different types of floral appendages on an ABF and the number of organs involved in a homeotic transformation are two independent phenomena.

The morphology and distribution of the structures were studied under light microscopy, and their chemical composition was analysed using thin-layer chromatography and high-performance liquid chromatography. To further investigate the function of these structures, several simple field experiments and observations were also conducted.

Three types of small structures on P. tomentosa were investigated: bowl-shaped organs, glandular hairs and dendritic trichomes. The bowl-shaped organs were densely aggregated on the leaves near flower buds and were determined to be extrafloral nectarines (EFNs) that secrete sugar and attract ants. Nectar production of these organs was increased by artificial damage to the leaves, suggesting an anti-herbivore function through symbiosis with ants. Glandular hairs were found on the surfaces of young and/or reproductive organs. Glandular hairs on leaves, stems and flowers secreted mucilage containing glycerides and trapped small insects. Secretions from glandular hairs on flowers and immature fruits contained flavonoids, which may provide protection against some herbivores. Yellow dendritic trichomes on the adaxial side of leaves also contained flavonoids identical to those secreted by the glandular hairs on fruits and flowers. Three special types of leaves, which differed from the standard leaves in shape, size and identity of small structures, developed near young shoot tips or young flower buds. The density of small structures on these leaf types was higher than on standard leaves, suggesting that these leaf types may be specialized to protect young leaves or reproductive organs. Changes in the small structures during leaf development suggested that leaves of P. tomentosa are primarily protected by glandular hairs and dendritic trichomes at young stages and by the EFNs at mature stages.

The results indicate that P. tomentosa protects young and/or reproductive organs from herbivores through the distribution and allocation of small structures, the nature of which depends on the developmental stage of leaves and shoots.

Data for 110 published cell cycle times of seedlings grown at temperatures between 20–25 °C were compared with DNA C-values (58 values for monocots and 52 for eudicots). Regression analyses were undertaken for all species, and separately for monocots and eudicots, diploids and polyploids, and annuals and perennials. Cell cycle times were plotted against the nuclear DNA C-values.

A positive relationship was observed between DNA C-value and cell cycle time for all species and for eudicots and monocots separately, regardless of the presence or absence of polyploid values. In this sample, among 52 eudicots the maximum cell cycle length was 18 h, whereas the 58 monocot values ranged from 8–120 h. There was a striking additional increase in cell cycle duration in perennial monocots with C-values greater than 25 pg. Indeed, the most powerful relationship between DNA C-value and cell cycle time and the widest range of cell cycle times was in perennials regardless of ploidy level.

DNA replication is identified as a rate limiting step in the cell cycle, the flexibility of DNA replication is explored, and we speculate on how the licensing of initiation points of DNA replication may be a responsive component of the positive nucleotypic effect of C-value on the duration of the mitotic cell cycle.

The results show that genome size is a strong predictor of phenotypic traits at the cellular level (guard cell length and epidermal cell area had significant positive relationships with genome size). Stomatal density decreased with increasing genome size, but this did not lead to decreased photosynthetic rate. At higher phenotypic scales, the predictive power of genome size generally diminishes (genome size had weak predictive power for both LMA and seed mass), except in the interesting case of maximum plant height (tree species tend to have small genomes). There was no relationship with wood density. The general observation that species with larger genome size have larger seed mass was supported; however, species with small genome size can also have large seed masses. All of these analyses involved robust comparative methods that incorporate the phylogenetic relationships of species.

Genome size correlations are quite strong at the cellular level but decrease in predictive power with increasing phenotypic scale. Our hope is that these results may lead to new mechanistic hypotheses about why genome size scaling exists at the cellular level, and why nucleotypic consequences diminish at higher phenotypic scales.

It is a sine qua non that Bs originate from the A chromosomes, in a variety of ways. We can study their modes of drive and ask how it is that chromosomes which apparently lack genes can have control over their own drive process which leads to their survival in natural populations. Molecular cytogenetic studies are opening up new avenues of investigation. Population equilibria for B frequencies are determined by a balance between accumulation and harmful effects. Bs are also subject to meiotic loss due to polysomy and to elimination at meiosis as univalents. These balancing forces can be seen in the context of host/parasite interaction, based on a dissection of the genetic elements in both As and Bs (in maize) which interact to bring about a stable equilibrium, at least for a snapshot in time.

Aside from their intrinsic enigmatic properties, B chromosomes make useful experimental tools to study genome organization. Thus far they have not been exploited for their applications, other than through the use of A-B translocations used for gene mapping in maize; but there are opportunities to use them to modulate the frequency and distribution of recombination, to diploidize allopolyploids, to study centromeres and to be developed as plant artificial chromosomes; given that they can be structurally modified and their inheritance stabilized.

DNA content was estimated as fluorescence of propidium iodide (PI)-stained nuclei of poinsettia and/or pea (Pisum sativum) using flow cytometry. Tissue was chopped, or two tissues co-chopped, in Galbraith buffer alone or with six concentrations of cyanidin-3-rutinoside (a cyanidin-3-rhamnoglucoside contributing to red coloration in poinsettia).

There were large differences in PI staining (35–70 %) between 2C nuclei from green leaf and red bract tissue in poinsettia. These largely disappeared when pea leaflets were co-chopped with poinsettia tissue as an internal standard. However, smaller (2·8–6·9 %) differences remained, and red bracts gave significantly lower 1C genome size estimates (1·69–1·76 pg) than green leaves (1·81 pg). Chopping pea or poinsettia tissue in buffer with 0–200 µm cyanidin-3-rutinoside showed that the effects of natural inhibitors in red bracts of poinsettia on PI staining were largely reproduced in a dose-dependent way by this anthocyanin.

Given their near-ubiquitous distribution, many suspected roles and known affects on DNA staining, anthocyanins are a potent, potential cause of significant error variation in genome size estimations for many plant tissues and taxa. This has important implications of wide practical and theoretical significance. When choosing genome size calibration standards it seems prudent to select materials producing little or no anthocyanin. Reviewing the literature identifies clear examples in which claims of intraspecific variation in genome size are probably artefacts caused by natural variation in anthocyanin levels or correlated with environmental factors known to induce variation in pigmentation.

The purpose of this article is to present a selective overview of the field of nuclear DNA content measurement, and C-values in particular, with a focus on the technical difficulties imposed by the characteristics of the biological material and with some comments on the photometrical aspects of the work. For over 20 years it has been known that plant polyphenols cause problems in Feulgen DNA cytophotometry, since they act as major staining inhibitors leading to unreliable results. However, little information is available about the chemical classes of plant metabolites capable of DNA staining interference and the mechanisms of their inhibition. Plant slimes are another source of concern.

In FCM research to uncover the effects of secondary metabolites on measurement results has begun only recently. In particular, the analysis of intraspecific genome size variation demands a stringent methodology which accounts for inhibitors. FCM tests for inhibitory effects of endogenous metabolites should become obligatory. The use of dry seeds for harvesting embryo and endosperm nuclei for FCM and Feulgen densitometry may often provide a means of circumventing staining inhibitors. The importance of internal standardization is highlighted. Our goal is a better understanding of phytochemical/cytochemical interactions in plant DNA photometry for the benefit of an ever-growing list of plant genome sizes.

Genome sizes were measured in 18 species of Nicotiana (nine diploids and nine polyploids) ranging in age from <200 000 years to approx. 4·5 Myr old, to determine the direction and extent of genome size change following polyploidy. These data were combined with data from genomic in situ hybridization and increasing amounts of information on sequence composition in Nicotiana to provide insights into the molecular basis of genome size changes.

By comparing the expected genome size of the polyploid (based on summing the genome size of species identified as either a parent or most closely related to the diploid progenitors) with the observed genome size, four polyploids showed genome downsizing and five showed increases. There was no discernable pattern in the direction of genome size change with age of polyploids, although with increasing age the amount of genome size change increased. In older polyploids (approx. 4·5 million years old) the increase in genome size was associated with loss of detectable genomic in situ hybridization signal, whereas some hybridization signal was still detected in species exhibiting genome downsizing. The possible significance of these results is discussed.

Allopolyploid species of Nicotiana are ideal for studying rDNA evolution because phylogenetic reconstruction of DNA sequences has revealed patterns of species divergence and their parents. From these studies we also know that polyploids formed over widely different timeframes (thousands to millions of years), enabling comparative and temporal studies of rDNA structure, activity and chromosomal distribution. In addition studies on synthetic polyploids enable the consequences of de novo polyploidy on rDNA activity to be determined.

We propose that rDNA epigenetic expression patterns established even in F₁ hybrids have a material influence on the likely patterns of divergence of rDNA. It is the active rDNA units that are vulnerable to homogenization, which probably acts to reduce mutational load across the active array. Those rDNA units that are epigenetically silenced may be less vulnerable to sequence homogenization. Selection cannot act on these silenced genes, and they are likely to accumulate mutations and eventually be eliminated from the genome. It is likely that whole silenced arrays will be deleted in polyploids of 1 million years of age and older.

This review summarizes recent genomics studies implemented in hexaploid and octoploid triticales and discusses the mechanisms of the changes and compares the major differences between genomic changes in triticale and other allopolyploid species.

Molecular studies have indicated extensive non-additive sequence changes or modifications in triticale, and the degree of variation appears to be higher than in other allopolyploid species. The data indicate that at least some sequence changes are non-random, and appear to be a function of genome relations, ploidy levels and sequence types. Specifically, the rye parental genome demonstrated a higher level of changes than the wheat genome. The frequency of lost parental bands was much higher than the frequency of gained novel bands, suggesting that sequence modification and/or elimination might be a major force causing genome variation in triticale. It was also shown that 68 % of the total changes occurred immediately following wide hybridization, but before chromosome doubling. Genome evolution following chromosome doubling occurred more slowly at a very low rate and the changes were mainly observed in the first five or so generations. The data suggest that cytoplasm and relationships between parental genomes are key factors in determining the direction, amount, timing and rate of genomic sequence variation that occurred during inter-generic allopolyploidization in this system.