Gene Regulatory Signaling Research Articles

Cardiac Progenitor Cells of the First and Second Heart Fields.

The heart forms from the first and second heart fields, which contribute to distinct regions of the myocardium. This is supported by clonal analyses, which identify corresponding first and second cardiac cell lineages in the heart. Progenitor cells of the second heart field and its sub-domains are controlled by a gene regulatory network and signaling pathways, which determine their behavior. Multipotent cells in this field can also contribute cardiac endothelial and smooth muscle cells. Furthermore, the skeletal muscles of the head and neck are clonally related to myocardial cells that form the arterial and venous poles of the heart. These lineage relationships, together with the genes that regulate the heart fields, have major implications for congenital heart disease.

In-Depth Quantification ofCell Division andElongation Dynamics at the Tip ofGrowing Arabidopsis Roots Using 4D Microscopy, AI-Assisted Image Processing andData Sonification.

One of the fundamental questions in plant developmental biology is how cell proliferation and cell expansion coordinately determine organ growth and morphology. An amenable system to address this question is the Arabidopsis root tip, where cell proliferation and elongation occur in spatially separated domains, and cell morphologies can easily be observed using a confocal microscope. While past studies revealed numerous elements of root growth regulation including gene regulatory networks, hormone transport and signaling, cell mechanics and environmental perception, how cells divide and elongate under possible constraints from cell lineages and neighboring cell files has not been analyzed quantitatively. This is mainly due to the technical difficulties in capturing cell division and elongation dynamics at the tip of growing roots, as well as an extremely labor-intensive task of tracing the lineages of frequently dividing cells. Here, we developed a motion-tracking confocal microscope and an Artificial Intelligence (AI)-assisted image-processing pipeline that enables semi-automated quantification of cell division and elongation dynamics at the tip of vertically growing Arabidopsis roots. We also implemented a data sonification tool that facilitates human recognition of cell division synchrony. Using these tools, we revealed previously unnoted lineage-constrained dynamics of cell division and elongation, and their contribution to the root zonation boundaries.

'Chip'-ing away at morphogenesis - application of organ-on-chip technologies to study tissue morphogenesis.

Emergent cell behaviors that drive tissue morphogenesis are the integrated product of instructions from gene regulatory networks, mechanics and signals from the local tissue microenvironment. How these discrete inputs intersect to coordinate diverse morphogenic events is a critical area of interest. Organ-on-chip technology has revolutionized the ability to construct and manipulate miniaturized human tissues with organotypic three-dimensional architectures in vitro. Applications of organ-on-chip platforms have increasingly transitioned from proof-of-concept tissue engineering to discovery biology, furthering our understanding of molecular and mechanical mechanisms that operate across biological scales to orchestrate tissue morphogenesis. Here, we provide the biological framework to harness organ-on-chip systems to study tissue morphogenesis, and we highlight recent examples where organ-on-chips and associated microphysiological systems have enabled new mechanistic insight in diverse morphogenic settings. We further highlight the use of organ-on-chip platforms as emerging test beds for cell and developmental biology.

Multi-omic analyses and network biology in cardiovascular disease.

Heart disease remains a leading cause of death in North America and worldwide. Despite advances in therapies, the chronic nature of cardiovascular diseases ultimately results in frequent hospitalizations and steady rates of mortality. Systems biology approaches have provided a new frontier toward unraveling the underlying mechanisms of cell, tissue, and organ dysfunction in disease. Mapping the complex networks of molecular functions across the genome, transcriptome, proteome, and metabolome has enormous potential to advance our understanding of cardiovascular disease, discover new disease biomarkers, and develop novel therapies. Computational workflows to interpret these data-intensive analyses as well as integration between different levels of interrogation remain important challenges in the advancement and application of systems biology-based analyses in cardiovascular research. This review will focus on summarizing the recent developments in network biology-level profiling in the heart, with particular emphasis on modeling of human heart failure. We will provide new perspectives on integration between different levels of large "omics" datasets, including integration of gene regulatory networks, protein-protein interactions, signaling networks, and metabolic networks in the heart.

Gene Regulatory Networks and Signaling Pathways in Palatogenesis and Cleft Palate: A Comprehensive Review.

