The Unreasonable Effectiveness of Reaction Diffusion in Vertebrate Skin Color Patterning.

In 1952, Alan Turing published the reaction-diffusion (RD) mathematical framework, laying the foundations of morphogenesis as a self-organized process emerging from physicochemical first principles. Regrettably, this approach has been widely doubted in the field of developmental biology. First, we summarize Turing's line of thoughts to alleviate the misconception that RD is an artificial mathematical construct. Second, we discuss why phenomenological RD models are particularly effective for understanding skin color patterning at the meso/macroscopic scales, without the need to parameterize the profusion of variables at lower scales. More specifically, we discuss how RD models (a) recapitulate the diversity of actual skin patterns, (b) capture the underlying dynamics of cellular interactions, (c) interact with tissue size and shape, (d) can lead to ordered sequential patterning, (e) generate cellular automaton dynamics in lizards and snakes, (f) predict actual patterns beyond their statistical features, and (g) are robust to model variations. Third, we discuss the utility of linear stability analysis and perform numerical simulations to demonstrate how deterministic RD emerges from the underlying chaotic microscopic agents.

A living mesoscopic cellular automaton made of skin scales.

In vertebrates, skin colour patterns emerge from nonlinear dynamical microscopic systems of cell interactions. Here we show that in ocellated lizards a quasi-hexagonal lattice of skin scales, rather than individual chromatophore cells, establishes a green and black labyrinthine pattern of skin colour. We analysed time series of lizard scale colour dynamics over four years of their development and demonstrate that this pattern is produced by a cellular automaton (a grid of elements whose states are iterated according to a set of rules based on the states of neighbouring elements) that dynamically computes the colour states of individual mesoscopic skin scales to produce the corresponding macroscopic colour pattern. Using numerical simulations and mathematical derivation, we identify how a discrete von Neumann cellular automaton emerges from a continuous Turing reaction-diffusion system. Skin thickness variation generated by three-dimensional morphogenesis of skin scales causes the underlying reaction-diffusion dynamics to separate into microscopic and mesoscopic spatial scales, the latter generating a cellular automaton. Our study indicates that cellular automata are not merely abstract computational systems, but can directly correspond to processes generated by biological evolution.

Genome mapping of a LYST mutation in corn snakes indicates that vertebrate chromatophore vesicles are lysosome-related organelles.

Reptiles exhibit a spectacular diversity of skin colors and patterns brought about by the interactions among three chromatophore types: black melanophores with melanin-packed melanosomes, red and yellow xanthophores with pteridine- and/or carotenoid-containing vesicles, and iridophores filled with light-reflecting platelets generating structural colors. Whereas the melanosome, the only color-producing endosome in mammals and birds, has been documented as a lysosome-related organelle, the maturation paths of xanthosomes and iridosomes are unknown. Here, we first use 10x Genomics linked-reads and optical mapping to assemble and annotate a nearly chromosome-quality genome of the corn snake Pantherophis guttatus The assembly is 1.71 Gb long, with an N50 of 16.8 Mb and L50 of 24. Second, we perform mapping-by-sequencing analyses and identify a 3.9-Mb genomic interval where the lavender variant resides. The lavender color morph in corn snakes is characterized by gray, rather than red, blotches on a pink, instead of orange, background. Third, our sequencing analyses reveal a single nucleotide polymorphism introducing a premature stop codon in the lysosomal trafficking regulator gene (LYST) that shortens the corresponding protein by 603 amino acids and removes evolutionary-conserved domains. Fourth, we use light and transmission electron microscopy comparative analyses of wild type versus lavender corn snakes and show that the color-producing endosomes of all chromatophores are substantially affected in the LYST mutant. Our work provides evidence characterizing xanthosomes in xanthophores and iridosomes in iridophores as lysosome-related organelles.

Lizard Skin Patterns and the Ising Model.

