Control of osteoblast regeneration by a train of Erk activity waves.

Regeneration is a complex chain of events that restores a tissue to its original size and shape. The tissue-wide coordination of cellular dynamics that is needed for proper morphogenesis is challenged by the large dimensions of regenerating body parts. Feedback mechanisms in biochemical pathways can provide effective communication across great distances1-5, but how they might regulate growth during tissue regeneration is unresolved6,7. Here we report that rhythmic travelling waves of Erk activity control the growth of bone in time and space in regenerating zebrafish scales, millimetre-sized discs of protective body armour. We find that waves of Erk activity travel across the osteoblast population as expanding concentric rings that are broadcast from a central source, inducing ring-like patterns of tissue growth. Using a combination of theoretical and experimental analyses, we show that Erk activity propagates as excitable trigger waves that are able to traverse the entire scale in approximately two days and that the frequency of wave generation controls the rate of scale regeneration. Furthermore, the periodic induction of synchronous, tissue-wide activation of Erk in place of travelling waves impairs tissue growth, which indicates that wave-distributed Erk activation is key to regeneration. Our findings reveal trigger waves as a regulatory strategy to coordinate cell behaviour and instruct tissue form during regeneration.

Modeling Edar expression reveals the hidden dynamics of tooth signaling center patterning.

When patterns are set during embryogenesis, it is expected that they are straightly established rather than subsequently modified. The patterning of the three mouse molars is, however, far from straight, likely as a result of mouse evolutionary history. The first-formed tooth signaling centers, called MS and R2, disappear before driving tooth formation and are thought to be vestiges of the premolars found in mouse ancestors. Moreover, the mature signaling center of the first molar (M1) is formed from the fusion of two signaling centers (R2 and early M1). Here, we report that broad activation of Edar expression precedes its spatial restriction to tooth signaling centers. This reveals a hidden two-step patterning process for tooth signaling centers, which was modeled with a single activator-inhibitor pair subject to reaction-diffusion (RD). The study of Edar expression also unveiled successive phases of signaling center formation, erasing, recovering, and fusion. Our model, in which R2 signaling center is not intrinsically defective but erased by the broad activation preceding M1 signaling center formation, predicted the surprising rescue of R2 in Edar mutant mice, where activation is reduced. The importance of this R2-M1 interaction was confirmed by ex vivo cultures showing that R2 is capable of forming a tooth. Finally, by introducing chemotaxis as a secondary process to RD, we recapitulated in silico different conditions in which R2 and M1 centers fuse or not. In conclusion, pattern formation in the mouse molar field relies on basic mechanisms whose dynamics produce embryonic patterns that are plastic objects rather than fixed end points.

Mathematical modeling of Erk activity waves in regenerating zebrafish scales.

Erk signaling regulates cellular decisions in many biological contexts. Recently, we have reported a series of Erk activity traveling waves that coordinate regeneration of osteoblast tissue in zebrafish scales. These waves originate from a central source region, propagate as expanding rings, and impart cell growth, thus controlling tissue morphogenesis. Here, we present a minimal reaction-diffusion model for Erk activity waves. The model considers three components: Erk, a diffusible Erk activator, and an Erk inhibitor. Erk stimulates both its activator and inhibitor, forming a positive and negative feedback loop, respectively. Our model shows that this system can be excitable and propagate Erk activity waves. Waves originate from a pulsatile source that is modeled by adding a localized basal production of the activator, which turns the source region from an excitable to an oscillatory state. As Erk activity periodically rises in the source, it can trigger an excitable wave that travels across the entire tissue. Analysis of the model finds that positive feedback controls the properties of the traveling wavefront and that negative feedback controls the duration of Erk activity peak and the period of Erk activity waves. The geometrical properties of the waves facilitate constraints on the effective diffusivity of the activator, indicating that waves are an efficient mechanism to transfer growth factor signaling rapidly across a large tissue.

Simultaneous CRISPR/Cas9-mediated editing of cassava eIF4E isoforms nCBP-1 and nCBP-2 reduces cassava brown streak disease symptom severity and incidence.

