Evo-Devo of skin appendages

abstract

Corn snakes are emerging models for animal colouration studies. Here, we focus on the Terrazzo morph, whose skin pattern is characterized by stripes rather than blotches. Using genome mapping, we discover a disruptive mutation in the coding region of the Premelanosome protein (PMEL) gene. Our transcriptomic analyses reveal that PMEL expression is significantly downregulated in Terrazzo embryonic tissues. We produce corn snake PMEL knockouts, which present a comparable colouration phenotype to Terrazzo and the subcellular structure of their melanosomes and xanthosomes is also similarly impacted. Our single-cell expression analyses of wild-type embryonic dorsal skin demonstrate that all chromatophore progenitors express PMEL at varying levels. Finally, we show that in wild-type embryos PMEL-expressing cells are initially uniformly spread before forming aggregates and eventually blotches, as seen in the adults. In Terrazzo embryos, the aggregates fail to form. Our results provide insights into the mechanisms governing colouration patterning in reptiles.

view more details on Pubmed

abstract

Reptilian skin coloration is spectacular and diverse, yet little is known about the ontogenetic processes that govern its establishment and the molecular signaling pathways that determine it. Here, we focus on the development of the banded pattern of leopard gecko hatchlings and the transition to black spots in the adult. With our histological analyses, we show that iridophores are present in the white and yellow bands of the hatchling and they gradually perish in the adult skin. Furthermore, we demonstrate that melanophores can autonomously form spots in the absence of the other chromatophores both on the regenerated skin of the tail and on the dorsal skin of the Mack Super Snow (MSS) leopard geckos. This color morph is characterized by uniform black coloration in hatchlings and black spots in adulthood; we establish that their skin is devoid of xanthophores and iridophores at both stages. Our genetic analyses identified a 13-nucleotide deletion in the PAX7 transcription factor of MSS geckos, affecting its protein coding sequence. With our single-cell transcriptomics analysis of embryonic skin, we confirm that PAX7 is expressed in iridophores and xanthophores, suggesting that it plays a key role in the differentiation of both chromatophores. Our in situ hybridizations on whole-mount embryos document the dynamics of the skin pattern formation and how it is impacted in the PAX7 mutants. We hypothesize that the melanophores-iridophores interactions give rise to the banded pattern of the hatchlings and black spot formation is an intrinsic capacity of melanophores in the postembryonic skin.

view more details on Pubmed

abstract

In some mammals, notably humans, recombination occurs almost exclusively where the protein PRDM9 binds, whereas in vertebrates lacking an intact , such as birds and canids, recombination rates are elevated near promoter-like features. To determine whether PRDM9 directs recombination in nonmammalian vertebrates, we focused on an exemplar species with a single, intact ortholog, the corn snake (). Analyzing historical recombination rates along the genome and crossovers in pedigrees, we found evidence that PRDM9 specifies the location of recombination events, but we also detected a separable effect of promoter-like features. These findings reveal that the uses of PRDM9 and promoter-like features need not be mutually exclusive and instead reflect a tug-of-war that is more even in some species than others.

view more details on Pubmed

abstract

Reptilian species, particularly snakes and lizards, are emerging models of animal coloration. Here, I focus on the role of the TFEC transcription factor in snake and lizard coloration based on a study on wild-type and piebald ball pythons. Genomic mapping previously identified a TFEC mutation linked to the piebald ball python phenotype. The association of TFEC with skin coloration was further supported by gene-editing experiments in the brown anole lizard. However, novel histological analyses presented here reveal discrepancies between the ball python and the anole TFEC mutants phenotype, cautioning against broad generalizations. Indeed, both wild-type and piebald ball pythons completely lack iridophores, whereas the TFEC anole lizard mutants lose their iridophores compared to the wild-type anole. Based on these findings, I discuss the potential role of the MiT/TFE family in skin pigmentation across vertebrate lineages and advocate the need for developmental analyses and additional gene-editing experiments to explore the reptilian coloration diversity.

view more details on Pubmed

abstract

