Evolution of a transcriptional repression domain in an insect Hox protein

Gene expression suggests conserved aspects of Hox gene regulation in arthropods and provides additional support for monophyletic Myriapoda

Antisense transcripts of Ultrabithorax (aUbx) in the millipede Glomeris and the centipede Lithobius are expressed in patterns complementary to that of the Ubx sense transcripts. A similar complementary expression pattern has been described for non-coding RNAs (ncRNAs) of the bithoraxoid (bxd) locus in Drosophila, in which the transcription of bxd ncRNAs represses Ubx via transcriptional interference. We discuss our findings in the context of possibly conserved mechanisms of Ubx regulation in myriapods and the fly.

Bicistronic transcription of Ubx and Antennapedia (Antp) has been reported previously for a myriapod and a number of crustaceans. In this paper, we show that Ubx/Antp bicistronic transcripts also occur in Glomeris and an onychophoran, suggesting further conserved mechanisms of Hox gene regulation in arthropods.

Myriapod monophyly is supported by the expression of aUbx in all investigated myriapods, whereas in other arthropod classes, including the Onychophora, aUbx is not expressed. Of the two splice variants of Ubx/Antp only one could be isolated from myriapods, representing a possible further synapomorphy of the Myriapoda.

The expression of wingless and Engrailed in developing embryos of the mayfly Ephoron leukon (Ephemeroptera: Polymitarcyidae)

The expression of the segment polarity genes wingless (wg) and engrailed (en) is highly conserved across arthropods, and these genes play a crucial role in patterning of the segmental body plan. Investigations of the expression and function of wg and en have focused primarily upon holometabolous insects, with the notable exception of recent detailed work in Oncopeltus (Hemiptera), Schistocerca, and Gryllus (Orthoptera). An increase in the phylogenetic breadth of our understanding of molecular patterning is crucial to ascertain the extent of conservation and divergence in molecular patterning mechanisms during insect embryogenesis. We examined the expression of wg mRNA transcripts and localization of En protein during embryogenesis in the mayfly Ephoron leukon (Ephemeroptera: Polymitarcyidae). These data represent one of the first embryonic gene expression pattern data for a mayfly, a lineage that may be the sister group to all other winged insects. Many aspects of wg and En expression are highly conserved, notably their expression in juxtaposed stripes in each parasegment, as well as expression domains in the procephalon, mouthparts, thoracic limbs, and nervous system. Future work in mayflies can be used to determine if conservation extends to other components of the segmentation hierarchy.

Comparing anterior and posterior Hox complex formation reveals guidelines for predicting cis-regulatory elements

Hox transcription factors specify numerous cell fates along the anterior-posterior axis by regulating the expression of downstream target genes. While expression analysis has uncovered large numbers of de-regulated genes in cells with altered Hox activity, determining which are direct versus indirect targets has remained a significant challenge. Here, we characterize the DNA binding activity of Hox transcription factor complexes on eight experimentally verified cis-regulatory elements. Hox factors regulate the activity of each element by forming protein complexes with two cofactor proteins, Extradenticle (Exd) and Homothorax (Hth). Using comparative DNA binding assays, we found that a number of flexible arrangements of Hox, Exd, and Hth binding sites mediate cooperative transcription factor complexes. Moreover, analysis of a Distal-less regulatory element (DMXR) that is repressed by abdominal Hox factors revealed that suboptimal binding sites can be combined to form high affinity transcription complexes. Lastly, we determined that the anterior Hox factors are more dependent upon Exd and Hth for complex formation than posterior Hox factors. Based upon these findings, we suggest a general set of guidelines to serve as a basis for designing bioinformatics algorithms aimed at identifying Hox regulatory elements using the wealth of recently sequenced genomes.

Natural selection drives the accumulation of amino acid tandem repeats in human proteins

Amino acid tandem repeats are found in a large number of eukaryotic proteins. They are often encoded by trinucleotide repeats and exhibit high intra- and interspecies size variability due to the high mutation rate associated with replication slippage. The extent to which natural selection is important in shaping amino acid repeat evolution is a matter of debate. On one hand, their high frequency may simply reflect their high probability of expansion by slippage, and they could essentially evolve in a neutral manner. On the other hand, there is experimental evidence that changes in repeat size can influence protein-protein interactions, transcriptional activity, or protein subcellular localization, indicating that repeats could be functionally relevant and thus shaped by selection. To gauge the relative contribution of neutral and selective forces in amino acid repeat evolution, we have performed a comparative analysis of amino acid repeat conservation in a large set of orthologous proteins from 12 vertebrate species. As a neutral model of repeat evolution we have used sequences with the same DNA triplet composition as the coding sequences--and thus expected to be subject to the same mutational forces--but located in syntenic noncoding genomic regions. The results strongly indicate that selection has played a more important role than previously suspected in amino acid tandem repeat evolution, by increasing the repeat retention rate and by modulating repeat size. The data obtained in this study have allowed us to identify a set of 92 repeats that are postulated to play important functional roles due to their strong selective signature, including five cases with direct experimental evidence.

Alternative splicing modulates Ubx protein function in Drosophila melanogaster

