mRNA diffusion explains protein gradients in Drosophila early development

We propose a new model describing the production and the establishment of the stable gradient of the Bicoid protein along the antero-posterior axis of the embryo of Drosophila. In this model, we consider that bicoid mRNA diffuses along the antero-posterior axis of the embryo and the protein is produced in the ribosomes localized near the syncytial nuclei. Bicoid protein stays localized near the syncytial nuclei as observed in experiments. We calibrate the parameters of the mathematical model with experimental data taken during the cleavage stages 11-14 of the developing embryo of Drosophila. We obtain good agreement between the experimental and the model gradients, with relative errors in the range 5-8%. The inferred diffusion coefficient of bicoid mRNA is in the range , in agreement with the theoretical predictions and experimental measurements for the diffusion of macromolecules in the cytoplasm. We show that the model based on the mRNA diffusion hypothesis is consistent with the known observational data, supporting the recent experimental findings of the gradient of bicoid mRNA in Drosophila [Spirov et al. (2009). Development 136, 605-614].     

A comparison of transplantable bicoid activity and partial bicoid homeobox sequences in several Drosophila and blowfly species (Calliphoridae)

In order to test for bicoid-like activity in insects other than Drosophila melanogaster, anterior egg cytoplasm from the following species was injected into cleavage stage embryos from mutant D. melanogaster lacking a functional bicoid (bcd) product: six other Drosophila species, the housefly, three blowfly species, the primitive cyclorrhaphic dipteran Megaselia, and the honeybee Apis mellifera; preliminary tests were made with four lower dipterans (Nematocera). Rescue effects were only observed with the drosophilids, housefly, and two of the three blowfly species. Rescue was stronger with the drosophilids than with the other flies as donors. Where checked (D. pseudoobscura), a positive correlation was found between the amount of cytoplasm injected and the number of pattern elements formed, suggesting threshold effects upon target genes as with the endogenous bcd product. By polymerase chain reaction, fragments from a bcd-orthologous homeobox were cloned from the three blowfly species. The derived sequence of 43 amino acids was identical in all blowflies and the housefly but differed at 4 positions from the orthologous D. melanogaster sequence. Localization of the mRNA recognized by the respective fragments in the blowflies Lucilia and Phormia resembled that known from D. melanogaster, while Calliphora — the blowfly species lacking rescue activity —showed remarkable differences of localization in both ovarian follicles and the deposited egg cell. This surprising divergence within a morphologically rather uniform family of cyclorrhaphic dipterans should be of interest from both functional and evolutionary points of view.     

Bicoid mRNA diffusion as a mechanism of morphogenesis in Drosophila early development

We show that mRNA diffusion is the main morphogenesis mechanism that consistently explains the establishment of Bicoid protein gradients in the embryo of Drosophila, contradicting the current view of protein diffusion. Moreover, we show that if diffusion for both bicoid mRNA and Bicoid protein were assumed, a steady distribution of Bicoid protein with a constant concentration along the embryo would result, contradicting observations.     

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

SUMMARY 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 homeo- domains 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.     

Pseudocleavage furrows restrict plasma membrane-associated PH domain in syncytial Drosophila embryos

Syncytial cells contain multiple nuclei and have local distribution and function of cellular components despite being synthesized in a common cytoplasm. The syncytial Drosophila blastoderm embryo shows reduced spread of organelle and plasma membrane-associated proteins between adjacent nucleo-cytoplasmic domains. Anchoring to the cytoarchitecture within a nucleo-cytoplasmic domain is likely to decrease the spread of molecules; however, its role in restricting this spread has not been assessed. In order to analyze the cellular mechanisms that regulate the rate of spread of plasma membrane-associated molecules in the syncytial Drosophila embryos, we express a pleckstrin homology (PH) domain in a localized manner at the anterior of the embryo by tagging it with the bicoid mRNA localization signal. Anteriorly expressed PH domain forms an exponential gradient in the anteroposterior axis with a longer length scale compared with Bicoid. Using a combination of experiments and theoretical modeling, we find that the characteristic distribution and length scale emerge due to plasma membrane sequestration and restriction within an energid. Loss of plasma membrane remodeling to form pseudocleavage furrows shows an enhanced spread of PH domain but not Bicoid. Modeling analysis suggests that the enhanced spread of the PH domain occurs due to the increased spread of the cytoplasmic population of the PH domain in pseudocleavage furrow mutants. Our analysis of cytoarchitecture interaction in regulating plasma membrane protein distribution and constraining its spread has implications on the mechanisms of spread of various molecules, such as morphogens in syncytial cells.     

