Robust target gene discovery through transcriptome perturbations and genome-wide enhancer predictions in Drosophila uncovers a regulatory basis for sensory specification

Prognostic iron-metabolism signature robustly stratifies single-cell characteristics of hepatocellular carcinoma

Cancer immunotherapy has shown to be a promising method in treating hepatocellular carcinoma (HCC), but suboptimal responses in patients are attributed to cellular and molecular heterogeneity. Iron metabolism-related genes (IRGs) are important in maintaining immune system homeostasis and have the potential to help develop new strategies for HCC treatment. Herein, we constructed and validated the iron-metabolism gene prognostic index (IPX) using univariate Cox proportional hazards regression and LASSO Cox regression analysis, successfully categorizing HCC patients into two groups with distinct survival risks. Then, we performed single-sample gene set enrichment analysis, weighted correlation network analysis, gene ontology enrichment analysis, cellular lineage analysis, and SCENIC analysis to reveal the key determinants underlying the ability of this model based on bulk and single-cell transcriptomic data. We identified several driver transcription factors specifically activated in specific malignant cell sub-populations to contribute to the adverse survival outcomes in the IPX-high subgroup. Within the tumor microenvironment (TME), T cells displayed significant diversity in their cellular characteristics and experienced changes in their developmental paths within distinct clusters identified by IPX. Interestingly, the proportion of Treg cells was increased in the high-risk group compared with the low-risk group. These results suggest that iron-metabolism could be involved in reshaping the TME, thereby disrupting the cell cycle of immune cells. This study utilized IRGs to construct a novel and reliable model, which can be used to assess the prognosis of patients with HCC and further clarify the molecular mechanisms of IRGs in HCC at single-cell resolution.

Cell-type-directed design of synthetic enhancers

Transcriptional enhancers act as docking stations for combinations of transcription factors and thereby regulate spatiotemporal activation of their target genes1. It has been a long-standing goal in the field to decode the regulatory logic of an enhancer and to understand the details of how spatiotemporal gene expression is encoded in an enhancer sequence. Here we show that deep learning models2-6, can be used to efficiently design synthetic, cell-type-specific enhancers, starting from random sequences, and that this optimization process allows detailed tracing of enhancer features at single-nucleotide resolution. We evaluate the function of fully synthetic enhancers to specifically target Kenyon cells or glial cells in the fruit fly brain using transgenic animals. We further exploit enhancer design to create 'dual-code' enhancers that target two cell types and minimal enhancers smaller than 50 base pairs that are fully functional. By examining the state space searches towards local optima, we characterize enhancer codes through the strength, combination and arrangement of transcription factor activator and transcription factor repressor motifs. Finally, we apply the same strategies to successfully design human enhancers, which adhere to enhancer rules similar to those of Drosophila enhancers. Enhancer design guided by deep learning leads to better understanding of how enhancers work and shows that their code can be exploited to manipulate cell states.

Single-Cell RNA-Sequencing Provides Insight into Skeletal Muscle Evolution during the Selection of Muscle Characteristics

Skeletal muscle comprises a large, heterogeneous assortment of cell populations that interact to maintain muscle homeostasis, but little is known about the mechanism that controls myogenic development in response to artificial selection. Different pig (Sus scrofa) breeds exhibit distinct muscle phenotypes resulting from domestication and selective breeding. Using unbiased single-cell transcriptomic sequencing analysis (scRNA-seq), the impact of artificial selection on cell profiles is investigated in neonatal skeletal muscle of pigs. This work provides panoramic muscle-resident cell profiles and identifies novel and breed-specific cells, mapping them on pseudotime trajectories. Artificial selection has elicited significant changes in muscle-resident cell profiles, while conserving signs of generational environmental challenges. These results suggest that fibro-adipogenic progenitors serve as a cellular interaction hub and that specific transcription factors identified here may serve as candidate target regulons for the pursuit of a specific muscle phenotype. Furthermore, a cross-species comparison of humans, mice, and pigs illustrates the conservation and divergence of mammalian muscle ontology. The findings of this study reveal shifts in cellular heterogeneity, novel cell subpopulations, and their interactions that may greatly facilitate the understanding of the mechanism underlying divergent muscle phenotypes arising from artificial selection.

Exploration of the Shared Molecular Mechanisms between COVID-19 and Neurodegenerative Diseases through Bioinformatic Analysis

The COVID-19 pandemic has caused millions of deaths and remains a major public health burden worldwide. Previous studies found that a large number of COVID-19 patients and survivors developed neurological symptoms and might be at high risk of neurodegenerative diseases, such as Alzheimer's disease (AD) and Parkinson's disease (PD). We aimed to explore the shared pathways between COVID-19, AD, and PD by using bioinformatic analysis to reveal potential mechanisms, which may explain the neurological symptoms and degeneration of brain that occur in COVID-19 patients, and to provide early intervention. In this study, gene expression datasets of the frontal cortex were employed to detect common differentially expressed genes (DEGs) of COVID-19, AD, and PD. A total of 52 common DEGs were then examined using functional annotation, protein-protein interaction (PPI) construction, candidate drug identification, and regulatory network analysis. We found that the involvement of the synaptic vesicle cycle and down-regulation of synapses were shared by these three diseases, suggesting that synaptic dysfunction might contribute to the onset and progress of neurodegenerative diseases caused by COVID-19. Five hub genes and one key module were obtained from the PPI network. Moreover, 5 drugs and 42 transcription factors (TFs) were also identified on the datasets. In conclusion, the results of our study provide new insights and directions for follow-up studies of the relationship between COVID-19 and neurodegenerative diseases. The hub genes and potential drugs we identified may provide promising treatment strategies to prevent COVID-19 patients from developing these disorders.

Gene expression variability across cells and species shapes the relationship between renal resident macrophages and infiltrated macrophages

BACKGROUND

Two main subclasses of macrophages are found in almost all solid tissues: embryo-derived resident tissue macrophages and bone marrow-derived infiltrated macrophages. These macrophage subtypes show transcriptional and functional divergence, and the programs that have shaped the evolution of renal macrophages and related signaling pathways remain poorly understood. To clarify these processes, we performed data analysis based on single-cell transcriptional profiling of renal tissue-resident and infiltrated macrophages in human, mouse and rat.

RESULTS

In this study, we (i) characterized the transcriptional divergence among species and (ii) illustrated variability in expression among cells of each subtype and (iii) compared the gene regulation network and (iv) ligand-receptor pairs in human and mouse. Using single-cell transcriptomics, we mapped the promoter architecture during homeostasis.

CONCLUSIONS

Transcriptionally divergent genes, such as the differentially TF-encoding genes expressed in resident and infiltrated macrophages across the three species, vary among cells and include distinct promoter structures. The gene regulatory network in infiltrated macrophages shows comparatively better species-wide consistency than resident macrophages. The conserved transcriptional gene regulatory network in infiltrated macrophages among species is uniquely enriched in pathways related to kinases, and TFs associated with largely conserved regulons among species are uniquely enriched in kinase-related pathways.