The ocellated lizard (Timon lepidus) exhibits an intricate skin color pattern made of monochromatic black and green skin scales, whose dynamics of color flipping are known to be well modeled by a stochastic cellular automaton. We show that the late-time probability distribution of the pattern corresponds to the canonical probability distribution of the antiferromagnetic Ising model and can be generated by dynamics different from the commonly-used Glauber. We comment on skin scale patterns generated by the Ising model on the triangular lattice in the low-temperature limit.

Feather arrays are patterned by interacting signalling and cell density waves.

Feathers are arranged in a precise pattern in avian skin. They first arise during development in a row along the dorsal midline, with rows of new feather buds added sequentially in a spreading wave. We show that the patterning of feathers relies on coupled fibroblast growth factor (FGF) and bone morphogenetic protein (BMP) signalling together with mesenchymal cell movement, acting in a coordinated reaction-diffusion-taxis system. This periodic patterning system is partly mechanochemical, with mechanical-chemical integration occurring through a positive feedback loop centred on FGF20, which induces cell aggregation, mechanically compressing the epidermis to rapidly intensify FGF20 expression. The travelling wave of feather formation is imposed by expanding expression of Ectodysplasin A (EDA), which initiates the expression of FGF20. The EDA wave spreads across a mesenchymal cell density gradient, triggering pattern formation by lowering the threshold of mesenchymal cells required to begin to form a feather bud. These waves, and the precise arrangement of feather primordia, are lost in the flightless emu and ostrich, though via different developmental routes. The ostrich retains the tract arrangement characteristic of birds in general but lays down feather primordia without a wave, akin to the process of hair follicle formation in mammalian embryos. The embryonic emu skin lacks sufficient cells to enact feather formation, causing failure of tract formation, and instead the entire skin gains feather primordia through a later process. This work shows that a reaction-diffusion-taxis system, integrated with mechanical processes, generates the feather array. In flighted birds, the key role of the EDA/Ectodysplasin A receptor (EDAR) pathway in vertebrate skin patterning has been recast to activate this process in a quasi-1-dimensional manner, imposing highly ordered pattern formation.

Translucent in air and iridescent in water: structural analysis of a salamander egg sac.

Females of some Asian salamanders of the genus Hynobius deposit in streams their eggs embedded in a translucent envelope called an 'egg sac'. The edges of the envelope exhibit a spectacular blue-to-yellow iridescent glow, which instantaneously disappears when the sac is removed from water. First, our scanning electron microscopy analyses reveal that the inner surface of the 100 µm-thick envelope displays striations (length scale of about 3 µm), which are themselves covered by much smaller (190 ± 30 nm) and quasi-periodic corrugations. The latter could constitute a surface diffraction grating generating iridescence by light interference. Second, our transmission electron microscopy and focused-ion-beam scanning electron microscopy analyses show that the bulk of the egg sac wall is composed of meandering fibres with a quasi-periodic modulation of 190 ± 60 nm along the thickness of the envelope, generating a photonic crystal. Third, Fourier power analyses of 450 electron microscopy images with varying incident angles indicate that changing the surrounding medium from water to air shifts most of the backscattered power spectrum to the ultraviolet range, hence, explaining that the egg sac loses visible iridescence when removed out of the water. Fourth, the results of our photography and optical spectroscopy experiments of submerged and emerged egg sacs rule out the possibility that the iridescence is due to a thin film or a multilayer, whereas the observed non-specular response is compatible with the backscattering expected from surface diffraction gratings and volumetric photonic crystals with spatial 1D modulation. Finally, although we mention several potential biological functions of the egg sac structural colours and iridescence, we emphasise that these optical properties might be the by-products of the envelope material internal structure selected during evolution for its mechanical properties.

Phylogenetic mapping of scale nanostructure diversity in snakes.

