Search

As colorectal cancer (CRC) mortality remains high due to metastatic spread in patients with advanced disease stages, we here aimed at identifying novel metastasis regulators in CRC and at understanding their mechanism of action. In CRC, tumor progression and metastasis formation are driven by cells that possess tumor-initiating cell (TIC) activity. We therefore employed a collection of primary patient-derived CRC tumor spheroid cultures (TSC) that are enriched for cells with TIC activity to perform a gene activating in vivo screen, intending to identify genes that enhance metastasis formation. To this end, TSC (n=4) were transduced with a lentiviral trapping vector that integrates into the genome and thereby drives the overexpression of genes located in the vicinity of the integration site (IS). Transduced TSC were then serially transplanted under the kidney capsule of immunodeficient mice (NSG) and monitored for the induction of metastasis formation. Analysis of lentiviral IS by linear amplification-mediated (LAM)-PCR revealed an enrichment of an IS within the gene locus of Estrogen Receptor Binding Site Associated Antigen 9 (EBAG9) in metastases. RNA sequencing confirmed the trapping vector mediated upregulation of a truncated version of the EBAG9 gene (exon 5-7), which strongly suggested a contribution of EBAG9 to metastatic spread. To validate these findings, both the truncated (TR) as well as the physiologically occurring full length (FL) EBAG9 gene were overexpressed in TSC (n=4) and serially transplanted into NSG mice. In addition to an augmentation in metastasis upon EBAG9 overexpression compared to controls (FL: 1.4-fold, TR: 2.4-fold), we also detected the occurrence of lymphogenous metastases when EBAG9 levels were increased, indicating that EBAG9 not only contributes to metastatic spread but also modulates tropism. To further characterize the underlying processes, global gene and protein expression analyses were performed in primary TSC as well as in the CRC cell line DLD1. Compared to control transduced cells, EBAG9 overexpression induced expression changes in cell adhesion and microtubule molecules, which led us to investigate cell migration and adhesion as potential contributing mechanisms to metastasis formation. We here observed an increase both in cell migration as well as in adhesion to the extracellular matrix protein collagen I upon EBAG9 overexpression in DLD1. Interestingly, DLD1 adhesion to lymphatic endothelial cells was significantly higher in EBAG9 FL (p=0.008) and EBAG9 TR (p=0.016) overexpressing cells, underscoring the role of EBAG9 in tropism modulation. In summary, we here identified EBAG9 as novel metastasis regulator in CRC which will not only contribute to a better mechanistic understanding of metastatic spread but also to the development of novel future strategies for preventing and targeting advanced disease stages. Citation Format: Karin Laaber, Taronish D. Dubash, Anna Maria Melzer, Tania Christiansen, Nati Ha, Benedikt Brors, Martin Schneider, Friederike Herbst, Hanno Glimm, Claudia R. Ball. A gene activating in vivo screen identifies EBAG9 as novel metastasis regulator in human colorectal cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2852.

Patient-derived cancer xenografts (PDX) are widely used to identify and evaluate novel therapeutic targets, and to test therapeutic approaches in preclinical mouse avatar trials. Despite their widespread use, potential caveats of PDX models remain considerably underappreciated. Here, we demonstrate that EBV-associated B-lymphoproliferations frequently develop following xenotransplantation of human colorectal and pancreatic carcinomas in highly immunodeficient NOD.Cg-Prkdcscid Il2rgtm1Wjl /SzJ (NSG) mice (18/47 and 4/37 mice, respectively), and in derived cell cultures in vitro. Strikingly, even PDX with carcinoma histology can host scarce EBV-infected B-lymphocytes that can fully overgrow carcinoma cells during serial passaging in vitro and in vivo. As serial xenografting is crucial to expand primary tumor tissue for biobanks and cohorts for preclinical mouse avatar trials, the emerging dominance of B-lymphoproliferations in serial PDX represents a serious confounding factor in these models. Consequently, repeated phenotypic assessments of serial PDX are mandatory at each expansion step to verify "bona fide" carcinoma xenografts.

