Patterns of gene expression during Drosophila mesoderm development

Whole-Mount RNA In Situ Hybridization of Zebrafish Embryos

A commonly employed technique in molecular biology to evaluate the temporal and spatial expression of a certain gene is in situ hybridization. This method is an effective strategy to construct synexpression groups, co-expressed genes acting in shared biological processes, and to find new members of genes engaged in the same signaling pathways to discover similar spatial and temporal expression patterns in zebrafish embryos. The major disadvantage of this method is that RNA probes can penetrate within 2 days of post-fertilization embryos, and therefore, in later developmental stages, the probe can only reach the surface tissues. Further application of the method in histological sections will be required for a complete and accurate gene expression investigation. However, this method is highly effective at late embryogenesis and early larval stages for observing gene expression in endodermal derivatives and sensory organs. RNA probes for in situ hybridization can be prepared through in vitro transcription from plasmids carrying specific promoter elements and mRNA-specific cDNA, or an alternative polymerase chain reaction (PCR) method can be used through PCR amplification. This chapter describes the procedures for detecting gene expression in zebrafish embryos using whole-mount RNA in situ hybridization.

Microarray Analysis to Determine Gene Expression Changes in Zebrafish Embryos

Chemical exposure in humans begins from the zygote stage and continues throughout the development of the embryo and the fetus. Zebrafish are one of the most powerful model organisms used in many research areas, including genetics, environmental toxicology, development, DNA damage and repair, cancer, and other diseases. Among the advantages that facilitate the use of zebrafish as a model for studies are features such as high homology with the human genome, small size, and high reproductive potential in short periods. The use of zebrafish embryos in research has increased rapidly due to their advantageous properties, including extrauterine development and the transparent feature of the embryos. However, there are thousands of genes that can be encountered in research, and in this case, the workforce is too much. This workload has been alleviated with the developed technologies. Microarray is one of these technologies. An important parameter in this assay is the RIN value. The RIN value ranges from 1 to 10, indicating mRNA degradation, and therefore helps to decide whether to continue the study. In this chapter, microarray analysis, which is one of the main techniques used in the determination of gene expression in zebrafish embryos, is described.

Dual mechanism of Let-7i in tumor progression

Let-7i regulates tumors primarily by binding to the 3' untranslated region (3' UTR) of mRNA, which indirectly regulates post-transcriptional gene expression. Let-7i also has an epigenetic function via modulating DNA methylation to directly regulate gene expression. Let-7i performs a dual role by inducing both the promotion and inhibition of various malignancies, depending on its target. The mechanism of Let-7i action involves cancer cell proliferation, migration, invasion, apoptosis, epithelial-mesenchymal transition, EV transmission, angiogenesis, autophagy, and drug resistance sensitization. Let-7i is closely related to cancer, and hence, is a potential biomarker for the diagnosis and prognosis of various cancers. Therapeutically, it can be used to promote an anti-cancer immune response by modifying exosomes, thus exerting a tumor-suppressive effect.

Quantitative proteomics and phosphoproteomics reveal insights into mechanisms of ocnus function in Drosophila testis development

BACKGROUND

Testis is the only organ supporting sperm production and with the largest number of proteins and tissue-specific proteins in animals. In our previous studies, we have found that knockdown of ocnus (ocn), a testis-specific gene, resulted in much smaller testis with no germ cells in Drosophila melanogaster. However, the molecular consequences of ocn knockdown in fly testes are unknown.

RESULTS

In this study, through iTRAQ quantitative proteomics sequencing, 606 proteins were identified from fly abdomens as having a significant and at least a 1.5-fold change in expression after ocn knockdown in fly testes, of which 85 were up-regulated and 521 were down-regulated. Among the differential expressed proteins (DEPs), apart from those proteins involved in spermatogenesis, the others extensively affected biological processes of generation of precursor metabolites and energy, metabolic process, and mitochondrial transport. Protein-protein interaction (PPI) analyses of DEPs showed that several kinases and/or phosphatases interacted with Ocn. Re-analyses of the transcriptome revealed 150 differential expressed genes (DEGs) appeared in the DEPs, and their changing trends in expressions after ocn knockdown were consistent. Many common down-regulated DEGs and DEPs were testis-specific or highly expressed in the testis of D. melanogaster. Quantitative RT-PCR (qRT-PCR) confirmed 12 genes appeared in both DEGs and DEPs were significantly down-regulated after ocn knockdown in fly testes. Furthermore, 153 differentially expressed phosphoproteins (DEPPs), including 72 up-regulated and 94 down-regulated phosphorylated proteins were also identified (13 phosphoproteins appeared in both up- and down-regulated groups due to having multiple phosphorylation sites). In addition to those DEPPs associated with spermatogenesis, the other DEPPs were enriched in actin filament-based process, protein folding, and mesoderm development. Some DEPs and DEPPs were involved in Notch, JAK/STAT, and cell death pathways.

CONCLUSIONS

Given the drastic effect of the ocn knockdown on tissue development and testis cells composition, the differences in protein abundance in the ocn knockdown flies might not necessarily be the direct result of differential gene regulation due to the inactivation of ocn. Nevertheless, our results suggest that the expression of ocn is essential for Drosophila testis development and that its down-regulation disturbs key signaling pathways related to cell survival and differentiation. These DEPs and DEPPs identified may provide significant candidate set for future studies on the mechanism of male reproduction of animals, including humans.

Actomyosin pulsing rescues embryonic tissue folding from disruption by myosin fluctuations

During early development, myosin II mechanically reshapes and folds embryo tissue. A much-studied example is ventral furrow formation in Drosophila , marking the onset of gastrulation. Furrowing is driven by contraction of actomyosin networks on apical cell surfaces, but how the myosin patterning encodes tissue shape is unclear, and elastic models failed to reproduce essential features of experimental cell contraction profiles. The myosin patterning exhibits substantial cell-to-cell fluctuations with pulsatile time-dependence, a striking but unexplained feature of morphogenesis in many organisms. Here, using biophysical modeling we find viscous forces offer the principle resistance to actomyosin-driven apical constriction. In consequence, tissue shape is encoded in the direction-dependent curvature of the myosin patterning which orients an anterior-posterior furrow. Tissue contraction is highly sensitive to cell-to-cell myosin fluctuations, explaining furrowing failure in genetically perturbed embryos whose fluctuations are temporally persistent. In wild-type embryos, this catastrophic outcome is averted by pulsatile myosin time-dependence, a time-averaging effect that rescues furrowing. This low pass filter mechanism may underlie the usage of actomyosin pulsing in diverse morphogenetic processes across many organisms.

Patched-Related Is Required for Proper Development of Embryonic Drosophila Nervous System

Patched-related (Ptr), classified primarily as a neuroectodermal gene, encodes a protein with predicted topology and domain organization closely related to those of Patched (Ptc), the canonical receptor of the Hedgehog (Hh) pathway. To investigate the physiological function of Ptr in the developing nervous system, Ptr null mutant embryos were immunolabeled and imaged under confocal microscopy. These embryos displayed severe alterations in the morphology of the primary axonal tracts, reduced number, and altered distribution of the Repo-positive glia as well as peripheral nervous system defects. Most of these alterations were recapitulated by downregulating Ptr expression, specifically in embryonic nerve cells. Because similar nervous system phenotypes have been observed in hh and ptc mutant embryos, we evaluated the Ptr participation in the Hh pathway by performing cell-based reporter assays. Clone-8 cells were transfected with Ptr-specific dsRNA or a Ptr DNA construct and assayed for changes in Hh-mediated induction of a luciferase reporter. The results obtained suggest that Ptr could act as a negative regulator of Hh signaling. Furthermore, co-immunoprecipitation assays from cell culture extracts premixed with a conditioned medium revealed a direct interaction between Ptr and Hh. Moreover, in vivo Ptr overexpression in the domain of the imaginal wing disc where Engrailed and Ptc coexist produced wing phenotypes at the A/P border. Thus, these results strongly suggest that Ptr plays a crucial role in nervous system development and appears to be a negative regulator of the Hh pathway.

