Synthesis of mathematical models of branching axons and dendrites

A synthesis was made of models of branching neuronal cable structures from a full set of standard basic models. The study aimed to produce an instrument of mathematical modelling making it possible to reflect true life morphological and electrophysiological characteristics of axons and dendrites, discarding some of the restrictions and simplifications characterizing existing models of the structures mentioned. Equivalent electrical circuits of branching axons and dendrites were set up with in-series and node connections of standard four-terminal networks corresponding to basic segments with active or passive membrane. Equations were obtained for electrical processes in branching neuronal neurites, generalized in the case of multiple binary branching with arbitrary symmetry and branching structure. A difference scheme common to the whole class of models contemplated was produced and the algorithm of a numerical solution to the difference equations thus obtained was elaborated. The instrument described makes it possible to synthesize diverse models of branching axons and dendrites, offering considerably greater opportunities for modelling the main electrophysiological processes developing in these structures of electrotonus, propagation of excitation, and interaction between these two factors.State University Commemorating Tricentenary of Russo-Ukrainian Union. Dnepropetrovsk. Translated from Neirofiziologiya, Vol. 20, No. 4, pp. 471–479, July–August, 1988.     

Predicting root biomass from branching patterns of Douglas-fir root systems

There are many examples of branching networks in nature, such as tree crowns, river systems, arteries and lungs. These networks have often been described as being self-similar, or following scale-invariant branching rules, and this property has been used to derive several scaling laws. In this paper we model root systems of Douglas-fir ( Pseudotsuga menziesii var. glauca (Beissn.) Franco) as branching networks following several simple branching rules. Our objective is to establish a relationship between trunk diameter and root biomass. We explore the effect of the self-similar branching assumption on this relationship. Using data collected from a mature stand in British Columbia, we find that branching asymmetry and the rate of root taper change with root size, thereby violating the assumption of self-similarity. However, the data are in general agreement with Leonardo da Vinci's area-preserving branching hypothesis. We use the field data to parameterize two models, one assuming self-similar branching and a second incorporating the measured size dependencies of branching parameters. The two models differ by only a small amount (≈8%) in their predictions. For both models, the predicted relationship between trunk diameter and root biomass is in good concordance with previously published empirical data. We conclude that the assumption of self-similar branching, although violated by the data, nevertheless provides a useful tool for predicting the allometric relationship between trunk diameter and root biomass. Finally, we use our models to show that the geometric properties of individual bifurcations fundamentally change the root biomass cost of different root topologies.     

Growth of branched actin networks against obstacles

A method for simulating the growth of branched actin networks against obstacles has been developed. The method is based on simple stochastic events, including addition or removal of monomers at filament ends, capping of filament ends, nucleation of branches from existing filaments, and detachment of branches; the network structure for several different models of the branching process has also been studied. The models differ with regard to their inclusion of effects such as preferred branch orientations, filament uncapping at the obstacle, and preferential branching at filament ends. The actin ultrastructure near the membrane in lamellipodia is reasonably well produced if preferential branching in the direction of the obstacle or barbed-end uncapping effects are included. Uncapping effects cause the structures to have a few very long filaments that are similar to those seen in pathogen-induced "actin tails." The dependence of the growth velocity, branch spacing, and network density on the rate parameters for the various processes is quite different among the branching models. An analytic theory of the growth velocity and branch spacing of the network is described. Experiments are suggested that could distinguish among some of the branching models.     

桃树修剪分枝模式的模拟

在不同修剪手法下,对栽培桃树(Prunuspersica(L.)Batsch)不同母枝上的分枝模式进行了比较研究.从分枝模式来看:修剪后的母枝基本由3个不同的区域组成,基部是不萌发的潜伏芽形成的未分枝区域;中部是延迟分枝和多次分枝组成的分枝区域(主要的枝条类型有短枝、长枝和多次枝);顶部是被剪除的部分.我们通过隐式半马尔可夫模型来模拟这一分枝模式,主要是定量描述1次枝和多次枝在母枝上的数量及其分布状况.在上述模型中,未分枝区、延迟分枝区和多次分枝区称为瞬时态,被剪除的部分称为吸收态.模拟的结果与观察的结果进行对比后发现,两者具有很好的一致性.这说明隐式半马尔可夫模型是模拟植物分枝过程的一种有效方法,尽管隐式半马尔可夫链模型只是一个描述性的模型,但仍能对其所描述的生物现象进行解释,在预测修剪手法对母枝分枝模式影响方面比传统的方法具有明显的优势.本研究结果是建立三维虚拟桃树树冠分枝结构的基础.     