Changes in gene expression represent an important source of phenotypic innovation. Yet how such changes emerge and impact the evolution of traits remains elusive. Here, we explore the molecular mechanisms associated with the development of masculinizing ovotestes in female moles. By performing integrative analyses of epigenetic and transcriptional data in mole and mouse, we identified the co-option of SALL1 expression in mole ovotestes formation. Chromosome conformation capture analyses highlight a striking conservation of the 3D organization at the SALL1 locus, but an evolutionary divergence of enhancer activity. Interspecies reporter assays support the capability of mole-specific enhancers to activate transcription in urogenital tissues. Through overexpression experiments in transgenic mice, we further demonstrate the capability of SALL1 to induce kidney-related gene programs, which are a signature of mole ovotestes. Our results highlight the co-option of gene expression, through changes in enhancer activity, as a plausible mechanism for the evolution of traits.

view more details on Pubmed

abstract

In vertebrates, there are two known mechanisms by which meiotic recombination is directed to the genome: in humans, mice, and other mammals, recombination occurs almost exclusively where the protein PRDM9 binds, while in species lacking an intact , such as birds and canids, recombination rates are elevated near promoter-like features. To test if PRDM9 also directs recombination in non-mammalian vertebrates, we focused on an exemplar species, the corn snake (). Unlike birds, this species possesses a single, intact ortholog. By inferring historical recombination rates along the genome from patterns of linkage disequilibrium and identifying crossovers in pedigrees, we found that PRDM9 specifies the location of recombination events outside of mammals. However, we also detected an independent effect of promoter-like features on recombination, which is more pronounced on macro- than microchromosomes. Thus, our findings reveal that the uses of PRDM9 and promoter-like features are not mutually-exclusive, and instead reflect a tug of war, which varies in strength along the genome and is more lopsided in some species than others.

view more details on Pubmed

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

abstract

Albinism and leucism are phenotypes resulting from impaired melanin pigmentation in the skin and skin appendages. However, melanin pigmentation of eyes remains unaffected in leucism. Here, using transmission electron microscopy, we show that the leucistic morph of the Texas rat snake (Pantherophis obsoletus lindheimeri) lacks both melanophores and xanthophores in its skin and exhibits a uniform ivory white color generated by iridophores and collagen fibers. In addition, we sequenced the full genome of a leucistic individual and obtained a highly-contiguous near-chromosome quality assembly of 1.69 Gb with an N50 of 14.5 Mb and an L50 of 29 sequences. Using a candidate-gene approach, we then identify in the leucistic genome a single-nucleotide deletion that generates a frameshift and a premature termination codon in the melanocyte inducing transcription factor (MITF) gene. This mutation shortens the translated protein from 574 to 286 amino acids, removing the helix-loop-helix DNA-binding domain that is highly conserved among vertebrates. Genotyping leucistic animals of independent lineages showed that not all leucistic individuals carry this single-nucleotide deletion. Subsequent gene expression analyses reveal that all leucistic individuals that we analyzed exhibit a significantly decreased expression of MITF. We thus suggest that mutations affecting the regulation and, in some cases, the coding sequence of MITF, the former probably predating the latter, could be associated with the leucistic phenotype in Texas rat snakes. MITF is involved in the development and survival of melanophores in vertebrates. In zebrafish, a classical model species for pigmentation that undergoes metamorphosis, larvae and adults of homozygous mitfa mutants lack melanophores, show an excess of iridophores and exhibit reduced yellow pigmentation. On the contrary, in the leucistic Texas rat snake, a non-metamorphic species, only iridophores persist. Our results suggest that fate determination of neural-crest derived melanophores and xanthophores, but not of iridophores, could require the expression of MITF during snake embryonic development.

see on external website

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

abstract

