Evolution of a core gene network for skeletogenesis in chordates

Extracellular Vesicles Carrying RUNX3 Promote Differentiation of Dental Pulp Stem Cells

BACKGROUND

This study aims to clarify the mechanism underlying dental pulp cells-extracellular vesicles (DPC-EVs) carrying runt-related transcription factor 3 (RUNX3) in mediating odontogenic differentiation of dental pulp stem cells (DPSCs) with the involvement of miR-30a-5p-regulated NOTCH1.

METHODS

Extracellular vesicles (EVs) were isolated from human DPSCs, and identified using transmission electron microscopy, and nanoparticle tracking analysis. PBS, EVs, or EV inhibitor GW4869 was added to DPSCs for co-culture, whilst odontogenic differentiation was assessed in terms of ratio of mineralized nodules and expression odontoblast differentiation markers. Dual luciferase reporter gene assay and chromatin immunoprecipitation for binding relation among RUNX3, miR-30a-5p and NOTCH1were employed to evaluate their roles in odontogenic differentiation was determined. Animal experiment was established to confirm the effect of DPC-EVs-loaded RUNX3 on dental pulp.

RESULTS

In vitro finding demonstrated that EVs delivered RUNX3 to DPSCs, thereby activated miR-30a-5p expression and inhibited NOTCH1 expression, which was reversed by addition of GW4869. RUNX3 upregulation promoted miR-30a-5p while miR-30a-5p targeted and inhibited NOTCH1. Silencing of RUNX3 in EVs decreased expression of those differentiation markers, downregulated miR-30a-5p and upregulated NOTCH1.

CONCLUSION

DPSC-EVs can carry RUNX3 to the DPSCs, promote the transcription of miR-30a-5p, and then inhibit the expression of NOTCH1, and finally promote the odontogenic differentiation of DPSCs.

Teeth outside the mouth: The evolution and development of shark denticles

Vertebrate skin appendages are incredibly diverse. This diversity, which includes structures such as scales, feathers, and hair, likely evolved from a shared anatomical placode, suggesting broad conservation of the early development of these organs. Some of the earliest known skin appendages are dentine and enamel-rich tooth-like structures, collectively known as odontodes. These appendages evolved over 450 million years ago. Elasmobranchs (sharks, skates, and rays) have retained these ancient skin appendages in the form of both dermal denticles (scales) and oral teeth. Despite our knowledge of denticle function in adult sharks, our understanding of their development and morphogenesis is less advanced. Even though denticles in sharks appear structurally similar to oral teeth, there has been limited data directly comparing the molecular development of these distinct elements. Here, we chart the development of denticles in the embryonic small-spotted catshark (Scyliorhinus canicula) and characterize the expression of conserved genes known to mediate dental development. We find that shark denticle development shares a vast gene expression signature with developing teeth. However, denticles have restricted regenerative potential, as they lack a sox2+ stem cell niche associated with the maintenance of a dental lamina, an essential requirement for continuous tooth replacement. We compare developing denticles to other skin appendages, including both sensory skin appendages and avian feathers. This reveals that denticles are not only tooth-like in structure, but that they also share an ancient developmental gene set that is likely common to all epidermal appendages.

The Evolution of Biomineralization through the Co-Option of Organic Scaffold Forming Networks

Biomineralization is the process in which organisms use minerals to generate hard structures like teeth, skeletons and shells. Biomineralization is proposed to have evolved independently in different phyla through the co-option of pre-existing developmental programs. Comparing the gene regulatory networks (GRNs) that drive biomineralization in different species could illuminate the molecular evolution of biomineralization. Skeletogenesis in the sea urchin embryo was extensively studied and the underlying GRN shows high conservation within echinoderms, larval and adult skeletogenesis. The organic scaffold in which the calcite skeletal elements form in echinoderms is a tubular compartment generated by the syncytial skeletogenic cells. This is strictly different than the organic cartilaginous scaffold that vertebrates mineralize with hydroxyapatite to make their bones. Here I compare the GRNs that drive biomineralization and tubulogenesis in echinoderms and in vertebrates. The GRN that drives skeletogenesis in the sea urchin embryo shows little similarity to the GRN that drives bone formation and high resemblance to the GRN that drives vertebrates’ vascular tubulogenesis. On the other hand, vertebrates’ bone-GRNs show high similarity to the GRNs that operate in the cells that generate the cartilage-like tissues of basal chordate and invertebrates that do not produce mineralized tissue. These comparisons suggest that biomineralization in deuterostomes evolved through the phylum specific co-option of GRNs that control distinct organic scaffolds to mineralization.

Mineralized Cartilage and Bone-Like Tissues in Chondrichthyans Offer Potential Insights Into the Evolution and Development of Mineralized Tissues in the Vertebrate Endoskeleton

The impregnation of biominerals into the extracellular matrix of living organisms, a process termed biomineralization, gives rise to diverse mineralized (or calcified) tissues in vertebrates. Preservation of mineralized tissues in the fossil record has provided insights into the evolutionary history of vertebrates and their skeletons. However, current understanding of the vertebrate skeleton and of the processes underlying its formation is biased towards biomedical models such as the tetrapods mouse and chick. Chondrichthyans (sharks, skates, rays, and chimaeras) and osteichthyans are the only vertebrate groups with extant (living) representatives that have a mineralized skeleton, but the basal phylogenetic position of chondrichthyans could potentially offer unique insights into skeletal evolution. For example, bone is a vertebrate novelty, but the internal supporting skeleton (endoskeleton) of extant chondrichthyans is commonly described as lacking bone. The molecular and developmental basis for this assertion is yet to be tested. Subperichondral tissues in the endoskeleton of some chondrichthyans display mineralization patterns and histological and molecular features of bone, thereby challenging the notion that extant chondrichthyans lack endoskeletal bone. Additionally, the chondrichthyan endoskeleton demonstrates some unique features and others that are potentially homologous with other vertebrates, including a polygonal mineralization pattern, a trabecular mineralization pattern, and an unconstricted perichordal sheath. Because of the basal phylogenetic position of chondrichthyans among all other extant vertebrates with a mineralized skeleton, developmental and molecular studies of chondrichthyans are critical to flesh out the evolution of vertebrate skeletal tissues, but only a handful of such studies have been carried out to date. This review discusses morphological and molecular features of chondrichthyan endoskeletal tissues and cell types, ultimately emphasizing how comparative embryology and transcriptomics can reveal homology of mineralized skeletal tissues (and their cell types) between chondrichthyans and other vertebrates.

The Genome of the "Sea Vomit" Didemnum vexillum