Cassava brown streak disease (CBSD) is a major constraint on cassava yields in East and Central Africa and threatens production in West Africa. CBSD is caused by two species of positive-sense RNA viruses belonging to the family Potyviridae, genus Ipomovirus: Cassava brown streak virus (CBSV) and Ugandan cassava brown streak virus (UCBSV). Diseases caused by the family Potyviridae require the interaction of viral genome-linked protein (VPg) and host eukaryotic translation initiation factor 4E (eIF4E) isoforms. Cassava encodes five eIF4E proteins: eIF4E, eIF(iso)4E-1, eIF(iso)4E-2, novel cap-binding protein-1 (nCBP-1), and nCBP-2. Protein-protein interaction experiments consistently found that VPg proteins associate with cassava nCBPs. CRISPR/Cas9-mediated genome editing was employed to generate ncbp-1, ncbp-2, and ncbp-1/ncbp-2 mutants in cassava cultivar 60444. Challenge with CBSV showed that ncbp-1/ncbp-2 mutants displayed delayed and attenuated CBSD aerial symptoms, as well as reduced severity and incidence of storage root necrosis. Suppressed disease symptoms were correlated with reduced virus titre in storage roots relative to wild-type controls. Our results demonstrate the ability to modify multiple genes simultaneously in cassava to achieve tolerance to CBSD. Future studies will investigate the contribution of remaining eIF4E isoforms on CBSD and translate this knowledge into an optimized strategy for protecting cassava from disease.

The first myriapod genome sequence reveals conservative arthropod gene content and genome organisation in the centipede Strigamia maritima.

Myriapods (e.g., centipedes and millipedes) display a simple homonomous body plan relative to other arthropods. All members of the class are terrestrial, but they attained terrestriality independently of insects. Myriapoda is the only arthropod class not represented by a sequenced genome. We present an analysis of the genome of the centipede Strigamia maritima. It retains a compact genome that has undergone less gene loss and shuffling than previously sequenced arthropods, and many orthologues of genes conserved from the bilaterian ancestor that have been lost in insects. Our analysis locates many genes in conserved macro-synteny contexts, and many small-scale examples of gene clustering. We describe several examples where S. maritima shows different solutions from insects to similar problems. The insect olfactory receptor gene family is absent from S. maritima, and olfaction in air is likely effected by expansion of other receptor gene families. For some genes S. maritima has evolved paralogues to generate coding sequence diversity, where insects use alternate splicing. This is most striking for the Dscam gene, which in Drosophila generates more than 100,000 alternate splice forms, but in S. maritima is encoded by over 100 paralogues. We see an intriguing linkage between the absence of any known photosensory proteins in a blind organism and the additional absence of canonical circadian clock genes. The phylogenetic position of myriapods allows us to identify where in arthropod phylogeny several particular molecular mechanisms and traits emerged. For example, we conclude that juvenile hormone signalling evolved with the emergence of the exoskeleton in the arthropods and that RR-1 containing cuticle proteins evolved in the lineage leading to Mandibulata. We also identify when various gene expansions and losses occurred. The genome of S. maritima offers us a unique glimpse into the ancestral arthropod genome, while also displaying many adaptations to its specific life history.

Noise Expands the Response Range of the Bacillus subtilis Competence Circuit.

Gene regulatory circuits must contend with intrinsic noise that arises due to finite numbers of proteins. While some circuits act to reduce this noise, others appear to exploit it. A striking example is the competence circuit in Bacillus subtilis, which exhibits much larger noise in the duration of its competence events than a synthetically constructed analog that performs the same function. Here, using stochastic modeling and fluorescence microscopy, we show that this larger noise allows cells to exit terminal phenotypic states, which expands the range of stress levels to which cells are responsive and leads to phenotypic heterogeneity at the population level. This is an important example of how noise confers a functional benefit in a genetic decision-making circuit.

Functional analysis of centipede development supports roles for Wnt genes in posterior development and segment generation.

The genes of the Wnt family play important and highly conserved roles in posterior growth and development in a wide range of animal taxa. Wnt genes also operate in arthropod segmentation, and there has been much recent debate regarding the relationship between arthropod and vertebrate segmentation mechanisms. Due to its phylogenetic position, body form, and possession of many (11) Wnt genes, the centipede Strigamia maritima is a useful system with which to examine these issues. This study takes a functional approach based on treatment with lithium chloride, which causes ubiquitous activation of canonical Wnt signalling. This is the first functional developmental study performed in any of the 15,000 species of the arthropod subphylum Myriapoda. The expression of all 11 Wnt genes in Strigamia was analyzed in relation to posterior development. Three of these genes, Wnt11, Wnt5, and WntA, were strongly expressed in the posterior region and, thus, may play important roles in posterior developmental processes. In support of this hypothesis, LiCl treatment of S. maritima embryos was observed to produce posterior developmental defects and perturbations in AbdB and Delta expression. The effects of LiCl differ depending on the developmental stage treated, with more severe effects elicited by treatment during germband formation than by treatment at later stages. These results support a role for Wnt signalling in conferring posterior identity in Strigamia. In addition, data from this study are consistent with the hypothesis of segmentation based on a "clock and wavefront" mechanism operating in this species.