PMEL is a pigment cell-specific protein that forms a functional amyloid matrix in melanosomes. The matrix consists of well-separated fibrillar sheets on which the pigment melanin is deposited. Using electron tomography, we demonstrate that this sheet architecture is governed by the PMEL repeat (RPT) domain, which associates with the amyloid as an accessory proteolytic fragment. Thus, the RPT domain is dispensable for amyloid formation as such but shapes the morphology of the matrix, probably in order to maximize the surface area available for pigment adsorption. Although the primary amino acid sequence of the RPT domain differs vastly among various vertebrates, we show that it is a functionally conserved, interchangeable module. RPT domains of all species are predicted to be very highly O-glycosylated, which is likely the common defining feature of this domain. O-glycosylation is indeed essential for RPT domain function and the establishment of the PMEL sheet architecture. Thus, O-glycosylation, not amino acid sequence, appears to be the major factor governing the characteristic PMEL amyloid morphology.

view more details on Pubmed

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

abstract

The study of chondrocrania has a long tradition with a focus on single specimens and stages. It revealed great interspecific diversity and a notion of intraspecific variation. As an embryonic structure, the chondrocranium is subject to major changes in ontogeny with resorption and ossification of different cartilaginous structures. The cupula nasi anterior is the anteriormost portion of the cartilaginous nasal capsule and is expected to mirror much of the animal's life history and lifestyle. Its diversity in mammals is reflected in the external nasal anatomy of newborns. Marsupials and placentals show marked differences, likely related to breathing and suckling behavior.

view more details on Pubmed

abstract

Cerebral cortex size differs dramatically between reptiles, birds, and mammals, owing to developmental differences in neuron production. In mammals, signaling pathways regulating neurogenesis have been identified, but genetic differences behind their evolution across amniotes remain unknown. We show that direct neurogenesis from radial glia cells, with limited neuron production, dominates the avian, reptilian, and mammalian paleocortex, whereas in the evolutionarily recent mammalian neocortex, most neurogenesis is indirect via basal progenitors. Gain- and loss-of-function experiments in mouse, chick, and snake embryos and in human cerebral organoids demonstrate that high Slit/Robo and low Dll1 signaling, via Jag1 and Jag2, are necessary and sufficient to drive direct neurogenesis. Attenuating Robo signaling and enhancing Dll1 in snakes and birds recapitulates the formation of basal progenitors and promotes indirect neurogenesis. Our study identifies modulation in activity levels of conserved signaling pathways as a primary mechanism driving the expansion and increased complexity of the mammalian neocortex during amniote evolution.

view more details on Pubmed

abstract

Multiple technologies and software are now available facilitating the de novo sequencing and assembly of any vertebrate genome. Yet the quality of most available sequenced genomes is substantially poorer than that of the golden standard in the field: the human genome. Here, we present a step-by-step protocol for the successful sequencing and assembly of a high-quality snake genome that can be applied to any other reptilian or avian species. We combine the great sequencing depth and accuracy of short reads with the use of different insert size libraries for extended scaffolding followed by optical mapping. We show that this procedure improved the corn snake scaffold N50 from 3.7 kbp to 1.4 Mbp, currently making it one of the snake genomes with the longest scaffolds. Short guidelines are also given on the extraction of long DNA molecules from reptilian blood and the necessary modifications in DNA extraction protocols. This chapter is accompanied by a website ( www.reptilomics.org/stepbystep.html ), where we provide links to the suggested software, examples of input and output files, and running parameters.

view more details on Pubmed

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

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

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

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

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

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

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

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

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

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

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

abstract

Until recently, facultative automictic parthenogenesis within the squamate reptiles exhibiting ZZ:ZW genetic sex determination has resulted in single reproductive events producing male (ZZ) or female (ZW) offspring. With the recent discovery of viable parthenogenetically produced female (WW) Boa constrictors, the existence of further parthenogenetic events resulting in WW females was questioned. Here, we provide genetic evidence for consecutive virgin births by a female Colombian rainbow boa (Epicrates maurus), resulting in the production of WW females likely through terminal fusion automixis. Samples were screened at 22 microsatellite loci with 12 amplifying unambiguous products. Of these, maternal heterozygosity was observed in 4, with the offspring differentially homozygous at each locus. This study documents the first record of parthenogenesis within the genus Epicrates, a second within the serpent lineage Boidae, and the third genetically confirmed case of consecutive virgin births of viable offspring within any vertebrate lineage. Unlike the recent record in Boa constrictors, the female described here was isolated from conspecifics from birth, demonstrating that males are not required to stimulate parthenogenetic reproduction in this species and possibly other Boas.

