Adaptive evolution of Hox-gene homeodomains after cluster duplications

Irreducible Complexity of Hox Gene: Path to the Canonical Function of the Hox Cluster

The evolution of major taxa is often associated with the emergence of new gene families. In all multicellular animals except sponges and comb jellies, the genomes contain Hox genes, which are crucial regulators of development. The canonical function of Hox genes involves colinear patterning of body parts in bilateral animals. This general function is implemented through complex, precisely coordinated mechanisms, not all of which are evolutionarily conserved and fully understood. We suggest that the emergence of this regulatory complexity was preceded by a stage of cooperation between more ancient morphogenetic programs or their individual elements. Footprints of these programs may be present in modern animals to execute non-canonical Hox functions. Non-canonical functions of Hox genes are involved in maintaining terminal nerve cell specificity, autophagy, oogenesis, pre-gastrulation embryogenesis, vertical signaling, and a number of general biological processes. These functions are realized by the basic properties of homeodomain protein and could have triggered the evolution of ParaHoxozoa and Nephrozoa subsequently. Some of these non-canonical Hox functions are discussed in our review.

Forkhead box O proteins: steering the course of stem cell fate

Stem cells are pivotal players in the intricate dance of embryonic development, tissue maintenance, and regeneration. Their behavior is delicately balanced between maintaining their pluripotency and differentiating as needed. Disruptions in this balance can lead to a spectrum of diseases, underscoring the importance of unraveling the complex molecular mechanisms that govern stem cell fate. Forkhead box O (FOXO) proteins, a family of transcription factors, are at the heart of this intricate regulation, influencing a myriad of cellular processes such as survival, metabolism, and DNA repair. Their multifaceted role in steering the destiny of stem cells is evident, as they wield influence over self-renewal, quiescence, and lineage-specific differentiation in both embryonic and adult stem cells. This review delves into the structural and regulatory intricacies of FOXO transcription factors, shedding light on their pivotal roles in shaping the fate of stem cells. By providing insights into the specific functions of FOXO in determining stem cell fate, this review aims to pave the way for targeted interventions that could modulate stem cell behavior and potentially revolutionize the treatment and prevention of diseases.

High-Throughput Discovery and Characterization of Human Transcriptional Effectors

Thousands of proteins localize to the nucleus; however, it remains unclear which contain transcriptional effectors. Here, we develop HT-recruit, a pooled assay where protein libraries are recruited to a reporter, and their transcriptional effects are measured by sequencing. Using this approach, we measure gene silencing and activation for thousands of domains. We find a relationship between repressor function and evolutionary age for the KRAB domains, discover that Homeodomain repressor strength is collinear with Hox genetic organization, and identify activities for several domains of unknown function. Deep mutational scanning of the CRISPRi KRAB maps the co-repressor binding surface and identifies substitutions that improve stability/silencing. By tiling 238 proteins, we find repressors as short as ten amino acids. Finally, we report new activator domains, including a divergent KRAB. These results provide a resource of 600 human proteins containing effectors and demonstrate a scalable strategy for assigning functions to protein domains.

Homeobox Gene Duplication and Divergence in Arachnids

Homeobox genes are key toolkit genes that regulate the development of metazoans and changes in their regulation and copy number have contributed to the evolution of phenotypic diversity. We recently identified a whole genome duplication (WGD) event that occurred in an ancestor of spiders and scorpions (Arachnopulmonata), and that many homeobox genes, including two Hox clusters, appear to have been retained in arachnopulmonates. To better understand the consequences of this ancient WGD and the evolution of arachnid homeobox genes, we have characterized and compared the homeobox repertoires in a range of arachnids. We found that many families and clusters of these genes are duplicated in all studied arachnopulmonates (Parasteatoda tepidariorum, Pholcus phalangioides, Centruroides sculpturatus, and Mesobuthus martensii) compared with nonarachnopulmonate arachnids (Phalangium opilio, Neobisium carcinoides, Hesperochernes sp., and Ixodes scapularis). To assess divergence in the roles of homeobox ohnologs, we analyzed the expression of P. tepidariorum homeobox genes during embryogenesis and found pervasive changes in the level and timing of their expression. Furthermore, we compared the spatial expression of a subset of P. tepidariorum ohnologs with their single copy orthologs in P. opilio embryos. We found evidence for likely subfunctionlization and neofunctionalization of these genes in the spider. Overall our results show a high level of retention of homeobox genes in spiders and scorpions post-WGD, which is likely to have made a major contribution to their developmental evolution and diversification through pervasive subfunctionlization and neofunctionalization, and paralleling the outcomes of WGD in vertebrates.

