Global Transcriptional Network Research Articles

Transcriptomics-based analysis of sex-differentiated mechanisms of hepatotoxicity in zebrafish after long-term exposure to bisphenol AF

Bisphenol AF (BPAF) is an emerging endocrine-disrupting chemical (EDC) prevalent in the environment as one of the main substitutes for bisphenol A. Sex-specific effects of EDCs have been commonly reported and closely linked to sexually dimorphic patterns of hormone metabolism and related gene expression during different exposure windows, but our understanding of these mechanisms is still limited. Here, following 28-day exposure of adult zebrafish to an environmentally relevant concentration of BPAF at 10 μg/L, the global transcriptional networks applying RNA sequencing (RNA-seq) and Ingenuity Pathway Analysis (IPA) were respectively investigated in the male and female fish liver, connecting the sex-dependent toxicity of the long-term exposure of BPAF to molecular responses. As a result, more differentially expressed genes (DEGs) were detected in males (811) than in females (195), and spermatogenesis was the most enriched Gene Ontology (GO) functional classification in males, while circadian regulation of gene expression was the most enriched GO term in females. The expression levels of selected DEGs were routinely verified using qRT-PCR, which showed consistent alterations with the transcriptional changes in RNA-seq data. The causal network analysis by IPA suggested that the adverse outcomes of BPAF in males including liver damage, apoptosis, inflammation of organ, and liver carcinoma, associated with the regulation of several key DEGs detected in RNA-seq, could be linked to the activation of upstream regulatory molecules ifnα, yap1, and ptger2; while, the inhibition of upstream regulators hif1α, ifng, and igf1, leading to the down-regulated expression of several key DEGs, might be involved in BPAF’s effects in females. Furthermore, BPAF exposure altered hepatic histological structure and inhibited antioxidant capability in both male and female livers. Overall, this study revealed different regulation networks involved in the sex-dependent effects of BPAF on the fish liver, and these detected DEGs upon BPAF exposure might be used as potential biomarkers for further assessing sex-specific hepatotoxicity following environmental EDC exposure.

Molecular basis for lethal cross-talk between two unrelated bacterial transcription factors -the regulatory protein of a restriction-modification system and the repressor of a defective prophage.

Bacterial gene expression depends on the efficient functioning of global transcriptional networks, however their interconnectivity and orchestration rely mainly on the action of individual DNA binding proteins called transcription factors (TFs). TFs interact not only with their specific target sites, but also with secondary (off-target) sites, and vary in their promiscuity. It is not clear yet what mechanisms govern the interactions with secondary sites, and how such rewiring affects the overall regulatory network, but this could clearly constrain horizontal gene transfer. Here, we show the molecular mechanism of one such off-target interaction between two unrelated TFs in Escherichia coli: the C regulatory protein of a Type II restriction-modification system, and the RacR repressor of a defective prophage. We reveal that the C protein interferes with RacR repressor expression, resulting in derepression of the toxic YdaT protein. These results also provide novel insights into regulation of the racR-ydaST operon. We mapped the C regulator interaction to a specific off-target site, and also visualized C protein dynamics, revealing intriguing differences in single molecule dynamics in different genetic contexts. Our results demonstrate an apparent example of horizontal gene transfer leading to adventitious TF cross-talk with negative effects on the recipient's viability. More broadly, this study represents an experimentally-accessible model of a regulatory constraint on horizontal gene transfer.

Expanding the genetic toolkit helps dissect a global stress response in the early-branching species Fusobacterium nucleatum

Fusobacterium nucleatum, long known as a common oral microbe, has recently garnered attention for its ability to colonize tissues and tumors elsewhere in the human body. Clinical and epidemiological research has now firmly established F. nucleatum as an oncomicrobe associated with several major cancer types. However, with the current research focus on host associations, little is known about gene regulation in F. nucleatum itself, including global stress-response pathways that typically ensure the survival of bacteria outside their primary niche. This is due to the phylogenetic distance of Fusobacteriota to most model bacteria, their limited genetic tractability, and paucity of known gene functions. Here, we characterize a global transcriptional stress-response network governed by the extracytoplasmic function sigma factor, σE. To this aim, we developed several genetic tools for this anaerobic bacterium, including four different fluorescent marker proteins, inducible gene expression, scarless gene deletion, and transcriptional and translational reporter systems. Using these tools, we identified a σE response partly reminiscent of phylogenetically distant Proteobacteria but induced by exposure to oxygen. Although F. nucleatum lacks canonical RNA chaperones, such as Hfq, we uncovered conservation of the noncoding arm of the σE response in form of the noncoding RNA FoxI. This regulatory small RNA acts as an mRNA repressor of several membrane proteins, thereby supporting the function of σE. In addition to the characterization of a global stress response in F. nucleatum, the genetic tools developed here will enable further discoveries and dissection of regulatory networks in this early-branching bacterium.

