Master regulators in development: Views from the Drosophila retinal determination and mammalian pluripotency gene networks

Fatecode enables cell fate regulator prediction using classification-supervised autoencoder perturbation

Cell reprogramming, which guides the conversion between cell states, is a promising technology for tissue repair and regeneration, with the ultimate goal of accelerating recovery from diseases or injuries. To accomplish this, regulators must be identified and manipulated to control cell fate. We propose Fatecode, a computational method that predicts cell fate regulators based only on single-cell RNA sequencing (scRNA-seq) data. Fatecode learns a latent representation of the scRNA-seq data using a deep learning-based classification-supervised autoencoder and then performs in silico perturbation experiments on the latent representation to predict genes that, when perturbed, would alter the original cell type distribution to increase or decrease the population size of a cell type of interest. We assessed Fatecode's performance using simulations from a mechanistic gene-regulatory network model and scRNA-seq data mapping blood and brain development of different organisms. Our results suggest that Fatecode can detect known cell fate regulators from single-cell transcriptomics datasets.

Identification of master regulator genes controlling pathogenic CD4+ T cell fate in inflammatory bowel disease through transcriptional network analysis

Inflammatory bowel diseases (IBD) are a group of chronic inflammatory conditions of the gastrointestinal tract associated with multiple pathogenic factors, including dysregulation of the immune response. Effector CD4+ T cells and regulatory CD4+ T cells (Treg) are central players in maintaining the balance between tolerance and inflammation. Interestingly, genetic modifications in these cells have been implicated in regulating the commitment of specific phenotypes and immune functions. However, the transcriptional program controlling the pathogenic behavior of T helper cells in IBD progression is still unknown. In this study, we aimed to find master transcription regulators controlling the pathogenic behavior of effector CD4+ T cells upon gut inflammation. To achieve this goal, we used an animal model of IBD induced by the transfer of naïve CD4+ T cells into recombination-activating gene 1 (Rag1) deficient mice, which are devoid of lymphocytes. As a control, a group of Rag1-/- mice received the transfer of the whole CD4+ T cells population, which includes both effector T cells and Treg. When gut inflammation progressed, we isolated CD4+ T cells from the colonic lamina propria and spleen tissue, and performed bulk RNA-seq. We identified differentially up- and down-regulated genes by comparing samples from both experimental groups. We found 532 differentially expressed genes (DEGs) in the colon and 30 DEGs in the spleen, mostly related to Th1 response, leukocyte migration, and response to cytokines in lamina propria T-cells. We integrated these data into Gene Regulatory Networks to identify Master Regulators, identifying four up-regulated master gene regulators (Lef1, Dnmt1, Mybl2, and Jup) and only one down-regulated master regulator (Foxo3). The altered expression of master regulators observed in the transcriptomic analysis was confirmed by qRT-PCR analysis and found an up-regulation of Lef1 and Mybl2, but without differences on Dnmt1, Jup, and Foxo3. These two master regulators have been involved in T cells function and cell cycle progression, respectively. We identified two master regulator genes associated with the pathogenic behavior of effector CD4+ T cells in an animal model of IBD. These findings provide two new potential molecular targets for treating IBD.

A Drosophila Model Reveals the Potential Role for mtt in Retinal Disease

Congenital stationary night blindness (CSNB) is a genetically heterogeneous inherited retinal disorder, caused by over 300 mutations in 17 different genes. While there are numerous fly models available for simulating ocular diseases, most are focused on mimicking retinitis pigmentosa (RP), with animal models specifically addressing CSNB limited to mammals. Here, we present a CSNB fly model associated with the mtt gene, utilizing RNA interference (RNAi) to silence the mtt gene in fly eyes (homologous to the mammalian GRM6 gene) and construct a CSNB model. Through this approach, we observed significant defects in the eye structure and function upon reducing mtt expression in fly eyes. This manifested as disruptions in the compound eye lens structure and reduced sensitivity to light responses. These results suggest a critical role for mtt in the function of fly adult eyes. Interestingly, we found that the mtt gene is not expressed in the photoreceptor neurons of adult flies but is localized to the inner lamina neurons. In summary, these results underscore the crucial involvement of mtt in fly retinal function, providing a framework for understanding the pathogenic mechanisms of CSNB and facilitating research into potential therapeutic interventions.

The multifaceted roles of Myb domain-containing proteins in apicomplexan parasites

Apicomplexan parasites are a large and diverse clade of protists responsible for significant diseases of humans and animals. Central to the ability of these parasites to colonize their host and evade immune responses is an expanded repertoire of gene-expression programs that requires the coordinated action of complex transcriptional networks. DNA-binding proteins and chromatin regulators are essential orchestrators of apicomplexan gene expression that often act in concert. Although apicomplexan genomes encode various families of putative DNA-binding proteins, most remain functionally and mechanistically unexplored. This review highlights the versatile role of myeloblastosis (Myb) domain-containing proteins in apicomplexan parasites as transcription factors and chromatin regulators. We explore the diversity of Myb domain structure and use phylogenetic analysis to identify common features across the phylum. This provides a framework to discuss functional heterogeneity and regulation of Myb domain-containing proteins particularly emphasizing their role in parasite differentiation.

Reversing pathological cell states: the road less travelled can extend the therapeutic horizon

Acquisition of omics data advances at a formidable pace. Yet, our ability to utilize these data to control cell phenotypes and design interventions that reverse pathological states lags behind. Here, we posit that cell states are determined by core networks that control cell-wide networks. To steer cell fate decisions, core networks connecting genotype to phenotype must be reconstructed and understood. A recent method, cell state transition assessment and regulation (cSTAR), applies perturbation biology to quantify causal connections and mechanistically models how core networks influence cell phenotypes. cSTAR models are akin to digital cell twins enabling us to purposefully convert pathological states back to physiologically normal states. While this capability has a range of applications, here we discuss reverting oncogenic transformation.

Transcriptional Signatures and Network-Based Approaches Identified Master Regulators Transcription Factors Involved in Experimental Periodontitis Pathogenesis

