Divergence of Hoxc8 early enhancer parallels diverged axial morphologies between mammals and fishes

Genomic code for Sox2 binding uncovers its regulatory role in Six3 activation in the forebrain

The SRY-related HMG box transcription factor Sox2 plays critical roles throughout embryogenesis. Haploinsufficiency for SOX2 results in human developmental defects including anophthalmia, microphthalmia and septo-optic dysplasia, a congenital forebrain defect. To understand how Sox2 plays a role in neurogenesis, we combined genomic and in vivo transgenic approaches to characterize genomic regions occupied by Sox2 in the developing forebrain. Six3, a homeobox gene associated with holoprosencephaly, a forebrain midline defect, was identified as a Sox2 transcriptional target. This study shows that Sox2 directly regulates a previously unidentified long-range forebrain enhancer to activate Six3 expression in the rostral diencephalon. Further biochemical and genetic evidences indicated a direct regulatory link between Sox2 and Six3 during forebrain development, providing a better understanding of a common molecular mechanism underlying these forebrain defects.

Changes in Cis-regulatory Elements during Morphological Evolution

How have animals evolved new body designs (morphological evolution)? This requires explanations both for simple morphological changes, such as differences in pigmentation and hair patterns between different Drosophila populations and species, and also for more complex changes, such as differences in the forelimbs of mice and bats, and the necks of amphibians and reptiles. The genetic changes and pathways involved in these evolutionary steps require identification. Many, though not all, of these events occur by changes in cis-regulatory (enhancer) elements within developmental genes. Enhancers are modular, each affecting expression in only one or a few tissues. Therefore it is possible to add, remove or alter an enhancer without producing changes in multiple tissues, and thereby avoid widespread (pleiotropic) deleterious effects. Ideally, for a given step in morphological evolution it is necessary to identify (i) the change in phenotype, (ii) the changes in gene expression, (iii) the DNA region, enhancer or otherwise, affected, (iv) the mutation involved, (v) the nature of the transcription or other factors that bind to this site. In practice these data are incomplete for most of the published studies upon morphological evolution. Here, the investigations are categorized according to how far these analyses have proceeded.

Ancient vertebrate conserved noncoding elements have been evolving rapidly in teleost fishes

Vertebrate genomes contain thousands of conserved noncoding elements (CNEs) that often function as tissue-specific enhancers. In this study, we have identified CNEs in human, dog, chicken, Xenopus, and four teleost fishes (zebrafish, stickleback, medaka, and fugu) using elephant shark, a cartilaginous vertebrate, as the base genome and investigated the evolution of these ancient vertebrate CNEs (aCNEs) in bony vertebrate lineages. Our analysis shows that aCNEs have been evolving at different rates in different bony vertebrate lineages. Although 78-83% of CNEs have diverged beyond recognition ("lost") in different teleost fishes, only 24% and 40% have been lost in the chicken and mammalian lineages, respectively. Relative rate tests of substitution rates in CNEs revealed that the teleost fish CNEs have been evolving at a significantly higher rate than those in other bony vertebrates. In the ray-finned fish lineage, 68% of aCNEs were lost before the divergence of the four teleosts. This implicates the "fish-specific" whole-genome duplication in the accelerated evolution and the loss of a large number of both copies of duplicated CNEs in teleost fishes. The aCNEs are rich in tissue-specific enhancers and thus many of them are likely to be evolutionarily constrained cis-regulatory elements. The rapid evolution of aCNEs might have affected the expression patterns driven by them. Transgenic zebrafish assay of some human CNE enhancers that have been lost in teleosts has indicated instances of conservation or changes in trans-acting factors between mammals and fishes.

cis-Regulatory remodeling of the SCL locus during vertebrate evolution

