Publications

Tissue morphogenesis requires the collective, coordinated motion and deformation of a large number of cells. Vertex model simulations for tissue mechanics have been developed to bridge the scales between force generation at the cellular level and tissue deformation and flows. We review here various formulations of vertex models that have been proposed for describing tissues in two and three dimensions. We discuss a generic formulation using a virtual work differential, and we review applications of vertex models to biological morphogenetic processes. We also highlight recent efforts to obtain continuum theories of tissue mechanics, which are effective, coarse-grained descriptions of vertex models.This article is part of the themed issue 'Systems morphodynamics: understanding the development of tissue hardware'.

Most mammals, birds, and reptiles are readily recognized by their hairs, feathers, and scales, respectively. However, the lack of fossil intermediate forms between scales and hairs and substantial differences in their morphogenesis and protein composition have fueled the controversy pertaining to their potential common ancestry for decades. Central to this debate is the apparent lack of an "anatomical placode" (that is, a local epidermal thickening characteristic of feathers' and hairs' early morphogenesis) in reptile scale development. Hence, scenarios have been proposed for the independent development of the anatomical placode in birds and mammals and parallel co-option of similar signaling pathways for their morphogenesis. Using histological and molecular techniques on developmental series of crocodiles and snakes, as well as of unique wild-type and EDA (ectodysplasin A)-deficient scaleless mutant lizards, we show for the first time that reptiles, including crocodiles and squamates, develop all the characteristics of an anatomical placode: columnar cells with reduced proliferation rate, as well as canonical spatial expression of placode and underlying dermal molecular markers. These results reveal a new evolutionary scenario where hairs, feathers, and scales of extant species are homologous structures inherited, with modification, from their shared reptilian ancestor's skin appendages already characterized by an anatomical placode and associated signaling molecules.

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.

Understanding how speciation relates to ecological divergence has long fascinated biologists. It is assumed that ecological divergence is essential to sympatric speciation, as a mechanism to avoid competition and eventually lead to reproductive isolation, while divergence in allopatry is not necessarily associated with niche differentiation. The impact of the spatial context of divergence on the evolutionary rates of abiotic dimensions of the ecological niche has rarely been explored for an entire clade. Here, we compare the magnitude of climatic niche shifts between sympatric versus allopatric divergence of lineages in butterflies. By combining next-generation sequencing, parametric biogeography and ecological niche analyses applied to a genus-wide phylogeny of Palaearctic butterflies, we compare evolutionary rates along eight climatic dimensions across sister lineages that diverged in large-scale sympatry versus allopatry. In order to examine the possible effects of the spatial scale at which sympatry is defined, we considered three sets of biogeographic assignments, ranging from narrow to broad definition. Our findings suggest higher rates of niche evolution along all climatic dimensions for sister lineages that diverge in sympatry, when using a narrow delineation of biogeographic areas. This result contrasts with significantly lower rates of climatic niche evolution found in cases of allopatric speciation, despite the biogeographic regions defined here being characterized by significantly different climates. Higher rates in allopatry are retrieved when biogeographic areas are too widely defined-in such a case allopatric events may be recorded as sympatric. Our results reveal the macro-evolutionary significance of abiotic niche differentiation involved in speciation processes within biogeographic regions, and illustrate the importance of the spatial scale chosen to define areas when applying parametric biogeographic analyses.

The transition from benthos to plankton requires multiple adaptations, yet so far it remains unclear how these are acquired in the course of the transition. To investigate this process, we analyzed the genetic diversity and distribution patterns of a group of foraminifera of the genus Bolivina with a tychopelagic mode of life (same species occurring both in benthos and plankton). We assembled a global sequence data set for this group from single-cell DNA extractions and occurrences in metabarcodes from pelagic environmental samples. The pelagic sequences all cluster within a single monophyletic clade within Bolivina. This clade harbors three distinct genetic lineages, which are associated with incipient morphological differentiation. All lineages occur in the plankton and benthos, but only one lineage exhibits no limit to offshore dispersal and has been shown to grow in the plankton. These observations indicate that the emergence of buoyancy regulation within the clade preceded the evolution of pelagic feeding and that the evolution of both traits was not channeled into a full transition into the plankton. We infer that in foraminifera, colonization of the planktonic niche may occur by sequential cooptation of independently acquired traits, with holoplanktonic species being recruited from tychopelagic ancestors.

Several attempts to include post-Devonian lungfish in phylogenies of dipnoans have been made, but are hampered by the poor preservation of most Mesozoic and Cenozoic lungfish, and by the paucity of the occurrences of these taxa. This contribution has made use of the few post-Devonian fossils that are known from cranial, dental and post-cranial remains, to compare them with Devonian material and with living lungfish. Characters have been chosen to cover the best preserved structures of the Mesozoic and Cenozoic fossils, and the resulting phylogeny has been discussed in relation to previous analyses. The post-Devonian phylogeny has been anchored to a published phylogeny of Devonian taxa and a phylogenetic diversity curve has been computed. Based on this phylogeny rates of origination and extinction have been calculated, and environmental transfers and trends in body size changes have been detected. Our analyses show that the Permian gnathorhizodontids and the extant lepidosirenids are closely related, that lungfishes have experienced two phases of taxic diversification, a marine one in the Devonian and a freshwater one in the Permian. Major negative size shifts occurred during the marine phase and major positive size shifts occurred during the freshwater phase. A new classification of post-Devonian dipnoan families is also presented.

