Functional evolution of the Ultrabithorax protein

Evidence of positive selection on six spider developmental genes

Spiders constitute more than 49,000 described species distributed all over the world, and all ecological environments. Their order, Araneae, is defined by a set of characteristics with no parallel among their arachnid counterparts (e.g., spinnerets, silk glands, chelicerae that inoculate venom, among others). Changes in developmental pathways often underlie the evolution of morphological synapomorphies, and as such spiders are a promising model to study the role of developmental genes in the origin of evolutionary novelties. With that in mind, we investigated changes in the evolutionary regime of a set of six developmental genes, using spiders as our model. The genes were mainly chosen for their roles in spinneret ontogeny, yet they are pleiotropic, and it is likely that the origins of other unique morphological phenotypes are also linked to changes in their sequences. Our results indicate no great differences in the selective pressures on those genes when comparing spiders to other arachnids, but a few site-specific positive selection evidence were found in the Araneae lineage. These findings lead us to new insights on spider evolution that are to be further tested.

Ultrabithorax is a key regulator for the dimorphism of wings, a main cause for the outbreak of planthoppers in rice

Rice planthoppers, the most devastating rice pests, occur in two wing forms: the short-wing form for rapid population growth and long-wing form for long-distance migration, which together create the mechanism for outbreak. Here we show that Ultrabithorax (Ubx) is a key regulator for switching between the long- and short-wing forms of rice planthoppers. Ubx is expressed in both forewing and hindwing pads, which is different from the canonical model of Ubx expression. In brown planthoppers, expression of Ubx (NlUbx) is regulated by nutritional status of the rice host. High-quality young plants induce NlUbx expression leading to the short-wing form; low-quality ripe plants reduce NlUbx expression resulting in long-wing form. We also showed that NlUbx is regulated by the insulin receptors NlInR1 and NlInR2. The default expression of NlInR1 inhibits NlUbx resulting in long-wings, while high-quality hosts induce NlInR2 expression, which represses NlInR1 thus promoting NlUbx expression to produce short-wings.

A comparative genomic analysis of targets of Hox protein Ultrabithorax amongst distant insect species

In the fruitfly Drosophila melanogaster, the differential development of wing and haltere is dependent on the function of the Hox protein Ultrabithorax (Ubx). Here we compare Ubx-mediated regulation of wing patterning genes between the honeybee, Apis mellifera, the silkmoth, Bombyx mori and Drosophila. Orthologues of Ubx are expressed in the third thoracic segment of Apis and Bombyx, although they make functional hindwings. When over-expressed in transgenic Drosophila, Ubx derived from Apis or Bombyx could suppress wing development, suggesting evolutionary changes at the level of co-factors and/or targets of Ubx. To gain further insights into such events, we identified direct targets of Ubx from Apis and Bombyx by ChIP-seq and compared them with those of Drosophila. While majority of the putative targets of Ubx are species-specific, a considerable number of wing-patterning genes are retained, over the past 300 millions years, as targets in all the three species. Interestingly, many of these are differentially expressed only between wing and haltere in Drosophila but not between forewing and hindwing in Apis or Bombyx. Detailed bioinformatics and experimental validation of enhancer sequences suggest that, perhaps along with other factors, changes in the cis-regulatory sequences of earlier targets contribute to diversity in Ubx function.

Mechanisms of Specificity for Hox Factor Activity

Metazoans encode clusters of paralogous Hox genes that are critical for proper development of the body plan. However, there are a number of unresolved issues regarding how paralogous Hox factors achieve specificity to control distinct cell fates. First, how do Hox paralogs, which have very similar DNA binding preferences in vitro, drive different transcriptional programs in vivo? Second, the number of potential Hox binding sites within the genome is vast compared to the number of sites bound. Hence, what determines where in the genome Hox factors bind? Third, what determines whether a Hox factor will activate or repress a specific target gene? Here, we review the current evidence that is beginning to shed light onto these questions. In particular, we highlight how cooperative interactions with other transcription factors (especially PBC and HMP proteins) and the sequences of cis-regulatory modules provide a basis for the mechanisms of Hox specificity. We conclude by integrating a number of the concepts described throughout the review in a case study of a highly interrogated Drosophila cis-regulatory module named “The Distal-less Conserved Regulatory Element” (DCRE).

Rethinking gene regulatory networks in light of alternative splicing, intrinsically disordered protein domains, and post-translational modifications

Models for genetic regulation and cell fate specification characteristically assume that gene regulatory networks (GRNs) are essentially deterministic and exhibit multiple stable states specifying alternative, but pre-figured cell fates. Mounting evidence shows, however, that most eukaryotic precursor RNAs undergo alternative splicing (AS) and that the majority of transcription factors contain intrinsically disordered protein (IDP) domains whose functionalities are context dependent as well as subject to post-translational modification (PTM). Consequently, many transcription factors do not have fixed cis-acting regulatory targets, and developmental determination by GRNs alone is untenable. Modeling these phenomena requires a multi-scale approach to explain how GRNs operationally interact with the intra- and intercellular environments. Evidence shows that AS, IDP, and PTM complicate gene expression and act synergistically to facilitate and promote time- and cell-specific protein modifications involved in cell signaling and cell fate specification and thereby disrupt a strict deterministic GRN-phenotype mapping. The combined effects of AS, IDP, and PTM give proteomes physiological plasticity, adaptive responsiveness, and developmental versatility without inefficiently expanding genome size. They also help us understand how protein functionalities can undergo major evolutionary changes by buffering mutational consequences.

