Influence of structural reinforcements on the twist-to-bend ratio of plant axes: a case study on Carex pendula

During biological evolution, plants have developed a wide variety of body plans and concepts that enable them to adapt to changing environmental conditions. The trade-off between flexural and torsional rigidity is an important example of sometimes conflicting mechanical requirements, the adaptation to which can be quantified by the dimensionless twist-to-bend ratio. Our study considers the triangular flower stalk of Carex pendula, which shows the highest twist-to-bend ratios ever measured for herbaceous plant axes. For an in-depth understanding of this peak value, we have developed geometric models reflecting the 2D setting of triangular cross-sections comprised of a parenchymatous matrix with vascular bundles surrounded by an epidermis. We analysed the mathematical models (using finite elements) to measure the effect of either reinforcements of the epidermal tissue or fibre reinforcements such as collenchyma and sclerenchyma on the twist-to-bend ratio. The change from an epidermis to a covering tissue of corky periderm increases both the flexural and the torsional rigidity and decreases the twist-to-bend ratio. Furthermore, additional individual fibre reinforcement strands located in the periphery of the cross-section and embedded in a parenchymatous ground tissue lead to a strong increase of the flexural and a weaker increase of the torsional rigidity and thus resulted in a marked increase of the twist-to-bend ratio. Within the developed model, a reinforcement by 49 sclerenchyma fibre strands or 24 collenchyma fibre strands is optimal in order to achieve high twist-to-bend ratios. Dependent on the mechanical quality of the fibres, the twist-to-bend ratio of collenchyma-reinforced axes is noticeably smaller, with collenchyma having an elastic modulus that is approximately 20 times smaller than that of sclerenchyma. Based on our mathematical models, we can thus draw conclusions regarding the influence of mechanical requirements on the development of plant axis geometry, in particular the placement of reinforcements.