Genetic and Molecular Interactions between HΔCT, a Novel Allele of the Notch Antagonist Hairless, and the Histone Chaperone Asf1 in Drosophila melanogaster

Cellular differentiation relies on the highly conserved Notch signaling pathway. Notch activity induces gene expression changes that are highly sensitive to chromatin landscape. We address Notch gene regulation using Drosophila as a model, focusing on the genetic and molecular interactions between the Notch antagonist Hairless and the histone chaperone Asf1. Earlier work implied that Asf1 promotes the silencing of Notch target genes via Hairless (H). Here, we generate a novel HΔCT allele by genome engineering. Phenotypically, HΔCT behaves as a Hairless gain of function allele in several developmental contexts, indicating that the conserved CT domain of H has an attenuator role under native biological contexts. Using several independent methods to assay protein-protein interactions, we define the sequences of the CT domain that are involved in Hairless-Asf1 binding. Based on previous models, where Asf1 promotes Notch repression via Hairless, a loss of Asf1 binding should reduce Hairless repressive activity. However, tissue-specific Asf1 overexpression phenotypes are increased, not rescued, in the HΔCT background. Counterintuitively, Hairless protein binding mitigates the repressive activity of Asf1 in the context of eye development. These findings highlight the complex connections of Notch repressors and chromatin modulators during Notch target-gene regulation and open the avenue for further investigations.

Single-cell RNA-seq analysis of testicular somatic cell development in pigs

Spermatogenesis is the process by which diploid male germ cells propagate and differentiate into haploid flagellated spermatozoa. This highly complicated process is dependent on testicular somatic cells maturation. While the role of these somatic cells in spermatogenesis is relatively well established, knowledge about their development and maturation, particularly at the molecular level, is limited. In this study, we profiled the testicular single-cell transcriptomes of Guanzhong black pigs at the age of 7, 30, 60, 90, and 150 days. Five types of Sertoli cells, five types of Leydig cells, and four types of peritubular myoid cells were identified. Histological analysis revealed the changes in proliferation levels and marker gene expressions, and the prion-like protein gene (PRND) was identified as a novel marker for Sertoli cells. Additionally, integrated analyses of porcine and human datasets revealed similarities between human and pig testicular somatic cells. Overall, the data obtained in this study contribute to the understanding of testicular development in pigs as a model species.

Identification and Validation of a DNA Damage Repair-Related Signature for Diffuse Large B-Cell Lymphoma

BACKGROUND

Diffuse large B-cell lymphoma (DLBCL) is the most common subtype of non-Hodgkin's lymphoma in adults, whose prognostic scoring system remains to be improved. Dysfunction of DNA repair genes is closely associated with the development and prognosis of diffuse large B-cell lymphoma. The aim of this study was to establish and validate a DNA repair-related gene signature associated with the prognosis of DLBCL and to investigate the clinical predictive value of this signature.

METHODS

DLBCL cases were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. One hundred ninety-nine DNA repair-related gene sets were retrieved from the GeneCards database. The LASSO Cox regression was used to generate the DNA repair-related gene signature. Subsequently, the level of immune cell infiltration and the correlation between the gene signature and immune cells were analyzed using the CIBERSORT algorithm. Based on the Genomics of Drug Sensitivity in Cancer (GDSC) database, the relationship between the signature and drug sensitivity was analyzed, and together with the nomogram and gene set variation analysis (GSVA), the value of the signature for clinical application was evaluated.

RESULTS

A total of 14 DNA repair genes were screened out and included in the final risk model. Subgroup analysis of the training and validation cohorts showed that the risk model accurately predicted overall survival of DLBCL patients, with patients in the high-risk group having a worse prognosis than patients in the low-risk group. Subsequently, the risk score was confirmed as an independent prognostic factor by multivariate analysis. Furthermore, by CIBERSORT analysis, we discovered that immune cells, such as regulatory T cells (Tregs), activated memory CD4+ T cells, and gamma delta T cells showed significant differences between the high- and low-risk groups. In addition, we found some interesting associations of our signature with immune checkpoint genes (CD96, TGFBR1, and TIGIT). By analyzing drug sensitivity data in the GDSC database, we were able to identify potential therapeutics for DLBCL patients stratified according to our signature.

CONCLUSIONS

Our study identified and validated a 14-DNA repair-related gene signature for stratification and prognostic prediction of DLBCL patients, which might guide clinical personalization of treatment.

Single-cell RNA sequencing of the Mongolia sheep testis reveals a conserved and divergent transcriptome landscape of mammalian spermatogenesis

Spermatogenesis is a highly coordinated and complex process, and is pivotal for transmitting genetic information between mammalian generations. In this study, we investigated the conservation, differences, and biological functions of homologous genes during spermatogenesis in Mongolia sheep, humans, cynomolgus monkey, and mice using single-cell RNA sequencing technology. We compared X chromosome meiotic inactivation events in Mongolia sheep, humans, cynomolgus monkey, and mice to uncover the concerted activity of X chromosome genes. Subsequently, we focused on the dynamics of gene expression, key biological functions, and signaling pathways at various stages of spermatogenesis in Mongolia sheep and humans. Additionally, the ligand-receptor networks of Mongolia sheep and humans in testicular somatic and germ cells at different developmental stages were mapped to reveal conserved germ cell-soma communication using single-cell resolution. These datasets provided novel information and insights to unravel the molecular regulatory mechanisms of Mongolia sheep spermatogenesis and highlight conservation in gene expression during spermatogenesis between Mongolia sheep and humans, providing a foundation for the establishment of a large mammalian disease model of male infertility.

Decoding gene regulation in the fly brain

The Drosophila brain is a frequently used model in neuroscience. Single-cell transcriptome analysis^1-6, three-dimensional morphological classification⁷ and electron microscopy mapping of the connectome^8,9 have revealed an immense diversity of neuronal and glial cell types that underlie an array of functional and behavioural traits in the fly. The identities of these cell types are controlled by gene regulatory networks (GRNs), involving combinations of transcription factors that bind to genomic enhancers to regulate their target genes. Here, to characterize GRNs at the cell-type level in the fly brain, we profiled the chromatin accessibility of 240,919 single cells spanning 9 developmental timepoints and integrated these data with single-cell transcriptomes. We identify more than 95,000 regulatory regions that are used in different neuronal cell types, of which 70,000 are linked to developmental trajectories involving neurogenesis, reprogramming and maturation. For 40 cell types, uniquely accessible regions were associated with their expressed transcription factors and downstream target genes through a combination of motif discovery, network inference and deep learning, creating enhancer GRNs. The enhancer architectures revealed by DeepFlyBrain lead to a better understanding of neuronal regulatory diversity and can be used to design genetic driver lines for cell types at specific timepoints, facilitating their characterization and manipulation.

Par3 cooperates with Sanpodo for the assembly of Notch clusters following asymmetric division of Drosophila sensory organ precursor cells

