Shared cis-regulatory modules control expression of the tandem paralogs midline and H15 in the follicular epithelium.

The posterior end of the follicular epithelium is patterned by midline (MID) and its paralog H15, the Drosophila homologs of the mammalian Tbx20 transcription factor. We have previously identified two cis-regulatory modules (CRMs) that recapitulate the endogenous pattern of mid in the follicular epithelium. Here, using CRISPR/Cas9 genome editing, we demonstrate redundant activity of these mid CRMs. Although the deletion of either CRM alone generated marginal change in mid expression, the deletion of both CRMs reduced expression by 60%. Unexpectedly, the deletion of the 5' proximal CRM of mid eliminated H15 expression. Interestingly, expression of these paralogs in other tissues remained unaffected in the CRM deletion backgrounds. These results suggest that the paralogs are regulated by a shared CRM that coordinates gene expression during posterior fate determination. The consistent overlapping expression of mid and H15 in various tissues may indicate that the paralogs could also be under shared regulation by other CRMs in these tissues.

Nebulous without white: annotated long-read genome assembly and CRISPR/Cas9 genome engineering in Drosophila nebulosa.

The diversity among Drosophila species presents an opportunity to study the molecular mechanisms underlying the evolution of biological phenomena. A challenge to investigating these species is that, unlike the plethora of molecular and genetics tools available for D. melanogaster research, many other species do not have sequenced genomes; a requirement for employing these tools. Selecting transgenic flies through white (w) complementation has been commonly practiced in numerous Drosophila species. While tolerated, the disruption of w is associated with impaired vision, among other effects in D. melanogaster. The D. nebulosa fly has a unique mating behavior which requires vision, and is thus unable to successfully mate in dark conditions. Here, we hypothesized that the disruption of w will impede mating success. As a first step, using PacBio long-read sequencing, we assembled a high-quality annotated genome of D. nebulosa. Using these data, we employed CRISPR/Cas9 to successfully disrupt the w gene. As expected, D. nebulosa males null for w did not court females, unlike several other mutant strains of Drosophila species whose w gene has been disrupted. In the absence of mating, no females became homozygous null for w. We conclude that gene disruption via CRISPR/Cas9 genome engineering is a successful tool in D. nebulosa, and that the w gene is necessary for mating. Thus, an alternative selectable marker unrelated to vision is desirable.

The ETS-transcription factor Pointed is sufficient to regulate the posterior fate of the follicular epithelium.

The Janus-kinase/signal transducer and activator of transcription (JAK/STAT) pathway regulates the anterior posterior axis of the Drosophila follicle cells. In the anterior, it activates the bone morphogenetic protein (BMP) signaling pathway through expression of the BMP ligand decapentaplegic (dpp). In the posterior, JAK/STAT works with the epidermal growth factor receptor (EGFR) pathway to express the T-box transcription factor midline (mid). Although MID is necessary for establishing the posterior fate of the egg chamber, we show that it is not sufficient to determine a posterior fate. The ETS-transcription factor pointed (pnt) is expressed in an overlapping domain to mid in the follicle cells. This study shows that pnt is upstream of mid and that it is sufficient to induce a posterior fate in the anterior end, which is characterized by the induction of mid, the prevention of the stretched cells formation and the abrogation of border cell migration. We demonstrate that the anterior BMP signaling is abolished by PNT through dpp repression. However, ectopic DPP cannot rescue the anterior fate formation, suggesting additional targets of PNT participate in the posterior fate determination.

Quantitative analyses of EGFR localization and trafficking dynamics in the follicular epithelium.

To bridge the gap between qualitative and quantitative analyses of the epidermal growth factor receptor (EGFR) in tissues, we generated an sfGFP-tagged EGF receptor (EGFR-sfGFP) in Drosophila The homozygous fly appears similar to wild type with EGFR expression and activation patterns that are consistent with previous reports in the ovary, early embryo, and imaginal discs. Using ELISA, we quantified an average of 1100, 6200 and 2500 receptors per follicle cell (FC) at stages 8/9, 10 and ≥11 of oogenesis, respectively. Interestingly, the spatial localization of the EGFR to the apical side of the FCs at early stages depended on the TGFα-like ligand Gurken. At later stages, EGFR localized to basolateral positions of the FCs. Finally, we followed the endosomal localization of EGFR in the FCs. The EGFR colocalized with the late endosome, but no significant colocalization of the receptor was found with the early endosome. The EGFR-sfGFP fly is an exciting new resource for studying cellular localization and regulation of EGFR in tissues.