Modelling of biological materials systems. Plants are multifunctional structures whose diverse functions are anchored to six hierarchical levels (molecule, organelle, cell, tissue, organ, plant). In this study, we focus on the tissue level and the modelling of the influence of the two-dimensional tissue distribution on the mechanical performance of the entire plant axis. Experiments and findings within the models go far beyond those on plants. For example, the number of fibre bundles can be increased or reduced, even though the amount of material remains the same if the dimensions (diameter) of the bundles change accordingly. Furthermore, in the simulation, lignified sclerenchyma fibres can be replaced by non-lignified collenchyma fibres, whereby the elastic modulus of the sclerenchyma fibres is one order of magnitude higher than that of the collenchyma fibres. In contrast, the non-lignified parenchyma, which often takes up the largest part of the plant cross-section, and the epidermis, a dermal tissue comprising a single layer of cells, have very low elastic moduli, which are on average one to three orders of magnitude smaller than those of the collenchyma fibres 4,8,[11][12][13] . The simulation also Scientific Reports | (2021) 11:21232 | https://doi.org/10.1038/s41598-021-00569-z www.nature.com/scientificreports/ allows us to change a non-lignified single-layered epidermis into a corky multi-layered peridermal tissue, which is at least 10-times as stiff as the epidermis (unpublished data). A cornerstone of plant biomechanics is the performance of mechanical tests under a variety of applied loads such as tension, bending, compression and torsion, both in the linear-elastic range and up to ultimate strength, i.e. at failure. Data of the mechanical tests and the corresponding geometrical properties in terms of size, shape and tissue arrangement are necessary to develop analytical and numerical models. The models not only enable a deeper understanding of the functional principles of the plant model, but are also an indispensable precondition for the transfer of knowledge to technical developments, because models represent a common language for natural scientists, mathematicians, materials scientists and engineers 14 . According to the biomimetic approaches of the "biology push process" or "bottom-up approach" and the "technology pull process" or "top-down approach", every successful biomimetic product has to go through a step of abstraction 15 .
Both analytical and numerical modelling have their intrinsic advantages and short comings when they are used to improve our understanding of the form-structure-function relationship of the biological models or to facilitate the transfer into biomimetic applications. Analytical models, for example, allow fast predictions about the influence of a variety of structural, morphological and/or anatomical changes in a plant organ on its mechanical properties and also enable the inclusion of predictions about non-existing intermediate or extreme forms in the analysis. Similar to the requirements of a "closed" mathematical description, as desired in analytical models, the number of variables and the boundary conditions have to be limited, and a careful simplification in the description of the biological model is typically necessary. The reduction to a few characterising variables allows, on the one hand, a (much) faster analysis of the biological models. On the other hand, this procedure unfortunately includes dangers of oversimplification and of decisive variables being overlooked. The potential of analytical models for an in-depth understanding of the biological model and the transfer to a technical application was demonstrated in the following studies: mechanically driven self-sealing function of the leaves of Delosperma cooperi 16 , tuyere surfaces of metals inspired by lotus leaves 17 , and modelling of stomatal density response to atmospheric CO 2 18 . Numerical models (often) permit the precise description of form and structure of a biological model, but frequently need a plethora of variables. The measurement of all of these variables to the accuracy needed for numerical models is mostly extremely time-consuming (and sometimes even impossible). Moreover, numerical models are only of limited value for the inclusion and prediction of non-existing intermediate or extreme forms in the analysis. The latter is especially true as far as changes in form and structure are involved, whereas variations in mechanical variables are typically easy to include in numerical models. The potential of numerical models for a deeper understanding of the biological model and the transfer to a technical application was demonstrated in the following studies: hydraulically driven self-sealing function of the leaves of Delosperma cooperi 19 , the biomimetic cellular actuator inspired by turgor driven plant movement 20 , and elastic systems in architecture transferred from plant movements 21 .
Now that the advantages and shortcomings of analytical and numerical modelling have been pointed out, it becomes obvious that a combination of both approaches can be considered as the "royal road" of modelling, combining the advantages and avoiding the disadvantages of both approaches. In this work, we have therefore used simplified materials models and have made some geometrical assumptions to be able to focus on the essential variables of the model. The resulting description is too complex to provide a closed form solution but, is simple enough that, by using a finite element approach, we can efficiently analyse a whole range of variable values. Detailed information on the chosen approach is given in section "Mathematical models".
Aim of the study. In the study presented, we selected the triangular flower stalk of Carex pendula as model system, because to our knowledge, it posses the highest twist-to-bend ratios ever measured for herbaceous plant axes. Therefore, it represents a prime example for testing our models as to their predictive strength. We aim to find answers to the following scientific question: "To what extent do individual tissues such as fibres, vascular bundles, epidermis and parenchyma contribute to the flexural rigidity, the torsional rigidity and thus to the twist-to-bend ratio of a triangular plant axis?" Mathematical calculations have therefore been carried out, based on mechanical and geometrical properties from the literature. In this context, the effect of additional reinforcements on the flexural rigidity, the torsional rigidity and the twist-to-bend ratio have been examined with regard to: (1) the formation of a periderm instead of an epidermis, (2) an increasing number of fibre strands up to an optimum, while keeping their total area in a cross-section constant and (3) the replacement of sclerenchyma fibres by collenchyma fibres, two fibre types that differ notably in their elastic moduli.

Plant data for modelling
Carex pendula with its peak values of twist-to-bend ratio is a very suitable model plant for studying the triangular geometry. The necessary data of geometrical and mechanical properties for the mathematical calculations presented below have been collected in previous studies of the Plant Biomechanics Group Freiburg 2,6,8 . Figure 1 presents a stained section and schematic drawing of its flower stalks showing the cross-sectional distribution of tissues. The triangular cross-section shows an outer epidermis (e). Individual lignified sclerenchyma fibre strands (sc) in the periphery and vascular bundles (vb) scattered within the outer half of the cross-section are embedded in a non-lignified ground tissue consisting largely of parenchyma (pa) with interspaced smaller regions of aerenchyma (ae) close to the periphery and a layer of chlorenchyma (ch) directly underneath the epidermis. Table 1 presents experimental data of the flower stalk of C. pendula. Mechanical properties from bending and torsional tests are provided, as are geometrical properties of the entire internode and of individual tissues.  Table 2 are also included in the mathematical analyses of the influence of fibre reinforcement of plant axes on the twist-to-bend ratio presented below.

