Probing the SELEX Process with Next-Generation Sequencing

Background SELEX is an iterative process in which highly diverse synthetic nucleic acid libraries are selected over many rounds to finally identify aptamers with desired properties. However, little is understood as how binders are enriched during the selection course. Next-generation sequencing offers the opportunity to open the black box and observe a large part of the population dynamics during the selection process. Methodology We have performed a semi-automated SELEX procedure on the model target streptavidin starting with a synthetic DNA oligonucleotide library and compared results obtained by the conventional analysis via cloning and Sanger sequencing with next-generation sequencing. In order to follow the population dynamics during the selection, pools from all selection rounds were barcoded and sequenced in parallel. Conclusions High affinity aptamers can be readily identified simply by copy number enrichment in the first selection rounds. Based on our results, we suggest a new selection scheme that avoids a high number of iterative selection rounds while reducing time, PCR bias, and artifacts

A Reporter Screen in a Human Haploid Cell Line Identifies CYLD as a Constitutive Inhibitor of NF-κB

The development of forward genetic screens in human haploid cells has the potential to transform our understanding of the genetic basis of cellular processes unique to man. So far, this approach has been limited mostly to the identification of genes that mediate cell death in response to a lethal agent, likely due to the ease with which this phenotype can be observed. Here, we perform the first reporter screen in the near-haploid KBM7 cell line to identify constitutive inhibitors of NF-κB. CYLD was the only currently known negative regulator of NF-κB to be identified, thus uniquely distinguishing this gene. Also identified were three genes with no previous known connection to NF-κB. Our results demonstrate that reporter screens in haploid human cells can be applied to investigate the many complex signaling pathways that converge upon transcription factors

FOXM1 binds directly to non-consensus sequences in the human genome.

BACKGROUND: The Forkhead (FKH) transcription factor FOXM1 is a key regulator of the cell cycle and is overexpressed in most types of cancer. FOXM1, similar to other FKH factors, binds to a canonical FKH motif in vitro. However, genome-wide mapping studies in different cell lines have shown a lack of enrichment of the FKH motif, suggesting an alternative mode of chromatin recruitment. We have investigated the role of direct versus indirect DNA binding in FOXM1 recruitment by performing ChIP-seq with wild-type and DNA binding deficient FOXM1. RESULTS: An in vitro fluorescence polarization assay identified point mutations in the DNA binding domain of FOXM1 that inhibit binding to a FKH consensus sequence. Cell lines expressing either wild-type or DNA binding deficient GFP-tagged FOXM1 were used for genome-wide mapping studies comparing the distribution of the DNA binding deficient protein to the wild-type. This shows that interaction of the FOXM1 DNA binding domain with target DNA is essential for recruitment. Moreover, analysis of the protein interactome of wild-type versus DNA binding deficient FOXM1 shows that the reduced recruitment is not due to inhibition of protein-protein interactions. CONCLUSIONS: A functional DNA binding domain is essential for FOXM1 chromatin recruitment. Even in FOXM1 mutants with almost complete loss of binding, the protein-protein interactions and pattern of phosphorylation are largely unaffected. These results strongly support a model whereby FOXM1 is specifically recruited to chromatin through co-factor interactions by binding directly to non-canonical DNA sequences.We would like to acknowledge the Genomics and bioinformatics core at the CRUK Research Institute for the Illumina sequencing and the Proteomics core for the LC/MS-MS protein analysis for the RIME experiments. We acknowledge the support from The University of Cambridge and Cancer Research UK. The Balasubramanian Laboratory is supported by core funding from Cancer Research UK (C14303/A17197). SB is a Wellcome Trust Principle Investigator.This is the final version of the article. It first appeared from BioMed Central via http://dx.doi.org/10.1186/s13059-015-0696-

A systematic, large-scale comparison of transcription factor binding site models

Background The modelling of gene regulation is a major challenge in biomedical research. This process is dominated by transcription factors (TFs) and mutations in their binding sites (TFBSs) may cause the misregulation of genes, eventually leading to disease. The consequences of DNA variants on TF binding are modelled in silico using binding matrices, but it remains unclear whether these are capable of accurately representing in vivo binding. In this study, we present a systematic comparison of binding models for 82 human TFs from three freely available sources: JASPAR matrices, HT-SELEX-generated models and matrices derived from protein binding microarrays (PBMs). We determined their ability to detect experimentally verified “real” in vivo TFBSs derived from ENCODE ChIP-seq data. As negative controls we chose random downstream exonic sequences, which are unlikely to harbour TFBS. All models were assessed by receiver operating characteristics (ROC) analysis. Results While the area- under-curve was low for most of the tested models with only 47 % reaching a score of 0.7 or higher, we noticed strong differences between the various position-specific scoring matrices with JASPAR and HT-SELEX models showing higher success rates than PBM-derived models. In addition, we found that while TFBS sequences showed a higher degree of conservation than randomly chosen sequences, there was a high variability between individual TFBSs. Conclusions Our results show that only few of the matrix-based models used to predict potential TFBS are able to reliably detect experimentally confirmed TFBS. We compiled our findings in a freely accessible web application called ePOSSUM (http:/mutationtaster.charite.de/ePOSSUM/) which uses a Bayes classifier to assess the impact of genetic alterations on TF binding in user-defined sequences. Additionally, ePOSSUM provides information on the reliability of the prediction using our test set of experimentally confirmed binding sites