The Drosophila Hox gene Ultrabithorax (Ubx) produces a family of protein isoforms through alternative splicing. Isoforms differ from one another by the presence of optional segments-encoded by individual exons-that modify the distance between the homeodomain and a cofactor-interaction module termed the "YPWM" motif. To investigate the functional implications of Ubx alternative splicing, here we analyze the in vivo effects of the individual Ubx isoforms on the activation of a natural Ubx molecular target, the decapentaplegic (dpp) gene, within the embryonic mesoderm. These experiments show that the Ubx isoforms differ in their abilities to activate dpp in mesodermal tissues during embryogenesis. Furthermore, using a Ubx mutant that reduces the full Ubx protein repertoire to just one single isoform, we obtain specific anomalies affecting the patterning of anterior abdominal muscles, demonstrating that Ubx isoforms are not functionally interchangeable during embryonic mesoderm development. Finally, a series of experiments in vitro reveals that Ubx isoforms also vary in their capacity to bind DNA in presence of the cofactor Extradenticle (Exd). Altogether, our results indicate that the structural changes produced by alternative splicing have functional implications for Ubx protein function in vivo and in vitro. Since other Hox genes also produce splicing isoforms affecting similar protein domains, we suggest that alternative splicing may represent an underestimated regulatory system modulating Hox gene specificity during fly development.

Regulation of Hox activity: insights from protein motifs

Deciphering the molecular bases of animal body plan construction is a central question in developmental and evolutionary biology. Genome analyses of a number of metazoans indicate that widely conserved regulatory molecules underlie the amazing diversity of animal body plans, suggesting that these molecules are reiteratively used for multiple purposes. Hox proteins constitute a good example of such molecules and provide the framework to address the mechanisms underlying transcriptional specificity and diversity in development and evolution. Here we examine the current knowledge of the molecular bases of Hox-mediated transcriptional control, focusing on how this control is encoded within protein sequences and structures. The survey suggests that the homeodomain is part of an extended multifunctional unit coordinating DNA binding and activity regulation and highlights the need for further advances in our understanding of Hox protein activity.

Functional analysis of Ultrabithorax in the silkworm, Bombyx mori, using RNAi

The formation of abdominal appendages in insects is suppressed by the Hox genes Ultrabithorax (Ubx) and abdominal-A (abd-A), but mechanisms of the suppression can differ among species. As the function of Ubx and abd-A has been described in only a few species, more data from various insects are necessary to elucidate the evolutionary transition of regulation on abdominal appendages. We examined the function of Ubx in the silkworm Bombyx mori (Bm-Ubx) by embryonic RNA interference (RNAi). This is the first case in which functional analysis for Ubx is performed in lepidopteran insects. Larvae treated with Bm-Ubx dsRNA displayed an additional pair of thoracic leg-like protuberances in A1, whereas the other abdominal segments had no transformation. Our results suggest that Bm-Ubx is a suppressor of leg development in A1.

Probing the evolution of appendage specialization by Hox gene misexpression in an emerging model crustacean

Changes in the expression of Hox genes have been widely linked to the evolution of animal body plans, but functional demonstrations of this relationship have been impeded by the lack of suitable model organisms. A classic case study involves the repeated evolution of specialized feeding appendages, called maxillipeds, from anterior thoracic legs, in many crustacean lineages. These leg-to-maxilliped transformations correlate with the loss of Ultrabithorax (Ubx) expression from corresponding segments, which is proposed to be the underlying genetic cause. To functionally test this hypothesis, we establish tools for conditional misexpression and use these to misexpress Ubx in the crustacean Parhyale hawaiensis. Ectopic Ubx leads to homeotic transformations of anterior appendages toward more posterior thoracic fates, including maxilliped-to-leg transformations, confirming the capacity of Ubx to control thoracic (leg) versus gnathal (feeding) segmental identities. We find that maxillipeds not only are specified in the absence of Ubx, but also can develop in the presence of low/transient Ubx expression. Our findings suggest a path for the gradual evolutionary transition from thoracic legs to maxillipeds, in which stepwise changes in Hox gene expression have brought about this striking morphological and functional transformation.

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

Rapid sequence evolution of transcription factors controlling neuron differentiation in Caenorhabditis

Whether phenotypic evolution proceeds predominantly through changes in regulatory sequences is a controversial issue in evolutionary genetics. Ample evidence indicates that the evolution of gene regulatory networks via changes in cis-regulatory sequences is an important determinant of phenotypic diversity. However, recent experimental work suggests that the role of transcription factor (TF) divergence in developmental evolution may be underestimated. In order to help understand what levels of constraints are acting on the coding sequence of developmental regulatory genes, evolutionary rates were investigated among 48 TFs required for neuronal development in Caenorhabditis elegans. Allelic variation was then sampled for 28 of these genes within a population of the related species Caenorhabditis remanei. Neuronal TFs are more divergent, both within and between species, than structural genes. TFs affecting different neuronal classes are under different levels of selective constraints. The regulatory genes controlling the differentiation of chemosensory neurons evolve particularly fast and exhibit higher levels of within- and between-species nucleotide variation than TFs required for the development of several neuronal classes and TFs required for motorneuron differentiation. The TFs affecting chemosensory neuron development are also more divergent than chemosensory genes expressed in the neurons they differentiate. These results illustrate that TFs are not as highly constrained as commonly thought and suggest that the role of divergence in developmental regulatory genes during the evolution of gene regulatory networks requires further attention.

Classification of sequence signatures: a guide to Hox protein function