Development progresses through a sequence of cellular identities which are determined by the activities of networks of transcription factor genes. Alterations in cis-regulatory elements of these genes play a major role in evolutionary change, but little is known about the mechanisms responsible for maintaining conserved patterns of gene expression. We have studied the evolution of cis-regulatory mechanisms controlling the SCL gene, which encodes a key transcriptional regulator of blood, vasculature, and brain development and exhibits conserved function and pattern of expression throughout vertebrate evolution. SCL cis-regulatory elements are conserved between frog and chicken but accrued alterations at an accelerated rate between 310 and 200 million years ago, with subsequent fixation of a new cis-regulatory pattern at the beginning of the mammalian radiation. As a consequence, orthologous elements shared by mammals and lower vertebrates exhibit functional differences and binding site turnover between widely separated cis-regulatory modules. However, the net effect of these alterations is constancy of overall regulatory inputs and of expression pattern. Our data demonstrate remarkable cis-regulatory remodelling across the SCL locus and indicate that stable patterns of expression can mask extensive regulatory change. These insights illuminate our understanding of vertebrate evolution.

Elephant shark (Callorhinchus milii) provides insights into the evolution of Hox gene clusters in gnathostomes

We have sequenced and analyzed Hox gene clusters from elephant shark, a holocephalian cartilaginous fish. Elephant shark possesses 4 Hox clusters with 45 Hox genes that include orthologs for a higher number of ancient gnathostome Hox genes than the 4 clusters in tetrapods and the supernumerary clusters in teleost fishes. Phylogenetic analysis of elephant shark Hox genes from 7 paralogous groups that contain all of the 4 members indicated an ((AB)(CD)) topology for the order of Hox cluster duplication, providing support for the 2R hypothesis (i.e., 2 rounds of whole-genome duplication during the early evolution of vertebrates). Comparisons of noncoding sequences of the elephant shark and human Hox clusters have identified a large number of conserved noncoding elements (CNEs), which represent putative cis-regulatory elements that may be involved in the regulation of Hox genes. Interestingly, in fugu more than 50% of these ancient CNEs have diverged beyond recognition in the duplicated (HoxA, HoxB, and HoxD) as well as the singleton (HoxC) Hox clusters. Furthermore, the b-paralogs of the duplicated fugu Hox clusters are virtually devoid of unique ancient CNEs. In contrast to fugu Hox clusters, elephant shark and human Hox clusters have lost fewer ancient CNEs. If these ancient CNEs are indeed enhancers directing tissue-specific expression of Hox genes, divergence of their sequences in vertebrate lineages might have led to altered expression patterns and presumably the functions of their associated Hox genes.

Conservation of core gene expression in vertebrate tissues

BACKGROUND

Vertebrates share the same general body plan and organs, possess related sets of genes, and rely on similar physiological mechanisms, yet show great diversity in morphology, habitat and behavior. Alteration of gene regulation is thought to be a major mechanism in phenotypic variation and evolution, but relatively little is known about the broad patterns of conservation in gene expression in non-mammalian vertebrates.

RESULTS

We measured expression of all known and predicted genes across twenty tissues in chicken, frog and pufferfish. By combining the results with human and mouse data and considering only ten common tissues, we have found evidence of conserved expression for more than a third of unique orthologous genes. We find that, on average, transcription factor gene expression is neither more nor less conserved than that of other genes. Strikingly, conservation of expression correlates poorly with the amount of conserved nonexonic sequence, even using a sequence alignment technique that accounts for non-collinearity in conserved elements. Many genes show conserved human/fish expression despite having almost no nonexonic conserved primary sequence.

CONCLUSIONS

There are clearly strong evolutionary constraints on tissue-specific gene expression. A major challenge will be to understand the precise mechanisms by which many gene expression patterns remain similar despite extensive cis-regulatory restructuring.

Comparison of diverged Hoxc8 early enhancer activities reveals modification of regulatory interactions at conserved cis-acting elements