In multiple cell lineages, Delta-Notch signalling regulates cell fate decisions owing to unidirectional signalling between daughter cells. In Drosophila pupal sensory organ lineage, Notch regulates the intra-lineage pIIa/pIIb fate decision at cytokinesis. Notch and Delta that localise apically and basally at the pIIa-pIIb interface are expressed at low levels and their residence time at the plasma membrane is in the order of minutes. How Delta can effectively interact with Notch to trigger signalling from a large plasma membrane area remains poorly understood. Here, we report that the signalling interface possesses a unique apico-basal polarity with Par3/Bazooka localising in the form of nano-clusters at the apical and basal level. Notch is preferentially targeted to the pIIa-pIIb interface, where it co-clusters with Bazooka and its cofactor Sanpodo. Clusters whose assembly relies on Bazooka and Sanpodo activities are also positive for Neuralized, the E3 ligase required for Delta activity. We propose that the nano-clusters act as snap buttons at the new pIIa-pIIb interface to allow efficient intra-lineage signalling.

Evaluating the Reproducibility of Single-Cell Gene Regulatory Network Inference Algorithms

Networks are powerful tools to represent and investigate biological systems. The development of algorithms inferring regulatory interactions from functional genomics data has been an active area of research. With the advent of single-cell RNA-seq data (scRNA-seq), numerous methods specifically designed to take advantage of single-cell datasets have been proposed. However, published benchmarks on single-cell network inference are mostly based on simulated data. Once applied to real data, these benchmarks take into account only a small set of genes and only compare the inferred networks with an imposed ground-truth. Here, we benchmark six single-cell network inference methods based on their reproducibility, i.e., their ability to infer similar networks when applied to two independent datasets for the same biological condition. We tested each of these methods on real data from three biological conditions: human retina, T-cells in colorectal cancer, and human hematopoiesis. Once taking into account networks with up to 100,000 links, GENIE3 results to be the most reproducible algorithm and, together with GRNBoost2, show higher intersection with ground-truth biological interactions. These results are independent from the single-cell sequencing platform, the cell type annotation system and the number of cells constituting the dataset. Finally, GRNBoost2 and CLR show more reproducible performance once a more stringent thresholding is applied to the networks (1,000–100 links). In order to ensure the reproducibility and ease extensions of this benchmark study, we implemented all the analyses in scNET, a Jupyter notebook available at https://github.com/ComputationalSystemsBiology/scNET.

A binding protein regulates myosin-7a dimerization and actin bundle assembly

Myosin-7a, despite being monomeric in isolation, plays roles in organizing actin-based cell protrusions such as filopodia, microvilli and stereocilia, as well as transporting cargoes within them. Here, we identify a binding protein for Drosophila myosin-7a termed M7BP, and describe how M7BP assembles myosin-7a into a motile complex that enables cargo translocation and actin cytoskeletal remodeling. M7BP binds to the autoinhibitory tail of myosin-7a, extending the molecule and activating its ATPase activity. Single-molecule reconstitution show that M7BP enables robust motility by complexing with myosin-7a as 2:2 translocation dimers in an actin-regulated manner. Meanwhile, M7BP tethers actin, enhancing complex’s processivity and driving actin-filament alignment during processive runs. Finally, we show that myosin-7a-M7BP complex assembles actin bundles and filopodia-like protrusions while migrating along them in living cells. Together, these findings provide insights into the mechanisms by which myosin-7a functions in actin protrusions.

Myosin-7a is found in actin bundles, microvilli and stereocilia, and plays conserved roles in hearing and vision. Here the authors identify M7BP, a myosin-7a binding protein that activates and dimerizes myosin-7a, enabling cargo transport and assembly of actin bundles and filopodia-like protrusions

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 Clathrin adaptor AP-1 and Stratum act in parallel pathways to control Notch activation in Drosophila sensory organ precursors cells

Drosophila sensory organ precursors divide asymmetrically to generate pIIa/pIIb cells, the identity of which relies on activation of Notch at cytokinesis. Although Notch is present apically and basally relative to the midbody at the pIIa-pIIb interface, the basal pool of Notch is reported to be the main contributor for Notch activation in the pIIa cell. Intra-lineage signalling requires appropriate apico-basal targeting of Notch, its ligand Delta and its trafficking partner Sanpodo. We have previously reported that AP-1 and Stratum regulate the trafficking of Notch and Sanpodo from the trans-Golgi network to the basolateral membrane. Loss of AP-1 or Stratum caused mild Notch gain-of-function phenotypes. Here, we report that their concomitant loss results in a penetrant Notch gain-of-function phenotype, indicating that they control parallel pathways. Although unequal partitioning of cell fate determinants and cell polarity were unaffected, we observed increased amounts of signalling-competent Notch as well as Delta and Sanpodo at the apical pIIa-pIIb interface, at the expense of the basal pool of Notch. We propose that AP-1 and Stratum operate in parallel pathways to localize Notch and control where receptor activation takes place.

Single-Cell RNA Sequencing of the Cynomolgus Macaque Testis Reveals Conserved Transcriptional Profiles during Mammalian Spermatogenesis

Spermatogenesis is highly orchestrated and involves the differentiation of diploid spermatogonia into haploid sperm. The process is driven by spermatogonial stem cells (SSCs). SSCs undergo mitotic self-renewal, whereas sub-populations undergo differentiation and later gain competence to initiate meiosis. Here, we describe a high-resolution single-cell RNA-seq atlas of cells derived from Cynomolgus macaque testis. We identify gene signatures that define spermatogonial populations and explore self-renewal versus differentiation dynamics. We detail transcriptional changes occurring over the entire process of spermatogenesis and highlight the concerted activity of DNA damage response (DDR) pathway genes, which have dual roles in maintaining genomic integrity and effecting meiotic sex chromosome inactivation (MSCI). We show remarkable similarities and differences in gene expression during spermatogenesis with two other eutherian mammals, i.e., mouse and humans. Sex chromosome expression in the male germline in all three species demonstrates conserved features of MSCI but divergent multicopy and ampliconic gene content.

Atonal bHLH transcription factor 1 is an important factor for maintaining the balance of cell proliferation and differentiation in tumorigenesis

Establishing the link between cellular processes and oncogenesis may aid the elucidation of targeted and effective therapies against tumor cell proliferation and metastasis. Previous studies have investigated the mechanisms involved in maintaining the balance between cell proliferation, differentiation and migration. There is increased interest in determining the conditions that allow cancer stem cells to differentiate as well as the identification of molecules that may serve as novel drug targets. Furthermore, the study of various genes, including transcription factors, which serve a crucial role in cellular processes, may present a promising direction for future therapy. The present review described the role of the transcription factor atonal bHLH transcription factor 1 (ATOH1) in signaling pathways in tumorigenesis, particularly in cerebellar tumor medulloblastoma and colorectal cancer, where ATOH1 serves as an oncogene or tumor suppressor, respectively. Additionally, the present review summarized the associated therapeutic interventions for these two types of tumors and discussed novel clinical targets and approaches.

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

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

Gene regulatory networks during the development of the Drosophila visual system