A mosaic of conserved and novel modes of gene expression and morphogenesis in mesoderm and muscle formation of a larval bivalve

The mesoderm gives rise to several key morphological features of bilaterian animals including endoskeletal elements and the musculature. A number of regulatory genes involved in mesoderm and/or muscle formation (e.g., Brachyury (Bra), even-skipped (eve), Mox, myosin II heavy chain (mhc)) have been identified chiefly from chordates and the ecdysozoans Drosophila and Caenorhabditis elegans, but data for non-model protostomes, especially those belonging to the ecdysozoan sister clade, Lophotrochozoa (e.g., flatworms, annelids, mollusks), are only beginning to emerge. Within the lophotrochozoans, Mollusca constitutes the most speciose and diverse phylum. Interestingly, however, information on the morphological and molecular underpinnings of key ontogenetic processes such as mesoderm formation and myogenesis remains scarce even for prominent molluscan sublineages such as the bivalves. Here, we investigated myogenesis and developmental expression of Bra, eve, Mox, and mhc in the quagga mussel Dreissena rostriformis, an invasive freshwater bivalve and an emerging model in invertebrate evodevo. We found that all four genes are expressed during mesoderm formation, but some show additional, individual sites of expression during ontogeny. While Mox and mhc are involved in early myogenesis, eve is also expressed in the embryonic shell field and Bra is additionally present in the foregut. Comparative analysis suggests that Mox has an ancestral role in mesoderm and possibly muscle formation in bilaterians, while Bra and eve are conserved regulators of mesoderm development of nephrozoans (protostomes and deuterostomes). The fully developed Dreissena veliger larva shows a highly complex muscular architecture, supporting a muscular ground pattern of autobranch bivalve larvae that includes at least a velum muscle ring, three or four pairs of velum retractors, one or two pairs of larval retractors, two pairs of foot retractors, a pedal plexus, possibly two pairs of mantle retractors, and the muscles of the pallial line, as well as an anterior and a posterior adductor. As is typical for their molluscan kin, remodelling and loss of prominent larval features such as the velum musculature and various retractor systems appear to be also common in bivalves.

Supplementary information

The online version contains supplementary material available at 10.1007/s13127-022-00569-5.

Bisdemethoxycurcumin Promotes Apoptosis and Inhibits the Epithelial-Mesenchymal Transition through the Inhibition of the G-Protein-Coupled Receptor 161/Mammalian Target of Rapamycin Signaling Pathway in Triple Negative Breast Cancer Cells

Triple negative breast cancer (TNBC) is one of the leading causes of cancer death in the world and lacks an effective targeted therapy. G-protein-coupled receptor 161 (GPR161) has been demonstrated to perform the functional regulations on TNBC progression and might be a potential new target for TNBC therapy. This study showed the effects of bisdemethoxycurcumin (BDMC) on GPR161 regulation, indicating that BDMC effectively inhibited GPR161 expression and downregulated GPR161-driven signaling. BDMC showed the potent inhibitory effects on TNBC proliferation through suppressing GPR161-mediated mammalian target of rapamycin (mTOR)/70 kDa ribosomal protein S6 kinase (p70S6K) activation. Besides, in this study, we discover the mechanism of GPR161-driven TNBC metastasis, linking to GPR161-mediated twist-related protein 1 (Twist1)/matrix metallopeptidase 9 (MMP9) contributing to the epithelial-mesenchymal transition (EMT). BDMC effectively repressed GPR161-mediated TNBC metastasis via inhibiting Twist1/MMP9-induced EMT. The three-dimensional invasion assay also showed that BDMC significantly inhibited TNBC invasion. The combination treatment of BDMC and rapamycin enhanced the inhibition of TNBC proliferation and metastasis through increasing the blockage of mTOR activation. Furthermore, this study also observed that BDMC activated the caspase 3/9 signaling pathway to induce TNBC apoptosis. Therefore, BDMC could be applicable to anticancer therapy, especially targeting on the GPR161-driven cancer type.

Extremely rapid and reversible optogenetic perturbation of nuclear proteins in living embryos

Many developmental regulators have complex and context-specific roles in different tissues and stages, making the dissection of their function extremely challenging. As regulatory processes often occur within minutes, perturbation methods that match these dynamics are needed. Here, we present the improved light-inducible nuclear export system (iLEXY), an optogenetic loss-of-function approach that triggers translocation of proteins from the nucleus to the cytoplasm. By introducing a series of mutations, we substantially increased LEXY's efficiency and generated variants with different recovery times. iLEXY enables rapid (t1/2 < 30 s), efficient, and reversible nuclear protein depletion in embryos, and is generalizable to proteins of diverse sizes and functions. Applying iLEXY to the Drosophila master regulator Twist, we phenocopy loss-of-function mutants, precisely map the Twist-sensitive embryonic stages, and investigate the effects of timed Twist depletions. Our results demonstrate the power of iLEXY to dissect the function of pleiotropic factors during embryogenesis with unprecedented temporal precision.

Combinatorial patterns of graded RhoA activation and uniform F-actin depletion promote tissue curvature

During development, gene expression regulates cell mechanics and shape to sculpt tissues. Epithelial folding proceeds through distinct cell shape changes that occur simultaneously in different regions of a tissue. Here, using quantitative imaging in Drosophila melanogaster, we investigate how patterned cell shape changes promote tissue bending during early embryogenesis. We find that the transcription factors Twist and Snail combinatorially regulate a multicellular pattern of lateral F-actin density that differs from the previously described Myosin-2 gradient. This F-actin pattern correlates with whether cells apically constrict, stretch or maintain their shape. We show that the Myosin-2 gradient and F-actin depletion do not depend on force transmission, suggesting that transcriptional activity is required to create these patterns. The Myosin-2 gradient width results from a gradient in RhoA activation that is refined through the balance between RhoGEF2 and the RhoGAP C-GAP. Our experimental results and simulations of a 3D elastic shell model show that tuning gradient width regulates tissue curvature.

Identification of genomic enhancers through spatial integration of single-cell transcriptomics and epigenomics

Single-cell technologies allow measuring chromatin accessibility and gene expression in each cell, but jointly utilizing both layers to map bona fide gene regulatory networks and enhancers remains challenging. Here, we generate independent single-cell RNA-seq and single-cell ATAC-seq atlases of the Drosophila eye-antennal disc and spatially integrate the data into a virtual latent space that mimics the organization of the 2D tissue using ScoMAP (Single-Cell Omics Mapping into spatial Axes using Pseudotime ordering). To validate spatially predicted enhancers, we use a large collection of enhancer-reporter lines and identify ~ 85% of enhancers in which chromatin accessibility and enhancer activity are coupled. Next, we infer enhancer-to-gene relationships in the virtual space, finding that genes are mostly regulated by multiple, often redundant, enhancers. Exploiting cell type-specific enhancers, we deconvolute cell type-specific effects of bulk-derived chromatin accessibility QTLs. Finally, we discover that Prospero drives neuronal differentiation through the binding of a GGG motif. In summary, we provide a comprehensive spatial characterization of gene regulation in a 2D tissue.

The Physical Mechanisms of Drosophila Gastrulation: Mesoderm and Endoderm Invagination

