Whole-genome ChIP-chip analysis of Dorsal, Twist, and Snail suggests integration of diverse patterning processes in the Drosophila embryo

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest that it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how "Coordinator," a long DNA motif composed of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines the regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, whereas HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in the shared regulation of genes involved in cell-type and positional identities and ultimately shapes facial morphology and evolution.

Tissue-specific RNA Polymerase II promoter-proximal pause release and burst kinetics in a Drosophila embryonic patterning network

BACKGROUND

Formation of tissue-specific transcriptional programs underlies multicellular development, including dorsoventral (DV) patterning of the Drosophila embryo. This involves interactions between transcriptional enhancers and promoters in a chromatin context, but how the chromatin landscape influences transcription is not fully understood.

RESULTS

Here we comprehensively resolve differential transcriptional and chromatin states during Drosophila DV patterning. We find that RNA Polymerase II pausing is established at DV promoters prior to zygotic genome activation (ZGA), that pausing persists irrespective of cell fate, but that release into productive elongation is tightly regulated and accompanied by tissue-specific P-TEFb recruitment. DV enhancers acquire distinct tissue-specific chromatin states through CBP-mediated histone acetylation that predict the transcriptional output of target genes, whereas promoter states are more tissue-invariant. Transcriptome-wide inference of burst kinetics in different cell types revealed that while DV genes are generally characterized by a high burst size, either burst size or frequency can differ between tissues.

CONCLUSIONS

The data suggest that pausing is established by pioneer transcription factors prior to ZGA and that release from pausing is imparted by enhancer chromatin state to regulate bursting in a tissue-specific manner in the early embryo. Our results uncover how developmental patterning is orchestrated by tissue-specific bursts of transcription from Pol II primed promoters in response to enhancer regulatory cues.

Dual role of Snail1 as transcriptional repressor and activator

Snail1 transcriptional factor plays a key role in the control of epithelial to mesenchymal transition, a process that remodels tumor cells increasing their invasion and chemo-resistance as well as reprograms their metabolism and provides stemness properties. During this transition, Snail1 acts as a transcriptional repressor and, as growing evidences have demonstrated, also as a direct activator of mesenchymal genes. In this review, I describe the different proteins that interact with Snail1 and are responsible for these two different functions on gene expression; I focus on the transcriptional factors that associate to Snail1 in their target promoters, both activated and repressed. I also present working models for Snail1 action both as repressor and activator and raise some issues that still need to be investigated.

Prediction and analysis of cis-regulatory elements in Dorsal and Ventral patterning genes of Tribolium castaneum and its comparison with Drosophila melanogaster

Insect embryonic development and morphology are characterized by their anterior-posterior and dorsal-ventral (DV) patterning. In Drosophila embryos, DV patterning is mediated by a dorsal protein gradient which activates twist and snail proteins, the important regulators of DV patterning. To activate or repress gene expression, some regulatory proteins bind in clusters to their target gene at sites known as cis-regulatory elements or enhancers. To understand how variations in gene expression in different lineages might lead to different phenotypes, it is necessary to understand enhancers and their evolution. Drosophila melanogaster has been widely studied to understand the interactions between transcription factors and the transcription factor binding sites. Tribolium castaneum is an upcoming model animal which is catching the interest of biologists and the research on the enhancer mechanisms in the insect's axes patterning is still in infancy. Therefore, the current study was designed to compare the enhancers of DV patterning in the two insect species. The sequences of ten proteins involved in DV patterning of D. melanogaster were obtained from Flybase. The protein sequences of T. castaneum orthologous to those obtained from D. melanogaster were acquired from NCBI BLAST, and these were then converted to DNA sequences which were modified by adding 20 kb sequences both upstream and downstream to the gene. These modified sequences were used for further analysis. Bioinformatics tools (Cluster-Buster and MCAST) were used to search for clusters of binding sites (enhancers) in the modified DV genes. The results obtained showed that the transcription factors in Drosophila melanogaster and Tribolium castaneum are nearly identical; however, the number of binding sites varies between the two species, indicating transcription factor binding site evolution, as predicted by two different computational tools. It was observed that dorsal, twist, snail, zelda, and Supressor of Hairless are the transcription factors responsible for the regulation of DV patterning in the two insect species.

Modulation of protein-DNA binding reveals mechanisms of spatiotemporal gene control in early Drosophila embryos

It is well known that enhancers regulate the spatiotemporal expression of their target genes by recruiting transcription factors (TFs) to the cognate binding sites in the region. However, the role of multiple binding sites for the same TFs and their specific spatial arrangement in determining the overall competency of the enhancer has yet to be fully understood. In this study, we utilized the MS2-MCP live imaging technique to quantitatively analyze the regulatory logic of the snail distal enhancer in early Drosophila embryos. Through systematic modulation of Dorsal and Twist binding motifs in this enhancer, we found that a mutation in any one of these binding sites causes a drastic reduction in transcriptional amplitude, resulting in a reduction in mRNA production of the target gene. We provide evidence of synergy, such that multiple binding sites with moderate affinities cooperatively recruit more TFs to drive stronger transcriptional activity than a single site. Moreover, a Hidden Markov-based stochastic model of transcription reveals that embryos with mutated binding sites have a higher probability of returning to the inactive promoter state. We propose that TF-DNA binding regulates spatial and temporal gene expression and drives robust pattern formation by modulating transcriptional kinetics and tuning bursting rates.

Chromatin accessibility in the Drosophila embryo is determined by transcription factor pioneering and enhancer activation

Chromatin accessibility is integral to the process by which transcription factors (TFs) read out cis-regulatory DNA sequences, but it is difficult to differentiate between TFs that drive accessibility and those that do not. Deep learning models that learn complex sequence rules provide an unprecedented opportunity to dissect this problem. Using zygotic genome activation in Drosophila as a model, we analyzed high-resolution TF binding and chromatin accessibility data with interpretable deep learning and performed genetic validation experiments. We identify a hierarchical relationship between the pioneer TF Zelda and the TFs involved in axis patterning. Zelda consistently pioneers chromatin accessibility proportional to motif affinity, whereas patterning TFs augment chromatin accessibility in sequence contexts where they mediate enhancer activation. We conclude that chromatin accessibility occurs in two tiers: one through pioneering, which makes enhancers accessible but not necessarily active, and the second when the correct combination of TFs leads to enhancer activation.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how 'Coordinator', a long DNA motif comprised of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, while HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in shared regulation of genes involved in cell-type and positional identities, and ultimately shapes facial morphology and evolution.

Cis-regulatory atlas of primary human CD4+ T cells