BACKGROUND: Many species of snakes exhibit epidermal surface nanostructures that form complex motifs conferring self-cleaning properties, and sometimes structural iridescence, to their skin. RESULTS: Using confocal microscopy, we show that these specialised cells can be greatly elongated along their left-right axis and that different types of nanostructures are generated by cell borders and cell surface. To characterise the complexity and diversity of these surface gratings, we analysed scanning electron microscopy images of skin sheds from 353 species spanning 19 of the 26 families of snakes and characterised the observed nanostructures with four characters. The full character matrix, as well as one representative SEM image of each of the corresponding species, is available as a MySQL relational database at https://snake-nanogratings.lanevol.org . We then performed continuous-time Markov phylogenetic mapping on the snake phylogeny, providing an evolutionary dynamical estimate for the different types of nanostructures. These analyses suggest that the presence of cell border digitations is the ancestral state for snake skin nanostructures which was subsequently and independently lost in multiple lineages. Our analyses also indicate that cell shape and cell border shape are co-dependent characters whereas we did not find correlation between a simple life habit classification and any specific nanomorphological character. CONCLUSIONS: These results, compatible with the fact that multiple types of nanostructures can generate hydrophobicity, suggest that the diversity and complexity of snake skin surface nano-morphology are dominated by phylogenetic rather than habitat-specific functional constraints. The present descriptive study opens the perspective of investigating the cellular self-organisational cytoskeletal processes controlling the patterning of different skin surface nanostructures in snakes and lizards.

Bifurcation Analysis of Reaction Diffusion Systems on Arbitrary Surfaces.

In this paper, we present computational techniques to investigate the effect of surface geometry on biological pattern formation. In particular, we study two-component, nonlinear reaction-diffusion (RD) systems on arbitrary surfaces. We build on standard techniques for linear and nonlinear analysis of RD systems and extend them to operate on large-scale meshes for arbitrary surfaces. In particular, we use spectral techniques for a linear stability analysis to characterise and directly compose patterns emerging from homogeneities. We develop an implementation using surface finite element methods and a numerical eigenanalysis of the Laplace-Beltrami operator on surface meshes. In addition, we describe a technique to explore solutions of the nonlinear RD equations using numerical continuation. Here, we present a multiresolution approach that allows us to trace solution branches of the nonlinear equations efficiently even for large-scale meshes. Finally, we demonstrate the working of our framework for two RD systems with applications in biological pattern formation: a Brusselator model that has been used to model pattern development on growing plant tips, and a chemotactic model for the formation of skin pigmentation patterns. While these models have been used previously on simple geometries, our framework allows us to study the impact of arbitrary geometries on emerging patterns.

Retroviral envelope syncytin capture in an ancestrally diverged mammalian clade for placentation in the primitive Afrotherian tenrecs.

Syncytins are fusogenic envelope (env) genes of retroviral origin that have been captured for a function in placentation. Syncytins have been identified in Euarchontoglires (primates, rodents, Leporidae) and Laurasiatheria (Carnivora, ruminants) placental mammals. Here, we searched for similar genes in species that retained characteristic features of primitive mammals, namely the Malagasy and mainland African Tenrecidae. They belong to the superorder Afrotheria, an early lineage that diverged from Euarchotonglires and Laurasiatheria 100 Mya, during the Cretaceous terrestrial revolution. An in silico search for env genes with full coding capacity within a Tenrecidae genome identified several candidates, with one displaying placenta-specific expression as revealed by RT-PCR analysis of a large panel of Setifer setosus tissues. Cloning of this endogenous retroviral env gene demonstrated fusogenicity in an ex vivo cell-cell fusion assay on a panel of mammalian cells. Refined analysis of placental architecture and ultrastructure combined with in situ hybridization demonstrated specific expression of the gene in multinucleate cellular masses and layers at the materno-fetal interface, consistent with a role in syncytium formation. This gene, which we named "syncytin-Ten1," is conserved among Tenrecidae, with evidence of purifying selection and conservation of fusogenic activity. To our knowledge, it is the first syncytin identified to date within the ancestrally diverged Afrotheria superorder.

Changes in Hox genes' structure and function during the evolution of the squamate body plan.