Mathematical models
In order to describe the influence of fibre reinforcement to the twist-to-bend ratio of a plant, we can consider a 2D-model determining the flexural and torsional rigidity of a beam its cross-section. In our setting, it does not suffice to simply consider the products of Young's and shear modulus with second and polar moment of area, respectively given by EI and GJ, as these approximations are only valid for nearly circular cross-sections of homogeneous tissue. As in our previous work 24 we use methods from linearized elasticity, which we repeat here for the readers convenience. As we are interested in investigating mechanical properties of the flower stalk of C. pendula we  www.nature.com/scientificreports/ describe a plant stem as a long thin elastic rod with domain B = � × (0, L) of length L and simply connected cross-section remaining constant along the longitudinal axis. We can now consider the domain as a composite of different materials bounded by a sufficiently regular boundary curve ∂� . Further, we assume L ≫ diam and isotropy for the materials involved. Anisotropic effects, viscosity and other time-dependent processes are neglected here since we are only interested in the influence of the cross-sectional geometry and the contained distribution of various materials on the mechanical properties of the stem. For our modelling, we consider B fixed at z = 0 and bending of B to be due to an outer normal force on at z = L. Mora and Müller 25 have mathematically rigorously derived the flexural (or bending) rigidity by considering the limit of a very slender and long rod, confirming the classical approaches used here. Following this classical theory, see Crandall et al. 26 , flexural rigidity can be deduced from the moment curvature relation where M y , M x denote the bending moments applied at the end of the beam ( z = L ), κ x , κ y denote the curvature in the direction of x and y respectively and the moments of inertia D x , D y and the product of inertia D xy are given by where we have considered the coordinate system (x,ŷ) which has its origin at the centroid of the cross-section , i.e, Because of the heterogeneity of the elastic modulus E(x, y) depends on the cross-sectional coordinates, being piece-wise constant. The maximal and minimal flexural rigidities D max and D min along the principal axes are then given by the maximal and minimal eigenvalue of the matrix in (1) leading to with D mean = D x +D y 2 . In the following we are concerned with the problem of generating cross-sections with high flexural rigidity. Therefore, we will confine ourselves to the computation of the minimal flexural rigidity D min . High values of D min then lead to high resistance against bending forces in any direction orthogonal to the flower stalk.
As for the flexural rigidity, the torsional rigidity for an elastic slender rod with domain B has also been mathematically rigorously derived by Mora and Müller 25 . In order to describe the torsional rigidity, we thus use St.Venant's theory of pure torsion of nonhomogeneous elastic beams, which has been employed, among others, by Ecsedi 27 . Torsion is assumed to be due to a moment T at the top of B and thus the torsional rigidity can be expressed by Prandtl's stress function φ(x, y) satisfying with the shear modulus G(x, y) depending on the cross-sectional coordinates, being piece-wise constant. By using the stress function φ , the torsional rigidity is then given by  (2,4) we have employed a P1 triangular finite element discretisation of the cross-section . To be precise, was chosen as an equilateral triangle and discretized using approximately 5 · 10 6 triangular elements, thus finely resolving the material heterogeneities due to fibre reinforcement (see, e.g., Fig. 6a). The implementation of this standard finite element method (C++-code) is available in the supplementary material.