Cis-regulatory elements (CRE) are critical for coordinating gene expression programs that dictate cell-specific differentiation and homeostasis. Recently developed self-transcribing active regulatory region sequencing (STARR-Seq) has allowed for genome-wide annotation of functional CREs. Despite this, STARR-Seq assays are only employed in cell lines, in part, due to difficulties in delivering reporter constructs. Herein, we implemented and validated a STARR-Seq-based screen in human CD4+ T cells using a non-integrating lentiviral transduction system. Lenti-STARR-Seq is the first example of a genome-wide assay of CRE function in human primary cells, identifying thousands of functional enhancers and negative regulatory elements (NREs) in human CD4+ T cells. We find an unexpected difference in nucleosome organization between enhancers and NRE: enhancers are located between nucleosomes, whereas NRE are occupied by nucleosomes in their endogenous locations. We also describe chromatin modification, eRNA production, and transcription factor binding at both enhancers and NREs. Our findings support the idea of silencer repurposing as enhancers in alternate cell types. Collectively, these data suggest that Lenti-STARR-Seq is a successful approach for CRE screening in primary human cell types, and provides an atlas of functional CREs in human CD4+ T cells.

The epithelial polarity genes frazzled and GUK-holder adjust morphogen gradients to coordinate changes in cell position with cell fate specification

Morphogenetic gradients specify distinct cell populations within tissues. Originally, morphogens were conceived as substances that act on a static field of cells, yet cells usually move during development. Thus, the way cell fates are defined in moving cells remains a significant and largely unsolved problem. Here, we investigated this issue using spatial referencing of cells and 3D spatial statistics in the Drosophila blastoderm to reveal how cell density responds to morphogenetic activity. We show that the morphogen decapentaplegic (DPP) attracts cells towards its peak levels in the dorsal midline, whereas dorsal (DL) stalls them ventrally. We identified frazzled and GUK-holder as the downstream effectors regulated by these morphogens that constrict cells and provide the mechanical force necessary to draw cells dorsally. Surprisingly, GUKH and FRA modulate the DL and DPP gradient levels and this regulation creates a very precise mechanism of coordinating cell movement and fate specification.

Modulation of protein-DNA binding reveals mechanisms of spatiotemporal gene control in early Drosophila embryos

It is well known that enhancers regulate the spatiotemporal expression of their target genes by recruiting transcription factors (TFs) to the cognate binding sites in the region. However, the role of multiple binding sites for the same TFs and their specific spatial arrangement in determining the overall competency of the enhancer has yet to be fully understood. In this study, we utilized the MS2-MCP live imaging technique to quantitatively analyze the regulatory logic of the snail distal enhancer in early Drosophila embryos. Through systematic modulation of Dorsal and Twist binding motifs in this enhancer, we found that a mutation in any one of these binding sites causes a drastic reduction in transcriptional amplitude, resulting in a reduction in total mRNA production of the target gene. We provide evidence of synergy, such that multiple binding sites with moderate affinities cooperatively recruit more TFs to drive stronger transcriptional activity than a single site. Moreover, a Hidden Markov-based stochastic model of transcription reveals that embryos with mutated binding sites have a higher probability of returning to the inactive promoter state. We propose that TF-DNA binding regulates spatial and temporal gene expression and drives robust pattern formation by modulating transcriptional kinetics and tuning bursting rates.

Prognostic significance of circulating tumor cell measurement in the peripheral blood of patients with nasopharyngeal carcinoma

OBJECTIVE

Nasopharyngeal Carcinoma (NPC) is lethal cancer. Typically, relapse and metastasis are the outcomes of most patients. Against this backdrop, this study aimed to investigate the correlation between Circulating Tumor Cell (CTC) profiles and clinicopathological features in patients with NPC.

PATIENTS AND METHODS

A total of 119 blood samples from 79 patients were collected from patients with NPC during treatment. CanPatrol^TM CTC enrichment and RNA In Situ Hybridization (RNA-ISH) were used to characterize CTCs, including epithelial, Mesenchymal (MCTCs), and epithelial/mesenchymal mixed types according to their surface markers.

RESULTS

The number of CTCs and MCTCs in the pre-treatment group was significantly higher than that in the post-treatment group (p < 0.05). The total number of CTCs and MCTCs cell numbers was significant correlation with Tumor-Node-Metastasis (TNM) staging (p < 0.05), Progression-Free Survival (PFS), and Overall Survival (OS). The PFS of patients with > 7 CTCs or > 5 MCTCs per 5 mL blood was significantly shorter PFS than those patients with ≤ 7 CTCs or ≤ 5 MCTCs (p < 0.05). Patients treated with targeted therapy combined with chemoradiotherapy had poorer PFS and OS rates than those treated with chemoradiotherapy (p < 0.05). The Kaplan-Meier survival analysis also demonstrated that patients with changes in CTC > 4 were strongly associated with PFS and OS rates (p < 0.05).

CONCLUSION

CTC and MCTC number detection in patients with NPC is a useful biomarker for predicting patient progress. Patients with more than 7 CTCs or 5 MCTCs in 5 mL of blood had shorter PFS and OS rates. CTC and MCTC count changes were also significantly associated with the patient's therapy.

Identification of sex-biased and neurodevelopment genes via brain transcriptome in Ostrinia furnacalis

Insect brains play important roles in the regulation of sex-biased behaviors such as mating and oviposition. The neural structure and function of brain differences between males and females have been identified, in which the antenna lobes (AL) showed the most discrepancy, however, the whole repertoire of the genes expressed in the brains and the molecular mechanism of neural signaling and structural development are still unclear. In this study, high-throughput transcriptome analysis of male and female brains was carried on in the Asia corn borer, Ostrinia furnacalis, and a total of 39.23 Gb data and 34,092 unigenes were obtained. Among them, 276 genes displayed sex-biased expression by DEG analysis, of which 125 genes were highly expressed in the males and 151 genes were highly expressed in the females. Besides, by homology analysis against genes that have been confirmed to be related to brain neurodevelopment, a total of 24 candidate genes were identified in O. furnacalis. In addition, to further screen the core genes that may be important for sex-biased nerve signaling and neurodevelopment, protein-protein interaction networks were constructed for the sex-biased genes and neurodevelopment genes. We identified 10 (Mhc, Mlc1, Mlc2, Prm, Mf, wupA, TpnC25D, fln, l(2)efl, and Act5C), 11 (PPO2, GNBP3, Spn77Ba, Ppn, yellow-d2, PGRP-LB, PGRP-SD, PGRP-SC2, Hml, Cg25C, and vkg) and 8 (dac, wg, hh, ci, run, Lim1, Rbp9, and Bx) core hub genes that may be related to brain neural development from male-biased, female-biased, and neurodevelopment gene groups. Our results provide a reference for further analysis of the dimorphism of male and female brain structures in agricultural pests.

SUMOylation of Dorsal attenuates Toll/NF-κB signaling

In Drosophila, Toll/NF-κB signaling plays key roles in both animal development and in host defense. The activation, intensity, and kinetics of Toll signaling are regulated by posttranslational modifications such as phosphorylation, SUMOylation, or ubiquitination that target multiple proteins in the Toll/NF-κB cascade. Here, we have generated a CRISPR-Cas9 edited Dorsal (DL) variant that is SUMO conjugation resistant. Intriguingly, embryos laid by dlSCR mothers overcome dl haploinsufficiency and complete the developmental program. This ability appears to be a result of higher transcriptional activation by DLSCR. In contrast, SUMOylation dampens DL transcriptional activation, ultimately conferring robustness to the dorso-ventral program. In the larval immune response, dlSCR animals show an increase in crystal cell numbers, stronger activation of humoral defense genes, and high cactus levels. A mathematical model that evaluates the contribution of the small fraction of SUMOylated DL (1-5%) suggests that it acts to block transcriptional activation, which is driven primarily by DL that is not SUMO conjugated. Our findings define SUMO conjugation as an important regulator of the Toll signaling cascade, in both development and host defense. Our results broadly suggest that SUMO attenuates DL at the level of transcriptional activation. Furthermore, we hypothesize that SUMO conjugation of DL may be part of a Ubc9-dependent mechanism that restrains Toll/NF-κB signaling.