Postgastrular zen expression is required to develop distinct amniotic and serosal epithelia in the scuttle fly Megaselia

The amnioserosa is an extraembryonic epithelium that evolved in higher cyclorrhaphan flies from distinct serosal and amniotic epithelia. The underlying genetic mechanism of this evolutionary transition is unknown. Amnioserosa development of Drosophila correlates with novel expression characteristics of the homeobox gene zerknüllt (zen), including a broad zen expression domain in the syncytial blastoderm and the complete absence of postgastrular zen expression. Here we examine the functional significance of these features by altering the activity profile of zen in Megaselia (a lower cyclorrhaphan fly with distinct serosal and amniotic epithelia) and Drosophila, and by examining in Megaselia the function of u-shaped group (ush-group) genes, which in Drosophila maintain the amnioserosa after gastrulation when zen is no longer expressed. In Megaselia, loss of postgastrular zen expression abrogates serosa development but allows amnion development. Ectopic expression of zen in early Megaselia embryos allows serosa formation but perturbs amnion development. Megaselia homologues of u-shaped group genes are not essential for serosa formation but mediate germband retraction and dorsal closure. Finally, ectopic postgastrular zen expression in Drosophila causes an enlargement of amnioserosa cells and interferes with the morphogenetic functions of the amnioserosa. Our results suggest that the origin of the amnioserosa involved the loss of postgastrular zen expression from extraembryonic tissue, that the early broad expression domain of Drosophila zen evolved afterwards, and that the ush-group genes ancestrally played a role in morphogenetic functions of the amnion.     

The Drosophila Dead end Arf-like3 GTPase controls vesicle trafficking during tracheal fusion cell morphogenesis

The Drosophila larval tracheal system consists of a highly branched tubular organ that becomes interconnected by migration-fusion events during embryonic development. Fusion cells at the tip of each branch guide migration, adhere, and then undergo extensive remodeling as the tracheal lumen extends between the two branches. The Drosophila dead end gene is expressed in fusion cells, and encodes an Arf-like3 GTPase. Analyses of dead end RNAi and mutant embryos reveal that the lumen fails to connect between the two branches. Expression of a constitutively active form of Dead end in S2 cells reveals that it influences the state of actin polymerization, and is present on particles that traffic along actin/microtubule-containing processes. Imaging experiments in vivo reveal that Dead end-containing vesicles are associated with recycling endosomes and the exocyst, and control exocyst localization in fusion cells. These results indicate that the Dead end GTPase plays an important role in trafficking membrane components involved in tracheal fusion cell morphogenesis and lumenal development.     

Strong Purifying Selection on the Odysseus Gene in Two Clades of Sibling Species of the Drosophila montium Species Subgroup

The Odysseus (OdsH) gene was duplicated from its ancestral neuron-expressed gene, unc-4, and then evolved very rapidly under strong positive Darwinian selection as a speciation gene causing hybrid-male sterility between closely related species of the Drosophila simulans clade. Has OdsH also experienced similar positive selection between Drosophila sibling species other than those of the simulans clade? We cloned and sequenced OdsH and unc-4 from two clades of the Drosophila montium species subgroup, the Drosophila lini and the Drosophila kikkawai clades. The ratios of Ka/Ks for OdsH were remarkably low between sibling species of these two clades, suggesting that OdsH has been subjected to strong purifying selection in these two clades.     

Stability and nuclear dynamics of the bicoid morphogen gradient