A critical juncture in early development is the partitioning of cells that will adopt different fates into three germ layers: the ectoderm, the mesoderm, and the endoderm. This step is achieved through the internalization of specified cells from the outermost surface layer, through a process called gastrulation. In Drosophila, gastrulation is achieved through cell shape changes (i.e., apical constriction) that change tissue curvature and lead to the folding of a surface epithelium. Folding of embryonic tissue results in mesoderm and endoderm invagination, not as individual cells, but as collective tissue units. The tractability of Drosophila as a model system is best exemplified by how much we know about Drosophila gastrulation, from the signals that pattern the embryo to the molecular components that generate force, and how these components are organized to promote cell and tissue shape changes. For mesoderm invagination, graded signaling by the morphogen, Spätzle, sets up a gradient in transcriptional activity that leads to the expression of a secreted ligand (Folded gastrulation) and a transmembrane protein (T48). Together with the GPCR Mist, which is expressed in the mesoderm, and the GPCR Smog, which is expressed uniformly, these signals activate heterotrimeric G-protein and small Rho-family G-protein signaling to promote apical contractility and changes in cell and tissue shape. A notable feature of this signaling pathway is its intricate organization in both space and time. At the cellular level, signaling components and the cytoskeleton exhibit striking polarity, not only along the apical-basal cell axis, but also within the apical domain. Furthermore, gene expression controls a highly choreographed chain of events, the dynamics of which are critical for primordium invagination; it does not simply throw the cytoskeletal "on" switch. Finally, studies of Drosophila gastrulation have provided insight into how global tissue mechanics and movements are intertwined as multiple tissues simultaneously change shape. Overall, these studies have contributed to the view that cells respond to forces that propagate over great distances, demonstrating that cellular decisions, and, ultimately, tissue shape changes, proceed by integrating cues across an entire embryo.

Genetic signatures of dominance hierarchies reveal conserved cis-regulatory and brain gene expression underlying aggression in a facultatively social bee

Agonistic interactions among individuals can result in the formation of dominance hierarches that can reinforce individual behavior and social status. Such dominance hierarches precede the establishment of reproductive dominance, division of labor and caste formation in highly social insect taxa. As such, deciphering the molecular basis of aggression is fundamental in understanding the mechanisms of social evolution. Assessing the proximate mechanisms of aggression in incipiently social bees can provide insights into the foundations of genomic mechanisms of social behavior. Here, we measured the effects of aggression on brain gene expression in the incipiently social bee, Ceratina australensis. We examine the brain transcriptomic differences between individuals who have experienced recurrent winning, losing, or a change in rank during repeated encounters. Using comparative analyses across taxa, we identify deeply conserved candidate genes, pathways, and regulatory networks for the formation of social hierarchies.

GLIS1-3 transcription factors: critical roles in the regulation of multiple physiological processes and diseases

Krüppel-like zinc finger proteins form one of the largest families of transcription factors. They function as key regulators of embryonic development and a wide range of other physiological processes, and are implicated in a variety of pathologies. GLI-similar 1-3 (GLIS1-3) constitute a subfamily of Krüppel-like zinc finger proteins that act either as activators or repressors of gene transcription. GLIS3 plays a critical role in the regulation of multiple biological processes and is a key regulator of pancreatic β cell generation and maturation, insulin gene expression, thyroid hormone biosynthesis, spermatogenesis, and the maintenance of normal kidney functions. Loss of GLIS3 function in humans and mice leads to the development of several pathologies, including neonatal diabetes and congenital hypothyroidism, polycystic kidney disease, and infertility. Single nucleotide polymorphisms in GLIS3 genes have been associated with increased risk of several diseases, including type 1 and type 2 diabetes, glaucoma, and neurological disorders. GLIS2 plays a critical role in the kidney and GLIS2 dysfunction leads to nephronophthisis, an end-stage, cystic renal disease. In addition, GLIS1-3 have regulatory functions in several stem/progenitor cell populations. GLIS1 and GLIS3 greatly enhance reprogramming efficiency of somatic cells into induced embryonic stem cells, while GLIS2 inhibits reprogramming. Recent studies have obtained substantial mechanistic insights into several physiological processes regulated by GLIS2 and GLIS3, while a little is still known about the physiological functions of GLIS1. The localization of some GLIS proteins to the primary cilium suggests that their activity may be regulated by a downstream primary cilium-associated signaling pathway. Insights into the upstream GLIS signaling pathway may provide opportunities for the development of new therapeutic strategies for diabetes, hypothyroidism, and other diseases.

The Insulator Protein CTCF Is Required for Correct Hox Gene Expression, but Not for Embryonic Development in Drosophila

Insulator binding proteins (IBPs) play an important role in regulating gene expression by binding to specific DNA sites to facilitate appropriate gene regulation. There are several IBPs in Drosophila, each defined by their ability to insulate target gene promoters in transgenic assays from the activating or silencing effects of neighboring regulatory elements. Of these, only CCCTC-binding factor (CTCF) has an obvious ortholog in mammals. CTCF is essential for mammalian cell viability and is an important regulator of genome architecture. In flies, CTCF is both maternally deposited and zygotically expressed. Flies lacking zygotic CTCF die as young adults with homeotic defects, suggesting that specific Hox genes are misexpressed in inappropriate body segments. The lack of any major embryonic defects was assumed to be due to the maternal supply of CTCF protein, as maternally contributed factors are often sufficient to progress through much of embryogenesis. Here, we definitively determined the requirement of CTCF for developmental progression in Drosophila We generated animals that completely lack both maternal and zygotic CTCF and found that, contrary to expectation, these mutants progress through embryogenesis and larval life. They develop to pharate adults, which fail to eclose from their pupal case. These mutants show exacerbated homeotic defects compared to zygotic mutants, misexpressing the Hox gene Abdominal-B outside of its normal expression domain early in development. Our results indicate that loss of Drosophila CTCF is not accompanied by widespread effects on gene expression, which may be due to redundant functions with other IBPs. Rather, CTCF is required for correct Hox gene expression patterns and for the viability of adult Drosophila.

The cis-regulatory dynamics of embryonic development at single-cell resolution

Understanding how gene regulatory networks control the progressive restriction of cell fates is a long-standing challenge. Recent advances in measuring single cell gene expression are providing new insights into lineage commitment. However, the regulatory events underlying these changes remain elusive. Here we investigate the dynamics of chromatin regulatory landscapes during embryogenesis at single cell resolution. Using single cell combinatorial indexing assay for transposase accessible chromatin (sci-ATAC-seq)1, we profiled chromatin accessibility in over 20,000 single nuclei from fixed Drosophila embryos spanning three landmark embryonic stages: 2-4 hours (hrs) after egg laying (predominantly stage 5 blastoderm nuclei), when each embryo comprises ~6,000 multipotent cells; 6-8hrs (predominantly stage 10-11), to capture a midpoint in embryonic development when major lineages in the mesoderm and ectoderm are specified; and 10-12hrs (predominantly stage 13), when each of the embryo’s >20,000 cells are undergoing terminal differentiation. Our results reveal spatial heterogeneity in the usage of the regulatory genome prior to gastrulation, a feature that aligns with future cell fate, and nuclei can be temporally ordered along developmental trajectories. During mid-embryogenesis, tissue granularity emerges such that individual cell types can be inferred by their chromatin accessibility, while maintaining a signature of their germ layer of origin. The data reveal overlapping usage of regulatory elements between cells of the endoderm and non-myogenic mesoderm, suggesting a common developmental program reminiscent of the mesendoderm lineage in other species2–4. Altogether, we identify over 30,000 distal regulatory elements exhibiting tissue-specific accessibility. We validated the germ layer specificity of a subset of these predicted enhancers in transgenic embryos, achieving 90% accuracy. Overall, our results demonstrate the power of shotgun single cell profiling of embryos to resolve dynamic changes in the chromatin landscape during development, and to uncover the cis-regulatory programs of metazoan germ layers and cell types.

Revisiting epithelial-mesenchymal transition in cancer metastasis: the connection between epithelial plasticity and stemness