Multi-transcriptome analysis following an acute skeletal muscle growth stimulus yields tools for discerning global and MYC regulatory networks

Myc is a powerful transcription factor implicated in epigenetic reprogramming, cellular plasticity, and rapid growth as well as tumorigenesis. Cancer in skeletal muscle is extremely rare despite marked and sustained Myc induction during loading-induced hypertrophy. Here, we investigated global, actively transcribed, stable, and myonucleus-specific transcriptomes following an acute hypertrophic stimulus in mouse plantaris. With these datasets, we define global and Myc-specific dynamics at the onset of mechanical overload-induced muscle fiber growth. Data collation across analyses reveals an under-appreciated role for the muscle fiber in extracellular matrix remodeling during adaptation, along with the contribution of mRNA stability to epigenetic-related transcript levels in muscle. We also identify Runx1 and Ankrd1 (Marp1) as abundant myonucleus-enriched loading-induced genes. We observed that a strong induction of cell cycle regulators including Myc occurs with mechanical overload in myonuclei. Additionally, in vivo Myc-controlled gene expression in the plantaris was defined using a genetic muscle fiber-specific doxycycline-inducible Myc-overexpression model. We determined Myc is implicated in numerous aspects of gene expression during early-phase muscle fiber growth. Specifically, brief induction of Myc protein in muscle represses Reverbα, Reverbβ, and Myh2 while increasing Rpl3, recapitulating gene expression in myonuclei during acute overload. Experimental, comparative, and in silico analyses place Myc at the center of a stable and actively transcribed, loading-responsive, muscle fiber–localized regulatory hub. Collectively, our experiments are a roadmap for understanding global and Myc-mediated transcriptional networks that regulate rapid remodeling in postmitotic cells. We provide open webtools for exploring the five RNA-seq datasets as a resource to the field.

Potassium perchlorate effects on primordial germ cells of developing medaka larvae

Perchlorate is a chemical compound commonly used in military artillery and equipment. It has been detected in drinking water, air, soil, and breast milk. Exposure of humans to perchlorate can occur in the theater of war and areas adjacent to military training grounds. A high concentration of perchlorate has been found to affect reproduction in vertebrates, including fish. However, whether environmental concentrations of perchlorate can affect primordial germ cells (PGCs), the founders of sperm and eggs, is not clearly understood. In the present study, we examined the effects of 0, 10, 100, and 1000 μg/L potassium perchlorate exposure on the embryonic development of medaka and their PGCs. Perchlorate exposure delayed hatching time, reduced heartbeat, inhibited migration of PGCs, and increased developmental deformities in the larvae. The 10 and 20 mg/L concentrations of perchlorate were lethal to embryos, whereas vitamin C co-treatment (1 mg/L) completely blocked perchlorate-induced mortality. RNA-seq analysis of isolated PGCs showed a non-linear pattern in expression profiles of differentially altered genes. Significantly upregulated genes were found in PGCs from the 10 and 1000 μg/L groups, whereas the 100 μg/L groups showed the highest number of significantly downregulated genes. Gene ontology analysis predicted differentially expressed genes to be involved in proteolysis, metabolic processes, peptides activity, hydrolase activity, and hormone activity. Among the cellular components, extracellular, intracellular, sarcoplasmic, and 6-phosphofructokinase and membrane-bounded processes were affected. Ingenuity Pathway Analysis of PGC transcriptomes revealed thyroid hormone signaling to be affected by all concentrations of perchlorate. The present results suggested that perchlorate affected the development of medaka larvae and vitamin C was able to ameliorate perchlorate-induced embryo mortality. Additionally, perchlorate altered the global transcriptional network in PGCs in a non-linear fashion suggesting its potential effects on developing germ cells and fertility.

