Comparison of bicoid-dependent regulation of hunchback between Musca domestica and Drosophila melanogaster

Modularity, criticality, and evolvability of a developmental gene regulatory network

The existence of discrete phenotypic traits suggests that the complex regulatory processes which produce them are functionally modular. These processes are usually represented by networks. Only modular networks can be partitioned into intelligible subcircuits able to evolve relatively independently. Traditionally, functional modularity is approximated by detection of modularity in network structure. However, the correlation between structure and function is loose. Many regulatory networks exhibit modular behaviour without structural modularity. Here we partition an experimentally tractable regulatory network-the gap gene system of dipteran insects-using an alternative approach. We show that this system, although not structurally modular, is composed of dynamical modules driving different aspects of whole-network behaviour. All these subcircuits share the same regulatory structure, but differ in components and sensitivity to regulatory interactions. Some subcircuits are in a state of criticality, while others are not, which explains the observed differential evolvability of the various expression features in the system.

The Functionality and Evolution of Eukaryotic Transcriptional Enhancers

Enhancers regulate precise spatial and temporal patterns of gene expression in eukaryotes and, moreover, evolutionary changes in these modular cis-regulatory elements may represent the predominant genetic basis for phenotypic evolution. Here, we review approaches to identify and functionally analyze enhancers and their transcription factor binding sites, including assay for transposable-accessible chromatin-sequencing (ATAC-Seq) and clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9, respectively. We also explore enhancer functionality, including how transcription factor binding sites combine to regulate transcription, as well as research on shadow and super enhancers, and how enhancers can act over great distances and even in trans. Finally, we discuss recent theoretical and empirical data on how transcription factor binding sites and enhancers evolve. This includes how the function of enhancers is maintained despite the turnover of transcription factor binding sites as well as reviewing studies where mutations in enhancers have been shown to underlie morphological change.

Emerging developmental genetic model systems in holometabolous insects

The number of insect species that are amenable to functional genetic studies is growing rapidly and provides many new research opportunities in developmental and evolutionary biology. The holometabolous insects represent a disproportionate percentage of animal diversity and are thus well positioned to provide model species for a wide variety of developmental processes. Here we discuss emerging holometabolous models, and review some recent breakthroughs. For example, flies and midges were found to use structurally unrelated long-range pattern organizers, butterflies and moths revealed extensive pattern formation during oogenesis, new imaging possibilities in the flour beetle Tribolium castaneum showed how embryos break free of their extraembryonic membranes, and the complex genetics governing interspecies difference in head shape were revealed in Nasonia wasps.

Sequential construction of a model for modular gene expression control, applied to spatial patterning of theDrosophilagenehunchback

Gene network simulations are increasingly used to quantify mutual gene regulation in biological tissues. These are generally based on linear interactions between single-entity regulatory and target genes. Biological genes, by contrast, commonly have multiple, partially independent, cis-regulatory modules (CRMs) for regulator binding, and can produce variant transcription and translation products. We present a modeling framework to address some of the gene regulatory dynamics implied by this biological complexity. Spatial patterning of the hunchback (hb) gene in Drosophila development involves control by three CRMs producing two distinct mRNA transcripts. We use this example to develop a differential equations model for transcription which takes into account the cis-regulatory architecture of the gene. Potential regulatory interactions are screened by a genetic algorithms (GAs) approach and compared to biological expression data.

The achaete-scute complex in Diptera: patterns of noncoding sequence evolution

The achaete‐scute complex (AS‐C) has been a useful paradigm for the study of pattern formation and its evolution. achaete‐scute genes have duplicated and evolved distinct expression patterns during the evolution of cyclorraphous Diptera. Are the expression patterns in different species driven by conserved regulatory elements? If so, when did such regulatory elements arise? Here, we have sequenced most of the AS‐C of the fly Calliphora vicina (including the genes achaete, scute and lethal of scute) to compare noncoding sequences with known cis‐regulatory sequences in Drosophila. The organization of the complex is conserved with respect to Drosophila species. There are numerous small stretches of conserved noncoding sequence that, in spite of high sequence turnover, display binding sites for known transcription factors. Synteny of the blocks of conserved noncoding sequences is maintained suggesting not only conservation of the position of regulatory elements but also an origin prior to the divergence between these two species. We propose that some of these enhancers originated by duplication with their target genes.