A 20 bp Duplication in Exon 2 of the Aristaless-Like Homeobox 4 Gene (ALX4) Is the Candidate Causative Mutation for Tibial Hemimelia Syndrome in Galloway Cattle

Aristaless-like homeobox 4 (ALX4) gene is an important transcription regulator in skull and limb development. In humans and mice ALX4 mutations or loss of function result in a number of skeletal and organ malformations, including polydactyly, tibial hemimelia, omphalocele, biparietal foramina, impaired mammary epithelial morphogenesis, alopecia, coronal craniosynostosis, hypertelorism, depressed nasal bridge and ridge, bifid nasal tip, hypogonadism, and body agenesis. Here we show that a complex skeletal malformation of the hind limb in Galloway cattle together with other developmental anomalies is a recessive autosomal disorder most likely caused by a duplication of 20 bp in exon 2 of the bovine ALX4 gene. A second duplication of 34 bp in exon 4 of the same gene has no known effect, although both duplications result in a frameshift and premature stop codon leading to a truncated protein. Genotyping of 1,688 Black/Red/Belted/Riggit Galloway (GA) and 289 White Galloway (WGA) cattle showed that the duplication in exon 2 has allele frequencies of 1% in GA and 6% in WGA and the duplication in exon 4 has frequencies of 23% in GA and 38% in WGA. Both duplications were not detected in 876 randomly selected German Holstein Friesian and 86 cattle of 21 other breeds. Hence, we have identified a candidate causative mutation for tibial hemimelia syndrome in Galloway cattle and selection against this mutation can be used to eliminate the mutant allele from the breed.

Adaptive evolution of Hoxc13 genes in the origin and diversification of the vertebrate integument

The problem of origination and diversification of integument derivatives in vertebrates is still a challenge. The homeobox (Hox) genes Hoxc13 control integument formation in vertebrate. Hoxc13 show strong expression in the integument development, are highly conserved across vertebrates, and show mutations that are associated with skin and appendages. To test whether the evolution of the integument is associated with positive selection or relaxation of Hoxc13, we obtained these genes in a wide range of vertebrates. In Hoxc13, we found evidence of diversifying selection after speciation during the origin of vertebrates. In addition, we found the glycine-rich regions in Hoxc13 protein in mammals, but not among non-mammalian taxa. Our results strongly implicate that Hoxc13 genes could have played an important role in the evolution of integument structure.

Quantifying protein modularity and evolvability: a comparison of different techniques

Modularity increases evolvability by reducing constraints on adaptation and by allowing preexisting parts to function in new contexts for novel uses. Protein evolution provides an excellent context to study the causes and consequences of biological modularity. In order to address such questions, however, an index for protein modularity is necessary. This paper proposes a simple index for protein modularity-"module density"-which is the number of evolutionarily independent modules that compose a protein divided by the number of amino acids in the protein. The decomposition of proteins into constituent modules can be accomplished by either of two classes of methods. The first class of methods relies on "suppositional" criteria to assign amino acids to modules, whereas the second class of methods relies on "coevolutionary" criteria for this task. One simple and practical method from the first class consists of approximating the number of modules in a protein as the number of regular secondary structure elements (i.e., helices and sheets). Methods based on coevolutionary criteria require more elaborate data, but they have the advantage of being able to specify modules without prior assumptions about why they exist. Given the increasing availability of datasets sampling protein mutational spectra (e.g., from comparative genomics, experimental evolution, and computational prediction), methods based on coevolutionary criteria will likely become more promising in the near future. The ability to meaningfully quantify protein modularity via simple indices has the potential to aid future efforts to understand protein evolutionary rate determinants, improve molecular evolution models and engineer novel proteins.