Long non protein coding RNAs (lncRNAs) have been identified in many different organisms and cell types. Emerging examples emphasize the biological importance of these RNA species but their regulation and functions remain poorly understood. In the filamentous fungus Neurospora crassa, the annotation and characterization of lncRNAs is incomplete.We have performed a comprehensive transcriptome analysis of Neurospora crassa by using ChIP-seq, RNA-seq and polysome fractionation datasets. We have annotated and characterized 1478 long intergenic noncoding RNAs (lincRNAs) and 1056 natural antisense transcripts, indicating that 20% of the RNA Polymerase II transcripts of Neurospora are not coding for protein. Both classes of lncRNAs accumulate at lower levels than protein-coding mRNAs and they are considerably shorter. Our analysis showed that the vast majority of lincRNAs and antisense transcripts do not contain introns and carry less H3K4me2 modifications than similarly expressed protein coding genes. In contrast, H3K27me3 modifications inversely correlate with transcription of protein coding and lincRNA genes. We show furthermore most lincRNA sequences evolve rapidly, even between phylogenetically close species.Our transcriptome analyses revealed distinct features of Neurospora lincRNAs and antisense transcripts in comparison to mRNAs and showed that the prevalence of noncoding transcripts in this organism is higher than previously anticipated. The study provides a broad repertoire and a resource for further studies of lncRNAs.

Marker genes are specifically expressed in a tissue, organ or time of development. Here we used a computational screen to identify marker genes of the root in Arabidopsis thaliana. We mined the existing transcriptome datasets for genes having high expression in roots while being low in all other organs under a wide range of growth conditions. We show that the root-specificity of these genes is conserved in the sister species Arabidopsis lyrata, indicating that their expression pattern is under selective pressure. We delineated the cis-regulatory elements responsible for root-specific expression and validated two third of those in planta as bona fide root-specific regulatory sequences. We identified three motifs over-represented in these sequences, which mutation resulted in alteration of root-specific expression, demonstrating that these motifs are functionally relevant. In addition, the three motifs are also over-represented in the cis-regulatory regions of the A. lyrata orthologs of our root-specific genes, and this despite an overall low degree of sequence conservation of these regions. Our results provide a resource to assess root-identity in the model genus Arabidopsis and shed light on the evolutionary history of gene regulation in plants.

Patient-derived cancer xenografts (PDX) are widely used to identify and evaluate novel therapeutic targets, and to test therapeutic approaches in preclinical mouse avatar trials. Despite their widespread use, potential caveats of PDX models remain considerably underappreciated. Here, we demonstrate that EBV-associated B-lymphoproliferations frequently develop following xenotransplantation of human colorectal and pancreatic carcinomas in highly immunodeficient NOD.Cg-Prkdc

Precise gene expression is a fundamental aspect of organismal function and depends on the combinatorial interplay of transcription factors (TFs) with cis-regulatory DNA elements. While much is known about TF function in general, our understanding of their cell type-specific activities is still poor. To address how widely expressed transcriptional regulators modulate downstream gene activity with high cellular specificity, we have identified binding regions for the Hox TF Deformed (Dfd) in the Drosophila genome. Our analysis of architectural features within Hox cis-regulatory response elements (HREs) shows that HRE structure is essential for cell type-specific gene expression. We also find that Dfd and Ultrabithorax (Ubx), another Hox TF specifying different morphological traits, interact with non-overlapping regions in vivo, despite their similar DNA binding preferences. While Dfd and Ubx HREs exhibit comparable design principles, their motif compositions and motif-pair associations are distinct, explaining the highly selective interaction of these Hox proteins with the regulatory environment. Thus, our results uncover the regulatory code imprinted in Hox enhancers and elucidate the mechanisms underlying functional specificity of TFs in vivo.