A coherent FOXO3-SNAI2 feed-forward loop in autophagy

Understanding autophagy regulation is instrumental in developing therapeutic interventions for autophagy-associated disease. Here, we identified SNAI2 as a regulator of autophagy from a genome-wide screen in HeLa cells. Upon energy stress, SNAI2 is transcriptionally activated by FOXO3 and interacts with FOXO3 to form a feed-forward regulatory loop to reinforce the expression of autophagy genes. Of note, SNAI2-increased FOXO3-DNA binding abrogates CRM1-dependent FOXO3 nuclear export, illuminating a pivotal role of DNA in the nuclear retention of nucleocytoplasmic shuttling proteins. Moreover, a dFoxO-Snail feed-forward loop regulates both autophagy and cell size in Drosophila, suggesting this evolutionarily conserved regulatory loop is engaged in more physiological activities.

Autophagy is a highly conserved programmed degradation process that regulates a variety of physiological and pathological activities in health, aging, and disease. To identify additional factors that modulate autophagy, we utilized serum-free starvation or Torin1 to induce autophagy in HeLa cells for unbiased mRNA-sequencing analysis and identified SNAI2, a crucial player in epithelial-to-mesenchymal transition and cancer progression, as a regulator of autophagy. Mechanistically, SNAI2 promotes autophagy by physically interacting with FOXO3 and enhancing FOXO3 binding affinity to its response elements in autophagy-related genes. Intriguingly, binding to the DNA targets appears necessary and sufficient for FOXO3 to antagonize its CRM1-dependent nuclear export, illustrating a critical role of DNA in regulating protein nuclear localization. Moreover, stress-elevated SNAI2 expression is mediated by FOXO3, which activates SNAI2 transcription by directly binding to its promoter. Herein, FOXO3 and SNAI2 form a coherent feed-forward regulatory loop to reinforce autophagy genes induction in response to energy stress. Strikingly, a dFoxO-Snail feed-forward circuit also regulates autophagy in Drosophila, suggesting this mechanism is evolutionarily conserved from fly to human.

BMP Signaling: Lighting up the Way for Embryonic Dorsoventral Patterning

One of the most significant events during early embryonic development is the establishment of a basic embryonic body plan, which is defined by anteroposterior, dorsoventral (DV), and left-right axes. It is well-known that the morphogen gradient created by BMP signaling activity is crucial for DV axis patterning across a diverse set of vertebrates. The regulation of BMP signaling during DV patterning has been strongly conserved across evolution. This is a remarkable regulatory and evolutionary feat, as the BMP gradient has been maintained despite the tremendous variation in embryonic size and shape across species. Interestingly, the embryonic DV axis exhibits robust stability, even in face of variations in BMP signaling. Multiple lines of genetic, molecular, and embryological evidence have suggested that numerous BMP signaling components and their attendant regulators act in concert to shape the developing DV axis. In this review, we summarize the current knowledge of the function and regulation of BMP signaling in DV patterning. Throughout, we focus specifically on popular model animals, such as Xenopus and zebrafish, highlighting the similarities and differences of the regulatory networks between species. We also review recent advances regarding the molecular nature of DV patterning, including the initiation of the DV axis, the formation of the BMP gradient, and the regulatory molecular mechanisms behind BMP signaling during the establishment of the DV axis. Collectively, this review will help clarify our current understanding of the molecular nature of DV axis formation.

Generation and timing of graded responses to morphogen gradients

Morphogen gradients are known to subdivide a naive cell field into distinct zones of gene expression. Here, we examine whether morphogens can also induce a graded response within such domains. To this end, we explore the role of the Dorsal protein nuclear gradient along the dorsoventral axis in defining the graded pattern of actomyosin constriction that initiates gastrulation in early Drosophila embryos. Two complementary mechanisms for graded accumulation of mRNAs of crucial zygotic Dorsal target genes were identified. First, activation of target-gene expression expands over time from the ventral-most region of high nuclear Dorsal to lateral regions, where the levels are lower, as a result of a Dorsal-dependent activation probability of transcription sites. Thus, sites that are activated earlier will exhibit more mRNA accumulation. Second, once the sites are activated, the rate of RNA Polymerase II loading is also dependent on Dorsal levels. Morphological restrictions require that translation of the graded mRNA be delayed until completion of embryonic cell formation. Such timing is achieved by large introns, which provide a delay in production of the mature mRNAs. Spatio-temporal regulation of key zygotic genes therefore shapes the pattern of gastrulation.

Chiffon triggers global histone H3 acetylation and expression of developmental genes in Drosophila embryos

The histone acetyltransferase Gcn5 is critical for gene expression and development. In Drosophila, Gcn5 is part of four complexes (SAGA, ATAC, CHAT and ADA) that are essential for fly viability and have key roles in regulating gene expression. Here, we show that although the SAGA, ADA and CHAT complexes play redundant roles in embryonic gene expression, the insect-specific CHAT complex uniquely regulates expression of a subset of developmental genes. We also identify a substantial decrease in histone acetylation in chiffon mutant embryos that exceeds that observed in Ada2b, suggesting broader roles for Chiffon in regulating histone acetylation outside of the Gcn5 complexes. The chiffon gene encodes two independent polypeptides that nucleate formation of either the CHAT or Dbf4-dependent kinase (DDK) complexes. DDK includes the cell cycle kinase Cdc7, which is necessary for maternally driven DNA replication in the embryo. We identify a temporal switch between the expression of these chiffon gene products during a short window during the early nuclear cycles in embryos that correlates with the onset of zygotic genome activation, suggesting a potential role for CHAT in this process.

This article has an associated First Person interview with the first author of the paper.

Summary: Chiffon within the Gcn5-containing CHAT complex plays key roles in embryonic gene expression and histone H3 acetylation in Drosophila.

Fitting thermodynamic-based models: Incorporating parameter sensitivity improves the performance of an evolutionary algorithm

A detailed comprehension of transcriptional regulation is critical to understanding the genetic control of development and disease across many different organisms. To more fully investigate the complex molecular interactions controlling the precise expression of genes, many groups have constructed mathematical models to complement their experimental approaches. A critical step in such studies is choosing the most appropriate parameter estimation algorithm to enable detailed analysis of the parameters that contribute to the models. In this study, we develop a novel set of evolutionary algorithms that use a pseudo-random Sobol Set to construct the initial population and incorporate parameter sensitivities into the adaptation of mutation rates, using local, global, and hybrid strategies. Comparison of the performance of these new algorithms to a number of current state-of-the-art global parameter estimation algorithms on a range of continuous test functions, as well as synthetic biological data representing models of gene regulatory systems, reveals improved performance of the new algorithms in terms of runtime, error and reproducibility. In addition, by analyzing the ability of these algorithms to fit datasets of varying quality, we provide the experimentalist with a guide to how the algorithms perform across a range of noisy data. These results demonstrate the improved performance of the new set of parameter estimation algorithms and facilitate meaningful integration of model parameters and predictions in our understanding of the molecular mechanisms of gene regulation.