The Drosophila visual system integrates input from 800 ommatidia and extracts different features in stereotypically connected optic ganglia. The development of the Drosophila visual system is controlled by gene regulatory networks that control the number of precursor cells, generate neuronal diversity by integrating spatial and temporal information, coordinate the timing of retinal and optic lobe cell differentiation, and determine distinct synaptic targets of each cell type. In this chapter, we describe the known gene regulatory networks involved in the development of the different parts of the visual system and explore general components in these gene networks. Finally, we discuss the advantages of the fly visual system as a model for gene regulatory network discovery in the era of single-cell transcriptomics.

TMT-Based Quantitative Proteomic Profiling of Overwintering Lissorhoptrus oryzophilus

Adaptations to low temperature play a critical role in restricting the geographical distribution of insects. Decreasing day lengths and temperatures trigger seasonal cold adaptations in insects. These adaptions include changes in expression at the miRNA, mRNA and protein levels. The rice water weevil (RWW), Lissorhoptrus oryzophilus, introduced from the Mississippi River, is a globally invasive pest of wetland rice that can survive at the northern border of China. To investigate the changes in expression at the protein level in overwintering female RWW adults, 6-plex tandem mass tags (TMTs) were used in overwintering and summer adults. By using a proteome database available for Curculionidae, 1077 proteins were quantified, 183 of which differed significantly between the overwintering and summer samples. To further understand these differentially expressed proteins (DEPs), bioinformatics analyses such as gene ontology (GO) enrichment analyses were performed. DEPs associated with the terms binding, structural molecule activity, catalytic activity, multicellular organismal process, extracellular region, chitin binding, metabolic process, intracellular part and organic cyclic compound binding were altered by selection during winter. The changes in the expression of these proteins suggest that the proteins are important for RWW survival in winter.

Altering the Temporal Regulation of One Transcription Factor Drives Evolutionary Trade-Offs between Head Sensory Organs

Size trade-offs of visual versus olfactory organs is a pervasive feature of animal evolution. This could result from genetic or functional constraints. We demonstrate that head sensory organ size trade-offs in Drosophila are genetically encoded and arise through differential subdivision of the head primordium into visual versus non-visual fields. We discover that changes in the temporal regulation of the highly conserved eyeless/Pax6 gene expression during development is a conserved mechanism for sensory trade-offs within and between Drosophila species. We identify a natural single nucleotide polymorphism in the cis-regulatory region of eyeless in a binding site of its repressor Cut that is sufficient to alter its temporal regulation and eye size. Because eyeless/Pax6 is a conserved regulator of head sensory placode subdivision, we propose that its temporal regulation is key to define the relative size of head sensory organs.

Retrograde BMP signaling activates neuronal gene expression through widespread deployment of a conserved BMP-responsive cis-regulatory activation element

Retrograde Bone Morphogenetic Protein (BMP) signaling in neurons is essential for the differentiation and synaptic function of many neuronal subtypes. BMP signaling regulates these processes via Smad transcription factor activity, yet the scope and nature of Smad-dependent gene regulation in neurons are mostly unknown. Here, we applied a computational approach to predict Smad-binding cis-regulatory BMP-Activating Elements (BMP-AEs) in Drosophila, followed by transgenic in vivo reporter analysis to test their neuronal subtype enhancer activity in the larval central nervous system (CNS). We identified 34 BMP-AE-containing genomic fragments that are responsive to BMP signaling in neurons, and showed that the embedded BMP-AEs are required for this activity. RNA-seq analysis identified BMP-responsive genes in the CNS and revealed that BMP-AEs selectively enrich near BMP-activated genes. These data suggest that functional BMP-AEs control nearby BMP-activated genes, which we validated experimentally. Finally, we demonstrated that the BMP-AE motif mediates a conserved Smad-responsive function in the Drosophila and vertebrate CNS. Our results provide evidence that BMP signaling controls neuronal function by directly coordinating the expression of a battery of genes through widespread deployment of a conserved Smad-responsive cis-regulatory motif.

Identification of key gene modules and transcription factors for human osteoarthritis by weighted gene co-expression network analysis

Osteoarthritis (OA) is one of the most prevalent causes of joint disease. However, the pathological mechanisms of OA have remained to be completely elucidated, and further investigation into the underlying mechanisms of OA development and the identification of novel therapeutic targets are urgently required. In the present study, the dataset GSE114007 was downloaded from the Gene Expression Omnibus database. Based on weighted gene co-expression network analysis (WGCNA) and the identification of differentially expressed genes (DEGs), the microarray data were further analyzed to identify hub genes, key transcription factors (TFs) and pivotal signaling pathways involved in the pathogenesis of OA. A total of 1,898 genes were identified to be differentially expressed between OA samples and normal samples. Based on WGCNA, the present study identified 5 hub modules closely associated with OA, and the potential key TFs for hub modules were further explored based on CisTargetX. The results demonstrated that B-Cell Lymphoma 6, Myelin Gene Expression Factor 2, Activating Transcription Factor 3, CCAAT Enhancer Binding Protein γ, Nuclear Factor Interleukin-3-Regulated, FOS Like Antigen-2, FOS-Like Antigen-1, Fos Proto-Oncogene, JunD Proto-Oncogene, Transcription Factor CP2 Like 1, RELA proto-oncogene NF-kB subunit, SRY-box transcription factor 3, V-Ets Avian Erythroblastosis Virus E26 Oncogene Homolog 2, Interferon Regulatory Factor 4 and REL proto-oncogene, NF-kB subunit were the potential key TFs. In addition, osteoclast differentiation, FoxO, MAPK and PI3K/Akt signaling pathways were revealed to be imperative for the pathogenesis of OA, as these 4 pivotal signaling pathways were observed to be tightly linked through 4 key TFs Fos Proto-Oncogene, JUN, JunD Proto-Oncogene and MYC, and 4 DEGs Vascular Endothelial Growth Factor A, Growth Arrest and DNA Damage Inducible α, Growth Arrest and DNA Damage Inducible β and Cyclin D1. The present study identified a set of potential key genes and signaling pathways, and provided an important opportunity to advance the current understanding of OA.

Systems Biology Reveals NR2F6 and TGFB1 as Key Regulators of Feed Efficiency in Beef Cattle

Systems biology approaches are used as strategy to uncover tissue-specific perturbations and regulatory genes related to complex phenotypes. We applied this approach to study feed efficiency (FE) in beef cattle, an important trait both economically and environmentally. Poly-A selected RNA of five tissues (adrenal gland, hypothalamus, liver, skeletal muscle and pituitary) of eighteen young bulls, selected for high and low FE, were sequenced (Illumina HiSeq 2500, 100 bp, pared-end). From the 17,354 expressed genes considering all tissues, 1,335 were prioritized by five selection categories (differentially expressed, harboring SNPs associated with FE, tissue-specific, secreted in plasma and key regulators) and used for network construction. NR2F6 and TGFB1 were identified and validated by motif discovery as key regulators of hepatic inflammatory response and muscle tissue development, respectively, two biological processes demonstrated to be associated with FE. Moreover, we indicated potential biomarkers of FE, which are related to hormonal control of metabolism and sexual maturity. By using robust methodologies and validation strategies, we confirmed the main biological processes related to FE in Bos indicus and indicated candidate genes as regulators or biomarkers of superior animals.