Tunicates are the sister group of vertebrates and thus occupy a key position for investigations into vertebrate innovations as well as into the consequences of the vertebrate-specific genome duplications. Nevertheless, tunicate genomes have not been studied extensively in the past, and comparative studies of tunicate genomes have remained scarce. The carpet sea squirt Didemnum vexillum, commonly known as "sea vomit", is a colonial tunicate considered an invasive species with substantial ecological and economical risk. We report the assembly of the D. vexillum genome using a hybrid approach that combines 28.5 Gb Illumina and 12.35 Gb of PacBio data. The new hybrid scaffolded assembly has a total size of 517.55 Mb that increases contig length about eightfold compared to previous, Illumina-only assembly. As a consequence of an unusually high genetic diversity of the colonies and the moderate length of the PacBio reads, presumably caused by the unusually acidic milieu of the tunic, the assembly is highly fragmented (L50 = 25,284, N50 = 6539). It is sufficient, however, for comprehensive annotations of both protein-coding genes and non-coding RNAs. Despite its shortcomings, the draft assembly of the "sea vomit" genome provides a valuable resource for comparative tunicate genomics and for the study of the specific properties of colonial ascidians.

Evolutionary developmental genetics of teeth and odontodes in jawed vertebrates: a perspective from the study of elasmobranchs

Most extant vertebrates display a high variety of tooth and tooth-like organs (odontodes) that vary in shape, position over the body and nature of composing tissues. The development of these structures is known to involve similar genetic cascades and teeth and odontodes are believed to share a common evolutionary history. Gene expression patterns have previously been compared between mammalian and teleost tooth development but we highlight how the comparative framework was not always properly defined to deal with different tooth types or tooth developmental stages. Larger-scale comparative analyses also included cartilaginous fishes: sharks display oral teeth and dermal scales for which the gene expression during development started to be investigated in the small-spotted catshark Scyliorhinus canicula during the past decade. We report several descriptive approaches to analyse the embryonic tooth and caudal scale gene expressions in S. canicula. We compare these expressions wih the ones reported in mouse molars and teleost oral and pharyngeal teeth and highlight contributions and biases that arise from these interspecific comparisons. We finally discuss the evolutionary processes that can explain the observed intra and interspecific similarities and divergences in the genetic cascades involved in tooth and odontode development in jawed vertebrates.

Evolution and expansion of the RUNX2 QA repeat corresponds with the emergence of vertebrate complexity

Runt-related transcription factor 2 (RUNX2) is critical for the development of the vertebrate bony skeleton. Unlike other RUNX family members, RUNX2 possesses a variable poly-glutamine, poly-alanine (QA) repeat domain. Natural variation within this repeat is able to alter the transactivation potential of RUNX2, acting as an evolutionary 'tuning knob' suggested to influence mammalian skull shape. However, the broader role of the RUNX2 QA repeat throughout vertebrate evolution is unknown. In this perspective, we examine the role of the RUNX2 QA repeat during skeletal development and discuss how its emergence and expansion may have facilitated the evolution of morphological novelty in vertebrates.

Evidence for a role of RUNX1 as recombinase cofactor for TCRβ rearrangements and pathological deletions in ETV6-RUNX1 ALL

T-cell receptor gene beta (TCRβ) gene rearrangement represents a complex, tightly regulated molecular mechanism involving excision, deletion and recombination of DNA during T-cell development. RUNX1, a well-known transcription factor for T-cell differentiation, has recently been described to act in addition as a recombinase cofactor for TCRδ gene rearrangements. In this work we employed a RUNX1 knock-out mouse model and demonstrate by deep TCRβ sequencing, immunostaining and chromatin immunoprecipitation that RUNX1 binds to the initiation site of TCRβ rearrangement and its homozygous inactivation induces severe structural changes of the rearranged TCRβ gene, whereas heterozygous inactivation has almost no impact. To compare the mouse model results to the situation in Acute Lymphoblastic Leukemia (ALL) we analyzed TCRβ gene rearrangements in T-ALL samples harboring heterozygous Runx1 mutations. Comparable to the Runx1+/- mouse model, heterozygous Runx1 mutations in T-ALL patients displayed no detectable impact on TCRβ rearrangements. Furthermore, we reanalyzed published sequence data from recurrent deletion borders of ALL patients carrying an ETV6-RUNX1 translocation. RUNX1 motifs were significantly overrepresented at the deletion ends arguing for a role of RUNX1 in the deletion mechanism. Collectively, our data imply a role of RUNX1 as recombinase cofactor for both physiological and aberrant deletions.

HMGA Genes and Proteins in Development and Evolution

HMGA (high mobility group A) (HMGA1 and HMGA2) are small non-histone proteins that can bind DNA and modify chromatin state, thus modulating the accessibility of regulatory factors to the DNA and contributing to the overall panorama of gene expression tuning. In general, they are abundantly expressed during embryogenesis, but are downregulated in the adult differentiated tissues. In the present review, we summarize some aspects of their role during development, also dealing with relevant studies that have shed light on their functioning in cell biology and with emerging possible involvement of HMGA1 and HMGA2 in evolutionary biology.

Evolution of dental tissue mineralization: an analysis of the jawed vertebrate SPARC and SPARC-L families

BACKGROUND

The molecular bases explaining the diversity of dental tissue mineralization across gnathostomes are still poorly understood. Odontodes, such as teeth and body denticles, are serial structures that develop through deployment of a gene regulatory network shared between all gnathostomes. Dentin, the inner odontode mineralized tissue, is produced by odontoblasts and appears well-conserved through evolution. In contrast, the odontode hypermineralized external layer (enamel or enameloid) produced by ameloblasts of epithelial origin, shows extensive structural variations. As EMP (Enamel Matrix Protein) genes are as yet only found in osteichthyans where they play a major role in the mineralization of teeth and others skeletal organs, our understanding of the molecular mechanisms leading to the mineralized odontode matrices in chondrichthyans remains virtually unknown.

RESULTS

We undertook a phylogenetic analysis of the SPARC/SPARC-L gene family, from which the EMPs are supposed to have arisen, and examined the expression patterns of its members and of major fibrillar collagens in the spotted catshark Scyliorhinus canicula, the thornback ray Raja clavata, and the clawed frog Xenopus tropicalis. Our phylogenetic analyses reveal that the single chondrichthyan SPARC-L gene is co-orthologous to the osteichthyan SPARC-L1 and SPARC-L2 paralogues. In all three species, odontoblasts co-express SPARC and collagens. In contrast, ameloblasts do not strongly express collagen genes but exhibit strikingly similar SPARC-L and EMP expression patterns at their maturation stage, in the examined chondrichthyan and osteichthyan species, respectively.

CONCLUSIONS

A well-conserved odontoblastic collagen/SPARC module across gnathostomes further confirms dentin homology. Members of the SPARC-L clade evolved faster than their SPARC paralogues, both in terms of protein sequence and gene duplication. We uncover an osteichthyan-specific duplication that produced SPARC-L1 (subsequently lost in pipidae frogs) and SPARC-L2 (independently lost in teleosts and tetrapods).Our results suggest the ameloblastic expression of the single chondrichthyan SPARC-L gene at the maturation stage reflects the ancestral gnathostome situation, and provide new evidence in favor of the homology of enamel and enameloids in all gnathostomes.