Mid-embryo patterning and precision in Drosophila segmentation: Krüppel dual regulation of hunchback

In early development, genes are expressed in spatial patterns which later define cellular identities and tissue locations. The mechanisms of such pattern formation have been studied extensively in early Drosophila (fruit fly) embryos. The gap gene hunchback (hb) is one of the earliest genes to be expressed in anterior-posterior (AP) body segmentation. As a transcriptional regulator for a number of downstream genes, the spatial precision of hb expression can have significant effects in the development of the body plan. To investigate the factors contributing to hb precision, we used fine spatial and temporal resolution data to develop a quantitative model for the regulation of hb expression in the mid-embryo. In particular, modelling hb pattern refinement in mid nuclear cleavage cycle 14 (NC14) reveals some of the regulatory contributions of simultaneously-expressed gap genes. Matching the model to recent data from wild-type (WT) embryos and mutants of the gap gene Krüppel (Kr) indicates that a mid-embryo Hb concentration peak important in thoracic development (at parasegment 4, PS4) is regulated in a dual manner by Kr, with low Kr concentration activating hb and high Kr concentration repressing hb. The processes of gene expression (transcription, translation, transport) are intrinsically random. We used stochastic simulations to characterize the noise generated in hb expression. We find that Kr regulation can limit the positional variability of the Hb mid-embryo border. This has been recently corroborated in experimental comparisons of WT and Kr- mutant embryos. Further, Kr regulation can decrease uncertainty in mid-embryo hb expression (i.e. contribute to a smooth Hb boundary) and decrease between-copy transcriptional variability within nuclei. Since many tissue boundaries are first established by interactions between neighbouring gene expression domains, these properties of Hb-Kr dynamics to diminish the effects of intrinsic expression noise may represent a general mechanism contributing to robustness in early development.

Maternal co-ordinate gene regulation and axis polarity in the scuttle fly Megaselia abdita

Axis specification and segment determination in dipteran insects are an excellent model system for comparative analyses of gene network evolution. Antero-posterior polarity of the embryo is established through systems of maternal morphogen gradients. In Drosophila melanogaster, the anterior system acts through opposing gradients of Bicoid (Bcd) and Caudal (Cad), while the posterior system involves Nanos (Nos) and Hunchback (Hb) protein. These systems act redundantly. Both Bcd and Hb need to be eliminated to cause a complete loss of polarity resulting in mirror-duplicated abdomens, so-called bicaudal phenotypes. In contrast, knock-down of bcd alone is sufficient to induce double abdomens in non-drosophilid cyclorrhaphan dipterans such as the hoverfly Episyrphus balteatus or the scuttle fly Megaselia abdita. We investigate conserved and divergent aspects of axis specification in the cyclorrhaphan lineage through a detailed study of the establishment and regulatory effect of maternal gradients in M. abdita. Our results show that the function of the anterior maternal system is highly conserved in this species, despite the loss of maternal cad expression. In contrast, hb does not activate gap genes in this species. The absence of this activatory role provides a precise genetic explanation for the loss of polarity upon bcd knock-down in M. abdita, and suggests a general scenario in which the posterior maternal system is increasingly replaced by the anterior one during the evolution of the cyclorrhaphan dipteran lineage.

The basic head-to-tail polarity of an animal is established very early in development. In dipteran insects (flies, midges, and mosquitoes), polarity is established with the help of so-called morphogen gradients. Morphogens are regulatory proteins that are distributed as a concentration gradient, often involving diffusion from a localised source. This graded distribution then leads to the concentration-dependent activation of different target genes along the embryo’s axis. We examine this process, which differs to a surprising extent between dipteran species, in the scuttle fly Megaselia abdita, and compare our results to the model organism Drosophila melanogaster. In this way, we not only gain insights into how the mechanisms that establish polarity function differently in different species, but also how the system has evolved since these two flies shared a common ancestor. Specifically, we pin down the main difference between Drosophila and Megaselia in the altered function of the maternal Hunchback morphogen gradient, which activates target genes in the former, but not the latter species, where it has been completely replaced by the Bicoid morphogen during evolution.