The centipede Strigamia maritima possesses a large complement of Wnt genes with diverse expression patterns.

The genes of the Wnt family play important roles in the development of many animals. In the arthropods, these genes are known to have multiple functions, including roles in posterior development and segmentation. Despite this, secondary loss of Wnt genes is common among the Arthropoda. Unlike many arthropods, Strigamia maritima, a geophilomorph centipede, possesses a large complement of Wnt ligands, with 11 Wnt genes present. In this study, the expression of each of these genes was examined across a range of stages during embryonic development. The expression of Wnt genes in Strigamia displays much variability. Most Wnt genes are expressed in segmental stripes in the trunk; near the proctodeum; and in the head region. However, despite this overall broad similarity, there are many differences between the various Wnt genes in their exact patterns of expression. These data should be considered in the context of different hypotheses regarding the functional relationships between the Wnt genes and the degree of redundancy present in this system. The findings of this study are consistent with one particular model of Wnt activity, the combinatorial model, whereby the combination of Wnt ligands present in a particular region defines its identity. These findings should also be useful in attempts to reconstruct the evolutionary history of Wnt signaling in arthropods.

Expression patterns of Wnt genes in the venom claws of centipedes.

The venom claws of centipedes, also known as forcipules, represent an evolutionary novelty that must have arisen in the centipede stem species, as they are not found in any other myriapods. The developmental-genetic changes that are involved in the origin of novelties are of considerable interest. It has previously been shown that centipede forcipules have a unique Hox code. However, this is a combinatorial code: no single Hox gene has a forcipule-specific expression. Here, we focus on Wnt genes. Two genes of this family show forcipule-specific expression in the "model centipede" Strigamia maritima: Wnt7 and Wnt11. For Wnt7, this forcipular expression zone seems to be a new one, which has arisen in evolution subsequently to other expression zones of the same gene. However, for Wnt11, the forcipule-specific expression probably arose by reduction of a more general pattern that originally included most or all of the limbs of an ancestral myriapod. Thus the developmental-genetic basis of the evolutionary change that turned the first pair of walking legs into venom claws is complex, involving different types of change in expression pattern. This sort of complexity is likely to be the case regarding evolutionary changes in morphology in general. Whether the origins of those features that can be considered as novelties are different in terms of their developmental-genetic basis from more routine evolutionary changes remains an open question.

Development of the venom ducts in the centipede Scolopendra: an example of recapitulation.

In contrast to previous claims that (a) there is a law of recapitulation and, conversely, (b) recapitulation never happens, the evolutionary repatterning of development can take many forms, of which recapitulation is one. Here, we add another example to the list of case studies of recapitulation. This example involves the development of the venom claws (forcipules) in the centipede Scolopendra subspinipes mutilans, and in particular the development of the duct through which venom flows from the gland that produces it (proximal) to the opening called the meatus (distal) through which it is injected into prey. Most of the information we present is from early postembryonic stages--these have been neglected in previous work on centipede development. We show that the venom ducts arise from sutures that are invaginations of the cuticle. In S. s. mutilans, the invagination in each forcipule forms into a tubular structure that detaches itself from the exoskeleton and moves toward the center of the forcipule. This is in contrast to extant Scutigera, and also, probably, Scolopendra's extinct Scutigera-like ancestors, where the duct remains attached to the cuticle of throughout development. Thus, S. s. mutilans exhibits a recapitulatory repatterning of development.

Manipulating the nature of embryonic mitotic waves.

Early embryogenesis is characterized by rapid and synchronous cleavage divisions, which are often controlled by wave-like patterns of Cdk1 activity. Two mechanisms have been proposed for mitotic waves: sweep and trigger waves.1,2 The two mechanisms give rise to different wave speeds, dependencies on physical and molecular parameters, and spatial profiles of Cdk1 activity: upward sweeping gradients versus traveling wavefronts. Both mechanisms hinge on the transient bistability governing the cell cycle and are differentiated by the speed of the cell-cycle progression: sweep and trigger waves arise for rapid and slow drives, respectively. Here, using quantitative imaging of Cdk1 activity and theory, we illustrate that sweep waves are the dominant mechanism in Drosophila embryos and test two fundamental predictions on the transition from sweep to trigger waves. We demonstrate that sweep waves can be turned into trigger waves if the cell cycle is slowed down genetically or if significant delays in the cell-cycle progression are introduced across the embryo by altering nuclear density. Our genetic experiments demonstrate that Polo kinase is a major rate-limiting regulator of the blastoderm divisions, and genetic perturbations reducing its activity can induce the transition from sweep to trigger waves. Furthermore, we show that changes in temperature cause an essentially uniform slowdown of interphase and mitosis. That results in sweep waves being observed across a wide temperature range despite the cell-cycle durations being significantly different. Collectively, our combination of theory and experiments elucidates the nature of mitotic waves in Drosophila embryogenesis, their control mechanisms, and their mutual transitions.