Animal Evolution: The Hard Problem of Cartilage Origins

Our skeletons evolved from cartilaginous tissue, but it remains a mystery how cartilage itself first arose in evolution. Characterization of cartilage in cuttlefish and horseshoe crabs reveals surprising commonalities with chordate chondrocytes, suggesting a common evolutionary origin.

Developmental basis of phenotypic integration in two Lake Malawi cichlids

Background

Cichlid fishes from the Rift Lakes of East Africa have undergone the most spectacular adaptive radiations in vertebrate history. Eco-morphological adaptations in lakes Victoria, Malawi and Tanganyika have resulted in a vast array of skull shapes and sizes, yet primary axes of morphological variation are conserved in all three radiations, prominently including the size of the preorbital region of the skull. This conserved pattern suggests that development may constrain the trajectories of cichlid head morphological evolution.

Results

Here, we (1) present a comparative analysis of adult head morphology in two sand-dweller cichlids from Lake Malawi with preorbital size differences representative of the main axis of variation among the three lakes and (2) analyze the ontogeny of shape and size differences by focusing on known developmental modules throughout the head. We find that (1) developmental differences between the two species correlate with known developmental modules; (2) differences in embryonic cartilage development result in phenotypically integrated changes among all bones derived from a single cartilage, while differences in dermal bone development tend to influence isolated regions within a bone; and lastly (3) species-specific morphologies appear in the embryo as subtle differences, which become progressively amplified throughout ontogeny. We propose that this amplification takes place at skeletal growth zones, the locations and shapes of which are patterned during embryogenesis.

Conclusions

This study is the most anatomically comprehensive analysis of the developmental differences underlying cichlid skull evolution in the Rift Lakes of East Africa. The scale of our analysis reveals previously unnoticed correlations between developmental modules and patterns of phenotypic integration. We propose that the primary axes of morphological variation among East African cichlid adaptive radiations are constrained by the hierarchical modularity of the teleost head skeleton.

On the evolutionary relationship between chondrocytes and osteoblasts

Vertebrates are the only animals that produce bone, but the molecular genetic basis for this evolutionary novelty remains obscure. Here, we synthesize information from traditional evolutionary and modern molecular genetic studies in order to generate a working hypothesis on the evolution of the gene regulatory network (GRN) underlying bone formation. Since transcription factors are often core components of GRNs (i.e., kernels), we focus our analyses on Sox9 and Runx2. Our argument centers on three skeletal tissues that comprise the majority of the vertebrate skeleton: immature cartilage, mature cartilage, and bone. Immature cartilage is produced during early stages of cartilage differentiation and can persist into adulthood, whereas mature cartilage undergoes additional stages of differentiation, including hypertrophy and mineralization. Functionally, histologically, and embryologically, these three skeletal tissues are very similar, yet unique, suggesting that one might have evolved from another. Traditional studies of the fossil record, comparative anatomy and embryology demonstrate clearly that immature cartilage evolved before mature cartilage or bone. Modern molecular approaches show that the GRNs regulating differentiation of these three skeletal cell fates are similar, yet unique, just like the functional and histological features of the tissues themselves. Intriguingly, the Sox9 GRN driving cartilage formation appears to be dominant to the Runx2 GRN of bone. Emphasizing an embryological and evolutionary transcriptomic view, we hypothesize that the Runx2 GRN underlying bone formation was co-opted from mature cartilage. We discuss how modern molecular genetic experiments, such as comparative transcriptomics, can test this hypothesis directly, meanwhile permitting levels of constraint and adaptation to be evaluated quantitatively. Therefore, comparative transcriptomics may revolutionize understanding of not only the clade-specific evolution of skeletal cells, but also the generation of evolutionary novelties, providing a modern paradigm for the evolutionary process.

Role of RUNX2 in Breast Carcinogenesis

RUNX2 is a transcription factor playing the major role in osteogenesis, but it can be involved in DNA damage response, which is crucial for cancer transformation. RUNX2 can interact with cell cycle regulators: cyclin-dependent kinases, pRB and p21Cip1 proteins, as well as the master regulator of the cell cycle, the p53 tumor suppressor. RUNX2 is involved in many signaling pathways, including those important for estrogen signaling, which, in turn, are significant for breast carcinogenesis. RUNX2 can promote breast cancer development through Wnt and Tgfβ signaling pathways, especially in estrogen receptor (ER)-negative cases. ERα interacts directly with RUNX2 and regulates its activity. Moreover, the ERα gene has a RUNX2 binding site within its promoter. RUNX2 stimulates the expression of aromatase, an estrogen producing enzyme, increasing the level of estrogens, which in turn stimulate cell proliferation and replication errors, which can be turned into carcinogenic mutations. Exploring the role of RUNX2 in the pathogenesis of breast cancer can lead to revealing new therapeutic targets.

Characterization of the Runx gene family in a jawless vertebrate, the Japanese lamprey (Lethenteron japonicum)

The cyclostomes (jawless vertebrates), comprising lampreys and hagfishes, are the sister group of jawed vertebrates (gnathostomes) and are hence an important group for the study of vertebrate evolution. In mammals, three Runx genes, Runx1, Runx2 and Runx3, encode transcription factors that are essential for cell proliferation and differentiation in major developmental pathways such as haematopoiesis, skeletogenesis and neurogenesis and are frequently associated with diseases. We describe here the characterization of Runx gene family members from a cyclostome, the Japanese lamprey (Lethenteron japonicum). The Japanese lamprey contains three Runx genes, RunxA, RunxB, and RunxC. However, phylogenetic and synteny analyses suggest that they are not one-to-one orthologs of gnathostome Runx1, Runx2 and Runx3. The major protein domains and motifs found in gnathostome Runx proteins are highly conserved in the lamprey Runx proteins. Although all gnathostome Runx genes each contain two alternative promoters, P1 (distal) and P2 (proximal), only lamprey RunxB possesses the alternative promoters; lamprey RunxA and RunxC contain only P2 and P1 promoter, respectively. Furthermore, the three lamprey Runx genes give rise to fewer alternative isoforms than the three gnathostome Runx genes. The promoters of the lamprey Runx genes lack the tandem Runx-binding motifs that are highly conserved among the P1 promoters of gnathostome Runx1, Runx2 and Runx3 genes; instead these promoters contain dispersed single Runx-binding motifs. The 3'UTR of lamprey RunxB contains binding sites for miR-27 and miR-130b/301ab, which are conserved in mammalian Runx1 and Runx3, respectively. Overall, the Runx genes in lamprey seem to have experienced a different evolutionary trajectory from that of gnathostome Runx genes which are highly conserved all the way from cartilaginous fishes to mammals.