Quantitative system drift compensates for altered maternal inputs to the gap gene network of the scuttle fly Megaselia abdita

The segmentation gene network in insects can produce equivalent phenotypic outputs despite differences in upstream regulatory inputs between species. We investigate the mechanistic basis of this phenomenon through a systems-level analysis of the gap gene network in the scuttle fly Megaselia abdita (Phoridae). It combines quantification of gene expression at high spatio-temporal resolution with systematic knock-downs by RNA interference (RNAi). Initiation and dynamics of gap gene expression differ markedly between M. abdita and Drosophila melanogaster, while the output of the system converges to equivalent patterns at the end of the blastoderm stage. Although the qualitative structure of the gap gene network is conserved, there are differences in the strength of regulatory interactions between species. We term such network rewiring 'quantitative system drift'. It provides a mechanistic explanation for the developmental hourglass model in the dipteran lineage. Quantitative system drift is likely to be a widespread mechanism for developmental evolution.

Comparative transcriptomics of early dipteran development

Background

Modern sequencing technologies have massively increased the amount of data available for comparative genomics. Whole-transcriptome shotgun sequencing (RNA-seq) provides a powerful basis for comparative studies. In particular, this approach holds great promise for emerging model species in fields such as evolutionary developmental biology (evo-devo).

Results

We have sequenced early embryonic transcriptomes of two non-drosophilid dipteran species: the moth midge Clogmia albipunctata, and the scuttle fly Megaselia abdita. Our analysis includes a third, published, transcriptome for the hoverfly Episyrphus balteatus. These emerging models for comparative developmental studies close an important phylogenetic gap between Drosophila melanogaster and other insect model systems. In this paper, we provide a comparative analysis of early embryonic transcriptomes across species, and use our data for a phylogenomic re-evaluation of dipteran phylogenetic relationships.

Conclusions

We show how comparative transcriptomics can be used to create useful resources for evo-devo, and to investigate phylogenetic relationships. Our results demonstrate that de novo assembly of short (Illumina) reads yields high-quality, high-coverage transcriptomic data sets. We use these data to investigate deep dipteran phylogenetic relationships. Our results, based on a concatenation of 160 orthologous genes, provide support for the traditional view of Clogmia being the sister group of Brachycera (Megaselia, Episyrphus, Drosophila), rather than that of Culicomorpha (which includes mosquitoes and blackflies).

Coevolution within and between regulatory loci can preserve promoter function despite evolutionary rate acceleration

Phenotypes that appear to be conserved could be maintained not only by strong purifying selection on the underlying genetic systems, but also by stabilizing selection acting via compensatory mutations with balanced effects. Such coevolution has been invoked to explain experimental results, but has rarely been the focus of study. Conserved expression driven by the unc-47 promoters of Caenorhabditis elegans and C. briggsae persists despite divergence within a cis-regulatory element and between this element and the trans-regulatory environment. Compensatory changes in cis and trans are revealed when these promoters are used to drive expression in the other species. Functional changes in the C. briggsae promoter, which has experienced accelerated sequence evolution, did not lead to alteration of gene expression in its endogenous environment. Coevolution among promoter elements suggests that complex epistatic interactions within cis-regulatory elements may facilitate their divergence. Our results offer a detailed picture of regulatory evolution in which subtle, lineage-specific, and compensatory modifications of interacting cis and trans regulators together maintain conserved gene expression patterns.

Some phenotypes, including gene expression patterns, are conserved between distantly related species. However, the molecular bases of those phenotypes are not necessarily conserved. Instead, regulatory DNA sequences and the proteins with which they interact can change over time with balanced effects, preserving expression patterns and concealing regulatory divergence. Coevolution between interacting molecules makes gene regulation highly species-specific, and it can be detected when the cis-regulatory DNA of one species is used to drive expression in another species. In this way, we identified regions of the C. elegans and C. briggsae unc-47 promoters that have coevolved with the lineage-specific trans-regulatory environments of these organisms. The C. briggsae promoter experienced accelerated sequence change relative to related species. All of this evolution occurred without changing the expression pattern driven by the promoter in its endogenous environment.

The gap gene network

