The human EDAR 370V/A polymorphism affects tooth root morphology potentially through the modification of a reaction-diffusion system.

Morphological variations in human teeth have long been recognized and, in particular, the spatial and temporal distribution of two patterns of dental features in Asia, i.e., Sinodonty and Sundadonty, have contributed to our understanding of the human migration history. However, the molecular mechanisms underlying such dental variations have not yet been completely elucidated. Recent studies have clarified that a nonsynonymous variant in the ectodysplasin A receptor gene (EDAR 370V/A; rs3827760) contributes to crown traits related to Sinodonty. In this study, we examined the association between the EDAR polymorphism and tooth root traits by using computed tomography images and identified that the effects of the EDAR variant on the number and shape of roots differed depending on the tooth type. In addition, to better understand tooth root morphogenesis, a computational analysis for patterns of tooth roots was performed, assuming a reaction-diffusion system. The computational study suggested that the complicated effects of the EDAR polymorphism could be explained when it is considered that EDAR modifies the syntheses of multiple related molecules working in the reaction-diffusion dynamics. In this study, we shed light on the molecular mechanisms of tooth root morphogenesis, which are less understood in comparison to those of tooth crown morphogenesis.

Assessment of Polygala paniculata (Polygalaceae) characteristics for evolutionary studies of legume-rhizobia symbiosis.

Root nodule (RN) symbiosis is a mutualistic interaction observed between nitrogen-fixing soil bacteria and nodulating plants, which are scattered in only four orders of angiosperms called nitrogen-fixing clade. Most of legumes engage in RN symbiosis with rhizobia. Molecular genetic analyses with legumes and non-leguminous nodulating plants revealed that RN symbiosis utilizes early signalling components that are required for symbiosis with arbuscular mycorrhizal (AM) fungi. However detailed evolutionary processes are still largely unknown. Comparative analyses with non-nodulating species phylogenetically related to legumes could be better strategies to study the evolution of RN symbiosis in legumes. Polygala paniculata is a non-leguminous species that belongs to a family different from legumes but that is classified into the same order, Fabales. It has appropriate characteristics for cultivation in laboratories: small body size, high fertility and short lifecycles. Therefore, we further assessed whether this species is suitable as a model species for comparative studies with legumes. We first validated that the plant we obtained in Palau was truly P. paniculata by molecular phylogenetic analysis using rbcL sequences. The estimated genome size of this species was less than those of two model legumes, Lotus japonicus and Medicago truncatula. We determined conditions for cultivation in vitro and for hairy root formation from P. paniculata seedlings. It would facilitate to investigate gene functions in this species. The ability of P. paniculata to interact with AM fungi was confirmed by inoculation with Rhizophagus irregularis, suggesting the presence of early signalling factors that might be involved in RN symbiosis. Unexpectedly, branching of root hairs was observed when inoculated with Mesorhizobium loti broad host range strain NZP2037, indicating that P. paniculata has the biological potential to respond to rhizobia. We propose that P. paniculata is used as a model plant for the evolutionary study of RN symbiosis.

Spatial regulation of resource allocation in response to nutritional availability.

It is critical for a living organism to appropriately allocate resources among its organs, or within a specific organ, because available resources are generally limited. For example, in response to the nutritional environment of their soil, plants regulate resource allocation in their roots in order to plastically change their root system architecture (RSA) for efficiently absorbing nutrients. However, it is still not understood why and how RSA is adaptively controlled. Therefore, we modeled and investigated the spatial regulation of resource allocation, focusing on RSA in response to nutrient availability, and provided analytical solutions to the optimal strategy in the case of simple fitness functions. We first showed that our model could explain the experimental evidence where root growth is maximized at the optimal nutrient concentration under the homogeneous condition. Next, we extended our model to incorporate the spatial heterogeneity of nutrient availability. This extended model revealed that growth suppression by systemic control is required for adapting to high nutrient conditions, whereas growth promotion by local control is sufficient for adaptation to low-nutrient environments. This evidence predicts that systemic control can be evolved in the presence of excessive amounts of nutrition, consistent with the 'N-supply' systemic signal that is observed experimentally. Furthermore, our model can also explain various experimental results using nitrogen nutrition. Our model provides a theoretical basis for understanding the spatial regulation of adaptive resource allocation in response to nutritional environment.