Computational biology approaches for mapping transcriptional regulatory networks

Transcriptional Regulatory Networks (TRNs) are mainly responsible for the cell-type- or cell-state-specific expression of gene sets from the same DNA sequence. However, so far there are no precise maps of TRNs available for each cell-type or cell-state, and no ideal tool to map those networks clearly and in full from biological samples. In this review, major approaches and tools to map TRNs from high-throughput data are presented, depending on the type of methods or data used to infer them, and their advantages and limitations are discussed. After summarizing the main principles defining the topology and structure–function relationships in TRNs, an overview of the extensive work done to map TRNs from bulk transcriptomic data will be presented by type of methodological approach. Most recent modellings of TRNs using other types of molecular data or integrating different data types, including single-cell RNA-sequencing and chromatin information, will then be discussed, before briefly concluding with improvements expected to come in the field.

Enhancer redundancy in development and disease

Shadow enhancers are seemingly redundant transcriptional cis-regulatory elements that regulate the same gene and drive overlapping expression patterns. Recent studies have shown that shadow enhancers are remarkably abundant and control most developmental gene expression in both invertebrates and vertebrates, including mammals. Shadow enhancers might provide an important mechanism for buffering gene expression against mutations in non-coding regulatory regions of genes implicated in human disease. Technological advances in genome editing and live imaging have shed light on how shadow enhancers establish precise gene expression patterns and confer phenotypic robustness. Shadow enhancers can interact in complex ways and may also help to drive the formation of transcriptional hubs within the nucleus. Despite their apparent redundancy, the prevalence and evolutionary conservation of shadow enhancers underscore their key role in emerging metazoan gene regulatory networks.

SCRT1 is a novel beta cell transcription factor with insulin regulatory properties

Here we show that scratch family transcriptional repressor 1 (SCRT1), a zinc finger transcriptional regulator, is a novel regulator of beta cell function. SCRT1 was found to be expressed in beta cells in rodent and human islets. In human islets, expression of SCRT1 correlated with insulin secretion capacity and the expression of the insulin (INS) gene. Furthermore, SCRT1 mRNA expression was lower in beta cells from T2D patients. siRNA-mediated Scrt1 silencing in INS-1832/13 cells, mouse- and human islets resulted in impaired glucose-stimulated insulin secretion and decreased expression of the insulin gene. This is most likely due to binding of SCRT1 to E-boxes of the Ins1 gene as shown with ChIP. Scrt1 silencing also reduced the expression of several key beta cell transcription factors. Moreover, Scrt1 mRNA expression was reduced by glucose and SCRT1 protein was found to translocate between the nucleus and the cytosol in a glucose-dependent fashion in INS-1832/13 cells as well as in a rodent model of T2D. SCRT1 was also regulated by a GSK3β-dependent SCRT1-serine phosphorylation. Taken together, SCRT1 is a novel beta cell transcription factor that regulates insulin secretion and is affected in T2D.

Dynamic patterning by morphogens illuminated by cis-regulatory studies

Morphogen concentration changes in space as well as over time during development. However, how these dynamics are interpreted by cells to specify fate is not well understood. Here, we focus on two morphogens: the maternal transcription factors Bicoid and Dorsal, which directly regulate target genes to pattern Drosophila embryos. The actions of these factors at enhancers has been thoroughly dissected and provides a rich platform for understanding direct input by morphogens and their changing roles over time. Importantly, Bicoid and Dorsal do not work alone; we also discuss additional inputs that work with morphogens to control spatiotemporal gene expression in embryos.

Independence of chromatin conformation and gene regulation during Drosophila dorsoventral patterning

The relationship between chromatin organization and gene regulation remains unclear. While disruption of chromatin domains and domain boundaries can lead to misexpression of developmental genes, acute depletion of regulators of genome organization has a relatively small effect on gene expression. It is therefore uncertain whether gene expression and chromatin state drive chromatin organization or whether changes in chromatin organization facilitate cell-type-specific activation of gene expression. Here, using the dorsoventral patterning of the Drosophila melanogaster embryo as a model system, we provide evidence for the independence of chromatin organization and dorsoventral gene expression. We define tissue-specific enhancers and link them to expression patterns using single-cell RNA-seq. Surprisingly, despite tissue-specific chromatin states and gene expression, chromatin organization is largely maintained across tissues. Our results indicate that tissue-specific chromatin conformation is not necessary for tissue-specific gene expression but rather acts as a scaffold facilitating gene expression when enhancers become active.

Emerging role of tumor cell plasticity in modifying therapeutic response

Resistance to cancer therapy is a major barrier to cancer management. Conventional views have proposed that acquisition of resistance may result from genetic mutations. However, accumulating evidence implicates a key role of non-mutational resistance mechanisms underlying drug tolerance, the latter of which is the focus that will be discussed here. Such non-mutational processes are largely driven by tumor cell plasticity, which renders tumor cells insusceptible to the drug-targeted pathway, thereby facilitating the tumor cell survival and growth. The concept of tumor cell plasticity highlights the significance of re-activation of developmental programs that are closely correlated with epithelial-mesenchymal transition, acquisition properties of cancer stem cells, and trans-differentiation potential during drug exposure. From observations in various cancers, this concept provides an opportunity for investigating the nature of anticancer drug resistance. Over the years, our understanding of the emerging role of phenotype switching in modifying therapeutic response has considerably increased. This expanded knowledge of tumor cell plasticity contributes to developing novel therapeutic strategies or combination therapy regimens using available anticancer drugs, which are likely to improve patient outcomes in clinical practice.

Functional Transcription Factor Target Networks Illuminate Control of Epithelial Remodelling

Cell identity is governed by gene expression, regulated by transcription factor (TF) binding at cis-regulatory modules. Decoding the relationship between TF binding patterns and gene regulation is nontrivial, remaining a fundamental limitation in understanding cell decision-making. We developed the NetNC software to predict functionally active regulation of TF targets; demonstrated on nine datasets for the TFs Snail, Twist, and modENCODE Highly Occupied Target (HOT) regions. Snail and Twist are canonical drivers of epithelial to mesenchymal transition (EMT), a cell programme important in development, tumour progression and fibrosis. Predicted "neutral" (non-functional) TF binding always accounted for the majority (50% to 95%) of candidate target genes from statistically significant peaks and HOT regions had higher functional binding than most of the Snail and Twist datasets examined. Our results illuminated conserved gene networks that control epithelial plasticity in development and disease. We identified new gene functions and network modules including crosstalk with notch signalling and regulation of chromatin organisation, evidencing networks that reshape Waddington's epigenetic landscape during epithelial remodelling. Expression of orthologous functional TF targets discriminated breast cancer molecular subtypes and predicted novel tumour biology, with implications for precision medicine. Predicted invasion roles were validated using a tractable cell model, supporting our approach.