Genome-wide linkage studies have been used to localize rare and highly penetrant prostate cancer (PRCA) susceptibility genes. Linkage studies performed in different ethnic backgrounds and populations have been somewhat disparate, resulting in multiple, often irreproducible signals because of genetic heterogeneity and high sporadic background of the disease. Our first genome-wide linkage study and subsequent fine-mapping study of Finnish hereditary prostate cancer (HPC) families gave evidence of linkage to one region. Here, we conducted subsequent scans with microsatellites and SNPs in a total of 69 Finnish HPC families. GENEHUNTER-PLUS was used for parametric and non-parametric analyses. Our microsatellite genome-wide linkage study provided evidence of linkage to 17q12-q23, with a heterogeneity LOD (HLOD) score of 3.14 in a total of 54 of the 69 families. Genome-wide SNP analysis of 59 of the 69 families gave a highest HLOD score of 3.40 at 2q37.3 under a dominant high penetrance model. Analyzing all 69 families by combining microsatellite and SNP maps also yielded HLOD scores of > 3.3 in two regions (2q37.3 and 17q12-q21.3). These significant linkage peaks on chromosome 2 and 17 confirm previous linkage evidence of a locus on 17q from other populations and provide a basis for continued research into genetic factors involved in PRCA. Fine-mapping analysis of these regions is ongoing and candidate genes at linked loci are currently under analysis.

Study Type – Aetiology (cohort) Level of Evidence 2b OBJECTIVE • To describe clinical and histopathological characteristics of Finnish familial prostate cancer (PCa) through a detailed analysis of cases in families. PATIENTS AND METHODS • In total, 202 Finnish families with 617 histopathologically confirmed PCa cases of confirmed genealogy were collected. • Complete clinical data, including age and prostate-specific antigen (PSA) at diagnosis, stage, grade and primary treatment, were gathered. The mean (range) number of affected men per family was 3 (2–8). • All the available diagnostic biopsy samples (n= 323) were collected and regraded by the same uropathologist. • A population-based cohort of 3011 hospital district Pirkanmaa PCa patients was used as a control group. RESULTS • The mean (range) year of diagnosis of PCa was 1993 (1962–2006) and the mean (range) age at diagnosis was 68 (43–98 years). • The median (range) primary PSA level was 12.0 (0.8–11 000) ng/mL. After regrading, the Gleason score was ≤6 in 38%, 7 in 37% and ≥8 in 25% of men. • The subset of familial PCa men diagnosed after 1995 had higher PSA levels (P= 9.9 × 10−6) and an earlier age of onset (P= 1.7 × 10−6) than men in the control group, although there were no differences in cancer-specific survival. CONCLUSIONS • We observed an earlier age of onset and higher PSA in familial PCa. • However, differences between sporadic and familial or hereditary PCa cannot be truly solved until genetic testing of high-risk genes in addition to family history is used to define PCa families. • We also emphasize that, when histological samples are collected over a longer study period, reanalysis of the samples by the same experienced uropathologist should be considered. What’s known on the subject? and What does the study add? Previous hospital- or population-based cross-sectional studies comparing the clinical and histopathological features of hereditary, familial and sporadic PCa either reported weak trends or no differences in features measured except the age of onset. In present study we observed higher PSA and earlier age of onset in the subset of 257 familial PCa menin Finnish PCa families.

Despite the independent evolution of multicellularity in plants and animals, the basic organization of their stem cell niches is remarkably similar. Here, we report the genome-wide regulatory potential of WUSCHEL, the key transcription factor for stem cell maintenance in the shoot apical meristem of the reference plant Arabidopsis thaliana. WUSCHEL acts by directly binding to at least two distinct DNA motifs in more than 100 target promoters and preferentially affects the expression of genes with roles in hormone signaling, metabolism, and development. Striking examples are the direct transcriptional repression of CLAVATA1, which is part of a negative feedback regulation of WUSCHEL, and the immediate regulation of transcriptional repressors of the TOPLESS family, which are involved in auxin signaling. Our results shed light on the complex transcriptional programs required for the maintenance of a dynamic and essential stem cell niche. Fil: Busch, Wolfang. Max Planck Institute For Developmental Biology; Alemania. University of Duke; Estados Unidos Fil: Miotk, Andrej. University of Heidelberg; Alemania Fil: Ariel, Federico Damian. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Instituto de Agrobiotecnología del Litoral. Universidad Nacional del Litoral. Instituto de Agrobiotecnología del Litoral; Argentina Fil: Zhao, Zhong. University of Heidelberg; Alemania Fil: Forner, Joachim. University of Heidelberg; Alemania Fil: Daum, Gabor. University of Heidelberg; Alemania Fil: Suzaki, Takuya. University of Heidelberg; Alemania Fil: Schuster, Christoph. University of Heidelberg; Alemania Fil: Schultheiss, Sebastian J.. Max Planck Institute For Developmental Biology; Alemania Fil: Leibfried, Andrea. Max Planck Institute For Developmental Biology; Alemania Fil: Haubeiß, Silke. Max Planck Institute For Developmental Biology; Alemania Fil: Ha, Nati. University of Heidelberg; Alemania Fil: Chan, Raquel Lia. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Santa Fe. Instituto de Agrobiotecnología del Litoral. Universidad Nacional del Litoral. Instituto de Agrobiotecnología del Litoral; Argentina Fil: Lohmann, Jan U.. Max Planck Institute For Developmental Biology; Alemania. University of Heidelberg; Alemania