Periodontitis is a chronic inflammatory disease characterized by the progressive and irreversible destruction of the periodontium. Its aetiopathogenesis lies in the constant challenge of the dysbiotic biofilm, which triggers a deregulated immune response responsible for the disease phenotype. Although the molecular mechanisms underlying periodontitis have been extensively studied, the regulatory mechanisms at the transcriptional level remain unclear. To generate transcriptomic data, we performed RNA shotgun sequencing of the oral mucosa of periodontitis-affected mice. Since genes are not expressed in isolation during pathological processes, we disclose here the complete repertoire of differentially expressed genes (DEG) and co-expressed modules to build Gene Regulatory Networks (GRNs) and identify the Master Transcriptional Regulators of periodontitis. The transcriptional changes revealed 366 protein-coding genes and 42 non-coding genes differentially expressed and enriched in the immune response. Furthermore, we found 13 co-expression modules with different representation degrees and gene expression levels. Our GRN comprises genes from 12 gene clusters, 166 nodes, of which 33 encode Transcription Factors, and 201 connections. Finally, using these strategies, 26 master regulators of periodontitis were identified. In conclusion, combining the transcriptomic analyses with the regulatory network construction represents a powerful and efficient strategy for identifying potential periodontitis-therapeutic targets.

Transcriptional and Histone Acetylation Changes Associated with CRE Elements Expose Key Factors Governing the Regulatory Circuit in the Early Stage of Huntington's Disease Models

Huntington's disease (HD) is a disorder caused by an abnormal expansion of trinucleotide CAG repeats within the huntingtin (Htt) gene. Under normal conditions, the CREB Binding Protein interacts with CREB elements and acetylates Lysine 27 of Histone 3 to direct the expression of several genes. However, mutant Htt causes depletion of CBP, which in turn induces altered histone acetylation patterns and transcriptional deregulation. Here, we have studied a differential expression analysis and H3K27ac variation in 4- and 6-week-old R6/2 mice as a model of juvenile HD. The analysis of differential gene expression and acetylation levels were integrated into Gene Regulatory Networks revealing key regulators involved in the altered transcription cascade. Our results show changes in acetylation and gene expression levels that are related to impaired neuronal development, and key regulators clearly defined in 6-week-old mice are proposed to drive the downstream regulatory cascade in HD. Here, we describe the first approach to determine the relationship among epigenetic changes in the early stages of HD. We determined the existence of changes in pre-symptomatic stages of HD as a starting point for early onset indicators of the progression of this disease.

A positive feedback loop controls Toxoplasma chronic differentiation

Successful infection strategies must balance pathogen amplification and persistence. In the obligate intracellular parasite Toxoplasma gondii this is accomplished through differentiation into dedicated cyst-forming chronic stages that avoid clearance by the host immune system. The transcription factor BFD1 is both necessary and sufficient for stage conversion; however, its regulation is not understood. In this study we examine five factors that are transcriptionally activated by BFD1. One of these is a cytosolic RNA-binding protein of the CCCH-type zinc-finger family, which we name bradyzoite formation deficient 2 (BFD2). Parasites lacking BFD2 fail to induce BFD1 and are consequently unable to fully differentiate in culture or in mice. BFD2 interacts with the BFD1 transcript under stress, and deletion of BFD2 reduces BFD1 protein levels but not messenger RNA abundance. The reciprocal effects on BFD2 transcription and BFD1 translation outline a positive feedback loop that enforces the chronic-stage gene-expression programme. Thus, our findings help explain how parasites both initiate and commit to chronic differentiation. This work provides new mechanistic insight into the regulation of T. gondii persistence, and can be exploited in the design of strategies to prevent and treat these key reservoirs of human infection.

A Mesp1-dependent developmental breakpoint in transcriptional and epigenomic specification of early cardiac precursors

Transcriptional networks governing cardiac precursor cell (CPC) specification are incompletely understood owing, in part, to limitations in distinguishing CPCs from non-cardiac mesoderm in early gastrulation. We leveraged detection of early cardiac lineage transgenes within a granular single-cell transcriptomic time course of mouse embryos to identify emerging CPCs and describe their transcriptional profiles. Mesp1, a transiently expressed mesodermal transcription factor, is canonically described as an early regulator of cardiac specification. However, we observed perdurance of CPC transgene-expressing cells in Mesp1 mutants, albeit mislocalized, prompting us to investigate the scope of the role of Mesp1 in CPC emergence and differentiation. Mesp1 mutant CPCs failed to robustly activate markers of cardiomyocyte maturity and crucial cardiac transcription factors, yet they exhibited transcriptional profiles resembling cardiac mesoderm progressing towards cardiomyocyte fates. Single-cell chromatin accessibility analysis defined a Mesp1-dependent developmental breakpoint in cardiac lineage progression at a shift from mesendoderm transcriptional networks to those necessary for cardiac patterning and morphogenesis. These results reveal Mesp1-independent aspects of early CPC specification and underscore a Mesp1-dependent regulatory landscape required for progression through cardiogenesis.

Functional Resilience of Mutually Repressing Motifs Embedded in Larger Networks

Elucidating the design principles of regulatory networks driving cellular decision-making has important implications for understanding cell differentiation and guiding the design of synthetic circuits. Mutually repressing feedback loops between 'master regulators' of cell fates can exhibit multistable dynamics enabling "single-positive" phenotypes: (high A, low B) and (low A, high B) for a toggle switch, and (high A, low B, low C), (low A, high B, low C) and (low A, low B, high C) for a toggle triad. However, the dynamics of these two motifs have been interrogated in isolation in silico, but in vitro and in vivo, they often operate while embedded in larger regulatory networks. Here, we embed these motifs in complex larger networks of varying sizes and connectivity to identify hallmarks under which these motifs maintain their canonical dynamical behavior. We show that an increased number of incoming edges onto a motif leads to a decay in their canonical stand-alone behaviors. We also show that this decay can be exacerbated by adding self-inhibition but not self-activation loops on the 'master regulators'. These observations offer insights into the design principles of biological networks containing these motifs and can help devise optimal strategies for the integration of these motifs into larger synthetic networks.

The regulatory role of pioneer factors during cardiovascular lineage specification – A mini review