Palatogenesis is a complex and intricate process involving the formation of the palate through various morphogenetic events highly dependent on the surrounding context. These events comprise outgrowth of palatal shelves from embryonic maxillary prominences, their elevation from a vertical to a horizontal position above the tongue, and their subsequent adhesion and fusion at the midline to separate oral and nasal cavities. Disruptions in any of these processes can result in cleft palate, a common congenital abnormality that significantly affects patient's quality of life, despite surgical intervention. Although many genes involved in palatogenesis have been identified through studies on genetically modified mice and human genetics, the precise roles of these genes and their products in signaling networks that regulate palatogenesis remain elusive. Recent investigations have revealed that palatal shelf growth, patterning, adhesion, and fusion are intricately regulated by numerous transcription factors and signaling pathways, including Sonic hedgehog (Shh), bone morphogenetic protein (Bmp), fibroblast growth factor (Fgf), transforming growth factor beta (Tgf-β), Wnt signaling, and others. These studies have also identified a significant number of genes that are essential for palate development. Integrated information from these studies offers novel insights into gene regulatory networks and dynamic cellular processes underlying palatal shelf elevation, contact, and fusion, deepening our understanding of palatogenesis, and facilitating the development of more efficacious treatments for cleft palate.

Lineage segregation in human pre-implantation embryos is specified by YAP1 and TEAD1.

Lineage segregation in human pre-implantation embryos is specified by YAP1 and TEAD1.

Survey of genomic and physiological characteristics for survival in lymphoma: The NCI genomic data portal

Lymphoma represents a myriad collection of neoplasms that impact lymphocytes. This cancer often involves disrupted cytokine, immune surveillance, and gene regulatory signaling, sometimes with expression of Epstein Barr Virus (EBV). We explored mutation patterns for People experiencing Lymphoma (PeL) in the National Cancer Institute (NCI) Genomic Data Commons (GDC), which contains detailed, deidentified genomic data on 86,046 people who have/had cancer with 2,730,388 distinctive mutations in 21,773 genes. The database included information on 536 (PeL), with the primary focal sample being the n = 30 who had complete mutational genomic data. We used correlations, independent samples t-tests, and linear regression to compare PeL demographics and vital status on mutation numbers, BMI, and mutation deleterious score across functional categories of 23 genes. PeL demonstrated varied patterns of mutated genes, consistent with most other cancer types. The primary PeL gene mutations clustered around five functional protein groups: transcriptional regulatory proteins, TNF/NFKB and cell signaling regulators, cytokine signaling proteins, cell cycle regulators, and immunoglobulins. Diagnosis Age, Birth Year, and BMI negatively (P < 0.05) correlated with Days to Death, and cell cycle mutations negatively correlated (P = 0.004) with survival days (R2 = 0.389). There were commonalities in some PeL for mutations across other cancer types based upon large sequence length, but also for 6 small cell lung cancer genes. Immunoglobulin mutations were prevalent but not for all cases. Research indicates a need for greater personalized genomics and multi-level systems analysis to evaluate facilitators and barriers for lymphoma survival.

MiRNA as a Potential Target for Multiple Myeloma Therapy-Current Knowledge and Perspectives.

Multiple myeloma (MM) is the second most common hematological malignancy. Despite the huge therapeutic progress thanks to the introduction of novel therapies, MM remains an incurable disease. Extensive research is currently ongoing to find new options. MicroRNAs (miRNAs) are small, non-coding RNA molecules that regulate gene expression at a post-transcriptional level. Aberrant expression of miRNAs in MM is common. Depending on their role in MM development, miRNAs have been reported as oncogenes and tumor suppressors. It was demonstrated that specific miRNA alterations using miRNA mimics or antagomirs can normalize the gene regulatory network and signaling pathways in the microenvironment and MM cells. These properties make miRNAs attractive targets in anti-myeloma therapy. However, only a few miRNA-based drugs have been entered into clinical trials. In this review, we discuss the role of the miRNAs in the pathogenesis of MM, their current status in preclinical/clinical trials, and the mechanisms by which miRNAs can theoretically achieve therapeutic benefit in MM treatment.

Human stem cell models to study placode development, function and pathology.

Placodes are embryonic structures originating from the rostral ectoderm that give rise to highly diverse organs and tissues, comprising the anterior pituitary gland, paired sense organs and cranial sensory ganglia. Their development, including the underlying gene regulatory networks and signalling pathways, have been for the most part characterised in animal models. In this Review, we describe how placode development can be recapitulated by the differentiation of human pluripotent stem cells towards placode progenitors and their derivatives, highlighting the value of this highly scalable platform as an optimal in vitro tool to study the development of human placodes, and identify human-specific mechanisms in their development, function and pathology.

