Combinatorial Maps, a New Framework to Model Agroforestry Systems.

Agroforestry systems are complex due to the diverse interactions between their elements, and they develop over several decades. Existing numerical models focus either on the structure or on the functions of agroforestry systems. However, both of these aspects are necessary, as function influences structure and vice versa. Here, we present a representation of agroforestry systems based on combinatorial maps (which are a type of multidimensional graphs), that allows conceptualizing the structure-function relationship at the agroecosystem scale. We show that such a model can represent the structure of agroforestry systems at multiple scales and its evolution through time. We propose an implementation of this framework, coded in Python, which is available on GitHub. In the future, this framework could be coupled with knowledge based or with biophysical simulation models to predict the production of ecosystem services. The code can also be integrated into visualization tools. Combinatorial maps seem promising to provide a unifying and generic description of agroforestry systems, including their structure, functions, and dynamics, with the possibility to translate to and from other representations.

Coupling Plant Growth Models and Pest and Disease Models: An Interaction Structure Proposal, MIMIC.

Coupling plant growth model with pests and diseases (P&D) models, with consideration for the long-term feedback that occurs after the interaction, is still a challenging task nowadays. While a number of studies have examined various methodologies, none of them provides a generic frame able to host existing models and their codes without updating deeply their architecture. We developed MIMIC (Mediation Interface for Model Inner Coupling), an open-access framework/tool for this objective. MIMIC allows to couple plant growth and P&D models in a variety of ways. Users can experiment with various interaction configurations, ranging from a weak coupling that is mediated by the direct exchange of inputs and outputs between models to an advanced coupling that utilizes a third-party tool if the models' data or operating cycles do not align. The users decide how the interactions operate, and the platform offers powerful tools to design key features of the interactions, mobilizing metaprogramming techniques. The proposed framework is demonstrated, implementing coffee berry borers' attacks on Coffea arabica fruits. Observations conducted in a field in Sumatra (Indonesia) assess the coupled interaction model. Finally, we highlight the user-centric implementation characteristics of MIMIC, as a practical and convenient tool that requires minimal coding knowledge to use.

Breakups of Chitosan microcapsules in extensional flow.

The controlled rupture of a core-shell capsule and the timely release of encapsulated materials are essential steps of the efficient design of such carriers. The mechanical and physico-chemical properties of their shells (or membranes) mainly govern the evolution of such systems under stress and notably the link between the dynamics of rupture and the mechanical properties. This issue is addressed considering weakly cohesive shells made by the interfacial complexation of Chitosan and PFacid in a planar extensional flow. Three regimes are observed, thanks to the two observational planes. Whatever the time of reaction in membrane assembly, there is no rupture in deformation as long as the hydrodynamic stress is below a critical value. At low times of complexation (weak shear elastic modulus), the rupture is reminiscent of the breakup of droplets: a dumbell or a waist. Fluorescent labelling of the membrane shows that this process is governed by continuous thinning of the membrane up to the destabilization. It is likely that the membrane shows a transition from a solid to liquid state. At longer times of complexation, the rupture has a feature of solid-like breakup (breakage) with a discontinuity of the membrane. The maximal internal constraint determined numerically marks the initial location of breakup as shown. The pattern becomes more complex as the elongation rate increases with several points of rupture. A phase diagram in the space parameters of the shear elastic modulus and the hydrodynamic stress is established.

Two decades of research with the GreenLab model in agronomy.

BACKGROUND: With up to 200 published contributions, the GreenLab mathematical model of plant growth, developed since 2000 under Sino-French co-operation for agronomic applications, is descended from the structural models developed in the AMAP unit that characterize the development of plants and encompass them in a conceptual mathematical framework. The model also incorporates widely recognized crop model concepts (thermal time, light use efficiency and light interception), adapting them to the level of the individual plant. SCOPE: Such long-term research work calls for an overview at some point. That is the objective of this review paper, which retraces the main history of the model's development and its current status, highlighting three aspects. (1) What are the key features of the GreenLab model? (2) How can the model be a guide for defining relevant measurement strategies and experimental protocols? (3) What kind of applications can such a model address? This last question is answered using case studies as illustrations, and through the Discussion. CONCLUSIONS: The results obtained over several decades illustrate a key feature of the GreenLab model: owing to its concise mathematical formulation based on the factorization of plant structure, it comes along with dedicated methods and experimental protocols for its parameter estimation, in the deterministic or stochastic cases, at single-plant or population levels. Besides providing a reliable statistical framework, this intense and long-term research effort has provided new insights into the internal trophic regulations of many plant species and new guidelines for genetic improvement or optimization of crop systems.