Testis-specific glyceraldehyde-3-phosphate dehydrogenase: origin and evolution

Background

Glyceraldehyde-3-phosphate dehydrogenase (GAPD) catalyses one of the glycolytic reactions and is also involved in a number of non-glycolytic processes, such as endocytosis, DNA excision repair, and induction of apoptosis. Mammals are known to possess two homologous GAPD isoenzymes: GAPD-1, a well-studied protein found in all somatic cells, and GAPD-2, which is expressed solely in testis. GAPD-2 supplies energy required for the movement of spermatozoa and is tightly bound to the sperm tail cytoskeleton by the additional N-terminal proline-rich domain absent in GAPD-1. In this study we investigate the evolutionary history of GAPD and gain some insights into specialization of GAPD-2 as a testis-specific protein.

Results

A dataset of GAPD sequences was assembled from public databases and used for phylogeny reconstruction by means of the Bayesian method. Since resolution in some clades of the obtained tree was too low, syntenic analysis was carried out to define the evolutionary history of GAPD more precisely. The performed selection tests showed that selective pressure varies across lineages and isoenzymes, as well as across different regions of the same sequences.

Conclusions

The obtained results suggest that GAPD-1 and GAPD-2 emerged after duplication during the early evolution of chordates. GAPD-2 was subsequently lost by most lineages except lizards, mammals, as well as cartilaginous and bony fishes. In reptilians and mammals, GAPD-2 specialized to a testis-specific protein and acquired the novel N-terminal proline-rich domain anchoring the protein in the sperm tail cytoskeleton. This domain is likely to have originated by exonization of a microsatellite genomic region. Recognition of the proline-rich domain by cytoskeletal proteins seems to be unspecific. Besides testis, GAPD-2 of lizards was also found in some regenerating tissues, but it lacks the proline-rich domain due to tissue-specific alternative splicing.

Molecular evolution of the EGF-CFC protein family

The epidermal growth factor-Cripto-1/FRL-1/Cryptic (EGF-CFC) proteins, characterized by the highly conserved EGF and CFC domains, are extracellular membrane associated growth factor-like glycoproteins. These proteins are essential components of the Nodal signaling pathway during early vertebrate embryogenesis. Homologs of the EGF-CFC family have also been implicated in tumorigenesis in humans. Yet, little is known about the mode of molecular evolution in this family. Here we investigate the origin, extent of conservation and evolutionary relationships of EGF-CFC proteins across the metazoa. The results suggest that the first appearance of the EGF-CFC gene occurred in the ancestor of the deuterostomes. Phylogenetic analysis supports the classification of the family into distinct subfamilies that appear to have evolved through lineage-specific duplication and divergence. Site-specific analyses of evolutionary rate shifts between the two major mammalian paralogous subfamilies, Cripto and Cryptic, reveal critical amino acid sites that may account for the observed functional divergence. Furthermore, estimates of functional divergence suggest that rapid change of evolutionary rates at sites located mainly in the CFC domain may contribute towards distinct functional properties of the two paralogs.

The evolution of vertebrate tetraspanins: gene loss, retention, and massive positive selection after whole genome duplications

BACKGROUND

The vertebrate tetraspanin family has many features which make it suitable for preserving the imprint of ancient sequence evolution and amenable for phylogenomic analysis. So we believe that an in-depth analysis of the tetraspanin evolution not only provides more complete understanding of tetraspanin biology, but offers new insights into the influence of the two rounds of whole genome duplication (2R-WGD) at the origin of vertebrates.

RESULTS