Gap genes are involved in segment determination during the early development of the fruit fly Drosophila melanogaster as well as in other insects. This review attempts to synthesize the current knowledge of the gap gene network through a comprehensive survey of the experimental literature. I focus on genetic and molecular evidence, which provides us with an almost-complete picture of the regulatory interactions responsible for trunk gap gene expression. I discuss the regulatory mechanisms involved, and highlight the remaining ambiguities and gaps in the evidence. This is followed by a brief discussion of molecular regulatory mechanisms for transcriptional regulation, as well as precision and size-regulation provided by the system. Finally, I discuss evidence on the evolution of gap gene expression from species other than Drosophila. My survey concludes that studies of the gap gene system continue to reveal interesting and important new insights into the role of gene regulatory networks in development and evolution.

Co-evolution of transcriptional silencing proteins and the DNA elements specifying their assembly

Co-evolution of transcriptional regulatory proteins and their sites of action has been often hypothesized but rarely demonstrated. Here we provide experimental evidence of such co-evolution in yeast silent chromatin, a finding that emerged from studies of hybrids formed between two closely related Saccharomyces species. A unidirectional silencing incompatibility between S. cerevisiae and S. bayanus led to a key discovery: asymmetrical complementation of divergent orthologs of the silent chromatin component Sir4. In S. cerevisiae/S. bayanus interspecies hybrids, ChIP-Seq analysis revealed a restriction against S. cerevisiae Sir4 associating with most S. bayanus silenced regions; in contrast, S. bayanus Sir4 associated with S. cerevisiae silenced loci to an even greater degree than did S. cerevisiae's own Sir4. Functional changes in silencer sequences paralleled changes in Sir4 sequence and a reduction in Sir1 family members in S. cerevisiae. Critically, species-specific silencing of the S. bayanus HMR locus could be reconstituted in S. cerevisiae by co-transfer of the S. bayanus Sir4 and Kos3 (the ancestral relative of Sir1) proteins. As Sir1/Kos3 and Sir4 bind conserved silencer-binding proteins, but not specific DNA sequences, these rapidly evolving proteins served to interpret differences in the two species' silencers presumably involving emergent features created by the regulatory proteins that bind sequences within silencers. The results presented here, and in particular the high resolution ChIP-Seq localization of the Sir4 protein, provided unanticipated insights into the mechanism of silent chromatin assembly in yeast.

Anterior development in the parthenogenetic and viviparous form of the pea aphid, Acyrthosiphon pisum: hunchback and orthodenticle expression

In the dipteran Drosophila, the genes bicoid and hunchback work synergistically to pattern the anterior blastoderm during embryogenesis. bicoid, however, appears to be an innovation of the higher Diptera. Hence, in some non-dipteran insects, anterior specification instead relies on a synergistic interaction between maternally transcribed hunchback and orthodenticle. Here we describe how orthologues of hunchback and orthodenticle are expressed during oogenesis and embryogenesis in the parthenogenetic and viviparous form of the pea aphid, Acyrthosiphon pisum. A. pisum hunchback (Aphb) mRNA is localized to the anterior pole in developing oocytes and early embryos prior to blastoderm formation - a pattern strongly reminiscent of bicoid localization in Drosophila. A. pisum orthodenticle (Apotd), on the other hand, is not expressed prior to gastrulation, suggesting that it is the asymmetric localization of Aphb, rather than synergy between Aphb and Apotd, that regulates anterior specification in asexual pea aphids.

Evolution of cis-regulatory sequence and function in Diptera

Cis-regulatory sequences direct patterns of gene expression essential for development and physiology. Evolutionary changes in these sequences contribute to phenotypic divergence. Despite their importance, cis-regulatory regions remain one of the most enigmatic features of the genome. Patterns of sequence evolution can be used to identify cis-regulatory elements, but the power of this approach depends upon the relationship between sequence and function. Comparative studies of gene regulation among Diptera reveal that divergent sequences can underlie conserved expression, and that expression differences can evolve despite largely similar sequences. This complex structure-function relationship is the primary impediment for computational identification and interpretation of cis-regulatory sequences. Biochemical characterization and in vivo assays of cis-regulatory sequences on a genomic-scale will relieve this barrier.