Tree Branching: Leonardo da Vinci's Rule versus Biomechanical Models

This study examined Leonardo da Vinci''s rule (i.e., the sum of the cross-sectional area of all tree branches above a branching point at any height is equal to the cross-sectional area of the trunk or the branch immediately below the branching point) using simulations based on two biomechanical models: the uniform stress and elastic similarity models. Model calculations of the daughter/mother ratio (i.e., the ratio of the total cross-sectional area of the daughter branches to the cross-sectional area of the mother branch at the branching point) showed that both biomechanical models agreed with da Vinci''s rule when the branching angles of daughter branches and the weights of lateral daughter branches were small; however, the models deviated from da Vinci''s rule as the weights and/or the branching angles of lateral daughter branches increased. The calculated values of the two models were largely similar but differed in some ways. Field measurements of Fagus crenata and Abies homolepis also fit this trend, wherein models deviated from da Vinci''s rule with increasing relative weights of lateral daughter branches. However, this deviation was small for a branching pattern in nature, where empirical measurements were taken under realistic measurement conditions; thus, da Vinci''s rule did not critically contradict the biomechanical models in the case of real branching patterns, though the model calculations described the contradiction between da Vinci''s rule and the biomechanical models. The field data for Fagus crenata fit the uniform stress model best, indicating that stress uniformity is the key constraint of branch morphology in Fagus crenata rather than elastic similarity or da Vinci''s rule. On the other hand, mechanical constraints are not necessarily significant in the morphology of Abies homolepis branches, depending on the number of daughter branches. Rather, these branches were often in agreement with da Vinci''s rule.     

Fractal modeling of pulmonary blood flow heterogeneity

The heterogeneity of pulmonary blood flow is not adequately described by gravitational forces alone. We investigated the flow distributions predicted by two fractally branching vascular models to determine how well such networks could explain the observed heterogeneity. The distribution of flow was modeled with a dichotomously branching tree in which the fraction of blood flow from the parent to the daughter branches was gamma and 1-gamma repeatedly at each generation. In one model gamma was held constant throughout the network, and in the other model gamma varied about a mean of 0.5 with a standard deviation of sigma. Both gamma and sigma were optimized in each model for the best fit to pulmonary blood flow data from experimental animals. The predicted relative dispersion of flow from the two model fractal networks produced an excellent fit to the observed data. These fractally branching models relate structure and function of the pulmonary vascular tree and provide a mechanism to describe the spatially correlated distribution of flow and the gravity-independent heterogeneity of blood flow.     

桃树修剪分枝模式的模拟

在不同修剪手法下,对栽培桃树(Prunuspersica(L.)Batsch)不同母枝上的分枝模式进行了比较研究。从分枝模式来看修剪后的母枝基本由3个不同的区域组成,基部是不萌发的潜伏芽形成的未分枝区域;中部是延迟分枝和多次分枝组成的分枝区域(主要的枝条类型有短枝、长枝和多次枝);顶部是被剪除的部分。我们通过隐式半马尔可夫模型来模拟这一分枝模式,主要是定量描述1次枝和多次枝在母枝上的数量及其分布状况。在上述模型中,未分枝区、延迟分枝区和多次分枝区称为瞬时态,被剪除的部分称为吸收态。模拟的结果与观察的结果进行对比后发现,两者具有很好的一致性。这说明隐式半马尔可夫模型是模拟植物分枝过程的一种有效方法,尽管隐式半马尔可夫链模型只是一个描述性的模型,但仍能对其所描述的生物现象进行解释,在预测修剪手法对母枝分枝模式影响方面比传统的方法具有明显的优势。本研究结果是建立三维虚拟桃树树冠分枝结构的基础。     

