# Cauliflowers or how the perseverance of a plant to make flowers produces an amazing fractal structure

Eugenio Azpeitia, François Parcy, Christophe Godin

► To cite this version:

Eugenio Azpeitia, François Parcy, Christophe Godin. Cauliflowers or how the perseverance of a plant to make flowers produces an amazing fractal structure. Comptes Rendus. Biologies, 2023, 346 (G1), pp.75-83. 10.5802/crbiol.120. hal-04242307

HAL Id: hal-04242307

https://inria.hal.science/hal-04242307v1

Submitted on 14 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
# Comptes Rendus
Biologies

Eugenio Azpeitia, François Parcy and Christophe Godin

Cauliflowers or how the perseverance of a plant to make flowers produces an
amazing fractal structure

Volume 346 (2023), p. 75-83

Published online: 14 June 2023

https://doi.org/10.5802/crbiol.120

<figure><img src="image_2.png" /></figure>

This article is licensed under the
CREATIVE COMMONS ATTRIBUTION 4.0 INTERNATIONAL LICENSE.
http://creativecommons.org/licenses/by/4.0/

<figure><img src="image_3.png" /></figure>

Les Comptes Rendus. Biologies sont membres du
Centre Mersenne pour l'édition scientifique ouverte
www.centre-mersenne.org
e-ISSN : 1768-3238
Review Article / *Article de revue*

# Cauliflowers or how the perseverance of a plant to make flowers produces an amazing fractal structure

*Le chou-fleur ou comment la persévérance d'une plante à vouloir faire des fleurs produit une étonnante structure fractale*

Eugenio Azpeitia<sup>①, a</sup>, François Parcy<sup>①, b</sup> and Christophe Godin<sup>①, c</sup>

<sup>a</sup> Centro de Ciencias Matemáticas, Universidad Nacional Autónoma de México, Morelia, México

<sup>b</sup> Laboratoire Physiologie Cellulaire et Végétale, Univ. Grenoble Alpes, CNRS, CEA, INRAE, IRIG-DBSCI-LPCV, F-38054 Grenoble, France

<sup>c</sup> Laboratoire de Reproduction et Développement des Plantes, Univ. Lyon, ENS de Lyon, UCB Lyon 1, CNRS, INRAE, Inria, F-69364 Lyon, France

*E-mails*: eazpeitia@matmor.unam.mx (E. Azpeitia), francois.parcy@cea.fr (F. Parcy), christophe.godin@inria.fr (C. Godin)

**Abstract.** Biological organisms have an immense diversity of forms. Some of them exhibit conspicuous and fascinating fractal structures that present self-similar patterns at all scales. How such structures are produced by biological processes is intriguing. In a recent publication, we used a multi-scale modelling approach to understand how gene activity can produce macroscopic cauliflower curds. Our work provides a plausible explanation for the appearance of fractal-like structures in plants, linking gene activity with development.

**Résumé.** Les organismes biologiques possèdent une immense diversité de formes. Certains d'entre eux arborent des structures fractales remarquables et fascinantes présentant des motifs similaires à toutes les échelles. La manière dont ces structures sont produites par des processus biologiques est intrigante. Dans une publication récente, nous avons utilisé une approche de modélisation multi-échelle pour comprendre comment l'activité de certains gènes peut engendrer des pommes de chou-fleur macroscopiques. Notre travail fournit une explication plausible quant à l'apparition de structures de type fractal chez les plantes, en reliant l'activité des gènes au développement.

**Keywords.** Biological fractals, Plant development, Multiscale modelling, Cauliflower curds, Flower, Meristems, Plant shapes.

**Mots-clés.** Fractales biologiques, Développement des plantes, Modélisation multi-échelle, Pommes de chou-fleur, Fleurs, Méristèmes, Formes végétales.

*Published online: 14 June 2023*

<footnote>* Corresponding authors.</footnote>
# 1. Introduction