Mathematical model studies of the comprehensive generation of major and minor phyllotactic patterns in plants with a predominant focus on orixate phyllotaxis.

Plant leaves are arranged around the stem in a beautiful geometry that is called phyllotaxis. In the majority of plants, phyllotaxis exhibits a distichous, Fibonacci spiral, decussate, or tricussate pattern. To explain the regularity and limited variety of phyllotactic patterns, many theoretical models have been proposed, mostly based on the notion that a repulsive interaction between leaf primordia determines the position of primordium initiation. Among them, particularly notable are the two models of Douady and Couder (alternate-specific form, DC1; more generalized form, DC2), the key assumptions of which are that each leaf primordium emits a constant power that inhibits new primordium formation and that this inhibitory effect decreases with distance. It was previously demonstrated by computer simulations that any major type of phyllotaxis can occur as a self-organizing stable pattern in the framework of DC models. However, several phyllotactic types remain unaddressed. An interesting example is orixate phyllotaxis, which has a tetrastichous alternate pattern with periodic repetition of a sequence of different divergence angles: 180°, 90°, -180°, and -90°. Although the term orixate phyllotaxis was derived from Orixa japonica, this type is observed in several distant taxa, suggesting that it may reflect some aspects of a common mechanism of phyllotactic patterning. Here we examined DC models regarding the ability to produce orixate phyllotaxis and found that model expansion via the introduction of primordial age-dependent changes of the inhibitory power is absolutely necessary for the establishment of orixate phyllotaxis. The orixate patterns generated by the expanded version of DC2 (EDC2) were shown to share morphological details with real orixate phyllotaxis. Furthermore, the simulation results obtained using EDC2 fitted better the natural distribution of phyllotactic patterns than did those obtained using the previous models. Our findings imply that changing the inhibitory power is generally an important component of the phyllotactic patterning mechanism.

Spatial regularity control of phyllotaxis pattern generated by the mutual interaction between auxin and PIN1.

Phyllotaxis, the arrangement of leaves on a plant stem, is well known because of its beautiful geometric configuration, which is derived from the constant spacing between leaf primordia. This phyllotaxis is established by mutual interaction between a diffusible plant hormone auxin and its efflux carrier PIN1, which cooperatively generate a regular pattern of auxin maxima, small regions with high auxin concentrations, leading to leaf primordia. However, the molecular mechanism of the regular pattern of auxin maxima is still largely unknown. To better understand how the phyllotaxis pattern is controlled, we investigated mathematical models based on the auxin-PIN1 interaction through linear stability analysis and numerical simulations, focusing on the spatial regularity control of auxin maxima. As in previous reports, we first confirmed that this spatial regularity can be reproduced by a highly simplified and abstract model. However, this model lacks the extracellular region and is not appropriate for considering the molecular mechanism. Thus, we investigated how auxin maxima patterns are affected under more realistic conditions. We found that the spatial regularity is eliminated by introducing the extracellular region, even in the presence of direct diffusion between cells or between extracellular spaces, and this strongly suggests the existence of an unknown molecular mechanism. To unravel this mechanism, we assumed a diffusible molecule to verify various feedback interactions with auxin-PIN1 dynamics. We revealed that regular patterns can be restored by a diffusible molecule that mediates the signaling from auxin to PIN1 polarization. Furthermore, as in the one-dimensional case, similar results are observed in the two-dimensional space. These results provide a great insight into the theoretical and molecular basis for understanding the phyllotaxis pattern. Our theoretical analysis strongly predicts a diffusible molecule that is pivotal for the phyllotaxis pattern but is yet to be determined experimentally.