Hox proteins are part of the conserved superfamily of homeodomain-containing transcription factors and play fundamental roles in shaping animal body plans in development and evolution. However, molecular mechanisms underlying their diverse and specific biological functions remain largely enigmatic. Here, we have analyzed Hox sequences from the main evolutionary branches of the Bilateria group. We have found that four classes of Hox protein signatures exist, which together provide sufficient support to explain how different Hox proteins differ in their control and function. The homeodomain and its surrounding sequences accumulate nearly all signatures, constituting an extended module where most of the information distinguishing Hox proteins is concentrated. Only a small fraction of these signatures has been investigated at the functional level, but these show that approaches relying on Hox protein alterations still have a large potential for deciphering molecular mechanisms of Hox differential control.

Internal regulatory interactions determine DNA binding specificity by a Hox transcription factor

In developing bilaterans, the Hox transcription factor family regulates batteries of downstream genes to diversify serially repeated units. Given Hox homeodomains bind a wider array of DNA binding sites in vitro than are regulated by the full-length protein in vivo, regions outside the homeodomain must aid DNA site selection. Indeed, we find affinity for disparate DNA sequences varies less than 3-fold for the homeodomain isolated from the Drosophila Hox protein Ultrabithorax Ia (UbxHD), whereas for the full-length protein (UbxIa) affinity differs by more than 10-fold. The rank order of preferred DNA sequences also differs, further demonstrating distinct DNA binding preferences. The increased specificity of UbxIa can be partially attributed to the I1 region, which lies adjacent to the homeodomain and directly impacts binding energetics. Each of three segments within I1-the Extradenticle-binding YPWM motif, the six amino acids immediately N-terminal to this motif, and the eight amino acids abutting the YPWM C-terminus-uniquely contribute to DNA specificity. Combination of these regions synergistically modifies DNA binding to further enhance specificity. Intriguingly, the presence of the YPWM motif in UbxIa inhibits DNA binding only to Ubx-Extradenticle heterodimer binding sites, potentially functioning in vivo to prevent Ubx monomers from binding and misregulating heterodimer target genes. However, removal of the surrounding region allows the YPWM motif to also inhibit binding to Hox-only recognition sequences. Despite a modular domain design for Hox proteins, these results suggest that multiple Hox protein regions form a network of regulatory interactions that coordinate context- and gene-specific responses. Since most nonhomeodomain regions are not conserved between Hox family members, these regulatory interactions have the potential to diversify binding by the highly homologous Hox homeodomains.

Domain shuffling and the evolution of vertebrates

The evolution of vertebrates has included a number of important events: the development of cartilage, the immune system, and complicated craniofacial structures. Here, we examine domain shuffling as one of the mechanisms that contributes novel genetic material required for vertebrate evolution. We mapped domain-shuffling events during the evolution of deuterostomes with a focus on how domain shuffling contributed to the evolution of vertebrate- and chordate-specific characteristics. We identified approximately 1000 new domain pairs in the vertebrate lineage, including approximately 100 that were shared by all seven of the vertebrate species examined. Some of these pairs occur in the protein components of vertebrate-specific structures, such as cartilage and the inner ear, suggesting that domain shuffling made a marked contribution to the evolution of vertebrate-specific characteristics. The evolutionary history of the domain pairs is traceable; for example, the Xlink domain of aggrecan, one of the major components of cartilage, was originally utilized as a functional domain of a surface molecule of blood cells in protochordate ancestors, and it was recruited by the protein of the matrix component of cartilage in the vertebrate ancestor. We also identified genes that were created as a result of domain shuffling in ancestral chordates. Some of these are involved in the functions of chordate structures, such as the endostyle, Reissner's fiber of the neural tube, and the notochord. Our analyses shed new light on the role of domain shuffling, especially in the evolution of vertebrates and chordates.

Wings, horns, and butterfly eyespots: how do complex traits evolve?

Do novel complex traits evolve when pre-existent complex developmental networks are recruited into new places in the body? Here we look closely at the genomic fingerprints that are produced as a result of gene network co-option.

Axial patterning in snakes and caecilians: evidence for an alternative interpretation of the Hox code

It is generally assumed that the characteristic deregionalized body plan of species with a snake-like morphology evolved through a corresponding homogenization of Hox gene expression domains along the primary axis. Here, we examine the expression of Hox genes in snake embryos and show that a collinear pattern of Hox expression is retained within the paraxial mesoderm of the trunk. Genes expressed at the anterior and most posterior, regionalized, parts of the skeleton correspond to the expected anatomical boundaries. Unexpectedly however, also the dorsal (thoracic), homogenous rib-bearing region of trunk, is regionalized by unconventional gradual anterior limits of Hox expression that are not obviously reflected in the skeletal anatomy. In the lateral plate mesoderm we also detect regionalized Hox expression yet the forelimb marker Tbx5 is not restricted to a rudimentary forelimb domain but is expressed throughout the entire flank region. Analysis of several Hox genes in a caecilian amphibian, which convergently evolved a deregionalized body plan, reveals a similar global collinear pattern of Hox expression. The differential expression of posterior, vertebra-modifying or even rib-suppressing Hox genes within the dorsal region is inconsistent with the homogeneity in vertebral identity. Our results suggest that the evolution of a deregionalized, snake-like body involved not only alterations in Hox gene cis-regulation but also a different downstream interpretation of the Hox code.

Analysis of the sequence and phenotype of Drosophila Sex combs reduced alleles reveals potential functions of conserved protein motifs of the Sex combs reduced protein