Fractals are intriguing by their mathematical and geometric beauty [1]. Their intricate shapes are made up by the repetition of a similar pattern at every scale. Despite their complexity, they can be easily simulated using intensive recursive computation on standard computers (see e.g. [2]). Biology provides a wealth of illustrations of fractal-like shapes. However, none of the computational algorithms used to generate fractals consider the complex chemical and physical processes responsible for the development of biological shapes, therefore giving little insight into their origins in nature, which thus remains largely unknown.

Cauliflower curds provide a unique opportunity to address this question. Their fractal-like shape derives from a modification of the flower development process. Strikingly, a cauliflower-like severe change of morphology is produced by mutations of the *APETALA1* (*API*) and *CAULIFLOWER* (*CAL*) genes in *Arabidopsis thaliana* [3, 4], a brassicaceae from the same family as the cauliflower and broccoli (Figure 1a, b). In a recent work [5], we provided a mechanistic explanation of how these mutations alter flower development, by recursively triggering lateral meristem development, which results in the fractal-like shape of the curds. Moreover, our work showed how different cauliflower morphologies can be caused by changes in the regulation of plant growth dynamics. In the following, we explain the main results and contributions of our work. In the first section, we describe how some mutations of the gene regulatory network (GRN) involved in flower development can open the way to the development of iterative structures necessary to produce fractals. Then, we explain how this GRN regulates the growth dynamics of plant axes and induces a chain reaction of growing axes leading to the emergence of fractal curds. We then describe how a simple model of plant architecture development was used to show that the changes in plant growth dynamics can be at the origin of the remarkable Romanesco cauliflower. Finally, we briefly discuss the relevance of our findings in the context of plant domestication.

# 2. The gene regulatory network

Aerial plant meristems are structures composed of undifferentiated cells that are located at the most

apical part of shoots. These meristems produce lateral meristems that can take different identities (e.g. leaf, inflorescence (IM) or flower meristems (FM)) depending on the regulatory genes they express (Figure 1c). The gene expression in meristems is regulated by interactions between master regulatory genes. In particular, in *A. thaliana*, shoot meristems express *TERMINAL FLOWER1* (*TFL1*), a shoot identity gene that prevents meristems from differentiation by repressing flower identity genes, such as *APETALA1* (*API*), *CAULIFLOWER* (*CAL*) and *LEAFY* (*LFY*) [6–11] (Figure 2a). During the vegetative phase of plant development, flower identity genes are inactive and lateral meristems take shoot and leaf identity. However, as plant ages, flowering signals induce the expression of *API*, *CAL* and *LFY* in new lateral meristems (Figure 2b), which turns them into flowers. *API* inhibits shoot identity by repressing *TFL1* in FM [7, 12, 13] and stabilises the expression of floral regulators by forming a positive feedback loop with *LFY* (Figure 2a) [6–10, 14–17]. In *A. thaliana*, the double mutant of *ap1 cal*, produces cauliflower-like structures instead of flowers [3, 4, 6]. Those structures correspond to meristems initiated as floral (with *LFY* expression) but that revert to shoot (with *TFL1* expression) because the *API/CAL-LFY* positive feedback is missing (Figure 2c).

The first attempts to model the floral GRN only based on *TFL1*, *API/CAL* (*API* and *CAL* were included as a single variable because they are homologous genes with redundant functions) and *LFY*, led us to realise that there were missing elements in the network. In particular, those inducing the expression of *TFL1* were lacking. We therefore looked for possible candidate genes and identified additional genes intricately linked to the original master regulators. These genes are the MADS-box genes *SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1*, *AGAMOUS-LIKE24* and *XAANTAL2* (collectively renamed the “SAX” genes). Interestingly, on the one hand, SAX proteins promote the expression of the flowering genes *API* and *LFY* [14, 18–20] and repress *TFL1* forming a complex with *API* [11]. On the other hand, work from our collaborators showed that SAX proteins are also able to activate *TFL1* expression. We updated the GRN with this information (Figure 3), formalised these interactions in a model of ordinary differential equations and analysed both its transitory and stationary properties. These analyses
<figure><img src="image_4.png" /><figcaption>**Figure 1.** WT vs *ap1 cal* morphology. (a) WT plant with an inflorescence meristem producing flower meristems. (b) The cauliflower morphology produced by the overproduction of lateral meristems in the *ap1 cal* mutant. (c) A microscopy image of an inflorescence meristem (IM), surrounded by flower meristems (FM). Photographs credits: Jan Trass (c) and Marie Le Masson (a and b).</figcaption></figure>