Runx family genes in a cartilaginous fish, the elephant shark (Callorhinchus milii)

The Runx family genes encode transcription factors that play key roles in hematopoiesis, skeletogenesis and neurogenesis and are often implicated in diseases. We describe here the cloning and characterization of Runx1, Runx2, Runx3 and Runxb genes in the elephant shark (Callorhinchus milii), a member of Chondrichthyes, the oldest living group of jawed vertebrates. Through the use of alternative promoters and/or alternative splicing, each of the elephant shark Runx genes expresses multiple isoforms similar to their orthologs in human and other bony vertebrates. The expression profiles of elephant shark Runx genes are similar to those of mammalian Runx genes. The syntenic blocks of genes at the elephant shark Runx gene loci are highly conserved in human, but represented by shorter conserved blocks in zebrafish indicating a higher degree of rearrangements in this teleost fish. Analysis of promoter regions revealed conservation of binding sites for transcription factors, including two tandem binding sites for Runx that are totally conserved in the distal promoter regions of elephant shark Runx1-3. Several conserved noncoding elements (CNEs), which are putative cis-regulatory elements, and miRNA binding sites were identified in the elephant shark and human Runx gene loci. Some of these CNEs and miRNA binding sites are absent in teleost fishes such as zebrafish and fugu. In summary, our analysis reveals that the genomic organization and expression profiles of Runx genes were already complex in the common ancestor of jawed vertebrates.

The vertebrate RCAN gene family: novel insights into evolution, structure and regulation

Recently there has been much interest in the Regulators of Calcineurin (RCAN) proteins which are important endogenous modulators of the calcineurin-NFATc signalling pathway. They have been shown to have a crucial role in cellular programmes such as the immune response, muscle fibre remodelling and memory, but also in pathological processes such as cardiac hypertrophy and neurodegenerative diseases. In vertebrates, the RCAN family form a functional subfamily of three members RCAN1, RCAN2 and RCAN3 whereas only one RCAN is present in the rest of Eukarya. In addition, RCAN genes have been shown to collocate with RUNX and CLIC genes in ACD clusters (ACD21, ACD6 and ACD1). How the RCAN genes and their clustering in ACDs evolved is still unknown. After analysing RCAN gene family evolution using bioinformatic tools, we propose that the three RCAN vertebrate genes within the ACD clusters, which evolved from single copy genes present in invertebrates and lower eukaryotes, are the result of two rounds of whole genome duplication, followed by a segmental duplication. This evolutionary scenario involves the loss or gain of some RCAN genes during evolution. In addition, we have analysed RCAN gene structure and identified the existence of several characteristic features that can be involved in RCAN evolution and gene expression regulation. These included: several transposable elements, CpG islands in the 5′ region of the genes, the existence of antisense transcripts (NAT) associated with the three human genes, and considerable evidence for bidirectional promoters that regulate RCAN gene expression. Furthermore, we show that the CpG island associated with the RCAN3 gene promoter is unmethylated and transcriptionally active. All these results provide timely new insights into the molecular mechanisms underlying RCAN function and a more in depth knowledge of this gene family whose members are obvious candidates for the development of future therapies.

Distinct global shifts in genomic binding profiles of limb malformation-associated HOXD13 mutations

Gene regulation by transcription factors (TFs) determines developmental programs and cell identity. Consequently, mutations in TFs can lead to dramatic phenotypes in humans by disrupting gene regulation. To date, the molecular mechanisms that actually cause these phenotypes have been difficult to address experimentally. ChIP-seq, which couples chromatin immunoprecipitation with high-throughput sequencing, allows TF function to be investigated on a genome-wide scale, enabling new approaches for the investigation of gene regulation. Here, we present the application of ChIP-seq to explore the effect of missense mutations in TFs on their genome-wide binding profile. Using a retroviral expression system in chicken mesenchymal stem cells, we elucidated the mechanism underlying a novel missense mutation in HOXD13 (Q317K) associated with a complex hand and foot malformation phenotype. The mutated glutamine (Q) is conserved in most homeodomains, a notable exception being bicoid-type homeodomains that have lysine (K) at this position. Our results show that the mutation results in a shift in the binding profile of the mutant toward a bicoid/PITX1 motif. Gene expression analysis and functional assays using in vivo overexpression studies confirm that the mutation results in a partial conversion of HOXD13 into a TF with bicoid/PITX1 properties. A similar shift was not observed with another mutation, Q317R, which is associated with brachysyndactyly, suggesting that the bicoid/PITX1-shift observed for Q317K might be related to the severe clinical phenotype. The methodology described can be used to investigate a wide spectrum of TFs and mutations that have not previously been amenable to ChIP-seq experiments.

Morpholino oligonucleotide knockdown of the extracellular calcium-sensing receptor impairs early skeletal development in zebrafish

The complex vertebrate skeleton depends on regulated cell activities to lay down protein matrix and mineral components of bone. As a distinctive vertebrate characteristic, bone is a storage site for physiologically-important calcium ion. The extracellular calcium-sensing receptor (CaSR) is linked to homeostatic regulation of calcium through its expression in endocrine glands that secrete calcium homeostatic hormones, in Ca(2+)- and ion-transporting epithelia, and in skeleton. Since CaSR is restricted in its presence to the chordate-vertebrate evolutionary lineage, we propose there to be important functional ties between CaSRs and vertebrate skeleton in the context of that group's characteristic form of calcium-mineralized skeleton. Since little is known about CaSR in the skeletal biology of non-mammalian vertebrates, reverse transcription-polymerase chain reaction (RT-PCR), in situ hybridization and immunohistochemistry were applied to adult and embryonic zebrafish to reveal CaSR transcript and protein expression in several tissues, including, among these, chondrocytes and developing bone and notochord as components in skeletal development. Morpholino oligonucleotide (MO) knockdown technique was used to probe CaSR role(s) in the zebrafish model system. By RT-PCR assessment, injection of a splice-inhibiting CaSR MO reduced normally-spliced Casr gene transcript expression measured at 2days postfertilization (dpf). Corresponding to the knockdown of normally-spliced mRNA by the CaSR MO, we observed a morphant phenotype characterized by stunted growth and disorganization of the notochord and axial skeleton by 1dpf. We conclude that, like its critically important role in normal bone development in mammals, CaSR is essential in skeletogenesis in fishes.

Expression of the amelogenin gene in the skin of Xenopus tropicalis