The transcription factor Grainy head primes epithelial enhancers for spatiotemporal activation by displacing nucleosomes

Transcriptional enhancers function as docking platforms for combinations of transcription factors (TFs) to control gene expression. How enhancer sequences determine nucleosome occupancy, TF recruitment and transcriptional activation in vivo remains unclear. Using ATAC-seq across a panel of Drosophila inbred strains, we found that SNPs affecting binding sites of the TF Grainy head (Grh) causally determine the accessibility of epithelial enhancers. We show that deletion and ectopic expression of Grh cause loss and gain of DNA accessibility, respectively. However, although Grh binding is necessary for enhancer accessibility, it is insufficient to activate enhancers. Finally, we show that human Grh homologs-GRHL1, GRHL2 and GRHL3-function similarly. We conclude that Grh binding is necessary and sufficient for the opening of epithelial enhancers but not for their activation. Our data support a model positing that complex spatiotemporal expression patterns are controlled by regulatory hierarchies in which pioneer factors, such as Grh, establish tissue-specific accessible chromatin landscapes upon which other factors can act.

Dynamic genome wide expression profiling of Drosophila head development reveals a novel role of Hunchback in retinal glia cell development and blood-brain barrier integrity

Drosophila melanogaster head development represents a valuable process to study the developmental control of various organs, such as the antennae, the dorsal ocelli and the compound eyes from a common precursor, the eye-antennal imaginal disc. While the gene regulatory network underlying compound eye development has been extensively studied, the key transcription factors regulating the formation of other head structures from the same imaginal disc are largely unknown. We obtained the developmental transcriptome of the eye-antennal discs covering late patterning processes at the late 2^nd larval instar stage to the onset and progression of differentiation at the end of larval development. We revealed the expression profiles of all genes expressed during eye-antennal disc development and we determined temporally co-expressed genes by hierarchical clustering. Since co-expressed genes may be regulated by common transcriptional regulators, we combined our transcriptome dataset with publicly available ChIP-seq data to identify central transcription factors that co-regulate genes during head development. Besides the identification of already known and well-described transcription factors, we show that the transcription factor Hunchback (Hb) regulates a significant number of genes that are expressed during late differentiation stages. We confirm that hb is expressed in two polyploid subperineurial glia cells (carpet cells) and a thorough functional analysis shows that loss of Hb function results in a loss of carpet cells in the eye-antennal disc. Additionally, we provide for the first time functional data indicating that carpet cells are an integral part of the blood-brain barrier. Eventually, we combined our expression data with a de novo Hb motif search to reveal stage specific putative target genes of which we find a significant number indeed expressed in carpet cells.

The development of different cell types must be tightly coordinated, and the eye-antennal imaginal discs of Drosophila melanogaster represent an excellent model to study the molecular mechanisms underlying this coordination. These imaginal discs contain the anlagen of nearly all adult head structures, such as the antennae, the head cuticle, the ocelli and the compound eyes. While large scale screens have been performed to unravel the gene regulatory network underlying compound eye development, a comprehensive understanding of genome wide expression dynamics throughout head development is still missing to date. We studied the genome wide gene expression dynamics during eye-antennal disc development in D. melanogaster to identify new central regulators of the underlying gene regulatory network. Expression based gene clustering and transcription factor motif enrichment analyses revealed a central regulatory role of the transcription factor Hunchback (Hb). We confirmed that hb is expressed in two polyploid retinal subperineurial glia cells (carpet cells). Our functional analysis shows that Hb is necessary for carpet cell development and we show for the first time that the carpet cells are an integral part of the blood-brain barrier.

SCENIC: single-cell regulatory network inference and clustering

Although single-cell RNA-seq is revolutionizing biology, data interpretation remains a challenge. We present SCENIC for the simultaneous reconstruction of gene regulatory networks and identification of cell states. We apply SCENIC to a compendium of single-cell data from tumors and brain, and demonstrate that the genomic regulatory code can be exploited to guide the identification of transcription factors and cell states. SCENIC provides critical biological insights into the mechanisms driving cellular heterogeneity.

Conserved Transcription Factors Steer Growth-Related Genomic Programs in Daphnia

Ecological genomics aims to understand the functional association between environmental gradients and the genes underlying adaptive traits. Many genes that are identified by genome-wide screening in ecologically relevant species lack functional annotations. Although gene functions can be inferred from sequence homology, such approaches have limited power. Here, we introduce ecological regulatory genomics by presenting an ontology-free gene prioritization method. Specifically, our method combines transcriptome profiling with high-throughput cis-regulatory sequence analysis in the water fleas Daphnia pulex and Daphnia magna. It screens coexpressed genes for overrepresented DNA motifs that serve as transcription factor binding sites, thereby providing insight into conserved transcription factors and gene regulatory networks shaping the expression profile. We first validated our method, called Daphnia-cisTarget, on a D. pulex heat shock data set, which revealed a network driven by the heat shock factor. Next, we performed RNA-Seq in D. magna exposed to the cyanobacterium Microcystis aeruginosa. Daphnia-cisTarget identified coregulated gene networks that associate with the moulting cycle and potentially regulate life history changes in growth rate and age at maturity. These networks are predicted to be regulated by evolutionary conserved transcription factors such as the homologues of Drosophila Shavenbaby and Grainyhead, nuclear receptors, and a GATA family member. In conclusion, our approach allows prioritising candidate genes in Daphnia without bias towards prior knowledge about functional gene annotation and represents an important step towards exploring the molecular mechanisms of ecological responses in organisms with poorly annotated genomes.

Identification of novel direct targets of Drosophila Sine oculis and Eyes absent by integration of genome-wide data sets

Drosophila eye development is a complex process that involves many transcription factors (TFs) and interactions with their cofactors and targets. The TF Sine oculis (So) and its cofactor Eyes absent (Eya) are highly conserved and are both necessary and sufficient for eye development. Despite their many important roles during development, the direct targets of So are still largely unknown. Therefore the So-dependent regulatory network governing eye determination and differentiation is poorly understood. In this study, we intersected gene expression profiles of so or eya mutant eye tissue prepared from three different developmental stages and identified 1731 differentially expressed genes across the Drosophila genome. A combination of co-expression analyses and motif discovery identified a set of twelve putative direct So targets, including three known and nine novel targets. We also used our previous So ChIP-seq data to assess motif predictions for So and identified a canonical So binding motif. Finally, we performed in vivo enhancer reporter assays to test predicted enhancers from six candidate target genes and find that at least one enhancer from each gene is expressed in the developing eye disc and that their expression patterns overlap with that of So. We furthermore confirmed that the expression level of predicted direct So targets, for which antibodies are available, are reduced in so or eya post-mitotic knockout eye discs. In summary, we expand the set of putative So targets and show for the first time that the combined use of expression profiling of so with its cofactor eya is an effective method to identify novel So targets. Moreover, since So is highly conserved throughout the metazoa, our results provide the basis for future functional studies in a wide variety of organisms.

The Capacity to Act in Trans Varies Among Drosophila Enhancers

