Laboratory of Artificial & Natural Evolution

Michel C Milinkovitch

Full Professor

  • T: +41 22 379 33 38
  • office 4024b (Sciences III)

The Laboratory of Artificial & Natural Evolution — LANE

In a nutshell, our highly multidisciplinary team of biologists, bioinformaticians, physicists, computer scientists and mathematicians investigates the developmental and evolutionary mechanisms generating Life’s complexity and diversity (or maybe should I say ‘Life’s Beauty’). We investigate multiple non-classical model species, mainly reptiles and ‘exotic’ mammals, that can inform us on yet unknown exciting biological and physical processes generating this complex and diverse living world.

Central to our reasoning is that a proper understanding of the complexity and diversity of organismal forms cannot be achieved without integrating the physical constrains acting on the developmental and Darwinian processes. Our research requires integrating data and methods from comparative genomics, molecular developmental genetics, as well as physical experiments, mathematical modelling and numerical simulations. More specifically, we investigate the interactions between physical (mechanics, reaction-diffusion) and biological (cell signalling, proliferation) mechanisms that generate and constrain the variety and complexity of skin appendages (scales, hairs, spines), skin colours (pigmentary and structural), and skin colour patterns in tetrapods (four-limbed vertebrates).

Much additional information is available on the LANE web site.