A detailed phylogeny of vertebrate tetraspanins was constructed by using multiple lines of information, including sequence-based phylogenetics, key structural features, intron configuration and genomic synteny. In particular, a total of 38 modern tetraspanin ortholog lineages in bony vertebrates have been identified and subsequently classified into 17 ancestral lineages existing before 2R-WGD. Based on this phylogeny, we found that the ohnolog retention rate of tetraspanins after 2R-WGD was three times as the average (a rate similar to those of transcription factors and protein kinases). This high rate didn't increase the tetrapanin family size, but changed the family composition, possibly by displacing vertebrate-specific gene lineages with the lineages conserved across deuterostomes. We also found that the period from 2R-WGD to recent time is controlled by gene losses. Meanwhile, positive selection has been detected on 80% of the branches right after 2R-WGDs, which declines significantly on both magnitude and extensity on the following speciation branches. Notably, the loss of mammalian RDS2 is accompanied by strong positive selection on mammalian ROM1, possibly due to gene loss-induced compensatory evolution.

CONCLUSIONS

First, different from transcription factors and kinases, high duplicate retention rate after 2R-WGD didn't increase the tetraspanin family size but just reshaped the family composition. Second, the evolution of tetraspanins right after 2R-WGD had been impacted by a massive wave of gene loss and positive selection on coding sequences. Third, the lingering effect of 2R-WGD on tetraspanin gene loss and positive selection might last for 300-400 million years.

Are the deuterostome posterior Hox genes a fast-evolving class?

There has been a great deal of interest in analysing the molecular evolution of the Hox cluster using both bioinformatic and experimental approaches. The posterior Hox genes have been of particular interest to both groups of biologists for a number of reasons: they appear to be associated with the evolution of a number of morphological novelties; the protostomes appear to be have lost a highly-conserved and functionally important amino acid motif (the hexapeptide motif) from their posterior Hox genes; and deuterostome posterior Hox genes seem to be evolving more quickly than all other Hox genes. In this chapter I will discuss the last of these points. The idea that Deuterostome posterior Hox genes were evolving more quickly than other Hox genes was first suggested by David Ferrier and colleagues. In this chapter, I start by introducing the posterior Hox genes--their distribution among the animal phyla and the likely sequence of duplications that led to this distribution. I then introduce the idea of 'deuterostome posterior flexibility' and examine this hypothesis in light of more recent phylogenetic and genomic work on the Hox cluster. Finally, I discuss some new approaches that could be used to test directly for differential rates of evolution among Hox genes and to assess what might underlie these differences.

The origin of conserved protein domains and amino acid repeats via adaptive competition for control over amino acid residues

Some proteins, such as homeodomain transcription factors, contain highly conserved regions of sequence. It has recently been suggested that multiple functional domains overlap in the homeodomain, together explaining this high conservation. However, the question remains why so many functional domains cluster together in one relatively small and constrained region of the protein. Here we have modeled an evolutionary mechanism that can produce this kind of clustering: conserved functional domains are displaced from the parts of the molecule that are undergoing adaptive evolution because novel functions generally out-compete conserved functions for control over the identity of amino acid residues. We call this model COAA, for Competition Over Amino Acids. We also studied the evolution of amino acid repeats (a.k.a. homopeptides), which are especially prevalent in transcription factors. Repeats that are encoded by non-homogenous mixtures of synonymous codons cannot be explained by replication slippage alone. Our model provides two explanations for their origin, maintenance, and over-representation in highly conserved proteins. We demonstrate that either competition between multiple functional domains for space within a sequence, or reuse of a sequence for many functions over time, can cause the evolution of amino acid repeats. Both of these processes are characteristic of multifunctional proteins such as homeodomain transcription factors. We conclude that the COAA model can explain two widely recognized features of transcription factor proteins: conserved domains and a tendency to accumulate homopeptides.

FoxO gene family evolution in vertebrates

Background

Forkhead box, class O (FoxO) belongs to the large family of forkhead transcription factors that are characterized by a conserved forkhead box DNA-binding domain. To date, the FoxO group has four mammalian members: FoxO1, FoxO3a, FoxO4 and FoxO6, which are orthologs of DAF16, an insulin-responsive transcription factor involved in regulating longevity of worms and flies. The degree of homology between these four members is high, especially in the forkhead domain, which contains the DNA-binding interface. Yet, mouse FoxO knockouts have revealed that each FoxO gene has its unique role in the physiological process. Whether the functional divergences are primarily due to adaptive selection pressure or relaxed selective constraint remains an open question. As such, this study aims to address the evolutionary mode of FoxO, which may lead to the functional divergence.