Ancestral Stem Cell Reprogramming Genes Active in Hemichordate Regeneration.

Hemichordate enteropneust worms regenerate extensively in a manner that resembles the regeneration for which planaria and hydra are well known. Although hemichordates are often classified as an extant phylogenetic group that may hold ancestral deuterostome body plans at the base of the deuterostome evolutionary line leading to chordates, mammals, and humans, extensive regeneration is not known in any of these more advanced groups. Here we investigated whether hemichordates deploy functional homologs of canonical Yamanaka stem cell reprogramming factors, Oct4, Sox2, Nanog, and Klf4, as they regenerate. These reprogramming factors are not expressed during regeneration of limbs, fins, eyes or other structures that represent the best examples of regeneration in chordates. We first examined Ptychodera flava EST libraries and identified Pf-Pou3, Pf-SoxB1, Pf-Msxlx, and Pf-Klf1/2/4 as most closely related to the Yamanaka factors, respectively. In situ hybridization analyses revealed that all these homologs are expressed in a distinct manner during head regeneration. Furthermore, Pf-Pou3 partially rescued the loss of endogenous Oct4 in mouse embryonic stem cells in maintaining the pluripotency gene expression program. Based on these results, we propose that hemichordates may have co-opted these reprogramming factors for their extensive regeneration or that chordates may have lost the ability to mobilize these factors in response to damage. The robustness of these pluripotency gene circuits in the inner cell mass and in formation of induced pluripotent stem cells from mammalian somatic cells shows that these programs are intact in humans and other mammals and that these circuits may respond to as yet unknown gene regulatory signals, mobilizing full regeneration in hemichordates.

Identification of MAD2L1 as a Potential Biomarker in Hepatocellular Carcinoma via Comprehensive Bioinformatics Analysis.

Background Hepatocellular carcinoma (HCC) is widely acknowledged as a malignant tumor with rapid progression, high recurrence rate, and poor prognosis. At present, there is a paucity of reliable biomarkers at the clinical level to guide the management of HCC and improve patient outcomes. Our research is aimed at assessing the prognostic value of MAD2L1 in HCC. Methods Four datasets, GSE121248, GSE101685, GSE85598, and GSE62232, were selected from the GEO database to analyze differentially expressed genes (DEGs) between HCC and normal liver tissues. After functional analysis, we constructed a protein-protein interaction network (PPI) for DEGs and identified core genes in this network with high connectivity with other genes. We assessed the relationship between core genes and the pathogenesis and prognosis of HCC. Finally, we explored the gene regulatory signaling mechanisms involved in HCC pathogenesis. Results 145 DEGs were screened from the intersection of the four GEO datasets. MAD2L1 was associated with most genes according to the PPI network and was selected as a candidate gene for further study. Survival analysis suggested that high MAD2L1 expression in HCC correlated with a worse prognosis. In addition, real-time quantitative PCR (RT-qPCR), western blot (WB), and immunohistochemistry (IHC) findings suggested that the expression of MAD2L1 was abnormally increased in HCC tissues and cells compared to paraneoplastic tissues and normal hepatocytes. Conclusion We found that high MAD2L1 expression in HCC was significantly associated with overall patient survival and clinical features. We also explored the potential biological properties of this gene.

Regulation of dhurrin pathway gene expression during Sorghum\xa0bicolor development

Main conclusionDevelopmental and organ-specific expression of genes in dhurrin biosynthesis, bio-activation, and recycling offers dynamic metabolic responses optimizing growth and defence responses in Sorghum.Plant defence models evaluate the costs and benefits of resource investments at different stages in the life cycle. Poor understanding of the molecular regulation of defence deployment and remobilization hampers accuracy of the predictions. Cyanogenic glucosides, such as dhurrin are phytoanticipins that release hydrogen cyanide upon bio-activation. In this study, RNA-seq was used to investigate the expression of genes involved in the biosynthesis, bio-activation and recycling of dhurrin in Sorghum bicolor. Genes involved in dhurrin biosynthesis were highly expressed in all young developing vegetative tissues (leaves, leaf sheath, roots, stems), tiller buds and imbibing seeds and showed gene specific peaks of expression in leaves during diel cycles. Genes involved in dhurrin bio-activation were expressed early in organ development with organ-specific expression patterns. Genes involved in recycling were expressed at similar levels in the different organ during development, although post-floral initiation when nutrients are remobilized for grain filling, expression of GSTL1 decreased > tenfold in leaves and NITB2 increased > tenfold in stems. Results are consistent with the establishment of a pre-emptive defence in young tissues and regulated recycling related to organ senescence and increased demand for nitrogen during grain filling. This detailed characterization of the transcriptional regulation of dhurrin biosynthesis, bioactivation and remobilization genes during organ and plant development will aid elucidation of gene regulatory networks and signalling pathways that modulate gene expression and dhurrin levels. In-depth knowledge of dhurrin metabolism could improve the yield, nitrogen use efficiency and stress resilience of Sorghum.