showed that the GRN was sufficient to produce the expected patterns of gene expression in a wild-type (WT) plant, in an *ap1 cal* mutant, and in the other single mutants of the different genes present in the network. In addition, the model predicted a high level of SAX expression in *ap1 cal* meristems that is not observed in WT plants and we experimentally confirmed this prediction with the *SOCI* gene [5].

Having found a minimal GRN able to recapitulate the different gene expression states and explain the shoot identity of new lateral meristems in the *ap1 cal* mutant, we wondered whether it was sufficient to explain the cauliflower inflorescence morphology.

To address this question, we first took into account that the plant architecture does not solely depend on the fate of its meristems, but also on additional developmental factors such as phyllotaxis of the shoots, the growth rate of the internodes, the time it takes for newly created meristems to produce new lateral organs, etc. Moreover, we also noticed that the GRN final states for *lfy* and *ap1 cal* mutant simulations are identical with a high level of TFL1 and SAX and a low level of AP1/CAL and LFY, indicating that the architecture of the plant cannot be solely explained by the GRN stationary states. To understand the origin of curd morphology, we were thus led to analyse in greater depth the dynamics of the different actors in the growing branching structure. For this, we designed a 3D computational model simulating the growth of a plant whose development dynamics

was regulated by the state of GRNs computed in each meristem. We call this model, the architectural model.

**3. The architectural model**

More precisely, our architectural model comprises two main components, meristems and internodes. Each meristem has its own identity based on the state of its GRN (Figure 4a). The identity of a meristem determines its behaviour. Vegetative and inflorescence meristems produce internodes and new lateral primordia (meristems) at each plastochron. Flower meristems turn into flowers, thus ending the axis growth. Vegetative meristems have low *TFL1*, *AP1* and *LFY* expression. Inflorescence meristems have high *TFL1* expression and low *AP1* and *LFY* expression. Finally, flower meristems have low *TFL1* expression and high *AP1* and *LFY* expression. At each step, the model inputs a level of flower inducing signal (that was assumed to increase in each meristem with aging) and a level of auxin that was set maximal in each new meristem and then exponentially decay. The state of the GRN was thus reevaluated in each meristem at each time step and was used to determine the identity of the meristem. In addition, since most meristems of order 3 remained inactive in *A. thaliana*, we added a rule to the model, specifying that meristems of order 3 should be kept dormant, irrespective of their identity (Figure 4b).
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Gene regulatory network and gene expression profiles involved in the flowering transition. (a) Interactions between the major floral regulators. Gene expression during late development in the shoot and lateral meristems of wild-type (WT) (b) and *ap1 cal* mutant (c) plants. SAM = shoot apical meristem. FM = Floral meristem (lateral meristem). In green and red, TFL1 and LFY expression, respectively.</figcaption></figure>

This model, calibrated for WT plant development, was able to recapitulate morphologies of mutated plants with the expected architecture (Figure 4b). For example, *tfl1* mutants terminate with the production of a few flowers due to AP1 and LFY expression in the shoot apical meristem [21], *lfy* mutant produced

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Gene regulatory network of the main genes involved in early floral development (Image from [5]). Arrows represent positive regulations, flat head arrows represent negative regulations and the merging arrows from SAX and AP1 towards TFL1 represent the formation of the SAX-AP1 complex. F and Auxin are inputs of the model. The equations of the model are available in [5].</figcaption></figure>

mostly secondary branches [22]. However, the model could not reproduce the *ap1 cal* morphology. First the *ap1 cal* inflorescence, contrary to *lfy*, is mostly devoid of leaves and no element was accounting for this in the model. Second, the cauliflower phenotype arises from the accumulation of meristems on the top of each other, reaching orders higher than 3 while we had imposed a limitation to the growing orders to model the wild-type plant properly (Figure 4c).