Anuran skin contains a calcified dermal layer, referred to as the Eberth-Kastschenko (EK) layer, which is found between the stratum spongiosum and the stratum compactum. Although it is established that some anuran species possess the EK layer, little is known about this layer from the standpoint of evolutionary and developmental biology. We conducted a morphological analysis by staining the dorsal skin from many species with alizarin red S to investigate the calcified layer. This layer was observed in all of the anurans tested, as well as in fishes and one species of caecilian with dermal scales, but not in urodeles, amniotes, or a scaleless caecilian. All of the investigated species with dermal scales exhibited a calcified layer in their dermis, while the anurans showed the EK layer, but no scales. We also analyzed the expression of genes related to scale formation (sparc, mmp9, and mmp2) in the dorsal skin of X. tropicalis. These genes were highly expressed at the metamorphic climax stage, which preceded the deposition of calcium. Furthermore, we examined the gene expression profile of amelogenin, the major protein found in the enamel matrix of the developing tooth. In X. tropicalis, amelogenin was upregulated in the skin at the climax stage and was expressed in the adult dermis at a high level. These data provide the first experimental evidence of the expression of amelogenin in the skin. These findings will lead to a better understanding of the developmental formation of the EK layer and the function of amelogenin.

Progenitors of the protochordate ocellus as an evolutionary origin of the neural crest

The neural crest represents a highly multipotent population of embryonic stem cells found only in vertebrate embryos. Acquisition of the neural crest during the evolution of vertebrates was a great advantage, providing Chordata animals with the first cellular cartilage, bone, dentition, advanced nervous system and other innovations. Today not much is known about the evolutionary origin of neural crest cells. Here we propose a novel scenario in which the neural crest originates from neuroectodermal progenitors of the pigmented ocelli in Amphioxus-like animals. We suggest that because of changes in photoreception needs, these multipotent progenitors of photoreceptors gained the ability to migrate outside of the central nervous system and subsequently started to give rise to neural, glial and pigmented progeny at the periphery.

Putative chitin synthases from Branchiostoma floridae show extracellular matrix-related domains and mosaic structures

The transition from unicellular to multicellular life forms requires the development of a specialized structural component, the extracellular matrix (ECM). In Metazoans, there are two main supportive systems, which are based on chitin and collagen/hyaluronan, respectively. Chitin is the major constituent of fungal cell walls and arthropod exoskeleton. However, presence of chitin/chitooligosaccharides has been reported in lower chordates and during specific stages of vertebrate development. In this study, the occurrence of chitin synthases (CHSs) was investigated with a bioinformatics approach in the cephalochordate Branchiostoma floridae, in which the presence of chitin was initially reported in the skeletal rods of the pharyngeal gill basket. Twelve genes coding for proteins containing conserved amino acid residues of processive glycosyltransferases from GT2 family were found and 10 of them display mosaic structures with novel domains never reported previously in a chitin synthase. In particular, the presence of a discoidin (DS) and a sterile alpha motif (SAM) domain was found in nine identified proteins. Sequence analyses and homology modelling suggest that these domains might interact with the extracellular matrix and mediate protein–protein interaction. The multi-domain putative chitin synthases from B. floridae constitute an emblematic example of the explosion of domain innovation and shuffling which predate Metazoans.

Evolutionary crossroads in developmental biology: cyclostomes (lamprey and hagfish)

Lampreys and hagfish, which together are known as the cyclostomes or ‘agnathans’, are the only surviving lineages of jawless fish. They diverged early in vertebrate evolution, before the origin of the hinged jaws that are characteristic of gnathostome (jawed) vertebrates and before the evolution of paired appendages. However, they do share numerous characteristics with jawed vertebrates. Studies of cyclostome development can thus help us to understand when, and how, key aspects of the vertebrate body evolved. Here, we summarise the development of cyclostomes, highlighting the key species studied and experimental methods available. We then discuss how studies of cyclostomes have provided important insight into the evolution of fins, jaws, skeleton and neural crest.

Evolution of the bone gene regulatory network

Current fossil, embryological and genetic data shed light on the evolution of the gene regulatory network (GRN) governing bone formation. The key proteins and genes involved in skeletogenesis are well accepted. We discuss when these essential components of the GRN evolved and propose that the Runx genes, master regulators of skeletogenesis, functioned in early cartilages well before they were co-opted to function in the making of bone. Two rounds of whole genome duplication, together with additional tandem gene duplications, created a genetic substrate for segregation of one GRN into several networks regulating the related tissues of cartilage, bone, enamel, and dentin. During this segregation, Runx2 assumed its position at the top of the bone GRN, and Sox9 was excluded from bone, retaining its ancient role in cartilage.

Evolution of the osteoblast: skeletogenesis in gar and zebrafish

BACKGROUND

Although the vertebrate skeleton arose in the sea 500 million years ago, our understanding of the molecular fingerprints of chondrocytes and osteoblasts may be biased because it is informed mainly by research on land animals. In fact, the molecular fingerprint of teleost osteoblasts differs in key ways from that of tetrapods, but we do not know the origin of these novel gene functions. They either arose as neofunctionalization events after the teleost genome duplication (TGD), or they represent preserved ancestral functions that pre-date the TGD. Here, we provide evolutionary perspective to the molecular fingerprints of skeletal cells and assess the role of genome duplication in generating novel gene functions. We compared the molecular fingerprints of skeletogenic cells in two ray-finned fish: zebrafish (Danio rerio)--a teleost--and the spotted gar (Lepisosteus oculatus)--a "living fossil" representative of a lineage that diverged from the teleost lineage prior to the TGD (i.e., the teleost sister group). We analyzed developing embryos for expression of the structural collagen genes col1a2, col2a1, col10a1, and col11a2 in well-formed cartilage and bone, and studied expression of skeletal regulators, including the transcription factor genes sox9 and runx2, during mesenchymal condensation.

RESULTS

Results provided no evidence for the evolution of novel functions among gene duplicates in zebrafish compared to the gar outgroup, but our findings shed light on the evolution of the osteoblast. Zebrafish and gar chondrocytes both expressed col10a1 as they matured, but both species' osteoblasts also expressed col10a1, which tetrapod osteoblasts do not express. This novel finding, along with sox9 and col2a1 expression in developing osteoblasts of both zebrafish and gar, demonstrates that osteoblasts of both a teleost and a basally diverging ray-fin fish express components of the supposed chondrocyte molecular fingerprint.

CONCLUSIONS

Our surprising finding that the "chondrogenic" transcription factor sox9 is expressed in developing osteoblasts of both zebrafish and gar can help explain the expression of chondrocyte genes in osteoblasts of ray-finned fish. More broadly, our data suggest that the molecular fingerprint of the osteoblast, which largely is constrained among land animals, was not fixed during early vertebrate evolution.

Gene regulatory network reconstruction using Bayesian networks, the Dantzig Selector, the Lasso and their meta-analysis

Modern technologies and especially next generation sequencing facilities are giving a cheaper access to genotype and genomic data measured on the same sample at once. This creates an ideal situation for multifactorial experiments designed to infer gene regulatory networks. The fifth "Dialogue for Reverse Engineering Assessments and Methods" (DREAM5) challenges are aimed at assessing methods and associated algorithms devoted to the inference of biological networks. Challenge 3 on "Systems Genetics" proposed to infer causal gene regulatory networks from different genetical genomics data sets. We investigated a wide panel of methods ranging from Bayesian networks to penalised linear regressions to analyse such data, and proposed a simple yet very powerful meta-analysis, which combines these inference methods. We present results of the Challenge as well as more in-depth analysis of predicted networks in terms of structure and reliability. The developed meta-analysis was ranked first among the 16 teams participating in Challenge 3A. It paves the way for future extensions of our inference method and more accurate gene network estimates in the context of genetical genomics.