The Drosophila Hox gene, Sex combs reduced (Scr), is required for patterning the larval and adult, labial and prothoracic segments. Fifteen Scr alleles were sequenced and the phenotypes analyzed in detail. Six null alleles were nonsense mutations (Scr(2), Scr(4), Scr(11), Scr(13), Scr(13A), and Scr(16)) and one was an intragenic deletion (Scr(17)). Five hypomorphic alleles were missense mutations (Scr(1), Scr(3), Scr(5), Scr(6), and Scr(8)) and one was a small protein deletion (Scr(15)). Protein sequence changes were found in four of the five highly conserved domains of SCR: the DYTQL motif (Scr(15)), YPWM motif (Scr(3)), Homeodomain (Scr(1)), and C-terminal domain (CTD) (Scr(6)), indicating importance for SCR function. Analysis of the pleiotropy of viable Scr alleles for the formation of pseudotracheae suggests that the DYTQL motif and the CTD mediate a genetic interaction with proboscipedia. One allele Scr(14), a missense allele in the conserved octapeptide, was an antimorphic allele that exhibited three interesting genetic properties. First, Scr(14)/Df had the same phenotype as Scr(+)/Df. Second, the ability of the Scr(14) allele to interact intragenetically with Scr alleles mapped to the first 82 amino acids of SCR, which contains the octapeptide motif. Third, Scr(6), which has two missense changes in the CTD, did not interact genetically with Scr(14).

The evolution of hierarchical gene regulatory networks

Comparative developmental evidence indicates that reorganizations in developmental gene regulatory networks (GRNs) underlie evolutionary changes in animal morphology, including body plans. We argue here that the nature of the evolutionary alterations that arise from regulatory changes depends on the hierarchical position of the change within a GRN. This concept cannot be accomodated by microevolutionary nor macroevolutionary theory. It will soon be possible to investigate these ideas experimentally, by assessing the effects of GRN changes on morphological evolution.

Hox specificity unique roles for cofactors and collaborators

Hox proteins are well known for executing highly specific functions in vivo, but our understanding of the molecular mechanisms underlying gene regulation by these fascinating proteins has lagged behind. The premise of this review is that an understanding of gene regulation-by any transcription factor-requires the dissection of the cis-regulatory elements that they act upon. With this goal in mind, we review the concepts and ideas regarding gene regulation by Hox proteins and apply them to a curated list of directly regulated Hox cis-regulatory elements that have been validated in the literature. Our analysis of the Hox-binding sites within these elements suggests several emerging generalizations. We distinguish between Hox cofactors, proteins that bind DNA cooperatively with Hox proteins and thereby help with DNA-binding site selection, and Hox collaborators, proteins that bind in parallel to Hox-targeted cis-regulatory elements and dictate the sign and strength of gene regulation. Finally, we summarize insights that come from examining five X-ray crystal structures of Hox-cofactor-DNA complexes. Together, these analyses reveal an enormous amount of flexibility into how Hox proteins function to regulate gene expression, perhaps providing an explanation for why these factors have been central players in the evolution of morphological diversity in the animal kingdom.

Evolution of an insect-specific GROUCHO-interaction motif in the ENGRAILED selector protein

Animal morphology evolves through alterations in the genetic regulatory networks that control development. Regulatory connections are commonly added, subtracted, or modified via mutations in cis-regulatory elements, but several cases are also known where transcription factors have gained or lost activity-modulating peptide motifs. In order to better assess the role of novel transcription factor peptide motifs in evolution, we searched for synapomorphic motifs in the homeotic selectors of Drosophila melanogaster and related insects. Here, we describe an evolutionarily novel GROUCHO (GRO)-interaction motif in the ENGRAILED (EN) selector protein. This "ehIFRPF" motif is not homologous to the previously characterized "engrailed homology 1" (eh1) GRO-interaction motif of EN. This second motif is an insect-specific "WRPW"-type motif that has been maintained by purifying selection in at least the dipteran/lepidopteran lineage. We demonstrate that this motif contributes to in vivo repression of the wingless (wg) target gene and to interaction with GRO in vitro. The acquisition and conservation of this auxiliary peptide motif shows how the number and activity of short peptide motifs can evolve in transcription factors while existing regulatory functions are maintained.

Genome adaptation to chemical stress: clues from comparative transcriptomics in Saccharomyces cerevisiae and Candida glabrata

Comparative transcriptomics of Saccharomyces cerevisiae and Candida glabrata revealed a remarkable conservation of response to drug-induced stress, despite underlying differences in the regulatory networks.

Background

Recent technical and methodological advances have placed microbial models at the forefront of evolutionary and environmental genomics. To better understand the logic of genetic network evolution, we combined comparative transcriptomics, a differential clustering algorithm and promoter analyses in a study of the evolution of transcriptional networks responding to an antifungal agent in two yeast species: the free-living model organism Saccharomyces cerevisiae and the human pathogen Candida glabrata.

Results

We found that although the gene expression patterns characterizing the response to drugs were remarkably conserved between the two species, part of the underlying regulatory networks differed. In particular, the roles of the oxidative stress response transcription factors ScYap1p (in S. cerevisiae) and Cgap1p (in C. glabrata) had diverged. The sets of genes whose benomyl response depends on these factors are significantly different. Also, the DNA motifs targeted by ScYap1p and Cgap1p are differently represented in the promoters of these genes, suggesting that the DNA binding properties of the two proteins are slightly different. Experimental assays of ScYap1p and Cgap1p activities in vivo were in accordance with this last observation.

Conclusions

Based on these results and recently published data, we suggest that the robustness of environmental stress responses among related species contrasts with the rapid evolution of regulatory sequences, and depends on both the coevolution of transcription factor binding properties and the versatility of regulatory associations within transcriptional networks.

High nucleotide divergence in developmental regulatory genes contrasts with the structural elements of olfactory pathways in caenorhabditis