Results

Sequence similarity searches have performed in genome and scaffold data to identify homologues of FoxO in vertebrates. Phylogenetic analysis was used to characterize the family evolutionary history by identifying two duplications early in vertebrate evolution. To determine the mode of evolution in vertebrates, we performed a rigorous statistical analysis with FoxO gene sequences, including relative rate ratio tests, branch-specific d_N/d_S ratio tests, site-specific d_N/d_S ratio tests, branch-site d_N/d_S ratio tests and clade level amino acid conservation/variation patterns analysis. Our results suggest that FoxO is constrained by strong purifying selection except four sites in FoxO6, which have undergone positive Darwinian selection. The functional divergence in this family is best explained by either relaxed purifying selection or positive selection.

Conclusion

We present a phylogeny describing the evolutionary history of the FoxO gene family and show that the genes have evolved through duplications followed by purifying selection except for four sites in FoxO6 fixed by positive selection lie mostly within the non-conserved optimal PKB motif in the C-terminal part. Relaxed selection may play important roles in the process of functional differentiation evolved through gene duplications as well.

Implications of duplicated cis-regulatory elements in the evolution of metazoans: the DDI model or how simplicity begets novelty

The discovery that most regulatory genes were conserved among animals from distant phyla challenged the ideas that gene duplication and divergence of homologous coding sequences were the basis for major morphological changes in metazoan evolution. In recent years, however, the interest for the roles, conservation and changes of non-coding sequences grew-up in parallel with genome sequencing projects. Presently, many independent studies are highlighting the importance that subtle changes in cis-regulatory regions had in the evolution of morphology trough the Animal Kingdom. Here we will show and discuss some of these studies, and underscore the future of cis-Evo-Devo research. Nevertheless, we would also explore how gene duplication, which includes duplication of regulatory regions, may have been critical for spatial or temporal co-option of new regulatory networks, causing the deployment of new transcriptome scenarios, and how these induced morphological changes were critical for the evolution of new forms. Forty years after Susumu Ohno famous sentence 'natural selection merely modifies, while redundancy creates', we suggest the alternative: 'natural selection modifies, while redundancy of cis-regulatory elements innovates', and propose the Duplication-Degeneration-Innovation model to explain the increased evolvability of duplicated cis-regulatory regions. Paradoxically, making regulation simpler by subfunctionalization paved the path for future complexity or, in other words, 'to make it simple to make it complex'.

Positive selection and ancient duplications in the evolution of class B floral homeotic genes of orchids and grasses

BACKGROUND

Positive selection is recognized as the prevalence of nonsynonymous over synonymous substitutions in a gene. Models of the functional evolution of duplicated genes consider neofunctionalization as key to the retention of paralogues. For instance, duplicate transcription factors are specifically retained in plant and animal genomes and both positive selection and transcriptional divergence appear to have played a role in their diversification. However, the relative impact of these two factors has not been systematically evaluated. Class B MADS-box genes, comprising DEF-like and GLO-like genes, encode developmental transcription factors essential for establishment of perianth and male organ identity in the flowers of angiosperms. Here, we contrast the role of positive selection and the known divergence in expression patterns of genes encoding class B-like MADS-box transcription factors from monocots, with emphasis on the family Orchidaceae and the order Poales. Although in the monocots these two groups are highly diverse and have a strongly canalized floral morphology, there is no information on the role of positive selection in the evolution of their distinctive flower morphologies. Published research shows that in Poales, class B-like genes are expressed in stamens and in lodicules, the perianth organs whose identity might also be specified by class B-like genes, like the identity of the inner tepals of their lily-like relatives. In orchids, however, the number and pattern of expression of class B-like genes have greatly diverged.

RESULTS