Epithelial renewal in the Drosophila intestine is orchestrated by Intestinal Stem Cells (ISCs). Following damage or stress the intestinal epithelium produces ligands that activate the epidermal growth factor receptor (EGFR) in ISCs. This promotes their growth and division and, thereby, epithelial regeneration. Here we demonstrate that the HMG-box transcriptional repressor, Capicua (Cic), mediates these functions of EGFR signaling. Depleting Cic in ISCs activated them for division, whereas overexpressed Cic inhibited ISC proliferation and midgut regeneration. Epistasis tests showed that Cic acted as an essential downstream effector of EGFR/Ras signaling, and immunofluorescence showed that Cic’s nuclear localization was regulated by EGFR signaling. ISC-specific mRNA expression profiling and DNA binding mapping using DamID indicated that Cic represses cell proliferation via direct targets including string (Cdc25), Cyclin E, and the ETS domain transcription factors Ets21C and Pointed (pnt). pnt was required for ISC over-proliferation following Cic depletion, and ectopic pnt restored ISC proliferation even in the presence of overexpressed dominant-active Cic. These studies identify Cic, Pnt, and Ets21C as critical downstream effectors of EGFR signaling in Drosophila ISCs., This work was supported by the DKFZ, DFG grant SFB 873, and ERC Advanced Grant 268515 to BAE. NH was supported by SFB638, SFB1036. MF and GJ were supported by research grants from the Spanish Government (BFU2011-23611) and Fundació La Marató de TV3 (20131730); GJ is an ICREA investigator

Background: Neuroblastoma (NB), a neural crest-derived tumor of the sympathetic nervous system, is the most common extracranial solid tumor in children. We have previously shown that genomic rearrangements activate proto-oncogenic telomerase by juxtaposing active enhancer elements to the TERT gene in a large fraction of high-risk NBs. In the present study, we applied a global approach integrating whole genome sequencing (WGS), Chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-seq) data of NB cells and tumors to identify further key oncogenes activated by enhancer-hijacking in NB. Methods: WGS was applied to search for structural rearrangements in 120 NB tumors and five NB cell lines. Corresponding RNA-seq data were used to discover mono-allelic and/or outlier expression of candidate genes potentially involved in enhancer hijacking events. ChIP-seq of 34 NB tumors and 17 NB cell lines was applied to identify active enhancer elements in NB. Circular chromatin conformation capture sequencing (4C-seq) was used to confirm physical promoter-enhancer interactions in NB cell lines. Results: WGS analyses revealed that chromosomal rearrangements are common events in NB tumors and cell lines and frequently affect regions harboring proto-oncogenes and lineage specific enhancers. ChIP-seq analyses of the chromatin mark histone 3 lysine 27 acetylation (H3K27ac), surrogate for enhancer activity, confirmed that these rearrangements recurrently juxtapose active enhancer elements to oncogenes including MYCN and MYC in NB. Intriguingly, quantification of H3K27ac ChIP-seq profiles uncovered that the enhancer elements translocated to MYC were among the most active ones within the respective epigenomes. 4C-seq analyses proofed physical interactions between translocated enhancer elements and promoters of the respective oncogenes, which is in line with their elevated expression in rearranged cases. Conclusions: Our study reveals that structural rearrangements in high-risk neuroblastoma frequently juxtapose strong enhancers to key oncogenes, including MYCN and MYC, leading to physical promoter-enhancer interactions which likely drive overexpression of the oncogenes observed in rearranged cases. The common mechanism of oncogene activation by enhancer-hijacking may open a therapeutic window for epigenetic drugs including BET or CDK7 inhibitors in high-risk NBs. Citation Format: Daniel Dreidax, Moritz Gartlgruber, Sebastian Steinhauser, Larisa Savelyeva, Ron Schwessinger, Umut Toprak, Nati Ha, Dilafruz Juraeva, Martin Peifer, Matthias Fischer, Stefan Gröschel, Kai-Oliver Henrich, Young-Gyu Park, Benedikt Brors, Matthias Schlesner, Carl Herrmann, Frank Westermann. Activation of proto-oncogenes by enhancer-hijacking in high-risk neuroblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 1023. doi:10.1158/1538-7445.AM2017-1023