Almost all organismal function is controlled by pathways composed of interacting genetic components. The relationship between pathway structure and the evolution of individual pathway components is not completely understood. For the nematode Caenorhabditis elegans, chemosensory pathways regulate critical aspects of an individual's life history and development. To help understand how olfaction evolves in Caenorhabditis and to examine patterns of gene evolution within transduction pathways in general, we analyzed nucleotide variation within and between species across two well-characterized olfactory pathways, including regulatory genes controlling the fate of the cells in which the pathways are expressed. In agreement with previous studies, we found much higher levels of polymorphism within C. remanei than within the related species C. elegans and C. briggsae. There are significant differences in the rates of nucleotide evolution for genes across the two pathways but no particular association between evolutionary rate and gene position, suggesting that the evolution of functional pathways must be considered within the context of broader gene network structure. However, developmental regulatory genes show both higher levels of divergence and polymorphism than the structural genes of the pathway. These results show that, contrary to the emerging paradigm in the evolution of development, important structural changes can accumulate in transcription factors.

Ultraconserved coding regions outside the homeobox of mammalian Hox genes

Background

All bilaterian animals share a general genetic framework that controls the formation of their body structures, although their forms are highly diversified. The Hox genes that encode transcription factors play a central role in this framework. All Hox proteins contain a highly conserved homeodomain encoded by the homeobox motif, but the other regions are generally assumed to be less conserved. In this study, we used comparative genomic methods to infer possible functional elements in the coding regions of mammalian Hox genes.

Results

We identified a set of ultraconserved coding regions (UCRs) outside the homeobox of mammalian Hox genes. Here a UCR is defined as a region of at least 120 nucleotides without synonymous and nonsynonymous nucleotide substitutions among different orders of mammals. Further analysis has indicated that these UCRs occur only in placental mammals and they evolved apparently after the split of placental mammals from marsupials. Analysis of human SNP data suggests that these UCRs are maintained by strong purifying selection.

Conclusion

Although mammalian genomes are known to contain ultraconserved non-coding elements (UNEs), this paper seems to be the first to report the UCRs in protein coding genes. The extremely high degree of sequence conservation in non-homeobox regions suggests that they might have important roles for the functions of Hox genes. We speculate that UCRs have some gene regulatory functions possibly in relation to the development of the intra-uterus child-bearing system.

Adaptive changes in the transcription factor HoxA-11 are essential for the evolution of pregnancy in mammals

Evolutionary change in gene regulation can result from changes in cis-regulatory elements, leading to differences in the temporal and spatial expression of genes or in the coding region of transcription factors leading to novel functions or both. Although there is a growing body of evidence supporting the importance of cis-regulatory evolution, examples of protein-mediated evolution of novel developmental pathways have not been demonstrated. Here, we investigate the evolution of prolactin (PRL) expression in endometrial cells, which is essential for placentation/pregnancy in eutherian mammals and is a direct regulatory target of the transcription factor HoxA-11. Here, we show that (i) endometrial PRL expression is a derived feature of placental mammals, (ii) the PRL regulatory gene HoxA-11 experienced a period of strong positive selection in the stem-lineage of eutherian mammals, and (iii) only HoxA-11 proteins from placental mammals, including the reconstructed ancestral eutherian gene, are able to up-regulate PRL from the promoter used in endometrial cells. In contrast, HoxA-11 from the reconstructed therian ancestor, opossum, platypus, and chicken are unable to up-regulate PRL expression. These results demonstrate that the evolution of novel gene expression domains is not only mediated by the evolution of cis-regulatory elements but can also require evolutionary changes of transcription factor proteins themselves.

Hox and senseless antagonism functions as a molecular switch to regulate EGF secretion in the Drosophila PNS

Hox factors are key regulators of distinct cells, tissues, and organs along the body plan. However, little is known about how Hox factors regulate cell-specific gene expression to pattern diverse tissues. Here, we show an unexpected Hox transcriptional mechanism: the permissive regulation of EGF secretion, and thereby cell specification, by antagonizing the Senseless transcription factor in the peripheral nervous system. rhomboid expression in a subset of sensory cells stimulates EGF secretion to induce hepatocyte-like cell development. We identified a rhomboid enhancer that is active in these cells and show that an abdominal Hox complex directly competes with Senseless for enhancer binding, with the transcriptional outcome dependent upon their relative binding activities. Thus, Hox-Senseless antagonism forms a molecular switch that integrates neural and anterior-posterior positional information. As the vertebrate senseless homolog is essential for neural development as well as hematopoiesis, we propose Hox-Senseless antagonism will broadly control cell fate decisions.

An empirical test for branch-specific positive selection

The use of phylogenetic analysis to predict positive selection specific to human genes is complicated by the very close evolutionary relationship with our nearest extant primate relatives, chimpanzees. To assess the power and limitations inherent in use of maximum-likelihood (ML) analysis of codon substitution patterns in such recently diverged species, a series of simulations was performed to assess the impact of several parameters of the evolutionary model on prediction of human-specific positive selection, including branch length and d(N)/d(S) ratio. Parameters were varied across a range of values observed in alignments of 175 transcription factor (TF) genes that were sequenced in 12 primate species. The ML method largely lacks the power to detect positive selection that has occurred since the most recent common ancestor between humans and chimpanzees. An alternative null model was developed on the basis of gene-specific evaluation of the empirical distribution of ML results, using simulated neutrally evolving sequences. This empirical test provides greater sensitivity to detect lineage-specific positive selection in the context of recent evolutionary divergence.