Cardiovascular disease (CVD) remains the number one cause of death worldwide. Ischemic heart disease contributes to heart failure and has considerable morbidity and mortality. Therefore, alternative therapeutic strategies are urgently needed. One class of epigenetic regulators known as pioneer factors has emerged as an important tool for the development of regenerative therapies for the treatment of CVD. Pioneer factors bind closed chromatin and remodel it to drive lineage specification. Here, we review pioneer factors within the cardiovascular lineage, particularly during development and reprogramming and highlight the implications this field of research has for the future development of cardiac specific regenerative therapies.

Regulation of Eye Determination and Regionalization in the Spider Parasteatoda tepidariorum

Animal visual systems are enormously diverse, but their development appears to be controlled by a set of conserved retinal determination genes (RDGs). Spiders are particular masters of visual system innovation, and offer an excellent opportunity to study the evolution of animal eyes. Several RDGs have been identified in spider eye primordia, but their interactions and regulation remain unclear. From our knowledge of RDG network regulation in Drosophila melanogaster, we hypothesize that orthologs of Pax6, eyegone, Wnt genes, hh, dpp, and atonal could play important roles in controlling eye development in spiders. We analyzed the expression of these genes in developing embryos of the spider Parasteatodatepidariorum, both independently and in relation to the eye primordia, marked using probes for the RDG sine oculis. Our results support conserved roles for Wnt genes in restricting the size and position of the eye field, as well as for atonal initiating photoreceptor differentiation. However, we found no strong evidence for an upstream role of Pax6 in eye development, despite its label as a master regulator of animal eye development; nor do eyg, hh or dpp compensate for the absence of Pax6. Conversely, our results indicate that hh may work with Wnt signaling to restrict eye growth, a role similar to that of Sonichedgehog (Shh) in vertebrates.

Linking Molecular Mechanisms and Evolutionary Consequences of Resource Polyphenism

Resource polyphenism—the occurrence of environmentally induced, discrete, and intraspecific morphs showing differential niche use—is taxonomically widespread and fundamental to the evolution of ecological function where it has arisen. Despite longstanding appreciation for the ecological and evolutionary significance of resource polyphenism, only recently have its proximate mechanisms begun to be uncovered. Polyphenism switches, especially those influencing and influenced by trophic interactions, offer a route to integrating proximate and ultimate causation in studies of plasticity, and its potential influence on evolution more generally. Here, we use the major events in generalized polyphenic development as a scaffold for linking the molecular mechanisms of polyphenic switching with potential evolutionary outcomes of polyphenism and for discussing challenges and opportunities at each step in this process. Not only does the study of resource polyphenism uncover interesting details of discrete plasticity, it also illuminates and informs general principles at the intersection of development, ecology, and evolution.

The timing of cell fate decisions is crucial for initiating pattern formation in the Drosophila eye

The eye-antennal disc of Drosophila is composed of three cell layers: a columnar epithelium called the disc proper (DP); an overlying sheet of squamous cells called the peripodial epithelium (PE); and a strip of cuboidal cells that joins the other two cellular sheets to each other and comprises the outer margin (M) of the disc. The M cells play an important role in patterning the eye because it is here that the Hedgehog (Hh), Decapentaplegic (Dpp) and JAK/STAT pathways function to initiate pattern formation. Dpp signaling is lost from the margin of eyes absent (eya) mutant discs and, as a result, the initiation of retinal patterning is blocked. Based on these observations, Eya has been proposed to control the initiation of the morphogenetic furrow via regulation of Dpp signaling within the M. We show that the failure in pattern formation surprisingly results from M cells prematurely adopting a head epidermis fate. This switch in fate normally takes place during pupal development after the eye has been patterned. Our results suggest that the timing of cell fate decisions is essential for correct eye development.

An Evolutionary Perspective on Hox Binding Site Preferences in Two Different Tissues

Transcription factor (TF) networks define the precise development of multicellular organisms. While many studies focused on TFs expressed in specific cell types to elucidate their contribution to cell specification and differentiation, it is less understood how broadly expressed TFs perform their precise functions in the different cellular contexts. To uncover differences that could explain tissue-specific functions of such TFs, we analyzed here genomic chromatin interactions of the broadly expressed Drosophila Hox TF Ultrabithorax (Ubx) in the mesodermal and neuronal tissues using bioinformatics. Our investigations showed that Ubx preferentially interacts with multiple yet tissue-specific chromatin sites in putative regulatory regions of genes in both tissues. Importantly, we found the classical Hox/Ubx DNA binding motif to be enriched only among the neuronal Ubx chromatin interactions, whereas a novel Ubx-like motif with rather low predicted Hox affinities was identified among the regions bound by Ubx in the mesoderm. Finally, our analysis revealed that tissues-specific Ubx chromatin sites are also different with regards to the distribution of active and repressive histone marks. Based on our data, we propose that the tissue-related differences in Ubx binding behavior could be a result of the emergence of the mesoderm as a new germ layer in triploblastic animals, which might have required the Hox TFs to relax their binding specificity.

Shaping an optical dome: The size and shape of the insect compound eye

The insect compound eye is the most abundant eye architecture on earth. It comes in a wide variety of shapes and sizes, which are exquisitely adapted to specific ecosystems. Here, we explore the organisational principles and pathways, from molecular to tissular, that underpin the building of this organ and highlight why it is an excellent model system to investigate the relationship between genes and tissue form. The compound eye offers wide fields of view, high sensitivity in motion detection and infinite depth of field. It is made of an array of visual units called ommatidia, which are precisely tiled in 3D to shape the retinal tissue as a dome-like structure. The eye starts off as a 2D epithelium, and it acquires its 3D organisation as ommatidia get into shape. Each ommatidium is made of a complement of retinal cells, including light-detecting photoreceptors and lens-secreting cells. The lens cells generate the typical hexagonal facet lens that lies atop the photoreceptors so that the eye surface consists of a quasi-crystalline array of these hexagonal facet-lenses. This array is curved to various degree, depending on the size and shape of the eye, and on the region of the retina. This curvature sets the resolution and visual field of the eye and is determined by i) the number and size of the facet lens - large ommatidial lenses can be used to generate flat, higher resolution areas, while smaller facets allow for stronger curvature of the eye, and ii) precise control of the inter facet-lens angle, which determines the optical axis of the each ommatidium. In this review we discuss how combinatorial variation in eye primordium shape, ommatidial number, facet lens size and inter facet-lens angle underpins the wide variety of insect eye shapes, and we explore what is known about the mechanisms that might control these parameters.