Internal trophic pressure, a regulator of plant development? Insights from a stochastic functional-structural plant growth model applied to Coffea trees.

BACKGROUND AND AIMS: Using internal trophic pressure as a regulating variable to model the complex interaction loops between organogenesis, production of assimilates and partitioning in functional-structural models of plant growth has attracted increasing interest in recent years. However, this approach is hampered by the fact that internal trophic pressure is a non-measurable quantity that can be assessed only through model parametric estimation, for which the methodology is not straightforward, especially when the model is stochastic. METHODS: A stochastic GreenLab model of plant growth (called 'GL4') is developed with a feedback effect of internal trophic competition, represented by the ratio of biomass supply to demand (Q/D), on organogenesis. A methodology for its parameter estimation is presented and applied to a dataset of 15 two-year-old Coffea canephora trees. Based on the fitting results, variations in Q/D are reconstructed and analysed in relation to the estimated variations in organogenesis parameters. KEY RESULTS: Our stochastic retroactive model was able to simulate realistically the progressive set-up of young plant architecture and the branch pruning effect. Parameter estimation using real data for Coffea trees provided access to the internal trophic dynamics. These dynamics correlated with the organogenesis probabilities during the establishment phase. CONCLUSIONS: The model can satisfactorily reproduce the measured data, thus opening up promising avenues for further applying this original procedure to other experimental data. The framework developed can serve as a model-based toolkit to reconstruct the hidden internal trophic dynamics of plant growth.

Realistic Procedural Plant Modeling from Multiple View Images.

In this paper, we describe a novel procedural modeling technique for generating realistic plant models from multi-view photographs. The realism is enhanced via visual and spatial information acquired from images. In contrast to previous approaches that heavily rely on user interaction to segment plants or recover branches in images, our method automatically estimates an accurate depth map of each image and extracts a 3D dense point cloud by exploiting an efficient stereophotogrammetry approach. Taking this point cloud as a soft constraint, we fit a parametric plant representation to simulate the plant growth progress. In this way, we are able to synthesize parametric plant models from real data provided by photos and 3D point clouds. We demonstrate the robustness of the proposed approach by modeling various plants with complex branching structures and significant self-occlusions. We also demonstrate that the proposed framework can be used to reconstruct ground-covering plants, such as bushes and shrubs which have been given little attention in the literature. The effectiveness of our approach is validated by visually and quantitatively comparing with the state-of-the-art approaches.

Estimating Sink Parameters of Stochastic Functional-Structural Plant Models Using Organic Series-Continuous and Rhythmic Development.

Functional-structural plant models (FSPMs) generally simulate plant development and growth at the level of individual organs (leaves, flowers, internodes, etc.). Parameters that are not directly measurable, such as the sink strength of organs, can be estimated inversely by fitting the weights of organs along an axis (organic series) with the corresponding model output. To accommodate intracanopy variability among individual plants, stochastic FSPMs have been built by introducing the randomness in plant development; this presents a challenge in comparing model output and experimental data in parameter estimation since the plant axis contains individual organs with different amounts and weights. To achieve model calibration, the interaction between plant development and growth is disentangled by first computing the occurrence probabilities of each potential site of phytomer, as defined in the developmental model (potential structure). On this basis, the mean organic series is computed analytically to fit the organ-level target data. This process is applied for plants with continuous and rhythmic development simulated with different development parameter sets. The results are verified by Monte-Carlo simulation. Calibration tests are performed both in silico and on real plants. The analytical organic series are obtained for both continuous and rhythmic cases, and they match well with the results from Monte-Carlo simulation, and vice versa. This fitting process works well for both the simulated and real data sets; thus, the proposed method can solve the source-sink functions of stochastic plant architectures through a simplified approach to plant sampling. This work presents a generic method for estimating the sink parameters of a stochastic FSPM using statistical organ-level data, and it provides a method for sampling stems. The current work breaks a bottleneck in the application of FSPMs to real plants, creating the opportunity for broad applications.

Stochastic modelling of tree architecture and biomass allocation: application to teak (Tectona grandis L. f.), a tree species with polycyclic growth and leaf neoformation.