In locus analysis of patterning evolution of the BMP type II receptor Wishful thinking.

Proper tissue patterning is an essential step during organ formation. During this process, genes are expressed in distinct patterns, defining boundaries for future functional domains. The bone morphogenetic protein (BMP) signaling pathway sets the anterior domain during eggshell patterning. Previously, the Drosophila melanogaster homolog of BMPR2, Wishful thinking (WIT), was shown to be required for BMP signaling and patterning during eggshell formation. Expressed in a conserved anterior pattern, the width of wit patterning in the follicular epithelium is evolutionarily divergent between Drosophila species. We used genome editing to demonstrate how the gene pattern divergence is controlled in cis within the wit locus of D. virilis Furthermore, unlike direct targets of BMP signaling, we demonstrate how one transcription factor binding site shapes the pattern of WIT in D. melanogaster by negative regulation. However, changes in this site are not sufficient to explain the evolution of wit patterning, suggesting that a positive regulatory element that controls pattern divergence remains to be discovered.

Simple Expression Domains Are Regulated by Discrete CRMs During Drosophila Oogenesis.

Eggshell patterning has been extensively studied in Drosophila melanogaster However, the cis-regulatory modules (CRMs), which control spatiotemporal expression of these patterns, are vastly unexplored. The FlyLight collection contains >7000 intergenic and intronic DNA fragments that, if containing CRMs, can drive the transcription factor GAL4. We cross-listed the 84 genes known to be expressed during D. melanogaster oogenesis with the ∼1200 listed genes of the FlyLight collection, and found 22 common genes that are represented by 281 FlyLight fly lines. Of these lines, 54 show expression patterns during oogenesis when crossed to an UAS-GFP reporter. Of the 54 lines, 16 recapitulate the full or partial pattern of the associated gene pattern. Interestingly, while the average DNA fragment size is ∼3 kb in length, the vast majority of fragments show one type of spatiotemporal pattern in oogenesis. Mapping the distribution of all 54 lines, we found a significant enrichment of CRMs in the first intron of the associated genes' model. In addition, we demonstrate the use of different anteriorly active FlyLight lines as tools to disrupt eggshell patterning in a targeted manner. Our screen provides further evidence that complex gene patterns are assembled combinatorially by different CRMs controlling the expression of genes in simple domains.

Gene regulation during Drosophila eggshell patterning.

A common path to the formation of complex 3D structures starts with an epithelial sheet that is patterned by inductive cues that control the spatiotemporal activities of transcription factors. These activities are then interpreted by the cis-regulatory regions of the genes involved in cell differentiation and tissue morphogenesis. Although this general strategy has been documented in multiple developmental contexts, the range of experimental models in which each of the steps can be examined in detail and evaluated in its effect on the final structure remains very limited. Studies of the Drosophila eggshell patterning provide unique insights into the multiscale mechanisms that connect gene regulation and 3D epithelial morphogenesis. Here we review the current understanding of this system, emphasizing how the recent identification of cis-regulatory regions of genes within the eggshell patterning network enables mechanistic analysis of its spatiotemporal dynamics and evolutionary diversification. It appears that cis-regulatory changes can account for only some aspects of the morphological diversity of Drosophila eggshells, such as the prominent differences in the number of the respiratory dorsal appendages. Other changes, such as the appearance of the respiratory eggshell ridges, are caused by changes in the spatial distribution of inductive signals. Both types of mechanisms are at play in this rapidly evolving system, which provides an excellent model of developmental patterning and morphogenesis.

Control of reaction-diffusion equations on time-evolving manifolds.

Among the main actors of organism development there are morphogens, which are signaling molecules diffusing in the developing organism and acting on cells to produce local responses. Growth is thus determined by the distribution of such signal. Meanwhile, the diffusion of the signal is itself affected by the changes in shape and size of the organism. In other words, there is a complete coupling between the diffusion of the signal and the change of the shapes. In this paper, we introduce a mathematical model to investigate such coupling. The shape is given by a manifold, that varies in time as the result of a deformation given by a transport equation. The signal is represented by a density, diffusing on the manifold via a diffusion equation. We show the non-commutativity of the transport and diffusion evolution by introducing a new concept of Lie bracket between the diffusion and the transport operator. We also provide numerical simulations showing this phenomenon.

Evolutionary changes in TGFα distribution underlie morphological diversity in eggshells from Drosophila species.