Rotaliida is the most diversified order of modern benthic foraminifera, comprising 73 extant families. The current classification of Rotaliida is based exclusively on morphological characters of their calcareous tests. Previous molecular phylogenetic studies partly confirmed morphological classifications but they were based on relatively limited taxon sampling and SSU rRNA gene sequences only. Here, we investigate phylogenetic relationships of rotaliid foraminifera by sequencing complete or partial SSU and LSU rRNA genes for 80 phylotypes representing 70 genera. The analysis of 87 sequences of concatenated genes enabled us to recognize seven rotaliid superfamilies, four of them being established ones (Planorbulinoidea, Discorboidea, Rotalioidea, Nummulitoidea) and three new ones (Glabratelloidea, Calcarinoidea, Serioidea). Although some of these superfamilies were not strongly supported in phylogenetic trees, due to a weak signal of ribosomal genes, they comprise genera and families of distinctive morphological features that have been often united in previous morphology-based classifications. A total of 34 families has been recognized, 10 of them being newly established ones. The placement of three families (Rubratelliidae, Cymbaloporidae, Murrayinelliidae), which occupied isolated branches among the seven superfamilies, could not be determined and they are considered here as incertae sedis. An assemblage of 22 genera representing several morphologically different families of small rotaliids (Clade 3) was not raised to superfamily level because of the lack of molecular and/or morphological support. The phylogenetic relations within this assemblage as well as within the superfamily Rotalioidea have been examined in detail by analyzing partial SSU rDNA sequences for large taxon sampling.

During embryonic development, mechanical forces are essential for cellular rearrangements driving tissue morphogenesis. Here, we show that in the early zebrafish embryo, friction forces are generated at the interface between anterior axial mesoderm (prechordal plate, ppl) progenitors migrating towards the animal pole and neurectoderm progenitors moving in the opposite direction towards the vegetal pole of the embryo. These friction forces lead to global rearrangement of cells within the neurectoderm and determine the position of the neural anlage. Using a combination of experiments and simulations, we show that this process depends on hydrodynamic coupling between neurectoderm and ppl as a result of E-cadherin-mediated adhesion between those tissues. Our data thus establish the emergence of friction forces at the interface between moving tissues as a critical force-generating process shaping the embryo.

We describe four new species of Flabellula, Leptomyxa and Rhizamoeba and publish new SSU rRNA gene and actin gene sequences of leptomyxids. Using these data we provide the most comprehensive SSU phylogeny of leptomyxids to date. Based on the analyses of morphological data and results of the SSU rRNA gene phylogeny we suggest changes in the systematics of the order Leptomyxida (Amoebozoa: Lobosa: Tubulinea). We propose to merge the genera Flabellula and Paraflabellula (the genus Flabellula remains valid by priority rule). The genus Rhizamoeba is evidently polyphyletic in all phylogenetic trees; we suggest retaining the generic name Rhizamoeba for the group unifying R. saxonica, R.matisi n. sp. and R. polyura, the latter remains the type species of the genus Rhizamoeba. Based on molecular and morphological evidence we move all remaining Rhizamoeba species to the genus Leptomyxa. New family Rhizamoebidae is established here in order to avoid paraphyly of the family Leptomyxidae. With the suggested changes both molecular and morphological systems of the order Leptomyxida are now fully congruent to each other.

Current biodiversity assessment and biomonitoring are largely based on the morphological identification of selected bioindicator taxa. Recently, several attempts have been made to use eDNA metabarcoding as an alternative tool. However, until now, most applied metabarcoding studies have been based on the taxonomic assignment of sequences that provides reference to morphospecies ecology. Usually, only a small portion of metabarcoding data can be used due to a limited reference database and a lack of phylogenetic resolution. Here, we investigate the possibility to overcome these limitations by using a taxonomy-free approach that allows the computing of a molecular index directly from eDNA data without any reference to morphotaxonomy. As a case study, we use the benthic diatoms index, commonly used for monitoring the biological quality of rivers and streams. We analysed 87 epilithic samples from Swiss rivers, the ecological status of which was established based on the microscopic identification of diatom species. We compared the diatom index derived from eDNA data obtained with or without taxonomic assignment. Our taxonomy-free approach yields promising results by providing a correct assessment for 77% of examined sites. The main advantage of this method is that almost 95% of OTUs could be used for index calculation, compared to 35% in the case of the taxonomic assignment approach. Its main limitations are under-sampling and the need to calibrate the index based on the microscopic assessment of diatoms communities. However, once calibrated, the taxonomy-free molecular index can be easily standardized and applied in routine biomonitoring, as a complementary tool allowing fast and cost-effective assessment of the biological quality of watercourses. This article is protected by copyright. All rights reserved.