The interphase nucleus is organized such that genomic segments interact in cis, on the same chromosome, and in trans, between different chromosomes. In Drosophila and other Dipterans, extensive interactions are observed between homologous chromosomes, which can permit enhancers and promoters to communicate in trans Enhancer action in trans has been observed for a handful of genes in Drosophila, but it is as yet unclear whether this is a general property of all enhancers or specific to a few. Here, we test a collection of well-characterized enhancers for the capacity to act in trans Specifically, we tested 18 enhancers that are active in either the eye or wing disc of third instar Drosophila larvae and, using two different assays, found evidence that each enhancer can act in trans However, the degree to which trans-action was supported varied greatly between enhancers. Quantitative analysis of enhancer activity supports a model wherein an enhancer's strength of transcriptional activation is a major determinant of its ability to act in trans, but that additional factors may also contribute to an enhancer's trans-activity. In sum, our data suggest that a capacity to activate a promoter on a paired chromosome is common among Drosophila enhancers.

The transcription factor Glass links eye field specification with photoreceptor differentiation in Drosophila

Eye development requires an evolutionarily conserved group of transcription factors, termed the retinal determination network (RDN). However, little is known about the molecular mechanism by which the RDN instructs cells to differentiate into photoreceptors. We show that photoreceptor cell identity in Drosophila is critically regulated by the transcription factor Glass, which is primarily expressed in photoreceptors and whose role in this process was previously unknown. Glass is both required and sufficient for the expression of phototransduction proteins. Our results demonstrate that the RDN member Sine oculis directly activates glass expression, and that Glass activates the expression of the transcription factors Hazy and Otd. We identified hazy as a direct target of Glass. Induced expression of Hazy in the retina partially rescues the glass mutant phenotype. Together, our results provide a transcriptional link between eye field specification and photoreceptor differentiation in Drosophila, placing Glass at a central position in this developmental process.

Regulation of transcription factors on sexual dimorphism of fig wasps

Fig wasps exhibit extreme intraspecific morphological divergence in the wings, compound eyes, antennae, body color, and size. Corresponding to this, behaviors and lifestyles between two sexes are also different: females can emerge from fig and fly to other fig tree to oviposit and pollinate, while males live inside fig for all their lifetime. Genetic regulation may drive these extreme intraspecific morphological and behavioral divergence. Transcription factors (TFs) involved in morphological development and physiological activity may exhibit sex-specific expressions. Herein, we detect 865 TFs by using genomic and transcriptomic data of the fig wasp Ceratosolen solmsi. Analyses of transcriptomic data indicated that up-regulated TFs in females show significant enrichment in development of the wing, eye and antenna in all stages, from larva to adult. Meanwhile, TFs related to the development of a variety of organs display sex-specific patterns of expression in the adults and these may contribute significantly to their sexual dimorphism. In addition, up-regulated TFs in adult males exhibit enrichment in genitalia development and circadian rhythm, which correspond with mating and protandry. This finding is consistent with their sex-specific behaviors. In conclusion, our results strongly indicate that TFs play important roles in the sexual dimorphism of fig wasps.

Identification of Lineage-Specific Cis-Regulatory Modules Associated with Variation in Transcription Factor Binding and Chromatin Activity Using Ornstein-Uhlenbeck Models

Scoring the impact of noncoding variation on the function of cis-regulatory regions, on their chromatin state, and on the qualitative and quantitative expression levels of target genes is a fundamental problem in evolutionary genomics. A particular challenge is how to model the divergence of quantitative traits and to identify relationships between the changes across the different levels of the genome, the chromatin activity landscape, and the transcriptome. Here, we examine the use of the Ornstein-Uhlenbeck (OU) model to infer selection at the level of predicted cis-regulatory modules (CRMs), and link these with changes in transcription factor binding and chromatin activity. Using publicly available cross-species ChIP-Seq and STARR-Seq data we show how OU can be applied genome-wide to identify candidate transcription factors for which binding site and CRM turnover is correlated with changes in regulatory activity. Next, we profile open chromatin in the developing eye across three Drosophila species. We identify the recognition motifs of the chromatin remodelers, Trithorax-like and Grainyhead as mostly correlating with species-specific changes in open chromatin. In conclusion, we show in this study that CRM scores can be used as quantitative traits and that motif discovery approaches can be extended towards more complex models of divergence.

Transcriptional selectors, masters, and combinatorial codes: regulatory principles of neural subtype specification