Light is an important environmental cue that affects physiology and development of Neurospora crassa. The light-sensing transcription factor (TF) WCC, which consists of the GATA-family TFs WC1 and WC2, is required for light-dependent transcription. SUB1, another GATA-family TF, is not a photoreceptor but has also been implicated in light-inducible gene expression. To assess regulation and organization of the network of light-inducible genes, we analyzed the roles of WCC and SUB1 in light-induced transcription and nucleosome remodeling. We show that SUB1 co-regulates a fraction of light-inducible genes together with the WCC. WCC induces nucleosome eviction at its binding sites. Chromatin remodeling is facilitated by SUB1 but SUB1 cannot activate light-inducible genes in the absence of WCC. We identified FF7, a TF with a putative O-acetyl transferase domain, as an interaction partner of SUB1 and show their cooperation in regulation of a fraction of light-inducible and a much larger number of non light-inducible genes. Our data suggest that WCC acts as a general switch for light-induced chromatin remodeling and gene expression. SUB1 and FF7 synergistically determine the extent of light-induction of target genes in common with WCC but have in addition a role in transcription regulation beyond light-induced gene expression., Author Summary In this study we have investigated the roles of the Neurospora transcription factors (TFs) WCC and SUB1 in light-activation of transcription. In principle TFs could exert identical functions for transcriptional activation and the extent of transcription will be determined by the sum of activity of the TFs. In this case however, we found that the activity of the main blue-light photoreceptor WCC is essential for the activation of light-inducible genes. SUB1 cooperates synergistically with the WCC to enhance expression of a subset of genes controlled directly by the light-activated WCC but cannot activate its light-inducible target genes in the absence of WCC. WCC evicts nucleosomes at its binding sites. This process is supported by SUB1 at a subset of common target genes. Light-dependent nucleosome loss generally correlates with but is not dependent on induction of transcription. Light-induced nucleosome eviction by the WCC/SUB1 could sensitize promoters for activation via endogenous and exogenous cues other than light, which may modulate the plasticity of the light-responsive transcriptome.

Background Circadian clocks control rhythmic expression of a large number of genes in coordination with the 24 hour day-night cycle. The mechanisms generating circadian rhythms, their amplitude and circadian phase are dependent on a transcriptional network of immense complexity. Moreover, the contribution of post-transcriptional mechanisms in generating rhythms in RNA abundance is not known. Results Here, we analyzed the clock-controlled transcriptome of Neurospora crassa together with temporal profiles of elongating RNA polymerase II. Our data indicate that transcription contributes to the rhythmic expression of the vast majority of clock-controlled genes (ccgs) in Neurospora. The ccgs accumulate in two main clusters with peak transcription and expression levels either at dawn or dusk. Dawn-phased genes are predominantly involved in catabolic and dusk-phased genes in anabolic processes, indicating a clock-controlled temporal separation of the physiology of Neurospora. Genes whose expression is strongly dependent on the core circadian activator WCC fall mainly into the dawn-phased cluster while rhythmic genes regulated by the glucose-dependent repressor CSP1 fall predominantly into the dusk-phased cluster. Surprisingly, the number of rhythmic transcripts increases about twofold in the absence of CSP1, indicating that rhythmic expression of many genes is attenuated by the activity of CSP1. Conclusions The data indicate that the vast majority of transcript rhythms in Neurospora are generated by dawn and dusk specific transcription. Our observations suggest a substantial plasticity of the circadian transcriptome with respect to the number of rhythmic genes as well as amplitude and phase of the expression rhythms and emphasize a major role of the circadian clock in the temporal organization of metabolism and physiology. Electronic supplementary material The online version of this article (doi:10.1186/s12915-015-0126-4) contains supplementary material, which is available to authorized users.