The homology of odontodes in gnathostomes: insights from Dlx gene expression in the dogfish, Scyliorhinus canicula

BACKGROUND

Teeth and tooth-like structures, together named odontodes, are repeated organs thought to share a common evolutionary origin. These structures can be found in gnathostomes at different locations along the body: oral teeth in the jaws, teeth and denticles in the oral-pharyngeal cavity, and dermal denticles on elasmobranch skin. We, and other colleagues, had previously shown that teeth in any location were serially homologous because: i) pharyngeal and oral teeth develop through a common developmental module; and ii) the expression patterns of the Dlx genes during odontogenesis were highly divergent between species but almost identical between oral and pharyngeal dentitions within the same species. Here we examine Dlx gene expression in oral teeth and dermal denticles in order to test the hypothesis of serial homology between these odontodes.

RESULTS

We present a detailed comparison of the first developing teeth and dermal denticles (caudal primary scales) of the dogfish (Scyliorhinus canicula) and show that both odontodes develop through identical stages that correspond to the common stages of oral and pharyngeal odontogenesis. We identified six Dlx paralogs in the dogfish and found that three showed strong transcription in teeth and dermal denticles (Dlx3, Dlx4 and Dlx5) whereas a weak expression was detected for Dlx1 in dermal denticles and teeth, and for Dlx2 in dermal denticles. Very few differences in Dlx expression patterns could be detected between tooth and dermal denticle development, except for the absence of Dlx2 expression in teeth.

CONCLUSIONS

Taken together, our histological and expression data strongly suggest that teeth and dermal denticles develop from the same developmental module and under the control of the same set of Dlx genes. Teeth and dermal denticles should therefore be considered as serial homologs developing through the initiation of a common gene regulatory network (GRN) at several body locations. This mechanism of heterotopy supports the 'inside and out' model that has been recently proposed for odontode evolution.

SoxE gene duplication and development of the lamprey branchial skeleton: Insights into development and evolution of the neural crest

SoxE genes are multifunctional transcriptional regulators that play key roles in specification and differentiation of neural crest. Three members (Sox8, Sox9, Sox10) are expressed in the neural crest and are thought to modulate the expression and activity of each other. In addition to regulating the expression of other early neural crest marker genes, SoxE genes are required for development of cartilage. Here we investigated the role of SoxE genes in development of the neural crest-derived branchial skeleton in the sea lamprey. Using a morpholino knockdown approach, we show that all three SoxE genes described in lamprey are required for branchial basket development. Our results suggest that SoxE1 and SoxE2 are required for specification of the chondrogenic neural crest. SoxE3 plays a morphogenetic role in patterning of the branchial basket and may be required for the development of mucocartilage, a tissue unique to larval lampreys. While the lamprey branchial basket develops primarily from an elastin-like major extracellular matrix protein that is specific to lampreys, fibrillar collagen is also expressed in developing branchial cartilage and may be regulated by the lamprey SoxE genes. Our data suggest that the regulation of Type II collagen by Sox9 might have been co-opted by the neural crest in development of the branchial skeleton following the divergence of agnathan and gnathostome vertebrates. Finally, our results also have implications for understanding the independent evolution of duplicated SoxE genes among agnathan and gnathostome vertebrates.

The evolution of extracellular matrix

We present a perspective on the molecular evolution of the extracellular matrix (ECM) in metazoa that draws on research publications and data from sequenced genomes and expressed sequence tag libraries. ECM components do not function in isolation, and the biological ECM system or “adhesome” also depends on posttranslational processing enzymes, cell surface receptors, and extracellular proteases. We focus principally on the adhesome of internal tissues and discuss its origins at the dawn of the metazoa and the expansion of complexity that occurred in the chordate lineage. The analyses demonstrate very high conservation of a core adhesome that apparently evolved in a major wave of innovation in conjunction with the origin of metazoa. Integrin, CD36, and certain domains predate the metazoa, and some ECM-related proteins are identified in choanoflagellates as predicted sequences. Modern deuterostomes and vertebrates have many novelties and elaborations of ECM as a result of domain shuffling, domain innovations and gene family expansions. Knowledge of the evolution of metazoan ECM is important for understanding how it is built as a system, its roles in normal tissues and disease processes, and has relevance for tissue engineering, the development of artificial organs, and the goals of synthetic biology.

Where did bone come from?

Bone is specific to vertebrates, and originated as mineralization around the basal membrane of the throat or skin, giving rise to tooth-like structures and protective shields in animals with a soft cartilage-like endoskeleton. A combination of fossil anatomy and genetic information from modern species has improved our understanding of the evolution of bone. Thus, even in man, there are still similarities in the molecular regulation of skin appendages and bone. This article gives a brief overview of the major milestones in skeletal evolution. Some molecular machineries involving members of core genetic networks and their interactions are described in the context of both old theories and modern genetic approaches.

A new mechanistic scenario for the origin and evolution of vertebrate cartilage

The appearance of cellular cartilage was a defining event in vertebrate evolution because it made possible the physical expansion of the vertebrate "new head". Despite its central role in vertebrate evolution, the origin of cellular cartilage has been difficult to understand. This is largely due to a lack of informative evolutionary intermediates linking vertebrate cellular cartilage to the acellular cartilage of invertebrate chordates. The basal jawless vertebrate, lamprey, has long been considered key to understanding the evolution of vertebrate cartilage. However, histological analyses of the lamprey head skeleton suggest it is composed of modern cellular cartilage and a putatively unrelated connective tissue called mucocartilage, with no obvious transitional tissue. Here we take a molecular approach to better understand the evolutionary relationships between lamprey cellular cartilage, gnathostome cellular cartilage, and lamprey mucocartilage. We find that despite overt histological similarity, lamprey and gnathostome cellular cartilage utilize divergent gene regulatory networks (GRNs). While the gnathostome cellular cartilage GRN broadly incorporates Runx, Barx, and Alx transcription factors, lamprey cellular cartilage does not express Runx or Barx, and only deploys Alx genes in certain regions. Furthermore, we find that lamprey mucocartilage, despite its distinctive mesenchymal morphology, deploys every component of the gnathostome cartilage GRN, albeit in different domains. Based on these findings, and previous work, we propose a stepwise model for the evolution of vertebrate cellular cartilage in which the appearance of a generic neural crest-derived skeletal tissue was followed by a phase of skeletal tissue diversification in early agnathans. In the gnathostome lineage, a single type of rigid cellular cartilage became dominant, replacing other skeletal tissues and evolving via gene cooption to become the definitive cellular cartilage of modern jawed vertebrates.