EPCO-08. ATRX DEFICIENCY INDUCES DYSFUNCTIONAL HETEROCHROMATIN ARCHITECTURE IN GLIOMAS AND ESTABLISHES DISEASE-DEFINING TRANSCRIPTIONAL NETWORKS

Abstract Loss of ATRX (Alpha Thalassemia/Mental Retardation Syndrome X, a member of SWI/SNF family chromatin regulator is altered in diffuse gliomas and defines molecular subtypes with aggressive behavior. Mechanistically, ATRX regulates incorporation of histone H3.3 into chromatin sites across the genome, maintains alternative lengthening of telomeres and establishes genomic distribution of polycomb responsive genes. We have recently reported Atrx deficiency induces glioma oncogenic features via widespread alterations in chromatin accessibility using mouse Neural Progenitor Cells (mNPCs- Tp53 -/-,Atrx -/-). Surprisingly, Atrx was found to be associated with transcription start site and enhancer regions, suggesting their strong association with epigenome architecture. In this background, we have recently performed ChIP-seq for histone marks that define active transcription, enhancers, repressors and gene bodies and Cohesion molecules on Atrx intact and deficient mNPCs. Our integrated analysis reports depletion of H3K9me3 loci’s with enrichment of H3K27me3 marks that coincidently enriched with Atrx binding sites and Lamina-Associated Domains (LADs). GSEA confirmed that the genes corresponding to “newly formed LADs” in mNPC-to-astrocyte differentiation are significantly enriched for genes down-regulated in Atrx deficient mNPCs and belongs to Gene Ontology categories such as cell cycle, chromosome organization and DNA damage. Alternatively, genes corresponding to decreased LAD association are enriched for up-regulated genes in Atrx deficient mNPCs and for regulation of differentiation, adhesion and cell death. Additionally, whole-genome bisulphite sequencing further demonstrated loss of methylation marks at H3K9me3 sites in Atrx deficient mNPCs, suggesting perturbations of heterochromatin regions leading to activation of canonical signals that define glioma phenotype and disease-state. To summarize, our data establishes tangible links between Atrx deficiency and dysregulated chromatin and heterochromatin architecture in gliomas and suggests the role of Atrx in establishing global chromatin features and transcriptional networks. Further, our data unravel novel therapeutic molecules/pathways for engineering potential.

The TOR-EIN2 axis mediates nuclear signalling to modulate plant growth.

The evolutionarily conserved target of rapamycin (TOR) kinase acts as a master regulator that coordinates cell proliferation and growth by integrating nutrient, energy, hormone and stress signals in all eukaryotes1,2. Research has focused mainly on TOR-regulated translation, but how TOR orchestrates the global transcriptional network remains unclear. Here we identify ethylene-insensitive protein2 (EIN2), a central integrator3-5 that shuttles between the cytoplasm and the nucleus, as a direct substrate of TOR in Arabidopsis thaliana. Glucose-activated TOR kinase directly phosphorylates EIN2 to prevent its nuclear localization. Notably, the rapid global transcriptional reprogramming that is directed by glucose-TOR signalling is largely compromised in the ein2-5 mutant, and EIN2 negatively regulates the expression of a wide range of target genes of glucose-activated TOR that are involved in DNA replication, cell wall and lipid synthesis and various secondary metabolic pathways. Chemical, cellular and genetic analyses reveal that cell elongation and proliferation processes that are controlled by the glucose-TOR-EIN2 axis are decoupled from canonical ethylene-CTR1-EIN2 signalling, and mediated by different phosphorylation sites. Our findings reveal a molecular mechanism by which a central signalling hub is shared but differentially modulated by diverse signalling pathways using distinct phosphorylation codes that can be specified by upstream protein kinases.

Tomato Spotted Wilt Virus Benefits Its Thrips Vector by Modulating Metabolic and Plant Defense Pathways in Tomato.

Several plant viruses modulate vector fitness and behavior in ways that may enhance virus transmission. Previous studies have documented indirect, plant-mediated effects of tomato spotted wilt virus (TSWV) infection on the fecundity, growth and survival of its principal thrips vector, Frankliniella occidentalis, the western flower thrips. We conducted thrips performance and preference experiments combined with plant gene expression, phytohormone and total free amino acid analyses to determine if systemically-infected tomato plants modulate primary metabolic and defense-related pathways to culminate into a more favorable environment for the vector. In a greenhouse setting, we documented a significant increase in the number of offspring produced by F. occidentalis on TSWV-infected tomato plants compared to mock-inoculated plants, and in choice test assays, females exhibited enhanced settling on TSWV-infected leaves. Microarray analysis combined with phytohormone signaling pathway analysis revealed reciprocal modulation of key phytohormone pathways under dual attack, possibly indicating a coordinated and dampening defense against the vector on infected plants. TSWV infection, alone or in combination with thrips, suppressed genes associated with photosynthesis and chloroplast function thereby significantly impacting primary metabolism of the host plant, and hierarchical cluster and network analyses revealed that many of these genes were co-regulated with phytohormone defense signaling genes. TSWV infection increased expression of genes related to protein synthesis and degradation which was reflected in the increased total free amino acid content in virus-infected plants that harbored higher thrips populations. These results suggest coordinated gene networks that regulate plant primary metabolism and defense responses rendering virus-infected plants more conducive for vector colonization, an outcome that is potentially beneficial to the vector and the virus when considered within the context of the complex transmission biology of TSWV. To our knowledge this is the first study to identify global transcriptional networks that underlie the TSWV-thrips interaction as compared to a single mechanistic approach. Findings of this study increase our fundamental knowledge of host plant-virus-vector interactions and identifies underlying mechanisms of induced host susceptibility to the insect vector.

Differential Allele-Specific Expression Uncovers Breast Cancer Genes Dysregulated by Cis Noncoding Mutations.

Identifying cancer-relevant mutations in noncoding regions is challenging due to the large numbers of such mutations, their low levels of recurrence, and difficulties in interpreting their functional impact. To uncover genes that are dysregulated due to somatic mutations in cis, we build upon the concept of differential allele-specific expression (ASE) and introduce methods to identify genes within an individual's cancer whose ASE differs from what is found in matched normal tissue. When applied to breast cancer tumor samples, our methods detect the known allele-specific effects of copy number variation and nonsense-mediated decay. Further, genes that are found to recurrently exhibit differential ASE across samples are cancer relevant. Genes with cis mutations are enriched for differential ASE, and we find 147 potentially functional noncoding mutations cis to genes that exhibit significant differential ASE. We conclude that differential ASE is a promising means for discovering gene dysregulation due to cis noncoding mutations.

SOX21 Ensures Rostral Forebrain Identity by Suppression of WNT8B during Neural Regionalization of Human Embryonic Stem Cells.

SummaryThe generation of brain region-specific progenitors from human embryonic stem cells (hESCs) is critical for their application. However, transcriptional regulation of neural regionalization in humans is poorly understood. Here, we applied a rostrocaudal patterning system from hESCs to dissect global transcriptional networks controlling early neural regionalization. We found that SOX21 is required for rostral forebrain fate specification. SOX21 knockout led to activation of Wnt signaling, resulting in caudalization of regional identity of rostral forebrain neural progenitor cells. Moreover, we identified WNT8B as a SOX21 direct target. Deletion of WNT8B or inhibition of Wnt signaling in SOX21 knockout neural progenitor cells restored rostral forebrain identity. Furthermore, SOX21 interacted with β-catenin, interfering with the binding of TCF4/β-catenin complex to the WNT8B enhancer. Collectively, these results unveil the unknown role of SOX21 and shed light on how a transcriptional factor modulates early neural regionalization through crosstalk with a key component of Wnt signaling.

Transcriptional changes after herpes simplex virus type 1 infection in human trabecular meshwork cells.

BackgroundHerpes simplex virus type 1 (HSV-1) is causative for hypertensive anterior uveitis. Trabecular meshwork (TM) cells, which are the key cells regulating intraocular pressure (IOP), is considered to be the site of inflammation. We explored the profiles of genes expressed in human TM primary cells upon HSV-1 infection.MethodsHuman TM cells were infected with HSV-1 and total RNA was isolated. The global transcriptional gene network analyses were performed in mock-infected and HSV-1 infected TM cells. Using ingenuity pathway analysis, we determined the key biological networks upon HSV-1 infection. The results of microarray analyses were validated using quantitative PCR.ResultsTM cells had a high susceptibility to HSV-1 infection. HSV-1 induced transcriptional suppression of many components related to fibrosis in TM cells. The top biological network related to the genes which were significantly altered upon HSV-1 infection was organismal injury and abnormalities involving TGF-β1 and PDGF-BB. The results of PCR analyses for selected molecules were found to be in good agreement with the microarray data. HSV-1-infected TM cells showed an 80-fold increase in the expression of PDGF-BB, which was further increased by treatment with TGF-β1. HSV-1 also induced a 4-fold increase in the expression of the monocyte chemoattractant protein (MCP)-1, the downstream molecules of PDGF-BB.ConclusionsIn human TM cells, HSV-1 induced transcriptional suppression of many components related to fibrosis and enhanced expression of both PDGF-BB and MCP-1. Our study may provide a novel mechanism for the pathogenesis of HSV-1 infection in TM cells.

Genotype-specific suppression of multiple defense pathways in apple root during infection by Pythium ultimum

The genotype-specific defense activation in the roots of perennial tree crops to soilborne necrotrophic pathogens remains largely unknown. A recent phenotyping study indicated that the apple rootstock genotypes B.9 and G.935 have contrasting resistance responses to infection by Pythium ultimum. In the current study, a comparative transcriptome analysis by Illumina Solexa HiSeq 3000 platform was carried out to identify the global transcriptional regulation networks between the susceptible B.9 and the resistant G.935 to P. ultimum infection. Thirty-six libraries were sequenced to cover three timepoints after pathogen inoculation, with three biological replicates for each sample. The transcriptomes in the roots of the susceptible genotype B.9 were reflected by overrepresented differentially expressed genes (DEGs) with downregulated patterns and systematic suppression of cellular processes at 48 h post inoculation (hpi). In contrast, DEGs with annotated functions, such as kinase receptors, MAPK signaling, JA biosynthesis enzymes, transcription factors, and transporters, were readily induced at 24 hpi and continued up-regulation at 48 hpi in G.935 roots. The earlier and stronger defense activation is likely associated with an effective inhibition of necrosis progression in G.935 roots. Lack of effector-triggered immunity or existence of a susceptibility gene could contribute to the severely disturbed transcriptome and susceptibility in B.9 roots. The identified DEGs constitute a valuable resource for hypothesis-driven studies to elucidate the resistance/tolerance mechanisms in apple roots and validating their potential association with resistance traits.

The Role of Fur in the Transcriptional and Iron Homeostatic Response of Enterococcus faecalis.

The ferric uptake regulator (Fur) plays a major role in controlling the expression of iron homeostasis genes in bacterial organisms. In this work, we fully characterized the capacity of Fur to reconfigure the global transcriptional network and influence iron homeostasis in Enterococcus faecalis. The characterization of the Fur regulon from E. faecalis indicated that this protein (Fur) regulated the expression of genes involved in iron uptake systems, conferring to the system a high level of efficiency and specificity to respond under different iron exposure conditions. An RNAseq assay coupled with a systems biology approach allowed us to identify the first global transcriptional network activated by different iron treatments (excess and limited), with and without the presence of Fur. The results showed that changes in iron availability activated a complex network of transcriptional factors in E. faecalis, among them global regulators such as LysR, ArgR, GalRS, and local regulators, LexA and CopY, which were also stimulated by copper and zinc treatments. The deletion of Fur impacted the expression of genes encoding for ABC transporters, energy production and [Fe-S] proteins, which optimized detoxification and iron uptake under iron excess and limitation, respectively. Finally, considering the close relationship between iron homeostasis and pathogenesis, our data showed that the absence of Fur increased the internal concentration of iron in the bacterium and also affected its ability to produce biofilm. These results open new alternatives in the field of infection mechanisms of E. faecalis.

A global transcriptional network connecting noncoding mutations to changes in tumor gene expression.

Although cancer genomes are replete with noncoding mutations, the effects of these mutations remain poorly characterized. Here we perform an integrative analysis of 930 tumor whole genomes and matched transcriptomes, identifying a network of 193 noncoding loci in which mutations disrupt target gene expression. These “somatic eQTLs” (expression Quantitative Trait Loci) are frequently mutated in specific cancer tissues, and the majority can be validated in an independent cohort of 3,382 tumors. Among these, we find that the effects of noncoding mutations on DAAM1, MTG2 and HYI transcription are recapitulated in multiple cancer cell lines, and that increasing DAAM1 expression leads to invasive cell migration. Collectively the noncoding loci converge on a set of core pathways, permitting a classification of tumors into pathway-based subtypes. The somatic eQTL network is disrupted in 88% of tumors, suggesting widespread impact of noncoding mutations in cancer.

Changes of DNA topology affect the global transcription landscape and allow rapid growth of a Bacillus subtilis mutant lacking carbon catabolite repression

Bacteria are able to prioritize preferred carbon sources from complex mixtures. This is achieved by the regulatory phenomenon of carbon catabolite repression. To allow the simultaneous utilization of multiple carbon sources and to prevent the time-consuming adaptation to each individual nutrient in biotechnological applications, mutants lacking carbon catabolite repression can be used. However, such mutants often exhibit pleiotropic growth defects. We have isolated and characterized mutations that overcome the growth defect of Bacillus subtilis ccpA mutants lacking the major regulator of catabolite repression, in particular their glutamate auxotrophy. Here we show, that distinct mutations affecting the essential DNA topoisomerase I (TopA) cause glutamate prototrophy of the ccpA mutant. These suppressing variants of the TopA enzyme exhibit increased activity resulting in enhanced relaxation of the DNA. Reduced DNA supercoiling results in enhanced expression of the gltAB operon encoding the biosynthetic glutamate synthase. This is achieved by a significant re-organization of the global transcription network accompanied by re-routing of metabolism, which results in inactivation of the glutamate dehydrogenase. Our results provide a link between DNA topology, the global transcriptional network, and glutamate metabolism and suggest that specific topA mutants may be well suited for biotechnological purposes.

Growth phase-dependent expression profiles of three vital H-NS family proteins encoded on the chromosome of Pseudomonas putida KT2440 and on the pCAR1 plasmid

BackgroundH-NS family proteins are nucleoid-associated proteins that form oligomers on DNA and function as global regulators. They are found in both bacterial chromosomes and plasmids, and were suggested to be candidate effectors of the interaction between them. TurA and TurB are the predominantly expressed H-NS family proteins encoded on the chromosome of Pseudomonas putida KT2440, while Pmr is encoded on the carbazole-degradative incompatibility group P-7 plasmid pCAR1. Previous transcriptome analyses suggested that they function cooperatively, but play different roles in the global transcriptional network. In addition to differences in protein interaction and DNA-binding functions, cell expression levels are important in clarifying the detailed underlying mechanisms. Here, we determined the precise protein amounts of TurA, TurB, and Pmr in KT2440 in the presence and absence of pCAR1.ResultsThe intracellular amounts of TurA and TurB in KT2440 and KT2440(pCAR1) were determined by quantitative western blot analysis using specific antibodies. The amount of TurA decreased from the log phase (~80,000 monomers per cell) to the stationary phase (~20,000 monomers per cell), while TurB was only detectable upon entry into the stationary phase (maximum 6000 monomers per cell). Protein amounts were not affected by pCAR1 carriage. KT2440(pCAR1pmrHis), where histidine-tagged Pmr is expressed under its original promotor, was used to determine the intracellular amount of Pmr, which was constant (~30,000 monomers per cell) during cell growth. Quantitative reverse transcription PCR demonstrated that the transcriptional levels of turA and turB were consistent with protein expression, though the transcriptional and translational profiles of Pmr differed.ConclusionThe amount of TurB increases as TurA decreases, and the amount of Pmr does not affect the amounts of TurA and TurB. This is consistent with our previous observation that TurA and TurB play complementary roles, whereas Pmr works relatively independently. This study provides insight into the molecular mechanisms underlying reconstitution of the transcriptional network in KT2440 by pCAR1 carriage.

Bacterial regulon modeling and prediction based on systematic cis regulatory motif analyses.

Regulons are the basic units of the response system in a bacterial cell, and each consists of a set of transcriptionally co-regulated operons. Regulon elucidation is the basis for studying the bacterial global transcriptional regulation network. In this study, we designed a novel co-regulation score between a pair of operons based on accurate operon identification and cis regulatory motif analyses, which can capture their co-regulation relationship much better than other scores. Taking full advantage of this discovery, we developed a new computational framework and built a novel graph model for regulon prediction. This model integrates the motif comparison and clustering and makes the regulon prediction problem substantially more solvable and accurate. To evaluate our prediction, a regulon coverage score was designed based on the documented regulons and their overlap with our prediction; and a modified Fisher Exact test was implemented to measure how well our predictions match the co-expressed modules derived from E. coli microarray gene-expression datasets collected under 466 conditions. The results indicate that our program consistently performed better than others in terms of the prediction accuracy. This suggests that our algorithms substantially improve the state-of-the-art, leading to a computational capability to reliably predict regulons for any bacteria.

Abstract A28: Expanding the TLX1 regulome in T-cell acute lymphoblastic leukemia towards long noncoding RNAs

Abstract Introduction: T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive cancer that results from the malignant transformation of T-cell precursors and affects children, adolescents and adults. In T-ALL, genetic lesions in several possible oncogenes and tumor suppressors have been shown to cooperatively contribute to leukemogenesis. The TLX1 (T-cell leukemia homeobox protein-1, HOX11) oncoprotein is aberrantly expressed in in 5-10% of pediatric patients and 30% of adult T-ALL patients due to chromosomal translocations. Although many downstream protein coding targets genes of TLX1 have been identified, the non-coding network downstream of TLX1 remains elusive. In this study we expand the TLX1 regulome towards long non-coding RNAs (lncRNAs). Experimental procedures: We measured the transcriptional response of all protein coding genes and lncRNAs following TLX1 knock down in the ALL-SIL cell line by polyA and total RNA-sequencing. In addition, similar mRNA-lncRNA expression profiles of 64 primary T-ALL patient samples were generated which included five TLX1+ cases. To establish the direct transcriptional TLX1 targets, we generated TLX1 and H3K27ac ChIP-sequencing data from ALL-SIL leukemic cells. Results: We confirm direct regulation of previously established protein coding gene targets and de novo TLX1 motif discovery also identified RUNX1 as an important mediator of the global TLX1 transcriptional network (Della-Gatta et al., Nature Medicine, 2012). Complementary to these data, our analysis for the first time establishes the TLX1 driven lncRNAome in thymocyte derived leukemic cells. Remarkably, the majority of TLX1 controlled lncRNAs were upregulated suggesting that they may be implicated in the TLX1 driven repression of protein coding gene expression. Notably, an important subset of these candidates is clearly associated with H3K27ac marked super-enhancer regions. Finally, pairwise mRNA-lncRNA correlation analysis allowed functional annotation of TLX1 targeted lncRNAs. To functionally interrogate candidate TLX1 regulated lncRNAs, LNA-mediated lncRNA knockdown experiments are currently performed as well as 4C-seq to explore the regulatory interactions in which these lncRNAs are involved. Conclusion: We present the first landscaping of the genome-wide binding pattern of TLX1 and provide evidence for a previously unestablished role of lncRNAs in the TLX1 regulatory network. Citation Format: Kaat Durinck, Wouter Van Loocke, Inge Van de Walle, Joni Van der Meulen, Pieter-Jan Volders, Nadine Van Roy, Yves Benoit, Bruce Poppe, Pieter Mestdagh, Jo Vandesompele, Pieter Rondou, Tom Taghon, Jean Soulier, Pieter Van Vlierberghe, Frank Speleman. Expanding the TLX1 regulome in T-cell acute lymphoblastic leukemia towards long noncoding RNAs. [abstract]. In: Proceedings of the AACR Special Conference on Chromatin and Epigenetics in Cancer; Sep 24-27, 2015; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2016;76(2 Suppl):Abstract nr A28.

Abstract B05: Transcriptional antagonism between the cooperative oncogenes TLX1 and NOTCH1 in T-cell acute lymphoblastic leukemia.

Abstract Introduction: Combined activation of specific oncogenes is a general feature of human cancer and suggests that co-occurrence of particular oncogenic factors provides a selective advantage during cellular transformation. However, the exact molecular mechanisms by which oncoproteins cooperate during malignant transformation often remains elusive. Here, we study the functional relationship between the cooperative oncogenes NOTCH1 and TLX1 in the context of T-cell acute lymphoblastic leukemia to better understand their cooperative mechanism of action during T cell transformation. Methods: In this study, we performed chromatin immunoprecipitation followed by deep sequencing (ChIP-seq) for the TLX1 homeobox oncoprotein in the T-ALL cell line ALL-SIL and analyzed the transcriptional response before and after TLX1 modulation using microarray based methods. We performed integration of TLX1 ChIPseq data with publically available transcription factor binding profiles in T-ALL and evaluated the immunophenotypic and transcriptional effects of ectopic TLX1 expression in thymus-derived CD34+ T-cell progenitors. Results: Integration of TLX1 ChIPseq data with gene expression profiles after TLX1 knockdown in the TLX1 positive T-ALL cell line ALL-SIL, confirmed the previously established role for TLX1 as transcriptional repressor in T-ALL biology. In line with previous reports (Della-Gatta et al., Nature Medicine, 2012), de novo TLX1 motif discovery identified RUNX1 and ETS1 as important mediators of the global TLX1 transcriptional network. Next, we used TLX1 ChIPseq data to define TLX1 bound super-enhancer including several loci critically involved in T-cell biology (e.g. T-cell receptor loci, RAG2, MYB). Furthermore, Gene set Enrichment Analysis (GSEA) showed that TLX1-defined super-enhancers were significantly affected by JQ1 treatment in ALL-SIL. Integration of our TLX1 ChIP-seq data with publically available ChIP-seq data for ICN1, RUNX1 and ETS1 in T-ALL cells (Wang et al., PNAS, 2013) showed a remarkable genome-wide overlap between the binding sites of these four transcription factors. Integration of these binding patterns with transcriptional read-out revealed an unprecedented transcriptional antagonism between TLX1 and NOTCH1, in which TLX1 suppresses the oncogenic NOTCH1 transcriptional program including IL7R, NOTCH3 and c-MYC. In line with this observation, ectopic TLX1 expression in CD34+ human thymic precursor T-cells broadly interfered with the normal T-cell differentiation program causing differentiation arrest, massive apoptosis and a significant downregulation of NOTCH1 target genes including reduced IL7R-alpha surface expression. Conclusion: In conclusion, our study reveals unexpected transcriptional antagonism between the cooperative oncogenes TLX1 and NOTCH1 in the biology of T-ALL. TLX1 mediated suppression of NOTCH signaling might be critically involved in the pre-leukemic phenotype (reduced thymus size and decreased cellularity) that has been observed in Lck-TLX1 transgenic mice (De Keersmaecker et al., Nature Medicine, 2010). Moreover, our results suggests that full malignant transformation of TLX1-driven leukemias might only be possible through acquisition of secondary NOTCH1 mutations that can overcome the initial TLX1 mediated suppression of NOTCH1 signaling, which is in line with the high frequency of gain-of-function NOTCH1 mutations in TLX1 positive human T-ALL. Citation Format: Kaat Durinck, Wouter Van Loocke, Joni Van der Meulen, Inge Van de Walle, Pieter Rondou, Charles E. De Bock, Bruce Poppe, Jan Cools, Jean Soulier, Tom Taghon, Frank Speleman, Pieter Van Vlierberghe. Transcriptional antagonism between the cooperative oncogenes TLX1 and NOTCH1 in T-cell acute lymphoblastic leukemia. [abstract]. In: Proceedings of the AACR Special Conference on Hematologic Malignancies: Translating Discoveries to Novel Therapies; Sep 20-23, 2014; Philadelphia, PA. Philadelphia (PA): AACR; Clin Cancer Res 2015;21(17 Suppl):Abstract nr B05.

Integrative analysis of FOXP1 function reveals a tumor-suppressive effect in prostate cancer.

The transcriptional network of the androgen receptor (AR), a key molecule of prostate cancer, is frequently modulated by interactions with other transcriptional factors such as forkhead box protein A1 (FOXA1). However, global regulatory mechanisms of AR signaling mediated by such factors have not been well investigated. Here we conducted a chromatin immunoprecipitation sequence analysis, which revealed that another FOX family, FOXP1, is specifically regulated by both AR and FOXA1. We also found that FOXP1 acts as a tumor suppressor in prostate cancer through inhibiting cell proliferation and migration. We generated an extensive global map of FOXP1 binding sites and found that FOXP1 is directly involved in AR-mediated transcription. We demonstrated that FOXP1 has a repressive effect on AR-induced transcriptional activity or histone modification in enhancer regions. Moreover, by a global analysis of androgen-mediated transcriptional networks, we observed enrichment of FOXP1 binding genes in the gene cluster negatively regulated by FOXP1. Evaluation of FOXP1 expression in clinical samples indicated that the decreased expression of FOXP1 is another prognostic factor of prostate cancer. Taken together, our results suggest a novel mechanism in which AR-induced FOXP1 functions as a direct modulator of the AR and FOXA1 centric global transcriptional network.