Evolutionary and functional analysis of the tailless enhancer in Musca domestica and Drosophila melanogaster

To further understand the evolutionary dynamics of the regulatory interactions underlying development, we expand on our previous analysis of hunchback and compare the structure and function of the tailless enhancer between Musca domestica and Drosophila melanogaster. Our analysis shows that although the expression patterns and functional protein domains of tll are conserved between Musca and Drosophila, the enhancer sequences are unalignable. Upon closer investigation, we find that these highly diverged enhancer sequences encode the same regulatory information necessary for Bicoid, Dorsal, and the terminal system to drive tll expression. The binding sites for these transcription factors differ in the sequence, number, spacing, and position between the Drosophila and Musca tll enhancers, and we were unable to establish homology between binding sites from each species. This implies that the Musca and Drosophila Bcd-binding sites have evolved de novo in the 100 million years since these species diverged. However, in transgenic Drosophila embryos the Musca tll enhancer is able to drive the same expression pattern as endogenous Drosophila tll. Therefore, during the rapid evolution of enhancer sequences individual binding sites are continually lost and gained, but the transcriptional output is maintained by compensatory mutations in cis and in trans.

The origins of eukaryotic gene structure

Most of the phenotypic diversity that we perceive in the natural world is directly attributable to the peculiar structure of the eukaryotic gene, which harbors numerous embellishments relative to the situation in prokaryotes. The most profound changes include introns that must be spliced out of precursor mRNAs, transcribed but untranslated leader and trailer sequences (untranslated regions), modular regulatory elements that drive patterns of gene expression, and expansive intergenic regions that harbor additional diffuse control mechanisms. Explaining the origins of these features is difficult because they each impose an intrinsic disadvantage by increasing the genic mutation rate to defective alleles. To address these issues, a general hypothesis for the emergence of eukaryotic gene structure is provided here. Extensive information on absolute population sizes, recombination rates, and mutation rates strongly supports the view that eukaryotes have reduced genetic effective population sizes relative to prokaryotes, with especially extreme reductions being the rule in multicellular lineages. The resultant increase in the power of random genetic drift appears to be sufficient to overwhelm the weak mutational disadvantages associated with most novel aspects of the eukaryotic gene, supporting the idea that most such changes are simple outcomes of semi-neutral processes rather than direct products of natural selection. However, by establishing an essentially permanent change in the population-genetic environment permissive to the genome-wide repatterning of gene structure, the eukaryotic condition also promoted a reliable resource from which natural selection could secondarily build novel forms of organismal complexity. Under this hypothesis, arguments based on molecular, cellular, and/or physiological constraints are insufficient to explain the disparities in gene, genomic, and phenotypic complexity between prokaryotes and eukaryotes.

The origin of subfunctions and modular gene regulation

Evolutionary explanations for the origin of modularity in genetic and developmental pathways generally assume that modularity confers a selective advantage. However, our results suggest that even in the absence of any direct selective advantage, genotypic modularity may increase through the formation of new subfunctions under near-neutral processes. Two subfunctions may be formed from a single ancestral subfunction by the process of fission. Subfunction fission occurs when multiple functions under unified genetic control become subdivided into more restricted functions under independent genetic control. Provided that population size is sufficiently small, random genetic drift and mutation can conspire to produce changes in the number of subfunctions in the genome of a species without necessarily altering the phenotype. Extensive genotypic modularity may then accrue in a near-neutral fashion in permissive population-genetic environments, potentially opening novel pathways to morphological evolution. Many aspects of gene complexity in multicellular eukaryotes may have arisen passively as population size reductions accompanied increases in organism size, with the adaptive exploitation of such complexity occurring secondarily.

How to get ahead: the origin, evolution and function ofbicoid

In Drosophila, a Bcd protein gradient orchestrates patterning along the anteroposterior embryonic axis. However, studies of basal flies and other insects have revealed that bcd is a derived Hox3 gene found only in higher dipterans. To understand how bcd acquired its role in flies and how anteroposterior patterning mechanisms have evolved, I first review key features of bcd function in Drosophila: anterior localization and transcriptional and translation control of gene expression. I then discuss investigations of bcd in other higher dipterans that have provided insight into the evolution of regulatory interactions and the Bcd gradient. Finally, I review studies of Drosophila and other insects that address the evolution of bcd function and integration of bcd into ancestral regulatory mechanisms. I suggest further comparative studies may allow us to identify the intermediate steps in bcd evolution. This will make bcd a paradigm for the origin and evolution of genes and regulatory networks.