Effect of reinforcements on the twist-to-bend ratio
During their ontogeny, plants react to increasing bending forces triggered by continuous growth of stem length and the formation of top loads such as flowers, seeds and fruits. An increase of their resistance against bending forces can be achieved by an increase of their flexural rigidity along the principal axes, i.e. increasing D min . For a better understanding of this effect, we introduce a simplified model that is related to C. pendula and that measures flexural and torsional rigidity. As a first approach, we describe the tissue arrangement in the cross-section of C. pendula in terms of a distribution of circles (vascular bundles) in an equilateral triangular reference domain additionally filled with parenchyma and surrounded by an epidermis (see Fig. 2a). The contribution to the crosssectional area of the vascular bundles ( ≈ 29% ) and the parenchyma ( ≈ 52% ) and the thickness of the epidermis is selected to match the distribution of mass in the total cross-sectional area of C. pendula.
For the computation of the rigidities, we normalise the elastic moduli ( E est ) of all contained materials with respect to the elastic modulus of the sclerenchyma. This means that we set E norm = 1 for the sclerenchyma and obtain the normalised elastic moduli of the other materials by scaling accordingly (see Table 1). Furthermore, we consider the ratio between normalised elastic moduli of the stiffest and the most elastic material contained in the cross-section. As torsional and flexural rigidity are almost exclusively determined by these two materials, they are referred to as the mechanically decisive materials. Accordingly, the ratio between mechanically decisive materials is denoted by µ in the following.
We assume a constant Poisson's ratio ν for the materials involved and compute the normalised torsional modulus G norm as for a given elastic modulus E norm . The assumption of a constant Poisson's ratio is reasonable as the value range ν ∈ [0.2, 0.5] is typical for many plant axes 28 and, thus, a change in ν among the materials is negligible for our model. In the following, we set ν = 0.35 for C. pendula. Reinforcement without fibres: from epidermis to periderm. In this section, we compare the effect of a single-layered epidermis and a multi-layered corky periderm on the twist-to-bend ratio of the triangular cross-section (see Fig. 2). In this configuration and with respect to the normalised elastic modulus of the epidermis ( E norm = 0.00111 ), the torsional rigidity ( D z ≈ 0.00013 ) and the flexural rigidity ( D min ≈ 0.00118 ) are comparatively low and so is the twist-to-bend ratio ( D min D z ≈ 9.37 ) (see Fig. 2a). The resistance against bending forces and, hence, the flexural rigidity increases by the formation of a ring of corky periderm ( E norm = 0.01111 ). This formation of a closed ring of strengthening tissue, however, simultaneously increases the flexural rigidity Reinforcement by sclerenchyma fibre strands. In addition to the above-mentioned reinforcement of the epidermal tissue, fibre reinforcement is extremely common in plant axes. The cross-section of C. pendula shows individual sclerenchyma strands in the periphery directly under the epidermis (Fig. 1). Therefore, we additionally incorporated sclerenchyma fibre strands ( E norm = 1 ) into the parenchyma ( E norm = 0.00044 ) of the cross-section from Fig. 2a. We have developed a model to describe the effect of fibre reinforcement. The mechanically decisive materials for this experiment are parenchyma and sclerenchyma. The ratio µ between their elastic moduli is given by µ = 0.00044 . We fix the proportion of the sclerenchyma ( ≈ 17% ; see Table 1) in the total cross-sectional area. Further, we consider a distribution of fibre bundles (sclerenchyma) around the inner boundary ∂� of with circular cross-sections centred at a fixed distance to ∂� . Starting with 6 fibre strands for each side of the triangle we refine the distribution, such that a higher number of fibre strands is used in each step of the refinement, see Fig. 3. The proportion of the fibre bundles in the total cross-sectional area is fixed during the whole process. The arrangement and the total cross-sectional area of the fibre strands are now determined by the structure of the Carex ground tissue, see Fig. 1, and, hence, the number of fibre bundles is the only free variable in this model. The procedure is stopped before single fibre strands become connected, as such a closed ring of sclerenchyma would immediately (markedly) increase the torsional rigidity and thus decrease the twist-to-bend ratio. For simplicity, circular cross-sectional geometries of the fibre bundles are taken, as circular cross-sections exhibit high torsional rigidity and, therefore, any other cross-sectional geometries of the fibre strands is likely to amplify the effect of decreasing torsional rigidity while increasing the number of fibre bundles.
The numerical experiments show that, by increasing the number of fibre strands, the torsional rigidity D z initially decreases when the flexural rigidity remains nearly bounded and, thus, the twist-to-bend ratio increases (see Fig. 4). A regression analysis for torsional rigidity up to 51 fibre strands shows that D z is given as a function of the number of fibre bundles N, with for constants a, C ∈ R . Therefore, increasing the number of fibre bundles initially decreases the torsional rigidity scaling as N −1 in dependence on the number of fibre bundles with asymptote C. This decline in the torsional rigidity is driven by a decreasing amplitude of Prandtl's stress function around fibre strands, see Fig. 6c and d. Because of the boundedness of the flexural rigidity, the twist-to-bend ratio increases showing a similar asymptotic behaviour and scaling as N in dependence on the number of fibre strands.
However, when the distance between neighbouring fibre strands becomes very small, even when the fibre strands are not yet connected, see (a) in Fig. 6, the ring-like structure of the fibre arrangement causes the torsional rigidity to increase again leading to a decrease in the twist-to-bend ratio. This occurs because the gradient of φ in the space between the fibre strands increases as the distance between fibre strands decreases, similar to the behaviour of a Neumann sieve 29 . Thus, the value of φ in the inner part of the cross-section is raised, see (e) and (f) in Fig. 6, resulting in an increase of the torsional rigidity when the number of fibres increases beyond 49 (see Fig. 4a) even if no closed ring is formed yet.
Since this increase is integrated over a large domain (the inner part of the triangle) the effect becomes more dominant for smaller distances of neighbouring fibre strands and therefore leads to an increase in the torsional rigidity, see Fig 4a. Thus in the setting described above, torsional rigidity and the twist-to-bend ratio reach an optimum for 49 fibre strands, see Fig. 4a,c. This characterises the number of fibre bundles N as a design variable in an optimisation problem in order to maximise the twist-to-bend ratio and simultaneously minimise the torsional rigidity.   investigate the dependency of the optimal number of fibre strands needed to increase the twist-to-bend ratio on the ratio µ between the elastic moduli of the mechanically decisive materials involved. To do so, we replace the sclerenchyma fibres ( E norm = 1 ) shown in Fig. 3 by collenchyma fibres ( E norm = 0.05555 ), see Table 2. Since now parenchyma ( E norm = 0.00044 ) and collenchyma are the decisive materials, this leads to a ratio µ ≈ 0.008 . We can again carry out the experiment shown in Fig. 3 but with collenchyma fibres instead of sclerenchyma fibres and with the normalised elastic modulus E norm = 0.05555 for the collenchyma fibres. The elastic moduli of the other materials involved (epidermis, parenchyma, vascular bundles) remain the same. Again, we find that up to a certain number of collenchyma strands, the torsional rigidity decreases linearly, whereas the flexural rigidity remains nearly constant, see Fig. 5a,b. However, now that the ratio µ is lower, the number of fibre strands that is optimal in order to decrease torsional rigidity and to increase the twist-to-bend ratio, is reduced noticeably to 24 and,-after a slight drop-the twist-to-bend ratio remains nearly constant up to 36. For more than 24 collenchyma strands, we can see the same effect as above, namely the torsional rigidity increases again and, hence, the twist-to-bend ratio decreases, see Fig. 5a,c.
This experiment, in which the sclerenchyma fibres are replaced by collenchyma fibres, illustrates the effect of the ratio µ between the elastic moduli of the mechanically decisive materials on the twist-to-bend ratio and on the optimal number of fibre strands included in fibre reinforcement. As a preview to the discussion, we can conjecture that the ratio µ plays a role in the sense that, for the higher ratio µ , the optimal number of fibre bundles in order to increase the twist-to-bend ratio is lower and vice versa. When collenchyma fibres are incorporated  . Reinforcement by hypothetically occurring collenchyma fibre strands. Trend of (a) torsional rigidity, (b) flexural rigidity and (c) twist-to-bend ratio. The ratio µ between the elastic moduli of the parenchyma and collenchyma is µ = 0.008 . The collenchyma and parenchyma are the mechanically decisive materials. By increasing the number of fibre strands used for reinforcement, the torsional rigidity first decreases nearly linearly, see (a), whereas the flexural rigidity remains nearly constant during the procedure, see (b). The twistto-bend ratio first increases nearly linearly reaching a maximum for N = 24 . For more than 24 strands, the torsional rigidity is caused to increase again, whereas the twist-to-bend ratio is reduced by this effect, see (a,c).