Lineage-instructive functions of the E26-transformation-specific-family transcription factor Spi-C in immune cell development and disease

Cell fate decisions during hematopoiesis are the consequence of a complex mixture of inputs from cell-intrinsic and cell-extrinsic factors. In rare cases, expression of a single transcription factor, or a few key factors, may be sufficient to dictate lineage differentiation in a precursor cell. The E26-transformation-specific-family transcription factor Spi-C has emerged as an example of a lineage-instructive factor involved in the generation of mature, specialized subsets of both myeloid and lymphoid cells. Spi-C can instruct differentiation of splenic precursors into red pulp macrophages responsible for phagocytosing senescent red blood cells. In the B cell compartment, Spi-C acts as a key regulator of cell fate decisions at the pro-B to pre-B cell stage and for plasma cell differentiation. Spi-C regulates key genes including Nfkb1, Bach2, Syk, and Blnk to regulate cell cycle entry and B cell differentiation. Here, we review the biology of the lineage-instructive transcription factor Spi-C and its contribution to mechanisms of disease in macrophages and B cells. This article is categorized under: Cancer > Molecular and Cellular Physiology Immune System Diseases > Molecular and Cellular Physiology Infectious Diseases > Genetics/Genomics/Epigenetics.

A Core Transcription Regulatory Circuitry Defining Microglia Cell Identity Inferred from the Reanalysis of Multiple Human Microglia Differentiation Protocols

Microglia, the immune cells in the brain involved in both homeostasis and injury/infection control, play a predominant role in neurodegenerative diseases. In vivo studies on microglia are limited due to the requirement of surgical intervention, which can lead to the destruction of the tissues. Over the last few years, multiple protocols-presenting a variety of strategies-have described microglia differentiation issued from human pluripotent stem cells. Herein, we have reanalyzed the transcriptomes released on six different microglia differentiation protocols and revealed a consensus core of master transcription regulatory circuitry defining microglia identity. Furthermore, we have discussed the major divergencies among the studied protocols and have provided suggestions to further enhance microglia differentiation assays.

In Silico Analysis to Explore Lineage-Independent and -Dependent Transcriptional Programs Associated with the Process of Endothelial and Neural Differentiation of Human Induced Pluripotent Stem Cells

Despite a major interest in understanding how the endothelial cell phenotype is established, the underlying molecular basis of this process is not yet fully understood. We have previously reported the generation of induced pluripotent stem cells (iPS) from human umbilical vein endothelial cells and differentiation of the resulting HiPS back to endothelial cells (Ec-Diff), as well as neural (Nn-Diff) cell lineage that contained both neurons and astrocytes. Furthermore, the identities of these cell lineages were established by gene array analysis. Here, we explored the same arrays to gain insight into the gene alteration processes that accompany the establishment of endothelial vs. non-endothelial neural cell phenotypes. We compared the expression of genes that code for transcription factors and epigenetic regulators when HiPS is differentiated into these endothelial and non-endothelial lineages. Our in silico analyses have identified cohorts of genes that are similarly up- or downregulated in both lineages, as well as those that exhibit lineage-specific alterations. Based on these results, we propose that genes that are similarly altered in both lineages participate in priming the stem cell for differentiation in a lineage-independent manner, whereas those that are differentially altered in endothelial compared to neural cells participate in a lineage-specific differentiation process. Specific GATA family members and their cofactors and epigenetic regulators (DNMT3B, PRDM14, HELLS) with a major role in regulating DNA methylation were among participants in priming HiPS for lineage-independent differentiation. In addition, we identified distinct cohorts of transcription factors and epigenetic regulators whose alterations correlated specifically with the establishment of endothelial vs. non-endothelial neural lineages.

Energy Metabolism Focused Analysis of Sexual Dimorphism in Biological Aging and Hypothesized Sex-specificity in Sirtuin Dependency

The process of biological aging or senescence refers to the gradual loss of homeostasis and subsequent loss of function - leading to higher chances of mortality. Many mechanisms and driving forces have been suggested to facilitate the evolution of a molecular circuit acting as a trade-off between survival and proliferation, resulting in senescence. A major observation on biological aging and longevity in humans and model organisms is the prevalence of significant sexual divergence in the onset, mechanisms and effects of aging associated processes. In the current account, we describe possible mechanisms by which aging, sex and reproduction are evolutionarily intertwined in order to maintain systemic energy homeostasis. We also interrogate existing literature on the sexual dimorphism of genetic, cellular, metabolic, endocrine and epigenetic processes driving cellular and systemic aging. Subsequently, based on available evidence, we propose a hypothetic model of sex-limited decoupling of female longevity from sirtuins, a major family of regulator proteins of the survival-proliferation trade-off. We also provide necessary considerations to be made in order to test the hypothesis and explore the physiological and therapeutic implications of this decoupling event in male and female longevity after reaching reproductive maturity. HYPOTHESIS STATEMENT: Sirtuins provide survival benefits in a sex-nonspecific manner but the dependency on sirtuins in driving metabolic networks after reaching reproductive maturity is evolutionarily decoupled from female longevity.

Interspecies chimeras as a platform for exogenic organ production and transplantation

Chronic diseases are associated with considerable morbidity and mortality. Therefore, new therapeutic strategies are warranted. Here, we provide a brief review outlining the rationale and feasibility for the generation of intraspecies and interspecies chimeras, which one day may serve as a platform for organ transplantation. These strategies are further associated with consideration of scientific and ethical issues.

Automated generation of context-specific gene regulatory networks with a weighted approach in Drosophila melanogaster