[Extreme divergence of a homeotic gene: the bicoid case]

Evolutionary developmental genetics (evo-devo) reveals that the plasticity of development is so important that every developmental biology project should carefully take this point into consideration. The example of bicoid, the first discovered morphogen, illustrates how an essential gene can change its function during evolution. The search for bicoid homologues showed that this gene is surprisingly specific to flies (cyclorraphan diptera) and absent in other insects. In fact, recent studies demonstrate that bicoid is a very derived Hox3 homeotic gene. During insect evolution, the ancestral Hox3 gene lost its homeotic function and acquired new roles in oocytes and embryonic annexes. Then, in the lineage leading to modern flies, a duplication of this new gene, followed by functional divergence, led to the formation of bicoid and zerknüllt. Both genes are located within the Drosophila Hox complex; however, they have no homeotic function. Thanks to the power of Drosophila genetics, it is possible to suggest that torso and hunchback may constitute the insect primitive anterior organizer. The bicoid evolutionary history reveals several fundamental mechanisms of the evolution of developmental genes, such as changes of gene regulation, modifications of protein sequences and gene duplication. It also shows the need for studying a wider range of model organisms before generalisations can be made from data obtained with one particular species.

Evolution of development in closely related species of flies and worms

One of the main challenges in evolutionary biology is to identify the molecular changes that underlie phenotypic differences that are of evolutionary significance. Comparative studies of early development have shown that changes in the spatio-temporal use of regulatory genes, as well as changes in the specificity of regulatory proteins, are correlated with important differences in morphology between phylogenetically distant species. However, it is not known how such changes take place in natural populations, and whether they result from a single, or many small, additive events. Extending this approach to the study of development of closely related species promises to enrich this debate.

Coevolution in bicoid-dependent promoters and the inception of regulatory incompatibilities among species of higher Diptera

To what extent and in what way do gene promoters and their transacting regulatory proteins coevolve? In this and in earlier publications we show that the Bicoid-dependent promoters of the segmentation genes hunchback and tailless in species of higher Diptera (Drosophila, Musca, Calliphora, and Lucilia) are different with respect to the copy number, spacing, sequence, and orientation of Bicoid binding sites. At the same time there are significant amino acid differences in the Bicoid homeodomain. To test these interspecific differences, we used a series of functional assays, starting with the analysis of Bicoid binding affinities of individual sites, through to transgene rescue experiments, to compare within-species with between-species mixtures of Bicoid homeodomains and hunchback or tailless promoters. We observed that components taken from different species interact with less efficiency compared with those taken from within the same species. Our interpretation is that such interspecific incompatibilities are a consequence of interactive genetic elements coevolving one with another, hence maintaining functional compatibility within each species. At the same time such a process allows differences to accumulate between species regarding the precise molecular basis whereby the common function is effected.

Rapid restructuring of bicoid-dependent hunchback promoters within and between Dipteran species: implications for molecular coevolution

Interacting genetic elements need to coevolve if their joint function is to be maintained; for example, the correct binding of transcriptional regulators to defined binding sites in gene promoters needs to be maintained during evolution to ensure proper function. As part of a wider investigation into the molecular coevolution of the Dipteran homeodomain-bearing regulator bicoid (bcd) and Bcd-dependent promoters, we present data on the functional, structural, and sequence differences between the promoters of the segmentation gene hunchback (hb), in several species of Cyclorrhaphan (higher) Diptera. The result of phenocopying hb mutations using RNA interference (RNAi) in Musca domestica shows broadly similar functions to the hb gene in Drosophila melanogaster. However, the Bcd-binding sites in the hb promoters of Drosophila, Musca, and the two blowfly species Lucilia sericata and Calliphora vicina differ in copy number, sequence, orientation, and spacing. Furthermore, all promoters are subject to rapid turnover by slippage-like processes leading to high densities of short repetitive motifs. A study of polymorphism among six strains of M. domestica reveals that turnover by slippage also occurs in the promoter, untranslated leader, and exonic coding sequences of hb, but to different extents. We discuss these results in terms of the known interspecific differences in bcdand the potential coevolution of selected compensatory mutations in trans and cis in response to continuous promoter restructuring.