MicroRNA as an Important Target for Anticancer Drug Development.

Cancer has become the second greatest cause of death worldwide. Although there are several different classes of anticancer drugs that are available in clinic, some tough issues like side-effects and low efficacy still need to dissolve. Therefore, there remains an urgent need to discover and develop more effective anticancer drugs. MicroRNAs (miRNAs) are a class of small endogenous non-coding RNAs that regulate gene expression by inhibiting mRNA translation or reducing the stability of mRNA. An abnormal miRNA expression profile was found to exist widely in cancer cell, which induces limitless replicative potential and evading apoptosis. MiRNAs function as oncogenes (oncomiRs) or tumor suppressors during tumor development and progression. It was shown that regulation of specific miRNA alterations using miRNA mimics or antagomirs can normalize the gene regulatory network and signaling pathways, and reverse the phenotypes in cancer cells. The miRNA hence provides an attractive target for anticancer drug development. In this review, we will summarize the latest publications on the role of miRNA in anticancer therapeutics and briefly describe the relationship between abnormal miRNAs and tumorigenesis. The potential of miRNA-based therapeutics for anticancer treatment has been critically discussed. And the current strategies in designing miRNA targeting therapeutics are described in detail. Finally, the current challenges and future perspectives of miRNA-based therapy are conferred.

Nuclear Receptors in Asthma: Empowering Classical Molecules Against a Contemporary Ailment.

The escalation in living standards and adoption of ‘Western lifestyle’ has an allied effect on the increased allergy and asthma burden in both developed and developing countries. Current scientific reports bespeak an association between allergic diseases and metabolic dysfunction; hinting toward the critical requirement of organized lifestyle and dietary habits. The ubiquitous nuclear receptors (NRs) translate metabolic stimuli into gene regulatory signals, integrating diet inflences to overall developmental and physiological processes. As a consequence of such promising attributes, nuclear receptors have historically been at the cutting edge of pharmacy world. This review discusses the recent findings that feature the cardinal importance of nuclear receptors and how they can be instrumental in modulating current asthma pharmacology. Further, it highlights a possible future employment of therapy involving dietary supplements and synthetic ligands that would engage NRs and aid in eliminating both asthma and linked comorbidities. Therefore, uncovering new and evolving roles through analysis of genomic changes would represent a feasible approach in both prevention and alleviation of asthma.

A hidden integral structure endows absolute concentration robust systems with resilience to dynamical concentration disturbances.

Biochemical systems that express certain chemical species of interest at the same level at any positive steady state are called 'absolute concentration robust' (ACR). These species behave in a stable, predictable way, in the sense that their expression is robust with respect to sudden changes in the species concentration, provided that the system reaches a (potentially new) positive steady state. Such a property has been proven to be of importance in certain gene regulatory networks and signaling systems. In the present paper, we mathematically prove that a well-known class of ACR systems studied by Shinar and Feinberg in 2010 hides an internal integral structure. This structure confers these systems with a higher degree of robustness than was previously known. In particular, disturbances much more general than sudden changes in the species concentrations can be rejected, and robust perfect adaptation is achieved. Significantly, we show that these properties are maintained when the system is interconnected with other chemical reaction networks. This key feature enables the design of insulator devices that are able to buffer the loading effect from downstream systems-a crucial requirement for modular circuit design in synthetic biology. We further note that while the best performance of the insulators are achieved when these act at a faster timescale than the upstream module (as typically required), it is not necessary for them to act on a faster timescale than the downstream module in our construction.

Genome-wide landscape establishes novel association signals for metabolic traits in the Arab population