Ras-activated Dsor1 promotes Wnt signaling in Drosophila development

Wnt/Wingless (Wg) and Ras/MAPK signaling both play fundamental roles in growth, cell-fate determination, and when dysregulated, can lead to tumorigenesis. Several conflicting modes of interaction between Ras/MAPK and Wnt signaling have been identified in specific cellular contexts, causing synergistic or antagonistic effects on target genes. We find novel evidence that the dual specificity kinase MEK, Downstream of Raf1 (Dsor1), is required for Wnt signaling. Knockdown of Dsor1 results in loss of Wg target gene expression, as well as reductions in stabilized Armadillo (Arm; Drosophila β-catenin). We have identified a close physical interaction between Dsor1 and Arm, and find that catalytically inactive Dsor1 causes a reduction inactive Arm. These results suggest that Dsor1 normally counteracts the Axin-mediated destruction of Arm. We find that Ras-Dsor1 activity is independent of upstream activation by EGFR, rather it appears to be activated by the insulin-like growth factor receptor to promote Wg signaling. Together our results suggest novel crosstalk between Insulin and Wg signaling via Dsor1.

The intrinsically disordered regions of the Drosophila melanogaster Hox protein ultrabithorax select interacting proteins based on partner topology

Interactions between structured proteins require a complementary topology and surface chemistry to form sufficient contacts for stable binding. However, approximately one third of protein interactions are estimated to involve intrinsically disordered regions of proteins. The dynamic nature of disordered regions before and, in some cases, after binding calls into question the role of partner topology in forming protein interactions. To understand how intrinsically disordered proteins identify the correct interacting partner proteins, we evaluated interactions formed by the Drosophila melanogaster Hox transcription factor Ultrabithorax (Ubx), which contains both structured and disordered regions. Ubx binding proteins are enriched in specific folds: 23 of its 39 partners include one of 7 folds, out of the 1195 folds recognized by SCOP. For the proteins harboring the two most populated folds, DNA-RNA binding 3-helical bundles and α-α superhelices, the regions of the partner proteins that exhibit these preferred folds are sufficient for Ubx binding. Three disorder-containing regions in Ubx are required to bind these partners. These regions are either alternatively spliced or multiply phosphorylated, providing a mechanism for cellular processes to regulate Ubx-partner interactions. Indeed, partner topology correlates with the ability of individual partner proteins to bind Ubx spliceoforms. Partners bind different disordered regions within Ubx to varying extents, creating the potential for competition between partners and cooperative binding by partners. The ability of partners to bind regions of Ubx that activate transcription and regulate DNA binding provides a mechanism for partners to modulate transcription regulation by Ubx, and suggests that one role of disorder in Ubx is to coordinate multiple molecular functions in response to tissue-specific cues.

Differential activity of Drosophila Hox genes induces myosin expression and can maintain compartment boundaries

Compartments are units of cell lineage that subdivide territories with different developmental potential. In Drosophila, the wing and haltere discs are subdivided into anterior and posterior (A/P) compartments, which require the activity of Hedgehog, and into dorsal and ventral (D/V) compartments, needing Notch signaling. There is enrichment in actomyosin proteins at the compartment boundaries, suggesting a role for these proteins in their maintenance. Compartments also develop in the mouse hindbrain rhombomeres, which are characterized by the expression of different Hox genes, a group of genes specifying different structures along their main axis of bilaterians. We show here that the Drosophila Hox gene Ultrabithorax can maintain the A/P and D/V compartment boundaries when Hedgehog or Notch signaling is compromised, and that the interaction of cells with and without Ultrabithorax expression induces high levels of non-muscle myosin II. In the absence of Ultrabithorax there is occasional mixing of cells from different segments. We also show a similar role in cell segregation for the Abdominal-B Hox gene. Our results suggest that the juxtaposition of cells with different Hox gene expression leads to their sorting out, probably through the accumulation of non-muscle myosin II at the boundary of the different cell territories. The increase in myosin expression seems to be a general mechanism used by Hox genes or signaling pathways to maintain the segregation of different groups of cells.

The phylotypic stage as a boundary of modular memory: non mechanistic perspective

The concept of the phylotypic stage has been strongly integrated into developmental biology, thanks mostly to drawings presented by Haeckel (Anthropogenie oder Entwicklungsgeschichte des Menschen, 1874). They are printed in every textbook as proof of the existence of the phylotypic stage and the fact of its conservation, albeit many times criticized as misleading and simplifying (Richardson in Develop Biol 172:412-421, 1995, Richardson et al. in Anat Embryo 196:91-106, 1997; Bininda-Emons et al. in Proc R Soc Lond 270:341-346, 2003). Although generally accepted by modern biology, doubt still exists concerning the very existence or the usefulness of the concept. What kind of evolutionary and developmental horizons does it open indeed? This article begins with the history of the concept, discusses its validity and draws this into connotation with the idea of a memory activated throughout the development. Barbieri (The organic codes. An introduction to semantic biology, 2003) considers the phylotypic stage to be a crucial boundary when the genetic program ceases to suffice for further development of the embryo, and supracellular memory of the body plan is activated. This moment clearly coincides with the commencing of the modular development of the embryo. In this article the nature of such putative memory will be discussed.