Molecular Framework of a Regulatory Circuit Initiating Two-Dimensional Spatial Patterning of Stomatal Lineage.

Stomata, valves on the plant epidermis, are critical for plant growth and survival, and the presence of stomata impacts the global water and carbon cycle. Although transcription factors and cell-cell signaling components regulating stomatal development have been identified, it remains unclear as to how their regulatory interactions are translated into two-dimensional patterns of stomatal initial cells. Using molecular genetics, imaging, and mathematical simulation, we report a regulatory circuit that initiates the stomatal cell-lineage. The circuit includes a positive feedback loop constituting self-activation of SCREAMs that requires SPEECHLESS. This transcription factor module directly binds to the promoters and activates a secreted signal, EPIDERMAL PATTERNING FACTOR2, and the receptor modifier TOO MANY MOUTHS, while the receptor ERECTA lies outside of this module. This in turn inhibits SPCH, and hence SCRMs, thus constituting a negative feedback loop. Our mathematical model accurately predicts all known stomatal phenotypes with the inclusion of two additional components to the circuit: an EPF2-independent negative-feedback loop and a signal that lies outside of the SPCH•SCRM module. Our work reveals the intricate molecular framework governing self-organizing two-dimensional patterning in the plant epidermis.

Oriented cell division shapes carnivorous pitcher leaves of Sarracenia purpurea.

Complex morphology is an evolutionary outcome of phenotypic diversification. In some carnivorous plants, the ancestral planar leaf has been modified to form a pitcher shape. However, how leaf development was altered during evolution remains unknown. Here we show that the pitcher leaves of Sarracenia purpurea develop through cell division patterns of adaxial tissues that are distinct from those in bifacial and peltate leaves, subsequent to standard expression of adaxial and abaxial marker genes. Differences in the orientation of cell divisions in the adaxial domain cause bifacial growth in the distal region and adaxial ridge protrusion in the middle region. These different growth patterns establish pitcher morphology. A computer simulation suggests that the cell division plane is critical for the pitcher morphogenesis. Our results imply that tissue-specific changes in the orientation of cell division underlie the development of a morphologically complex leaf.

Missed massive morel-lavallee lesion.

A Morel-Lavallee lesion (MLL) involves posttraumatic fluid collection around the greater trochanter. Many cases of MLL are missed at the initial evaluation, and the treatment of MLL is not well established. We present two cases in which MLL was missed at the initial evaluation. Case 1. A 65-year-old man was run over by a parade float. There was subcutaneous hematoma around the left greater trochanter, and no fracture was found. We diagnosed this injury as MLL on the 7th day after the trauma. Although we performed percutaneous drainage, the injured area was infected. Case 2. A 57-year-old man was hit by a train in a factory. There was an iliac wing fracture, but an MLL was not initially recognized. On the 6th day after the trauma, when performing open reduction and internal fixation for the iliac fracture, we recognized the lesion and performed percutaneous drainage simultaneously. This lesion also became infected. In these two cases, the wounds finally healed after a long duration of treatment. We suggest that it is important to keep this injury in mind and debride the lesion early and completely in the treatment course.

Evolutionary dynamics of nitrogen fixation in the legume-rhizobia symbiosis.

The stabilization of host-symbiont mutualism against the emergence of parasitic individuals is pivotal to the evolution of cooperation. One of the most famous symbioses occurs between legumes and their colonizing rhizobia, in which rhizobia extract nutrients (or benefits) from legume plants while supplying them with nitrogen resources produced by nitrogen fixation (or costs). Natural environments, however, are widely populated by ineffective rhizobia that extract benefits without paying costs and thus proliferate more efficiently than nitrogen-fixing cooperators. How and why this mutualism becomes stabilized and evolutionarily persists has been extensively discussed. To better understand the evolutionary dynamics of this symbiosis system, we construct a simple model based on the continuous snowdrift game with multiple interacting players. We investigate the model using adaptive dynamics and numerical simulations. We find that symbiotic evolution depends on the cost-benefit balance, and that cheaters widely emerge when the cost and benefit are similar in strength. In this scenario, the persistence of the symbiotic system is compatible with the presence of cheaters. This result suggests that the symbiotic relationship is robust to the emergence of cheaters, and may explain the prevalence of cheating rhizobia in nature. In addition, various stabilizing mechanisms, such as partner fidelity feedback, partner choice, and host sanction, can reinforce the symbiotic relationship by affecting the fitness of symbionts in various ways. This result suggests that the symbiotic relationship is cooperatively stabilized by various mechanisms. In addition, mixed nodule populations are thought to encourage cheater emergence, but our model predicts that, in certain situations, cheaters can disappear from such populations. These findings provide a theoretical basis of the evolutionary dynamics of legume-rhizobia symbioses, which is extendable to other single-host, multiple-colonizer systems.