Formation, interpretation, and regulation of the Drosophila Dorsal/NF-κB gradient

The morphogen gradient of the transcription factor Dorsal in the early Drosophila embryo has become one of the most widely studied tissue patterning systems. Dorsal is a Drosophila homolog of mammalian NF-κB and patterns the dorsal-ventral axis of the blastoderm embryo into several tissue types by spatially regulating upwards of 100 zygotic genes. Recent studies using fluorescence microscopy and live imaging have quantified the Dorsal gradient and its target genes, which has paved the way for mechanistic modeling of the gradient. In this review, we describe the mechanisms behind the initiation of the Dorsal gradient and its regulation of target genes. The main focus of the review is a discussion of quantitative and computational studies of the Dl gradient system, including regulation of the Dl gradient. We conclude with a discussion of potential future directions.

Exploring a Drosophila Transcription Factor Interaction Network to Identify Cis-Regulatory Modules

Multiple transcription factors (TFs) bind to specific sites in the genome and interact among themselves to form the cis-regulatory modules (CRMs). They are essential in modulating the expression of genes, and it is important to study this interplay to understand gene regulation. In the present study, we integrated experimentally identified TF binding sites collected from published studies with computationally predicted TF binding sites to identify Drosophila CRMs. Along with the detection of the previously known CRMs, this approach identified novel protein combinations. We determined high-occupancy target sites, where a large number of TFs bind. Investigating these sites revealed that Giant, Dichaete, and Knirp are highly enriched in these locations. A common TAG team motif was observed at these sites, which might play a role in recruiting other TFs. While comparing the binding sites at distal and proximal promoters, we found that certain regulatory TFs, such as Zelda, were highly enriched in enhancers. Our study has shown that, from the information available concerning the TF binding sites, the real CRMs could be predicted accurately and efficiently. Although we only may claim co-occurrence of these proteins in this study, it may actually point to their interaction (as known interaction proteins typically co-occur together). Such an integrative approach can, therefore, help us to provide a better understanding of the interplay among the factors, even though further experimental verification is required.

Combinatorial chromatin dynamics foster accurate cardiopharyngeal fate choices

During embryogenesis, chromatin accessibility profiles control lineage-specific gene expression by modulating transcription, thus impacting multipotent progenitor states and subsequent fate choices. Subsets of cardiac and pharyngeal/head muscles share a common origin in the cardiopharyngeal mesoderm, but the chromatin landscapes that govern multipotent progenitors competence and early fate choices remain largely elusive. Here, we leveraged the simplicity of the chordate model Ciona to profile chromatin accessibility through stereotyped transitions from naive Mesp+ mesoderm to distinct fate-restricted heart and pharyngeal muscle precursors. An FGF-Foxf pathway acts in multipotent progenitors to establish cardiopharyngeal-specific patterns of accessibility, which govern later heart vs. pharyngeal muscle-specific expression profiles, demonstrating extensive spatiotemporal decoupling between early cardiopharyngeal enhancer accessibility and late cell-type-specific activity. We found that multiple cis-regulatory elements, with distinct chromatin accessibility profiles and motif compositions, are required to activate Ebf and Tbx1/10, two key determinants of cardiopharyngeal fate choices. We propose that these 'combined enhancers' foster spatially and temporally accurate fate choices, by increasing the repertoire of regulatory inputs that control gene expression, through either accessibility and/or activity.

The genome-wide multi-layered architecture of chromosome pairing in early Drosophila embryos

Genome organization involves cis and trans chromosomal interactions, both implicated in gene regulation, development, and disease. Here, we focus on trans interactions in Drosophila, where homologous chromosomes are paired in somatic cells from embryogenesis through adulthood. We first address long-standing questions regarding the structure of embryonic homolog pairing and, to this end, develop a haplotype-resolved Hi-C approach to minimize homolog misassignment and thus robustly distinguish trans-homolog from cis contacts. This computational approach, which we call Ohm, reveals pairing to be surprisingly structured genome-wide, with trans-homolog domains, compartments, and interaction peaks, many coinciding with analogous cis features. We also find a significant genome-wide correlation between pairing, transcription during zygotic genome activation, and binding of the pioneer factor Zelda. Our findings reveal a complex, highly structured organization underlying homolog pairing, first discovered a century ago in Drosophila. Finally, we demonstrate the versatility of our haplotype-resolved approach by applying it to mammalian embryos.

Functional divergence caused by mutations in an energetic hotspot in ERK2

The most frequent extracellular signal-regulated kinase 2 (ERK2) mutation occurring in cancers is E322K (E-K). ERK2 E-K reverses a buried charge in the ERK2 common docking (CD) site, a region that binds activators, inhibitors, and substrates. Little is known about the cellular consequences associated with this mutation, other than apparent increases in tumor resistance to pathway inhibitors. ERK2 E-K, like the mutation of the preceding aspartate (ERK2 D321N [D-N]) known as the sevenmaker mutation, causes increased activity in cells and evades inactivation by dual-specificity phosphatases. As opposed to findings in cancer cells, in developmental assays in Drosophila, only ERK2 D-N displays a significant gain of function, revealing mutation-specific phenotypes. The crystal structure of ERK2 D-N is indistinguishable from that of wild-type protein, yet this mutant displays increased thermal stability. In contrast, the crystal structure of ERK2 E-K reveals profound structural changes, including disorder in the CD site and exposure of the activation loop phosphorylation sites, which likely account for the decreased thermal stability of the protein. These contiguous mutations in the CD site of ERK2 are both required for docking interactions but lead to unpredictably different functional outcomes. Our results suggest that the CD site is in an energetically strained configuration, and this helps drive conformational changes at distal sites on ERK2 during docking interactions.

ATAC-seq reveals regional differences in enhancer accessibility during the establishment of spatial coordinates in the Drosophila blastoderm

Establishment of spatial coordinates during Drosophila embryogenesis relies on differential regulatory activity of axis patterning enhancers. Concentration gradients of activator and repressor transcription factors (TFs) provide positional information to each enhancer, which in turn promotes transcription of a target gene in a specific spatial pattern. However, the interplay between an enhancer regulatory activity and its accessibility as determined by local chromatin organization is not well understood. We profiled chromatin accessibility with ATAC-seq in narrow, genetically tagged domains along the antero-posterior axis in the Drosophila blastoderm. We demonstrate that one-quarter of the accessible genome displays significant regional variation in its ATAC-seq signal immediately after zygotic genome activation. Axis patterning enhancers are enriched among the most variable intervals, and their accessibility changes correlate with their regulatory activity. In an embryonic domain where an enhancer receives a net activating TF input and promotes transcription, it displays elevated accessibility in comparison to a domain where it receives a net repressive input. We propose that differential accessibility is a signature of patterning cis-regulatory elements in the Drosophila blastoderm and discuss potential mechanisms by which accessibility of enhancers may be modulated by activator and repressor TFs.

Differential transcriptional regulation by alternatively designed mechanisms: A mathematical modeling approach