Cullin-5 mutants reveal collective sensing of the nucleocytoplasmic ratio in Drosophila embryogenesis.

In most metazoans, early embryonic development is characterized by rapid division cycles that pause before gastrulation at the midblastula transition (MBT).1 These cleavage divisions are accompanied by cytoskeletal rearrangements that ensure proper nuclear positioning. However, the molecular mechanisms controlling nuclear positioning are not fully elucidated. In Drosophila, early embryogenesis unfolds in a multinucleated syncytium. Nuclei rapidly move across the anterior-posterior (AP) axis at cell cycles 4-6 in a process driven by actomyosin contractility and cytoplasmic flows.2,3 In shackleton (shkl) mutants, this axial spreading is impaired.4 Here, we show that shkl mutants carry mutations in the cullin-5 (cul-5) gene. Live imaging experiments show that Cul-5 is downstream of the cell cycle but is required for cortical actomyosin contractility. The nuclear spreading phenotype of cul-5 mutants can be rescued by reducing Src activity, suggesting that a major target of cul-5 is Src kinase. cul-5 mutants display gradients of nuclear density across the AP axis that we exploit to study cell-cycle control as a function of the N/C ratio. We found that the N/C ratio is sensed collectively in neighborhoods of about 100 µm, and such collective sensing is required for a precise MBT, in which all the nuclei in the embryo pause their division cycle. Moreover, we found that the response to the N/C ratio is slightly graded along the AP axis. These two features can be linked to Cdk1 dynamics. Collectively, we reveal a new pathway controlling nuclear positioning and provide a dissection of how nuclear cycles respond to the N/C ratio.

Developmental variability channels mouse molar evolution.

Do developmental systems preferentially produce certain types of variation that orient phenotypic evolution along preferred directions? At different scales, from the intra-population to the interspecific, the murine first upper molar shows repeated anterior elongation. Using a novel quantitative approach to compare the development of two mouse strains with short or long molars, we identified temporal, spatial and functional differences in tooth signaling center activity, that arise from differential tuning of the activation-inhibition mechanisms underlying tooth patterning. By tracing their fate, we could explain why only the upper first molar reacts via elongation of its anterior part. Despite a lack of genetic variation, individuals of the elongated strain varied in tooth length and the temporal dynamics of their signaling centers, highlighting the intrinsic instability of the upper molar developmental system. Collectively, these results reveal the variational properties of murine molar development that drive morphological evolution along a line of least resistance.

Over time species develop random mutations in their genetic sequence that causes their form to change. If this new form increases the survival of a species it will become favored through natural selection and is more likely to get passed on to future generations. But, the evolution of these new traits also depends on what happens during development. Developmental mechanisms control how an embryo progresses from a single cell to an adult organism made of many cells. Mutations that alter these processes can influence the physical outcome of development, and cause a new trait to form. This means that if many different mutations alter development in a similar way, this can lead to the same physical change, making it 'easy' for a new trait to repeatedly occur. Most of the research has focused on finding the mutations that underlie repeated evolution, but rarely on identifying the role of the underlying developmental mechanisms. To bridge this gap, Hayden et al. investigated how changes during development influence the shape and size of molar teeth in mice. In some wild species of mice, the front part of the first upper molar is longer than in other species. This elongation, which is repeatedly found in mice from different islands, likely came from developmental mechanisms. Tooth development in mice has been well-studied in the laboratory, and Hayden et al. started by identifying two strains of laboratory mice that mimic the teeth seen in their wild cousins, one with elongated upper first molars and another with short ones. Comparing how these two strains of mice developed their elongated or short teeth revealed key differences in the embryonic structures that form the upper molar and cause it to elongate. Further work showed that variations in these embryonic structures can even cause mice that are genetically identical to have longer or shorter upper first molars. These findings show how early differences during development can lead to small variations in form between adult species of mice. This study highlights how studying developmental differences as well as genetic sequences can further our understanding of how different species evolved.