Surprising flexibility in a conserved Hox transcription factor over 550 million years of evolution

Although metazoan body plans are remarkably diverse, the structure and function of many embryonic regulatory genes are conserved because large changes would be detrimental to development. However, the fushi tarazu (ftz) gene has changed dramatically during arthropod evolution from Hox-like to a pair-rule segmentation gene in Drosophila. Changes in both expression and protein sequence contributed to this new function: ftz expression switched from Hox-like to stripes and changes in Ftz cofactor interaction motifs led to loss of homeotic and gain of segmentation potential. Here, we reconstructed ftz changes in a rigorous phylogenetic context. We found that ftz did not simply switch from Hox-like to segmentation function; rather, ftz is remarkably labile, having undergone multiple changes in sequence and expression. The segmentation LXXLL motif was stably acquired in holometabolous insects after the appearance of striped expression in early insect lineages. The homeotic YPWM motif independently degenerated multiple times. These "degen-YPWMs" showed varying degrees of homeotic potential when expressed in Drosophila, suggesting variable loss of Hox function in different arthropods. Finally, the intensity of ftz Hox-like expression decreased to marginal levels in some crustaceans. We propose that decreased expression levels permitted ftz variants to arise and persist in populations without disadvantaging organismal development. This process, in turn, allowed evolutionary transitions in protein function, as weakly expressed "hopeful gene variants" were coopted into alternative developmental pathways. Our findings show that variation of a pleiotropic transcription factor is more extensive than previously imagined, suggesting that evolutionary plasticity may be widespread among regulatory genes.

Use with caution: developmental systems divergence and potential pitfalls of animal models

Transgenic animal models have played an important role in elucidating gene functions and the molecular basis development, physiology, behavior, and pathogenesis. Transgenic models have been so successful that they have become a standard tool in molecular genetics and biomedical studies and are being used to fulfill one of the main goals of the post-genomic era: to assign functions to each gene in the genome. However, the assumption that gene functions and genetic systems are conserved between models and humans is taken for granted, often in spite of evidence that gene functions and networks diverge during evolution. In this review, I discuss some mechanisms that generate functional divergence and highlight recent examples demonstrating that gene functions and regulatory networks diverge through time. These examples suggest that annotation of gene functions based solely on mutant phenotypes in animal models, as well as assumptions of conserved functions between species, can be wrong. Therefore, animal models of gene function and human disease may not provide appropriate information, particularly for rapidly evolving genes and systems.

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.

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.

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.

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.

The road to modularity

A network of interactions is called modular if it is subdivided into relatively autonomous, internally highly connected components. Modularity has emerged as a rallying point for research in developmental and evolutionary biology (and specifically evo-devo), as well as in molecular systems biology. Here we review the evidence for modularity and models about its origin. Although there is an emerging agreement that organisms have a modular organization, the main open problem is the question of whether modules arise through the action of natural selection or because of biased mutational mechanisms.

The tardigrade Hypsibius dujardini, a new model for studying the evolution of development

Studying development in diverse taxa can address a central issue in evolutionary biology: how morphological diversity arises through the evolution of developmental mechanisms. Two of the best-studied developmental model organisms, the arthropod Drosophila and the nematode Caenorhabditis elegans, have been found to belong to a single protostome superclade, the Ecdysozoa. This finding suggests that a closely related ecdysozoan phylum could serve as a valuable model for studying how developmental mechanisms evolve in ways that can produce diverse body plans. Tardigrades, also called water bears, make up a phylum of microscopic ecdysozoan animals. Tardigrades share many characteristics with C. elegans and Drosophila that could make them useful laboratory models, but long-term culturing of tardigrades historically has been a challenge, and there have been few studies of tardigrade development. Here, we show that the tardigrade Hypsibius dujardini can be cultured continuously for decades and can be cryopreserved. We report that H. dujardini has a compact genome, a little smaller than that of C. elegans or Drosophila, and that sequence evolution has occurred at a typical rate. H. dujardini has a short generation time, 13-14 days at room temperature. We have found that the embryos of H. dujardini have a stereotyped cleavage pattern with asymmetric cell divisions, nuclear migrations, and cell migrations occurring in reproducible patterns. We present a cell lineage of the early embryo and an embryonic staging series. We expect that these data can serve as a platform for using H. dujardini as a model for studying the evolution of developmental mechanisms.

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.

Lessons from leeches: a call for DNA barcoding in the lab

Many evolution of development labs study organisms that must be periodically collected from the wild. Whenever this is the case, there is the risk that different field collections will recover genetically different strains or cryptic species. Ignoring this potential for genetic variation may introduce an uncontrolled source of experimental variability, leading to confusion or misinterpretation of the results. Leeches in the genus Helobdella have been a workhorse of annelid developmental biology for 30 years. Nearly all early Helobdella research was based on a single isolate, but in recent years isolates from multiple field collections and multiple sites across the country have been used. To assess the genetic distinctness of different isolates, we obtained specimens from most Helobdella laboratory cultures currently or recently in use and from some of their source field sites. From these samples, we sequenced part of the mitochondrial gene cytochrome oxidase I (COI). Sequence divergences and phylogenetic analyses reveal that, collectively, the Helobdella development community has worked on five distinct species from two major clades. Morphologically similar isolates that were thought to represent the same species (H. robusta) actually represent three species, two of which coexist at the same locality. Another isolate represents part of a species complex (the "H. triserialis" complex), and yet another is an invasive species (H. europaea). We caution researchers similarly working on multiple wild-collected isolates to preserve voucher specimens and to obtain from these a molecular "barcode," such as a COI gene sequence, to reveal genetic variation in animals used for research.