Divergent structure and function of the bicoid gene in Muscoidea fly species

We have investigated the evolution of the bicoid (bcd) gene in fly species of the Muscoidea Superfamily. We obtained the complete bcd sequence from the housefly Musca domestica and found polymorphism in the coding region among Musca strains. In addition to Musca, we cloned most of the bcd coding sequences from two blowfly species Calliphora vicina and Lucilia sericata. The 5' and 3' regulatory regions flanking the Musca bcd gene are widely diverged in sequence from Drosophila; however, some important sequence motifs identified in Drosophila bcd are present. The predicted RNA secondary structures of the 3' UTRs are similar, despite sequence divergence. Comparison of Bicoid (Bcd) proteins shows a serine-rich domain of unknown function is present in the Muscoidea species, but is absent in other species. The in vivo function of bcd in Musca was tested by RNAi to mimic loss of function phenotype. We obtained a head defect phenotype similar to weak bcd alleles of Drosophila. Although our comparisons initially suggest functional conservation between species, closer inspection reveals significant differences. Divergence of structural motifs, such as regulatory elements in flanking regions and conservation of protein domains in some species but not in others, points to functional divergence between species. We suggest that the larger embryonic size in Muscoidea species restricts the morphogenetic activity of a weak Bcd activator, which has evolved a more specialized role in head determination and lost some functions in thoracic development.

How genomic and developmental dynamics affect evolutionary processes

Evolutionary genetics is concerned with natural selection and neutral drift, to the virtual exclusion of almost everything else. In its current focus on DNA variation, it reduces phenotypes to symbols. Varying phenotypes, however, are the units of evolution, and, if we want a comprehensive theory of evolution, we need to consider both the internal and external evolutionary forces that shape the development of phenotypes. Genetic systems are redundant, modular and subject to a variety of genomic mechanisms of "turnover" (transposition, gene conversion, unequal crossingover, slippage and so on). As such the construction and spread of novel combinations of modules by turnover, in particular within gene promoters, contributes significantly to the evolution of phenotypes. Furthermore, redundancy, turnover and modularity lead to ever more complex networks of genetic interactions and ever more functions for a given module. The significant interaction between genomic turnover and natural selection leads to a molecular coevolution between interacting modules and hence facilitates the establishment of biological novelties.

Target selectivity of bicoid is dependent on nonconsensus site recognition and protein-protein interaction

We describe experiments to compare the activities of two Drosophila homeodomain proteins, Bicoid (Bcd) and an altered-specificity mutant of Fushi tarazu, Ftz(Q50K). Although the homeodomains of these proteins share a virtually indistinguishable ability to recognize a consensus Bcd site, only Bcd can activate transcription from natural enhancer elements when assayed in both yeast and Drosophila Schneider S2 cells. Our analysis of chimeric proteins suggests that both the homeodomain of Bcd and sequences outside the homeodomain contribute to its ability to recognize natural enhancer elements. We further show that, unlike the Bcd homeodomain, the Ftz(Q50K) homeodomain fails to recognize nonconsensus sites found in natural enhancer elements. The defect of a chimeric protein containing the homeodomain of Ftz(Q50K) in place of that of Bcd can be preferentially restored by converting the nonconsensus sites in natural enhancer elements to consensus sites. Our experiments suggest that the biological specificity of Bcd is determined by combinatorial contributions of two important mechanisms: the nonconsensus site recognition function conferred by the homeodomain and the cooperativity function conferred primarily by sequences outside the homeodomain. A systematic comparison of different assay methods and enhancer elements further suggests a fluid nature of the requirements for these two Bcd functions in target selection.

Function of bicoid and hunchback homologs in the basal cyclorrhaphan fly Megaselia (Phoridae)