Discussion
When embryophytic plants colonised the land and lost the buoyancy of surrounding water around 470 million year ago during the mid-Ordovician, they faced entirely different mechanical constraints. These concerned, in particular, their anchorage (root system) and the mechanical loads on their upright aerial parts (stems with branches and leaves). Stems, roots and leaves had evolved as early as the mid Devonian (ca. 390 million years BP) and the first forests of tall trees existed by the late Devonian (ca. 370 million years BP) 30,31 . Since then, evolutionary processes have shaped stem and root form and internal structure enabling plants to cope with the various mechanical loads that act upon them in their diverse environments. A plethora of evolutionary adaptations can be observed in aerial stems, especially with regard to bending and torsional loading, which represent the predominant load cases in plants with self-supporting upright stems. In addition to experimental analyses, theoretical considerations including analytical and numerical simulations can help to decipher the complex interplay between form, structure and mechanical properties of the involved plant tissues 24,[32][33][34][35] .
In this article, we present our finite element analysis involving a triangular cross-section consisting of parenchyma with scattered embedded vascular bundles surrounded by an epidermis (see Fig. 2a). In contrast to real experiments, computer experiments allow one particular variable to be altered and the resulting effect to be quantified precisely. By replacing the epidermis with a periderm or by including additional fibre strands with variable numbers and/or mechanical quality, we can quantify the mechanical effect of structural reinforcement on the flexural and torsional rigidity and thus the twist-to-bend ratio of plant axes. The biological model that we have used for our simulation is the triangular cross-section of the flower stalk of C. pendula, which shows the so far highest twist-to-bend ratio measured in plant stems. The long, drooping stem attains lengths of up to 2.5 m and bears apically pendulous inflorescences that develop into a relatively heavy multiple fruit. The presented approach gives us the possibility to give answers to our scientific question "To what extent do individual tissues such as fibres, vascular bundles, epidermis and parenchyma contribute to the flexural rigidity, the torsional rigidity and thus to the twist-to-bend ratio of a triangular plant axis?" Table 3 gives an overview on the results, which will be discussed in detail in the following.
From epidermis to periderm. The outermost single cell layer of plants is called the epidermis, the main function of which is to protect the plant. This primary covering with an elastic modulus of E ≈ 50 MPa can support tensile loads. Woody stems can also produce a secondary protective covering, called the periderm, which replaces the epidermis locally or globally. The periderm has an elastic modulus of E ≈ 500 MPa (unpublished data). A triangular cross-section comprising of parenchyma and vascular bundles surrounded by an epidermis (see Fig. 2a) shows relatively low values of the flexural and torsional rigidity. If the epidermis is replaced by a periderm, which forms a closed ring (see Fig. 2b), the flexural rigidity increases by ≈ 14% and the torsional rigidity increases by ≈ 46% . This results in a decrease of the twist-to-bend ratio of ≈ 24% in cases with the formation of periderm (see Table 3).
The cactus O. ficus-indica shows pronounced formation of periderm at its stem base at which the greatest bending moments occur, because of the long lever arm, and high torsional forces can appear because of the firm clamping by rooting in the soil. In addition, increased periderm formation can be found at the narrowed crosssections of the transitions between its branches, at which point an increased risk of breakage through bending or torsion forces by wind loads and passing animals exists 7 . If exposed to moderate wind, tall and shrubby cacti do not become streamlined, but start to oscillate. Under these conditions, high torsional rigidity together with high bending rigidity is advantageous, because it prevents damage and breakage. We therefore hypothesise that periderm formation takes place at mechanically highly stressed or damage-prone areas (e.g. at the basis of the entire plant, the transition between the branches and the main stem and the branches). Periderm formation might be triggered by tiny injuries arising from overcritical local mechanical stress. If the cacti are exposed to heavy wind, single branches are torn off or the whole plant is uprooted.
The formation of an outer secondary protection tissues (secondary cortex = periderm) covering the epidermis in older ontogenetic stages was probably an evolutionary important step allowing an increase in stem girth Table 3. Summary of the simulation results. The mechanical effect of various changes of the structural reinforcement is given as a percentage increase or decrease with reference to the initial configuration (= reference value). Figure 2a shows the initial configuration in terms of a triangular cross-section comprising of parenchyma with embedded vascular bundles surrounded by an epidermis. µ is the ratio of estimated elastic moduli (E est ) of the decisive materials (see Table 2). The percentage values inherent to the fibre reinforcements are computed with respect to the optimal number of fibre strands. namely, *49 sclerenchyma strands and **24 collenchyma strands (see Figs. 4,5).  This because, as the distance between the fibre strands becomes smaller, the gradient of φ between two fibre strands is increased and, therefore, the value of φ in the inner part of the cross-section is raised. This effect becomes dominant from a certain number of fibre strands and, thus, increases the torsional rigidity. Further, the impact of these effects on the torsional rigidity depends on the ratio µ between the elastic moduli of the two mechanically decisive materials. For high ratios µ , the impact of the effect in (e) and (f) is more pronounced, whereas for low ratios, the impact of the effect in (c) and (d) is more intense. The fineness of the triangulation is sufficiently high, so that even small gaps between fibre strands can be resolved, see (a). www.nature.com/scientificreports/ because of the production of secondary xylem and phloem by secondary vascular growth. Some of the first fossil records of secondary cortex originate from the late Devonian/early Carboniferous members of the progymnosperms, lycopsids and pteridosperms 31,[34][35][36] . Periderm formation not only allows secondary stem growth, which increases girth, through the simultaneous growth of an outer secondary protection layer, but also may take over mechanical functions [36][37][38] and ensure repair after injuries 7 . Whereas, in extant perennial arborescent growth forms, the mechanical properties and stability are governed by secondary wood (= secondary xylem), this was not the case in a group of important forest-forming plants during the late Devonian and Carboniferous: the arborescent lycopsids. In this group with scale trees (genus: Lepidodendron) and seal trees (genus: Sigillaria), the periderm, which contributed over 90% of the stem volume, was by far the mechanically dominant tissue 31,32,36,39 . Peripheral stiffening structures were the first specialised stabilising tissues evolved in early land plants. After the colonisation of the land by early plant genera with turgorstabilised stems, specialised collenchymatous or sclerenchymatous stiffening tissues, the so-called hypodermal steromes lying directly underneath the epidermis, evolved as early as the lower Devonian. These hypodermal steromes allowed not only the colonisation of dryer habitats, but also a marked increase in plant height 30-33,36 . Reinforcement by sclerenchyma fibre strands. As documented in Table 2, sclerenchyma fibres are the stiffest material in our simulation with an elastic modulus of E ≈ 45 GPa. Thus, we selected it as our reference material ( E norm = 1 ). We added fibre strands to the initial configuration, with the total area of the strands remaining constant. In this simulation, we investigated the effect of an increasing number of individual peripheral sclerenchyma strands in a parenchymatous matrix with vascular bundles being surrounded by an epidermis (see Fig. 3). The flexural rigidity remained almost constant independent of the number of sclerenchyma strands, namely between 18 and 63 strands (see Fig. 4b). In contrast, the torsional rigidity decreased nearly linearly from 18 to 49 strands. From 50 to 63 strands, it increased again to ≈ 140% of the initial value (see Fig. 4a). Since the flexural rigidity remained almost constant, the twist-to-bend ratio initially increased almost linearly up to an optimum at 49 strands and then decreased to ≈ 2200% of the initial value (see Fig. 4c). Compared with the initial configuration, the addition of 49 sclerenchyma strands led to enormous increases in the flexural and torsional rigidity and the twist-to-bend ratio (see Table 3).
This optimum of N = 49 strands found in the simulation corresponds well with the average value of 49.40 ± 7.83 , which is calculated from the numbers of strands of the apical internode ( 42.80 ± 4.31 ) and the more basal internode ( 56.00 ± 3.87 , see Table 1) of C. pendula 2 . The clear optimum of the U-shaped curve (see Fig. 4a) is of interest in so far that sclerenchyma fibres are energetically highly costly for plants because of their extremely thick secondary cell walls that are impregnated with the macromolecule lignin. For annual plants, in particular, this is an expensive investment in a strengthening tissue that is, however, highly rigid. Since sclerenchyma fibres are dead cells, no further investment is required once they have been formed. This is especially advantageous for perennial plants.
With twist-to-bend ratios of up to 403, the internodes of C. pendula show the highest values ever measured in plant axes 2 . The high flexural rigidity and low torsional rigidity of the flower stalk are particularly advantageous under dynamic wind loads. The high torsional flexibility allows the streamlining of the stalks together with the apical pendulous inflorescences or multiple fruits. In summer, the flower stalks have to bear these additional heavy fruits without bending down to the ground. Our simulations show that flexural rigidity cannot be achieved by the formation of additional sclerenchyma strands, as the number is already in the range of the optimum. However, a further increase in flexural rigidity is ensured by the rigid leaf sheaths that enclose large parts of the flower stalk. Since the leaf sheaths are not firmly attached to the stem, they are unlikely to have a strong negative effect on torsional flexibility. Reinforcement by collenchyma fibre strands. Collenchyma fibres are living cells that possess a vacuole and a thick non-lignified primary cell wall. Therefore, their elastic modulus is turgor-dependent. With an elastic modulus of E ≈ 2.5 GPa, collenchyma fibres belong to the strengthening tissues ( E norm = 0.05555 ) and are often found in still-growing shoots and leaves. From an energy point of view, collenchyma fibres are much cheaper to build than sclerenchyma fibres, because neither thick secondary cell walls nor the macromolecule lignin need to be formed in the former. Furthermore, collenchyma fibres have the advantage that they are still able to grow or to be stretched. On the other hand, as living cells, they always consume physiological energy and hence, are predominantly found in young growing tissues and in annual plants or plant organs.
The square flower stalks of the perennial L. cardiaca show pronounced strands of collenchyma fibres at the four corners and small strands in the middle of the four sides. Dependent on the height above ground and increasing age, the cross-sectional percentage area of collenchyma decreases from 8.61% to 6.18% . With increasing age and mass attributable to the formation of flowers and heavy fruits, the peripheral lignified vascular tissues with a percentage cross-sectional area increasing from 15.55% to 21.00% increasingly take over the mechanical support of the stem. Analyses have demonstrated that the area sum of the vascular tissues and parenchyma exhibit moderate positive allometric scaling, whereas the collenchyma shows clear negative scaling 8 .
In the circular leaf stalks of C. bicolor 'Candyland' , the cross-sectional percentage area of the ≈ 40 peripheral collenchyma strands has a value of 3.9 ± 0.7% . The cross-sectional area of the parenchyma and the lignified elements are 89.0 ± 3.3% and 1.2 ± 0.4% , respectively. Parenchyma and collenchyma are turgor-dependent. Together, they form a doubly secured mechanical system that is sensitive to drought stress. The decrease of flexural rigidity and, thus, the wilting of the leaf stalk are the result of a turgor-loss-induced decrease of the elastic moduli of both the collenchyma fibres and the parenchyma cells 5 . As a withered leaf stalk cannot be restored to its healthy positioning, even with sufficient water support, the evolution of a redundant mechanical system to maintain the flexural rigidity of the plant, in particular, is of great advantage for selection.