Adaptive evolution of Hox-gene homeodomains after cluster duplications

BACKGROUND

Hox genes code for homeodomain-containing transcription factors that function in cell fate determination and embryonic development. Hox genes are arranged in clusters with up to 14 genes. This archetypical chordate cluster has duplicated several times in vertebrates, once at the origin of vertebrates and once at the origin of gnathostoms, an additional duplication event is associated with the origin of teleosts and the agnanths, suggesting that duplicated Hox cluster genes are involved in the genetic mechanisms behind the diversification of vertebrate body plans, and the origin of morphological novelties. Preservation of duplicate genes is promoted by functional divergence of paralogs, either by subfunction partitioning among paralogs or the acquisition of a novel function by one paralog. But for Hox genes the mechanisms of paralog divergence is unknown, leaving open the role of Hox gene duplication in morphological evolution.

RESULTS

Here, we use several complementary methods, including branch-specific dN/dS ratio tests, branch-site dN/dS ratio tests, clade level amino acid conservation/variation patterns, and relative rate ratio tests, to show that the homeodomain of Hox genes was under positive Darwinian selection after cluster duplications.

CONCLUSION

Our results suggest that positive selection acted on the homeodomain immediately after Hox clusters duplications. The location of sites under positive selection in the homeodomain suggests that they are involved in protein-protein interactions. These results further suggest that adaptive evolution actively contributed to Hox-gene homeodomain functions.

Identification of quantitative trait loci for larval morphological traits in interspecific hybrids of Ochlerotatus triseriatus and Ochlerotatus hendersoni (Diptera: Culicidae)

Ochlerotatus triseriatus is the natural vector of La Crosse virus, a common cause of pediatric encephalitis in the United States; the closely related Ochlerotatus hendersoni transmits this virus at low frequency. Adults of these mosquito species are difficult to distinguish morphologically; however, the larval stages show species-specific differences in several characters. We identified genomic regions contributing to the differences between the larvae of these species through interspecific hybridizations. Quantitative trait loci (QTL) were identified by standard interval mapping techniques and by univariate marker association analyses. We examined 159 F(2) progeny from an Oc. hendersoni female by Oc. triseriatus male interspecific cross for variation in the number of saddle and siphon hair branches, attachment of the acus, and morphology of the anal papillae. At least one putative QTL was identified for each of the phenotypes examined. QTL most commonly mapped to linkage group (LG) III, although QTL were identified on LGI and LGII for three phenotypes each. Several of these QTL, and particularly those on LGIII, also map to genome regions controlling adult female body size and ability to orally transmit La Crosse virus. Further studies are required to elucidate the relationships among these traits and the impact they may have had on the ecological specialization and speciation of these mosquitoes.

Modulation of Decapentaplegic gradient during haltere specification in Drosophila

Suppression of wing fate and specification of haltere fate in Drosophila by Ultrabithorax is a classical example of Hox regulation of serial homology. However, the mechanism of Ultrabithorax function in specifying haltere size and shape is not well understood. Here we show that Decapentaplegic signaling, which controls wing growth and shape, is a target of Ultrabithorax function during haltere specification. The Decapentaplegic signaling is down-regulated in haltere discs due to a combination of reduced levels of the Dpp, its trapping at the A/P boundary by increased levels of its receptor Thick-vein and its inability to diffuse in the absence of Dally. Results presented here suggest a complex mechanism adopted by Ultrabithorax to modulate Decapentaplegic signaling. We discuss how this complexity may regulate the final form of the adult haltere in the fly, without compromising features such as cell survival, which is also dependent on Decapentaplegic signaling.

Modulation of AP and DV signaling pathways by the homeotic gene Ultrabithorax during haltere development in Drosophila

Suppression of wing fate and specification of haltere fate in Drosophila by the homeotic gene Ultrabithorax is a classical example of Hox regulation of serial homology (Lewis, E.B. 1978. Nature 276, 565-570) and has served as a paradigm for understanding homeotic gene function. We have used DNA microarray analyses to identify potential targets of Ultrabithorax function during haltere specification. Expression patterns of 18 validated target genes and functional analyses of a subset of these genes suggest that down-regulation of both anterior-posterior and dorso-ventral signaling is critical for haltere fate specification. This is further confirmed by the observation that combined over-expression of Decapentaplegic and Vestigial is sufficient to override the effect of Ubx and cause dramatic haltere-to-wing transformations. Our results also demonstrate that analysis of the differential development of wing and haltere is a good assay system to identify novel regulators of key signaling pathways.

Evolutionarily conserved domains required for activation and repression functions of the Drosophila Hox protein Ultrabithorax

While testing the functions of deletion mutants in the Hox protein Ultrabithorax (Ubx), we found that the embryonic repression function of Ubx on Distal-less transcription in limb primordia is highly concentration dependent. The steep sigmoidal relationship between in vivo Ubx concentration and Distal-less repression is dependent on the Ubx YPWM motif. This suggests that Ubx cooperatively assembles a multi-protein repression complex on Distal-less regulatory DNA with the YPWM motif as a key protein-protein interface in this complex. Our deletion mutants also provide evidence for a transcriptional activation domain in the N-terminal 19 amino acids of Ubx. This proposed activation domain contains a variant of the SSYF motif that is found at the N termini of many Hox proteins, and is conserved in the activation domain of another Hox protein, Sex combs reduced. These results suggest that the N-terminal region containing the SSYF motif has been conserved in many Hox proteins for its role in transcriptional activation.