Epithelial‐mesenchymal transition (EMT) is an important process in embryonic development, fibrosis, and cancer metastasis. During cancer progression, the activation of EMT permits cancer cells to acquire migratory, invasive, and stem‐like properties. A growing body of evidence supports the critical link between EMT and cancer stemness. However, contradictory results have indicated that the inhibition of EMT also promotes cancer stemness, and that mesenchymal‐epithelial transition, the reverse process of EMT, is associated with the tumor‐initiating ability required for metastatic colonization. The concept of ‘intermediate‐state EMT’ provides a possible explanation for this conflicting evidence. In addition, recent studies have indicated that the appearance of ‘hybrid’ epithelial‐mesenchymal cells is favorable for the establishment of metastasis. In summary, dynamic changes or plasticity between the epithelial and the mesenchymal states rather than a fixed phenotype is more likely to occur in tumors in the clinical setting. Further studies aimed at validating and consolidating the concept of intermediate‐state EMT and hybrid tumors are needed for the establishment of a comprehensive profile of cancer metastasis.

The Potential Role of Krüppel-Like Zinc-Finger Protein Glis3 in Genetic Diseases and Cancers

Gli-similar 3 (Glis3) belongs to a Glis subfamily of Krüppel-like zinc-finger transcription factors characterized to regulate a set of downstream targets essential for cellular functions, including pancreatic development, β-cell maturation and maintenance, and insulin production. Examination of the DNA-binding domain of Glis3 reveals that this domain contains a repeated cysteine 2/histidine 2 (Cys2/His2) zinc-finger motif in the central region where the recognized DNA sequence binds. The loss of the production of pancreatic hormones, such as insulin 1 and 2, is linked to the down-regulation of β cells-related genes and promotes the apoptotic death of β cells found in mutant Glis3. Although accumulating studies converge on the Glis3 functioning in β cells, recently, there have been developments in the field of Glis3 using knockdown/mutant mice to better understand the role of Glis3 in diseases. The Glis3 mutant mice have been characterized for their propensity to develop congenital hypothyroidism, polycystic kidney disease, and some types of cancer. In this review, we attempt to comprehensively summarize the knowledge of Glis3, including its structure and general function in cells. We also collected and organized the academic achievements related to the possible mechanisms of Glis3-related diseases.

NF-kappaB Is Involved in the Regulation of EMT Genes in Breast Cancer Cells

The metastatic process in breast cancer is related to the expression of the epithelial-to-mesenchymal transition transcription factors (EMT-TFs) SNAIL, SLUG, SIP1 and TWIST1. EMT-TFs and nuclear factor-κB (NF-κB) activation have been associated with aggressiveness and metastatic potential in carcinomas. Here, we sought to examine the role of NF-κB in the aggressive properties and regulation of EMT-TFs in human breast cancer cells. Blocking NF-κB/p65 activity by reducing its transcript and protein levels (through siRNA-strategy and dehydroxymethylepoxyquinomicin [DHMEQ] treatment) in the aggressive MDA-MB-231 and HCC-1954 cell lines resulted in decreased invasiveness and migration, a downregulation of SLUG, SIP1, TWIST1, MMP11 and N-cadherin transcripts and an upregulation of E-cadherin transcripts. No significant changes were observed in the less aggressive cell line MCF-7. Bioinformatics tools identified several NF-κB binding sites along the promoters of SNAIL, SLUG, SIP1 and TWIST1 genes. Through chromatin immunoprecipitation and luciferase reporter assays, the NF-κB/p65 binding on TWIST1, SLUG and SIP1 promoter regions was confirmed. Thus, we suggest that NF-κB directly regulates the transcription of EMT-TF genes in breast cancer. Our findings may contribute to a greater understanding of the metastatic process of this neoplasia and highlight NF-κB as a potential target for breast cancer treatment.

Revisiting the role of the Gcm transcription factor, from master regulator to Swiss army knife

Master genes are known to induce the differentiation of a multipotent cell into a specific cell type. These molecules are often transcription factors that switch on the regulatory cascade that triggers cell specification. Gcm was first described as the master gene of the glial fate in Drosophila as it induces the differentiation of neuroblasts into glia in the developing nervous system. Later on, Gcm was also shown to regulate the differentiation of blood, tendon and peritracheal cells as well as that of neuronal subsets. Thus, the glial master gene is used in at least 4 additional systems to promote differentiation. To understand the numerous roles of Gcm, we recently reported a genome-wide screen of Gcm direct targets in the Drosophila embryo. This screen provided new insight into the role and mode of action of this powerful transcription factor, notably on the interactions between Gcm and major differentiation pathways such as the Hedgehog, Notch and JAK/STAT. Here, we discuss the mode of action of Gcm in the different systems, we present new tissues that require Gcm and we revise the concept of 'master gene'.

Novel Molecular Markers for Breast Cancer

The use of molecular biomarkers assures that breast cancer (BC) patients receive optimal treatment. Established biomarkers, such as estrogen receptor, progesterone receptor, HER2, and Ki67, have been playing significant roles in the subcategorization of BC to predict the prognosis and decide the specific therapy to each patient. Antihormonal therapy using 4-hydroxytamoxifen or aromatase inhibitors have been employed in patients whose tumor cells express hormone receptors, while monoclonal antibody to HER2 has been administered to HER2-positive BCs. Although new therapeutic agents have been developed in the past few decades, many patients still die of the disease due to relapse; thus, novel molecular markers that predict therapeutic failure and those that can be targets for specific therapy are expected. We have chosen four of such molecules by reviewing recent publications, which are cyclin E, B-Myb, Twist, and DMP1β. The oncogenicity of these molecules has been demonstrated in vivo and/or in vitro through studies using transgenic mice or siRNAs, and their expressions have been shown to be associated with shortened overall or disease-free survival of BC patients. The former three molecules have been shown to accelerate epithelial-mesenchymal transition that is often associated with cancer stem cell-ness and metastasis; all these four can be novel therapeutic targets as well. Thus, large prospective studies employing immunohistochemistry will be needed to establish the predictive values of these molecules in patients with BC.

Twist1 is involved in trophoblast syncytialization by regulating GCM1

INTRODUCTION

The multinucleated syncytiotrophoblast (STB) is maintained and regenerated by the fusion of underlying cytotrophoblast cells (CTBs) and is responsible for a number of functions in the human placenta. Deficiencies in this structure may result in pregnancy-associated diseases. However, the detailed mechanisms underlying trophoblast syncytialization await further investigation.

METHODS

The location of the transcription factor Twist1 in human placental tissues was identified by immunohistochemistry. The expression of Twist1 and glial cells missing-1 (GCM1) was evaluated by qPCR or western blotting in two cell-fusion models including forskolin-induced fusion of BeWo cells and spontaneous syncytialization of CTBs. The key role of Twist1 in trophoblast differentiation was identified using BeWo cells transfected with Twist1-specific siRNA. We investigated the effect of hypoxia on the expression of Twist1 and GCM1 in primary CTBs cultured with 2% oxygen. The Twist1 binding region in the GCM1 gene was detected by chromatin-immunoprecipitation.

RESULTS

Twist1 was expressed in human placental tissues, and the expression of Twist1 and GCM1 increased in a time-dependent manner during spontaneous syncytialization of primary CTBs and forskolin-induced fusion of BeWo cells. A reduction in Twist1 and GCM1 expression was observed under hypoxic conditions and was accompanied by inhibition of trophoblast syncytialization. Moreover, siRNA-mediated silencing of Twist1 resulted in inhibition of BeWo cells fusion and down-regulation of GCM1 expression. Furthermore, Twist1 was found to bind to the E-box-enriched region in intron 2 of the GCM1 gene during forskolin-induced fusion of BeWo cells.

DISCUSSION