Scientific Reports
| (2021) 11:21232 | https://doi.org/10.1038/s41598-021-00569-z www.nature.com/scientificreports/ In the case of C. pendula, the replacement of sclerenchyma fibres by collenchyma fibres showed an increase of flexural rigidity by only ≈ 315% instead of ≈ 5600% (see Table 3). Reinforcement by collenchyma fibre strands would be much too weak mechanically to support the stalk. Moreover, even the addition of an optimal number of sclerenchyma strands is probably insufficient. During autumn, in particular, when the flower stalk has a heavy top load of fruits, additional leaf sheaths increase the flexural rigidity of the overall system.

Conclusion
We have shown here that the ratio of the elastic moduli of the materials, which are decisive for the mechanical performance of the entire plant axis ( µ ), plays a crucial role in plant stems. Reinforcements generally increase the flexural and torsional rigidities in a triangular cross-section composed of a parenchymatous matrix with embedded scattered vascular bundles surrounded by an epidermis. Closed ring-shaped reinforcement of the epidermal tissue (e.g. a periderm) leads to a considerable increase in torsional rigidity and a moderate increase in flexural rigidity in areas of the plants that are mechanically heavily loaded or at risk of damage. Therefore, epidermal reinforcements decrease the twist-to-bend-ratio. Fibre reinforcement noticeably increases the flexural rigidity and moderately increases the torsional rigidity of the entire plant axis. The flexural rigidity is almost independent of the number of fibre strands, whereas the torsional rigidity and thus the twist-to-bend ratio is a function of strand numbers. Obviously, torsional rigidity is the key factor for changing the twist-to-bend ratio through structural reinforcements. The evolution of structural reinforcements including cells with walls strengthened by lignin was a prerequisite for land plants to be able to colonise terrestrial habitats and successfully enabled them to face notably different mechanical constraints compared with those experienced in an aquatic environment. As outlined above, the primary (e.g. hypodermal sterome) and secondary (e.g. wood, cortex, periderm) strengthening tissues allowed plants not only to increase markedly in height, but also to colonise more and more hostile (dryer) habitats leading to the plethora of plant life forms that we know from extant and fossil flora [30][31][32]36 .

Data availibility
This work includes no experimental data. A C++-implementation including the computation of rigidities is available in the supplementary material.