Gain, loss and divergence in primate zinc-finger genes: a rich resource for evolution of gene regulatory differences between species

The molecular changes underlying major phenotypic differences between humans and other primates are not well understood, but alterations in gene regulation are likely to play a major role. Here we performed a thorough evolutionary analysis of the largest family of primate transcription factors, the Krüppel-type zinc finger (KZNF) gene family. We identified and curated gene and pseudogene models for KZNFs in three primate species, chimpanzee, orangutan and rhesus macaque, to allow for a comparison with the curated set of human KZNFs. We show that the recent evolutionary history of primate KZNFs has been complex, including many lineage-specific duplications and deletions. We found 213 species-specific KZNFs, among them 7 human-specific and 23 chimpanzee-specific genes. Two human-specific genes were validated experimentally. Ten genes have been lost in humans and 13 in chimpanzees, either through deletion or pseudogenization. We also identified 30 KZNF orthologs with human-specific and 42 with chimpanzee-specific sequence changes that are predicted to affect DNA binding properties of the proteins. Eleven of these genes show signatures of accelerated evolution, suggesting positive selection between humans and chimpanzees. During primate evolution the most extensive re-shaping of the KZNF repertoire, including most gene additions, pseudogenizations, and structural changes occurred within the subfamily homininae. Using zinc finger (ZNF) binding predictions, we suggest potential impact these changes have had on human gene regulatory networks. The large species differences in this family of TFs stands in stark contrast to the overall high conservation of primate genomes and potentially represents a potent driver of primate evolution.

Tenascins and the importance of adhesion modulation

Tenascins are a family of extracellular matrix proteins that evolved in early chordates. There are four family members: tenascin-X, tenascin-R, tenascin-W, and tenascin-C. Tenascin-X associates with type I collagen, and its absence can cause Ehlers-Danlos Syndrome. In contrast, tenascin-R is concentrated in perineuronal nets. The expression of tenascin-C and tenascin-W is developmentally regulated, and both are expressed during disease (e.g., both are associated with cancer stroma and tumor blood vessels). In addition, tenascin-C is highly induced by infections and inflammation. Accordingly, the tenascin-C knockout mouse has a reduced inflammatory response. All tenascins have the potential to modify cell adhesion either directly or through interaction with fibronectin, and cell-tenascin interactions typically lead to increased cell motility. In the case of tenascin-C, there is a correlation between elevated expression and increased metastasis in several types of tumors.

Evolution of the parathyroid hormone family and skeletal formation pathways

Bone is considered to be a feature of higher vertebrates and one of the features that was required for the movement from water onto land. But there are a number of evolutionarily important species that have cartilaginous skeletons, including sharks. Both bony and cartilaginous fish are believed to have a common ancestor who had a bony skeleton. A number of factors and pathways have been shown to be involved in the development and maintenance of bony skeleton including the Wnt pathway and the parathyroid hormone gene family. The study of these pathways and factors in cartilaginous animals may shed light on the evolution of the vertebrate skeleton.

Evolution of developmental pattern for vertebrate dentitions: an oro-pharyngeal specific mechanism

Classically the oral dentition with teeth regulated into a successional iterative order was thought to have evolved from the superficial skin denticles migrating into the mouth at the stage when jaws evolved. The canonical view is that the initiation of a pattern order for teeth at the mouth margin required development of a sub-epithelial, permanent dental lamina. This provided regulated tooth production in advance of functional need, as exemplified by the Chondrichthyes. It had been assumed that teeth in the Osteichthyes form in this way as in tetrapods. However, this has been shown not to be true for many osteichthyan fish where a dental lamina of this kind does not form, but teeth are regularly patterned and replaced. We question the evolutionary origin of pattern information for the dentition driven by new morphological data on spatial initiation of skin denticles in the catshark. We review recent gene expression data for spatio-temporal order of tooth initiation for Scyliorhinus canicula, selected teleosts in both oral and pharyngeal dentitions, and Neoceratodus forsteri. Although denticles in the chondrichthyan skin appear not to follow a strict pattern order in space and time, tooth replacement in a functional system occurs with precise timing and spatial order. We suggest that the patterning mechanism observed for the oral and pharyngeal dentition is unique to the vertebrate oro-pharynx and independent of the skin system. Therefore, co-option of a successional iterative pattern occurred in evolution not from the skin but from mechanisms existing in the oro-pharynx of now extinct agnathans.

Two families of Xenopus tropicalis skeletal genes display well-conserved expression patterns with mammals in spite of their highly divergent regulatory regions

The origin of bone and cartilage, and their subsequent diversification in specific vertebrate lineages, is intimately linked to the precise transcriptional control of genes involved in matrix mineralization. It is not yet clear, however, to which extent the osteoblasts, osteocytes, and chondrocytes of each of the major vertebrate groups express similar sets of genes. In this study we have focused on the evolution of two independent families of genes that code for extracellular matrix components of the skeleton and that include secreted protein, acidic, cysteine-rich (SPARC), bone sialoprotein (BSP) and dentin matrix protein 1 (DMP1) paralogues, and the osteocalcin (OC) and matrix gla protein (MGP) paralogues. Analyzing developing Xenopus tropicalis skeletal elements, we show that the expression patterns of these genes are well conserved with mammals. The fact that only a few osteoblasts express DMP1, while only some osteocytes express SPARC and BSP, reveals a significant degree of molecular heterogeneity for these two populations of X. tropicalis cells, similarly to what has been described in mouse. Although the cis-regulatory modules (CRM) of the mammalian OC, DMP1, and BSP orthologs have been functionally characterized, we found no evidence of sequence similarity between these regions and the X. tropicalis genome. Furthermore, these regulatory elements evolve rapidly, as they are only poorly conserved between human and rodents. Therefore, the SPARC/DMP1/BSP and the OC/MGP families provide a good paradigm to study how transcriptional output can be maintained in skeletal cells despite extensive sequence divergence of CRM.

Origin and genetic evolution of the vertebrate skeleton

The current understanding of the origin and evolution of the genetic cassette for the vertebrate skeletal system is reviewed. Molecular phylogenetic analyses of fibrillar collagen genes, which encode the main component of both cartilage and mineralized bone, suggest that genome duplications in vertebrate ancestors were essential for producing distinct collagen fibers for cartilage and mineralized bone. Several data Indicate co-expression of the ancestral copy of fibrillar collagen with the SoxE and Runx transcription factors. Therefore, the genetic cassette may have already existed in protochordate ancestors, and may operate in the development of the pharyngeal gill skeleton. Accompanied by genome duplications in vertebrate ancestors, this genetic cassette may have also been duplicated and co-opted for cartilage and bone. Subsequently, the genetic cassette for cartilage recruited novel genetic material via domain shuffling. Aggrecan, acquired by means of domain shuffling, performs an essential role in cartilage as a shock absorber. In contrast, the cassette for bone recruited new genetic material produced by tandem duplication of the SPARC/osteonectin genes. Some of the duplicated copies of SPARC/osteonectin became secretory Cabinding phosphoproteins (SCPPs) performing a central role in mineralization by regulating the calcium phosphate concentration. Comparative genome analysis revealed similar molecular evolutionary histories for the genetic cassettes for cartilage and bone, namely duplication of the ancestral genetic cassette and recruitment of novel genetic material.