Hox genes are central to the specification of structures along the anterior-posterior body axis, and modifications in their expression have paralleled the emergence of diversity in vertebrate body plans. Here we describe the genomic organization of Hox clusters in different reptiles and show that squamates have accumulated unusually large numbers of transposable elements at these loci, reflecting extensive genomic rearrangements of coding and non-coding regulatory regions. Comparative expression analyses between two species showing different axial skeletons, the corn snake and the whiptail lizard, revealed major alterations in Hox13 and Hox10 expression features during snake somitogenesis, in line with the expansion of both caudal and thoracic regions. Variations in both protein sequences and regulatory modalities of posterior Hox genes suggest how this genetic system has dealt with its intrinsic collinear constraint to accompany the substantial morphological radiation observed in this group.

Phylogeography and support vector machine classification of colour variation in panther chameleons.

Lizards and snakes exhibit colour variation of adaptive value for thermoregulation, camouflage, predator avoidance, sexual selection and speciation. Furcifer pardalis, the panther chameleon, is one of the most spectacular reptilian endemic species in Madagascar, with pronounced sexual dimorphism and exceptionally large intraspecific variation in male coloration. We perform here an integrative analysis of molecular phylogeography and colour variation after collecting high-resolution colour photographs and blood samples from 324 F. pardalis individuals in locations spanning the whole species distribution. First, mitochondrial and nuclear DNA sequence analyses uncover strong genetic structure among geographically restricted haplogroups, revealing limited gene flow among populations. Bayesian coalescent modelling suggests that most of the mitochondrial haplogroups could be considered as separate species. Second, using a supervised multiclass support vector machine approach on five anatomical components, we identify patterns in 3D colour space that efficiently predict assignment of male individuals to mitochondrial haplogroups. We converted the results of this analysis into a simple visual classification key that can assist trade managers to avoid local population overharvesting.

Precise colocalization of interacting structural and pigmentary elements generates extensive color pattern variation in Phelsuma lizards.

BACKGROUND: Color traits in animals play crucial roles in thermoregulation, photoprotection, camouflage, and visual communication, and are amenable to objective quantification and modeling. However, the extensive variation in non-melanic pigments and structural colors in squamate reptiles has been largely disregarded. Here, we used an integrated approach to investigate the morphological basis and physical mechanisms generating variation in color traits in tropical day geckos of the genus Phelsuma. RESULTS: Combining histology, optics, mass spectrometry, and UV and Raman spectroscopy, we found that the extensive variation in color patterns within and among Phelsuma species is generated by complex interactions between, on the one hand, chromatophores containing yellow/red pteridine pigments and, on the other hand, iridophores producing structural color by constructive interference of light with guanine nanocrystals. More specifically, we show that 1) the hue of the vivid dorsolateral skin is modulated both by variation in geometry of structural, highly ordered narrowband reflectors, and by the presence of yellow pigments, and 2) that the reflectivity of the white belly and of dorsolateral pigmentary red marks, is increased by underlying structural disorganized broadband reflectors. Most importantly, these interactions require precise colocalization of yellow and red chromatophores with different types of iridophores, characterized by ordered and disordered nanocrystals, respectively. We validated these results through numerical simulations combining pigmentary components with a multilayer interferential optical model. Finally, we show that melanophores form dark lateral patterns but do not significantly contribute to variation in blue/green or red coloration, and that changes in the pH or redox state of pigments provide yet another source of color variation in squamates. CONCLUSIONS: Precisely colocalized interacting pigmentary and structural elements generate extensive variation in lizard color patterns. Our results indicate the need to identify the developmental mechanisms responsible for the control of the size, shape, and orientation of nanocrystals, and the superposition of specific chromatophore types. This study opens up new perspectives on Phelsuma lizards as models in evolutionary developmental biology.

Development and embryonic staging in non-model organisms: the case of an afrotherian mammal.