Background and aims: For a given genotype, the observed variability of tree forms results from the stochasticity of meristem functioning and from changing and heterogeneous environmental factors affecting biomass formation and allocation. In response to climate change, trees adapt their architecture by adjusting growth processes such as pre- and neoformation, as well as polycyclic growth. This is the case for the teak tree. The aim of this work was to adapt the plant model, GreenLab, in order to take into consideration both these processes using existing data on this tree species. Methods: This work adopted GreenLab formalism based on source-sink relationships at organ level that drive biomass production and partitioning within the whole plant over time. The stochastic aspect of phytomer production can be modelled by a Bernoulli process. The teak model was designed, parameterized and analysed using the architectural data from 2- to 5-year-old teak trees in open field stands. Key results: Growth and development parameters were identified, fitting the observed compound organic series with the theoretical series, using generalized least squares methods. Phytomer distributions of growth units and branching pattern varied depending on their axis category, i.e. their physiological age. These emerging properties were in accordance with the observed growth patterns and biomass allocation dynamics during a growing season marked by a short dry season. Conclusions: Annual growth patterns observed on teak, including shoot pre- and neoformation and polycyclism, were reproduced by the new version of the GreenLab model. However, further updating is discussed in order to ensure better consideration of radial variation in basic specific gravity of wood. Such upgrading of the model will enable teak ideotypes to be defined for improving wood production in terms of both volume and quality.

Spatiotemporal organization of cilia drives multiscale mucus swirls in model human bronchial epithelium.

Mucociliary clearance is a biomechanical mechanism of airway protection. It consists of the active transport along the bronchial tree of the mucus, a fluid propelled by the coordinated beating of a myriad of cilia on the epithelial surface of the respiratory tract. The physics of mucus transport is poorly understood because it involves complex phenomena such as long-range hydrodynamic interactions, active collective ciliary motion, and the complex rheology of mucus. We propose a quantitative physical analysis of the ciliary activity and mucus transport on a large panel of human bronchial cultures from control subjects, patients with asthma and chronic obstructive pulmonary disease obtained from endobronchial biopsies. Here we report on the existence of multiple ciliary domains with sizes ranging from the tens of a micron to the centimeter, where ciliary beats present a circular orientational order. These domains are associated with circular mucus flow patterns, whose size scales with the average cilia density. In these domains, we find that the radial increase of the ciliated cell density coupled with the increase in the orientational order of ciliary beats result in a net local force proportional to the mucus velocity. We propose a phenomenological physical model that supports our results.

Interfacial rheological properties of self-assembling biopolymer microcapsules.

Tuning the mechanical properties of microcapsules through a cost-efficient route of fabrication is still a challenge. The traditional method of layer-by-layer assembly of microcapsules allows building a tailored composite multi-layer membrane but is technically complex as it requires numerous steps. The objective of this article is to characterize the interfacial rheological properties of self-assembling biopolymer microcapsules that were obtained in one single facile step. This thorough study provides new insights into the mechanics of these weakly cohesive membranes. Firstly, suspensions of water-in-oil microcapsules were formed in microfluidic junctions by self-assembly of two oppositely charged polyelectrolytes, namely chitosan (water soluble) and phosphatidic fatty acid (oil soluble). In this way, composite membranes of tunable thickness (between 40 and 900 nm measured by AFM) were formed at water/oil interfaces in a single step by changing the composition. Secondly, microcapsules were mechanically characterized by stretching them up to break-up in an extensional flow chamber which extends the relevance and convenience of the hydrodynamic method to weakly cohesive membranes. Finally, we show that the design of microcapsules can be 'engineered' in an extensive way since they present a wealth of interfacial rheological properties in terms of elasticity, plasticity and yield stress whose magnitudes can be controlled by the composition. These behaviors are explained by the variation of the membrane thickness with the physico-chemical parameters of the process.

How the Cobra Got Its Flesh-Eating Venom: Cytotoxicity as a Defensive Innovation and Its Co-Evolution with Hooding, Aposematic Marking, and Spitting.