Cells maintain cellular homeostasis employing different regulatory mechanisms to respond external stimuli. We study two groups of signal-dependent transcriptional regulatory mechanisms. In the first group, we assume that repressor and activator proteins compete for binding to the same regulatory site on DNA (competitive mechanisms). In the second group, they can bind to different regulatory regions in a noncompetitive fashion (noncompetitive mechanisms). For both competitive and noncompetitive mechanisms, we studied the gene expression dynamics by increasing the repressor or decreasing the activator abundance (inhibition mechanisms), or by decreasing the repressor or increasing the activator abundance (activation mechanisms). We employed delay differential equation models. Our simulation results show that the competitive and noncompetitive inhibition mechanisms exhibit comparable repression effectiveness. However, response time is fastest in the noncompetitive inhibition mechanism due to increased repressor abundance, and slowest in the competitive inhibition mechanism by increased repressor level. The competitive and noncompetitive inhibition mechanisms through decreased activator abundance show comparable and moderate response times, while the competitive and noncompetitive activation mechanisms by increased activator protein level display more effective and faster response. Our study exemplifies the importance of mathematical modeling and computer simulation in the analysis of gene expression dynamics.

Redundant and Cryptic Enhancer Activities of the Drosophila yellow Gene

Cis-regulatory sequences known as enhancers play a key role in regulating gene expression. Evolutionary changes in these DNA sequences contribute to phenotypic evolution. The Drosophila yellow gene, which is required for pigmentation, has emerged as a model system for understanding how cis-regulatory sequences evolve, providing some of the most detailed insights available into how activities of orthologous enhancers have diverged between species. Here, we examine the evolution of yellow cis-regulatory sequences on a broader scale, by comparing the distribution and function of yellow enhancer activities throughout the 5' intergenic and intronic sequences of Drosophila melanogaster, D. pseudoobscura, and D. willistoni We find that cis-regulatory sequences driving expression in a particular tissue are not as modular as previously described, but rather have many redundant and cryptic enhancer activities distributed throughout the regions surveyed. Interestingly, cryptic enhancer activities of sequences from one species often drove patterns of expression observed in other species, suggesting that the frequent evolutionary changes in yellow expression observed among Drosophila species may be facilitated by gaining and losing repression of preexisting cis-regulatory sequences.

Signal Integration by Shadow Enhancers and Enhancer Duplications Varies across the Drosophila Embryo

Transcription of developmental genes is controlled by multiple enhancers. Frequently, more than one enhancer can activate transcription from the same promoter in the same cells. How is regulatory information from multiple enhancers combined to determine the overall expression output? We measure nascent transcription driven by a pair of shadow enhancers, each enhancer of the pair separately, and each duplicated, using live imaging in Drosophila embryos. This set of constructs allows us to quantify the input-output function describing signal integration by two enhancers. We show that signal integration performed by these shadow enhancers and duplications varies across the expression pattern, implying that how their activities are combined depends on the transcriptional regulators bound to the enhancers in different parts of the embryo. Characterizing signal integration by multiple enhancers is a critical step in developing conceptual and computational models of gene expression at the locus level, where multiple enhancers control transcription together.

Inference of Transcription Factor Regulation Patterns Using Gene Expression Covariation in Natural Populations of Drosophila melanogaster

Gene regulatory networks control the complex programs that drive development. Deciphering the connections between transcription factors (TFs) and target genes is challenging, in part because TFs bind to thousands of places in the genome but control expression through a subset of these binding events. We hypothesize that we can combine natural variation of expression levels and predictions of TF binding sites to identify TF targets. We gather RNA-seq data from 71 genetically distinct F1 Drosophila melanogaster embryos and calculate the correlations between TF and potential target genes' expression levels, which we call "regulatory strength." To separate direct and indirect TF targets, we hypothesize that direct TF targets will have a preponderance of binding sites in their upstream regions. Using 14 TFs active during embryogenesis, we find that 12 TFs showed a significant correlation between their binding strength and regulatory strength on downstream targets, and 10 TFs showed a significant correlation between the number of binding sites and the regulatory effect on target genes. The general roles, e.g. bicoid's role as an activator, and the particular interactions we observed between our TFs, e.g. twist's role as a repressor of sloppy paired and odd paired, generally coincide with the literature.

Insect Transcription Factors: A Landscape of Their Structures and Biological Functions in Drosophila and beyond

Transcription factors (TFs) play essential roles in the transcriptional regulation of functional genes, and are involved in diverse physiological processes in living organisms. The fruit fly Drosophila melanogaster, a simple and easily manipulated organismal model, has been extensively applied to study the biological functions of TFs and their related transcriptional regulation mechanisms. It is noteworthy that with the development of genetic tools such as CRISPR/Cas9 and the next-generation genome sequencing techniques in recent years, identification and dissection the complex genetic regulatory networks of TFs have also made great progress in other insects beyond Drosophila. However, unfortunately, there is no comprehensive review that systematically summarizes the structures and biological functions of TFs in both model and non-model insects. Here, we spend extensive effort in collecting vast related studies, and attempt to provide an impartial overview of the progress of the structure and biological functions of current documented TFs in insects, as well as the classical and emerging research methods for studying their regulatory functions. Consequently, considering the importance of versatile TFs in orchestrating diverse insect physiological processes, this review will assist a growing number of entomologists to interrogate this understudied field, and to propel the progress of their contributions to pest control and even human health.

A feed-forward relay integrates the regulatory activities of Bicoid and Orthodenticle via sequential binding to suboptimal sites

Datta et al. define three major classes of enhancers that are differentially sensitive to binding and transcriptional activation by Bicoid (Bcd) and Orthodenticle (Otd). The specific activities of enhancers in each class are mediated by DNA motif variants preferentially bound by Bcd or Otd and the presence or absence of sites for cofactors that interact with these proteins.

The K50 (lysine at amino acid position 50) homeodomain (HD) protein Orthodenticle (Otd) is critical for anterior patterning and brain and eye development in most metazoans. In Drosophila melanogaster, another K50HD protein, Bicoid (Bcd), has evolved to replace Otd's ancestral function in embryo patterning. Bcd is distributed as a long-range maternal gradient and activates transcription of a large number of target genes, including otd. Otd and Bcd bind similar DNA sequences in vitro, but how their transcriptional activities are integrated to pattern anterior regions of the embryo is unknown. Here we define three major classes of enhancers that are differentially sensitive to binding and transcriptional activation by Bcd and Otd. Class 1 enhancers are initially activated by Bcd, and activation is transferred to Otd via a feed-forward relay (FFR) that involves sequential binding of the two proteins to the same DNA motif. Class 2 enhancers are activated by Bcd and maintained by an Otd-independent mechanism. Class 3 enhancers are never bound by Bcd, but Otd binds and activates them in a second wave of zygotic transcription. The specific activities of enhancers in each class are mediated by DNA motif variants preferentially bound by Bcd or Otd and the presence or absence of sites for cofactors that interact with these proteins. Our results define specific patterning roles for Bcd and Otd and provide mechanisms for coordinating the precise timing of gene expression patterns during embryonic development.

Twist promotes tumor metastasis in basal-like breast cancer by transcriptionally upregulating ROR1