Pattern dynamics in adaxial-abaxial specific gene expression are modulated by a plastid retrograde signal during Arabidopsis thaliana leaf development.

The maintenance and reformation of gene expression domains are the basis for the morphogenic processes of multicellular systems. In a leaf primordium of Arabidopsis thaliana, the expression of FILAMENTOUS FLOWER (FIL) and the activity of the microRNA miR165/166 are specific to the abaxial side. This miR165/166 activity restricts the target gene expression to the adaxial side. The adaxial and abaxial specific gene expressions are crucial for the wide expansion of leaf lamina. The FIL-expression and the miR165/166-free domains are almost mutually exclusive, and they have been considered to be maintained during leaf development. However, we found here that the position of the boundary between the two domains gradually shifts from the adaxial side to the abaxial side. The cell lineage analysis revealed that this boundary shifting was associated with a sequential gene expression switch from the FIL-expressing (miR165/166 active) to the miR165/166-free (non-FIL-expressing) states. Our genetic analyses using the enlarged fil expression domain2 (enf2) mutant and chemical treatment experiments revealed that impairment in the plastid (chloroplast) gene expression machinery retards this boundary shifting and inhibits the lamina expansion. Furthermore, these developmental effects caused by the abnormal plastids were not observed in the genomes uncoupled1 (gun1) mutant background. This study characterizes the dynamic nature of the adaxial-abaxial specification process in leaf primordia and reveals that the dynamic process is affected by the GUN1-dependent retrograde signal in response to the failure of plastid gene expression. These findings advance our understanding on the molecular mechanism linking the plastid function to the leaf morphogenic processes.

Pattern formation by two-layer Turing system with complementary synthesis.

Many multicellular organisms have a layered structure. The interaction between these layers plays an essential role in many developmental processes, and key molecules involved in these processes are often expressed in a layer-specific manner. On the other hand, pattern formation of organisms has been frequently discussed in connection with the Turing system. However, the Turing system has so far been studied mainly in single-layered space. In this paper, we thus investigate a two-layer Turing system with complementary synthesis, in which two interacting molecules are exclusively synthesized in different layers. From a linear stability analysis, we determine the Turing condition of the complementary system, and show that this condition requires stronger regulatory interactions of the molecules than that of the system with usual ubiquitous synthesis. We then confirm that this complementary system affects pattern types in fixed and expanding two-dimensional spaces in a similar way to the system with ubiquitous synthesis. In addition, the two-layer system includes two types of diffusion, lateral and transversal, and these have distinct effects on pattern formation with lateral diffusion mainly determining the periodicity of patterns generated and transversal diffusion affecting pattern type. These results suggest that the transversal diffusion functions as a time delay in the two-layer system. Finally, we apply this complementary system to explain pattern formation of the shoot apical meristem of plants. These findings provide an understanding of pattern formation caused by the interaction between cell layers in multicellular organisms.

Strategy for shoot meristem proliferation in plants.