Drosophila eggshells display remarkable morphological diversity among species; however, the molecular origin of this structural diversification is mostly unknown. Here, we analyzed the dorsal ridge (DR), a lumen-like structure along the dorsal side of eggshells, from numerous Drosophila species. This structure varies in length and width across species, and is absent from D. melanogaster eggshells. We associated DR formation with distinct spatiotemporal changes in epidermal growth factor receptor (EGFR) activation, which acts as a key receptor in eggshell patterning. We show that changes in the distribution of the TGFα-like ligand Gurken (GRK), a crucial ligand for axis formation, underlies EGFR activation and DR formation in D. willistoni. Furthermore, we demonstrate that GRK from D. willistoni rescues a grk-null D. melanogaster fly and, remarkably, it is also sufficient to generate a DR-like structure on its eggshell.

Chorion patterning: a window into gene regulation and Drosophila species' relatedness.

Changes in gene regulation are associated with the evolution of morphologies. However, the specific sequence information controlling gene expression is largely unknown and discovery is time and labor consuming. We use the intricate patterning of follicle cells to probe species' relatedness in the absence of sequence information. We focus on one of the major families of genes that pattern the Drosophila eggshell, the Chorion protein (Cp). Systematically screening for the spatiotemporal patterning of all nine Cp genes in three species (Drosophila melanogaster, D. nebulosa, and D. willistoni), we found that most genes are expressed dynamically during mid and late stages of oogenesis. Applying an annotation code, we transformed the data into binary matrices that capture the complexity of gene expression. Gene patterning is sufficient to predict species' relatedness, consistent with their phylogeny. Surprisingly, we found that expression domains of most genes are different among species, suggesting that Cp regulation is rapidly evolving. In addition, we found a morphological novelty along the dorsalmost side of the eggshell, the dorsal ridge. Our matrix analysis placed the dorsal ridge domain in a cluster of epidermal growth factor receptor associated domains, which was validated through genetic and chemical perturbations. Expression domains are regulated cooperatively or independently by signaling pathways, supporting that complex patterns are combinatorially assembled from simple domains.

The Drosophila BMPRII, wishful thinking, is required for eggshell patterning.

The Drosophila eggshell is an elaborate structure that is derived from a monolayer of follicular epithelium surrounding the developing oocyte within the female ovary. The bone morphogenetic protein (BMP) signaling pathway is essential for controlling the patterning and morphogenesis of the eggshell. During oogenesis, the roles of patterning and morphogenesis by the BMP type I receptor thickveins (tkv) have been studied extensively. However, signaling through this pathway requires both type I and II receptors, and the latter has yet to be established in oogenesis. We focus on wishful thinking (wit), the Drosophila homolog to the mammalian BMP type II receptor, BMPRII. We found that wit is expressed dynamically in the FCs of D. melanogaster in an evolutionary conserved pattern. The expression patterns are highly correlated with the dynamics of the BMP signaling, which is consistent with our finding that wit is a target of BMP signaling. Furthermore, we established that WIT is necessary for BMP signaling, and loss of WIT is associated with cell autonomous loss of BMP responses. Of importance, we demonstrated that perturbations in WIT led to changes in eggshell morphologies in domains that are patterned by BMP signaling. Previous studies have shown a role for WIT in BMP signaling during neurogenesis; however, our results reveal a role for WIT in epithelial cells' development.

Evolution of BMP signaling in Drosophila oogenesis: a receptor-based mechanism.

The bone morphogenetic protein (BMP) signaling pathway is a conserved regulator of cellular and developmental processes in animals. The mechanisms underlying BMP signaling activation differ among tissues and mostly reflect changes in the expression of pathway components. BMP signaling is one of the major pathways responsible for the patterning of the Drosophila eggshell, a complex structure derived from a layer of follicle cells (FCs) surrounding the developing oocyte. Activation of BMP signaling in the FCs is dynamic. Initially, signaling is along the anterior-posterior (A/P) axis; later, signaling acquires dorsal-ventral (D/V) polarity. These dynamics are regulated by changes in the expression pattern of the type I BMP receptor thickveins (tkv). We recently found that signaling dynamics and TKV patterning are highly correlated in the FCs of multiple Drosophila species. In addition, we showed that signaling patterns are spatially different among species. Here, we use a mathematical model to simulate the dynamics and differences of BMP signaling in numerous species. This model predicts that qualitative and quantitative changes in receptor expression can lead to differences in the spatial pattern of BMP signaling. We tested these predications experimentally in three different Drosophila species and through genetic perturbations of BMP signaling in D. melanogaster. On the basis of our results, we concluded that changes in tkv patterning can account for the experimentally observed differences in the patterns of BMP signaling in multiple Drosophila species.