While the Arabian population has a high prevalence of metabolic disorders, it has not been included in global studies that identify genetic risk loci for metabolic traits. Determining the transferability of such largely Euro-centric established risk loci is essential to transfer the research tools/resources, and drug targets generated by global studies to a broad range of ethnic populations. Further, consideration of populations such as Arabs, that are characterized by consanguinity and a high level of inbreeding, can lead to identification of novel risk loci. We imputed published GWAS data from two Kuwaiti Arab cohorts (n = 1434 and 1298) to the 1000 Genomes Project haplotypes and performed meta-analysis for associations with 13 metabolic traits. We compared the observed association signals with those established for metabolic traits. Our study highlighted 70 variants from 9 different genes, some of which have established links to metabolic disorders. By relaxing the genome-wide significance threshold, we identified ‘novel’ risk variants from 11 genes for metabolic traits. Many novel risk variant association signals were observed at or borderline to genome-wide significance. Furthermore, 349 previously established variants from 187 genes were validated in our study. Pleiotropic effect of risk variants on multiple metabolic traits were observed. Fine-mapping illuminated rs7838666/CSMD1 rs1864163/CETP and rs112861901/[INTS10,LPL] as candidate causal variants influencing fasting plasma glucose and high-density lipoprotein levels. Computational functional analysis identified a variety of gene regulatory signals around several variants. This study enlarges the population ancestry diversity of available GWAS and elucidates new variants in an ethnic group burdened with metabolic disorders.

Dynamic organization of intracellular organelle networks.

Intracellular organelles are membrane-bound and biochemically distinct compartments constructed to serve specialized functions in eukaryotic cells. Through extensive interactions, they form networks to coordinate and integrate their specialized functions for cell physiology. A fundamental property of these organelle networks is that they constantly undergo dynamic organization via membrane fusion and fission to remodel their internal connections and to mediate direct material exchange between compartments. The dynamic organization not only enables them to serve critical physiological functions adaptively but also differentiates them from many other biological networks such as gene regulatory networks and cell signaling networks. This review examines this fundamental property of the organelle networks from a systems point of view. The focus is exclusively on homotypic networks formed by mitochondria, lysosomes, endosomes, and the endoplasmic reticulum, respectively. First, key mechanisms that drive the dynamic organization of these networks are summarized. Then, several distinct organizational properties of these networks are highlighted. Next, spatial properties of the dynamic organization of these networks are emphasized, and their functional implications are examined. Finally, some representative molecular machineries that mediate the dynamic organization of these networks are surveyed. Overall, the dynamic organization of intracellular organelle networks is emerging as a fundamental and unifying paradigm in the internal organization of eukaryotic cells. This article is categorized under: Metabolic Diseases > Molecular and Cellular Physiology.

Abstract 4409: Towards cell-type-specific gene regulation in heterogeneous cancer cells