Patterning in multicellular organisms results from spatial gradients in morphogen concentration, but the dynamics of these gradients remain largely unexplored. We characterize, through in vivo optical imaging, the development and stability of the Bicoid morphogen gradient in Drosophila embryos that express a Bicoid-eGFP fusion protein. The gradient is established rapidly (approximately 1 hr after fertilization), with nuclear Bicoid concentration rising and falling during mitosis. Interphase levels result from a rapid equilibrium between Bicoid uptake and removal. Initial interphase concentration in nuclei in successive cycles is constant (+/-10%), demonstrating a form of gradient stability, but it subsequently decays by approximately 30%. Both direct photobleaching measurements and indirect estimates of Bicoid-eGFP diffusion constants (D < or = 1 microm(2)/s) provide a consistent picture of Bicoid transport on short ( approximately min) time scales but challenge traditional models of long-range gradient formation. A new model is presented emphasizing the possible role of nuclear dynamics in shaping and scaling the gradient.     

The roles of wingless and decapentaplegic in axis and appendage development in the red flour beetle, Tribolium castaneum

Axis patterning and appendage development have been well studied in Drosophila melanogaster, a species in which both limb and segment morphogenesis are derived. In Drosophila, positional information from genes important in anteroposterior and dorsoventral axis formation, including wingless (wg) and decapentaplegic (dpp), is required for allocating and patterning the appendage primordia. We used RNA interference to characterize the functions of wg and dpp in the red flour beetle, Tribolium castaneum, which retains more ancestral modes of limb and segment morphogenesis. We also characterized the expression of potential targets of the WG and DPP signaling pathways in these embryos. Tribolium embryos in which dpp had been downregulated had defects in the dorsalmost body wall, but did not appear to have been globally repatterned and had normal appendages. Downregulation of wg led to the loss of segment boundaries, gnathal and thoracic appendages, and lateral head lobes, and to changes in the expression of dpp, Distal-less, and Engrailed. The functions of wg varied along both the anteroposterior and dorsoventral axes of the embryo. Phylogenetic comparisons indicate that the role of WNT signaling in segment boundary formation is evolutionarily old, but that its role in appendage allocation originated in the common ancestor of holometabolous insects.     

A Two-Dimensional Simulation Model of the Bicoid Gradient in Drosophila

Background

Bicoid (Bcd) is a Drosophila morphogenetic protein responsible for patterning the anterior structures in embryos. Recent experimental studies have revealed important insights into the behavior of this morphogen gradient, making it necessary to develop a model that can recapitulate the biological features of the system, including its dynamic and scaling properties.

Methodology/Principal Findings

We present a biologically realistic 2-D model of the dynamics of the Bcd gradient in Drosophila embryos. This model is based on equilibrium binding of Bcd molecules to non-specific, low affinity DNA sites throughout the Drosophila genome. It considers both the diffusion media within which the Bcd gradient is formed and the dynamic and other relevant properties of bcd mRNA from which Bcd protein is produced. Our model recapitulates key features of the Bcd protein gradient observed experimentally, including its scaling properties and the stability of its nuclear concentrations during development. Our simulation model also allows us to evaluate the effects of other biological activities on Bcd gradient formation, including the dynamic redistribution of bcd mRNA in early embryos. Our simulation results suggest that, in our model, Bcd protein diffusion is important for the formation of an exponential gradient in embryos.

Conclusions/Significance

The 2-D model described in this report is a simple and versatile simulation procedure, providing a quantitative evaluation of the Bcd gradient system. Our results suggest an important role of Bcd binding to non-specific, low-affinity DNA sites in proper formation of the Bcd gradient in our model. They demonstrate that highly complex biological systems can be effectively modeled with relatively few parameters.     

Evolutionary divergence of the paralogs Methoprene tolerant (Met) and germ cell expressed (gce) within the genus Drosophila