The regulation of gene expression is a key factor in the development and maintenance of life in all organisms. Even so, little is known at whole genome scale for most genes and contexts. We propose a method, Tool for Weighted Epigenomic Networks in Drosophila melanogaster (Fly T-WEoN), to generate context-specific gene regulatory networks starting from a reference network that contains all known gene regulations in the fly. Unlikely regulations are removed by applying a series of knowledge-based filters. Each of these filters is implemented as an independent module that considers a type of experimental evidence, including DNA methylation, chromatin accessibility, histone modifications and gene expression. Fly T-WEoN is based on heuristic rules that reflect current knowledge on gene regulation in D. melanogaster obtained from the literature. Experimental data files can be generated with several standard procedures and used solely when and if available. Fly T-WEoN is available as a Cytoscape application that permits integration with other tools and facilitates downstream network analysis. In this work, we first demonstrate the reliability of our method to then provide a relevant application case of our tool: early development of D. melanogaster. Fly T-WEoN together with its step-by-step guide is available at https://weon.readthedocs.io.

Efficient conversion of human induced pluripotent stem cells into microglia by defined transcription factors

Microglia, the immune cells of the central nervous system, play critical roles in brain physiology and pathology. We report a novel approach that produces, within 10 days, the differentiation of human induced pluripotent stem cells (hiPSCs) into microglia (iMG) by forced expression of both SPI1 and CEBPA. High-level expression of the main microglial markers and the purity of the iMG cells were confirmed by RT-qPCR, immunostaining, and flow cytometry analyses. Whole-transcriptome analysis demonstrated that these iMGs resemble human fetal/adult microglia but not human monocytes. Moreover, these iMGs exhibited appropriate physiological functions, including various inflammatory responses, ADP/ATP-evoked migration, and phagocytic ability. When co-cultured with hiPSC-derived neurons, the iMGs respond and migrate toward injured neurons. This study has established a protocol for the rapid conversion of hiPSCs into functional iMGs, which should facilitate functional studies of human microglia using different disease models and also help with drug discovery.

•

Efficient generation of human iMGs from iPSCs by forced expression of SPI1 and CEBPA

•

The transcriptome profile of iMGs resembles that of human primary microglia

•

The iMG cells possess appropriate physiological functioning

•

An iN-iMG co-culture model is established for studying neuron-microglia interactions

Master regulators as order parameters of gene expression states

Cell type-specific gene expression patterns are represented as memory states of a Hopfield neural network model. It is shown that order parameters of this model can be interpreted as concentrations of master transcription regulators that form concurrent positive feedback loops with a large number of downstream regulated genes. The order parameter free energy then defines an epigenetic landscape in which local minima correspond to stable cell states. The model is applied to gene expression data in the context of hematopoiesis.

Brachyury controls Ciona notochord fate as part of a feed-forward network

The notochord is a defining feature of the chordates. The transcription factor Brachyury (Bra) is a key regulator of notochord fate but here we show that it is not a unitary master regulator in the model chordate Ciona Ectopic Bra expression only partially reprograms other cell types to a notochord-like transcriptional profile and a subset of notochord-enriched genes is unaffected by CRISPR Bra disruption. We identify Foxa.a and Mnx as potential co-regulators, and find that combinatorial cocktails are more effective at reprogramming other cell types than Bra alone. We reassess the network relationships between Bra, Foxa.a and other components of the notochord gene regulatory network, and find that Foxa.a expression in the notochord is regulated by vegetal FGF signaling. It is a direct activator of Bra expression and has a binding motif that is significantly enriched in the regulatory regions of notochord-enriched genes. These and other results indicate that Bra and Foxa.a act together in a regulatory network dominated by positive feed-forward interactions, with neither being a classically defined master regulator.

Regulatory network analysis of Paneth cell and goblet cell enriched gut organoids using transcriptomics approaches

The epithelial lining of the small intestine consists of multiple cell types, including Paneth cells and goblet cells, that work in cohort to maintain gut health. 3D in vitro cultures of human primary epithelial cells, called organoids, have become a key model to study the functions of Paneth cells and goblet cells in normal and diseased conditions. Advances in these models include the ability to skew differentiation to particular lineages, providing a useful tool to study cell type specific function/dysfunction in the context of the epithelium. Here, we use comprehensive profiling of mRNA, microRNA and long non-coding RNA expression to confirm that Paneth cell and goblet cell enrichment of murine small intestinal organoids (enteroids) establishes a physiologically accurate model. We employ network analysis to infer the regulatory landscape altered by skewing differentiation, and using knowledge of cell type specific markers, we predict key regulators of cell type specific functions: Cebpa, Jun, Nr1d1 and Rxra specific to Paneth cells, Gfi1b and Myc specific for goblet cells and Ets1, Nr3c1 and Vdr shared between them. Links identified between these regulators and cellular phenotypes of inflammatory bowel disease (IBD) suggest that global regulatory rewiring during or after differentiation of Paneth cells and goblet cells could contribute to IBD aetiology. Future application of cell type enriched enteroids combined with the presented computational workflow can be used to disentangle multifactorial mechanisms of these cell types and propose regulators whose pharmacological targeting could be advantageous in treating IBD patients with Crohn's disease or ulcerative colitis.

New regulators of Drosophila eye development identified from temporal transcriptome changes

In the last larval instar, uncommitted progenitor cells in the Drosophila eye primordium start to adopt individual retinal cell fates, arrest their growth and proliferation, and initiate terminal differentiation into photoreceptor neurons and other retinal cell types. To explore the regulation of these processes, we have performed mRNA-Seq studies of the larval eye and antennal primordial at multiple developmental stages. A total of 10,893 fly genes were expressed during these stages and could be adaptively clustered into gene groups, some of whose expression increases or decreases in parallel with the cessation of proliferation and onset of differentiation. Using in situ hybridization of a sample of 98 genes to verify spatial and temporal expression patterns, we estimate that 534 genes or more are transcriptionally upregulated during retinal differentiation, and 1367 or more downregulated as progenitor cells differentiate. Each group of co-expressed genes is enriched for regulatory motifs recognized by co-expressed transcription factors, suggesting that they represent coherent transcriptional regulatory programs. Using available mutant strains, we describe novel roles for the transcription factors SoxNeuro (SoxN), H6-like homeobox (Hmx), CG10253, without children (woc), Structure specific recognition protein (Ssrp), and multisex combs (mxc).

The Hox Transcription Factor Ubx Ensures Somatic Myogenesis by Suppressing the Mesodermal Master Regulator Twist