We then realised that, despite the fact that the *lfy* mutant and the *ap1 cal* mutant have identical stationary states in the inflorescence lateral meristems, a critical difference between the two mutants lies in the dynamics to reach that state. While in the *lfy* mutant, *LFY* expression never rises up, a transient peak of *LFY* expression is produced in lateral meristems of *ap1 cal* mutants before they finally take a shoot identity. This key difference was enough to explain both previous issues. Indeed, it was already well documented that the presence of *LFY* is responsible for repressing bracts in primordia [22]. Second, previous work from our group [23] had shown that ectopic expression of *LFY* was able to allow high order meristems to grow. The transient rise of *LFY* in *ap1 cal* was thus essential to free meristem outgrowth at all orders and induce a chain reaction.

In conclusion, our work showed that the cauliflower phenotype is mainly due to the fact that in the *ap1 cal* mutant there is an unsuccessful flower development that produces undifferentiated meristems instead of flowers. The transient peak of
<figure><img src="image_7.png" /></figure>

**Figure 4.** The architectural plant model. (a) Schematic description of the architectural plant model. Green circles represent inflorescence/vegetative meristems. Red circles represent flower meristems. Red lines represent internodes. The simulation starts with a vegetative meristem. Vegetative and inflorescence meristems were able to produce lateral meristem and internode, therefore driving plant growth with time. Each meristem contains its own GRN. The state of the GRN determines its identity. (b) WT morphology produced by the model vs a photograph of a WT *A. thaliana*. The arrows indicate the meristems order. In (c) the left panel shows that the original model produced cauliflowers with meristems of maximal order 3. In the middle, a photograph showing an *ap1 cal* plant with a meristem order higher than 3. On the right panel, a simulation assuming that LFY expression releases the production of meristems of orders higher than 3.

LFY expression in the *ap1 cal* mutant is responsible for the activation of all meristems and causes the recursive production of new meristems at all orders. In other words, the fractal-like cauliflower morphology is the result of an alteration of two fundamental processes shaping plant architecture, namely meristem differentiation and meristem outgrowth repression.

**4. The Romanesco hypothesis**

Brassicas curds have different fractal morphologies. Cauliflowers exhibit a curd with a round dome. Their fractal nature can be observed with a curd

morphology observable at many scales, each curd being formed of smaller curds with similar round domes. This fractal pattern is reinforced by the typical spiral arrangement of smaller curds forming each larger one. The fractal structure of the Romanesco in nested spiralling curds is similar, but the conic shape of curds dramatically highlights its fractal geometry (Figure 5a, b).

We thought of taking advantage of our architectural model to better understand what regulates the development of the different curd morphologies. A previous comparison between cauliflower and Romanesco curds identified several features that are
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Cauliflower and Romanesco morphologies produced by the model and experimentally. (a) Cauliflower and (b) Romanesco morphologies. Curd morphologies produced by the model with (c) constant plastochron vs (d) accelerating plastochron. (e) *ap1 cal* curd morphology vs (f) *clv3 ap1 cal* curd morphology. Photograph credits: Marie Le Masson.</figcaption></figure>

constantly increasing during development in Ro-
manesco curds, but that are constant in cauliflow-
ers [24]. These properties are: (1) the production rate
of new lateral meristems; (2) the number of phyl-
lotaxis spirals in clockwise and counter-clockwise
directions; (3) the number of lateral primordia
produced by a meristem before the first one starts
producing its own lateral primordial (meristem onset

delay); and (4) the size of the meristems. Whether
some of these parameters are causal to the Ro-
manesco phenotype had remained elusive. However,
we reasoned that according to phyllotaxis stud-
ies [25–28], these properties must be linked to each
other. Consequently, we wondered if the morpho-
logical curd variations could also be explained by
the difference in any parameter controlling these
properties. To check this, we constructed a simplified version of our architectural model, making it possible to test the consequences of changing these parameters on the final curd morphology. Strikingly, we found that an increase in the meristem production rate was sufficient to produce conspicuous Romanesco-like shape whereas constant values of this parameter produce a normal cauliflower morphology (Figure 5c, d). This acceleration in the meristem production rate induces the insertion of more internodes in the meristem's axis. Assuming that this process does not affect the elongation law of the internodes (internode elongation was kept unchanged in the simple model), this acceleration of the meristem production rate indirectly induces an acceleration of overall stem elongation relative to its lateral axes, which results in conspicuous conical structures with very visible spirals themselves made of cones and so on.