A flow cell reactor for the study of growth kinetics of single hairy roots

Summary A flow cell reactor designed for the study of growth and branching of individual roots under defined fluid flow conditions is described and preliminary results with the hairy root clone Tagetes erecta T3 are presented. Observations of root growth in this reactor may help to formulate and verify mathematical models root growth kinetics.     

Relevance of the lariat model to the splicing of the H-2 gene family

We examined the splicing of the H-2 gene family, taking the H-2Kd as a prototype, in the framework of the lariat model. We mainly investigated the mechanism described by Konarska et al. [Nature (Lond.) 313, 552-557 (1985)] who propose a direct interaction between the 5' splicing site and the branching region. We also checked each of the H-2 introns for the presence of patterns resembling the published consensus for the branching region. The known splicing events in the H-2 gene family are not always consistent with the current models, and our results indicate that slightly different mechanisms govern the splicing of different introns. A tentative explanation of the alternative splicing of the first and last intron, previously described, is given. The removal of the large third intron is not easily rationalized unless new rules for an additional multistep processing are postulated.     

Growth velocities of branched actin networks

The growth of an actin network against an obstacle that stimulates branching locally is studied using several variants of a kinetic rate model based on the orientation-dependent number density of filaments. The model emphasizes the effects of branching and capping on the density of free filament ends. The variants differ in their treatment of side versus end branching and dimensionality, and assume that new branches are generated by existing branches (autocatalytic behavior) or independently of existing branches (nucleation behavior). In autocatalytic models, the network growth velocity is rigorously independent of the opposing force exerted by the obstacle, and the network density is proportional to the force. The dependence of the growth velocity on the branching and capping rates is evaluated by a numerical solution of the rate equations. In side-branching models, the growth velocity drops gradually to zero with decreasing branching rate, while in end-branching models the drop is abrupt. As the capping rate goes to zero, it is found that the behavior of the velocity is sensitive to the thickness of the branching region. Experiments are proposed for using these results to shed light on the nature of the branching process.     

Modeling of branching structures of plants

Previous studies of branching structures generally focused on arteries. Four cost models minimizing total surface area, total volume, total drag and total power losses at a junction point have been proposed to study branching structures. In this paper, we highlight the branching structures of plants and examine which model fits data of branching structures of plants the best. Though the effect of light (e.g. phototropism) and other possible factors are not included in these cost models, a simple cost model with physiological significance, needs to be verified before further research on modeling of branching structures is conducted. Therefore, data are analysed in this paper to determine the best cost model. Branching structures of plants are studied by measuring branching angles and diameters of 234 junctions from four species of plants. The sample includes small junctions, large junctions, two- and three-dimensional junctions, junctions with three branches joining at a point and those with four branches joining at a point. First, junction exponents (x) were determined. Second, log-log plots indicate that model of volume minimization fits data better than other models. Third, one-sided t -tests were used to compare the fitness of four models. It is found that model of volume minimization fits data better than other cost models.     

Distribution of end-points of a branching network with decaying branch length

The geometry of the human bronchial tree has been described as a network formed by successive dichotomous branching with constant branching angles and geometrically decaying branch lengths. Models having these properties and with randomly distributed branching planes are constructed. The distribution of the end points of the model networks are described by computing the variance of the distributions in the direction of the axis of the network and in the transverse directions. It is found that, for a given decay ratio, there is a branching angle for which the volume filled by the end points is a maximum. The advantages of the network with the decay ratio and branching angle of the human bronchial tree are discussed.     

Reaction-diffusion models of growing plant tips: bifurcations on hemispheres

We study two chemical models for pattern formation in growing plant tips. For hemisphere radius and parameter values together optimal for spherical surface harmonic patterns of index l = 3, the Brusselator model gives an 84% probability of dichotomous branching pattern and 16% of annular pattern, while the hyperchirality model gives 88% probability of dichotomous branching and 12% of annular pattern. The models are two-morphogen reaction-diffusion systems on the surface of a hemispherical shell, with Dirichlet boundary conditions. Bifurcation analysis shows that both models give possible mechanisms for dichotomous branching of the growing tips. Symmetries of the models are used in the analysis.     