Rationale: Twist is a key transcription factor for induction of epithelial-mesenchymal transition (EMT), which promotes cell migration, invasion, and cancer metastasis, confers cancer cells with stem cell-like characteristics, and provides therapeutic resistance. However, the functional roles and targeted genes of Twist in EMT and cancer progression remain elusive. Methods: The potential targeted genes of Twist were identified from the global transcriptomes of T47D/Twist cells by microarray analysis. EMT phenotype was detected by western blotting and immunofluorescence of marker proteins. The dual-luciferase reporter and chromatin immunoprecipitation assays were employed to observe the direct transcriptional induction of ROR1 by Twist. A lung metastasis model was used to study the pro-metastatic role of Twist and ROR1 by injecting MDA-MB-231 cells into tail vein of nude mice. Bio-informatics analysis was utilized to measure the metastasis-free survival of breast cancer patients. Results: Twist protein was proved to directly activate the transcription of ROR1 gene, a receptor of Wnt5a in non-canonical WNT signaling pathway. Silencing of ROR1 inhibited EMT process, cell migration, invasion, and cancer metastasis of basal-like breast cancer (BLBC) cells. Knockdown of ROR1 also ameliorated the pro-metastatic effect of Twist. Furthermore, analyses of clinical specimens indicated that high expression of both ROR1 and Twist tightly correlates with poor metastasis-free survival of breast cancer patients. Conclusion: ROR1 is a targeted gene of Twist. Twist/ROR1 signaling is critical for invasion and metastasis of BLBC cells.

Enhancer talk

Enhancers are short noncoding segments of DNA (100-1000 bp) that control the temporal and spatial activity of genes in an orientation-independent manner. They can be separated from their target genes by large distances and are thus known as distal regulatory elements. One consequence of the variability in the distance separating enhancers and their target promoters is that it is difficult to determine which elements are involved in the regulation of a particular gene. Moreover, enhancers can be found in clusters in which multiple regulatory elements control expression of the same target gene. However, little is known about how the individual elements contribute to gene expression. Here, we describe how chromatin conformation promotes and constraints enhancer activity. Further, we discuss enhancer clusters and what is known about the contribution of individual elements to the regulation of target genes. Finally, we examine the reliability of different methods used to identify enhancers.

Nutrient-Driven O-GlcNAcylation at Promoters Impacts Genome-Wide RNA Pol II Distribution

Nutrient-driven O-GlcNAcylation has been linked to epigenetic regulation of gene expression in metazoans. In C. elegans, O-GlcNAc marks the promoters of over 800 developmental, metabolic, and stress-related genes; these O-GlcNAc marked genes show a strong 5', promoter-proximal bias in the distribution of RNA Polymerase II (Pol II). In response to starvation or feeding, the steady state distribution of O-GlcNAc at promoters remain nearly constant presumably due to dynamic cycling mediated by the transferase OGT-1 and the O-GlcNAcase OGA-1. However, in viable mutants lacking either of these enzymes of O-GlcNAc metabolism, the nutrient-responsive GlcNAcylation of promoters is dramatically altered. Blocked O-GlcNAc cycling leads to a striking nutrient-dependent accumulation of O-GlcNAc on RNA Pol II. O-GlcNAc cycling mutants also show an exaggerated, nutrient-responsive redistribution of promoter-proximal RNA Pol II isoforms and extensive transcriptional deregulation. Our findings suggest a complex interplay between the O-GlcNAc modification at promoters, the kinase-dependent "CTD-code," and co-factors regulating RNA Pol II dynamics. Nutrient-responsive O-GlcNAc cycling may buffer the transcriptional apparatus from dramatic swings in nutrient availability by modulating promoter activity to meet metabolic and developmental needs.

Enhancer adoption caused by genomic insertion elicits interdigital Shh expression and syndactyly in mouse

In this study, we reexamined an old mouse mutant named Hammer toe (Hm), which arose spontaneously almost a half century ago and exhibits a limb phenotype with webbing. We revealed that a 150-kb noncoding genomic fragment that was originally located in chromosome 14 has been inserted into a genomic region proximal to Sonic hedgehog (Shh), located in chromosome 5. This inserted fragment possesses enhancer activity to induce Shh expression in the interdigital regions in Hm, which in turn down-regulates bone morphogenetic protein signaling and eventually results in syndactyly and web formation. Since the donor fragment residing in chromosome 14 has enhancer activity to induce interdigital gene expression, the Hm mutation appears to be an archetypal case of enhancer adoption.

Acquisition of new cis-regulatory elements (CREs) can cause alteration of developmental gene regulation and may introduce morphological novelty in evolution. Although structural variation in the genome generated by chromosomal rearrangement is one possible source of new CREs, only a few examples are known, except for cases of retrotransposition. In this study, we show the acquisition of novel regulatory sequences as a result of large genomic insertion in the spontaneous mouse mutation Hammer toe (Hm). Hm mice exhibit syndactyly with webbing, due to suppression of interdigital cell death in limb development. We reveal that, in the Hm genome, a 150-kb noncoding DNA fragment from chromosome 14 is inserted into the region upstream of the Sonic hedgehog (Shh) promoter in chromosome 5. Phenotyping of mouse embryos with a series of CRISPR/Cas9-aided partial deletion of the 150-kb insert clearly indicated that two different regions are necessary for the syndactyly phenotype of Hm. We found that each of the two regions contains at least one enhancer for interdigital regulation. These results show that a set of enhancers brought by the large genomic insertion elicits the interdigital Shh expression and the Hm phenotype. Transcriptome analysis indicates that ectopic expression of Shh up-regulates Chordin (Chrd) that antagonizes bone morphogenetic protein signaling in the interdigital region. Indeed, Chrd-overexpressing transgenic mice recapitulated syndactyly with webbing. Thus, the Hm mutation provides an insight into enhancer acquisition as a source of creation of novel gene regulation.

Caenorhabditis elegans CES-1 Snail Represses pig-1 MELK Expression To Control Asymmetric Cell Division

Snail-like transcription factors affect stem cell function through mechanisms that are incompletely understood. In the Caenorhabditis elegans neurosecretory motor neuron (NSM) neuroblast lineage, CES-1 Snail coordinates cell cycle progression and cell polarity to ensure the asymmetric division of the NSM neuroblast and the generation of two daughter cells of different sizes and fates. We have previously shown that CES-1 Snail controls cell cycle progression by repressing the expression of cdc-25.2 CDC25. However, the mechanism through which CES-1 Snail affects cell polarity has been elusive. Here, we systematically searched for direct targets of CES-1 Snail by genome-wide profiling of CES-1 Snail binding sites and identified >3000 potential CES-1 Snail target genes, including pig-1, the ortholog of the oncogene maternal embryonic leucine zipper kinase (MELK). Furthermore, we show that CES-1 Snail represses pig-1 MELK transcription in the NSM neuroblast lineage and that pig-1 MELK acts downstream of ces-1 Snail to cause the NSM neuroblast to divide asymmetrically by size and along the correct cell division axis. Based on our results we propose that by regulating the expression of the MELK gene, Snail-like transcription factors affect the ability of stem cells to divide asymmetrically and, hence, to self-renew. Furthermore, we speculate that the deregulation of MELK contributes to tumorigenesis by causing cells that normally divide asymmetrically to divide symmetrically instead.