The above results suggest that Twist1 is required during trophoblast syncytialization. Twist1 may promote trophoblast syncytialization by regulating the expression of GCM1.

Fear-of-intimacy mediated zinc transport controls the function of Zn-finger transcription factors involved in myogenesis

Zinc is a component of one tenth of all human proteins. Its cellular concentration is tightly regulated because its dyshomeostasis has catastrophic health consequences. Two families of zinc transporters control zinc homeostasis in organisms, but there is little information about their specific developmental roles. We show that the ZIP transporter fear-of-intimacy (foi) is necessary for the formation of Drosophila muscles. In foi mutants, myoblasts segregate normally, but their specification is affected, leading to the formation of a misshapen muscle pattern and distorted midgut. The observed phenotypes could be ascribed to the inactivation of specific zing-finger transcription factors (ZFTFs), supporting the hypothesis that they a consequence of a zinc intracellular depletion. Accordingly, foi phenotypes can be rescued by mesodermal expression of other ZIP members with similar subcellular localization. We propose that Foi acts mostly as a transporter to regulate zinc intracellular homeostasis, thereby impacting on the activity of ZFTFs that control specific developmental processes. Our results additionally suggest a possible explanation for the presence of large numbers of zinc transporters in organisms based on differences in ion transport specificity and/or degrees of activity among transporters.

Origin and diversification of wings: Insights from a neopteran insect

Winged insects underwent an unparalleled evolutionary radiation, but mechanisms underlying the origin and diversification of wings in basal insects are sparsely known compared with more derived holometabolous insects. In the neopteran species Oncopeltus fasciatus, we manipulated wing specification genes and used RNA-seq to obtain both functional and genomic perspectives. Combined with previous studies, our results suggest the following key steps in wing origin and diversification. First, a set of dorsally derived outgrowths evolved along a number of body segments including the first thoracic segment (T1). Homeotic genes were subsequently co-opted to suppress growth of some dorsal flaps in the thorax and abdomen. In T1 this suppression was accomplished by Sex combs reduced, that when experimentally removed, results in an ectopic T1 flap similar to prothoracic winglets present in fossil hemipteroids and other early insects. Global gene-expression differences in ectopic T1 vs. T2/T3 wings suggest that the transition from flaps to wings required ventrally originating cells, homologous with those in ancestral arthropod gill flaps/epipods, to migrate dorsally and fuse with the dorsal flap tissue thereby bringing new functional gene networks; these presumably enabled the T2/T3 wing's increased size and functionality. Third, "fused" wings became both the wing blade and surrounding regions of the dorsal thorax cuticle, providing tissue for subsequent modifications including wing folding and the fit of folded wings. Finally, Ultrabithorax was co-opted to uncouple the morphology of T2 and T3 wings and to act as a general modifier of hindwings, which in turn governed the subsequent diversification of lineage-specific wing forms.

Transcriptome profiles of Penaeus (Marsupenaeus) japonicus animal and vegetal half-embryos: identification of sex determination, germ line, mesoderm, and other developmental genes

There is virtually no knowledge of the molecular events controlling early embryogenesis in Penaeid shrimp. A combination of controlled spawning environment, shrimp embryo micro-dissection techniques, and next-generation sequencing was used to produce transcriptome EST datasets of Penaeus japonicus animal and vegetal half-embryos. Embryos were collected immediately after spawning, and then blastomeres were separated at the two-cell stage and allowed to develop to late gastrulation, then pooled for RNA isolation and cDNA synthesis. Ion Torrent sequencing of cDNA from approximately 500 pooled animal and vegetal half-embryos from multiple spawnings resulted in 560,516 and 493,703 reads, respectively. Reads from each library were assembled and Gene Ontogeny analysis produced 3479 annotated animal contigs and 4173 annotated vegetal contigs, with 159/139 hits for developmental processes in the animal/vegetal contigs, respectively. Contigs were subject to BLAST for selected developmental toolbox genes. Some of the genes found included the sex determination genes sex-lethal and transformer; the germ line genes argonaute 1, boule, germ cell-less, gustavus, maelstrom, mex-3, par-1, pumilio, SmB, staufen, and tudor; the mesoderm genes brachyury, mef2, snail, and twist; the axis determination/segmentation genes β-catenin, deformed, distal-less, engrailed, giant, hairy, hunchback, kruppel, orthodenticle, patched, tailless, and wingless/wnt-8c; and a number of cell-cycle regulators. Animal and vegetal contigs were computationally subtracted from each other to produce sets unique to either half-embryo library. Genes expressed only in the animal half included bmp1, kruppel, maelstrom, and orthodenticle. Genes expressed only in the vegetal half included boule, brachyury, deformed, dorsal, engrailed, hunchback, spalt, twist, and wingless/wnt-8c.

Mechanisms involved in breast cancer liver metastasis

Liver metastasis is a frequent occurrence in patients with breast cancer; however, the available treatments are limited and ineffective. While liver-specific homing of breast cancer cells is an important feature of metastasis, the formation of liver metastases is not random. Indeed, breast cancer cell factors contribute to the liver microenvironment. Major breakthroughs have been achieved recently in understanding breast cancer liver metastasis (BCLM). The process of liver metastasis consists of multiple steps and involves various factors from breast cancer cells and the liver microenvironment. A further understanding of the roles of breast cancer cells and the liver microenvironment is crucial to guide future work in clinical treatments. In this review we discuss the contribution of breast cancer cells and the liver microenvironment to liver metastasis, with the aim to improve therapeutic efficacy for patients with BCLM.

Twist1 induces endothelial differentiation of tumour cells through the Jagged1-KLF4 axis

The mechanisms controlling tumour-induced angiogenesis are presently not clear. In principle, angiogenesis can be achieved through the activation of endothelial cells in existing vessels or by transdifferentiation of tumour cells into endothelial cells. However, whether tumour cells can go through a prior epithelial-mesenchymal transition and further differentiate into endothelial cells remains unknown. Here we show that overexpression of Twist1, a transcriptional regulator that induces and promotes cancer metastasis, leads to endothelial differentiation in head and neck cancer (HNC) cells. Induction of Jagged1 expression by Twist1 is essential for Twist1-induced endothelial differentiation. The Jagged1/Notch signalling subsequently activates KLF4, inducing stem-like properties in HNC cells and conferring them with drug resistance. Our results indicate that the Twist1-Jagged1/KLF4 axis is essential both for transdifferentiation of tumour cells into endothelial cells and for chemoresistance acquisition.

Coordinated repression and activation of two transcriptional programs stabilizes cell fate during myogenesis

Molecular models of cell fate specification typically focus on the activation of specific lineage programs. However, the concurrent repression of unwanted transcriptional networks is also essential to stabilize certain cellular identities, as shown in a number of diverse systems and phyla. Here, we demonstrate that this dual requirement also holds true in the context of Drosophila myogenesis. By integrating genetics and genomics, we identified a new role for the pleiotropic transcriptional repressor Tramtrack69 in myoblast specification. Drosophila muscles are formed through the fusion of two discrete cell types: founder cells (FCs) and fusion-competent myoblasts (FCMs). When tramtrack69 is removed, FCMs appear to adopt an alternative muscle FC-like fate. Conversely, ectopic expression of this repressor phenocopies muscle defects seen in loss-of-function lame duck mutants, a transcription factor specific to FCMs. This occurs through Tramtrack69-mediated repression in FCMs, whereas Lame duck activates a largely distinct transcriptional program in the same cells. Lineage-specific factors are therefore not sufficient to maintain FCM identity. Instead, their identity appears more plastic, requiring the combination of instructive repressive and activating programs to stabilize cell fate.

Enhancer loops appear stable during development and are associated with paused polymerase