Evolution of the interaction between Runx2 and VDR, two transcription factors involved in osteoblastogenesis

Background

The mineralized skeleton is a major evolutionary novelty that has contributed to the impressive morphological diversifications of the vertebrates. Essential to bone biology is the solidified extracellular matrix secreted by highly specialized cells, the osteoblasts. We now have a rather complete view of the events underlying osteogenesis, from a cellular, molecular, genetic, and epigenetic perspective. Because this knowledge is still largely restricted to mammals, it is difficult, if not impossible, to deduce the evolutionary history of the regulatory network involved in osteoblasts specification and differentiation. In this study, we focused on the transcriptional regulators Runx2 and VDR (the Vitamin D Receptor) that, in mammals, directly interact together and stabilize complexes of co-activators and chromatin remodellers, thereby allowing the transcriptional activation of target genes involved in extracellular matrix mineralization. Using a combination of functional, biochemical, and histological approaches, we have asked if the interaction observed between Runx2 and VDR represents a recent mammalian innovation, or if it results from more ancient changes that have occurred deep in the vertebrate lineage.

Results

Using immunohistochemistry and in situ hybridization in developing embryos of chick, frog and teleost fishes, we have revealed that the co-expression of Runx2 and VDR in skeletal elements has been particularly strengthened in the lineage leading to amniotes. We show that the teleost Runx2 orthologue as well as the three mammalian Runx1, Runx2 and Runx3 paralogues are able to co-immunoprecipitate with the VDR protein present in nuclear extracts of rat osteoblasts stimulated with 1α,25-dihydroxyvitamin D₃. In addition, the teleost Runx2 can activate the transcription of the mammalian osteocalcin promoter in transfection experiments, and this response can be further enhanced by 1α,25-dihydroxyvitamin D₃. Finally, using pull-down experiments between recombinant proteins, we show that the VDR homologue from teleosts, but not from ascidians, is able to directly interact with the mammalian Runx2 homologue.

Conclusions

We propose an evolutionary scenario for the assembly of the molecular machinery involving Runx2 and VDR in vertebrates. In the last common ancestor of actinopterygians and sacropterygians, the three Runx paralogues possessed the potential to physically and functionally interact with the VDR protein. Therefore, 1α,25-dihydroxyvitamin D₃ might have been able to modulate the transcriptional activity of Runx1, Runx2 or Runx3 in the tissues expressing VDR. After the split from amphibians, in the lineage leading to amniotes, Runx2 and VDR became robustly co-expressed in developing skeletal elements, and their regulatory interaction was incorporated in the genetic program involved in the specification and differentiation of osteoblasts.

Lineage-specific transcription factors and the evolution of gene regulatory networks

Nature is replete with examples of diverse cell types, tissues and body plans, forming very different creatures from genomes with similar gene complements. However, while the genes and the structures of proteins they encode can be highly conserved, the production of those proteins in specific cell types and at specific developmental time points might differ considerably between species. A full understanding of the factors that orchestrate gene expression will be essential to fully understand evolutionary variety. Transcription factor (TF) proteins, which form gene regulatory networks (GRNs) to act in cooperative or competitive partnerships to regulate gene expression, are key components of these unique regulatory programs. Although many TFs are conserved in structure and function, certain classes of TFs display extensive levels of species diversity. In this review, we highlight families of TFs that have expanded through gene duplication events to create species-unique repertoires in different evolutionary lineages. We discuss how the hierarchical structures of GRNs allow for flexible small to large-scale phenotypic changes. We survey evidence that explains how newly evolved TFs may be integrated into an existing GRN and how molecular changes in TFs might impact the GRNs. Finally, we review examples of traits that evolved due to lineage-specific TFs and species differences in GRNs.

Pigmentation pathway evolution after whole-genome duplication in fish

Whole-genome duplications (WGDs) have occurred repeatedly in the vertebrate lineage, but their evolutionary significance for phenotypic evolution remains elusive. Here, we have investigated the impact of the fish-specific genome duplication (FSGD) on the evolution of pigmentation pathways in teleost fishes. Pigmentation and color patterning are among the most diverse traits in teleosts, and their pigmentary system is the most complex of all vertebrate groups.

Using a comparative genomic approach including phylogenetic and synteny analyses, the evolution of 128 vertebrate pigmentation genes in five teleost genomes following the FSGD has been reconstructed. We show that pigmentation genes have been preferentially retained in duplicate after the FSGD, so that teleosts have 30% more pigmentation genes compared with tetrapods. This is significantly higher than genome-wide estimates of FSGD gene duplicate retention in teleosts. Large parts of the melanocyte regulatory network have been retained in two copies after the FSGD. Duplicated pigmentation genes follow general evolutionary patterns such as the preservation of protein complex stoichiometries and the overrepresentation of developmental genes among retained duplicates. These results suggest that the FSGD has made an important contribution to the evolution of teleost-specific features of pigmentation, which include novel pigment cell types or the division of existing pigment cell types into distinct subtypes. Furthermore, we have observed species-specific differences in duplicate retention and evolution that might contribute to pigmentary diversity among teleosts.

Our study therefore strongly supports the hypothesis that WGDs have promoted the increase of complexity and diversity during vertebrate phenotypic evolution.

Expression of Sox and fibrillar collagen genes in lamprey larval chondrogenesis with implications for the evolution of vertebrate cartilage

Lampreys possess unique types of cartilage in which elastin-like proteins are the dominant matrix component, whereas gnathostome cartilage is mainly composed of fibrillar collagen. Despite the differences in protein composition, the Sox-col2a1 genetic cascade was suggested to be conserved between lamprey pharyngeal cartilage and gnathostome cartilage. We examined whether the cascade is conserved in another type of lamprey cartilage, the trabecular cartilage. We found that SoxD and SoxE are expressed in both trabecular and pharyngeal cartilages. However, trabecular cartilage shows no clade A fibrillar collagen gene expression, including genes expressed in pharyngeal cartilage of this animal. On the basis of these observations, we propose that lampreys possess an ancestral type of cartilage that is similar to amphioxus gill cartilage, and in this respect, gnathostome cartilage can be regarded as derived for the loss of elastin-like protein as a cartilage component and recruitment of fibrillar collagen, which is included as a minor component in the ancestral cartilage, as the main component.