According to our model, increasing the production rate of lateral meristems should transform a cauliflower into a Romanesco curd. To test this hypothesis in planta, we used the fact that, according to phyllotaxis models, meristem production rate may in principle be augmented by an increase in the size of the central zone. In A. thaliana, the clavata (clv) mutants (e.g. clv1 and clv3) are known to produce bigger meristems [29, 30]. Thus, we crossed a clv3 mutant with the ap1 cal double mutant. Indeed, the clv3 ap1 cal triple mutant has a modified curd shape, with a morphology being more pyramidal than for ap1 cal curds (Figure 5e, f). This result suggests that the meristem production rate is a fundamental regulator of the final curd morphology. Interestingly, CLV3 and a couple of other related genes seem to have a lower expression in Romanesco curds than in cauliflowerers [5], but the causal link between this observation and the curd phenotype needs to be studied.

This mechanism provides a first evidence that the strikingly different morphologies observed in curds could be due to simple change in a single parameter, such as the size of the meristem central zone. This changes meristem production rate, which is then propagated throughout the whole plant as a change in the relative elongation rates of the different axes, resulting in a dramatic change of the curd global morphology.

Altogether, our results provide a plausible explanation for the appearance of fractal-like structures in

plants and for the possibility to significantly modify their geometry based on variations of only a limited number of growth regulation factors.

**5. Discussion**

The pervasive presence of fractal-like structures in biology is intriguing. Until now, the developmental origin of most biological fractal-like structures was still unknown. However, deciphering their origins is a non-trivial task due to the high interconnectivity and non-linearity of developmental processes, which take place at multiple scales. In our work, we used modelling and experimental tools to study the origin of fractal-like cauliflower structures. Our results showed that the fractal structure of curds is very likely the consequence of three main mechanisms. First, the mutation in the *ap1 cal* genes produces an unsuccessful flowering transition. As a consequence, all meristems take a vegetative or inflorescence identity. The *ap1 cal* mutation also releases the production of lateral meristems at all orders. This is due to the transitory expression of *LFY* in the mutant. Finally, a Romanesco-like form of the curds can be induced when extra mutations induce an accelerated production of new lateral meristems. Our work provides a first plausible explanation of how fractal-like structures develop linking gene activity with macroscopic structures.

Our work also provides a framework to better understand fundamental aspects of plant development and plant domestication. The unexpected discovery that a transient exposure to LFY can have a profound effect on how the meristem responds to general cues such as those limiting the growth of high-order structures is not understood at the molecular level and it would be interesting to identify which LFY target genes exert this function.

Moreover, even if our work pointed towards flower and inflorescence regulators that can be responsible for the development of curds, the mutations that were selected during the long process of cabbage, broccoli, cauliflower and Romanesco domestication still need to be identified [31] and validated using genetic and functional analyses. Genome comparisons are a powerful tool, but yield many possible candidates that until now, except for mutations affecting the *CAL* gene, have not yet been validated [32]. We believe that our models can provide a very useful framework for their identification.
# Conflicts of interest

Authors have no conflict of interest to declare.

# References

[1] B. Mandelbrot, "Is nature fractal?", *Science* **279** (1998), no. 5352, p. 783.

[2] H. O. Peitgen, H. Jurgens, D. Saupe, *Chaos and Fractals: New Frontiers of Science*, vol. 2(864), Springer, New York, 1992, 62-63 pages.

[3] S. A. Kempin, B. Savidge, M. F. Yanofsky, "Molecular basis of the cauliflower phenotype in Arabidopsis", *Science* **267** (1995), no. 5197, p. 522-525.

[4] C. Ferrándiz, Q. Gu, R. Martienssen, M. F. Yanofsky, "Redundant regulation of meristem identity and plant architecture by FRUITFULL, APETALA1 and CAULIFLOWER", *Development* **127** (2000), no. 4, p. 725-734.