Studies of evolutionary developmental biology commonly use 'model organisms' such as fruit flies or mice, and questions are often functional or epigenetic. Phylogenetic investigations, in contrast, typically use species that are less common and mostly deal with broad scale analyses in the tree of life. However, important evolutionary transformations have taken place at all taxonomic levels, resulting in such diverse forms as elephants and shrews. To understand the mechanisms underlying morphological diversification, broader sampling and comparative approaches are paramount. Using a simple, standardized protocol, we describe for the first time the development of soft tissues and some parts of the skeleton, using µCT-imaging of developmental series of Echinops telfairi and Tenrec ecaudatus, two tenrecid afrotherian mammals. The developmental timing of soft tissue and skeletal characters described for the tenrecids is briefly compared with that of other mammals, including mouse, echidna, and the opossum. We found relatively few heterochronic differences in development in the armadillo vs. tenrec, consistent with a close relationship of Xenarthra and Afrotheria. Ossification in T. ecaudatus continues well into the second half of overall gestation, resembling the pattern seen in other small mammals and differing markedly from the advanced state of ossification evident early in the gestation of elephants, sheep, and humans.

Transient agonism of the sonic hedgehog pathway triggers a permanent transition of skin appendage fate in the chicken embryo.

Vertebrate skin appendage early development is mediated by conserved molecular signaling composing a dynamical reaction-diffusion-like system. Variations to such systems contribute to the remarkable diversity of skin appendage forms within and among species. Here, we demonstrate that stage-specific transient agonism of sonic hedgehog (Shh) pathway signaling in chicken triggers a complete and permanent transition from reticulate scales to feathers on the ventral surfaces of the foot and digits. Resulting ectopic feathers are developmentally comparable to feathers adorning the body, with down-type feathers transitioning into regenerative, bilaterally symmetric contour feathers in adult chickens. Crucially, this spectacular transition of skin appendage fate (from nodular reticulate scales to bona fide adult feathers) does not require sustained treatment. Our RNA sequencing analyses confirm that smoothened agonist treatment specifically promotes the expression of key Shh pathway-associated genes. These results indicate that variations in Shh pathway signaling likely contribute to the natural diversity and regionalization of avian integumentary appendages.

High-resolution confocal and light-sheet imaging of collagen 3D network architecture in very large samples.

Although notoriously difficult, imaging collagen network architecture, a key element affecting tissue mechanical properties, is of paramount importance in developmental and cancer biology. Here, we introduce a simple and robust method of whole-mount collagen staining with the 'Fast Green' dye that provides unmatched visualization of collagen 3D network architecture, via confocal or light-sheet microscopy, compatible with solvent-based tissue clearing and immunostaining.

Protocol for the rapid intravenous in ovo injection of developing amniote embryos.

We present a technique for precise drug delivery into the vascular system of developing amniote embryos via injection into chorioallantoic veins underlying the eggshell membrane. We describe steps for incubating and candling eggs, removing the shell to expose underlying veins, and precise intravenous injection. In addition to chicken embryos, this protocol is applicable to other amniote species that lay hard-shell eggs, including crocodiles and tortoises. This technique is rapid, is reproducible, is of low cost, and will provide an important resource for developmental biologists. For complete details on the use and execution of this protocol, please refer to Cooper & Milinkovitch.1.

Somitic positional information guides self-organized patterning of snake scales.

Two influential concepts in tissue patterning are Wolpert's positional information and Turing's self-organized reaction-diffusion (RD). The latter establishes the patterning of hair and feathers. Here, our morphological, genetic, and functional-by CRISPR-Cas9-mediated gene disruption-characterization of wild-type versus "scaleless" snakes reveals that the near-perfect hexagonal pattern of snake scales is established through interactions between RD in the skin and somitic positional information. First, we show that ventral scale development is guided by hypaxial somites and, second, that ventral scales and epaxial somites guide the sequential RD patterning of the dorsolateral scales. The RD intrinsic length scale evolved to match somite periodicity, ensuring the alignment of ribs and scales, both of which play a critical role in snake locomotion.

Phylogenomics of unusual histone H2A Variants in Bdelloid rotifers.