Juvenile hormone (JH) signaling underpins both regulatory and developmental pathways in insects. However, the JH receptor is poorly understood. Methoprene tolerant (Met) and germ cell expressed (gce) have been implicated in JH signaling in Drosophila. We investigated the evolution of Met and gce across 12 Drosophila species and found that these paralogs are conserved across at least 63 million years of dipteran evolution. Distinct patterns of selection found using estimates of dN/dS ratios across Drosophila Met and gce coding sequences, along with their incongruent temporal expression profiles in embryonic Drosophila melanogaster, illustrate avenues through which these genes have diverged within the Diptera. Additionally, we demonstrate that the annotated gene CG15032 is the 5′ terminus of gce.In mosquitoes and beetles, a single Met-like homolog displays structural similarity to both Met and gce, and the intron locations are conserved with those of gce. We found that Tribolium and mosquito Met orthologs are assembled from Met- and gce-specific domains in a modular fashion. Our results suggest that Drosophila Met and gce experienced divergent evolutionary pressures following the duplication of an ancestral gce-like gene found in less derived holometabolous insects.     

Stable,Precise, and Reproducible Patterning of Bicoid and Hunchback Molecules in the Early Drosophila Embryo

Precise patterning of morphogen molecules and their accurate reading out are of key importance in embryonic development. Recent experiments have visualized distributions of proteins in developing embryos and shown that the gradient of concentration of Bicoid morphogen in Drosophila embryos is established rapidly after fertilization and remains stable through syncytial mitoses. This stable Bicoid gradient is read out in a precise way to distribute Hunchback with small fluctuations in each embryo and in a reproducible way, with small embryo-to-embryo fluctuation. The mechanisms of such stable, precise, and reproducible patterning through noisy cellular processes, however, still remain mysterious. To address these issues, here we develop the one- and three-dimensional stochastic models of the early Drosophila embryo. The simulated results show that the fluctuation in expression of the hunchback gene is dominated by the random arrival of Bicoid at the hunchback enhancer. Slow diffusion of Hunchback protein, however, averages out this intense fluctuation, leading to the precise patterning of distribution of Hunchback without loss of sharpness of the boundary of its distribution. The coordinated rates of diffusion and transport of input Bicoid and output Hunchback play decisive roles in suppressing fluctuations arising from the dynamical structure change in embryos and those arising from the random diffusion of molecules, and give rise to the stable, precise, and reproducible patterning of Bicoid and Hunchback distributions.     

Antenna and all gnathal appendages are similarly transformed by homothorax knock-down in the cricket Gryllus bimaculatus

Our understanding of the developmental mechanisms underlying the vast diversity of arthropod appendages largely rests on the peculiar case of the dipteran Drosophila melanogaster. In this insect, homothorax (hth) and extradenticle (exd) together play a pivotal role in appendage patterning and identity. We investigated the role of the hth homologue in the cricket Gryllus bimaculatus by parental RNA interference. This species has a more generalized morphology than Oncopeltus fasciatus, the one other insect besides Drosophila where homothorax function has been investigated. The Gryllus head appendages represent the morphologically primitive state including insect-typical mandibles, maxillae and labium, structures highly modified or missing in Oncopeltus and Drosophila. We depleted Gb’hth function through parental RNAi to investigate its requirement for proper regulation of other appendage genes (Gb’wingless, Gb’dachshund, Gb’aristaless and Gb’Distalless) and analyzed the terminal phenotype of Gryllus nymphs. Gb’hth RNAi nymphs display homeotic and segmentation defects similar to hth mutants or loss-of-function clones in Drosophila. Intriguingly, however, we find that in Gb’hth RNAi nymphs not only the antennae but also all gnathal appendages are homeotically transformed, such that all head appendages differentiate distally as legs and proximally as antennae. Hence, Gb’hth is not specifically required for antennal fate, but fulfills a similar role in the specification of all head appendages. This suggests that the role of hth in the insect antenna is not fundamentally different from its function as cofactor of segment-specific homeotic genes in more posterior segments.     

Two ligands signal through the Drosophila PDGF/VEGF receptor to ensure proper salivary gland positioning