The Hoxc8 early enhancer that controls the initiation and establishment of Hoxc8 expression in the developing mouse embryo is found in different vertebrate lineages including mammals, birds and fish. Mouse and Fugu Hoxc8 early enhancers (200 bp) have diverged in the composition of elements located towards the 3' region. However, they share cis-acting elements A-E located in the 5' region. Mutations at these elements in the context of the mouse Hoxc8 early enhancer affect reporter gene expression in the posterior neural tube, somites and lateral plate mesoderm of day 9.5 mouse embryos. Here, we demonstrate that mutations introduced at the same elements but in the context of the Fugu Hoxc8 early enhancer had different consequences on the reporter gene expression in transgenic mouse embryos. Furthermore, in contrast to the mouse enhancer the Fugu enhancer does not utilize elements D and E in achieving posterior neural tube and somite expression. These results suggest that the diverged sequences prevent regulatory interactions at conserved cis-acting elements. We propose that divergent sequences modify regulatory interactions at conserved elements by providing a "contextual change". Our finding that the enhancer elements do not act in a unitary fashion but function in the context of the surrounding sequence brings a new dimension to the study of cis-regulatory evolution.

Variable gene expression in eukaryotes: a network perspective

Changes in gene expression underlie phenotypic plasticity, variation within species, and phenotypic divergence between species. These expression differences arise from modulation of regulatory networks. To understand the source of expression differences, networks of interactions among genes and gene products that orchestrate gene expression must be considered. Here I review the basic structure of eukaryotic regulatory networks and discuss selected case studies that provide insight into how these networks are altered to create expression differences within and between species.

Identification of regulatory elements by gene family footprinting and in vivo analysis

Gene families of recently duplicated but subsequently diverged genes provide an unique opportunity for comparative analysis of regulatory elements. We have studied the human SPRR gene family of small proline rich proteins involved in barrier function of stratified squamous epithelia. These genes are all expressed in normal human keratinocytes, but respond differently to environmental insults. Comparisons of the functional promoter regions allows the rapid identification of both conserved and of novel regulatory elements that appeared after gene duplication. Competitive electrophoretic mobility shift assays can be used to confirm their presence. Here we show the power of gene family footprinting by the identification of two novel elements in the SPRR3 promoter, not present in SPRR1A and SPRR2A. One of these elements binds a protein similar to GAAP-1, a pro-apoptotic activator of IRF-1 and p53. In vivo analysis shows that this element functions as an inhibitor of SPRR3 transcription. The second novel element functions as an activator of promoter activity and is characterized by its A/T rich sequence. The latter interacting protein indeed binds through contacts in the minor groove, and strikingly, depends on the presence of calcium for DNA interaction.

Duplicated Abd-B class genes in medaka hoxAa and hoxAb clusters exhibit differential expression patterns in pectoral fin buds

Hox genes form clusters. Invertebrates and Amphioxus have only one hox cluster, but in vertebrates, they are multiple, i.e., four in the basal teleost fish Polyodon and tetrapods (HoxA, B, C, D), but seven or eight in common teleosts. We earlier completely sequenced the entire hox gene loci in medaka fish, showing a total of 46 hox genes to be encoded in seven clusters (hoxAa, Ab, Ba, Bb, Ca, Da, Db). Among them, hoxAa, hoxAb and hoxDa clusters are presumed to be important for fin-to-limb evolution because of their key role in forelimb and pectoral fin development. In the present study, we compared genome organization and nucleotide sequences of the hoxAa and hoxAb clusters to these of tetrapod HoxA clusters, and found greater similarity in hoxAa case. We then analyzed expression of Abd-B family genes in the clusters. In the trunk, those from the hoxAa cluster, i.e., hoxA9a, hoxA10a, hoxA11a and hoxA13a, were expressed in a manner keeping the colinearity rule of the hox expression as those of tetrapods, while those from the hoxAb cluster, i.e., hoxA9b, hoxA10b, hoxA11b and hoxA13b, were not. In the pectoral fins, the hoxAa cluster was expressed in split domains and did not obey the rule. By contrast, those from the hoxAb and hoxDa clusters were expressed in a manner keeping the rule, i.e., an ancestral pattern similar to those of tetrapods. It is plausible that this differential expression of the two clusters is caused by changes occurred in global control regions after cluster duplications.