Rotifers of Class Bdelloidea are remarkable in having evolved for millions of years, apparently without males and meiosis. In addition, they are unusually resistant to desiccation and ionizing radiation and are able to repair hundreds of radiation-induced DNA double-strand breaks per genome with little effect on viability or reproduction. Because specific histone H2A variants are involved in DSB repair and certain meiotic processes in other eukaryotes, we investigated the histone H2A genes and proteins of two bdelloid species. Genomic libraries were built and probed to identify histone H2A genes in Adineta vaga and Philodina roseola, species representing two different bdelloid families. The expressed H2A proteins were visualized on SDS-PAGE gels and identified by tandem mass spectrometry. We find that neither the core histone H2A, present in nearly all other eukaryotes, nor the H2AX variant, a ubiquitous component of the eukaryotic DSB repair machinery, are present in bdelloid rotifers. Instead, they are replaced by unusual histone H2A variants of higher mass. In contrast, a species of rotifer belonging to the facultatively sexual, desiccation- and radiation-intolerant sister class of bdelloid rotifers, the monogononts, contains a canonical core histone H2A and appears to lack the bdelloid H2A variant genes. Applying phylogenetic tools, we demonstrate that the bdelloid-specific H2A variants arose as distinct lineages from canonical H2A separate from those leading to the H2AX and H2AZ variants. The replacement of core H2A and H2AX in bdelloid rotifers by previously uncharacterized H2A variants with extended carboxy-terminal tails is further evidence for evolutionary diversity within this class of histone H2A genes and may represent adaptation to unusual features specific to bdelloid rotifers.

Modeling convergent scale-by-scale skin color patterning in multiple species of lizards.

Skin color patterning in vertebrates emerges at the macroscale from microscopic cell-cell interactions among chromatophores. Taking advantage of the convergent scale-by-scale skin color patterning dynamics in five divergent species of lizards, we quantify the respective efficiencies of stochastic (Lenz-Ising and cellular automata, sCA) and deterministic reaction-diffusion (RD) models to predict individual patterns and their statistical attributes. First, we show that all models capture the underlying microscopic system well enough to predict, with similar efficiencies, neighborhood statistics of adult patterns. Second, we show that RD robustly generates, in all species, a substantial gain in scale-by-scale predictability of individual adult patterns without the need to parametrize the system down to its many cellular and molecular variables. Third, using 3D numerical simulations and Lyapunov spectrum analyses, we quantitatively demonstrate that, given the non-linearity of the dynamical system, uncertainties in color measurements at the juvenile stage and in skin geometry variation explain most, if not all, of the residual unpredictability of adult individual scale-by-scale patterns. We suggest that the efficiency of RD is due to its intrinsic ability to exploit mesoscopic information such as continuous scale colors and the relations among growth, scales geometries, and the pattern length scale. Our results indicate that convergent evolution of CA patterning dynamics, leading to dissimilar macroscopic patterns in different species, is facilitated by their spontaneous emergence under a large range of RD parameters, as long as a Turing instability occurs in a skin domain with periodic thickness. VIDEO ABSTRACT.

Alien eggs in duck nests: brood parasitism or a help from Grandma?

Intraspecific brood parasitism (IBP) is a remarkable phenomenon by which parasitic females can increase their reproductive output by laying eggs in conspecific females' nests in addition to incubating eggs in their own nest. Kin selection could explain the tolerance, or even the selective advantage, of IBP, but different models of IBP based on game theory yield contradicting predictions. Our analyses of seven polymorphic autosomal microsatellites in two eider duck colonies indicate that relatedness between host and parasitizing females is significantly higher than the background relatedness within the colony. This result is unlikely to be a by-product of relatives nesting in close vicinity, as nest distance and genetic identity are not correlated. For eider females that had been ring-marked during the decades prior to our study, our analyses indicate that (i) the average age of parasitized females is higher than the age of nonparasitized females, (ii) the percentage of nests with alien eggs increases with the age of nesting females, (iii) the level of IBP increases with the host females' age, and (iv) the number of own eggs in the nest of parasitized females significantly decreases with age. IBP may allow those older females unable to produce as many eggs as they can incubate to gain indirect fitness without impairing their direct fitness: genetically related females specialize in their energy allocation, with young females producing more eggs than they can incubate and entrusting these to their older relatives. Intraspecific brood parasitism in ducks may constitute cooperation among generations of closely related females.