The DEF-like genes from Orchidaceae form four well-supported, ancient clades of orthologues. In contrast, orchid GLO-like genes form a single clade of ancient orthologues and recent paralogues. DEF-like genes from orchid clade 2 (OMADS3-like genes) are under less stringent purifying selection than the other orchid DEF-like and GLO-like genes. In comparison with orchids, purifying selection was less stringent in DEF-like and GLO-like genes from Poales. Most importantly, positive selection took place before the major organ reduction and losses in the floral axis that eventually yielded the zygomorphic grass floret.

CONCLUSION

In DEF-like genes of Poales, positive selection on the region mediating interactions with other proteins or DNA could have triggered the evolution of the regulatory mechanisms behind the development of grass-specific reproductive structures. Orchidaceae show a different trend, where gene duplication and transcriptional divergence appear to have played a major role in the canalization and modularization of perianth development.

A Study of the Distribution of Phylogenetically Conserved Blocks within Clusters of Mammalian Homeobox Genes

Genome sequencing efforts of the last decade have produced a large amount of data, which has enabled whole-genome comparative analyses in order to locate potentially functional elements and study the overall patterns of phylogenetic conservation. In this paper we present a statistically based method for the characterization of these patterns in mammalian DNA sequences. We have applied this approach to the study of exceptionally well conserved homeobox gene clusters (Hox), based on the alignment of six species, and we have constructed a map of Hox cataloguing the conserved fragments, along with their locations in relation to the genes and other landmarks, sometimes showing unexpected layouts.

Multiple chromosomal rearrangements structured the ancestral vertebrate Hox-bearing protochromosomes

While the proposal that large-scale genome expansions occurred early in vertebrate evolution is widely accepted, the exact mechanisms of the expansion--such as a single or multiple rounds of whole genome duplication, bloc chromosome duplications, large-scale individual gene duplications, or some combination of these--is unclear. Gene families with a single invertebrate member but four vertebrate members, such as the Hox clusters, provided early support for Ohno's hypothesis that two rounds of genome duplication (the 2R-model) occurred in the stem lineage of extant vertebrates. However, despite extensive study, the duplication history of the Hox clusters has remained unclear, calling into question its usefulness in resolving the role of large-scale gene or genome duplications in early vertebrates. Here, we present a phylogenetic analysis of the vertebrate Hox clusters and several linked genes (the Hox "paralogon") and show that different phylogenies are obtained for Dlx and Col genes than for Hox and ErbB genes. We show that these results are robust to errors in phylogenetic inference and suggest that these competing phylogenies can be resolved if two chromosomal crossover events occurred in the ancestral vertebrate. These results resolve conflicting data on the order of Hox gene duplications and the role of genome duplication in vertebrate evolution and suggest that a period of genome reorganization occurred after genome duplications in early vertebrates.

Adaptive evolution of 5'HoxD genes in the origin and diversification of the cetacean flipper

The homeobox (Hox) genes Hoxd12 and Hoxd13 control digit patterning and limb formation in tetrapods. Both show strong expression in the limb bud during embryonic development, are highly conserved across vertebrates, and show mutations that are associated with carpal, metacarpal, and phalangeal deformities. The most dramatic evolutionary reorganization of the mammalian limb has occurred in cetaceans (whales, dolphins, and porpoises), in which the hind limbs have been lost and the forelimbs have evolved into paddle-shaped flippers. We reconstructed the phylogeny of digit patterning in mammals and inferred that digit number has changed twice in the evolution of the cetacean forelimb. First, the divergence of the early cetaceans from their even-toed relatives coincided with the reacquisition of the pentadactyl forelimb, whereas the ancestors of tetradactyl baleen whales (Mysticeti) later lost a digit again. To test whether the evolution of the cetacean forelimb is associated with positive selection or relaxation of Hoxd12 and Hoxd13, we sequenced these genes in a wide range of mammals. In Hoxd12, we found evidence of Darwinian selection associated with both episodes of cetacean forelimb reorganization. In Hoxd13, we found a novel expansion of a polyalanine tract in cetaceans compared with other mammals (17/18 residues vs. 14/15 residues, respectively), lengthening of which has previously been shown to be linked to synpolydactyly in humans and mice. Both genes also show much greater sequence variation among cetaceans than across other mammalian lineages. Our results strongly implicate 5'HoxD genes in the modulation of digit number, web forming, and the high morphological diversity of the cetacean manus.