Early lineage-specific master regulators are essential for the specification of cell types. However, once cells are committed to a specific fate, it is critical to restrict the activity of such factors to enable differentiation. To date, it remains unclear how these factors are silenced. Using the Drosophila mesoderm as a model and a comparative genomic approach, we identify the Hox transcription factor Ultrabithorax (Ubx) to be critical for the repression of the master regulator Twist. Mesoderm-specific Ubx loss-of-function experiments using CRISPR-Cas9 and overexpression studies demonstrate that Ubx majorly impacts twist transcription. A mechanistic analysis reveals that Ubx requires the NK-homeodomain protein Tinman to bind to the twist promoter. Furthermore, we find these factor interactions to be critical for silencing by recruiting the Polycomb DNA binding protein Pleiohomeotic. Altogether, our data reveal that Ubx is a critical player in mediating the silencing of Twist, which is crucial for coordinated muscle differentiation.

Scale free topology as an effective feedback system

Biological networks are often heterogeneous in their connectivity pattern, with degree distributions featuring a heavy tail of highly connected hubs. The implications of this heterogeneity on dynamical properties are a topic of much interest. Here we show that interpreting topology as a feedback circuit can provide novel insights on dynamics. Based on the observation that in finite networks a small number of hubs have a disproportionate effect on the entire system, we construct an approximation by lumping these nodes into a single effective hub, which acts as a feedback loop with the rest of the nodes. We use this approximation to study dynamics of networks with scale-free degree distributions, focusing on their probability of convergence to fixed points. We find that the approximation preserves convergence statistics over a wide range of settings. Our mapping provides a parametrization of scale free topology which is predictive at the ensemble level and also retains properties of individual realizations. Specifically, outgoing hubs have an organizing role that can drive the network to convergence, in analogy to suppression of chaos by an external drive. In contrast, incoming hubs have no such property, resulting in a marked difference between the behavior of networks with outgoing vs. incoming scale free degree distribution. Combining feedback analysis with mean field theory predicts a transition between convergent and divergent dynamics which is corroborated by numerical simulations. Furthermore, they highlight the effect of a handful of outlying hubs, rather than of the connectivity distribution law as a whole, on network dynamics.

Nature abounds with complex networks of interacting elements—from the proteins in our cells, through neural networks in our brains, to species interacting in ecosystems. In all of these fields, the relation between network structure and dynamics is an important research question. A recurring feature of natural networks is their heterogeneous structure: individual elements exhibit a huge diversity of connectivity patterns, which complicates the understanding of network dynamics. To address this problem, we devised a simplified approximation for complex structured networks which captures their dynamical properties. Separating out the largest “hubs”—a small number of nodes with disproportionately high connectivity—we represent them by a single node linked to the rest of the network. This enables us to borrow concepts from control theory, where a system’s output is linked back to itself forming a feedback loop. In this analogy, hubs in heterogeneous networks implement a feedback circuit with the rest of the network. The analogy reveals how these hubs can coordinate the network and drive it more easily towards stable states. Our approach enables analyzing dynamical properties of heterogeneous networks, which is difficult to achieve with existing techniques. It is potentially applicable to many fields where heterogeneous networks are important.

The PAX-SIX-EYA-DACH network modulates GATA-FOG function in fly hematopoiesis and human erythropoiesis

The GATA and PAX-SIX-EYA-DACH transcriptional networks (PSEDNs) are essential for proper development across taxa. Here, we demonstrate novel PSEDN roles in vivo in Drosophila hematopoiesis and in human erythropoiesis in vitro Using Drosophila genetics, we show that PSEDN members function with GATA to block lamellocyte differentiation and maintain the prohemocyte pool. Overexpression of human SIX1 stimulated erythroid differentiation of human erythroleukemia TF1 cells and primary hematopoietic stem-progenitor cells. Conversely, SIX1 knockout impaired erythropoiesis in both cell types. SIX1 stimulation of erythropoiesis required GATA1, as SIX1 overexpression failed to drive erythroid phenotypes and gene expression patterns in GATA1 knockout cells. SIX1 can associate with GATA1 and stimulate GATA1-mediated gene transcription, suggesting that SIX1-GATA1 physical interactions contribute to the observed functional interactions. In addition, both fly and human SIX proteins regulated GATA protein levels. Collectively, our findings demonstrate that SIX proteins enhance GATA function at multiple levels, and reveal evolutionarily conserved cooperation between the GATA and PSEDN networks that may regulate developmental processes beyond hematopoiesis.

Widespread activation of developmental gene expression characterized by PRC1-dependent chromatin looping

Polycomb repressive complexes 1 and 2 have been historically described as transcriptional repressors, but recent reports suggest that PRC1 might also support activation, although the underlying mechanisms remain elusive. Here, we show that stage-specific PRC1 binding at a subset of active promoters and enhancers during Drosophila development coincides with the formation of three-dimensional (3D) loops, an increase in expression during development and repression in PRC1 mutants. Dissection of the dachshund locus indicates that PRC1-anchored loops are versatile architectural platforms that persist when surrounding genes are transcriptionally active and fine-tune their expression. The analysis of RING1B binding profiles and 3D contacts during neural differentiation in mice suggests that this role is conserved in mammals.

The β-oxidation pathway is downregulated during diapause termination in Calanus copepods

Calanus copepods are keystone species in marine ecosystems, mainly due to their high lipid content, which is a nutritious food source for e.g. juvenile fish. Accumulated lipids are catabolized to meet energy requirements during dormancy (diapause), which occurs during the last copepodite stage (C5). The current knowledge of lipid degradation pathways during diapause termination is limited. We characterized changes in lipid fullness and generated transcriptional profiles in C5s during termination of diapause and progression towards adulthood. Lipid fullness of C5s declined linearly during developmental progression, but more β-oxidation genes were upregulated in early C5s compared to late C5s and adults. We identified four possible master regulators of energy metabolism, which all were generally upregulated in early C5s, compared to late C5s and adults. We discovered that one of two enzymes in the carnitine shuttle is absent from the calanoid copepod lineage. Based on the geographical location of the sampling site, the field-samples were initially presumed to consist of C. finmarchicus. However, the identification of C. glacialis in some samples underlines the need for performing molecular analyses to reliably identify Calanus species. Our findings contributes to a better understanding of molecular events occurring during diapause and diapause termination in calanoid copepods.