Daily rhythms of the sleep-wake cycle

The amount and timing of sleep and sleep architecture (sleep stages) are determined by several factors, important among which are the environment, circadian rhythms and time awake. Separating the roles played by these factors requires specific protocols, including the constant routine and altered sleep-wake schedules. Results from such protocols have led to the discovery of the factors that determine the amounts and distribution of slow wave and rapid eye movement sleep as well as to the development of models to determine the amount and timing of sleep. One successful model postulates two processes. The first is process S, which is due to sleep pressure (and increases with time awake) and is attributed to a 'sleep homeostat'. Process S reverses during slow wave sleep (when it is called process S'). The second is process C, which shows a daily rhythm that is parallel to the rhythm of core temperature. Processes S and C combine approximately additively to determine the times of sleep onset and waking. The model has proved useful in describing normal sleep in adults. Current work aims to identify the detailed nature of processes S and C. The model can also be applied to circumstances when the sleep-wake cycle is different from the norm in some way. These circumstances include: those who are poor sleepers or short sleepers; the role an individual's chronotype (a measure of how the timing of the individual's preferred sleep-wake cycle compares with the average for a population); and changes in the sleep-wake cycle with age, particularly in adolescence and aging, since individuals tend to prefer to go to sleep later during adolescence and earlier in old age. In all circumstances, the evidence that sleep times and architecture are altered and the possible causes of these changes (including altered S, S' and C processes) are examined

Modeling the Basal Dynamics of P53 System

The tumor suppressor p53 has become one of most investigated genes. Once activated by stress, p53 leads to cellular responses such as cell cycle arrest and apoptosis.Most previous models have ignored the basal dynamics of p53 under nonstressed conditions. To explore the basal dynamics of p53, we constructed a stochastic delay model by incorporating two negative feedback loops. We found that protein distribution of p53 under nonstressed condition is highly skewed with a fraction of cells showing high p53 levels comparable to those observed under stressed conditions. Under nonstressed conditions, asynchronous and spontaneous p53 pulses are triggered by basal DNA double strand breaks produced during normal cell cycle progression. The first peaking times show a predominant G1 distribution while the second ones are more widely distributed. The spontaneous pulses are triggered by an excitable mechanism. Once initiated, the amplitude and duration of pulses remain unchanged. Furthermore, the spontaneous pulses are filtered by ataxia telangiectasia mutated protein mediated posttranslational modifications and do not result in substantial p21 transcription. If challenged by externally severe DNA damage, cells generate synchronous p53 pulses and induce significantly high levels of p21. The high expression of p21 can also be partially induced by lowering the deacetylation rate.Our results demonstrated that the dynamics of p53 under nonstressed conditions is initiated by an excitable mechanism and cells become fully responsive only when cells are confronted with severe damage. These findings advance our understanding of the mechanism of p53 pulses and unlock many opportunities to p53-based therapy

SalMotifDB:a tool for analyzing putative transcription factor binding sites in salmonid genomes

publishedVersio

Expression of circadian clock genes and proteins in urothelial cancer is related to cancer-associated genes

Background: The purpose of this study was to evaluate invasive and metastatic potential of urothelial cancer by investigating differential expression of various clock genes/proteins participating in the 24 h circadian rhythms and to compare these gene expressions with transcription of other cancer-associated genes. Methods: Twenty seven paired samples of tumour and benign tissue collected from patients who underwent cystectomy were analysed and compared to 15 samples of normal bladder tissue taken from patients who underwent cystoscopy for benign prostate hyperplasia (unrelated donors). Immunohistochemical analyses were made for clock and clock-related proteins. In addition, the gene-expression levels of 22 genes (clock genes, casein kinases, oncogenes, tumour suppressor genes and cytokeratins) were analysed by real-time quantitative PCR (qPCR). Results: Considerable up- or down-regulation and altered cellular distribution of different clock proteins, a reduction of casein kinase1A1 (CSNK1A1) and increase of casein kinase alpha 1 E (CSNK1E) were found. The pattern was significantly correlated with simultaneous up-regulation of stimulatory tumour markers, and a down-regulation of several suppressor genes. The pattern was mainly seen in aneuploid high-grade cancers. Considerable alterations were also found in the neighbouring bladder mucosa. Conclusions: The close correlation between altered expression of various clock genes and common tumour markers in urothelial cancer indicates that disturbed function in the cellular clock work may be an important additional mechanism contributing to cancer progression and malignant behaviour