Pattern formation by a moving morphogen source.

During Drosophila melanogaster oogenesis, the follicular epithelium that envelops the germline cyst gives rise to an elaborate eggshell, which houses the future embryo and mediates its interaction with the environment. A prominent feature of the eggshell is a pair of dorsal appendages, which are needed for embryo respiration. Morphogenesis of this structure depends on broad, a zinc-finger transcription factor, regulated by the EGFR pathway. While much has been learned about the mechanisms of broad regulation by EGFR, current understanding of processes that shape the spatial pattern of broad expression is incomplete. We propose that this pattern is defined by two different phases of EGFR activation: an early, posterior-to-anterior gradient of EGFR signaling sets the posterior boundary of broad expression, while the anterior boundary is set by a later phase of EGFR signaling, distributed in a dorsoventral gradient. This model can explain the wild-type pattern of broad in D. melanogaster, predicts how this pattern responds to genetic perturbations, and provides insight into the mechanisms driving diversification of eggshell patterning. The proposed model of the broad expression pattern can be used as a starting point for the quantitative analysis of a large number of gene expression patterns in Drosophila oogenesis.

BMP signaling dynamics in the follicle cells of multiple Drosophila species.

The dorsal anterior region of the follicle cells (FCs) in the developing Drosophila egg gives rise to the respiratory eggshell appendages. These tubular structures display a wide range of qualitative and quantitative variations across Drosophila species, providing a remarkable example of a rapidly evolving morphology. In D. melanogaster, the bone morphogenetic protein (BMP) signaling pathway is an important regulator of FCs patterning and dorsal appendages morphology. To explore the mechanisms underlying the diversification of eggshell patterning, we analyzed BMP signaling in the FCs of 16 Drosophila species that span 45 million years of evolution. We found that the spatial patterns of BMP signaling in the FCs are dynamic and exhibit a range of interspecies' variations. In most of the species examined, the dynamics of BMP signaling correlate with the expression of the type I BMP receptor thickveins (tkv). This correlation suggests that interspecies' variations of tkv expression are responsible for the diversification of BMP signaling during oogenesis. This model was supported by genetic manipulations of tkv expression in the FCs of D. melanogaster that successfully recapitulated the signaling diversities found in the other species. Our results suggest that regulation of receptor expression mediates spatial diversification of BMP signaling in Drosophila oogenesis, and they provide insight into a mechanism underlying the evolution of eggshell patterning.

Pattern formation by dynamically interacting network motifs.

Systematic validation of pattern formation mechanisms revealed by molecular studies of development is essentially impossible without mathematical models. Models can provide a compact summary of a large number of experiments that led to mechanism formulation and guide future studies of pattern formation. Here, we realize this program by analyzing a mathematical model of epithelial patterning by the highly conserved EGFR and BMP signaling pathways in Drosophila oogenesis. The model accounts for the dynamic interaction of the feedforward and feedback network motifs that control the expression of Broad, a zinc finger transcription factor expressed in the cells that form the upper part of the respiratory eggshell appendages. Based on the combination of computational analysis and genetic experiments, we show that the model accounts for the key features of wild-type pattern formation, correctly predicts patterning defects in multiple mutants, and guides the identification of additional regulatory links in a complex pattern formation mechanism.

Expression patterns of cadherin genes in Drosophila oogenesis.

In Drosophila oogenesis, the follicular epithelium that envelops the oocyte is patterned by a small set of inductive signals and gives rise to an elaborate three-dimensional eggshell. Several eggshell structures provide sensitive readouts of the patterning signals, but the formation of these structures is still poorly understood. In other systems, epithelial morphogenesis is guided by the spatial patterning of cell adhesion and cytoskeleton genes. As a step towards developing a comprehensive description of patterning events leading to eggshell morphogenesis, we report the expression of Drosophila cadherins, calcium-dependent adhesion molecules that are repeatedly used throughout development. We found that 9/17 of Drosophila cadherins are expressed in the follicular epithelium in dynamic patterns during oogenesis. In late oogenesis, the expression patterns of cadherin genes in the main body follicle cells is summarized using a compact set of simple geometric shapes, reflecting the integration of the EGFR and DPP inductive signals. The multi-layered composite patterning of the cadherins is hypothesized to play a key role in the formation of the eggshell. Of particular note is the complex patterning of the region of the follicular epithelium that gives rise to the dorsal appendages, which are tubular structures that serve as respiratory organs for the developing embryo.