The anuran Bauplan: a review of the adaptive, developmental, and genetic underpinnings of frog and tadpole morphology

Anurans (frogs, toads, and their larvae) are among the most morphologically derived of vertebrates. While tightly conserved across the order, the anuran Bauplan (body plan) diverges widely from that of other vertebrates, particularly with respect to the skeleton. Here we address the adaptive, ontogenetic, and genetic bases of three such hallmark anuran features: (1) the absence of discrete caudal vertebrae, (2) a truncated axial skeleton, and (3) elongate hind limbs. We review the functional significance of each as it relates to the anuran lifestyle, which includes locomotor adaptations to both aquatic and terrestrial environments. We then shift our focus to the proximal origins of each feature, namely, ontogeny and its molecular regulation. Drawing on relatively limited data, we detail the development of each character and then, by extrapolating from comparative vertebrate data, propose molecular bases for these processes. Cast in this light, the divergent morphology of anurans emerges as a product of evolutionary modulation of the generalised vertebrate developmental machinery. Specifically, we hypothesise that: (1) the formation of caudal vertebrae is precluded due to a failure of sclerotomes to form cartilaginous condensations, perhaps resulting from altered expression of a suite of genes, including Pax1, Pax9, Msx1, Uncx-4.1, Sonic hedgehog, and noggin; (2) anteriorised Hox gene expression in the paraxial mesoderm has led to a rostral shift of morphological boundaries of the vertebral column; and, (3) spatial and temporal shifts in Hox expression may underlie the expanded tarsal elements of the anuran hind limb. Technology is currently in place to investigate each of these scenarios in the African clawed frog Xenopus. Experimental corroboration will further our understanding of the molecular regulation of the anuran Bauplan and provide insight into the origin of vertebrate morphological diversity as well as the role of development in evolution.

Cooperation of sonic hedgehog enhancers in midline expression

In zebrafish, as in other vertebrates, the secreted signalling molecule Sonic hedgehog (Shh) is expressed in organiser regions such as the embryonic midline and the zona limitans intrathalamica (zli). To investigate the regulatory mechanisms underlying the pattern of shh expression, we carried out a systematic analysis of the intronic regulatory sequences of zebrafish shh using stable transgenesis. Deletion analysis identified the modules responsible for expression in the embryonic shield, the hypothalamus and the zli and confirmed the activities of previously identified notochord and floor plate enhancers. We detected a strong synergism between regulatory regions. The degree of synergy varied over time in the hypothalamus suggesting different mechanisms for initiation and maintenance of expression. Our data show that the pattern of shh expression in the embryonic central nervous system involves an intricate crosstalk of at least 4 different regulatory regions. When compared to the enhancer activities of the mouse Shh gene, we observed a remarkable divergence of function of structurally conserved enhancer sequences. The activating region ar-C (61% identical to SFPE2 in mouse Shh), for example, mediates floor plate expression in the mouse embryo while it directs expression in the forebrain and the notochord and only weakly in the floor plate in the zebrafish embryo. This raises doubts on the predictive power of phylogenetic footprinting and indicates a stunning divergence of function of structurally conserved regulatory modules during vertebrate evolution.

Evolutionary constraint on Otx2 neuroectoderm enhancers-deep conservation from skate to mouse and unique divergence in teleost

Otx2 is a paired type homeobox gene that plays essential roles in each step and site of head development in vertebrates. In the mouse, Otx2 expression in the anterior neuroectoderm is regulated primarily by two distinct enhancers: anterior neuroectoderm (AN) and forebrain/midbrain (FM) enhancers at 92 kb and 75 kb 5'of the Otx2 locus, respectively. The AN enhancer has activity in the entire anterior neuroectoderm at headfold and early somite stages, whereas the FM enhancer is subsequently active in the future caudal forebrain and midbrain ectoderm. In tetrapods, both AN and FM enhancers are conserved, whereas the AN region is missing in teleosts, despite overt Otx2 expression in the anterior neuroectoderm. Here, we show that zebrafish and fugu FM regions drive expression not only in the forebrain and midbrain but also in the anterior neuroectoderm at headfold stage. The analysis of coelacanth and skate genomic Otx2 orthologues suggests that the utilization of the two enhancers, AN and FM, is an ancestral condition. In contrast, the AN enhancer has been specifically lost in the teleost lineage with a compensatory establishment of AN activity within the FM enhancer. Furthermore, the AN activity in the fish FM enhancer was established by recruiting upstream factors different from those that direct the tetrapod AN enhancer, yet zebrafish FM enhancer is active in both mouse and zebrafish anterior neuroectoderm at the headfold stage.