Unexpected dynamic gene family evolution in algal actins

Actin is a conserved cytoskeletal protein that is well studied in model organisms although much less is known about actin molecular evolution in taxonomically diverse algae. Here, we analyzed 107 novel partial algal actin sequences and report some unexpected results. First, monophyletic actin gene families in multiple, phylogenetically distantly related algal taxa contain two distinct clades of sequences. One of these clades contains highly conserved sequences, whereas the second has multiple members with a significantly elevated substitution rate. This rate difference is associated with an excess of synonymous substitutions, strongly suggesting that both isoforms are active. These results paint a novel picture of actin gene evolution in algae showing it to be a remarkably dynamic system with duplication, homogenization, and potential functional diversification occurring independently in distantly related lineages.

Resurrecting the role of transcription factor change in developmental evolution

A long-standing question in evolutionary and developmental biology concerns the relative contribution of cis-regulatory and protein changes to developmental evolution. Central to this argument is which mutations generate evolutionarily relevant phenotypic variation? A review of the growing body of evolutionary and developmental literature supports the notion that many developmentally relevant differences occur in the cis-regulatory regions of protein-coding genes, generally to the exclusion of changes in the protein-coding region of genes. However, accumulating experimental evidence demonstrates that many of the arguments against a role for proteins in the evolution of gene regulation, and the developmental evolution in general, are no longer supported and there is an increasing number of cases in which transcription factor protein changes have been demonstrated in evolution. Here, we review the evidence that cis-regulatory evolution is an important driver of phenotypic evolution and provide examples of protein-mediated developmental evolution. Finally, we present an argument that the evolution of proteins may play a more substantial, but thus far underestimated, role in developmental evolution.

Unexpected novel relational links uncovered by extensive developmental profiling of nuclear receptor expression

Nuclear receptors (NRs) are transcription factors that are implicated in several biological processes such as embryonic development, homeostasis, and metabolic diseases. To study the role of NRs in development, it is critically important to know when and where individual genes are expressed. Although systematic expression studies using reverse transcriptase PCR and/or DNA microarrays have been performed in classical model systems such as Drosophila and mouse, no systematic atlas describing NR involvement during embryonic development on a global scale has been assembled. Adopting a systems biology approach, we conducted a systematic analysis of the dynamic spatiotemporal expression of all NR genes as well as their main transcriptional coregulators during zebrafish development (101 genes) using whole-mount in situ hybridization. This extensive dataset establishes overlapping expression patterns among NRs and coregulators, indicating hierarchical transcriptional networks. This complete developmental profiling provides an unprecedented examination of expression of NRs during embryogenesis, uncovering their potential function during central nervous system and retina formation. Moreover, our study reveals that tissue specificity of hormone action is conferred more by the receptors than by their coregulators. Finally, further evolutionary analyses of this global resource led us to propose that neofunctionalization of duplicated genes occurs at the levels of both protein sequence and RNA expression patterns. Altogether, this expression database of NRs provides novel routes for leading investigation into the biological function of each individual NR as well as for the study of their combinatorial regulatory circuitry within the superfamily.

NRs are key molecules controlling development, metabolism, and reproduction in metazoans. Since NRs are implicated in many human diseases such as cancer, metabolic syndrome, and hormone resistance, they are important pharmaceutical targets and are under intense scrutiny to better understand their biological functions. In the present study, we determined the expression patterns of all NR genes as well as their main transcriptional coregulators during zebrafish development. We used zebrafish because the transparency of its embryo allows us to perform whole-mount in situ hybridization from early development to late organogenesis. This complete developmental profiling offers an unprecedented view of NR expression during embryogenesis, uncovering their potential function during central nervous system and retina formation. We observed that in contrast to NR genes, only a few coregulators exhibit a restricted expression pattern, suggesting that tissue specificity of hormone action is conferred more by the receptors than by their coregulators. Lastly, by evolutionary analysis of expression pattern divergence of duplicated genes, we observed that neofunctionalization occurs at the levels of both protein sequence and mRNA expression patterns. Taken together, our data provide the starting point for functional analysis of an entire gene family during development and call for the study of the intersection between metabolism and development.