Developmental enhancers initiate transcription and are fundamental to our understanding of developmental networks, evolution and disease. Despite their importance, the properties governing enhancer-promoter interactions and their dynamics during embryogenesis remain unclear. At the β-globin locus, enhancer-promoter interactions appear dynamic and cell-type specific, whereas at the HoxD locus they are stable and ubiquitous, being present in tissues where the target genes are not expressed. The extent to which preformed enhancer-promoter conformations exist at other, more typical, loci and how transcription is eventually triggered is unclear. Here we generated a high-resolution map of enhancer three-dimensional contacts during Drosophila embryogenesis, covering two developmental stages and tissue contexts, at unprecedented resolution. Although local regulatory interactions are common, long-range interactions are highly prevalent within the compact Drosophila genome. Each enhancer contacts multiple enhancers, and promoters with similar expression, suggesting a role in their co-regulation. Notably, most interactions appear unchanged between tissue context and across development, arising before gene activation, and are frequently associated with paused RNA polymerase. Our results indicate that the general topology governing enhancer contacts is conserved from flies to humans and suggest that transcription initiates from preformed enhancer-promoter loops through release of paused polymerase.

Subtle changes in motif positioning cause tissue-specific effects on robustness of an enhancer's activity

Deciphering the specific contribution of individual motifs within cis-regulatory modules (CRMs) is crucial to understanding how gene expression is regulated and how this process is affected by sequence variation. But despite vast improvements in the ability to identify where transcription factors (TFs) bind throughout the genome, we are limited in our ability to relate information on motif occupancy to function from sequence alone. Here, we engineered 63 synthetic CRMs to systematically assess the relationship between variation in the content and spacing of motifs within CRMs to CRM activity during development using Drosophila transgenic embryos. In over half the cases, very simple elements containing only one or two types of TF binding motifs were capable of driving specific spatio-temporal patterns during development. Different motif organizations provide different degrees of robustness to enhancer activity, ranging from binary on-off responses to more subtle effects including embryo-to-embryo and within-embryo variation. By quantifying the effects of subtle changes in motif organization, we were able to model biophysical rules that explain CRM behavior and may contribute to the spatial positioning of CRM activity in vivo. For the same enhancer, the effects of small differences in motif positions varied in developmentally related tissues, suggesting that gene expression may be more susceptible to sequence variation in one tissue compared to another. This result has important implications for human eQTL studies in which many associated mutations are found in cis-regulatory regions, though the mechanism for how they affect tissue-specific gene expression is often not understood.

Comparative gene expression analysis of Dtg, a novel target gene of Dpp signaling pathway in the early Drosophila melanogaster embryo

In the early Drosophila melanogaster embryo, Dpp, a secreted molecule that belongs to the TGF-β superfamily of growth factors, activates a set of downstream genes to subdivide the dorsal region into amnioserosa and dorsal epidermis. Here, we examined the expression pattern and transcriptional regulation of Dtg, a new target gene of Dpp signaling pathway that is required for proper amnioserosa differentiation. We showed that the expression of Dtg was controlled by Dpp and characterized a 524-bp enhancer that mediated expression in the dorsal midline, as well as, in the differentiated amnioserosa in transgenic reporter embryos. This enhancer contained a highly conserved region of 48-bp in which bioinformatic predictions and in vitro assays identified three Mad binding motifs. Mutational analysis revealed that these three motifs were necessary for proper expression of a reporter gene in transgenic embryos, suggesting that short and highly conserved genomic sequences may be indicative of functional regulatory regions in D. melanogaster genes. Dtg orthologs were not detected in basal lineages of Dipterans, which unlike D. melanogaster develop two extra-embryonic membranes, amnion and serosa, nevertheless Dtg orthologs were identified in the transcriptome of Musca domestica, in which dorsal ectoderm patterning leads to the formation of a single extra-embryonic membrane. These results suggest that Dtg was recruited as a new component of the network that controls dorsal ectoderm patterning in the lineage leading to higher Cyclorrhaphan flies, such as D. melanogaster and M. domestica.

Gα73B is a downstream effector of JAK/STAT signalling and a regulator of Rho1 in Drosophila haematopoiesis

JAK/STAT signalling regulates many essential developmental processes including cell proliferation and haematopoiesis, whereas its inappropriate activation is associated with the majority of myeloproliferative neoplasias and numerous cancers. Furthermore, high levels of JAK/STAT pathway signalling have also been associated with enhanced metastatic invasion by cancerous cells. Strikingly, gain-of-function mutations in the single Drosophila JAK homologue, Hopscotch, result in haemocyte neoplasia, inappropriate differentiation and the formation of melanised haemocyte-derived 'tumour' masses; phenotypes that are partly orthologous to human gain-of-function JAK2-associated pathologies. Here we show that Gα73B, a novel JAK/STAT pathway target gene, is necessary for JAK/STAT-mediated tumour formation in flies. In addition, although Gα73B does not affect haemocyte differentiation, it does regulate haemocyte morphology and motility under non-pathological conditions. We show that Gα73B is required for constitutive, but not injury-induced, activation of Rho1 and for the localisation of Rho1 into filopodia upon haemocyte activation. Consistent with these results, we also show that Rho1 interacts genetically with JAK/STAT signalling, and that wild-type levels of Rho1 are necessary for tumour formation. Our findings link JAK/STAT transcriptional outputs, Gα73B activity and Rho1-dependent cytoskeletal rearrangements and cell motility, therefore connecting a pathway associated with cancer with a marker indicative of invasiveness. As such, we suggest a mechanism by which JAK/STAT pathway signalling may promote metastasis.

Microarray comparison of anterior and posterior Drosophila wing imaginal disc cells identifies novel wing genes

Signaling between cells in the anterior (A) and posterior (P) compartments directs Drosophila wing disc development and is dependent on expression of the homeodomain transcription factor Engrailed (En) in P cells. Downstream of en, posteriorly expressed Hedgehog (Hh) protein signals across the A/P border to establish a developmental organizer that directs pattern formation and growth throughout the wing primordium. Here we extend investigations of the processes downstream of en by using expression array analysis to compare A and P cells. A total of 102 candidate genes were identified that express differentially in the A and P compartments; four were characterized: Stubble (Sb) expression is restricted to A cells due to repression by en. CG15905, CG16884; CG10200/hase und igel (hui) are expressed in A cells downstream of Hh signaling; and RNA interference for hui, Stubble, and CG16884 revealed that each is essential to wing development.

Highly parallel assays of tissue-specific enhancers in whole Drosophila embryos

Transcriptional enhancers are a primary mechanism by which tissue-specific gene expression is achieved. Despite the importance of these regulatory elements in development, responses to environmental stresses and disease, testing enhancer activity in animals remains tedious, with a minority of enhancers having been characterized. Here we describe 'enhancer-FACS-seq' (eFS) for highly parallel identification of active, tissue-specific enhancers in Drosophila melanogaster embryos. Analysis of enhancers identified by eFS as being active in mesodermal tissues revealed enriched DNA binding site motifs of known and putative, previously uncharacterized mesodermal transcription factors. Naive Bayes classifiers using transcription factor binding site motifs accurately predicted mesodermal enhancer activity. Application of eFS to other cell types and organisms should accelerate the cataloging of enhancers and understanding how transcriptional regulation is encoded in them.

Jeb/Alk signalling regulates the Lame duck GLI family transcription factor in the Drosophila visceral mesoderm