The cytotoxicity of the venom of 25 species of Old World elapid snake was tested and compared with the morphological and behavioural adaptations of hooding and spitting. We determined that, contrary to previous assumptions, the venoms of spitting species are not consistently more cytotoxic than those of closely related non-spitting species. While this correlation between spitting and non-spitting was found among African cobras, it was not present among Asian cobras. On the other hand, a consistent positive correlation was observed between cytotoxicity and utilisation of the defensive hooding display that cobras are famous for. Hooding and spitting are widely regarded as defensive adaptations, but it has hitherto been uncertain whether cytotoxicity serves a defensive purpose or is somehow useful in prey subjugation. The results of this study suggest that cytotoxicity evolved primarily as a defensive innovation and that it has co-evolved twice alongside hooding behavior: once in the Hemachatus + Naja and again independently in the king cobras (Ophiophagus). There was a significant increase of cytotoxicity in the Asian Naja linked to the evolution of bold aposematic hood markings, reinforcing the link between hooding and the evolution of defensive cytotoxic venoms. In parallel, lineages with increased cytotoxicity but lacking bold hood patterns evolved aposematic markers in the form of high contrast body banding. The results also indicate that, secondary to the evolution of venom rich in cytotoxins, spitting has evolved three times independently: once within the African Naja, once within the Asian Naja, and once in the Hemachatus genus. The evolution of cytotoxic venom thus appears to facilitate the evolution of defensive spitting behaviour. In contrast, a secondary loss of cytotoxicity and reduction of the hood occurred in the water cobra Naja annulata, which possesses streamlined neurotoxic venom similar to that of other aquatic elapid snakes (e.g., hydrophiine sea snakes). The results of this study make an important contribution to our growing understanding of the selection pressures shaping the evolution of snake venom and its constituent toxins. The data also aid in elucidating the relationship between these selection pressures and the medical impact of human snakebite in the developing world, as cytotoxic cobras cause considerable morbidity including loss-of-function injuries that result in economic and social burdens in the tropics of Asia and sub-Saharan Africa.

Viper Venom Botox: The Molecular Origin and Evolution of the Waglerin Peptides Used in Anti-Wrinkle Skin Cream.

The molecular origin of waglerin peptides has remained enigmatic despite their industrial application in skin cream products to paralyse facial muscles and thus reduce the incidence of wrinkles. Here we show that these neurotoxic peptides are the result of de novo evolution within the prepro region of the C-type natriuretic peptide gene in Tropidolaemus venoms, at a site distinct from the domain encoding for the natriuretic peptide. It is the same region that yielded the azemiopsin peptides from Azemiops feae, indicative of a close relationship of this toxin gene between these two genera. The precursor region for the molecular evolution is a biodiversity hotspot that has yielded other novel bioactive peptides with novel activities. We detail the diversity of components in this and other species in order to explore what characteristics enable it to be such a biodiscovery treasure trove. The unusual function of Tropidolaemus venoms may have been selected for due to evolutionary pressures brought about by a high likelihood of prey escape.

Automatic identification and characterization of radial files in light microscopy images of wood.

BACKGROUND AND AIMS: Analysis of anatomical sections of wood provides important information for understanding the secondary growth and development of plants. This study reports on a new method for the automatic detection and characterization of cell files in wood images obtained by light microscopy. To facilitate interpretation of the results, reliability coefficients have been determined, which characterize the files, their cells and their respective measurements. METHODS: Histological sections and blocks of the gymnosperms Pinus canariensis, P. nigra and Abies alba were used, together with histological sections of the angiosperm mahogany (Swietenia spp.). Samples were scanned microscopically and mosaic images were built up. After initial processing to reduce noise and enhance contrast, cells were identified using a 'watershed' algorithm and then cell files were built up by the successive aggregation of cells taken from progressively enlarged neighbouring regions. Cell characteristics such as thickness and size were calculated, and a method was developed to determine the reliability of the measurements relative to manual methods. KEY RESULTS: Image analysis using this method can be performed in less than 20 s, which compares with a time of approx. 40 min to produce the same results manually. The results are accompanied by a reliability indicator that can highlight specific configurations of cells and also potentially erroneous data. CONCLUSIONS: The method provides a fast, economical and reliable tool for the identification of cell files. The reliability indicator characterizing the files permits quick filtering of data for statistical analysis while also highlighting particular biological configurations present in the wood sections.

Sedimentation-induced tether on a settling vesicle.

Destabilization of soft interfaces into thin cylindrical filaments under external stresses is ubiquitous and is generally the first step toward breakup. We show that such filaments, called tethers, emerge from a vesicle subjected to gravity. Contrary to the pendant drop experiment, we demonstrate that the bending rigidity, a specific membrane property of vesicles, ensures the tethers reach a stationary state. Moreover, unlike point-like force experiments, we show that the family of shapes is continuous.

Plant growth and architectural modelling and its applications. Preface.