[5] E. Azpeitia, G. Tichtinsky, M. Le Masson, A. Serrano-Mislata, J. Lucas, V. Gregis *et al.*, "Cauliflower fractal forms arise from perturbations of floral gene networks", *Science* **373** (2021), p. 192-197.

[6] J. L. Bowman, J. Alvarez, D. Weigel, E. M. Meyerowitz, D. R. Smyth, "Control of flower development in Arabidopsis thaliana by APETALA1 and interacting genes", *Development* **119** (1993), no. 3, p. 721-743.

[7] O. J. Ratcliffe, D. J. Bradley, E. S. Coen, "Separation of shoot and floral identity in Arabidopsis", *Development* **126** (1999), no. 6, p. 1109-1120.

[8] S. J. Liljegren, C. Gustafson-Brown, A. Pinyopich, G. S. Ditta, M. F. Yanofsky, "Interactions among APETALA1, LEAFY, and TERMINAL FLOWER1 specify meristem fate", *Plant Cell* **11** (1999), no. 6, p. 1007-1018.

[9] C. M. Winter, R. S. Austin, S. Blanvillain-Baufumé, M. A. Reback, M. Monniaux, M. F. Wu, Y. Sang, A. Yamaguchi, N. Yamaguchi, J. E. Parker, F. Parcy, S. T. Jensen, H. Li, D. Wagner, "LEAFY target genes reveal floral regulatory logic, cis motifs, and a link to biotic stimulus response", *Dev. Cell* **20** (2011), no. 4, p. 430-443.

[10] K. Kaufmann, F. Wellmer, J. M. Muiño, T. Ferrier, S. E. Wuest, "Orchestration of floral initiation by APETALA1", *Science* **328** (2010), no. 5974, p. 85-89.

[11] C. Liu, Z. W. Teo, Y. Bi, S. Song, W. Xi, X. Yang, Z. Yin, H. Yu, "A conserved genetic pathway determines inflorescence architecture in Arabidopsis and rice", *Dev. Cell* **24** (2013), no. 6, p. 612-622.

[12] O. J. Ratcliffe, I. Amaya, C. A. Vincent, S. Rothstein, R. Carpenter, E. S. Coen *et al.*, "A common mechanism controls the life cycle and architecture of plants", *Development* **125** (1998), no. 9, p. 1609-1615.

[13] D. Goretti, M. Silvestre, S. Collani, T. Langenecker, C. Méndez, F. Madueno *et al.*, "TERMINAL FLOWER1 functions as a mobile transcriptional cofactor in the shoot apical meristem", *Plant Physiol.* **182** (2020), no. 4, p. 2081-2095.

[14] V. Grandi, V. Gregis, M. M. Kater, "Uncovering genetic and molecular interactions among floral meristem identity genes in Arabidopsis thaliana", *Plant J.* **69** (2012), no. 5, p. 881-893.

[15] E. Moyroud, E. G. Minguet, F. Ott, L. Yant, D. Posé, M. Monniaux *et al.*, "Prediction of regulatory interactions from genome sequences using a biophysical model for the *Arabidopsis* LEAFY transcription factor", *Plant Cell* **23** (2011), no. 4, p. 1293-1306.

[16] K. Goslin, B. Zheng, A. Serrano-Mislata, L. Rae, P. T. Ryan, K. Kwaśniewska, B. Thomson, D. S. Ó'Maoiléidigh, F. Madueño, F. Wellmer, E. Graciet, "Transcription factor interplay between LEAFY and APETALA1/CAULIFLOWER during floral Initiation", *Plant Physiol.* **174** (2017), no. 2, p. 1097-1109.

[17] F. Parcy, O. Nilsson, M. A. Busch, I. Lee, D. Weigel, "A genetic framework for floral patterning", *Nature* **395** (1998), no. 6702, p. 561-566.

[18] C. Liu, H. Chen, H. L. Er, H. M. Soo, P. P. Kumar, J. H. Han, Y. C. Liou, H. Yu, "Direct interaction of AGL24 and SOC1 integrates flowering signals in Arabidopsis", *Development* **135** (2008), no. 8, p. 1481-1491.