Physical and genetic interactions link hox function with diverse transcription factors and cell signaling proteins

Positional information provided by Hox homeotic transcription factors is integrated with other transcription factors and cell signaling cascades in specific combinations to dictate context- and gene-specific Hox activity. Protein-protein interactions between these groups have long been hypothesized to modulate Hox functions, yielding a context-specific function. However, difficulties in applying interaction screens to potent transcription factors have limited partner identification. A yeast two-hybrid screen using transcription activation-deficient mutants of the Drosophila melanogaster Hox protein Ultrabithorax IB identified an array of interacting proteins, consisting primarily of transcription factors and components of cell signaling pathways. Interactions were confirmed with wild-type Ultrabithorax (UBX) in phage display experiments and by immunoprecipitation for a subset of partners. In vivo assays demonstrated that two Ultrabithorax IB partners, Armadillo, regulated by Wingless/WNT signaling, and the homeodomain protein Aristaless, inhibit UBX-dependent haltere development from the default wing development pathway. Therefore, transcription factors and cell signaling proteins that subdivide Hox-specified tissues can both alter Hox function in vivo and interact with the corresponding Hox protein in vitro. UBX may also modulate partner function: the pupal death phenotype induced by ectopic expression of the UBX partner Hairy required the presence of UBX. Thus, Hox.transcription factor complexes may integrate a variety of positional cues, generating the specificity and versatility required for context-dependent Hox function.

Molecular evolution of evolutionary novelties: the vagina and uterus of therian mammals

Innovations are an integral part of the evolutionary process if we accept the fact that more complex organisms derived from anatomically simple ones. All major taxa are distinguished not only by their closer genealogical relatedness relative to other species but also by the possession of novel anatomical and physiological features. The question is whether the origin of these novel characters can be simply understood as adaptations, like all other phenotypic differences that arise by natural selection, or whether the origin of these characters requires more profound genetic changes. In this paper, we argue that innovations constitute a distinct class of evolutionary processes that require a research program complementary to the study of adaptation. The distinguishing feature of innovations is the origin of novel organ identity gene functions specific to the novel character. By implication, research into the origin of novel characters has to identify the developmental regulatory links that were involved in the evolution of these characters. We suggest that novel regulatory links will include the evolution of cis-regulatory elements as well as novel protein-protein interactions among transcription factor proteins. The latter hypothesis suggests that innovations should leave a trace in the evolution of the protein coding regions of transcription factor genes. We illustrate this idea with results on the evolution of HoxA-11 and HoxA-13 in the stem lineage of placental mammals. These genes are essential for female reproductive tract development and function. We show that, as predicted, these genes experience strong directional selection in the stem lineage of placental mammals and that these amino acid substitutions affect residues at the surface of the protein, consistent with their expected role in protein-protein interactions. We conclude that a careful analysis of sequence variation in developmental genes can aid in testing which developmental changes were instrumental in the origin of novel morphological characters.

Pleiotropic functions of a conserved insect-specific Hox peptide motif

The proteins that regulate developmental processes in animals have generally been well conserved during evolution. A few cases are known where protein activities have functionally evolved. These rare examples raise the issue of how highly conserved regulatory proteins with many roles evolve new functions while maintaining old functions. We have investigated this by analyzing the function of the ;QA' peptide motif of the Hox protein Ultrabithorax (Ubx), a motif that has been conserved throughout insect evolution since its establishment early in the lineage. We precisely deleted the QA motif at the endogenous locus via allelic replacement in Drosophila melanogaster. Although the QA motif was originally characterized as involved in the repression of limb formation, we have found that it is highly pleiotropic. Curiously, deleting the QA motif had strong effects in some tissues while barely affecting others, suggesting that QA function is preferentially required for a subset of Ubx target genes. QA deletion homozygotes had a normal complement of limbs, but, at reduced doses of Ubx and the abdominal-A (abd-A) Hox gene, ectopic limb primordia and adult abdominal limbs formed when the QA motif was absent. These results show that redundancy and the additive contributions of activity-regulating peptide motifs play important roles in moderating the phenotypic consequences of Hox protein evolution, and that pleiotropic peptide motifs that contribute quantitatively to several functions are subject to intense purifying selection.

Wing-to-Leg homeosis by spineless causes apoptosis regulated by Fish-lips, a novel leucine-rich repeat transmembrane protein

Growth, patterning, and apoptosis are mutually interactive during development. For example, cells that select an abnormal fate in a developing field are frequently removed by apoptosis. An important issue in this process that needs to be resolved is the mechanism used by cells to discern their correct fate from an abnormal fate. In order to examine this issue, we developed an animal model that expresses the dioxin receptor homolog Spineless (Ss) ectopically in the Drosophila wing. The presence of mosaic clones ectopically expressing ss results in a local transformation of organ identity, homeosis, from wing into a leg or antenna. The cells with misspecified fates subsequently activate c-Jun N-terminal kinase to undergo apoptosis in an autonomous or nonautonomous manner depending on their position within the wing, suggesting that a cell-cell interaction is, at least in some cases, involved in the detection of misspecified cells. Similar position dependence is commonly observed when various homeotic genes controlling the body segments are ectopically expressed. The autonomous and nonautonomous apoptosis caused by ss is regulated by a novel leucine-rich repeat family transmembrane protein, Fish-lips (Fili) that interacts with surrounding normal cells. These data support a mechanism in which the lack of some membrane proteins helps to recognize the presence of different cell types and direct these cells to an apoptotic fate in order to exclude them from the normal developing field.