view more details on Pubmed

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

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

abstract

Vertebrate skin is characterized by its patterned array of appendages, whether feathers, hairs, or scales. In avian skin the distribution of feathers occurs on two distinct spatial levels. Grouping of feathers within discrete tracts, with bare skin lying between the tracts, is termed the macropattern, while the smaller scale periodic spacing between individual feathers is referred to as the micropattern. The degree of integration between the patterning mechanisms that operate on these two scales during development and the mechanisms underlying the remarkable evolvability of skin macropatterns are unknown. A striking example of macropattern variation is the convergent loss of neck feathering in multiple species, a trait associated with heat tolerance in both wild and domestic birds. In chicken, a mutation called Naked neck is characterized by a reduction of body feathering and completely bare neck. Here we perform genetic fine mapping of the causative region and identify a large insertion associated with the Naked neck trait. A strong candidate gene in the critical interval, BMP12/GDF7, displays markedly elevated expression in Naked neck embryonic skin due to a cis-regulatory effect of the causative mutation. BMP family members inhibit embryonic feather formation by acting in a reaction-diffusion mechanism, and we find that selective production of retinoic acid by neck skin potentiates BMP signaling, making neck skin more sensitive than body skin to suppression of feather development. This selective production of retinoic acid by neck skin constitutes a cryptic pattern as its effects on feathering are not revealed until gross BMP levels are altered. This developmental modularity of neck and body skin allows simple quantitative changes in BMP levels to produce a sparsely feathered or bare neck while maintaining robust feather patterning on the body.

view more details on Pubmed

abstract

Vertebrate skin is characterized by its patterned array of appendages, whether feathers, hairs, or scales. In avian skin the distribution of feathers occurs on two distinct spatial levels. Grouping of feathers within discrete tracts, with bare skin lying between the tracts, is termed the macropattern, while the smaller scale periodic spacing between individual feathers is referred to as the micropattern. The degree of integration between the patterning mechanisms that operate on these two scales during development and the mechanisms underlying the remarkable evolvability of skin macropatterns are unknown. A striking example of macropattern variation is the convergent loss of neck feathering in multiple species, a trait associated with heat tolerance in both wild and domestic birds. In chicken, a mutation called Naked neck is characterized by a reduction of body feathering and completely bare neck. Here we perform genetic fine mapping of the causative region and identify a large insertion associated with the Naked neck trait. A strong candidate gene in the critical interval, BMP12/GDF7, displays markedly elevated expression in Naked neck embryonic skin due to a cis-regulatory effect of the causative mutation. BMP family members inhibit embryonic feather formation by acting in a reaction-diffusion mechanism, and we find that selective production of retinoic acid by neck skin potentiates BMP signaling, making neck skin more sensitive than body skin to suppression of feather development. This selective production of retinoic acid by neck skin constitutes a cryptic pattern as its effects on feathering are not revealed until gross BMP levels are altered. This developmental modularity of neck and body skin allows simple quantitative changes in BMP levels to produce a sparsely feathered or bare neck while maintaining robust feather patterning on the body.

see on external website

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

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

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

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

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

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

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

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

abstract

see on external website

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

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

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

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

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

abstract

Two species of the brine shrimp, namely Artemia franciscana Kellogg and A. persimilis Piccinelli and Prosdocimi, inhabit Chile. Most studies so far have shown that A. franciscana is the most widely distributed species in Chile, with A. persimilis present only in Chilean Patagonia. In general, there is good agreement between morphological and genetic comparisons of Chilean populations with respect to species discrimination. However, a number of results indicate an overlap with some populations tending to diverge from A. franciscana and/or resembling A. persimilis. Prior to the mid 90's the use of DNA markers in Artemia was rather limited, despite their successful application in numerous other species. In this study, we investigate whether the conclusions drawn from traditional comparative tools are congruent with the pattern of genetic divergence depicted by DNA analysis at the mitochondrial level.

see on external website