Aniridia and Axenfeld-Rieger Syndrome: Clinical presentations, molecular genetics and current/emerging therapies

Aniridia and Axenfeld-Rieger Syndrome are related, human ocular disorders that are typically inherited in an autosomal dominant manner. Both result from incorrect development of the eye and have, as their most serious consequences, elevated risk to develop the blinding condition glaucoma. This review will focus on describing the clinical presentations of Aniridia and Axenfeld-Rieger Syndrome as well as the molecular genetics and current and emerging therapies used to treat patients.

Novel Strategies for the Generation of Neuronal Diversity: Lessons From the Fly Visual System

Among all organs of an adult animal, the central nervous system stands out because of its vast complexity and morphological diversity. During early development, the entire central nervous system develops from an apparently homogenous group of progenitors that differentiate into all neural cell types. Therefore, understanding the molecular and genetic mechanisms that give rise to the cellular and anatomical diversity of the brain is a key goal of the developmental neurobiology field. With this aim in mind, the development of the central nervous system of model organisms has been extensively studied. From more than a century, the mechanisms of neurogenesis have been studied in the fruit fly Drosophila melanogaster. The visual system comprises one of the major structures of the Drosophila brain. The visual information is collected by the eye-retina photoreceptors and then processed by the four optic lobe ganglia: the lamina, medulla, lobula and lobula plate. The molecular mechanisms that originate neuronal diversity in the optic lobe have been unveiled in the past decade. In this article, we describe the early development and differentiation of the lobula plate ganglion, from the formation of the optic placode and the inner proliferation center to the specification of motion detection neurons. We focused specifically on how the precise combination of signaling pathways and cell-specific transcription factors patterns the pool of neural stem cells that generates the different neurons of the motion detection system.

Two modes of transvection at the eyes absent gene of Drosophila demonstrate plasticity in transcriptional regulatory interactions in cis and in trans

For many genes, proper gene expression requires coordinated and dynamic interactions between multiple regulatory elements, each of which can either promote or silence transcription. In Drosophila, the complexity of the regulatory landscape is further complicated by the tight physical pairing of homologous chromosomes, which can permit regulatory elements to interact in trans, a phenomenon known as transvection. To better understand how gene expression can be programmed through cis- and trans-regulatory interactions, we analyzed transvection effects for a collection of alleles of the eyes absent (eya) gene. We find that trans-activation of a promoter by the eya eye-specific enhancers is broadly supported in many allelic backgrounds, and that the availability of an enhancer to act in trans can be predicted based on the molecular lesion of an eya allele. Furthermore, by manipulating promoter availability in cis and in trans, we demonstrate that the eye-specific enhancers of eya show plasticity in their promoter preference between two different transcriptional start sites, which depends on promoter competition between the two potential targets. Finally, we show that certain alleles of eya demonstrate pairing-sensitive silencing resulting from trans-interactions between Polycomb Response Elements (PREs), and genetic and genomic data support a general role for PcG proteins in mediating transcriptional silencing at eya. Overall, our data highlight how eya gene regulation relies upon a complex but plastic interplay between multiple enhancers, promoters, and PREs.

Gene regulation requires interactions between regions of DNA known as regulatory elements, which, in combination, determine where and when a gene will be active or silenced. Some genes use just a few regulatory elements, whereas others rely on highly complex interactions between many different elements that are poorly understood. While we typically imagine regulatory elements interacting with one another along the length of a single chromosome, in a curious phenomenon called transvection, elements can communicate between two different chromosomes that are held in close proximity. Here, we use the study of transvection to better understand how different regulatory elements contribute to the expression of eyes absent (eya), a gene required for proper eye development in Drosophila. Our data show that a class of elements that initiate eya gene expression, called promoters, will compete with one another for activation by eya’s enhancers, a second class of regulatory element, with the promoter that is closest to the enhancers being the favored target for activation. Furthermore, our study of transvection uncovers an important role for a silencing element, called a PRE, in opposing eya gene expression. Overall, our study sheds new light on how different elements combine to produce patterned expression of eya.

Drosophila melanogaster: A Valuable Genetic Model Organism to Elucidate the Biology of Retinitis Pigmentosa

Retinitis pigmentosa (RP) is a complex inherited disease. It is associated with mutations in a wide variety of genes with many different functions. These mutations impact the integrity of rod photoreceptors and ultimately result in the progressive degeneration of rods and cone photoreceptors in the retina, leading to complete blindness. A hallmark of this disease is the variable degree to which symptoms are manifest in patients. This is indicative of the influence of the environment, and/or of the distinct genetic makeup of the individual.The fruit fly, Drosophila melanogaster, has effectively proven to be a great model system to better understand interconnected genetic networks. Unraveling genetic interactions and thereby different cellular processes is relatively easy because more than a century of research on flies has enabled the creation of sophisticated genetic tools to perturb gene function. A remarkable conservation of disease genes across evolution and the similarity of the general organization of the fly and vertebrate photoreceptor cell had prompted research on fly retinal degeneration. To date six fly models for RP, including RP4, RP11, RP12, RP14, RP25, and RP26, have been established, and have provided useful information on RP disease biology. In this chapter, an outline of approaches and experimental specifications are described to enable utilizing or developing new fly models of RP.

Dachshund Depletion Disrupts Mammary Gland Development and Diverts the Composition of the Mammary Gland Progenitor Pool

DACH1 abundance is reduced in human malignancies, including breast cancer. Herein DACH1 was detected among multipotent fetal mammary stem cells in the embryo, among mixed lineage precursors, and in adult basal cells and (ERα⁺) luminal progenitors. Dach1 gene deletion at 6 weeks in transgenic mice reduced ductal branching, reduced the proportion of mammary basal cells (Lin⁻ CD24^med CD29^high) and reduced abundance of basal cytokeratin 5, whereas DACH1 overexpression induced ductal branching, increased Gata3 and Notch1, and expanded mammosphere formation in LA-7 breast cells. Mammary gland-transforming growth factor β (TGF-β) activity, known to reduce ductal branching and to reduce the basal cell population, increased upon Dach1 deletion, associated with increased SMAD phosphorylation. Association of the scaffold protein Smad anchor for receptor activation with Smad2/3, which facilitates TGF-β activation, was reduced by endogenous DACH1. DACH1 increases basal cells, enhances ductal formation and restrains TGF-β activity in vivo.