The exact probabilities of branching patterns under terminal and segmental growth hypotheses

Information about the way of branching of dendritic arborizations may be obtained by comparing the frequency distributions of observed branching patterns with theoretical distributions based on well-defined growth models. Two models usually get much attention in geomorphological and (neuro)biological studies, viz. terminal growth and segmental growth. Formulae to construct the exact probability distributions for both growth models are presented. It is shown that ranking and lumping of the individual branching patterns enable the analysis of very large arborizations with relatively few data. The application of the Kolmogorov goodness-of-fit test for discrete distributions to the analysis is discussed.     

Temperature jump study of charge translocation during the bacteriorhodopsin photocycle

Temperature jump experiments were carried out on purple membranes oriented and fixed in polyacrylamide gel. With green background illumination a relaxation of the photocurrent after an infrared laser pulse could be observed. To simulate the temperature jump signals different models of the bacteriorhodopsin photocycle were tested. The parameters of these models were obtained by measuring absorbance changes and photocurrent after excitation with a 575-nm laser flash.

A model with a temperature-dependent branching before the M state turned out to be satisfying. Other models, especially those with a late branching or without branching, could not reproduce the temperature jump measurements.

     

The Probability of Extinction of Infectious Salmon Anemia Virus in One and Two Patches

Single-type and multitype branching processes have been used to study the dynamics of a variety of stochastic birth–death type phenomena in biology and physics. Their use in epidemiology goes back to Whittle’s study of a susceptible–infected–recovered (SIR) model in the 1950s. In the case of an SIR model, the presence of only one infectious class allows for the use of single-type branching processes. Multitype branching processes allow for multiple infectious classes and have latterly been used to study metapopulation models of disease. In this article, we develop a continuous time Markov chain (CTMC) model of infectious salmon anemia virus in two patches, two CTMC models in one patch and companion multitype branching process (MTBP) models. The CTMC models are related to deterministic models which inform the choice of parameters. The probability of extinction is computed for the CTMC via numerical methods and approximated by the MTBP in the supercritical regime. The stochastic models are treated as toy models, and the parameter choices are made to highlight regions of the parameter space where CTMC and MTBP agree or disagree, without regard to biological significance. Partial extinction events are defined and their relevance discussed. A case is made for calculating the probability of such events, noting that MTBPs are not suitable for making these calculations.     

A simplified method of casting the macroscopic airways of lungs

A technique for preparing casts of the macroscopic airways of mammalian lungs, which is both simplified and inexpensive in comparison with previous techniques, is described. The models are accurate, durable and flexible, and clearly demonstrate the orientation and branching pattern of the bronchial tree. The nature of the procedure also extends the availability of casts to laboratories or individuals with limited instrumentation and/or funding. Preliminary results using this technique to inject the lungs and certain air sacs of birds are also discussed.     

Adaptive dynamics in diploid, sexual populations and the evolution of reproductive isolation

Evolutionary branching is the process whereby an initially monomorphic population evolves to a point where it undergoes disruptive selection and splits up into two phenotypically diverging lineages. We studied evolutionary branching in three models that are ecologically identical but that have different genetic systems. The first model is clonal, the second is sexual diploid with additive genetics on a single locus and the third is like the second but with an additional locus for mate choice. Evolutionary branching occurred under exactly the same ecological circumstances in all three models. After branching the evolutionary dynamics may be qualitatively different. In particular, in the diploid, sexual models there can be multiple evolutionary outcomes whereas in the corresponding clonal model there is only one. We showed that evolutionary branching favours the evolution of (partial) assortative mating and that this in turn effectively restores the results from the clonal model by rendering the alternative outcomes unreachable except for the one that also occurs in the clonal model. The evolution of assortative mating during evolutionary branching can be interpreted as the initial phase of sympatric speciation with phenotypic divergence and partial reproductive isolation.     

羊草种群无性系生长格局的研究

 本文根据野外实测数据和采用数学分析方法,在研究羊草种群无性系营养扩散和生长过程的基础上,建立了符合描述羊草种群无性系生长格局的相关随机走动模型和简单扩散模型。对模型进行检验,结果表明,所建模型可以定量地描述羊草种群无性系的生长格局。