Subunits

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

    Annu Rev Cell Dev Biol 2023 Oct;39():145-174. 10.1146/annurev-cellbio-120319-024414.

    abstract

    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 () recapitulate the diversity of actual skin patterns, () capture the underlying dynamics of cellular interactions, () interact with tissue size and shape, () can lead to ordered sequential patterning, () generate cellular automaton dynamics in lizards and snakes, () predict actual patterns beyond their statistical features, and () 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.

    view more details on Pubmed

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

    Sci Adv 2023 Jun;9(24):eadf8834. PMC10266723. 10.1126/sciadv.adf8834.

    abstract

    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.

    view more details on Pubmed

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

    STAR Protoc 2023 May;4(2):102324. PMC10209871. 10.1016/j.xpro.2023.102324. S2666-1667(23)00291-5.

    abstract

    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..

    view more details on Pubmed

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

    Sci Adv 2023 May;9(20):eadg9619. 10.1126/sciadv.adg9619. PMC10191425.

    abstract

    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.

    view more details on Pubmed

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

    iScience 2023 Apr;26(4):106452. PMC10067766. 10.1016/j.isci.2023.106452. S2589-0042(23)00529-1.

    abstract

    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.

    view more details on Pubmed

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

    Curr Biol 2022 Nov;():. S0960-9822(22)01692-X. 10.1016/j.cub.2022.10.044.

    abstract

    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.

    view more details on Pubmed

  • Coupling Chemotaxis and Growth Poromechanics for the Modelling of Feather Primordia Patterning

    Mathematics. 2022; 10(21):4096. https://doi.org/10.3390/math10214096

    abstract

    In this paper we propose a new mathematical model for describing the complex interplay between skin cell populations with fibroblast growth factor and bone morphogenetic protein, occurring within deformable porous media describing feather primordia patterning. Tissue growth, in turn, modifies the transport of morphogens (described by reaction-diffusion equations) through diverse mechanisms such as advection from the solid velocity generated by mechanical stress, and mass supply. By performing an asymptotic linear stability analysis on the coupled poromechanical-chemotaxis system (assuming rheological properties of the skin cell aggregates that reside in the regime of infinitesimal strains and where the porous structure is fully saturated with interstitial fluid and encoding the coupling mechanisms through active stress) we obtain the conditions on the parameters—especially those encoding coupling mechanisms—under which the system will give rise to spatially heterogeneous solutions. We also extend the mechanical model to the case of incompressible poro-hyperelasticity and include the mechanisms of anisotropic solid growth and feedback by means of standard Lee decompositions of the tensor gradient of deformation. Because the model in question involves the coupling of several nonlinear PDEs, we cannot straightforwardly obtain closed-form solutions. We therefore design a suitable numerical method that employs backward Euler time discretisation, linearisation of the semidiscrete problem through Newton–Raphson’s method, a seven-field finite element formulation for the spatial discretisation, and we also advocate the construction and efficient implementation of tailored robust solvers. We present a few illustrative computational examples in 2D and 3D, briefly discussing different spatio-temporal patterns of growth factors as well as the associated solid response scenario depending on the specific poromechanical regime. Our findings confirm the theoretically predicted behaviour of spatio-temporal patterns, and the produced results reveal a qualitative agreement with respect to the expected experimental behaviour. We stress that the present study provides insight on several biomechanical properties of primordia patterning.

    see on external website

  • Lizard Skin Patterns and the Ising Model.

    Phys Rev Lett 2022 Jan;128(4):048102. 10.1103/PhysRevLett.128.048102.

    abstract

    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.

    view more details on Pubmed

  • Elephants evolved strategies reducing the biomechanical complexity of their trunk.

    Curr Biol 2021 Aug;():. S0960-9822(21)01133-7. 10.1016/j.cub.2021.08.029.

    abstract

    The elephant proboscis (trunk), which functions as a muscular hydrostat with a virtually infinite number of degrees of freedom, is a spectacular organ for delicate to heavy object manipulation as well as social and sensory functions. Using high-resolution motion capture and functional morphology analyses, we show here that elephants evolved strategies that reduce the biomechanical complexity of their trunk. Indeed, our behavioral experiments with objects of various shapes, sizes, and weights indicate that (1) complex behaviors emerge from the combination of a finite set of basic movements; (2) curvature, torsion, and strain provide an appropriate kinematic representation, allowing us to extract motion primitives from the trunk trajectories; (3) transport of objects involves the proximal propagation of an inward curvature front initiated at the tip; (4) the trunk can also form pseudo-joints for point-to-point motion; and (5) the trunk tip velocity obeys a power law with its path curvature, similar to human hand drawing movements. We also reveal with unprecedented precision the functional anatomy of the African and Asian elephant trunks using medical imaging and macro-scale serial sectioning, thus drawing strong connections between motion primitives and muscular synergies. Our study is the first combined quantitative analysis of the mechanical performance, kinematic strategies, and functional morphology of the largest animal muscular hydrostat on Earth. It provides data for developing innovative "soft-robotic" manipulators devoid of articulations, replicating the high compliance, flexibility, and strength of the elephant trunk. VIDEO ABSTRACT.

    view more details on Pubmed

  • Reaction-diffusion in a growing 3D domain of skin scales generates a discrete cellular automaton.

    Nat Commun 2021 04;12(1):2433. 10.1038/s41467-021-22525-1. 10.1038/s41467-021-22525-1.

    abstract

    We previously showed that the adult ocellated lizard skin colour pattern is effectively generated by a stochastic cellular automaton (CA) of skin scales. We additionally suggested that the canonical continuous 2D reaction-diffusion (RD) process of colour pattern development is transformed into this discrete CA by reduced diffusion coefficients at the borders of scales (justified by the corresponding thinning of the skin). Here, we use RD numerical simulations in 3D on realistic lizard skin geometries and demonstrate that skin thickness variation on its own is sufficient to cause scale-by-scale coloration and CA dynamics during RD patterning. In addition, we show that this phenomenon is robust to RD model variation. Finally, using dimensionality-reduction approaches on large networks of skin scales, we show that animal growth affects the scale-colour flipping dynamics by causing a substantial decrease of the relative length scale of the labyrinthine colour pattern of the lizard skin.

    view more details on Pubmed

  • Emergence of self-organizational patterning at the mesoscopic scale.

    Dev Cell 2021 Mar;56(6):719-721. S1534-5807(21)00200-8. 10.1016/j.devcel.2021.03.003.

    abstract

    In this issue of Developmental Cell, Chuyen et al. (2021) suggest that the Scf/Kit pathway controls mutual repulsion of multiciliated cells and their affinity for epidermal cell junctions through soluble and membrane-associated Scf ligands, respectively. Effective self-organizational patterning emerges at the mesoscopic scale as a small set of effective behaviors.

    view more details on Pubmed

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

    Proc Natl Acad Sci U S A 2020 Oct;():. 2003724117. 10.1073/pnas.2003724117.

    abstract

    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 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 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 () 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 mutant. Our work provides evidence characterizing xanthosomes in xanthophores and iridosomes in iridophores as lysosome-related organelles.

    view more details on Pubmed

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

    Soft Matter 2020 (16) Feb. 10.1039/c9sm02151e.

    abstract

    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.

    view more details on Pubmed

  • Elastic instability during branchial ectoderm development causes folding of the Chlamydosaurus erectile frill.

    Elife. 2019 Jun 25;8. pii: e44455. doi: 10.7554/eLife.44455

    abstract

    We study the morphogenesis and evolutionary origin of the spectacular erectile ruff of the frilled dragon (Chlamydosaurus kingii). Our comparative developmental analyses of multiple species suggest that the ancestor of Episquamata reptiles developed a neck fold from the hyoid branchial arch by preventing it to fully fuse with posterior arches. We also show that the Chlamydosaurus embryonic neck fold dramatically enlarges and its anterior surface wrinkles, establishing three convex ridges on each lobe of the frill. We suggest that this robust folding pattern is not due to localised increased growth at the positions of the ridges, but emerges from an elastic instability during homogeneous growth of the frill skin frustrated by its attachment to adjacent tissues. Our physical analog experiments and 3D computational simulations, using realistic embryonic tissue growth, thickness and stiffness values, recapitulate the transition from two to three ridges observed during embryonic development of the dragon's frill.

    view more details on Pubmed

  • Numerical approximation of a 3D mechanochemical interface model for skin patterning.

    Journal of Computational Physics 384 (2019) 383–404

    abstract

    We introduce a model for the mass transfer of molecular activators and inhibitors in two media separated by an interface, and study its interaction with the deformations exhibited by the two-layer skin tissue where they occur. The mathematical model results in a system of nonlinear advection-diffusion–reaction equations including cross-diffusion, and coupled with an interface elasticity problem. We propose a Galerkin method for the discretisation of the set of governing equations, involving also a suitable Newton linearisation, partitioned techniques, non-overlapping Schwarz alternating schemes, and high-order adaptive time stepping algorithms. The experimental accuracy and robustness of the proposed partitioned numerical methods is assessed, and some illustrating tests in 2D and 3D are provided to exemplify the coupling effects between the mechanical properties and the advection-diffusion–reaction interactions involving the two separate layers.

    see on external website

  • Phylogenetic mapping of scale nanostructure diversity in snakes.

    BMC Evol. Biol. 2019 04;19(1):91. 10.1186/s12862-019-1411-6. 10.1186/s12862-019-1411-6. PMC6469093.

    abstract

    Many species of snakes exhibit epidermal surface nanostructures that form complex motifs conferring self-cleaning properties, and sometimes structural iridescence, to their skin.

    view more details on Pubmed

  • Phylogenetic mapping of scale nanostructure diversity in snakes.

    BMC Evolutionary Biology2019 19:91 https://doi.org/10.1186/s12862-019-1411-6

    abstract

    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.

    see on external website

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

    PLoS Biol. 2019 Feb;17(2):e3000132. 10.1371/journal.pbio.3000132. PBIOLOGY-D-18-00761.

    abstract

    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.

    view more details on Pubmed

  • Locally-curved geometry generates bending cracks in the African elephant skin.

    Nat Commun 2018 Oct;9(1):3865. 10.1038/s41467-018-06257-3. 10.1038/s41467-018-06257-3. PMC6168576.

    abstract

    An intricate network of crevices adorns the skin surface of the African bush elephant, Loxodonta africana. These micrometre-wide channels enhance the effectiveness of thermal regulation (by water retention) as well as protection against parasites and intense solar radiation (by mud adherence). While the adaptive value of these structures is well established, their morphological characterisation and generative mechanism are unknown. Using microscopy, computed tomography and a custom physics-based lattice model, we show that African elephant skin channels are fractures of the animal brittle and desquamation-deficient skin outermost layer. We suggest that the progressive thickening of the hyperkeratinised stratum corneum causes its fracture due to local bending mechanical stress in the troughs of a lattice of skin millimetric elevations. The African elephant skin channels are therefore generated by thickening of a brittle material on a locally-curved substrate rather than by a canonical tensile cracking process caused by frustrated shrinkage.

    view more details on Pubmed

  • The anatomical placode in reptile scale morphogenesis indicates shared ancestry among skin appendages in amniotes.

    Sci Adv 2016 Jun;2(6):e1600708. 10.1126/sciadv.1600708. 1600708. PMC5392058.

    abstract

    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.

    view more details on Pubmed

  • A living mesoscopic cellular automaton made of skin scales.

    Nature 2017 Apr;544(7649):173-179. nature22031. 10.1038/nature22031.

    abstract

    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.

    view more details on Pubmed

  • Bifurcation Analysis of Reaction Diffusion Systems on Arbitrary Surfaces.

    Bull. Math. Biol. 2017 Feb;():. 10.1007/s11538-017-0255-8. 10.1007/s11538-017-0255-8.

    abstract

    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.

    view more details on Pubmed

  • A Step-by-step Guide to Assemble a Reptilian Genome

    Chapter in "Methods in Molecular Biology - Avian and Reptilian Developmental Biology".

    abstract

    This volume discusses recent advances in avian and reptilian biology that have caused this diverse field to re-emerge. The chapters in this book are divided into 4 parts: genomics and transcriptomics, genetic manipulation, stem cells, and new model systems. Part I details how to perform genomic and transcriptomic analyses in birds and reptiles; Part II highlights technological advancements in avian genetic manipulation; Part III focuses on methods to handle pluripotent cells; and part IV looks at the emerging models in avian and reptilian developmental biology. Written in the highly successful Methods in Molecular Biology series format, chapters include introductions to their respective topics, lists of the necessary materials and reagents, step-by-step, readily reproducible laboratory protocols, and tips on troubleshooting and avoiding known pitfalls. Cutting-edge and comprehensive, Avian and Reptilian Developmental Biology: Methods and Protocols explores a variety of approaches and different sauropsid models that will help facilitate communication and collaboration among researchers, which in turn will progress this field forward.

    see on external website

  • Reorganisation of Hoxd regulatory landscapes during the evolution of a snake-like body plan.

    Elife 2016 ;5():. 10.7554/eLife.16087. PMC4969037.

    abstract

    Within land vertebrate species, snakes display extreme variations in their body plan, characterized by the absence of limbs and an elongated morphology. Such a particular interpretation of the basic vertebrate body architecture has often been associated with changes in the function or regulation of Hox genes. Here, we use an interspecies comparative approach to investigate different regulatory aspects at the snake HoxD locus. We report that, unlike in other vertebrates, snake mesoderm-specific enhancers are mostly located within the HoxD cluster itself rather than outside. In addition, despite both the absence of limbs and an altered Hoxd gene regulation in external genitalia, the limb-associated bimodal HoxD chromatin structure is maintained at the snake locus. Finally, we show that snake and mouse orthologous enhancer sequences can display distinct expression specificities. These results show that vertebrate morphological evolution likely involved extensive reorganisation at Hox loci, yet within a generally conserved regulatory framework.

    view more details on Pubmed

  • The anatomical placode in reptile scale morphogenesis indicates shared ancestry among skin appendages in amniotes

    Science Advances 24 Jun 2016: Vol. 2, no. 6, e1600708; DOI: 10.1126/sciadv.1600708

    abstract

    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.

    see on external website

  • The SIB Swiss Institute of Bioinformatics' resources: focus on curated databases.

    Nucleic Acids Res. 2016 Jan 4;44(D1):D27-37. doi: 10.1093/nar/gkv1310. Epub 2015 Nov 28

    abstract

    The SIB Swiss Institute of Bioinformatics (www.isb-sib.ch) provides world-class bioinformatics databases, software tools, services and training to the international life science community in academia and industry. These solutions allow life scientists to turn the exponentially growing amount of data into knowledge. Here, we provide an overview of SIB's resources and competence areas, with a strong focus on curated databases and SIB's most popular and widely used resources. In particular, SIB's Bioinformatics resource portal ExPASy features over 150 resources, including UniProtKB/Swiss-Prot, ENZYME, PROSITE, neXtProt, STRING, UniCarbKB, SugarBindDB, SwissRegulon, EPD, arrayMap, Bgee, SWISS-MODEL Repository, OMA, OrthoDB and other databases, which are briefly described in this article.

    view more details on Pubmed

  • Amelanism in the corn snake is associated with the insertion of an LTR-retrotransposon in the OCA2 gene.

    Sci Rep 2015 ;5():17118. srep17118. 10.1038/srep17118. PMC4657000.

    abstract

    The corn snake (Pantherophis guttatus) is a new model species particularly appropriate for investigating the processes generating colours in reptiles because numerous colour and pattern mutants have been isolated in the last five decades. Using our captive-bred colony of corn snakes, transcriptomic and genomic next-generation sequencing, exome assembly, and genotyping of SNPs in multiple families, we delimit the genomic interval bearing the causal mutation of amelanism, the oldest colour variant observed in that species. Proceeding with sequencing the candidate gene OCA2 in the uncovered genomic interval, we identify that the insertion of an LTR-retrotransposon in its 11(th) intron results in a considerable truncation of the p protein and likely constitutes the causal mutation of amelanism in corn snakes. As amelanistic snakes exhibit white, instead of black, borders around an otherwise normal pattern of dorsal orange saddles and lateral blotches, our results indicate that melanocytes lacking melanin are able to participate to the normal patterning of other colours in the skin. In combination with research in the zebrafish, this work opens the perspective of using corn snake colour and pattern variants to investigate the generative processes of skin colour patterning shared among major vertebrate lineages.

    view more details on Pubmed

  • Reptilian Transcriptomes v2.0: An Extensive Resource for Sauropsida Genomics and Transcriptomics.

    Genome Biol Evol 2015 Jun;7(6):1827-41. evv106. 10.1093/gbe/evv106. PMC4494049.

    abstract

    Despite the availability of deep-sequencing techniques, genomic and transcriptomic data remain unevenly distributed across phylogenetic groups. For example, reptiles are poorly represented in sequence databases, hindering functional evolutionary and developmental studies in these lineages substantially more diverse than mammals. In addition, different studies use different assembly and annotation protocols, inhibiting meaningful comparisons. Here, we present the "Reptilian Transcriptomes Database 2.0," which provides extensive annotation of transcriptomes and genomes from species covering the major reptilian lineages. To this end, we sequenced normalized complementary DNA libraries of multiple adult tissues and various embryonic stages of the leopard gecko and the corn snake and gathered published reptilian sequence data sets from representatives of the four extant orders of reptiles: Squamata (snakes and lizards), the tuatara, crocodiles, and turtles. The LANE runner 2.0 software was implemented to annotate all assemblies within a single integrated pipeline. We show that this approach increases the annotation completeness of the assembled transcriptomes/genomes. We then built large concatenated protein alignments of single-copy genes and inferred phylogenetic trees that support the positions of turtles and the tuatara as sister groups of Archosauria and Squamata, respectively. The Reptilian Transcriptomes Database 2.0 resource will be updated to include selected new data sets as they become available, thus making it a reference for differential expression studies, comparative genomics and transcriptomics, linkage mapping, molecular ecology, and phylogenomic analyses involving reptiles. The database is available at www.reptilian-transcriptomes.org and can be enquired using a wwwblast server installed at the University of Geneva.

    view more details on Pubmed

  • R(2)OBBIE-3D, a Fast Robotic High-Resolution System for Quantitative Phenotyping of Surface Geometry and Colour-Texture.

    PLoS ONE 2015 ;10(6):e0126740. 10.1371/journal.pone.0126740. PONE-D-14-42512. PMC4454658.

    abstract

    While recent imaging techniques provide insights into biological processes from the molecular to the cellular scale, phenotypes at larger scales remain poorly amenable to quantitative analyses. For example, investigations of the biophysical mechanisms generating skin morphological complexity and diversity would greatly benefit from 3D geometry and colour-texture reconstructions. Here, we report on R(2)OBBIE-3D, an integrated system that combines a robotic arm, a high-resolution digital colour camera, an illumination basket of high-intensity light-emitting diodes and state-of-the-art 3D-reconstruction approaches. We demonstrate that R(2)OBBIE generates accurate 3D models of biological objects between 1 and 100 cm, makes multiview photometric stereo scanning possible in practical processing times, and enables the capture of colour-texture and geometric resolutions better than 15 μm without the use of magnifying lenses. R(2)OBBIE has the potential to greatly improve quantitative analyses of phenotypes in addition to providing multiple new applications in, e.g., biomedical science.

    view more details on Pubmed

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

    Mol. Ecol. 2015 Jul;24(13):3455-66. 10.1111/mec.13241. PMC4755148.

    abstract

    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.

    view more details on Pubmed

  • Bifurcation Analysis of Reaction Diffusion Systems on Arbitrary Surfaces

    arXiv:1605.01583 [cs.CE]

    abstract

    In this paper we present computational techniques to investigate the solutions of 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 characterize 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.

    see on external website

  • Photonic crystals cause active colour change in chameleons.

    Nat Commun 2015 ;6():6368. ncomms7368. 10.1038/ncomms7368. PMC4366488.

    abstract

    Many chameleons, and panther chameleons in particular, have the remarkable ability to exhibit complex and rapid colour changes during social interactions such as male contests or courtship. It is generally interpreted that these changes are due to dispersion/aggregation of pigment-containing organelles within dermal chromatophores. Here, combining microscopy, photometric videography and photonic band-gap modelling, we show that chameleons shift colour through active tuning of a lattice of guanine nanocrystals within a superficial thick layer of dermal iridophores. In addition, we show that a deeper population of iridophores with larger crystals reflects a substantial proportion of sunlight especially in the near-infrared range. The organization of iridophores into two superposed layers constitutes an evolutionary novelty for chameleons, which allows some species to combine efficient camouflage with spectacular display, while potentially providing passive thermal protection.

    view more details on Pubmed

  • Two waves of anisotropic growth generate enlarged follicles in the spiny mouse.

    Evodevo 2014 ;5():33. 10.1186/2041-9139-5-33. 2041-9139-5-33. PMC4335386.

    abstract

    Mammals exhibit a remarkable variety of phenotypes and comparative studies using novel model species are needed to uncover the evolutionary developmental mechanisms generating this diversity. Here, we undertake a developmental biology and numerical modeling approach to investigate the development of skin appendages in the spiny mouse, Acomys dimidiatus.

    view more details on Pubmed

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

    Proc. Natl. Acad. Sci. U.S.A. 2014 Oct;111(41):E4332-41. 1412268111. 10.1073/pnas.1412268111. PMC4205661.

    abstract

    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.

    view more details on Pubmed

  • Interactive Diffraction from Biological Nanostructures

    COMPUTER GRAPHICS forum Volume 33 (2014), number 8 pp. 177–188, DOI: 10.1111/cgf.12425

    abstract

    We describe a technique for interactive rendering of diffraction effects produced by biological nanostructures, such as snake skin surface gratings. Our approach uses imagery from atomic force microscopy that accurately captures the geometry of the nanostructures responsible for structural colouration, that is, colouration due to wave interference, in a variety of animals. We develop a rendering technique that constructs bidirectional reflection distribution functions (BRDFs) directly from the measured data and leverages pre-computation to achieve interactive performance. We demonstrate results of our approach using various shapes of the surface grating nanostructures. Finally, we evaluate the accuracy of our pre-computation-based technique and compare to a reference BRDF construction technique

    see on external website

  • The genome sequence of the corn snake (Pantherophis guttatus), a valuable resource for EvoDevo studies in squamates.

    Int. J. Dev. Biol. 2014 ;58(10-12):881-8. 150060at. 10.1387/ijdb.150060at.

    abstract

    Squamates (snakes and lizards) exhibit a striking variety of phenotypes, with little known on their generative mechanisms. Studies aiming to understand the genetic basis of this wide diversity in morphology, physiology and ecology will greatly benefit from whole genome sequencing initiatives, as they provide the foundation for comparative analyses and improve our understanding of the evolution, development and diversification of traits. Here, we present the first draft genome of the corn snake Pantherophis guttatus, an oviparous snake that we promote as a particularly appropriate model species for evolutionary developmental studies in squamates. We sequenced 100-base paired-end reads from multiple individuals of a single family (parents and offspring) that produced a genome assembly of 1.53 gigabases (Gb), roughly covering 75% of the expected total genome size, and 297,768 scaffolds >1 Kb. We were able to fully retrieve 86, and partially another 106, of the 248 CEGMA core genes, indicating that a high genome completeness was achieved, even though the assembly is fragmented. Using MAKER2, we annotated 10,917 genes with high confidence (Annotation Edit Distance (AED)<1) and an additional 5,263 predicted genes matched with the species' transcriptome. Numerous colour and colour pattern morphs exist in P. guttatus, making it an ideal model to study the genetic determinism, development, and evolution of adaptive colour traits in reptiles. Using our draft genome and a Single-Nucleotide Polymorphism (SNP) calling approach, we confirmed the interval with the causative mutation for the amelanistic phenotype, a result supported by a parallel exome-based study.

    view more details on Pubmed

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

    BMC Biol. 2013 ;11():105. 1741-7007-11-105. 10.1186/1741-7007-11-105. PMC4021644.

    abstract

    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.

    view more details on Pubmed

  • Multiple enhancers regulate Hoxd genes and the Hotdog LncRNA during cecum budding.

    Cell Rep 2013 Oct;5(1):137-50. S2211-1247(13)00500-7. 10.1016/j.celrep.2013.09.002.

    abstract

    Hox genes are required for the development of the intestinal cecum, a major organ of plant-eating species. We have analyzed the transcriptional regulation of Hoxd genes in cecal buds and show that they are controlled by a series of enhancers located in a gene desert flanking the HoxD cluster. The start site of two opposite long noncoding RNAs (lncRNAs), Hotdog and Twin of Hotdog, selectively contacts the expressed Hoxd genes in the framework of a topological domain, coinciding with robust transcription of these genes during cecum budding. Both lncRNAs are specifically transcribed in the cecum, albeit bearing no detectable function in trans. Hedgehogs have kept this regulatory potential despite the absence of the cecum, suggesting that these mechanisms are used in other developmental situations. In this context, we discuss the implementation of a common "budding toolkit" between the cecum and the limbs.

    view more details on Pubmed

  • Crocodylians evolved scattered multi-sensory micro-organs.

    Evodevo 2013 ;4(1):19. 2041-9139-4-19. 10.1186/2041-9139-4-19. PMC3711810.

    abstract

    During their evolution towards a complete life cycle on land, stem reptiles developed both an impermeable multi-layered keratinized epidermis and skin appendages (scales) providing mechanical, thermal, and chemical protection. Previous studies have demonstrated that, despite the presence of a particularly armored skin, crocodylians have exquisite mechanosensory abilities thanks to the presence of small integumentary sensory organs (ISOs) distributed on postcranial and/or cranial scales.

    view more details on Pubmed

  • Patterns of ossification in southern versus northern placental mammals.

    Evolution 2013 Jul;67(7):1994-2010. 10.1111/evo.12071.

    abstract

    Consensus on placental mammal phylogeny is fairly recent compared to that for vertebrates as a whole. A stable phylogenetic hypothesis enables investigation into the possibility that placental clades differ from one another in terms of their development. Here, we focus on the sequence of skeletal ossification as a possible source of developmental distinctiveness in "northern" (Laurasiatheria and Euarchontoglires) versus "southern" (Afrotheria and Xenarthra) placental clades. We contribute data on cranial and postcranial ossification events during growth in Afrotheria, including elephants, hyraxes, golden moles, tenrecs, sengis, and aardvarks. We use three different techniques to quantify sequence heterochrony: continuous method, sequence-ANOVA (analysis of variance) and event-paring/Parsimov. We show that afrotherians significantly differ from other placentals by an early ossification of the orbitosphenoid and caudal vertebrae. Our analysis also suggests that both southern placental groups show a greater degree of developmental variability; however, they rarely seem to vary in the same direction, especially regarding the shifts that differ statistically. The latter observation is inconsistent with the Atlantogenata hypothesis in which afrotherians are considered as the sister clade of xenarthrans. Interestingly, ancestral nodes for Laurasiatheria and Euarchontoglires show very similar trends and our results suggest that developmental homogeneity in some ossification sequences may be restricted to northern placental mammals (Boreoeutheria).

    view more details on Pubmed

  • Recovery of a nearly extinct Galápagos tortoise despite minimal genetic variation.

    Evol Appl 2013 Feb;6(2):377-83. 10.1111/eva.12014. PMC3586625.

    abstract

    A species of Galápagos tortoise endemic to Española Island was reduced to just 12 females and three males that have been bred in captivity since 1971 and have produced over 1700 offspring now repatriated to the island. Our molecular genetic analyses of juveniles repatriated to and surviving on the island indicate that none of the tortoises sampled in 1994 had hatched on the island versus 3% in 2004 and 24% in 2007, which demonstrates substantial and increasing reproduction in situ once again. This recovery occurred despite the parental population having an estimated effective population size <8 due to a combination of unequal reproductive success of the breeders and nonrandom mating in captivity. These results provide guidelines for adapting breeding regimes in the parental captive population and decreasing inbreeding in the repatriated population. Using simple morphological data scored on the sampled animals, we also show that a strongly heterogeneous distribution of tortoise sizes on Española Island observed today is due to a large variance in the number of animals included in yearly repatriation events performed in the last 40 years. Our study reveals that, at least in the short run, some endangered species can recover dramatically despite a lack of genetic variation and irregular repatriation efforts.

    view more details on Pubmed

  • Contrasted evolution of the vomeronasal receptor repertoires in mammals and squamate reptiles.

    Genome Biol Evol 2013 ;5(2):389-401. evt013. 10.1093/gbe/evt013. PMC3590772.

    abstract

    The vomeronasal organ (VNO) is an olfactory structure that detects pheromones and environmental cues. It consists of sensory neurons that express evolutionary unrelated groups of transmembrane chemoreceptors. The predominant V1R and V2R receptor repertoires are believed to detect airborne and water-soluble molecules, respectively. It has been suggested that the shift in habitat of early tetrapods from water to land is reflected by an increase in the ratio of V1R/V2R genes. Snakes, which have a very large VNO associated with a sophisticated tongue delivery system, are missing from this analysis. Here, we use RNA-seq and RNA in situ hybridization to study the diversity, evolution, and expression pattern of the corn snake vomeronasal receptor repertoires. Our analyses indicate that snakes and lizards retain an extremely limited number of V1R genes but exhibit a large number of V2R genes, including multiple lineages of reptile-specific and snake-specific expansions. We finally show that the peculiar bigenic pattern of V2R vomeronasal receptor gene transcription observed in mammals is conserved in squamate reptiles, hinting at an important but unknown functional role played by this expression strategy. Our results do not support the hypothesis that the shift to a vomeronasal receptor repertoire dominated by V1Rs in mammals reflects the evolutionary transition of early tetrapods from water to land. This study sheds light on the evolutionary dynamics of the vomeronasal receptor families in vertebrates and reveals how mammals and squamates differentially adapted the same ancestral vomeronasal repertoire to succeed in a terrestrial environment.

    view more details on Pubmed

  • Crocodile head scales are not developmental units but emerge from physical cracking.

    Science 2013 Jan;339(6115):78-81. science.1226265. 10.1126/science.1226265.

    abstract

    Various lineages of amniotes display keratinized skin appendages (feathers, hairs, and scales) that differentiate in the embryo from genetically controlled developmental units whose spatial organization is patterned by reaction-diffusion mechanisms (RDMs). We show that, contrary to skin appendages in other amniotes (as well as body scales in crocodiles), face and jaws scales of crocodiles are random polygonal domains of highly keratinized skin, rather than genetically controlled elements, and emerge from a physical self-organizing stochastic process distinct from RDMs: cracking of the developing skin in a stress field. We suggest that the rapid growth of the crocodile embryonic facial and jaw skeleton, combined with the development of a very keratinized skin, generates the mechanical stress that causes cracking.

    view more details on Pubmed

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

    J. Anat. 2013 Jan;222(1):2-18. 10.1111/j.1469-7580.2012.01509.x. PMC3552411.

    abstract

    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.

    view more details on Pubmed

  • Reptilian-transcriptome v1.0, a glimpse in the brain transcriptome of five divergent Sauropsida lineages and the phylogenetic position of turtles.

    Evodevo 2011 ;2(1):19. 2041-9139-2-19. 10.1186/2041-9139-2-19. PMC3192992.

    abstract

    Reptiles are largely under-represented in comparative genomics despite the fact that they are substantially more diverse in many respects than mammals. Given the high divergence of reptiles from classical model species, next-generation sequencing of their transcriptomes is an approach of choice for gene identification and annotation.

    view more details on Pubmed

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

    Mol. Ecol. 2011 Aug;20(15):3237-50. 10.1111/j.1365-294X.2011.05158.x.

    abstract

    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.

    view more details on Pubmed

  • Analysis of global and local population stratification of Finless Porpoises Neophocaena phocaenoides in Chinese waters

    Marine Biology, 158(8):1791-1804.

    abstract

    The existence of three distinct populations is widely accepted for the finless porpoise (Neophocaena phocaenoides) in Chinese waters: the Yellow Sea, Yangtze River, and South China Sea populations. Here, we use nine species-specific microsatellite loci, the complete mitochondrial DNA control region (912 bp), and the complete mitochondrial cytochrome b gene (1,140 bp) to further investigate potential population stratification in the Yellow Sea using 147 finless porpoise samples from the Bohai Sea and adjacent northern Yellow Sea, two regions that were largely underrepresented in previous genetic studies. Our F-statistics analyses confirm the previously described three populations, but further demonstrate significant genetic differentiation between the [Bohai + northern Yellow] Sea and the southern Yellow Sea. On the other hand, median-joining network analyses do not exhibit well-differentiated haplotype groups among different geographic populations, suggesting the existence of shared ancestral haplotypes. Levels of microsatellite diversity are moderate to high (mean HE = 0.794) among the 147 [Bohai + northern Yellow] Sea finless porpoises and no recent bottleneck was detected, whereas mtDNA control region and cytochrome b gene diversity is low to moderate. The microsatellite genotypic and mtDNA haplotypic data also confirm the presence of mother-calf pairs in single-net bycatch cases. The results presented here highlight the necessity to treat the [Bohai + northern Yellow] Sea population (highly impacted by anthropogenic threats) as a separate Management Unit.

    see on external website

  • Development of a multiplex PCR assay for fine-scale population genetic analysis of the Komodo monitor Varanus komodoensis based on 18 polymorphic microsatellite loci.

    Mol Ecol Resour 2011 May;11(3):550-6. 10.1111/j.1755-0998.2011.02986.x.

    abstract

    Multiplex PCR assays for the coamplification of microsatellite loci allow rapid and cost-effective genetic analyses and the production of efficient screening protocols for international breeding programs. We constructed a partial genomic library enriched for di-nucleotide repeats and characterized 14 new microsatellite loci for the Komodo monitor (or Komodo dragon, Varanus komodoensis). Using these novel microsatellites and four previously described loci, we developed multiplex PCR assays that may be loaded on a genetic analyser in three separate panels. We tested the novel set of microsatellites for polymorphism using 69 individuals from three island populations and evaluated the resolving power of the entire panel of 18 loci by conducting (i) a preliminary assignment test to determine population(s) of origin and (ii) a parentage analysis for 43 captive Komodo monitors. This panel of polymorphic loci proved useful for both purposes and thus can be exploited for fine-scale population genetic analyses and as part of international captive breeding programs directed at maintaining genetically viable ex situ populations and reintroductions.

    view more details on Pubmed

  • Considerations for the inclusion of 2x mammalian genomes in phylogenetic analyses (a response to Vilella et al.’s correspondence)

    Genome Biology, 12: 401.

    abstract

    a response to Vilella et al.’s correspondence (A response to 2x genomes - depth does matter by MC Milinkovitch, R Helaers, E Depiereux, AC Tzika and T Gabaldón. Genome Biol 2010, 11:R16.)

    see on external website

  • Molecular analyses of small cetacean samples from Peruvian fish markets

    Conservation Genetics, 11(6):2207-2218.

    abstract

    In the last 60 years, incidental entanglement in fishing gears (so called by-catch) became the main cause of mortality worldwide for small cetaceans and is pushing several populations and species to the verge of extinction. Thus, monitoring and quantifying by-catches is an important step towards proper and sustainable management of cetacean populations. Continuous studies indicated that by-catches and directed takes of small cetaceans in Peru greatly increased since 1985. Legal measures banning cetacean takes, enforced in 1994 and 1996, ironically made monitoring highly problematic as fishers continue catching these animals but utilize or dispose of carcasses clandes- tinely. Hence, in locations where cetaceans are landed covertly or already butchered, molecular genetic methods can provide the only means of identification of the species, sex, and sometimes the population of each sample. Here, we generate and analyse a fragment of the mitochondrial DNA cytochrome b gene and 5 nuclear microsatellite markers from 182 meat and skin samples of unidentified small cetaceans collected at three Peruvian markets between July 2006 and April 2007. Our results, compared to past surveys, indicate that Lagenorhynchus obscurus, Phocoena spinipinnis, Tursiops truncatus, Delphinus capensis, and D. delphis continue to be caught and marketed, but that the relative incidence of P. spinipinnis is highly reduced, pos- sibly because of population depletion. The small number of possible sampling duplicates demonstrates that a high monitoring frequency is required for a thorough evaluation of incidental catches in the area. A wide public debate on by-catch mitigation measures is greatly warranted in Peru.

    see on external website

  • MetaPIGA v2.0: maximum likelihood large phylogeny estimation using the metapopulation genetic algorithm and other stochastic heuristics.

    BMC Bioinformatics 2010 ;11():379. 1471-2105-11-379. 10.1186/1471-2105-11-379. PMC2912891.

    abstract

    The development, in the last decade, of stochastic heuristics implemented in robust application softwares has made large phylogeny inference a key step in most comparative studies involving molecular sequences. Still, the choice of a phylogeny inference software is often dictated by a combination of parameters not related to the raw performance of the implemented algorithm(s) but rather by practical issues such as ergonomics and/or the availability of specific functionalities.

    view more details on Pubmed

  • Historical constraints on vertebrate genome evolution.

    Genome Biol Evol 2010 ;2():13-8. 10.1093/gbe/evp052. PMC2839353.

    abstract

    Recent analyses indicated that genes with larger effect of knockout or mutation and with larger probability to revert to single copy after whole genome duplication are expressed earlier in development. Here, we further investigate whether tissue specificity of gene expression is constrained by the age of origin of the corresponding genes. We use 38 metazoan genomes and a comparative genomic application system to integrate inference of gene duplication with expression data from 17,503 human genes into a strictly phylogenetic framework. We show that the number of anatomical systems in which genes are expressed decreases steadily with decreased age of the genes' first appearance in the phylogeny: the oldest genes are expressed, on average, in twice as many anatomical systems than the genes gained recently in evolution. These results are robust to different sources of expression data, to different levels of the anatomical system hierarchy, and to the use of gene families rather than duplication events. Finally, we show that the rate of increase in gene tissue specificity correlates with the relative rate of increase in the maximum number of cell types in the corresponding taxa. Although subfunctionalization and increase in cell type number throughout evolution could constitute, respectively, the proximal and ultimate causes of this correlation, the two phenomena are intermingled. Our analyses identify a striking historical constraint in gene expression: the number of cell types in existence at the time of a gene appearance (through duplication or de novo origination) tends to determine its level of tissue specificity for tens or hundreds of millions of years.

    view more details on Pubmed

  • Preliminary Genetic Status of the Spotted Seal Phoca largha in Liaodong Bay Based on Microsatellite and Mitochondrial DNA Analyses

    Trends in Evolutionary Biology, 2(6):33-38.

    abstract

    The Liaodong Bay spotted seal (Phoca largha) population experienced several drastic declines in the last 80 years. Recent studies are contradictory regarding the level of genetic diversity and population structure of P. largha, possibly because of the use of non-species-specific nuclear markers. Here, we report on i) the first isolation and characterization of 10 species-specific polymorphic microsatellite loci for the spotted seal, ii) sequences of a 572 bp mtDNA fragment in 25 Liaodong Bay individuals that we analyzed together with all published haplotypes from Liaodong Bay and Japan. Intermediate genetic diversity in microsatellite loci was found in the Liaodong Bay population and the effective population size estimates were 41.8 to 86.8 individuals. Low mtDNA genetic variability, especially nucleotide diversity, in the Liaodong Bay population was detected, but Bayesian skyline plots did not show any evidence of recent bottleneck. Both F-statistics and the haplotypic network indicate a clear differentiation between the Liaodong Bay and Japanese populations separated by a fixed mutation. Analysis of mtDNA data indicates that Liaodong Bay female seals show fidelity to their breeding site, and breeding time data suggest that this population is reproductively isolated from populations in other breeding areas. The observed low genetic diversity in mtDNA and the intermediate levels of nuclear microsatellite diversity, combined with the potential genetic isolation, suggest that the Liaodong Bay population might be at risk and that further investigation of the population genetics of spotted seals across their whole range is warranted for proper management of the species.

    see on external website

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

    Nature 2010 Mar;464(7285):99-103. nature08789. 10.1038/nature08789.

    abstract

    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.

    view more details on Pubmed

  • 2x genomes — depth does matter.

    Genome Biol. 2010 ;11(2):R16. gb-2010-11-2-r16. 10.1186/gb-2010-11-2-r16. PMC2872876.

    abstract

    Given the availability of full genome sequences, mapping gene gains, duplications, and losses during evolution should theoretically be straightforward. However, this endeavor suffers from overemphasis on detecting conserved genome features, which in turn has led to sequencing multiple eutherian genomes with low coverage rather than fewer genomes with high-coverage and more even distribution in the phylogeny. Although limitations associated with analysis of low coverage genomes are recognized, they have not been quantified.

    view more details on Pubmed

  • Phylogenomics of unusual histone H2A Variants in Bdelloid rotifers.

    PLoS Genet. 2009 Mar;5(3):e1000401. 10.1371/journal.pgen.1000401. PMC2642717.

    abstract

    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.

    view more details on Pubmed

  • Molecular genetic analysis of a captive-breeding program: the vulnerable endemic Jamaican yellow boa

    Conservation Genetics, 10 (1):69-77.

    abstract

    The endemic Jamaican boa (or “yellow boa”, Epicrates subflavus) is a vulnerable species of the Caribbean biodiversity hotspot whose natural populations greatly declined mainly due to predation by introduced species, human persecution, and habitat destruction. A captive breeding program was initiated in 1976 and rationalized in 2002 by the establishment of a European Endangered Species Program. During the last 30 years, more than 600 offspring, of which 80 are still alive today, have been produced and distributed among European host institutions and privates. Here, using nine nuclear microsatellite loci and a fragment of the mitochondrial cytochrome b gene, we (i) determine the natural population from which the founders originate, (ii) identify parental allocation errors and ambiguities in the studbook, and (iii) assess the genetic diversity and estimate levels of inbreeding of the current captive population based on loss of alleles, variance in reproductive success, and relatedness among individuals. Combining measures of relatedness derived from multilocus genotypes with practical parameters such as age of animals and localization of host institutions, we propose mating groups that would maximize genetic diversity in the captive population of the Jamaican boa. Our analyses provide guidance for a more efficient breeding program that, in turn, could be used as the starting point of a repatriation program to increase the probability of the species long-term survival.

    see on external website

  • Population genetics of Galápagos land iguana (genus Conolophus) remnant populations.

    Mol. Ecol. 2008 Dec;17(23):4943-52. MEC3967. 10.1111/j.1365-294X.2008.03967.x.

    abstract

    The Galápagos land iguanas (genus Conolophus) have faced significant anthropogenic disturbances since the 17th century, leading to severe reduction of some populations and the extinction of others. Conservation activities, including the repatriation of captive-bred animals to depleted areas, have been ongoing since the late 1970s, but genetic information has not been extensively incorporated. Here we use nine species-specific microsatellite loci of 703 land iguanas from the six islands where the species occur today to characterize the genetic diversity within, and the levels of genetic differentiation among, current populations as well as test previous hypotheses about accidental translocations associated with early conservation efforts. Our analyses indicate that (i) five populations of iguanas represent distinct conservation units (one of them being the recently discovered rosada form) and could warrant species status, (ii) some individuals from North Seymour previously assumed to be from the natural Baltra population appear related to both Isabela and Santa Cruz populations, and (iii) the five different management units exhibit considerably different levels of intrapopulation genetic diversity, with the Plaza Sur and Santa Fe populations particularly low. Although the initial captive breeding programmes, coupled with intensive efforts to eradicate introduced species, saved several land iguana populations from extinction, our molecular results provide objective data for improving continuing in situ species survival plans and population management for this spectacular and emblematic reptile.

    view more details on Pubmed

  • Development of 10 highly-polymorphic microsatellite markers in the vulnerable Galápagos land iguanas (genus Conolophus).

    Mol Ecol Resour 2009 Jan;9(1):376-9. 10.1111/j.1755-0998.2008.02457.x.

    abstract

    The two species of Galápagos land iguanas (Conolophus subcristatus and C. pallidus) are listed as 'vulnerable' species by the International Union for the Conservation of Nature (IUCN Red List; http://www.iucnredlist.org). Here, we report on the isolation and characterization of 10 microsatellite markers using 562 individuals sampled on all Galápagos islands where Conolophus species occur today. We show that these 10 loci are highly polymorphic and display diagnostic alleles for five out of the six island populations. These markers will be useful for Conolophus population genetic analyses as well as for guiding ongoing captive breeding programmes.

    view more details on Pubmed

  • Characterization of microsatellite loci in the European pond turtle Emys orbicularis.

    Mol Ecol Resour 2009 Jan;9(1):189-91. 10.1111/j.1755-0998.2008.02205.x.

    abstract

    A set of eight highly polymorphic microsatellite markers was isolated and characterized from a genomic library enriched for dinucleotide repeats in the European pond turtle, Emys orbicularis. The markers were tested for polymorphism in a total of 33 turtles sampled in two natural ponds in the nature reserve of Kerkini, northern Greece. Number of alleles varied from 10 to 18, and expected heterozygosity ranged between 0.738 and 0.921. This novel set of loci will be particularly useful to assess fine-scale population structure and for parentage analysis in E. orbicularis.

    view more details on Pubmed

  • Mapping gene gains and losses among metazoan full genomes using an integrated phylogenetic framework

    Chapter 9 in “Evolution after Gene Duplication” Editors Dittmar K. and Liberles D. A., Wiley-Blackwell.

    abstract

    see on external website

  • Controlling population evolution in the laboratory to evaluate methods of historical inference.

    PLoS ONE 2008 ;3(8):e2960. 10.1371/journal.pone.0002960. PMC2491900.

    abstract

    Natural populations of known detailed past demographic history are extremely valuable to evaluate methods of historical inference, yet are extremely rare. As an alternative approach, we have generated multiple replicate microsatellite data sets from laboratory-cultured populations of a gonochoric free-living nematode, Caenorhabditis remanei, that were constrained to pre-defined demographic histories featuring different levels of migration among populations or bottleneck events of different magnitudes. These data sets were then used to evaluate the performances of two recently developed population genetics methods, BayesAss+, that estimates recent migration rates among populations, and Bottleneck, that detects the occurrence of recent bottlenecks. Migration rates inferred by BayesAss+ were generally over-estimates, although these were often included within the confidence interval. Analyses of data sets simulated in-silico, using a model mimicking the laboratory experiments, produced less biased estimates of the migration rates, and showed increased efficiency of the program when the number of loci and sampled genotypes per population was higher. In the replicates for which the pre-bottleneck laboratory-cultured populations did not significantly depart from a mutation/drift equilibrium, an important assumption of the program Bottleneck, only a portion of the bottleneck events were detected. This result was confirmed by in-silico simulations mirroring the laboratory bottleneck experiments. More generally, our study demonstrates the feasibility, and highlights some of the limits, of the approach that consists in generating molecular genetic data sets by controlling the evolution of laboratory-reared nematode populations, for the purpose of validating methods inferring population history.

    view more details on Pubmed

  • Automated discovery and phylogenetic analysis of new toxin-antitoxin systems.

    BMC Microbiol. 2008 ;8():104. 1471-2180-8-104. 10.1186/1471-2180-8-104. PMC2446400.

    abstract

    Although often viewed as elements "at the service of" bacteria, plasmids exhibit replication and maintenance mechanisms that make them purely "selfish DNA" candidates. Toxin-antitoxin (TA) systems are a spectacular example of such mechanisms: a gene coding for a cytotoxic stable protein is preceded by a gene coding for an unstable antitoxin. The toxin being more stable than the antitoxin, absence of the operon causes a reduction of the amount of the latter relative to the amount of the former. Thus, a cell exhibiting a TA system on a plasmid is 'condemned' either not to loose it or to die.

    view more details on Pubmed

  • Population structure of an endemic vulnerable species, the Jamaican boa (Epicrates subflavus).

    Mol. Ecol. 2008 Jan;17(2):533-44. MEC3588. 10.1111/j.1365-294X.2007.03588.x.

    abstract

    The Jamaican boa (Epicrates subflavus; also called Yellow boa) is an endemic species whose natural populations greatly and constantly declined since the late 19th century, mainly because of predation by introduced species, human persecution, and habitat destruction. In-situ conservation of the Jamaican boa is seriously hindered by the lack of information on demographic and ecological parameters as well as by a poor understanding of the population structure and species distribution in the wild. Here, using nine nuclear microsatellite loci and a fragment of the mitochondrial cytochrome b gene from 87 wild-born individuals, we present the first molecular genetic analyses focusing on the diversity and structure of the natural populations of the Jamaican boa. A model-based clustering analysis of multilocus microsatellite genotypes identifies three groups that are also significantly differentiated on the basis of F-statistics. Similarly, haplotypic network reconstruction methods applied on the cytochrome b haplotypes isolated here identify two well-differentiated haplogroups separated by four to six fixed mutations. Bayesian and metaGA analyses of the mitochondrial data set combined with sequences from other Boidae species indicate that rooting of the haplotypic network occurs most likely between the two defined haplogroups. Both analyses (based on nuclear and mitochondrial markers) underline an Eastern vs. (Western + Central) pattern of differentiation in agreement with geological data and patterns of differentiation uncovered in other vertebrate and invertebrate Jamaican species. Our results provide important insights for improving management of ex-situ captive populations and for guiding the development of proper in-situ species survival and habitat management plans for this spectacular, yet poorly known and vulnerable, snake.

    view more details on Pubmed

  • A Pragmatic Approach for Selecting Evo-Devo Model Species in Amniotes

    Chapter 7 in “Evolving Pathways: Key Themes in Evolutionary Developmental Biology.

    abstract

    One major classical justification of using a model metazoan species for experimentation has been that discoveries of biological phenomena in that species could be extrapolated to other multicellular species. Because the chances that this extrapolation is valid in humans depend on the phylogenetic distance between humans and the model species, many researchers have somewhat sacrificed the major benefits of small size, short generation time and ease of manipulation that characterise some invertebrates in order to use species that humans can more readily relate to, such as the laboratory mouse (Mus musculus). However, the community of biologists has continued to use additional model species because each of the selected taxa have specific features that make experimental manipulation easier (e.g. easy-to-score morphological variation and giant polytene chromosomes in Drosophila melanogaster, or accurate description of the largely invariant complete cell lineage and full neural connectivity in the roundworm Caenorhabditis elegans).

    see on external website

  • MANTIS: a phylogenetic framework for multi-species genome comparisons.

    Bioinformatics 2008 Jan;24(2):151-7. btm567. 10.1093/bioinformatics/btm567.

    abstract

    Practitioners of comparative genomics face huge analytical challenges as whole genome sequences and functional/expression data accumulate. Furthermore, the field would greatly benefit from a better integration of this wealth of data with evolutionary concepts.

    view more details on Pubmed

  • An ant colony optimization algorithm for phylogenetic estimation under the minimum evolution principle.

    BMC Evol. Biol. 2007 ;7():228. 1471-2148-7-228. 10.1186/1471-2148-7-228. PMC2211314.

    abstract

    Distance matrix methods constitute a major family of phylogenetic estimation methods, and the minimum evolution (ME) principle (aiming at recovering the phylogeny with shortest length) is one of the most commonly used optimality criteria for estimating phylogenetic trees. The major difficulty for its application is that the number of possible phylogenies grows exponentially with the number of taxa analyzed and the minimum evolution principle is known to belong to the NRho- hard class of problems.

    view more details on Pubmed

  • Conceptual bases for quantifying the role of the environment on gene evolution: the participation of positive selection and neutral evolution.

    Biol Rev Camb Philos Soc 2007 Nov;82(4):551-72. BRV24. 10.1111/j.1469-185X.2007.00024.x.

    abstract

    To demonstrate that a given change in the environment has contributed to the emergence of a given genotypic and phenotypic shift during the course of evolution, one should ask to what extent such shifts would have occurred without environmental change. Of course, such tests are rarely practical but phenotypic novelties can still be correlated to genomic shifts in response to environmental changes if enough information is available. We surveyed and re-evaluated the published data in order to estimate the role of environmental changes on the course of species and genomic evolution. Only a few published examples clearly demonstrate a causal link between a given environmental change and the fixation of a genomic variant resulting in functional modification (gain, loss or alteration of function). Many others suggested a link between a given phenotypic shift and a given environmental change but failed to identify the underlying genomic determinant(s) and/or the associated functional consequence(s). The proportion of genotypic and phenotypic variation that is fixed concomitantly with environmental changes is often considered adaptive and hence, the result of positive selection, even though alternative causes, such as genetic drift, are rarely investigated. Therefore, the second aim herein is to review evidence for the mechanisms leading to fixation.

    view more details on Pubmed

  • Escaping the mouse trap: the selection of new Evo-Devo model species.

    J. Exp. Zool. B Mol. Dev. Evol. 2007 Jul;308(4):337-46. 10.1002/jez.b.21180.

    abstract

    Among the many, sometimes contradictory, criteria that have been used for promoting model species, the most prominent has probably been their relevance for understanding human biology. Recently however, the debate has partly shifted from the search for evolutionary conservation (medicine-driven models) to a better understanding of the generative mechanisms underlying biological diversity (Evo-Devo-driven models). Integration of multiple disciplines, beyond developmental genetics and evolutionary molecular genetics, as well as of innovative technologies will help biologists to open the massive realm of living species to genome manipulation and phenotypic investigation. However, a consensual list of model species must still be reached for optimizing the interplay between in silico analyses and in vivo experiments, and we claim that the Evo-Devo community should play a more energetic role in this endeavor. We discuss here a few criteria and limitations of major relevance to the choice of model species for Evo-Devo studies, and promote the use of a pragmatic approach. Finally, given the difficulties related to manipulating and breeding model species, we suggest the development of Evo-Devo virtual zoos maintaining breeding colonies of a selected set of species and from which eggs or staged embryos are available on order.

    view more details on Pubmed

  • Mitochondrial genome and nuclear sequence data support myzostomida as part of the annelid radiation.

    Mol. Biol. Evol. 2007 Aug;24(8):1690-701. msm086. 10.1093/molbev/msm086.

    abstract

    The echinoderm symbionts Myzostomida are marine worms that show an enigmatic lophotrochozoan body plan. Historically, their phylogenetic origins were obscured due to disagreement about which morphological features are evolutionarily conserved, but now most morphological evidence points to annelid origins. In contrast, recent phylogenetic analyses using different molecular markers produced variable results regarding the position of myzostomids, but all suggested these worms are not derived annelids. To reexamine this issue, we analyzed data from nuclear genes (18S rDNA, 28S rDNA, Myosin II, and Elongation Factor-1alpha), and a nearly complete myzostomid mitochondrial genome. Here, we show that the molecular data are in agreement with the morphological evidence that myzostomids are part of the annelid radiation. This result is robustly supported by mitochondrial (gene order and sequence data) and nuclear data, as well as by recent ultrastructural investigations. Using Bayes factor comparison, alternative hypotheses are shown to lack support. Thus, myzostomids probably evolved from a segmented ancestor and gained a derived anatomy during their long evolutionary history as echinoderm symbionts.

    view more details on Pubmed

  • Giant Galápagos tortoises; molecular genetic analyses identify a trans-island hybrid in a repatriation program of an endangered taxon.

    BMC Ecol. 2007 ;7():2. 1472-6785-7-2. 10.1186/1472-6785-7-2. PMC1820773.

    abstract

    Giant Galápagos tortoises on the island of Española have been the focus of an intensive captive breeding-repatriation programme for over 35 years that saved the taxon from extinction. However, analysis of 118 samples from released individuals indicated that the bias sex ratio and large variance in reproductive success among the 15 breeders has severely reduced the effective population size (Ne).

    view more details on Pubmed

  • Phylogeny and biogeography of a cosmopolitan frog radiation: Late cretaceous diversification resulted in continent-scale endemism in the family ranidae.

    Syst. Biol. 2006 Aug;55(4):579-94. R5051H424T28H706. 10.1080/10635150600812551.

    abstract

    Ranidae is a large anuran group with a nearly cosmopolitan distribution. We investigated the phylogenetic relationships and early biogeographic history of ranid frogs, using 104 representatives of all subfamilies and families, sampled from throughout their distribution. Analyses of approximately 1570 bp of nuclear gene fragments (Rag-1, rhod, Tyr) and approximately 2100 bp of the mitochondrial genome (12S rRNA, tRNAVAL, 16S rRNA) indicate that the monophyly of several taxa can be rejected with high confidence. Our tree is characterized by a clear historical association of each major clade with one Gondwanan plate. This prevalence of continent-scale endemism suggests that plate tectonics has played a major role in the distribution of ranid frogs. We performed dispersal-vicariance analyses, as well as analyses constrained by paleogeographic data, to estimate ancestral distributions during early ranid diversification. Additionally, we used molecular clock analyses to evaluate whether these scenarios fit the temporal framework of continental breakup. Our analyses suggest that a scenario in which the ancestors of several clades (Rhacophorinae, Dicroglossinae, Raninae) reached Eurasia via the Indian subcontinent, and the ancestor of Ceratobatrachinae entered via the Australia-New Guinea plate, best fits the paleogeographic models and requires the fewest number of dispersal/vicariance events. However, several alternatives, in which part of the ranid fauna colonized Laurasia from Africa, are not significantly worse. Most importantly, all hypotheses make clear predictions as to where to expect key fossils and where to sample other living ranids, and thus constitute a strong basis for further research.

    view more details on Pubmed

  • Molecular phylogenetic analyses indicate multiple independent emergences of parasitism in Myzostomida (Protostomia).

    Syst. Biol. 2006 Apr;55(2):208-27. T127351LPM06186Q. 10.1080/10635150500481317.

    abstract

    The fossil record indicates that Myzostomida, an enigmatic group of marine worms, traditionally considered as annelids, have exhibited a symbiotic relationship with echinoderms, especially crinoids, for nearly 350 million years. All known extant myzostomids are associated with echinoderms and infest their integument, gonads, celom, or digestive system. Using nuclear (18S rDNA) and mitochondrial (16S and COI) DNA sequence data from 37 myzostomid species representing nine genera, we report here the first molecular phylogeny of the Myzostomida and investigate the evolution of their various symbiotic associations. Our analyses indicate that the two orders Proboscidea and Pharyngidea do not constitute natural groupings. Character reconstruction analyses strongly suggest that (1) the ancestor of all extant myzostomids was an ectocommensal that first infested crinoids, and then asteroids and ophiuroids, and (2) parasitism in myzostomids emerged multiple times independently.

    view more details on Pubmed

  • A non-linear optimization procedure to estimate distances and instantaneous substitution rate matrices under the GTR model.

    Bioinformatics 2006 Mar;22(6):708-15. btk001. 10.1093/bioinformatics/btk001.

    abstract

    The general-time-reversible (GTR) model is one of the most popular models of nucleotide substitution because it constitutes a good trade-off between mathematical tractability and biological reality. However, when it is applied for inferring evolutionary distances and/or instantaneous rate matrices, the GTR model seems more prone to inapplicability than more restrictive time-reversible models. Although it has been previously noted that the causes for intractability are caused by the impossibility of computing the logarithm of a matrix characterised by negative eigenvalues, the issue has not been investigated further.

    view more details on Pubmed

  • Assessing the applicability of the GTR nucleotide substitution model through simulations.

    Evol. Bioinform. Online 2006 ;2():145-55. PMC2674669.

    abstract

    The General Time Reversible (GTR) model of nucleotide substitution is at the core of many distance-based and character-based phylogeny inference methods. The procedure described by Waddell and Steel (1997), for estimating distances and instantaneous substitution rate matrices, R, under the GTR model, is known to be inapplicable under some conditions, ie, it leads to the inapplicability of the GTR model. Here, we simulate the evolution of DNA sequences along 12 trees characterized by different combinations of tree length, (non-)homogeneity of the substitution rate matrix R, and sequence length. We then evaluate both the frequency of the GTR model inapplicability for estimating distances and the accuracy of inferred alignments. Our results indicate that, inapplicability of the Waddel and Steel's procedure can be considered a real practical issue, and illustrate that the probability of this inapplicability is a function of substitution rates and sequence length.We also discuss the implications of our results on the current implementations of maximum likelihood and Bayesian methods.

    view more details on Pubmed

  • OligoFaktory: a visual tool for interactive oligonucleotide design.

    Bioinformatics 2006 Jan;22(1):115-6. bti728. 10.1093/bioinformatics/bti728.

    abstract

    The OligoFaktory is a set of tools for the design, on an arbitrary number of target sequences, of high-quality long oligonucleotide for micro-array, of primer pair for PCR, of siRNA and more. The user-centered interface exists in two flavours: a web portal and a standalone software for Mac OS X Tiger. A unified presentation of results provides overviews with distribution charts and relative location bar graphs, as well as detailed features for each oligonucleotide. Input and output files conform to a common XML interchange file format to allow both automatic generation of input data, archiving, and post-processing of results. The design pipeline can use BLAST servers to evaluate specificity of selected oligonucleotides.

    view more details on Pubmed

  • Evaluating intraspecific "network" construction methods using simulated sequence data: do existing algorithms outperform the global maximum parsimony approach?

    Syst. Biol. 2005 Jun;54(3):363-72. Q512W89524N36102. 10.1080/10635150590945377.

    abstract

    In intraspecific studies, reticulated graphs are valuable tools for visualization, within a single figure, of alternative genealogical pathways among haplotypes. As available software packages implementing the global maximum parsimony (MP) approach only give the possibility to merge resulting topologies into less-resolved consensus trees, MP has often been neglected as an alternative approach to purely algorithmic (i.e., methods defined solely on the basis of an algorithm) "network" construction methods. Here, we propose to search tree space using the MP criterion and present a new algorithm for uniting all equally most parsimonious trees into a single (possibly reticulated) graph. Using simulated sequence data, we compare our method with three purely algorithmic and widely used graph construction approaches (minimum-spanning network, statistical parsimony, and median-joining network). We demonstrate that the combination of MP trees into a single graph provides a good estimate of the true genealogy. Moreover, our analyses indicate that, when internal node haplotypes are not sampled, the median-joining and MP methods provide the best estimate of the true genealogy whereas the minimum-spanning algorithm shows very poor performances.

    view more details on Pubmed

  • Separate-component-stabilization system for protein and DNA production without the use of antibiotics.

    BioTechniques 2005 May;38(5):775-81. 05385RR02.

    abstract

    Plasmid instability is a significant concern in the industrial utilization of microorganisms for protein or DNA production. Here we report on the development of a new and highly effective stabilization system based on the use of the ccd antidote/poison genes. For the first time, we separated the antidote gene from the poison gene: localizing the former in the plasmid and integrating the latter in the bacterial chromosome. We show that this separate-component-stabilization (SCS) strategy: (i) allows for perfect stabilization without the use of antibiotics; (ii) increases three to five times the recombinant protein production levels; and (iii) does not require any specific modification of the protein production process or culture medium. We illustrate that point by using the classical T7 promotor (i.e., used in most expression systems). Finally, we demonstrate that the SCS system increases by five the yield in DNA production, a result especially important for the design and production of gene therapy constructs void of any antibiotic resistance gene.

    view more details on Pubmed

  • Inferring contemporary levels of gene flow and demographic history in a local population of the leaf beetle Gonioctena olivacea from mitochondrial DNA sequence variation.

    Mol. Ecol. 2005 May;14(6):1641-53. MEC2537. 10.1111/j.1365-294X.2005.02537.x.

    abstract

    We have studied mitochondrial DNA variation in a local population of the leaf beetle species Gonioctena olivacea, to check whether its apparent low dispersal behaviour affects its pattern of genetic variation at a small geographical scale. We have sampled 10 populations of G. olivacea within a rectangle of 5 x 2 km in the Belgian Ardennes, as well as five populations located approximately along a straight line of 30 km and separated by distances of 3-12 km. For each sampled individual (8-19 per population), a fragment of the mtDNA control region was polymerase chain reaction-amplified and sequenced. Sequence data were analysed to test whether significant genetic differentiation could be detected among populations separated by such relatively short distances. The reconstructed genealogy of the mitochondrial haplotypes was also used to investigate the demographic history of these populations. Computer simulations of the evolution of populations were conducted to assess the minimum amount of gene flow that is necessary to explain the observed pattern of variation in the samples. Results show that migration among populations included in the rectangle of 5 x 2 km is substantial, and probably involves the occurrence of dispersal flights. This appears difficult to reconcile with the results of a previous ecological field study that concluded that most of this species dispersal occurs by walking. While sufficient migration to homogenize genetic diversity occurs among populations separated by distances of a few hundred metres to a few kilometres, distances greater than 5 km results in contrast in strong differentiation among populations, suggesting that migration is drastically reduced on such distances. Finally, the results of coalescent simulations suggest that the star-like genealogy inferred from the mtDNA sequence data is fully compatible with a past demographic expansion. However, a metapopulation structure alone (without the need to invoke a population expansion event) cannot be dismissed as the cause of this star shape.

    view more details on Pubmed

  • Evidence for male dispersal along the coasts but no migration in pelagic waters in dusky dolphins (Lagenorhynchus obscurus).

    Mol. Ecol. 2005 Jan;14(1):107-21. MEC2407. 10.1111/j.1365-294X.2004.02407.x.

    abstract

    Using nine nuclear species-specific microsatellite loci and two mitochondrial gene fragments (cytochrome b and control region), we investigated the processes that have shaped the geographical distribution of genetic diversity exhibited by contemporary dusky dolphin (Lagenorhynchus obscurus) populations. A total of 221 individuals from four locations (Peru, Argentina, southern Africa, and New Zealand) were assayed, covering most of the species' distribution range. Although our analyses identify a general demographic decline in the Peruvian dusky dolphin stock (recently affected by high natural and human-induced mortality levels), comparison between the different molecular markers hint at an ancient bottleneck that predates recent El Niño oscillations and human exploitation. Moreover, we find evidence of a difference in dispersal behaviour of dusky dolphins along the South American coast and across the Atlantic. While data in Peruvian and Argentine waters are best explained by male-specific gene flow between these two populations, our analyses suggest that dusky dolphins from Argentina and southern Africa recently separated from an ancestral Atlantic population and, since then, diverged without considerable gene flow. The inclusion of a few New Zealand samples further confirms the low levels of genetic differentiation among most dusky dolphin populations. Only the Peruvian dusky dolphin stock is highly differentiated, especially at mitochondrial loci, suggesting that major fluctuations in its population size have led to an increased rate of genetic drift.

    view more details on Pubmed

  • Local endemism within the Western Ghats-sri Lanka biodiversity hotspot.

    Science 2004 Oct;306(5695):479-81. 306/5695/479. 10.1126/science.1100167.

    abstract

    The apparent biotic affinities between the mainland and the island in the Western Ghats-Sri Lanka biodiversity hotspot have been interpreted as the result of frequent migrations during recent periods of low sea level. We show, using molecular phylogenies of two invertebrate and four vertebrate groups, that biotic interchange between these areas has been much more limited than hitherto assumed. Despite several extended periods of land connection during the past 500,000 years, Sri Lanka has maintained a fauna that is largely distinct from that of the Indian mainland. Future conservation programs for the subcontinent should take into account such patterns of local endemism at the finest scale at which they may occur.

    view more details on Pubmed

  • Genetic analysis of a successful repatriation programme: giant Galápagos tortoises.

    Proc. Biol. Sci. 2004 Feb;271(1537):341-5. 10.1098/rspb.2003.2607. PMC1691607.

    abstract

    As natural populations of endangered species dwindle to precarious levels, remaining members are sometimes brought into captivity, allowed to breed and their offspring returned to the natural habitat. One goal of such repatriation programmes is to retain as much of the genetic variation of the species as possible. A taxon of giant Galápagos tortoises on the island of Española has been the subject of a captive breeding-repatriation programme for 33 years. Core breeders, consisting of 12 females and three males, have produced more than 1200 offspring that have been released on Española where in situ reproduction has recently been observed. Using microsatellite DNA markers, we have determined the maternity and paternity of 132 repatriated offspring. Contributions of the breeders are highly skewed. This has led to a further loss of genetic variation that is detrimental to the long-term survival of the population. Modifications to the breeding programme could alleviate this problem.

    view more details on Pubmed

  • Molecular phylogenetic analyses indicate a wide and ancient radiation of African hepatitis delta virus, suggesting a deltavirus genus of at least seven major clades.

    J. Virol. 2004 Mar;78(5):2537-44. PMC369207.

    abstract

    Hepatitis D virus (HDV) is a satellite of hepatitis B virus (HBV) for transmission and propagation and infects nearly 20 million people worldwide. The HDV genome is a compact circular single-stranded RNA genome with extensive intramolecular complementarity. Despite its different epidemiological and pathological patterns, the variability and geographical distribution of HDV are limited to three genotypes and two subtypes that have been characterized to date. Phylogenetic reconstructions based on the delta antigen gene and full-length genome sequence data show an extensive and probably ancient radiation of African lineages, suggesting that the genetic variability of HDV is much more complex than was previously thought, with evidence of additional clades. These results relate the geographic distribution of HDV more closely to the genetic variability of its helper HBV.

    view more details on Pubmed

  • A hidden Markov model for progressive multiple alignment.

    Bioinformatics 2003 Aug;19(12):1505-13.

    abstract

    Progressive algorithms are widely used heuristics for the production of alignments among multiple nucleic-acid or protein sequences. Probabilistic approaches providing measures of global and/or local reliability of individual solutions would constitute valuable developments.

    view more details on Pubmed

  • Structure and evolution of the mitochondrial control region of leaf beetles (Coleoptera: Chrysomelidae): a hierarchical analysis of nucleotide sequence variation.

    J. Mol. Evol. 2003 Jan;56(1):38-45. 10.1007/s00239-002-2378-1.

    abstract

    To assess the levels of variation at different evolutionary scales in the mitochondrial (mt) control region of leaf beetles, we sequenced and compared the full mt control region in two genera ( Chrysomela and Gonioctena), in two species within a genus ( Gonioctena olivacea and G. pallida), in individuals from distant populations of these species in Europe, and in individuals from populations separated by moderate (10- to 100-km) to short (<5-km) distances. In all individuals, a highly repetitive section consisting of the tandem repetition of 12 to 17 imperfect copies of a 107- to 159-bp-long core sequence was observed. This repetitive fragment accounts for roughly 50% of the full control-region length. The sequence variability among repeated elements within the control region of a given individual depends on the species considered: the variability within any G. olivacea individual is much higher than that within G. pallida individuals. Comparisons of the repeated elements, in a phylogenetic framework, within and among individuals of G. olivacea and G. pallida suggests that the repetitive section of the control region experienced recurrent duplications/deletions, leading to some degree of concerted evolution. Comparisons between Chrysomela and Gonioctena control regions revealed virtually no significant sequence similarity, except for two long stretches of A's and several [T(T)A(A)] repeats, all found in the control region of other insect orders. Our analyses allowed us to identify portions of the control region with enough variation for population genetic or phylogeographic studies.

    view more details on Pubmed

  • Microsatellite analysis of genetic divergence among populations of giant Galápagos tortoises.

    Mol. Ecol. 2002 Nov;11(11):2265-83. 1617.

    abstract

    Giant Galápagos tortoises represent an interesting model for the study of patterns of genetic divergence and adaptive differentiation related to island colonization events. Recent mitochondrial DNA work elucidated the evolutionary history of the species and helped to clarify aspects of nomenclature. We used 10 microsatellite loci to assess levels of genetic divergence among and within island populations. In particular, we described the genetic structure of tortoises on the island of Isabela, where discrimination of different taxa is still subject of debate. Individual island populations were all genetically distinct. The island of Santa Cruz harboured two distinct populations. On Isabela, populations of Volcan Wolf, Darwin and Alcedo were significantly different from each other. On the other hand, Volcan Wolf showed allelic similarity with the island of Santiago. On Southern Isabela, lower genetic divergence was found between Northeast Sierra Negra and Volcan Alcedo, while patterns of gene flow were recorded among tortoises of Cerro Azul and Southeast Sierra Negra. These tortoises have endured heavy exploitation during the last three centuries and recently attracted much concern due to the current number of stochastic and deterministic threats to extant populations. Our study complements previous investigation based on mtDNA diversity and provides further information that may help devising tortoise management plans.

    view more details on Pubmed

  • Investigation of the population genetic structure and mating system in the ant Pheidole pallidula.

    Mol. Ecol. 2002 Sep;11(9):1805-14. 1573.

    abstract

    The origin of eusociality in haplo-diploid organisms such as Hymenoptera has been mostly explained by kin selection. However, several studies have uncovered decreased relatedness values within colonies, resulting primarily from multiple queen matings (polyandry) and/or from the presence of more than one functional queen (polygyny). Here, we report on the use of microsatellite data for the investigation of sociogenetic parameters, such as relatedness, and levels of polygyny and polyandry, in the ant Pheidole pallidula. We demonstrate, through analysis of mother-offspring combinations and the use of direct sperm typing, that each queen is inseminated by a single male. The inbreeding coefficient within colonies and the levels of relatedness between the queens and their mate are not significantly different from zero, indicating that matings occur between unrelated individuals. Analyses of worker genotypes demonstrate that 38% of the colonies are polygynous with 2-4 functional queens, and suggest the existence of reproductive skew, i.e. unequal respective contribution of queens to reproduction. Finally, our analyses indicate that colonies are genetically differentiated and form a population exhibiting significant isolation-by-distance, suggesting that some colonies originate through budding.

    view more details on Pubmed

  • The metapopulation genetic algorithm: An efficient solution for the problem of large phylogeny estimation.

    Proc. Natl. Acad. Sci. U.S.A. 2002 Aug;99(16):10516-21. 10.1073/pnas.162224399. 162224399. PMC124960.

    abstract

    Large phylogeny estimation is a combinatorial optimization problem that no future computer will ever be able to solve exactly in practical computing time. The difficulty of the problem is amplified by the need to use complex evolutionary models and large taxon samplings. Hence, many heuristic approaches have been developed, with varying degrees of success. Here, we report on a heuristic approach, the metapopulation genetic algorithm, involving several populations of trees that are forced to cooperate in the search for the optimal tree. Within each population, trees are subjected to evaluation, selection, and mutation events, which are directed by using inter-population consensus information. The method proves to be both very accurate and vastly faster than existing heuristics, such that data sets comprised of hundreds of taxa can be analyzed in practical computing times under complex models of maximum-likelihood evolution. Branch support values produced by the metapopulation genetic algorithm might closely approximate the posterior probabilities of the corresponding branches.

    view more details on Pubmed

  • Dual chemical sequestration: a key mechanism in transitions among ecological specialization.

    Proc. Biol. Sci. 2002 Jan;269(1486):1-6. 10.1098/rspb.2001.1859. PMC1690864.

    abstract

    Platyphora leaf beetles form a vast group of tropical species each feeding on a restricted set of host plants and exhibiting bright coloration warning predators against their chemical protection. These beetles offer an exceptional opportunity for understanding the evolution of phytochemical sequestration. Indeed, qualitative studies of defensive secretions indicate that Platyphora species acquire toxicity via sequestration of plant secondary metabolites. All produce pentacyclic triterpene saponins from sequestered plant amyrins, but our analyses also indicate that many Platyphora species produce a dual chemical defence, that is, they are additionally protected by lycopsamine-type pyrolyzidine alkaloids that they also sequester from their host. This paper reports on the evolution of chemical defence and host affiliation in Platyphora leaf beetles as reconstructed on a molecular phylogeny of mitochondrial and nuclear DNA sequences. The analyses indicate that dual sequestration could be the key mechanistic means by which transitions among ecological specializations (i.e. restricted host-plant affiliations) are made possible.

    view more details on Pubmed


Subunits