The Jelly belly (Jeb)/Anaplastic Lymphoma Kinase (Alk) signalling pathway regulates myoblast fusion in the circular visceral mesoderm (VM) of Drosophila embryos via specification of founder cells. However, only a limited number of target molecules for this pathway are described. We have investigated the role of the Lame Duck (Lmd) transcription factor in VM development in relationship to Jeb/Alk signal transduction. We show that Alk signalling negatively regulates Lmd activity post-transcriptionally through the MEK/MAPK (ERK) cascade resulting in a relocalisation of Lmd protein from the nucleus to cytoplasm. It has previously been shown that downregulation of Lmd protein is necessary for the correct specification of founder cells. In the visceral mesoderm of lmd mutant embryos, fusion-competent myoblasts seem to be converted to 'founder-like' cells that are still able to build a gut musculature even in the absence of fusion. The ability of Alk signalling to downregulate Lmd protein requires the N-terminal 140 amino acids, as a Lmd(141-866) mutant remains nuclear in the presence of active ALK and is able to drive robust expression of the Lmd downstream target Vrp1 in the developing VM. Our results suggest that Lmd is a target of Jeb/Alk signalling in the VM of Drosophila embryos.

Functional knowledge transfer for high-accuracy prediction of under-studied biological processes

A key challenge in genetics is identifying the functional roles of genes in pathways. Numerous functional genomics techniques (e.g. machine learning) that predict protein function have been developed to address this question. These methods generally build from existing annotations of genes to pathways and thus are often unable to identify additional genes participating in processes that are not already well studied. Many of these processes are well studied in some organism, but not necessarily in an investigator's organism of interest. Sequence-based search methods (e.g. BLAST) have been used to transfer such annotation information between organisms. We demonstrate that functional genomics can complement traditional sequence similarity to improve the transfer of gene annotations between organisms. Our method transfers annotations only when functionally appropriate as determined by genomic data and can be used with any prediction algorithm to combine transferred gene function knowledge with organism-specific high-throughput data to enable accurate function prediction.

We show that diverse state-of-art machine learning algorithms leveraging functional knowledge transfer (FKT) dramatically improve their accuracy in predicting gene-pathway membership, particularly for processes with little experimental knowledge in an organism. We also show that our method compares favorably to annotation transfer by sequence similarity. Next, we deploy FKT with state-of-the-art SVM classifier to predict novel genes to 11,000 biological processes across six diverse organisms and expand the coverage of accurate function predictions to processes that are often ignored because of a dearth of annotated genes in an organism. Finally, we perform in vivo experimental investigation in Danio rerio and confirm the regulatory role of our top predicted novel gene, wnt5b, in leftward cell migration during heart development. FKT is immediately applicable to many bioinformatics techniques and will help biologists systematically integrate prior knowledge from diverse systems to direct targeted experiments in their organism of study.

Due to technical and ethical challenges many human diseases or biological processes are studied in model organisms. Discoveries in these organisms are then transferred back to human or other model organisms. Traditional methods for transferring novel gene function annotations have relied on finding genes with high sequence similarity believed to share evolutionary ancestry. However, sequence similarity does not guarantee a shared functional role in molecular pathways. In this study, we show that functional genomics can complement traditional sequence similarity measures to improve the transfer of gene annotations between organisms. We coupled our knowledge transfer method with current state-of-the-art machine learning algorithms and predicted gene function for 11,000 biological processes across six organisms. We experimentally validated our prediction of wnt5b's involvement in the determination of left-right heart asymmetry in zebrafish. Our results show that functional knowledge transfer can improve the coverage and accuracy of machine learning methods used for gene function prediction in a diverse set of organisms. Such an approach can be applied to additional organisms, and will be especially beneficial in organisms that have high-throughput genomic data with sparse annotations.

Analysis of Snail1 function and regulation by Twist1 in palatal fusion

Palatal fusion is a tightly controlled process which comprises multiple cellular events, including cell movement and differentiation. Midline epithelial seam (MES) degradation is essential to palatal fusion. In this study, we analyzed the function of Snail1 during the degradation of the MES. We also analyzed the mechanism regulating the expression of the Snail1 gene in palatal shelves. Palatal explants treated with Snail1 siRNA did not degrade the MES and E-cadherin was not repressed leading to failure of palatal fusion. Transforming growth factor beta 3 (Tgfβ3) regulated Snail1 mRNA, as Snail1 expression decreased in response to Tgfβ3 neutralizing antibody and a PI-3 kinase (PI3K) inhibitor. Twist1, in collaboration with E2A factors, regulated the expression of Snail1. Twist1/E47 dimers bond to the Snail1 promoter to activate expression. Without E47, Twist1 repressed Snail1 expression. These results support the hypothesis that Tgfβ3 may signal through Twist1 and then Snail1 to downregulate E-cadherin expression during palatal fusion.

Interpreting the regulatory genome: the genomics of transcription factor function in Drosophila melanogaster

Researchers have now had access to the fully sequenced Drosophila melanogaster genome for over a decade, and the sequenced genomes of 11 additional Drosophila species have been available for almost 5 years, with more species' genomes becoming available every year [Adams MD, Celniker SE, Holt RA, et al. The genome sequence of Drosophila melanogaster. Science 2000;287:2185-95; Clark AG, Eisen MB, Smith DR, et al. Evolution of genes and genomes on the Drosophila phylogeny. Nature 2007;450:203-18]. Although the best studied of the D. melanogaster transcription factors (TFs) were cloned before sequencing of the genome, the availability of sequence data promised to transform our understanding of TFs and gene regulatory networks. Sequenced genomes have allowed researchers to generate tools for high-throughput characterization of gene expression levels, genome-wide TF localization and analyses of evolutionary constraints on DNA elements across multiple species. With an estimated 700 DNA-binding proteins in the Drosophila genome, it will be many years before each potential sequence-specific TF is studied in detail, yet the last decade of functional genomics research has already impacted our view of gene regulatory networks and TF DNA recognition.

Aberrant expression of hypoxia-inducible factor 1α, TWIST and E-cadherin is associated with aggressive tumor phenotypes in endometrioid endometrial carcinoma

OBJECTIVE

Intratumoral hypoxia promotes angiogenesis, invasion and epithelial-mesenchymal transition, a pivotal event in tumor metastasis. TWIST is a master regulator of multiple developmental processes and has recently been shown to be the key factor responsible for cancer metastasis via the inhibition of E-cadherin expression, a hallmark of epithelial-mesenchymal transition. This study aimed to determine the expression of hypoxia-inducible factor 1α, TWIST and E-cadherin in patients with endometrioid endometrial carcinoma and to examine their clinical significance in endometrioid endometrial carcinoma progression.

METHODS

Using immunohistochemical and tissue microarray approaches, we evaluated the expression of hypoxia-inducible factor 1α, TWIST and E-cadherin in normal endometrial (n = 35), atypical hyperplasia (n = 28) and endometrioid endometrial carcinoma samples (n = 124). Furthermore, we statistically analyzed the association between these markers, as well as their correlation with clinicopathologic variables.

RESULTS

The expression of hypoxia-inducible factor 1α and TWIST were markedly increased, whereas E-cadherin was decreased, as lesions progressed from normal endometrium to atypical hyperplasia to carcinoma (P < 0.01). Among various clinical parameters, the expression of hypoxia-inducible factor 1α and TWIST was strikingly elevated with aggressive tumor characteristics, including higher pathologic grade, deep myometrial invasion and lymph node involvement (P < 0.05). More importantly, overexpression of hypoxia-inducible factor 1α positively correlated with enhanced TWIST expression in endometrioid endometrial carcinoma samples (r = 0.249, P < 0.01); however, statistical analysis showed a negative relationship between TWIST upregulation and E-cadherin downregulation (r = -0.183, P = 0.042).

CONCLUSIONS

These results demonstrated for the first time that the hypoxia-inducible factor 1α/TWIST/E-cadherin pathway may play a critical role in invasion and metastasis of endometrioid endometrial carcinoma. The combined evaluation of these markers may be useful in predicting aggressive phenotypes and thus prognosis in patients with endometrioid endometrial carcinoma.

Simultaneous detection of duplex DNA oligonucleotides using a SERS-based micro-network gradient chip