Adaptive evolution of HoxA-11 and HoxA-13 at the origin of the uterus in mammals

The evolution of morphological characters is mediated by the evolution of developmental genes. Evolutionary changes can either affect cis-regulatory elements, leading to differences in their temporal and spatial regulation, or affect the coding region. Although there is ample evidence for the importance of cis-regulatory evolution, it has only recently been shown that transcription factors do not remain functionally equivalent during evolution. These results suggest that the evolution of transcription factors may play an active role in the evolution of development. To test this idea we investigated the molecular evolution of two genes essential for the development and function of the mammalian female reproductive organs, HoxA-11 and HoxA-13. We predicted that if coding-region evolution plays an active role in developmental evolution, then these genes should have experienced adaptive evolution at the origin of the mammalian female reproductive system. We report the sequences of HoxA-11 from basal mammalian and amniote taxa and analyse HoxA-11 and HoxA-13 for signatures of adaptive molecular evolution. The data demonstrate that these genes were under strong positive (directional) selection in the stem lineage of therian and eutherian mammals, coincident with the evolution of the uterus and vagina. These results support the idea that adaptive evolution of transcription factors can be an integral part in the evolution of novel structures.

HOXA13 regulates the expression of bone morphogenetic proteins 2 and 7 to control distal limb morphogenesis

In humans and mice, loss of HOXA13 function causes defects in the growth and patterning of the digits and interdigital tissues. Analysis of Hoxa13 expression reveals a pattern of localization overlapping with sites of reduced Bmp2 and Bmp7 expression in Hoxa13mutant limbs. Biochemical analyses identified a novel series of Bmp2and Bmp7 enhancer regions that directly interact with the HOXA13 DNA-binding domain and activate gene expression in the presence of HOXA13. Immunoprecipitation of HOXA13-Bmp2 and HOXA13-Bmp7 enhancer complexes from the developing autopod confirm that endogenous HOXA13 associates with these regions. Exogenous application of BMP2 or BMP7 partially rescues the Hoxa13 mutant limb phenotype, suggesting that decreased BMP signaling contributes to the malformations present in these tissues. Together, these results provide conclusive evidence that HOXA13 regulates Bmp2 and Bmp7 expression, providing a mechanistic link between HOXA13, its target genes and the specific developmental processes affected by loss of HOXA13 function.

Gene Network Polymorphism Is the Raw Material of Natural Selection: The Selfish Gene Network Hypothesis

Population genetics, the mathematical theory of modern evolutionary biology, defines evolution as the alteration of the frequency of distinct gene variants (alleles) differing in fitness over the time. The major problem with this view is that in gene and protein sequences we can find little evidence concerning the molecular basis of phenotypic variance, especially those that would confer adaptive benefit to the bearers. Some novel data, however, suggest that a large amount of genetic variation exists in the regulatory region of genes within populations. In addition, comparison of homologous DNA sequences of various species shows that evolution appears to depend more strongly on gene expression than on the genes themselves. Furthermore, it has been demonstrated in several systems that genes form functional networks, whose products exhibit interrelated expression profiles. Finally, it has been found that regulatory circuits of development behave as evolutionary units. These data demonstrate that our view of evolution calls for a new synthesis. In this article I propose a novel concept, termed the selfish gene network hypothesis, which is based on an overall consideration of the above findings. The major statements of this hypothesis are as follows. (1) Instead of individual genes, gene networks (GNs) are responsible for the determination of traits and behaviors. (2) The primary source of microevolution is the intraspecific polymorphism in GNs and not the allelic variation in either the coding or the regulatory sequences of individual genes. (3) GN polymorphism is generated by the variation in the regulatory regions of the component genes and not by the variance in their coding sequences. (4) Evolution proceeds through continuous restructuring of the composition of GNs rather than fixing of specific alleles or GN variants.

Hox Transcription Factor Ultrabithorax Ib Physically and Genetically Interacts with Disconnected Interacting Protein 1, a Double-stranded RNA-binding Protein

The Hox protein family consists of homeodomain-containing transcription factors that are primary determinants of cell fate during animal development. Specific Hox function appears to rely on protein-protein interactions; however, the partners involved in these interactions and their function are largely unknown. Disconnected Interacting Protein 1 (DIP1) was isolated in a yeast two-hybrid screen of a 0-12-h Drosophila embryo library designed to identify proteins that interact with Ultrabithorax (Ubx), a Drosophila Hox protein. The Ubx.DIP1 physical interaction was confirmed using phage display, immunoprecipitation, pull-down assays, and gel retardation analysis. Ectopic expression of DIP1 in wing and haltere imaginal discs malforms the adult structures and enhances a decreased Ubx expression phenotype, establishing a genetic interaction. Ubx can generate a ternary complex by simultaneously binding its target DNA and DIP1. A large region of Ubx, including the repression domain, is required for interaction with DIP1. These more variable sequences may be key to the differential Hox function observed in vivo. The Ubx.DIP1 interaction prevents transcriptional activation by Ubx in a modified yeast one-hybrid assay, suggesting that DIP1 may modulate transcriptional regulation by Ubx. The DIP1 sequence contains two dsRNA-binding domains, and DIP1 binds double-stranded RNA with a 1000-fold higher affinity than either single-stranded RNA or double-stranded DNA. The strong interaction of Ubx with an RNA-binding protein suggests a wider range of proteins may influence Ubx function than previously appreciated.