Over the last decade, a growing number of scientists around the world have invested in research on plant growth and architectural modelling and applications (often abbreviated to plant modelling and applications, PMA). By combining physical and biological processes, spatially explicit models have shown their ability to help in understanding plant­environment interactions. This Special Issue on plant growth modelling presents new information within this topic, which are summarized in this preface. Research results for a variety of plant species growing in the field, in greenhouses and in natural environments are presented. Various models and simulation platforms are developed in this field of research, opening new features to a wider community of researchers and end users. New modelling technologies relating to the structure and function of plant shoots and root systems are explored from the cellular to the whole-plant and plant-community levels.

Three-dimensional distribution of vessels, passage cells and lateral roots along the root axis of winter wheat (Triticum aestivum).

BACKGROUND AND AIMS: The capacity of a plant to absorb and transport water and nutrients depends on anatomical structures within the roots and their co-ordination. However, most descriptions of root anatomical structure are limited to 2-D cross-sections, providing little information on 3-D spatial relationships and hardly anything on their temporal evolution. Three-dimensional reconstruction and visualization of root anatomical structures can illustrate spatial co-ordination among cells and tissues and provide new insights and understanding of the interrelation between structure and function. METHODS: Classical paraffin serial-section methods, image processing, computer-aided 3-D reconstruction and 3-D visualization techniques were combined to analyse spatial relationships among metaxylem vessels, passage cells and lateral roots in nodal roots of winter wheat (Triticum aestivum). KEY RESULTS: 3-D reconstruction demonstrated that metaxylem vessels were neither parallel, nor did they run directly along the root axis from the root base to the root tip; rather they underwent substitution and transition. Most vessels were connected to pre-existent or newly formed vessels by pits on their lateral walls. The spatial distributions of both passage cells and lateral roots exhibited similar position-dependent patterns. In the transverse plane, the passage cells occurred opposite the poles of the protoxylem and the lateral roots opposite those of the protophloem. Along the axis of a young root segment, the passage cells were arranged in short and discontinuous longitudinal files, thus as the tissues mature, the sequence in which the passage cells lose their transport function is not basipetal. In older segments, passage cells decreased drastically in number and coexisted with lateral roots. The spatial distribution of lateral roots was similar to that of the passage cells, mirroring their similar functions as lateral pathways for water and nutrient transport to the stele. CONCLUSIONS: With the 3-D reconstruction and visualization techniques developed here, the spatial relationships between vessels, passage cells and lateral roots and the temporal evolution of these relationships can be described. The technique helps to illustrate synchronization and spatial co-ordination among the root's radial and axial pathways for water and nutrient transport and the interdependence of structure and function in the root.

Dynamics of vesicle unbinding under axisymmetric flow.

The competition between adhesion and external flow to unbind settled vesicles from substrates is investigated. An experimental setup is developed to apply a hydrodynamic pulling force in the range of a few piconewtons to a vesicle with retained axisymmetry. In the limit of a small excess of membrane area, vesicles are found to transit during unbinding from a process of fluid film thickening at constant contact area to a finite-time process of contact radius drop to zero with an exponent 1/2. Both characteristic times vary linearly with the inverse flow rate. On the contrary, deflated vesicles under a moderate pulling force exhibit a decrease of contact area at a constant film thickness before a film thickening.

The first skull of the earliest giant panda.

Fossils of the giant panda Ailuropoda (Order Carnivora, Family Ursidae) are largely isolated teeth, mandibles, and a few rare skulls, known from the late Pliocene to late Pleistocene in China and Southeast Asia. Much of this material represents a Pleistocene chronospecies, Ailuropoda baconi, an animal larger than the living giant panda, Ailuropoda melanoleuca. The earliest certain record of Ailuropoda is the late Pliocene chronospecies, Ailuropoda microta, smaller than either A. baconi or A. melanoleuca, and previously known only from teeth and a few mandibles from karst caves in south China. Here, we report the discovery of the first skull of A. microta, establishing its cranial anatomy and demonstrating that the specialized cranial and dental adaptations of Ailuropoda for durophagous feeding behavior centered on bamboo were already evident in this late Pliocene species. The skull from Jinyin cave (Guangxi) and dental remains from other karst localities in southeastern China show that Ailuropoda microta occupied south China from approximately 2 to 2.4 Myr ago after a marked global climatic deterioration. Dental and basicranial anatomy indicate a less specialized morphology early in the history of the lineage and support derivation of the giant panda from the Miocene Asian ursid Ailurarctos.