Modeling the cis-regulatory modules of genes expressed in developmental stages of Drosophila melanogaster

Because transcription is the first step in the regulation of gene expression, understanding how transcription factors bind to their DNA binding motifs has become absolutely necessary. It has been shown that the promoters of genes with similar expression profiles share common structural patterns. This paper presents an extensive study of the regulatory regions of genes expressed in 24 developmental stages of Drosophila melanogaster. It proposes the use of a combination of structural features, such as positioning of individual motifs relative to the transcription start site, orientation, pairwise distance between motifs, and presence of motifs anywhere in the promoter for predicting gene expression from structural features of promoter sequences. RNA-sequencing data was utilized to create and validate the 24 models. When genes with high-scoring promoters were compared to those identified by RNA-seq samples, 19 (79.2%) statistically significant models, a number that exceeds previous studies, were obtained. Each model yielded a set of highly informative features, which were used to search for genes with similar biological functions.

Molecular control of gut formation in the spider Parasteatoda tepidariorum

The development of a digestive system is an essential feature of bilaterians. Studies of the molecular control of gut formation in arthropods have been studied in detail in the fruit fly Drosophila melanogaster. However, little is known in other arthropods, especially in noninsect arthropods. To better understand the evolution of arthropod alimentary system, we investigate the molecular control of gut development in the spider Parasteatoda tepidariorum (Pt), the primary chelicerate model species for developmental studies. Orthologs of the ectodermal genes Pt-wingless (Pt-wg) and Pt-hedgehog (Pt-hh), of the endodermal genes, Pt-serpent (Pt-srp) and Pt-hepatocyte-nuclear factor-4 (Pt-hnf4) and of the mesodermal gene Pt-twist (Pt-twi) are expressed in the same germ layers during spider gut development as in D. melanogaster. Thus, our expression data suggest that the downstream molecular components involved in gut development in arthropods are conserved. However, Pt-forkhead (Pt-fkh) expression and function in spiders is considerably different from its D. melanogaster ortholog. Pt-fkh is expressed before gastrulation in a cell population that gives rise to endodermal and mesodermal precursors, suggesting a possible role for this factor in specification of both germ layers. To test this hypothesis, we knocked down Pt-fkh via RNA interference. Pt-fkh RNAi embryos not only fail to develop a proper gut, but also lack the mesodermal Pt-twi expressing cells. Thus, in spiders Pt-fkh specifies endodermal and mesodermal germ layers. We discuss the implications of these findings for the evolution and development of gut formation in Ecdysozoans.

Mechanisms of Groucho-mediated repression revealed by genome-wide analysis of Groucho binding and activity

Background

The transcriptional corepressor Groucho (Gro) is required for the function of many developmentally regulated DNA binding repressors, thus helping to define the gene expression profile of each cell during development. The ability of Gro to repress transcription at a distance together with its ability to oligomerize and bind to histones has led to the suggestion that Gro may spread along chromatin. However, much is unknown about the mechanism of Gro-mediated repression and about the dynamics of Gro targeting.

Results

Our chromatin immunoprecipitation sequencing analysis of temporally staged Drosophila embryos shows that Gro binds in a highly dynamic manner primarily to clusters of discrete (<1 kb) segments. Consistent with the idea that Gro may facilitate communication between silencers and promoters, Gro binding is enriched at both cis-regulatory modules, as well as within the promotors of potential target genes. While this Gro-recruitment is required for repression, our data show that it is not sufficient for repression. Integration of Gro binding data with transcriptomic analysis suggests that, contrary to what has been observed for another Gro family member, Drosophila Gro is probably a dedicated repressor. This analysis also allows us to define a set of high confidence Gro repression targets. Using publically available data regarding the physical and genetic interactions between these targets, we are able to place them in the regulatory network controlling development. Through analysis of chromatin associated pre-mRNA levels at these targets, we find that genes regulated by Gro in the embryo are enriched for characteristics of promoter proximal paused RNA polymerase II.

Conclusions

Our findings are inconsistent with a one-dimensional spreading model for long-range repression and suggest that Gro-mediated repression must be regulated at a post-recruitment step. They also show that Gro is likely a dedicated repressor that sits at a prominent highly interconnected regulatory hub in the developmental network. Furthermore, our findings suggest a role for RNA polymerase II pausing in Gro-mediated repression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3589-6) contains supplementary material, which is available to authorized users.

PLCE1 Promotes Esophageal Cancer Cell Progression by Maintaining the Transcriptional Activity of Snail

Esophageal cancer is among the most deadly malignant diseases. However, the genetic factors contributing to its occurrence are poorly understood. Multiple studies with large clinic-based cohorts revealed that variations of the phospholipase C epsilon (PLCE1) gene were associated with esophageal cancer susceptibility. However, the causative role of PLCE1 in esophageal cancer is not clear. We inactivated the functional alleles of PLCE1 by CRISPR/Cas9 genome editing technology. The resultant PLCE1 inactivated cells were analyzed both in vitro and in vivo. Our results showed that loss of PLCE1 dramatically decreased the invasion and proliferation capacity of esophageal carcinoma cells in vitro. Moreover, such PLCE1 inactivated tumor grafts exhibited significantly decreased tumor size in mice. We found that PLCE1 was required to maintain protein level of snail a key transcription factor responsible for invasion. Our further transcriptomic data revealed that deficient cells were significantly decreased in expression of genes enriched as targets of Snail. Strikingly, recovery of Snail protein at least partially rescued the invasion and proliferation capacity in PLCE1 inactivated cells. In ESCC clinical specimens, PLCE1 was correlated with tumor stage (P<.0001). Interestingly, PLCE1 expression was positively correlated Snail by immunohistochemistry in such specimens (P<.0001). Therefore, our functional experiments showed the essential roles of PLCE1 in esophageal carcinoma cells and provided evidences that targeting PLCE1 and its downstream molecules could be effective therapies for esophageal cancer.

Transcriptional activity of the short gastrulation primary enhancer in the ventral midline requires its early activity in the presumptive neurogenic ectoderm

The short gastrulation (sog) shadow enhancer directs early and late sog expression in the neurogenic ectoderm and the ventral midline of the developing Drosophila embryo, respectively. Here, evidence is presented that the sog primary enhancer also has both activities, with the late enhancer activity dependent on the early activity. Computational analyses showed that the sog primary enhancer contains five Dorsal (Dl)-, four Zelda (Zld)-, three Bicoid (Bcd)-, and no Single-minded (Sim)-binding sites. In contrast to many ventral midline enhancers, the primary enhancer can direct lacZ expression in the ventral midline as well as in the neurogenic ectoderm without a canonical Simbinding site. Intriguingly, the impaired transcriptional synergy between Dl and either Zld or Bcd led to aberrant and abolished lacZ expression in the neurogenic ectoderm and in the ventral midline, respectively. These findings suggest that the two enhancer activities of the sog primary enhancer are functionally consolidated and geographically inseparable. [BMB Reports 2016; 49(10): 572-577]