SNPs and Hox gene mapping in Ciona intestinalis

BACKGROUND

The tunicate Ciona intestinalis (Enterogona, Ascidiacea), a major model system for evolutionary and developmental genetics of chordates, harbours two cryptic species. To assess the degree of intra- and inter-specific genetic variability, we report the identification and analysis of C. intestinalis SNP (Single Nucleotide Polymorphism) markers. A SNP subset was used to determine the genetic distance between Hox-5 and -10 genes.

RESULTS

DNA fragments were amplified from 12 regions of C. intestinalis sp. A. In total, 128 SNPs and 32 one bp indels have been identified within 8 Kb DNA. SNPs in coding regions cause 4 synonymous and 12 non-synonymous substitutions. The highest SNP frequency was detected in the Hox5 and Hox10 intragenic regions. In C. intestinalis, these two genes have lost their archetypal topology within the cluster, such that Hox10 is located between Hox4 and Hox5. A subset of the above primers was used to perform successful amplification in C. intestinalis sp. B. In this cryptic species, 62 SNPs were identified within 3614 bp: 41 in non-coding and 21 in coding regions. The genetic distance of the Hox-5 and -10 loci, computed combining a classical backcross approach with the application of SNP markers, was found to be 8.4 cM (Haldane's function). Based on the physical distance, 1 cM corresponds to 39.5 Kb. Linkage disequilibrium between the aforementioned loci was calculated in the backcross generation.

CONCLUSION

SNPs here described allow analysis and comparisons within and between C. intestinalis cryptic species. We provide the first reliable computation of genetic distance in this important model chordate. This latter result represents an important platform for future studies on Hox genes showing deviations from the archetypal topology.

The developmental genetics of homology

Homology is an essential idea of biology, referring to the historical continuity of characters, but it is also conceptually highly elusive. The main difficulty is the apparently loose relationship between morphological characters and their genetic basis. Here I propose that it is the historical continuity of gene regulatory networks rather than the expression of individual homologous genes that underlies the homology of morphological characters. These networks, here referred to as 'character identity networks', enable the execution of a character-specific developmental programme.

Duplicated Hox genes in the spider Cupiennius salei

Background

Hox genes are expressed in specific domains along the anterior posterior body axis and define the regional identity. In most animals these genes are organized in a single cluster in the genome and the order of the genes in the cluster is correlated with the anterior to posterior expression of the genes in the embryo. The conserved order of the various Hox gene orthologs in the cluster among most bilaterians implies that such a Hox cluster was present in their last common ancestor. Vertebrates are the only metazoans so far that have been shown to contain duplicated Hox clusters, while all other bilaterians seem to possess only a single cluster.

Results

We here show that at least three Hox genes of the spider Cupiennius salei are present as two copies in this spider. In addition to the previously described duplicated Ultrabithorax gene, we here present sequence and expression data of a second Deformed gene, and of two Sex comb reduced genes. In addition, we describe the sequence and expression of the Cupiennius proboscipedia gene. The spider Cupiennius salei is the first chelicerate for which orthologs of all ten classes of arthropod Hox genes have been described. The posterior expression boundary of all anterior Hox genes is at the tagma border of the prosoma and opisthosoma, while the posterior boundary of the posterior Hox genes is at the posterior end of the embryo.

Conclusion

The presence of at least three duplicated Hox genes points to a major duplication event in the lineage to this spider, perhaps even of the complete Hox cluster as has taken place in the lineage to the vertebrates. The combined data of all Cupiennius Hox genes reveal the existence of two distinct posterior expression boundaries that correspond to morphological tagmata boundaries.