Abstract The treatment of cancer suffers from diverse drug responsiveness mainly due to high heterogeneity among cell types/states. Sophisticated gene regulatory mechanisms define and maintain cancer types, and in turn, underlie different cancer evolution trends. Within a global gene regulatory system, a regulon represents a maximal group of genes co-regulated by the same transcription regulator (TR, including transcription factors and long non-coding RNAs). Empowered by single-cell RNA-Seq, substantial gene expression analyses have been carried out to reveal regulon activities and construct gene regulatory networks in melanoma, mammary tumors, etc. A clear and profound understanding of these regulons provides an effective way to characterize the heterogeneous regulatory mechanisms, as it allows us to pinpoint crucial regulators and their targeting genes encoded in diverse cell states along with the development of cancer. Hence, for the first time, we define an alternative regulon (AR) as a group of chromatin accessible genes that co-regulated by the same TR in a cell type. Not surprisingly, a TR can regulate multiple heterogeneous ARs with diverse gene availability and activity in different cell types, making the computational identification of ARs a tremendous challengeable problem. Here, we build a first-of-its-kind web server, CeRIS, to predict cell-type-specific ARs based on single-cell RNA-Seq data. Five steps are included: (i) pre-processing, including gene imputation (DrImpute), normalization (scran), and quality control (Seurat3.0); (ii) representing gene regulatory signals by the in-house left-truncated mixture Gaussian model; (iii) clustering cell types using Seurat3.0 and identifying active gene modules in each cell type based on the in-house biclustering method (QUBIC2); (iv) finding conserved cis-regulatory motif patterns for active gene modules using the in-house tool DMINDA2.0 and predicting ARs using TR-motif databases; and (v) ARs with component genes being uniquely accessible or highly expressed in one specific cell type than other cell types are recognized as cell-type-specific ARs. The performance of CeRIS was evaluated using &amp;gt;30 datasets and it outperformed existing tools in terms of the percentage of CTSRs in predicted regulons, marker gene coverage, and pathway enrichment. An applicative study on multiple myeloma suggested that CeRIS can identify novel onco-CTSRs including oncogenes and driver transcription factors in multiple myeloma. Another integrative study of single-cell multi-omics data from human lung indicated the alternative regulons across cell types that can be used to generate heterogeneous landscapes for gene regulatory mechanisms. We believe such novel characteristics of a regulon can be critically useful for cell reprogramming as it can show the transition mechanisms between cell types influenced by the specific dynamic changes of gene regulations. Citation Format: Qin Ma, Anjun Ma, Cankun Wang, Yang Li, Chi Zhang. Towards cell-type-specific gene regulation in heterogeneous cancer cells [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 4409.

Unraveling the Genomic-Epigenomic Interaction Landscape in Triple Negative and Non-Triple Negative Breast Cancer.

Background: The recent surge of next generation sequencing of breast cancer genomes has enabled development of comprehensive catalogues of somatic mutations and expanded the molecular classification of subtypes of breast cancer. However, somatic mutations and gene expression data have not been leveraged and integrated with epigenomic data to unravel the genomic-epigenomic interaction landscape of triple negative breast cancer (TNBC) and non-triple negative breast cancer (non-TNBC). Methods: We performed integrative data analysis combining somatic mutation, epigenomic and gene expression data from The Cancer Genome Atlas (TCGA) to unravel the possible oncogenic interactions between genomic and epigenomic variation in TNBC and non-TNBC. We hypothesized that within breast cancers, there are differences in somatic mutation, DNA methylation and gene expression signatures between TNBC and non-TNBC. We further hypothesized that genomic and epigenomic alterations affect gene regulatory networks and signaling pathways driving the two types of breast cancer. Results: The investigation revealed somatic mutated, epigenomic and gene expression signatures unique to TNBC and non-TNBC and signatures distinguishing the two types of breast cancer. In addition, the investigation revealed molecular networks and signaling pathways enriched for somatic mutations and epigenomic changes unique to each type of breast cancer. The most significant pathways for TNBC were: retinal biosynthesis, BAG2, LXR/RXR, EIF2 and P2Y purigenic receptor signaling pathways. The most significant pathways for non-TNBC were: UVB-induced MAPK, PCP, Apelin endothelial, Endoplasmatic reticulum stress and mechanisms of viral exit from host signaling Pathways. Conclusion: The investigation revealed integrated genomic, epigenomic and gene expression signatures and signing pathways unique to TNBC and non-TNBC, and a gene signature distinguishing the two types of breast cancer. The study demonstrates that integrative analysis of multi-omics data is a powerful approach for unravelling the genomic-epigenomic interaction landscape in TNBC and non-TNBC.

Transcriptome analysis reveals wingless regulates neural development and signaling genes in the region of wing pigmentation of a polka-dotted fruit fly.

How evolutionary novelties have arisen is one of the central questions in evolutionary biology. Preexisting gene regulatory networks or signaling pathways have been shown to be co-opted for building novel traits in several organisms. However, the structure of entire gene regulatory networks and evolutionary events of gene co-option for emergence of a novel trait are poorly understood. In this study, to explore the genetic and molecular bases of the novel wing pigmentation pattern of a polka-dotted fruit fly (Drosophilaguttifera), we performed denovo genome sequencing and transcriptome analyses. As a result, we comprehensively identified the genes associated with the pigmentation pattern. Furthermore, we revealed that 151 of these associated genes were positively or negatively regulated by wingless, a master regulator of wing pigmentation. Genes for neural development, Wnt signaling, Dpp signaling, and effectors (such as enzymes) for melanin pigmentation were included among these 151 genes. None of the known regulatory genes that regulate pigmentation pattern formation in other fruit fly species were included. Our results suggest that the novel pigmentation pattern of a polka-dotted fruit fly might have emerged through multistep co-options of multiple gene regulatory networks, signaling pathways, and effector genes, rather than recruitment of one large gene circuit.