Multiple intrinsically disordered sequences alter DNA binding by the homeodomain of the Drosophila hox protein ultrabithorax

During animal development, distinct tissues, organs, and appendages are specified through differential gene transcription by Hox transcription factors. However, the conserved Hox homeodomains bind DNA with high affinity yet low specificity. We have therefore explored the structure of the Drosophila melanogaster Hox protein Ultrabithorax and the impact of its nonhomeodomain regions on DNA binding properties. Computational and experimental approaches identified several conserved, intrinsically disordered regions outside the homeodomain of Ultrabithorax that impact DNA binding by the homeodomain. Full-length Ultrabithorax bound to target DNA 2.5-fold weaker than its isolated homeodomain. Using N-terminal and C-terminal deletion mutants, we demonstrate that the YPWM region and the disordered microexons (termed the I1 region) inhibit DNA binding approximately 2-fold, whereas the disordered I2 region inhibits homeodomain-DNA interaction a further approximately 40-fold. Binding is restored almost to homeodomain affinity by the mostly disordered N-terminal 174 amino acids (R region) in a length-dependent manner. Both the I2 and R regions contain portions of the activation domain, functionally linking DNA binding and transcription regulation. Given that (i) the I1 region and a portion of the R region alter homeodomain-DNA binding as a function of pH and (ii) an internal deletion within I1 increases Ultrabithorax-DNA affinity, I1 must directly impact homeodomain-DNA interaction energetics. However, I2 appears to indirectly affect DNA binding in a manner countered by the N terminus. The amino acid sequences of I2 and much of the I1 and R regions vary significantly among Ultrabithorax orthologues, potentially diversifying Hox-DNA interactions.

Evo-devo and an expanding evolutionary synthesis: a genetic theory of morphological evolution

Biologists have long sought to understand which genes and what kinds of changes in their sequences are responsible for the evolution of morphological diversity. Here, I outline eight principles derived from molecular and evolutionary developmental biology and review recent studies of species divergence that have led to a genetic theory of morphological evolution, which states that (1) form evolves largely by altering the expression of functionally conserved proteins, and (2) such changes largely occur through mutations in the cis-regulatory sequences of pleiotropic developmental regulatory loci and of the target genes within the vast networks they control.

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.

Context-dependent regulation of Hox protein functions by CK2 phosphorylation sites

Variations in Hox protein sequences and functions have been proposed to contribute to evolutionary changes in appendage shape and number in crustaceans and insects. One model is that insect Hox proteins of the Ultrabithorax (UBX) ortholog class evolved increased abilities to repress Distal-less (Dll) transcription and appendage development in part through the loss of serine and threonine residues in casein kinase 2 (CK2) phosphorylation sites. To explore this possibility, we constructed and tested the appendage repression function of chimeric proteins with insertions of different CK2 consensus sites or phosphomimetics of CK2 sites in C-terminal regions of Drosophila melanogaster UBX. Our results indicate that CK2 sites C-terminal to the homeodomain can inhibit the appendage repression functions of UBX proteins, but only in the context of specific amino acid sequences. Our results, combined with previous findings on evolutionary changes in Hox protein, suggest how intra-protein regulatory changes can diversify Hox protein function, and thus animal morphology.

Evolution of eukaryotic transcription circuits

The gradual modification of transcription circuits over evolutionary time scales is an important source of the diversity of life. Over the past decade, studies in animals have shown how seemingly small molecular changes in gene regulation can have large effects on morphology and physiology and how selective pressures can act on these changes. More recently, genome-wide studies, particularly those in single-cell yeasts, have uncovered evidence of extensive transcriptional rewiring, indicating that even closely related organisms regulate their genes using markedly different circuitries.

Heterotypy in the N-terminal region of growth/differentiation factor 5 (GDF5) mature protein during teleost evolution

Heterotypy is now recognized as a generative force in the formation of new proteins through modification of existing proteins. We report that heterotypy in the N-terminal region of the mature growth/differentiation factor 5 (GDF5) protein occurred during evolution of teleosts. N-terminal length variation of GDF5 was found among teleost interfamilies and interorders but not within teleost families or among tetrapods. We further show that increase of proline and glutamine to the N-terminal region of mature GDF5 occurred in Eurypterygii, the higher lineage of teleosts. Because the basic amino acids, believed to control diffusion, are conserved in this region across all species examined, we suggest that the N-terminal elongation of the mature GDF5 protein during evolution has altered the protein diffusion in Eurypterygii, leading to high concentrations of the protein in the joint of the pharyngeal skeleton, the location of cartilage formation during development.

The correlated evolution of Runx2 tandem repeats, transcriptional activity, and facial length in carnivora

To assess the ability of protein-coding mutations to contribute to subtle, inter-specific morphologic evolution, here, we test the hypothesis that mutations within the protein-coding region of runt-related transcription factor 2 (Runx2) have played a role in facial evolution in 30 species from a naturally evolving group, the mammalian order Carnivora. Consistent with this hypothesis, we find significant correlations between changes in Runx2 glutamine-alanine tandem-repeat ratio, and both Runx2 transcriptional activity and carnivoran facial length. Furthermore, we identify a potential evolutionary mechanism for the correlation between Runx2 tandem repeat ratio and facial length. Specifically, our results are consistent with the Runx2 tandem repeat system providing a flexible genetic mechanism to rapidly change the timing of ossification. These heterochronic changes, in turn, potentially act on existing allometric variation in carnivoran facial length to generate the disparity in adult facial lengths observed among carnivoran species. Our results suggest that despite potentially great pleiotropic effects, changes to the protein-coding regions of genes such as Runx2 do occur and have the potential to affect subtle morphologic evolution across a diverse array of species in naturally evolving lineages.