The Drosophila gene bicoid functions at the beginning of a gene cascade that specifies anterior structures in the embryo. Its transcripts are localized at the anterior pole of the oocyte, giving rise to a Bicoid protein gradient, which regulates the spatially restricted expression of target genes along the anterior-posterior axis of the embryo in a concentration-dependent manner. The morphogen function of Bicoid requires the coactivity of the zinc finger transcription factor Hunchback, which is expressed in a Bicoid-dependent fashion in the anterior half of the embryo. Whereas hunchback is conserved throughout insects, bicoid homologs are known only from cyclorrhaphan flies. Thus far, identification of hunchback and bicoid homologs rests only on sequence comparison. In this study, we used double-stranded RNA interference (RNAi) to address the function of bicoid and hunchback homologs in embryos of the lower cyclorrhaphan fly Megaselia abdita (Phoridae). Megaselia-hunchback RNAi causes hunchback-like phenotypes as observed in Drosophila, but Megaselia-bicoid RNAi causes phenotypes different from corresponding RNAi experiments in Drosophila and bicoid mutant embryos. Megaselia-bicoid is required not only for the head and thorax but also for the development of four abdominal segments. This difference between Megaselia and Drosophila suggests that the range of functional bicoid activity has been reduced in higher flies.

Genetic analysis of bristle loss in hybrids between Drosophila melanogaster and D. simulans provides evidence for divergence of cis-regulatory sequences in the achaete-scute gene complex

The two closely related species of Drosophila, D. melanogaster and D. simulans, display an identical bristle pattern on the notum, but hybrids between the two are lacking a variable number of bristles. We show that the loss is temperature-dependent and provide evidence for two periods of temperature sensitivity. A first period of heat sensitivity occurs during larval development and corresponds to the time when the prepattern of expression of genes whose products activate achaete-scute in the proneural clusters preceding bristle precursor formation is established. A second period of cold sensitivity corresponds to the time of emergence of the bristle precursor cells and the maintenance of their neural fate, a process requiring high levels of Achaete-Scute. Expression of achaete-scute at these two critical periods depends on cis-regulatory elements of the achaete-scute complex (AS-C). The differences between males, which have only one copy of the X-linked AS-C from D. simulans, and females, which have copies from both parental species, are compared, together with the effects of crossing in different rearrangements of the D. melanogaster AS-C that delete regulatory and/or coding sequences. We provide evidence that bristle loss in the hybrids may result from a decrease in the level of transcription at the AS-C and argue that interaction between trans-acting factors and cis-regulatory elements within the AS-C has diverged between the two species.

Regulation of the Tribolium homologues of caudal and hunchback in Drosophila: evidence for maternal gradient systems in a short germ embryo

In short germ embryos, the germ rudiment forms at the posterior ventral side of the egg, while the anterior-dorsal region becomes the extra-embryonic serosa. It is difficult to see how an anterior gradient like that of bicoid in Drosophila could in these embryos be directly involved in patterning of the germ rudiment. Moreover, since it has not yet been possible to recover a bicoid homologue from any species outside the diptera, it has been speculated that the anterior bicoid gradient could be a late addition during insect evolution. We addressed this question by analysing the regulation of potential target genes of bicoid in the short germ embryo of Tribolium castaneum. We demonstrate that homologues of caudal and hunchback from Tribolium are regulated by Drosophila bicoid. In Drosophila, maternal caudal RNA is translationally repressed by bicoid. We find that Tribolium caudal RNA is also translationally repressed by bicoid, when it is transferred into Drosophila embryos under a maternal promoter. This strongly suggests that a functional bicoid homologue must exist in Tribolium. The second target gene, hunchback, is transcriptionally activated by bicoid in Drosophila. Transfer of the regulatory region of Tribolium hunchback into Drosophila also results in regulation by early maternal factors, including bicoid, but in a pattern that is more reminiscent of Tribolium hunchback expression, namely in two early blastoderm domains. Using enhancer mapping constructs and footprinting, we show that caudal activates the posterior of these domains via a specific promoter. Our experiments suggest that a major event in the evolutionary transition from short to long germ embryogenesis was the switch from activation of the hunchback gap domain by caudal to direct activation by bicoid. This regulatory switch can explain how this domain shifted from a posterior location in short germ embryos to its anterior position in long germ insects, and it also suggest how an anterior gradient can pattern the germ rudiment in short germ embryos, i.e. by regulating the expression of caudal.