Differential expression patterns of the hox gene are associated with differential growth of insect hind legs

Diversification of leg appendages is one of the hallmarks of morphological evolution in insects. In particular, insect hind (T3) legs exhibit a whole spectrum of morphological diversification, ranging from uniform to extremely modified. To elucidate the developmental basis of T3 leg evolution, we have examined the expression patterns of two homeotic genes, Ultrabithorax and abdominal-A (collectively referred to as UbdA), in a broad range of species. First, our results show that UbdA expression in hemimetabolous insects is localized only in specific T3 leg segments undergoing differential growth (compared to their foreleg counterparts). In contrast, in basal hexapod and insect lineages, the absence of the UbdA signal coincides with uniform leg morphology. The same situation exists in first instar larvae of holometabolous insects, in which absence of UbdA expression in the embryonic T3 legs is associated with the lack of larval T3 leg diversification. Second, there is a clear difference in the timing of expression between species with greatly enlarged T3 leg, such as crickets and grasshoppers, and species that exhibit more moderate enlargement of hind legs, such as mantids and cockroaches. In the former, the UbdA expression starts much earlier, coinciding with the elongation of T3 limb buds. In the latter, however, the UbdA expression starts at much later stages of development, coinciding with the establishment of distinct leg segments. These results suggest that diversification of insect hind legs was influenced by changes in both the spatial and temporal regulation of the UbdA expression.

Bapx1 regulates patterning in the middle ear: altered regulatory role in the transition from the proximal jaw during vertebrate evolution

The middle ear apparatus is composed of three endochondrial ossicles (the stapes, incus and malleus) and two membranous bones, the tympanic ring and the gonium, which act as structural components to anchor the ossicles to the skull. Except for the stapes, these skeletal elements are unique to mammals and are derived from the first and second branchial arches. We show that, in combination with goosecoid (Gsc), the Bapx1 gene defines the structural components of the murine middle ear.

During embryogenesis, Bapx1 is expressed in a discrete domain within the mandibular component of the first branchial arch and later in the primordia of middle ear-associated bones, the gonium and tympanic ring. Consistent with the expression pattern of Bapx1, mouse embryos deficient for Bapx1 lack a gonium and display hypoplasia of the anterior end of the tympanic ring. At E10.5, expression of Bapx1partially overlaps that of Gsc and although Gsc is required for development of the entire tympanic ring, the role of Bapx1 is restricted to the specification of the gonium and the anterior tympanic ring. Thus, simple overlapping expression of these two genes appears to account for the patterning of the elements that compose the structural components of the middle ear and suggests that they act in concert.

In addition, Bapx1 is expressed both within and surrounding the incus and the malleus. Examination of the malleus shows that the width, but not the length, of this ossicle is decreased in the mutant mice. In non-mammalian jawed vertebrates, the bones homologous to the mammalian middle ear ossicles compose the proximal jaw bones that form the jaw articulation(primary jaw joint). In fish, Bapx1 is responsible for the formation of the joint between the quadrate and articular (homologues of the malleus and incus, respectively) enabling an evolutionary comparison of the role of a regulatory gene in the transition of the proximal jawbones to middle ear ossicles. Contrary to expectations, murine Bapx1 does not affect the articulation of the malleus and incus. We show that this change in role of Bapx1 following the transition to the mammalian ossicle configuration is not due to a change in expression pattern but results from an inability to regulate Gdf5 and Gdf6, two genes predicted to be essential in joint formation.

Molecular insights into Cumacean family relationships (Crustacea, Cumacea)

Cumaceans are a diverse order of small, benthic marine crustaceans. Phylogenetic hypotheses for the eight currently recognized cumacean families have not been formally proposed. However, based on external morphological traits and Linnean classification, a few conflicting hypotheses of relatedness have been proposed. Family definitions rely on morphological characters that often overlap and diagnoses are based on a combination of non-unique characters. Morphological analysis does not provide a well-resolved phylogeny. In the present study, we use amino acid sequences from the mitochondrial cytochrome oxidase I gene to produce a molecular phylogenetic hypothesis for the families of Cumacea. Phylogenetic analyses at the amino acid level were performed under Bayesian, likelihood, and parsimony methods. Results strongly suggest that families lacking an articulated telson form a monophyletic group. This pleotelson clade, composed of the families Bodotriidae, Leuconidae, and Nannastacidae, is the most derived within the Cumacea. Within this group, the Bodotriidae resolve paraphyletically, with Leuconidae and Nannastacidae embedded within it. Comparison of the molecular phylogeny with that based on morphology suggests that many "diagnostic" characters are homoplasious.

Insect allatotropins belong to a family of structurally-related myoactive peptides present in several invertebrate phyla