Biomimetic evolutionary analysis: testing the adaptive value of vertebrate tail stiffness in autonomous swimming robots

For early vertebrates, a long-standing hypothesis is that vertebrae evolved as a locomotor adaptation, stiffening the body axis and enhancing swimming performance. While supported by biomechanical data, this hypothesis has not been tested using an evolutionary approach. We did so by extending biomimetic evolutionary analysis (BEA), which builds physical simulations of extinct systems, to include use of autonomous robots as proxies of early vertebrates competing in a forage navigation task. Modeled after free-swimming larvae of sea squirts (Chordata, Urochordata), three robotic tadpoles (`Tadros'), each with a propulsive tail bearing a biomimetic notochord of variable spring stiffness, k (N m-1), searched for, oriented to, and orbited in two dimensions around a light source. Within each of ten generations, we selected for increased swimming speed, U (m s-1) and decreased time to the light source, t (s),average distance from the source, R (m) and wobble maneuvering, W (rad s-2). In software simulation, we coded two quantitative trait loci (QTL) that determine k: bending modulus, E (Nm-2) and length, L (m). Both QTL were mutated during replication, independently assorted during meiosis and, as haploid gametes, entered into the gene pool in proportion to parental fitness. After random mating created three new diploid genotypes, we fabricated three new offspring tails. In the presence of both selection and chance events(mutation, genetic drift), the phenotypic means of this small population evolved. The classic hypothesis was supported in that k was positively correlated (r2=0.40) with navigational prowess, NP, the dimensionless ratio of U to the product of R, t and W. However, the plausible adaptive scenario, even in this simplified system, is more complex, since the remaining variance in NP was correlated with the residuals of R and U taken with respect to k, suggesting that changes in k alone are insufficient to explain the evolution of NP.

Conserved sequences and the evolution of gene regulatory signals

Studies of evolutionary conservation of gene regulatory signals have led to a paradox: extensive sequence similarity implies functional conservation in non-coding regions across mammalian species; however, this stands in contrast to our understanding of transcriptional regulatory sites composed of degenerate recognition sequences for transcription factors that can maintain functional equivalence despite considerable sequence divergence. The latter observation provides an explanation for the rapid evolution of new traits through the gain and loss of transcription factor binding sites that bring new genes under the control of an existing genetic regulatory network. The former observation might point to novel mechanisms of gene regulation and/or chromosome function that are currently unappreciated. Recent comparative genome analysis has highlighted extensive conserved sequences in mammalian genomes that are beginning to be functionally characterized.

Anterior epidermis-specific expression of the cuticle gene EDG84A is controlled by many cis-regulatory elements in Drosophila melanogaster

During insect metamorphosis, a pulse of ecdysteroids induces many different morphological changes depending on different parts of the body. In Drosophila, although a number of transcription factors are expressed in a stage-specific manner in response to an ecdysteroid pulse, little is known on the regulatory mechanism for space-specific gene expression during metamorphosis. The EDG84A gene encoding pupal cuticle protein is one of the targets of ecdysteroid-inducible transcription factor betaFTZ-F1 and is expressed only in anterior epidermis of the body during mid- to late prepupal period, whereas betaFTZ-F1 is expressed in almost all tissues. To address the regulatory mechanism of the tissue-specific expression of the EDG84A gene, we established transgenic fly lines which carry various upstream regions of the gene fused to the LacZ gene and examined the expression pattern of the reporter gene. Results of the transgenic fly reporter assays showed that the space-specific expression is controlled by at least four positive and two negative elements within a 263-bp region near the transcription start site, and at least three of them showed space-specific effects to the anterior body trunk. These results suggest that both high expression level and differential expression are achieved through many cis-regulatory elements.

Hoxc8 early enhancer of the Indonesian coelacanth, Latimeria menadoensis

Hoxc8 early enhancer controls the initiation and establishment phase of Hoxc8 expression in the mouse. Comparative studies indicate the presence of Hoxc8 early enhancer sequences in different vertebrate clades including mammals, birds and fish. Previous studies have shown differences between teleost and mammalian Hoxc8 early enhancers with respect to sequence and organization of protein binding elements. This raises the question of when the Hoxc8 early enhancer arose and how it has become modified in different vertebrate lineages. Here, we describe Hoxc8 early enhancer from the Indonesian coelacanth, Latimeria menadoensis. Coelacanths are the only extant lobefinned fish whose genome is tractable to genome analysis. The Latimeria Hoxc8 early enhancer sequence more closely resembles that of the mouse than that of Fugu or zebrafish. When assayed for enhancer activity by reporter gene analysis in transgenic mouse embryos, Latimeria Hoxc8 early enhancer directs expression to the posterior neural tube and mesoderm similar to that of the mouse enhancer. These observations support a close relationship between coelacanths and tetrapods and place the origin of a common Hoxc8 early enhancer sequence within the sarcopterygian lineage. The divergence of teleost (actinopterygii) Hoxc8 early enhancer may reflect a case of relaxed selection or other forms of instability induced by genome duplication events.

Comparative cis-regulatory analyses identify new elements of the mouse Hoxc8 early enhancer

The Hoxc8 early enhancer is a 200 bp region that controls the early phase of Hoxc8 expression during mouse embryonic development. This enhancer defines the domain of Hoxc8 expression in the neural tube and mesoderm of the posterior regions of the developing embryo. Five distinct cis-acting elements, A-E, were previously shown to govern early phase Hoxc8 expression. Significant divergence between mammalian and fish Hoxc8 early enhancer sequences and activities suggested additional cis-acting elements. Here we describe four additional cis-acting elements (F-I) within the 200 bp Hoxc8 early enhancer region identified by comparative regulatory analysis and transgene-mutation studies. These elements affect posterior neural tube and mesoderm expression of the reporter gene, either singly or in combination. Surprisingly, these new elements are missing from the zebrafish and Fugu Hoxc8 early enhancer sequences. Considering that fish enhancers direct robust reporter expression in transgenic mouse embryos, it is tempting to postulate that fish and mammalian Hoxc8 early enhancers utilize different sets of elements to direct Hoxc8 early expression. These observations reveal a remarkable plasticity in the Hoxc8 early enhancer, suggesting different modes of initiation and establishment of Hoxc8 expression in different species. We postulate that extensive restructuring and remodeling of Hox cis-regulatory regions occurring in different taxa lead to relatively different Hox expression patterns, which in turn may act as a driving force in generating diverse axial morphologies.

Interpreting mammalian evolution using Fugu genome comparisons

Recently, it has been shown that a significant number of evolutionarily conserved human-Fugu noncoding elements function as tissue-specific transcriptional enhancers in vivo, suggesting that distant comparisons are capable of identifying a particular class of regulatory elements. We therefore hypothesized that by juxtaposing human/Fugu and human/mouse conservation patterns we can define conservation criteria for discovering transcriptional regulatory elements specific to mammals. Genome-scale comparisons of noncoding human/Fugu evolutionary conserved elements (ECRs) and their humans/mouse counterparts revealed a particular signature common to human/mouse ECRs (>or=350 bp long, >or=77% identity) that are also conserved in fishes. This newly defined threshold identifies 90% of all human/Fugu noncoding ECRs without the assistance of human-Fugu genome alignments and provides a very efficient filter for identifying functional human/mouse ECRs.