Regulatory evolution in proteins by turnover and lineage-specific changes of cyclin-dependent kinase consensus sites

Evolutionary change in gene regulation is a key mechanism underlying the genetic component of organismal diversity. Here, we study evolution of regulation at the posttranslational level by examining the evolution of cyclin-dependent kinase (CDK) consensus phosphorylation sites in the protein subunits of the pre-replicative complex (RC). The pre-RC, an assembly of proteins formed during an early stage of DNA replication, is believed to be regulated by CDKs throughout the animals and fungi. Interestingly, although orthologous pre-RC components often contain clusters of CDK consensus sites, the positions and numbers of sites do not seem conserved. By analyzing protein sequences from both distantly and closely related species, we confirm that consensus sites can turn over rapidly even when the local cluster of sites is preserved, consistent with the notion that precise positioning of phosphorylation events is not required for regulation. We also identify evolutionary changes in the clusters of sites and further examine one replication protein, Mcm3, where a cluster of consensus sites near a nucleocytoplasmic transport signal is confined to a specific lineage. We show that the presence or absence of the cluster of sites in different species is associated with differential regulation of the transport signal. These findings suggest that the CDK regulation of MCM nuclear localization was acquired in the lineage leading to Saccharomyces cerevisiae after the divergence with Candida albicans. Our results begin to explore the dynamics of regulatory evolution at the posttranslational level and show interesting similarities to recent observations of regulatory evolution at the level of transcription.

A unique Extradenticle recruitment mode in the Drosophila Hox protein Ultrabithorax

Hox transcription factors are essential for shaping body morphology in development and evolution. The control of Hox protein activity in part arises from interaction with the PBC class of partners, pre-B cell transcription factor (Pbx) proteins in vertebrates and Extradenticle (Exd) in Drosophila. Characterized interactions occur through a single mode, involving a short hexapeptide motif in the Hox protein. This apparent uniqueness in Hox-PBC interaction provides little mechanistic insight in how the same cofactors endow Hox proteins with specific and diverse activities. Here, we identify in the Drosophila Ultrabithorax (Ubx) protein a short motif responsible for an alternative mode of Exd recruitment. Together with previous reports, this finding highlights that the Hox protein Ubx has multiple ways to interact with the Exd cofactor and suggests that flexibility in Hox-PBC contacts contributes to specify and diversify Hox protein function.

Collaboration between Smads and a Hox protein in target gene repression

Hox proteins control the differentiation of serially iterated structures in arthropods and chordates by differentially regulating many target genes. It is yet unclear to what extent Hox target gene selection is dependent upon other regulatory factors and how these interactions might affect target gene activation or repression. We find that two Smad proteins, effectors of the Drosophila Dpp/TGF-β pathway, that are genetically required for the activation of the spalt (sal) gene in the wing,collaborate with the Hox protein Ultrabithorax (Ubx) to directly repress sal in the haltere. The repression of sal is integrated by a cis-regulatory element (CRE) through a remarkably conserved set of Smad binding sites flanked by Ubx binding sites. If the Ubx binding sites are relocated at a distance from the Smad binding sites, the proteins no longer collaborate to repress gene expression. These results support an emerging view of Hox proteins acting in collaboration with a much more diverse set of transcription factors than has generally been appreciated.

Ubx regulates differential enlargement and diversification of insect hind legs

Differential enlargement of hind (T3) legs represents one of the hallmarks of insect evolution. However, the actual mechanism(s) responsible are yet to be determined. To address this issue, we have now studied the molecular basis of T3 leg enlargement in Oncopeltus fasciatus (milkweed bug) and Acheta domesticus (house cricket). In Oncopeltus, the T3 tibia displays a moderate increase in size, whereas in Acheta, the T3 femur, tibia, and tarsus are all greatly enlarged. Here, we show that the hox gene Ultrabithorax (Ubx) is expressed in the enlarged segments of hind legs. Furthermore, we demonstrate that depletion of Ubx during embryogenesis has a primary effect in T3 legs and causes shortening of leg segments that are enlarged in a wild type. This result shows that Ubx is regulating the differential growth and enlargement of T3 legs in both Oncopeltus and Acheta. The emerging view suggests that Ubx was co-opted for a novel role in regulating leg growth and that the transcriptional modification of its expression may be a universal mechanism for the evolutionary diversification of insect hind legs.

Genome clashes in hybrids: insights from gene expression

In interspecific hybrids, novel phenotypes often emerge from the interaction of two divergent genomes. Interactions between the two transcriptional networks are assumed to contribute to these unpredicted new phenotypes by inducing novel patterns of gene expression. Here we provide a review of the recent literature on the accumulation of regulatory incompatibilities. We review specific examples of regulatory incompatibilities reported at particular loci as well as genome-scale surveys of gene expression in interspecific hybrids. Finally, we consider and preview novel technologies that could help decipher how divergent transcriptional networks interact in hybrids between species.

The new mutation theory of phenotypic evolution