The Drosophila embryonic salivary gland is a migrating tissue that undergoes a stereotypic pattern of migration into the embryo. We demonstrate that the migratory path of the salivary gland requires the PDGF/VEGF pathway. The PDGF/VEGF receptor, Pvr, is strongly expressed in the salivary glands, and Pvr mutations cause abnormal ventral curving of the glands, suggesting that Pvr is involved in gland migration. Although the Pvr ligands, Pvf1 and Pvf2, have distinct expression patterns in the Drosophila embryo, mutations for either one of the ligands result in salivary gland migration defects similar to those seen in embryos that lack Pvr. Rescue experiments indicate that the PDGF/VEGF pathway functions autonomously in the salivary gland. The results of this study demonstrate that the Drosophila PDGF/VEGF pathway is essential for proper positioning of the salivary glands.     

A precise Bicoid gradient is nonessential during cycles 11-13 for precise patterning in the Drosophila blastoderm

Background

During development, embryos decode maternal morphogen inputs into highly precise zygotic gene expression. The discovery of the morphogen Bicoid and its profound effect on developmental programming in the Drosophila embryo has been a cornerstone in understanding the decoding of maternal inputs. Bicoid has been described as a classical morphogen that forms a concentration gradient along the antero-posterior axis of the embryo by diffusion and initiates expression of target genes in a concentration-dependent manner in the syncytial blastoderm. Recent work has emphasized the stability of the Bicoid gradient as a function of egg length and the role of nuclear dynamics in maintaining the Bicoid gradient. Bicoid and nuclear dynamics were observed but not modulated under the ideal conditions used previously. Therefore, it has not been tested explicitly whether a temporally stable Bicoid gradient prior to cellularization is required for precise patterning.

Principal Findings

Here, we modulate both nuclear dynamics and the Bicoid gradient using laminar flows of different temperature in a microfluidic device to determine if stability of the Bicoid gradient prior to cellularization is essential for precise patterning. Dramatic motion of both cytoplasm and nuclei was observed prior to cellularization, and the Bicoid gradient was disrupted by nuclear motion and was highly abnormal as a function of egg length. Despite an abnormal Bicoid gradient during cycles 11–13, Even-skipped patterning in these embryos remained precise.

Conclusions

These results indicate that the stability of the Bicoid gradient as a function of egg length is nonessential during syncytial blastoderm stages. Further, presumably no gradient formed by simple diffusion on the scale of egg length could be responsible for the robust antero-posterior patterning observed, as severe cytoplasmic and nuclear motion would disrupt such a gradient. Additional mechanisms for how the embryo could sense its dimensions and interpret the Bicoid gradient are discussed.     

Larval cells become imaginal cells under the control of homothorax prior to metamorphosis in the Drosophila tracheal system

In Drosophila melanogaster, one of the most derived species among holometabolous insects, undifferentiated imaginal cells that are set-aside during larval development are thought to proliferate and replace terminally differentiated larval cells to constitute adult structures. Essentially all tissues that undergo extensive proliferation and drastic morphological changes during metamorphosis are thought to derive from these imaginal cells and not from differentiated larval cells. The results of studies on metamorphosis of the Drosophila tracheal system suggested that large larval tracheal cells that are thought to be terminally differentiated may be eliminated via apoptosis and rapidly replaced by small imaginal cells that go on to form the adult tracheal system. However, the origin of the small imaginal tracheal cells has not been clear. Here, we show that large larval cells in tracheal metamere 2 (Tr2) divide and produce small imaginal cells prior to metamorphosis. In the absence of homothorax gene activity, larval cells in Tr2 become non-proliferative and small imaginal cells are not produced, indicating that homothorax is necessary for proliferation of Tr2 larval cells. These unexpected results suggest that larval cells can become imaginal cells and directly contribute to the adult tissue in the Drosophila tracheal system. During metamorphosis of less derived species of holometabolous insects, adult structures are known to be formed via cells constituting larval structures. Thus, the Drosophila tracheal system may utilize ancestral mode of metamorphosis.