Originally named for its ability to stimulate juvenile hormone production by lepidopteran corpora allata, allatotropin has emerged as a neuropeptide with multiple neural, endocrine and myoactive roles. This paper describes the experimental evidence for allatotropin action, its localization in several species of insects, and its multiple effects on a variety of different tissues that lead to increased hemolymph circulation and gut motility. The overall physiological effects may also include species-specific effects such as the regulation of nutrient absorption, modulation of the circadian cycle and migratory preparedness. In addition, we present evidence suggesting that allatotropins are members of a family of myoactive peptides found in several invertebrate phyla. Finally, we speculate that the myoactive properties of allatotropins are basal and it is likely that the stimulatory action of allatotropins on juvenile hormone synthesis evolved secondarily.

Functional specificity of the nematode Hox gene mab-5

Hox genes encode evolutionarily conserved transcription factors involved in morphological specification along the anteroposterior body axis of animals. The two most striking features of Hox genes are colinearity and the strong sequence conservation. Among all animals studied so far, the nematode Caenorhabditis elegans contains one of the most divergent Hox clusters. The core cluster contains only four members, which in part deviate from the colinearity rule. In addition, orthologous and paralogous nematode Hox sequences diverged substantially. Given these nematode-specific features, we asked how these Hox proteins evolved and how they provide functional specificity. We investigated the role of MAB-5 during ray formation and established an in vivo assay using Cel-mab-5 regulatory elements to express orthologous, paralogous and chimeric cDNAs in a Cel-mab-5 mutant background. We show that the MAB-5 ortholog from Pristionchus pacificus, but not the C. elegans paralogous Hox proteins can rescue Cel-mab-5. Experiments with chimeric, truncated and mutagenized Hox proteins suggest the specificity to be conferred by the N-terminal arm and helix I, but not helix II of the homeodomain.

Hox genes and the evolution of the arthropod body plan

In recent years researchers have analyzed the expression patterns of the Hox genes in a multitude of arthropod species, with the hope of understanding the mechanisms at work in the evolution of the arthropod body plan. Now, with Hox expression data representing all four major groups of arthropods (chelicerates, myriapods, crustaceans, and insects), it seems appropriate to summarize the results and take stock of what has been learned. In this review we summarize the expression and functional data regarding the 10 arthropod Hox genes: labial proboscipedia, Hox3/zen, Deformed, Sex combs reduced, fushi tarazu, Antennapedia, Ultrabithorax, abdominal-A, and Abdominal-B. In addition, we discuss mechanisms of developmental evolutionary change thought to be important for the emergence of novel morphological features within the arthropods.

Specificity of Distalless repression and limb primordia development by abdominal Hox proteins

In Drosophila, differences between segments, such as the presence or absence of appendages, are controlled by Hox transcription factors. The Hox protein Ultrabithorax (Ubx) suppresses limb formation in the abdomen by repressing the leg selector gene Distalless, whereas Antennapedia (Antp), a thoracic Hox protein, does not repress Distalless. We show that the Hox cofactors Extradenticle and Homothorax selectively enhance Ubx, but not Antp, binding to a Distalless regulatory sequence. A C-terminal peptide in Ubx stimulates binding to this site. However, DNA binding is not sufficient for Distalless repression. Instead, an additional alternatively spliced domain in Ubx is required for Distalless repression but not DNA binding. Thus, the functional specificities of Hox proteins depend on both DNA binding-dependent and -independent mechanisms.

Drosophila atonal fully rescues the phenotype of Math1 null mice: new functions evolve in new cellular contexts

Many genes share sequence similarity between species, but their properties often change significantly during evolution. For example, the Drosophila genes engrailed and orthodenticle and the onychophoran gene Ultrabithorax only partially substitute for their mouse or Drosophila homologs. We have been analyzing the relationship between atonal (ato) in the fruit fly and its mouse homolog, Math1. In flies, ato acts as a proneural gene that governs the development of chordotonal organs (CHOs), which serve as stretch receptors in the body wall and joints and as auditory organs in the antennae. In the fly CNS, ato is important not for specification but for axonal arborization. Math1, in contrast, is required for the specification of cells in both the CNS and the PNS. Furthermore, Math1 serves a role in the development of secretory lineage cells in the gut, a function that does not parallel any known to be served by ato. We wondered whether ato and Math1 might be more functionally homologous than they appear, so we expressed Math1 in ato mutant flies and ato in Math1 null mice. To our surprise, the two proteins are functionally interchangeable.

Functional evolution of Hox proteins in arthropods

It is presumed that the evolution of morphological diversity in animals and plants is driven by changes in the developmental processes that govern morphology, hence basically by changes in the function and/or expression of a defined set of genes that control these processes. A large body of evidence has suggested that changes in developmental gene regulation are the predominant mechanisms that sustain morphological evolution, being much more important than the evolution of the primary sequences and functions of proteins. Recent reports challenge this idea by highlighting functional evolution of Hox proteins during the evolutionary history of arthropods.

Arthropod Hox genes: insights on the evolutionary forces that shape gene functions

Comparative studies suggest that gene duplication, changes in cis-regulatory elements and changes in protein sequence all contribute to the evolution of Hox gene functions, but the evolutionary dynamics of these changes are probably different. It seems likely that gene duplications arise as neutral changes and acquire an adaptive significance later on. By contrast, some changes in regulatory and protein-coding sequences can have immediate consequences in morphological evolution.