We report the development of a programmable surface-enhanced Raman scattering (SERS)-based micro-network gradient platform to simultaneously detect two different types of DNA oligomer mixtures. The utility of this platform was demonstrated by quantitative analysis of two breast cancer-related (BRCA1) DNA oligomer mixtures. To generate on-demand concentration gradients, the microfluidic circuit was designed using an electric-hydraulic analogy. Then a multi-gradient microfluidic channel was fabricated based on the theoretical design of the concentration control module. These micro-network structures automatically produce a series of different concentration gradients by continuously mixing Cy3-labeled DNA oligomers (BRAC1-Mutation) with TAMRA-labeled DNA oligomer (BRAC1-Wild). The SERS signals for different ratios of duplex DNA oligomer mixtures, adsorbed on the surface of silver nanoparticles, were measured under flowing conditions. Total analysis time from serial mixing to SERS detection takes less than 10 min because all experimental conditions are automatically controlled inside the exquisitely designed microfluidic channel. This novel SERS-based DNA sensing technology in a micro-network gradient channel is expected to be a powerful analytical tool to simultaneously detect multiple DNA oligomer mixtures.

Poised RNA polymerase II changes over developmental time and prepares genes for future expression

Poised RNA polymerase II (Pol II) is predominantly found at developmental control genes and is thought to allow their rapid and synchronous induction in response to extracellular signals. How the recruitment of poised RNA Pol II is regulated during development is not known. By isolating muscle tissue from Drosophila embryos at five stages of differentiation, we show that the recruitment of poised Pol II occurs at many genes de novo and this makes them permissive for future gene expression. A comparison with other tissues shows that these changes are stage specific and not tissue specific. In contrast, Polycomb group repression is tissue specific, and in combination with Pol II (the balanced state) marks genes with highly dynamic expression. This suggests that poised Pol II is temporally regulated and is held in check in a tissue-specific fashion. We compare our data with findings in mammalian embryonic stem cells and discuss a framework for predicting developmental programs on the basis of the chromatin state.

Predicting spatial and temporal gene expression using an integrative model of transcription factor occupancy and chromatin state

Precise patterns of spatial and temporal gene expression are central to metazoan complexity and act as a driving force for embryonic development. While there has been substantial progress in dissecting and predicting cis-regulatory activity, our understanding of how information from multiple enhancer elements converge to regulate a gene's expression remains elusive. This is in large part due to the number of different biological processes involved in mediating regulation as well as limited availability of experimental measurements for many of them. Here, we used a Bayesian approach to model diverse experimental regulatory data, leading to accurate predictions of both spatial and temporal aspects of gene expression. We integrated whole-embryo information on transcription factor recruitment to multiple cis-regulatory modules, insulator binding and histone modification status in the vicinity of individual gene loci, at a genome-wide scale during Drosophila development. The model uses Bayesian networks to represent the relation between transcription factor occupancy and enhancer activity in specific tissues and stages. All parameters are optimized in an Expectation Maximization procedure providing a model capable of predicting tissue- and stage-specific activity of new, previously unassayed genes. Performing the optimization with subsets of input data demonstrated that neither enhancer occupancy nor chromatin state alone can explain all gene expression patterns, but taken together allow for accurate predictions of spatio-temporal activity. Model predictions were validated using the expression patterns of more than 600 genes recently made available by the BDGP consortium, demonstrating an average 15-fold enrichment of genes expressed in the predicted tissue over a naïve model. We further validated the model by experimentally testing the expression of 20 predicted target genes of unknown expression, resulting in an accuracy of 95% for temporal predictions and 50% for spatial. While this is, to our knowledge, the first genome-wide approach to predict tissue-specific gene expression in metazoan development, our results suggest that integrative models of this type will become more prevalent in the future.

Development is a complex process in which a single cell gives rise to a multi-cellular organism comprised of diverse cell types and well-organized tissues. This transformation requires tightly coordinated expression, both spatially and temporally, of hundreds to thousands of genes specific to any given tissue. To orchestrate these patterns, gene expression is regulated at multiple steps, from TF binding to cis-regulatory modules, general transcription factor and RNA polymerase II recruitment to promoters, chromatin remodeling, and three-dimensional looping interactions. Despite this level of complexity, the regulation of gene expression is typically modeled in the context of transcription factor binding and a single enhancer's activity as this is where the majority of experimental data is available. Recent advances in the measurement of chromatin modifications and insulator binding during embryogenesis provide new datasets that can be used for modeling gene expression. Here we use a Bayesian approach to integrate all three levels of information to combine the activity of multiple regulatory elements into a single model of a gene's expression, implementing an expectation maximization strategy to overcome the problem of missing data. Importantly, while the data for histone modifications and insulator binding represents merged signals from all cells in the embryo, the model can extract cell type specific and stage-specific predictions on gene expression for hundreds of genes of unknown expression.

Linkage between Twist1 and Bmi1: molecular mechanism of cancer metastasis/stemness and clinical implications

Cancer metastasis is the major cause of cancer-related death despite significant improvements in multimodal cancer therapy. Epithelial-mesenchymal transition (EMT), a major mechanism of cancer metastasis, is a process that generates cells with stem cell-like properties (cancer stemness). Cancer stemness is a concept that describes a minor population of cells (cancer stem cells) residing within a tumour that are able to self-renew and are resistant to conventional therapy. The mechanisms delineating the generation of cancer stemness and its connection to cancer metastasis remain largely unknown. Twist1 is an EMT regulator and increased Twist1 expression, which has prognostic significance in various human cancers, has been widely reported. Bmi1 is a critical component of polycomb repressive complex (PRC) 1, which maintains self-renewal and stemness. Bmi1 is frequently overexpressed in different types of human cancers and can induce drug resistance (Table 2). Recent studies have shown that Twist1 directly activates Bmi1 expression and that these two molecules function together to mediate cancer stemness and EMT. These results present a unique mechanism of EMT-induced cancer metastasis and stemness. Further investigation of the mechanisms of EMT-mediated cancer metastasis and stemness will contribute to the management and treatment of metastatic cancers.

Integrative analysis of the zinc finger transcription factor Lame duck in the Drosophila myogenic gene regulatory network

Contemporary high-throughput technologies permit the rapid identification of transcription factor (TF) target genes on a genome-wide scale, yet the functional significance of TFs requires knowledge of target gene expression patterns, cooperating TFs, and cis-regulatory element (CRE) structures. Here we investigated the myogenic regulatory network downstream of the Drosophila zinc finger TF Lame duck (Lmd) by combining both previously published and newly performed genomic data sets, including ChIP sequencing (ChIP-seq), genome-wide mRNA profiling, cell-specific expression patterns of putative transcriptional targets, analysis of histone mark signatures, studies of TF cooccupancy by additional mesodermal regulators, TF binding site determination using protein binding microarrays (PBMs), and machine learning of candidate CRE motif compositions. Our findings suggest that Lmd orchestrates an extensive myogenic regulatory network, a conclusion supported by the identification of Lmd-dependent genes, histone signatures of Lmd-bound genomic regions, and the relationship of these features to cell-specific gene expression patterns. The heterogeneous cooccupancy of Lmd-bound regions with additional mesodermal regulators revealed that different transcriptional inputs are used to mediate similar myogenic gene expression patterns. Machine learning further demonstrated diverse combinatorial motif patterns within tissue-specific Lmd-bound regions. PBM analysis established the complete spectrum of Lmd DNA binding specificities, and site-directed mutagenesis of Lmd and additional newly discovered motifs in known enhancers demonstrated the critical role of these TF binding sites in supporting full enhancer activity. Collectively, these findings provide insights into the transcriptional codes regulating muscle gene expression and offer a generalizable approach for similar studies in other systems.