The broad range of tissue and cellular diversity of animals is generated to a large extent by the hierarchical deployment of sequence-specific transcription factors and co-factors (collectively referred to as TF's herein) during development. Our understanding of these developmental processes has been facilitated by the recognition that the activities of many TF's can be meaningfully described by a few functional categories that usefully convey a sense for how the TF's function, and also provides a sense for the regulatory organization of the developmental processes in which they participate. Here, we draw on examples from studies in Caenorhabditis elegans, Drosophila melanogaster, and vertebrates to discuss how the terms spatial selector, temporal selector, tissue/cell type selector, terminal selector and combinatorial code may be usefully applied to categorize the activities of TF's at critical steps of nervous system construction. While we believe that these functional categories are useful for understanding the organizational principles by which TF's direct nervous system construction, we however caution against the assumption that a TF's function can be solely or fully defined by any single functional category. Indeed, most TF's play diverse roles within different functional categories, and their roles can blur the lines we draw between these categories. Regardless, it is our belief that the concepts discussed here are helpful in clarifying the regulatory complexities of nervous system development, and hope they prove useful when interpreting mutant phenotypes, designing future experiments, and programming specific neuronal cell types for use in therapies. WIREs Dev Biol 2015, 4:505–528. doi: 10.1002/wdev.191

For further resources related to this article, please visit the WIREs website.

Identifying transcriptional cis-regulatory modules in animal genomes

Gene expression is regulated through the activity of transcription factors (TFs) and chromatin-modifying proteins acting on specific DNA sequences, referred to as cis-regulatory elements. These include promoters, located at the transcription initiation sites of genes, and a variety of distal cis-regulatory modules (CRMs), the most common of which are transcriptional enhancers. Because regulated gene expression is fundamental to cell differentiation and acquisition of new cell fates, identifying, characterizing, and understanding the mechanisms of action of CRMs is critical for understanding development. CRM discovery has historically been challenging, as CRMs can be located far from the genes they regulate, have few readily identifiable sequence characteristics, and for many years were not amenable to high-throughput discovery methods. However, the recent availability of complete genome sequences and the development of next-generation sequencing methods have led to an explosion of both computational and empirical methods for CRM discovery in model and nonmodel organisms alike. Experimentally, CRMs can be identified through chromatin immunoprecipitation directed against TFs or histone post-translational modifications, identification of nucleosome-depleted 'open' chromatin regions, or sequencing-based high-throughput functional screening. Computational methods include comparative genomics, clustering of known or predicted TF-binding sites, and supervised machine-learning approaches trained on known CRMs. All of these methods have proven effective for CRM discovery, but each has its own considerations and limitations, and each is subject to a greater or lesser number of false-positive identifications. Experimental confirmation of predictions is essential, although shortcomings in current methods suggest that additional means of validation need to be developed. For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.

Hindsight regulates photoreceptor axon targeting through transcriptional control of jitterbug/Filamin and multiple genes involved in axon guidance in Drosophila

During axon targeting, a stereotyped pattern of connectivity is achieved by the integration of intrinsic genetic programs and the response to extrinsic long and short-range directional cues. How this coordination occurs is the subject of intense study. Transcription factors play a central role due to their ability to regulate the expression of multiple genes required to sense and respond to these cues during development. Here we show that the transcription factor HNT regulates layer-specific photoreceptor axon targeting in Drosophila through transcriptional control of jbug/Filamin and multiple genes involved in axon guidance and cytoskeleton organization.Using a microarray analysis we identified 235 genes whose expression levels were changed by HNT overexpression in the eye primordia. We analyzed nine candidate genes involved in cytoskeleton regulation and axon guidance, six of which displayed significantly altered gene expression levels in hnt mutant retinas. Functional analysis confirmed the role of OTK/PTK7 in photoreceptor axon targeting and uncovered Tiggrin, an integrin ligand, and Jbug/Filamin, a conserved actin- binding protein, as new factors that participate of photoreceptor axon targeting. Moreover, we provided in silico and molecular evidence that supports jbug/Filamin as a direct transcriptional target of HNT and that HNT acts partially through Jbug/Filamin in vivo to regulate axon guidance. Our work broadens the understanding of how HNT regulates the coordinated expression of a group of genes to achieve the correct connectivity pattern in the Drosophila visual system. © 2015 Wiley Periodicals, Inc. Develop Neurobiol 75: 1018-1032, 2015.

Neural precursor-specific expression of multiple Drosophila genes is driven by dual enhancer modules with overlapping function

Transcriptional cis-regulatory modules (CRMs), or enhancers, are responsible for directing gene expression in specific territories and cell types during development. In some instances, the same gene may be served by two or more enhancers with similar specificities. Here we show that the utilization of dual, or "shadow", enhancers is a common feature of genes that are active specifically in neural precursor (NP) cells in Drosophila. By genome-wide computational discovery of statistically significant clusters of binding motifs for both proneural activator (P) proteins and basic helix-loop-helix (bHLH) repressor (R) factors (a "P+R" regulatory code), we have identified NP-specific enhancer modules associated with multiple genes expressed in this cell type. These CRMs are distinct from those previously identified for the corresponding gene, establishing the existence of a dual-enhancer arrangement in which both modules reside close to the gene they serve. Using wild-type and mutant reporter gene constructs in vivo, we show that P sites in these modules mediate activation by proneural factors in "proneural cluster" territories, whereas R sites mediate repression by bHLH repressors, which serves to restrict expression specifically to NP cells. To our knowledge, our results identify the first direct targets of these bHLH repressors. Finally, using genomic rescue constructs for neuralized (neur), we demonstrate that each of the gene's two NP-specific enhancers is sufficient to rescue neur function in the lateral inhibition process by which adult sensory organ precursor (SOP) cells are specified, but that deletion of both enhancers results in failure of this event.

A fluorescent tagging approach in Drosophila reveals late endosomal trafficking of Notch and Sanpodo

Signaling and endocytosis are highly integrated processes that regulate cell fate. In the Drosophila melanogaster sensory bristle lineages, Numb inhibits the recycling of Notch and its trafficking partner Sanpodo (Spdo) to regulate cell fate after asymmetric cell division. In this paper, we have used a dual GFP/Cherry tagging approach to study the distribution and endosomal sorting of Notch and Spdo in living pupae. The specific properties of GFP, i.e., quenching at low pH, and Cherry, i.e., slow maturation time, revealed distinct pools of Notch and Spdo: cargoes exhibiting high GFP/low Cherry fluorescence intensities localized mostly at the plasma membrane and early/sorting endosomes, whereas low GFP/high Cherry cargoes accumulated in late acidic endosomes. These properties were used to show that Spdo is sorted toward late endosomes in a Numb-dependent manner. This dual-tagging approach should be generally applicable to study the trafficking dynamics of membrane proteins in living cells and tissues.

Current directions and future perspectives from the third Nematostella research conference

The third Nematostella vectensis Research Conference took place in December 2013 in Eilat, Israel, as a satellite to the 8th International Conference on Coelenterate Biology. The starlet sea anemone, N. vectensis, has emerged as a powerful cnidarian model, in large part due to the extensive genomic and transcriptomic resources and molecular approaches that are becoming available for Nematostella, which were the focus of several presentations. In addition, research was presented highlighting the broader utility of this species for studies of development, circadian rhythms, signal transduction, and gene-environment interactions.

Nemo promotes Notch-mediated lateral inhibition downstream of proneural factors

During neurogenesis, conserved tissue-specific proneural factors establish a cell's competence to take on neural fate from within a field of unspecified cells. Proneural genes encode basic helix-loop-helix transcription factors that promote the expression of 'core' and subtype-specific target genes. Target genes include both pan-neuronal genes and genes that aid in the process of refinement, known as lateral inhibition. In this process, proneural gene expression is increased in the neural progenitor while simultaneously down-regulated in the surrounding cells, in a Notch signalling-dependent manner. Here, we identify nemo (nmo) as a target of members of both Drosophila Atonal and Achaete-Scute proneural factor families and find that mammalian proneural homologs induce Nemo-like-kinase (Nlk) expression in cell culture. We find that nmo loss of function leads to reduced expression of Notch targets and to perturbations in Notch-mediated lateral inhibition. Furthermore, Notch hyperactivity can compensate for nmo loss in the Drosophila eye. Thus nmo promotes Notch-mediated lateral inhibition downstream of proneural factors during neurogenesis.

iRegulon: from a gene list to a gene regulatory network using large motif and track collections

Identifying master regulators of biological processes and mapping their downstream gene networks are key challenges in systems biology. We developed a computational method, called iRegulon, to reverse-engineer the transcriptional regulatory network underlying a co-expressed gene set using cis-regulatory sequence analysis. iRegulon implements a genome-wide ranking-and-recovery approach to detect enriched transcription factor motifs and their optimal sets of direct targets. We increase the accuracy of network inference by using very large motif collections of up to ten thousand position weight matrices collected from various species, and linking these to candidate human TFs via a motif2TF procedure. We validate iRegulon on gene sets derived from ENCODE ChIP-seq data with increasing levels of noise, and we compare iRegulon with existing motif discovery methods. Next, we use iRegulon on more challenging types of gene lists, including microRNA target sets, protein-protein interaction networks, and genetic perturbation data. In particular, we over-activate p53 in breast cancer cells, followed by RNA-seq and ChIP-seq, and could identify an extensive up-regulated network controlled directly by p53. Similarly we map a repressive network with no indication of direct p53 regulation but rather an indirect effect via E2F and NFY. Finally, we generalize our computational framework to include regulatory tracks such as ChIP-seq data and show how motif and track discovery can be combined to map functional regulatory interactions among co-expressed genes. iRegulon is available as a Cytoscape plugin from http://iregulon.aertslab.org.

Gene regulatory networks control developmental, homeostatic, and disease processes by governing precise levels and spatio-temporal patterns of gene expression. Determining their topology can provide mechanistic insight into these processes. Gene regulatory networks consist of interactions between transcription factors and their direct target genes. Each regulatory interaction represents the binding of the transcription factor to a specific DNA binding site near its target gene. Here we present a computational method, called iRegulon, to identify master regulators and direct target genes in a human gene signature, i.e. a set of co-expressed genes. iRegulon relies on the analysis of the regulatory sequences around each gene in the gene set to detect enriched TF motifs or ChIP-seq peaks, using databases of nearly 10.000 TF motifs and 1000 ChIP-seq data sets or “tracks”. Next, it associates enriched motifs and tracks with candidate transcription factors and determines the optimal subset of direct target genes. We validate iRegulon on ENCODE data, and use it in combination with RNA-seq and ChIP-seq data to map a p53 downstream network with new predicted co-factors and targets. iRegulon is available as a Cytoscape plugin, supporting human, mouse, and Drosophila genes, and provides access to hundreds of cancer-related TF-target subnetworks or “regulons”.

Identification of cis-regulatory modules encoding temporal dynamics during development

Background

Developmental transcriptional regulatory networks are circuits of transcription factors (TFs) and cis-acting DNA elements (Cis Regulatory Modules, CRMs) that dynamically control expression of downstream genes. Comprehensive knowledge of these networks is an essential step towards our understanding of developmental processes. However, this knowledge is mostly based on genome-wide mapping of transcription factor binding sites, and therefore requires prior knowledge regarding the TFs involved in the network.

Results

Focusing on how temporal control of gene expression is integrated within a developmental network, we applied an in silico approach to discover regulatory motifs and CRMs of co-expressed genes, with no prior knowledge about the involved TFs. Our aim was to identify regulatory motifs and potential trans-acting factors which regulate the temporal expression of co-expressed gene sets during a particular process of organogenesis, namely adult heart formation in Drosophila. Starting from whole genome tissue specific expression dynamics, we used an in silico method, cisTargetX, to predict TF binding motifs and CRMs. Potential Nuclear Receptor (NR) binding motifs were predicted to control the temporal expression profile of a gene set with increased expression levels during mid metamorphosis. The predicted CRMs and NR motifs were validated in vivo by reporter gene essays. In addition, we provide evidence that three NRs modulate CRM activity and behave as temporal regulators of target enhancers.

Conclusions

Our approach was successful in identifying CRMs and potential TFs acting on the temporal regulation of target genes. In addition, our results suggest a modular architecture of the regulatory machinery, in which the temporal and spatial regulation can be uncoupled and encoded by distinct CRMs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-534) contains supplementary material, which is available to authorized users.