Genes are often transcribed in random bursts followed by long periods of inactivity. Here we employ the light-activatable white collar complex (WCC) of Neurospora to study the transcriptional bursting with a population approach. Activation of WCC by a light pulse triggers a synchronized wave of transcription from the frequency promoter followed by an extended period (∼1 h) during which the promoter is refractory towards restimulation. When challenged by a second light pulse, the newly activated WCC binds to refractory promoters and has the potential to recruit RNA polymerase II (Pol II). However, accumulation of Pol II and phosphorylation of its C-terminal domain repeats at serine 5 are impaired. Our results suggest that refractory promoters carry a physical memory of their recent transcription history. Genome-wide analysis of light-induced transcription suggests that refractoriness is rather widespread and a property of promoter architecture.

The white-collar complex (WCC), the core transcription factor of the circadian clock of Neurospora, activates morning-specific expression of the transcription repressor CSP1. Newly synthesized CSP1 exists in a transient complex with the corepressor RCM1/RCO1 and the ubiquitin ligase UBR1. CSP1 is rapidly hyperphosphorylated and degraded via UBR1 and its ubiquitin conjugase RAD6. Genes controlled by CSP1 are rhythmically expressed and peak in the evening (i.e., in antiphase to morning-specific genes directly controlled by WCC). Rhythmic expression of these second-tier genes depends crucially on phosphorylation and rapid turnover of CSP1, which ensures tight coupling of CSP1 abundance and function to the circadian activity of WCC. Negative feedback of CSP1 on its own transcription buffers the amplitude of CSP1-dependent oscillations against fluctuations of WCC activity. CSP1 predominantly regulates genes involved in metabolism. It controls ergosterol synthesis and fatty acid desaturases and thereby modulates the lipid composition of membranes.

SummaryProper niche architecture is critical for stem cell function, yet only few upstream regulators are known. Here, we report that the Hox transcription factor Abdominal-B (Abd-B), active in premeiotic spermatocytes of Drosophila testes, is essential for positioning the niche to the testis anterior by regulating integrin in neighboring somatic cyst cells. Abd-B also non-cell-autonomously controls critical features within the niche, including centrosome orientation and division rates of germline stem cells. By using genome-wide binding studies, we find that Abd-B mediates its effects on integrin localization by directly controlling at multiple levels the signaling activity of the Sev ligand Boss via its direct targets src42A and sec63, two genes involved in protein trafficking and recycling. Our data show that Abd-B, through local signaling between adjucent cell types, provides positional cues for integrin localization, which is critical for placement of the distant stem cell niche and stem cell activity.

In multi-cellular organisms, spatiotemporal activity of cis-regulatory DNA elements depends on their occupancy by different transcription factors (TFs). In recent years, genome-wide ChIP-on-Chip, ChIP-Seq and DamID assays have been extensively used to unravel the combinatorial interaction of TFs with cis-regulatory modules (CRMs) in the genome. Even though genome-wide binding profiles are increasingly becoming available for different TFs, single TF binding profiles are in most cases not sufficient for dissecting complex regulatory networks. Thus, potent computational tools detecting statistically significant and biologically relevant TF-motif co-occurrences in genome-wide datasets are essential for analyzing context-dependent transcriptional regulation. We have developed COPS (Co-Occurrence Pattern Search), a new bioinformatics tool based on a combination of association rules and Markov chain models, which detects co-occurring TF binding sites (BSs) on genomic regions of interest. COPS scans DNA sequences for frequent motif patterns using a Frequent-Pattern tree based data mining approach, which allows efficient performance of the software with respect to both data structure and implementation speed, in particular when mining large datasets. Since transcriptional gene regulation very often relies on the formation of regulatory protein complexes mediated by closely adjoining TF binding sites on CRMs, COPS additionally detects preferred short distance between co-occurring TF motifs. The performance of our software with respect to biological significance was evaluated using three published datasets containing genomic regions that are independently bound by several TFs involved in a defined biological process. In sum, COPS is a fast, efficient and user-friendly tool mining statistically and biologically significant TFBS co-occurrences and therefore allows the identification of TFs that combinatorially regulate gene expression.