Shoot apical meristem (SAM) of plants harbors stem cells capable of generating the aerial tissues including reproductive organs. Therefore, it is very important for plants to control SAM proliferation and its density as a survival strategy. The SAM is regulated by the dynamics of a specific gene network, such as the WUS-CLV interaction of A. thaliana. By using a mathematical model, we previously proposed six possible SAM patterns in terms of the manner and frequency of stem cell proliferation. Two of these SAM patterns are predicted to generate either dichotomous or axillary shoot branch. Dichotomous shoot branches caused by this mechanism are characteristic of the earliest vascular plants, such as Cooksonia and Rhynia, but are observed in only a small minority of plant species of the present day. On the other hand, axillary branches are observed in the majority of plant species and are induced by a different dynamics of the feedback regulation between auxin and the asymmetric distribution of PIN auxin efflux carriers. During evolution, some plants may have adopted this auxin-PIN system to more strictly control SAM proliferation.

Reaction-diffusion pattern in shoot apical meristem of plants.

A fundamental question in developmental biology is how spatial patterns are self-organized from homogeneous structures. In 1952, Turing proposed the reaction-diffusion model in order to explain this issue. Experimental evidence of reaction-diffusion patterns in living organisms was first provided by the pigmentation pattern on the skin of fishes in 1995. However, whether or not this mechanism plays an essential role in developmental events of living organisms remains elusive. Here we show that a reaction-diffusion model can successfully explain the shoot apical meristem (SAM) development of plants. SAM of plants resides in the top of each shoot and consists of a central zone (CZ) and a surrounding peripheral zone (PZ). SAM contains stem cells and continuously produces new organs throughout the lifespan. Molecular genetic studies using Arabidopsis thaliana revealed that the formation and maintenance of the SAM are essentially regulated by the feedback interaction between WUSHCEL (WUS) and CLAVATA (CLV). We developed a mathematical model of the SAM based on a reaction-diffusion dynamics of the WUS-CLV interaction, incorporating cell division and the spatial restriction of the dynamics. Our model explains the various SAM patterns observed in plants, for example, homeostatic control of SAM size in the wild type, enlarged or fasciated SAM in clv mutants, and initiation of ectopic secondary meristems from an initial flattened SAM in wus mutant. In addition, the model is supported by comparing its prediction with the expression pattern of WUS in the wus mutant. Furthermore, the model can account for many experimental results including reorganization processes caused by the CZ ablation and by incision through the meristem center. We thus conclude that the reaction-diffusion dynamics is probably indispensable for the SAM development of plants.

Arabidopsis PIS1 encodes the ABCG37 transporter of auxinic compounds including the auxin precursor indole-3-butyric acid.

Differential distribution of the plant hormone auxin within tissues mediates a variety of developmental processes. Cellular auxin levels are determined by metabolic processes including synthesis, degradation, and (de)conjugation, as well as by auxin transport across the plasma membrane. Whereas transport of free auxins such as naturally occurring indole-3-acetic acid (IAA) is well characterized, little is known about the transport of auxin precursors and metabolites. Here, we identify a mutation in the ABCG37 gene of Arabidopsis that causes the polar auxin transport inhibitor sensitive1 (pis1) phenotype manifested by hypersensitivity to auxinic compounds. ABCG37 encodes the pleiotropic drug resistance transporter that transports a range of synthetic auxinic compounds as well as the endogenous auxin precursor indole-3-butyric acid (IBA), but not free IAA. ABCG37 and its homolog ABCG36 act redundantly at outermost root plasma membranes and, unlike established IAA transporters from the PIN and ABCB families, transport IBA out of the cells. Our findings explore possible novel modes of regulating auxin homeostasis and plant development by means of directional transport of the auxin precursor IBA and presumably also other auxin metabolites.

Morphological effects of sinefungin, an inhibitor of S-adenosylmethionine-dependent methyltransferases, on Anabaena sp. PCC 7120.