The Little Fly that Could: Wizardry and Artistry of Drosophila Genomics

For more than 100 years now, the fruit fly Drosophila melanogaster has been at the forefront of our endeavors to unlock the secrets of the genome. From the pioneering studies of chromosomes and heredity by Morgan and his colleagues, to the generation of fly models for human disease, Drosophila research has been at the forefront of genetics and genomics. We present a broad overview of some of the most powerful genomics tools that keep Drosophila research at the cutting edge of modern biomedical research.

Dissection of thousands of cell type-specific enhancers identifies dinucleotide repeat motifs as general enhancer features

Gene expression is determined by genomic elements called enhancers, which contain short motifs bound by different transcription factors (TFs). However, how enhancer sequences and TF motifs relate to enhancer activity is unknown, and general sequence requirements for enhancers or comprehensive sets of important enhancer sequence elements have remained elusive. Here, we computationally dissect thousands of functional enhancer sequences from three different Drosophila cell lines. We find that the enhancers display distinct cis-regulatory sequence signatures, which are predictive of the enhancers’ cell type-specific or broad activities. These signatures contain transcription factor motifs and a novel class of enhancer sequence elements, dinucleotide repeat motifs (DRMs). DRMs are highly enriched in enhancers, particularly in enhancers that are broadly active across different cell types. We experimentally validate the importance of the identified TF motifs and DRMs for enhancer function and show that they can be sufficient to create an active enhancer de novo from a nonfunctional sequence. The function of DRMs as a novel class of general enhancer features that are also enriched in human regulatory regions might explain their implication in several diseases and provides important insights into gene regulation.

Proper connectivity of Drosophila motion detector neurons requires Atonal function in progenitor cells

Background

Vertebrates and invertebrates obtain visual motion information by channeling moving visual cues perceived by the retina through specific motion sensitive synaptic relays in the brain. In Drosophila, the series of synaptic relays forming the optic lobe are known as the lamina, medulla, lobula and lobula plate neuropiles. The fly’s motion detection output neurons, called the T4 and T5 cells, reside in the lobula plate. Adult optic lobe neurons are derived from larval neural progenitors in two proliferating compartments known as the outer and inner proliferation centers (OPC and IPC). Important insight has been gained into molecular mechanisms involved in the development of the lamina and medulla from the OPC, though less is known about the development of the lobula and lobula plate.

Results

Here we show that the proneural gene Atonal is expressed in a subset of IPC progenitors that give rise to the higher order motion detection neurons, T4 and T5, of the lobula plate. We also show that Atonal does not act as a proneural gene in this context. Rather, it is required specifically in IPC neural progenitors to regulate neurite outgrowth in the neuronal progeny.

Conclusions

Our findings reveal that a proneural gene is expressed in progenitors but is required for neurite development of their progeny neurons. This suggests that transcriptional programs initiated specifically in progenitors are necessary for subsequent neuronal morphogenesis.

Genome-wide analyses of Shavenbaby target genes reveals distinct features of enhancer organization

Background

Developmental programs are implemented by regulatory interactions between Transcription Factors (TFs) and their target genes, which remain poorly understood. While recent studies have focused on regulatory cascades of TFs that govern early development, little is known about how the ultimate effectors of cell differentiation are selected and controlled. We addressed this question during late Drosophila embryogenesis, when the finely tuned expression of the TF Ovo/Shavenbaby (Svb) triggers the morphological differentiation of epidermal trichomes.

Results

We defined a sizeable set of genes downstream of Svb and used in vivo assays to delineate 14 enhancers driving their specific expression in trichome cells. Coupling computational modeling to functional dissection, we investigated the regulatory logic of these enhancers. Extending the repertoire of epidermal effectors using genome-wide approaches showed that the regulatory models learned from this first sample are representative of the whole set of trichome enhancers. These enhancers harbor remarkable features with respect to their functional architectures, including a weak or non-existent clustering of Svb binding sites. The in vivo function of each site relies on its intimate context, notably the flanking nucleotides. Two additional cis-regulatory motifs, present in a broad diversity of composition and positioning among trichome enhancers, critically contribute to enhancer activity.

Conclusions

Our results show that Svb directly regulates a large set of terminal effectors of the remodeling of epidermal cells. Further, these data reveal that trichome formation is underpinned by unexpectedly diverse modes of regulation, providing fresh insights into the functional architecture of enhancers governing a terminal differentiation program.