Recent studies of developmental biology have shown that the genes controlling phenotypic characters expressed in the early stage of development are highly conserved and that recent evolutionary changes have occurred primarily in the characters expressed in later stages of development. Even the genes controlling the latter characters are generally conserved, but there is a large component of neutral or nearly neutral genetic variation within and between closely related species. Phenotypic evolution occurs primarily by mutation of genes that interact with one another in the developmental process. The enormous amount of phenotypic diversity among different phyla or classes of organisms is a product of accumulation of novel mutations and their conservation that have facilitated adaptation to different environments. Novel mutations may be incorporated into the genome by natural selection (elimination of preexisting genotypes) or by random processes such as genetic and genomic drift. However, once the mutations are incorporated into the genome, they may generate developmental constraints that will affect the future direction of phenotypic evolution. It appears that the driving force of phenotypic evolution is mutation, and natural selection is of secondary importance.

Polyglutamine variation in a flowering time protein correlates with island age in a Hawaiian plant radiation

Background

A controversial topic in evolutionary developmental biology is whether morphological diversification in natural populations can be driven by expansions and contractions of amino acid repeats in proteins. To promote adaptation, selection on protein length variation must overcome deleterious effects of multiple correlated traits (pleiotropy). Thus far, systems that demonstrate this capacity include only ancient or artificial morphological diversifications. The Hawaiian Islands, with their linear geological sequence, present a unique environment to study recent, natural radiations. We have focused our research on the Hawaiian endemic mints (Lamiaceae), a large and diverse lineage with paradoxically low genetic variation, in order to test whether a direct relationship between coding-sequence repeat diversity and morphological change can be observed in an actively evolving system.

Results

Here we show that in the Hawaiian mints, extensive polyglutamine (CAG codon repeat) polymorphism within a homolog of the pleiotropic flowering time protein and abscisic acid receptor FCA tracks the natural environmental cline of the island chain, consequent with island age, across a period of 5 million years. CAG expansions, perhaps following their natural tendency to elongate, are more frequent in colonists of recently-formed, nutrient-rich islands than in their forebears on older, nutrient-poor islands. Values for several quantitative morphological variables related to reproductive investment, known from Arabidopsis fca mutant studies, weakly though positively correlate with increasing glutamine tract length. Together with protein modeling of FCA, which indicates that longer polyglutamine tracts could induce suboptimally mobile functional domains, we suggest that CAG expansions may form slightly deleterious alleles (with respect to protein function) that become fixed in founder populations.

Conclusion

In the Hawaiian mint FCA system, we infer that contraction of slightly deleterious CAG repeats occurred because of competition for resources along the natural environmental cline of the island chain. The observed geographical structure of FCA variation and its correlation with morphologies expected from Arabidopsis mutant studies may indicate that developmental pleiotropy played a role in the diversification of the mints. This discovery is important in that it concurs with other suggestions that repetitive amino acid motifs might provide a mechanism for driving morphological evolution, and that variation at such motifs might permit rapid tuning to environmental change.

Direct interaction between Teashirt and Sex combs reduced proteins, via Tsh's acidic domain, is essential for specifying the identity of the prothorax in Drosophila

teashirt (tsh) encodes a zinc-finger protein that is thought to be part of a network that contributes to regionalization of the Drosophila embryo and establishes the domains of Hox protein function. tsh and the Hox gene Sex combs reduced (Scr) are essential to establish the identity of the first thoracic segment. We used the development of the first thoracic segment as a paradigm for Scr dependent regional morphological distinctions. In this specific context, we asked whether Tsh protein could have a direct influence on Scr activity. Here we present evidence that Tsh interacts directly with Scr and this interaction depends in part on the presence of a short domain located in the N-terminal half of Teashirt called "acidic domain". In vivo, expression of full length Tsh can rescue the tsh null phenotype throughout the trunk whereas Tsh lacking the Scr interacting domain rescues all the trunk defects except in the prothorax. We suggest this provides insights into the mechanism by which Tsh, in concert with Scr, specifies the prothoracic identity.

The cnidarian-bilaterian ancestor possessed at least 56 homeoboxes: evidence from the starlet sea anemone, Nematostella vectensis

BACKGROUND

Homeodomain transcription factors are key components in the developmental toolkits of animals. While this gene superclass predates the evolutionary split between animals, plants, and fungi, many homeobox genes appear unique to animals. The origin of particular homeobox genes may, therefore, be associated with the evolution of particular animal traits. Here we report the first near-complete set of homeodomains from a basal (diploblastic) animal.

RESULTS

Phylogenetic analyses were performed on 130 homeodomains from the sequenced genome of the sea anemone Nematostella vectensis along with 228 homeodomains from human and 97 homeodomains from Drosophila. The Nematostella homeodomains appear to be distributed among established homeodomain classes in the following fashion: 72 ANTP class; one HNF class; four LIM class; five POU class; 33 PRD class; five SINE class; and six TALE class. For four of the Nematostella homeodomains, there is disagreement between neighbor-joining and Bayesian trees regarding their class membership. A putative Nematostella CUT class gene is also identified.

CONCLUSION

The homeodomain superclass underwent extensive radiations prior to the evolutionary split between Cnidaria and Bilateria. Fifty-six homeodomain families found in human and/or fruit fly are also found in Nematostella, though seventeen families shared by human and fly appear absent in Nematostella. Homeodomain loss is also apparent in the bilaterian taxa: eight homeodomain families shared by Drosophila and Nematostella appear absent from human (CG13424, EMXLX, HOMEOBRAIN, MSXLX, NK7, REPO, ROUGH, and UNC4), and six homeodomain families shared by human and Nematostella appear absent from fruit fly (ALX, DMBX, DUX, HNF, POU1, and VAX).