Anabaena cells develop regular one-dimensional filaments through cell division in planes parallel to each other. A gcvP mutant displayed morphological defects such as filaments with sharp bends and/or branching and irregular cell clumps. The defects probably result from depletion of S-adenosylmethionine (AdoMet), because they were rescued by the application of methionine, an AdoMet precursor, and because sinefungin, a strong inhibitor of AdoMet-dependent methyltransferases, caused morphological abnormalities in wild-type Anabaena similar to those of the mutant. AdoMet-dependent methylation is involved in the spatial regulation of cell polarity in Anabaena.

The origin of the diversity of leaf venation pattern.

The leaf venation pattern of plants shows remarkable diversity and species-specificity. However, the mechanism underlying the pattern formation and pattern diversity remains unclear. We developed a mathematical model that is based on the positive feedback regulation between plant hormone auxin and its efflux carrier. This system can generate auxin flow pathways by self-organization from an almost homogeneous state. This result explains a well-known experimental phenomenon referred as to "polar auxin transport." The model can produce diverse leaf venation patterns with spatial regularity under similar conditions to those of leaf development, that is, in the presence of leaf expansion and auxin sink. Final venation patterns are strikingly affected by leaf shape and leaf expansion. These results indicate that the positive feedback regulation between auxin and its efflux carrier is a central dynamic in leaf venation pattern formation. The diversity of leaf venation patterns in plant species is probably due to the differences of leaf shape and leaf expansion pattern.

Pattern formation of leaf veins by the positive feedback regulation between auxin flow and auxin efflux carrier.

The generation of vascular pattern formation in plants is an interesting process of pattern formation in organisms. It is well known that the plant hormone auxin is involved in plant vascular differentiation and that the PIN1 protein, an auxin efflux carrier, localizes to one side of the cell membrane. Several hypotheses have been proposed to explain the formation of leaf venation. One is the canalization hypothesis that is based on the assumption that a positive feedback regulation exists between the flow of a signal molecule and the capacity of its flow. Here, we attempted to integrate the canalization hypothesis and experimental data. We investigated models of the positive feedback regulation between the auxin flow and PIN1 localization. Model 1, with conserved PIN1 amount in each cell, can generate a branching pattern similar to that of plant leaf venation. We introduced the diffusible enhancer "e" into the model as unknown factor. The obtained patterns show a quasi-periodic distribution of auxin flow paths, when the model dynamics includes domain growth. In order to understand the early initiation process that generates an inhomogeneity from an almost homogeneous distribution, we introduced model 2, a simplified version of model 1. Model 2 can generate inhomogeneity with a parameter dependency similar to that of model 1. To analyse parameter condition required for pattern development, approximated equations are obtained from model 2. The isocline analysis of the equations without spatial structure shows that the inhomogeneous distribution occurs from an almost homogeneous distribution. This parameter condition for generating inhomogeneity is consistent with the results of models 1 and 2.

A Lotus basic leucine zipper protein with a RING-finger motif negatively regulates the developmental program of nodulation.

The developmental program of nodulation is regulated systemically in leguminous host species. A mutant astray (Ljsym77) in Lotus japonicus has lost some sort of its ability to regulate this symtem, and shows enhanced and early nodulation. In the absence of rhizobia, this mutant exhibits characteristics associated with defects in light and gravity responses. These nonsymbiotic phenotypes of astray are very similar to those observed in photomorphogenic Arabidopsis mutant hy5. Based on this evidence, we predicted that astray might contain a mutation in the HY5 homologue of L. japonicus. The homologue, named LjBzf, encodes a basic leucine zipper protein in the C-terminal half that shows the highest level of identity with HY5 of all Arabidopsis proteins. It also encodes legume-characteristic combination of motifs, including a RING-finger motif and an acidic region in the N-terminal half. The astray phenotypes were cosegregated with LjBzf, and the failure to splice the intron was detected. Nonsymbiotic and symbiotic phenotypes of astray were complemented by introduction of CaMV35SLjBzf. It is noteworthy that although Arabidopsis hy5 showed an enhancement of lateral root initiation, Lotus astray showed an enhancement of nodule initiation but not of lateral root initiation. Legume-characteristic combination of motifs of ASTRAY may play specific roles in the regulation of nodule development.