[19] J. Lee, M. Oh, H. Park, I. Lee, "SOC1 translocated to the nucleus by interaction with AGL24 directly regulates leafy", *Plant J.* **55** (2008), no. 5, p. 832-843.

[20] R. V. Pérez-Ruiz, B. García-Ponce, N. Marsch-Martínez, Y. Ugartechea-Chirino, M. Villajuana-Bonequi, S. de Folter, E. Azpeitia, J. Dávila-Velderrain, D. Cruz-Sánchez, A. Garay-Arroyo, L. Sánchez Mde, J. M. Estévez-Palmas, E. R. Álvarez Buylla, "XAANTAL2 (AGL14) is an important component of the complex gene regulatory network that underlies Arabidopsis shoot apical meristem transitions", *Mol. Plant* **8** (2015), no. 5, p. 796-813.

[21] S. Shannon, D. R. Meeks-Wagner, "A mutation in the arabidopsis TFL1 gene affects inflorescence meristem development", *Plant Cell* **3** (1991), p. 877-892.

[22] D. Weigel, J. Alvarez, D. R. Smyth, M. F. Yanofsky, E. M. Meyerowitz, "LEAFY controls floral meristem identity in Arabidopsis", *Cell* **69** (1992), p. 843-859.

[23] H. Chahtane, G. Vachon, M. Le Masson, E. Thévenon, S. Périgon, N. Mihajlovic, A. Kalinina, R. Michard, E. Moyroud, M. Monniaux, C. Sayou, V. Grbic, F. Parcy, G. Tichtinsky, "A variant of LEAFY reveals its capacity to stimulate meristem development by inducing RAX1", *Plant J.* **74** (2013), p. 678-689.

[24] M. Kieffer, M. Fuller, A. Jellings, "Explaining curd and spear geometry in broccoli, cauliflower and 'romanesco': quantitative variation in activity of primary meristems", *Planta* **206** (1998), p. 34-43.

[25] F. T. Richards, "Phyllotaxis: its quantitative expression and relation to growth in the apex", *Philos. Trans. R. Soc. Lond. B Biol. Sci.* **235** (1951), p. 509-564.

[26] S. Douady, Y. Couder, "Phyllotaxis as a dynamical self organizing process Part I: The spiral modes resulting from time-periodic iterations", *J. Theor. Biol.* **178** (1996), p. 255-273.

[27] S. Douady, Y. Couder, "Phyllotaxis as a dynamical self organizing process Part II: The spontaneous formation of a periodicity and the coexistence of spiral and whorled patterns", *J. Theor. Biol.* **178** (1996), p. 275-294.

[28] S. Douady, Y. Couder, "Phyllotaxis as a dynamical self organizing process Part III: The simulation of the transient regimes of ontogeny", *J. Theor. Biol.* **178** (1996), p. 295-312.

[29] J. C. Fletcher, U. Brand, M. P. Running, R. Simon, E. M.
Meyerowitz, "Signaling of cell fate decisions by CLAVATA3 in Arabidopsis shoot meristems", *Science* **283** (1999), no. 5409, p. 1911-1914.

[30] M. Kitagawa, D. Jackson, "Control of meristem size", *Annu. Rev. Plant Biol.* **70** (2019), p. 269-291.

[31] C. Cai, J. Bucher, F. T. Bakker, G. Bonnema, "Evidence for two domestication lineages supporting a middle-eastern origin for Brassica oleracea crops from diversified kale populations", *Hortic Res.* **9** (2022), article no. uhac033.

[32] N. Guo, S. Wang, L. Gao, Y. Liu, X. Wang, E. Lai, M. Duan, G. Wang, J. Li, M. Yang, M. Zong, S. Han, Y. Pei, T. Borm, H. Sun, L. Miao, D. Liu, F. Yu, W. Zhang, H. Ji, C. Zhu, Y. Xu, G. Bonnema, J. Li, Z. Fei, F. Liu, "Genome sequencing sheds light on the contribution of structural variants to Brassica oleracea diversification", *BMC Biol.* **19** (2021), no. 1, article no. 93.