A combinatorial code for pattern formation in Drosophila oogenesis.

Two-dimensional patterning of the follicular epithelium in Drosophila oogenesis is required for the formation of three-dimensional eggshell structures. Our analysis of a large number of published gene expression patterns in the follicle cells suggests that they follow a simple combinatorial code based on six spatial building blocks and the operations of union, difference, intersection, and addition. The building blocks are related to the distribution of inductive signals, provided by the highly conserved epidermal growth factor receptor and bone morphogenetic protein signaling pathways. We demonstrate the validity of the code by testing it against a set of patterns obtained in a large-scale transcriptional profiling experiment. Using the proposed code, we distinguish 36 distinct patterns for 81 genes expressed in the follicular epithelium and characterize their joint dynamics over four stages of oogenesis. The proposed combinatorial framework allows systematic analysis of the diversity and dynamics of two-dimensional transcriptional patterns and guides future studies of gene regulation.

Spatial regulation of BMP signaling by patterned receptor expression.

Local delivery of TGF-beta/BMP ligands is commonly used as a tissue engineering strategy for the spatial regulation of cell growth and differentiation. While the location and the dose of ligand are the only parameters that influence the spatial distribution and biological effects of the ligand in vitro, in vivo genetic studies of development reveal that spatial control of TGF-beta/BMP signaling can be accomplished at multiple levels, from ligand release to signal interpretation. Here we focus on spatial control of BMP signaling by patterned receptor expression. Motivated by our recent experimental analysis of the two-dimensional BMP signaling patterns in the developing Drosophila egg, we formulate one- and two-dimensional models of ligand diffusion and internalization in the presence of patterned receptor expression. Our analysis of these models shows that they can capture the quantitative features of the experimentally observed pattern of phosphorylated SMAD in Drosophila oogenesis and shows that patterned receptor expression provides versatile control of BMP signaling in developing tissues. Quantitative understanding of the mechanisms of spatiotemporal control of signaling pathways in development is essential for successful harnessing of these pathways in tissue engineering.

Cad74A is regulated by BR and is required for robust dorsal appendage formation in Drosophila oogenesis.

Drosophila egg development is an established model for studying epithelial patterning and morphogenesis, but the connection between signaling pathways and egg morphology is still incompletely understood. We have identified a non-classical cadherin, Cad74A, as a putative adhesion gene that bridges epithelial patterning and morphogenesis in the follicle cells. Starting in mid-oogenesis, Cad74A is expressed in the follicle cells that contact the oocyte, including the border cells and most of the columnar follicle cells. However, Cad74A is repressed in two dorsolateral patches of follicle cells, which participate in the formation of tubular respiratory appendages. We show genetically that Cad74A is downstream of the EGFR and BMP signaling pathways and is repressed by the Zn-finger transcription factor Broad. The correlation of Cad74A repression in the cells that bend out of the plane of the follicular epithelium is preserved across Drosophila species and mutant backgrounds exhibiting a range of eggshell phenotypes. Complete removal of Cad74A from the follicle cells causes defects in dorsal appendage formation. Ectopic expression of Cad74A in the roof cells results in shortened, flattened appendages due to the hindered migration of the roof cells. Based on these results, we propose that Cad74A is part of the adhesive machinery that enables robust dorsal appendage formation, and as such provides a link between the patterning of the follicle cells and eggshell morphogenesis.

Drosophila eggshell is patterned by sequential action of feedforward and feedback loops.

During Drosophila oogenesis, patterning activities of the EGFR and Dpp pathways specify several subpopulations of the follicle cells that give rise to dorsal eggshell structures. The roof of dorsal eggshell appendages is formed by the follicle cells that express Broad (Br), a zinc-finger transcription factor regulated by both pathways. EGFR induces Br in the dorsal follicle cells. This inductive signal is overridden in the dorsal midline cells, which are exposed to high levels of EGFR activation, and in the anterior cells, by Dpp signaling. We show that the resulting changes in the Br pattern affect the expression of Dpp receptor thickveins (tkv), which is essential for Dpp signaling. By controlling tkv, Br controls Dpp signaling in late stages of oogenesis and, as a result, regulates its own repression in a negative-feedback loop. We synthesize these observations into a model, whereby the dynamics of Br expression are driven by the sequential action of feedforward and feedback loops. The feedforward loop controls the spatial pattern of Br expression, while the feedback loop modulates this pattern in time. This mechanism demonstrates how complex patterns of gene expression can emerge from simple inputs, through the interaction of regulatory network motifs.