•

Dach1 is expressed in mammary gland fetal stem cells and adult luminal cells

•

Dach1 expands mammary gland basal/myoepithelial cells

•

Dach1 induces post-natal mammary gland ductal formation

•

Dach1 retrains TGF-β activity in the mammary gland in vivo

Human eye conditions: insights from the fly eye

The fruit fly Drosophila melanogaster has served as an excellent model to study and understand the genetics of many human diseases from cancer to neurodegeneration. Studying the regulation of growth, determination and differentiation of the compound eyes of this fly, in particular, have provided key insights into a wide range of diseases. Here we review the regulation of the development of fly eyes in light of shared aspects with human eye development. We also show how understanding conserved regulatory pathways in eye development together with the application of tools for genetic screening and functional analyses makes Drosophila a powerful model to diagnose and characterize the genetics underlying many human eye conditions, such as aniridia and retinitis pigmentosa. This further emphasizes the importance and vast potential of basic research to underpin applied research including identifying and treating the genetic basis of human diseases.

Glass promotes the differentiation of neuronal and non-neuronal cell types in the Drosophila eye

Transcriptional regulators can specify different cell types from a pool of equivalent progenitors by activating distinct developmental programs. The Glass transcription factor is expressed in all progenitors in the developing Drosophila eye, and is maintained in both neuronal and non-neuronal cell types. Glass is required for neuronal progenitors to differentiate as photoreceptors, but its role in non-neuronal cone and pigment cells is unknown. To determine whether Glass activity is limited to neuronal lineages, we compared the effects of misexpressing it in neuroblasts of the larval brain and in epithelial cells of the wing disc. Glass activated overlapping but distinct sets of genes in these neuronal and non-neuronal contexts, including markers of photoreceptors, cone cells and pigment cells. Coexpression of other transcription factors such as Pax2, Eyes absent, Lozenge and Escargot enabled Glass to induce additional genes characteristic of the non-neuronal cell types. Cell type-specific glass mutations generated in cone or pigment cells using somatic CRISPR revealed autonomous developmental defects, and expressing Glass specifically in these cells partially rescued glass mutant phenotypes. These results indicate that Glass is a determinant of organ identity that acts in both neuronal and non-neuronal cells to promote their differentiation into functional components of the eye.

Pleiotropy in Drosophila organogenesis: Mechanistic insights from Combgap and the retinal determination gene network

Master regulatory transcription factors cooperate in networks to shepherd cells through organogenesis. In the Drosophila eye, a collection of master control proteins known as the retinal determination gene network (RDGN) switches the direction and targets of its output to choreograph developmental transitions, but the molecular partners that enable such regulatory flexibility are not known. We recently showed that two RDGN members, Eyes absent (Eya) and Sine oculis (So), promote exit from the terminal cell cycle known as the second mitotic wave (SMW) to permit differentiation. A search for co-factors identified the ubiquitously expressed Combgap (Cg) as a novel transcriptional partner that impedes cell cycle exit and interferes with Eya-So activity specifically in this context. Here, we argue that Cg acts as a flexible transcriptional platform that contributes to numerous gene expression outcomes by a variety of mechanisms. For example, Cg provides repressive activities that dampen Eya-So output, but not by recruiting Polycomb chromatin-remodeling complexes as it does in other contexts. We propose that master regulators depend on both specifically expressed co-factors that assemble the combinatorial code and broadly expressed partners like Cg that recruit the diverse molecular activities needed to appropriately regulate their target enhancers.

Master Regulatory Transcription Factors in Plant Development: A Blooming Perspective

Transcription factors that trigger major developmental decisions in plants and animals are termed "master regulators". Such master regulators are classically seen as acting on the top of a regulatory hierarchy that determines a complete developmental program, and they usually encode transcription factors. Here, we introduce master regulators of flowering time and flower development as examples to show how analysis of molecular interactions and gene-regulatory networks in plants has changed our view on the molecular mechanisms by which these factors control developmental processes. A picture has emerged that emphasizes a complex combinatorial interplay in determining cell-type transcriptional programs, and a high level of feedback control. The expression of master regulators themselves is usually regulated by multiple factors integrating environmental and endogenous spatiotemporal cues. Master regulatory transcription factors regulate gene expression by different mechanisms, including modifications in chromatin status in the bound regions. A poorly understood phenomenon is how developmental master regulators exert functions in different cell- and organ types. This is especially relevant for those factors that have important functions in several developmental processes.

The rewiring of transcription circuits in evolution

The binding of transcription regulators to cis-regulatory sequences is a key step through which all cells regulate expression of their genes. Due to gains and losses of cis-regulatory sequences and changes in the transcription regulators themselves, the binding connections between regulators and their target genes rapidly change over evolutionary time and constitute a major source of biological novelty. This review covers recent work, carried out in a wide range of species, that addresses the overall extent of these evolutionary changes, their consequences, and some of the molecular mechanisms that lie behind them.

Mutations that impair Eyes absent tyrosine phosphatase activity in vitro reduce robustness of retinal determination gene network output in Drosophila

A limited collection of signaling networks and transcriptional effectors directs the full spectrum of cellular behaviors that comprise development. One mechanism to diversify regulatory potential is to combine multiple biochemical activities into the same protein. Exemplifying this principle of modularity, Eyes absent (Eya), originally identified as a transcriptional co-activator within the retinal determination gene network (RDGN), also harbors tyrosine and threonine phosphatase activities. Although mounting evidence argues for the importance of Eya’s phosphatase activities to mammalian biology, genetic rescue experiments in Drosophila have shown that the tyrosine phosphatase function is dispensable for normal development. In this study, we repeated these rescue experiments in genetically sensitized backgrounds in which the dose of one or more RDGN factor was reduced. Heterozygosity for sine oculis or dachshund, both core RDGN members, compromised the ability of phosphatase-dead eya, but not of the control wild type eya transgene, to rescue the retinal defects and reduced viability associated with eya loss. We speculate that Eya’s tyrosine phosphatase activity, although non-essential, confers robustness to RDGN output.