Tiling Array Analysis Research Articles

Global regulation of gene expression by OxyR in an important human opportunistic pathogen

Most bacteria control oxidative stress through the H2O2-responsive transactivator OxyR, a member of the LTTR family (LysR Type Transcriptional Regulators), which activates the expression of defensive genes such as those encoding catalases, alkyl hydroperoxide reductases and superoxide dismutases. In the human opportunistic pathogen Pseudomonas aeruginosa, OxyR positively regulates expression of the oxidative stress response genes katA, katB, ahpB and ahpCF. To identify additional targets of OxyR in P. aeruginosa PAO1, we performed chromatin immunoprecipitation in combination with whole genome tiling array analyses (ChIP-chip). We detected 56 genes including all the previously identified defensive genes and a battery of novel direct targets of OxyR. Electrophoretic mobility shift assays (EMSAs) for selected newly identified targets indicated that ∼70% of those were bound by purified oxidized OxyR and their regulation was confirmed by quantitative real-time polymerase chain reaction. Furthermore, a thioredoxin system was identified to enzymatically reduce OxyR under oxidative stress. Functional classification analysis showed that OxyR controls a core regulon of oxidative stress defensive genes, and other genes involved in regulation of iron homeostasis (pvdS), quorum-sensing (rsaL), protein synthesis (rpsL) and oxidative phosphorylation (cyoA and snr1). Collectively, our results indicate that OxyR is involved in oxidative stress defense and regulates other aspects of cellular metabolism as well.

Genome-wide analyses of Nkx2-1 binding to transcriptional target genes uncover novel regulatory patterns conserved in lung development and tumors.

The homeodomain transcription factor Nkx2-1 is essential for normal lung development and homeostasis. In lung tumors, it is considered a lineage survival oncogene and prognostic factor depending on its expression levels. The target genes directly bound by Nkx2-1, that could be the primary effectors of its functions in the different cellular contexts where it is expressed, are mostly unknown. In embryonic day 11.5 (E11.5) mouse lung, epithelial cells expressing Nkx2-1 are predominantly expanding, and in E19.5 prenatal lungs, Nkx2-1-expressing cells are predominantly differentiating in preparation for birth. To evaluate Nkx2-1 regulated networks in these two cell contexts, we analyzed genome-wide binding of Nkx2-1 to DNA regulatory regions by chromatin immunoprecipitation followed by tiling array analysis, and intersected these data to expression data sets. We further determined expression patterns of Nkx2-1 developmental target genes in human lung tumors and correlated their expression levels to that of endogenous NKX2-1. In these studies we uncovered differential Nkx2-1 regulated networks in early and late lung development, and a direct function of Nkx2-1 in regulation of the cell cycle by controlling the expression of proliferation-related genes. New targets, validated in Nkx2-1 shRNA transduced cell lines, include E2f3, Cyclin B1, Cyclin B2, and c-Met. Expression levels of Nkx2-1 direct target genes identified in mouse development significantly correlate or anti-correlate to the levels of endogenous NKX2-1 in a dosage-dependent manner in multiple human lung tumor expression data sets, supporting alternative roles for Nkx2-1 as a transcriptional activator or repressor, and direct regulator of cell cycle progression in development and tumors.

Methylated Retroviral Integration Mutagenesis (MRIM) in Mice and Comparative Analysis in Human AML Identify Reduced PTP4A3 Expression As a Prognostic Indicator

Abstract 746Acute myeloid leukemia (AML) is a complex disease driven by multiple cytogenetic abnormalities, aberrant expression levels and mutations of different genes. Frequent occurrence of chromosomal deletions suggests that haplo-insufficiencies contribute to the pathogenesis of AML. However, because deleted regions often harbor numerous genes, it remains difficult to pin point critical haplo-insufficient genes (HIGs) involved in the pathogenesis of AML. Retroviral insertion mutagenesis in mouse models has been used to discover novel genes involved in the development of different types of cancer. Only a small minority of identified genes in these screens has been classified as tumor suppressor genes or HIGs, based on disruption of coding sequences by proviral integrations. Gene therapy studies using murine leukemia virus based vectors have shown that epigenetic changes of long terminal repeats (LTRs) of integrated proviruses often result in silencing of therapeutic genes, and that preventing methylation of the CpG islands within LTRs overcomes this problem. Based on these observations, we hypothesized that methylation of viral sequences not only results in silencing of retroviral genes themselves but may also affect host genes located proximal to proviral integrations. Methylated LTRs located in proximity of promoter regions may thus identify genes that are deregulated leading to haplo-insufficiency.To identify novel HIGs, genomic DNA of 6 murine leukemias induced by the Graffi 1.4 murine leukemia virus, was isolated and digested. Subsequently, methylated fragments were immunoprecipitated using methylated DNA immunoprecipitation (MeDIP), followed by inverse PCR to amplify fragments flanking methylated viral integrations (mVIS). Amplified fragments were labeled, fragmented and hybridized on murine promoter DNA arrays. DNA methylated viral integration sites were detected using a novel algorithm, hypergeometric analysis of tiling arrays (HAT) [1]. All integrations were confirmed by directed PCR and Sanger sequencing. Next, their methylation status was confirmed with methylation sensitive restriction analysis. We refer to this strategy as Methylated Virus Integration Mutagenesis (MRIM). In the murine leukemias, expression levels of genes identified by MRIM, relative to Tbp (TATA box binding protein) expression, were determined by quantitative PCR (qPCR). Survival analysis was performed in a cohort of 454 de novo AML patients under the age of 60 using a Cox-regression model. For this analysis, MAS5 normalised, log2 transformed PTP4A3 gene expression levels were used, determined by gene expression profiling.Using MRIM, 6 genes that were flanked by a methylated LTR were identified. Of these genes, Ptp4a3 was transcriptionally down regulated relative to normal bone marrow. Ptp4a3 encodes a tyrosine phosphatase that is implicated in Flt3 signaling. Studies in a human AML patient cohort under the age of 60 (n=454) revealed that PTP4A3 expression is negatively correlated with prognostic outcome, both for overall survival (p-value <0.0001, HR = 1.269) and event-free survival (p-value <0.0001, HR = 1.261). Additional analysis showed that this negative correlation with survival was independent of other prognostic parameters as age, white blood cell count, cytogenetic risk and mutation status of NPM1, CEBPA and internal tandem duplication of FLT3 (FLT3ITD).In conclusion, by mapping DNA methylated viral integration sites in murine leukemias induced by retroviral integration mutagenesis (MRIM), followed by expression level analysis, we identified Ptp4a3 as a candidate HIG. PTP4A3 has recently been proposed to have a role in drug-resistance in AMLs with FLT3ITD. This finding, together with the observation that high PTP4A3 expression negatively correlates with prognostic outcome, marks PTP4A3 as a potential therapeutic target in AML.1. E. Taskesen et al., HAT: hypergeometric analysis of tiling-arrays with application to promoter-GeneChip data. BMC Bioinformatics, 2010. 11: p. 275. Disclosures:No relevant conflicts of interest to declare.

Discovering transcription factor binding sites in highly repetitive regions of genomes with multi-read analysis of ChIP-Seq data.

Chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq) is rapidly replacing chromatin immunoprecipitation combined with genome-wide tiling array analysis (ChIP-chip) as the preferred approach for mapping transcription-factor binding sites and chromatin modifications. The state of the art for analyzing ChIP-seq data relies on using only reads that map uniquely to a relevant reference genome (uni-reads). This can lead to the omission of up to 30% of alignable reads. We describe a general approach for utilizing reads that map to multiple locations on the reference genome (multi-reads). Our approach is based on allocating multi-reads as fractional counts using a weighted alignment scheme. Using human STAT1 and mouse GATA1 ChIP-seq datasets, we illustrate that incorporation of multi-reads significantly increases sequencing depths, leads to detection of novel peaks that are not otherwise identifiable with uni-reads, and improves detection of peaks in mappable regions. We investigate various genome-wide characteristics of peaks detected only by utilization of multi-reads via computational experiments. Overall, peaks from multi-read analysis have similar characteristics to peaks that are identified by uni-reads except that the majority of them reside in segmental duplications. We further validate a number of GATA1 multi-read only peaks by independent quantitative real-time ChIP analysis and identify novel target genes of GATA1. These computational and experimental results establish that multi-reads can be of critical importance for studying transcription factor binding in highly repetitive regions of genomes with ChIP-seq experiments.

Histone modifications influence mediator interactions with chromatin

The Mediator complex transmits activation signals from DNA bound transcription factors to the core transcription machinery. Genome wide localization studies have demonstrated that Mediator occupancy not only correlates with high levels of transcription, but that the complex also is present at transcriptionally silenced locations. We provide evidence that Mediator localization is guided by an interaction with histone tails, and that this interaction is regulated by their post-translational modifications. A quantitative, high-density genetic interaction map revealed links between Mediator components and factors affecting chromatin structure, especially histone deacetylases. Peptide binding assays demonstrated that pure wild-type Mediator forms stable complexes with the tails of Histone H3 and H4. These binding assays also showed Mediator—histone H4 peptide interactions are specifically inhibited by acetylation of the histone H4 lysine 16, a residue critical in transcriptional silencing. Finally, these findings were validated by tiling array analysis that revealed a broad correlation between Mediator and nucleosome occupancy in vivo, but a negative correlation between Mediator and nucleosomes acetylated at histone H4 lysine 16. Our studies show that chromatin structure and the acetylation state of histones are intimately connected to Mediator localization.

Genome-Wide Analysis Reveals PADI4 Cooperates with Elk-1 to Activate c-Fos Expression in Breast Cancer Cells

Peptidylarginine deiminase IV (PADI4) catalyzes the conversion of positively charged arginine and methylarginine residues to neutrally charged citrulline, and this activity has been linked to the repression of a limited number of target genes. To broaden our knowledge of the regulatory potential of PADI4, we utilized chromatin immunoprecipitation coupled with promoter tiling array (ChIP-chip) analysis to more comprehensively investigate the range of PADI4 target genes across the genome in MCF-7 breast cancer cells. Results showed that PADI4 is enriched in gene promoter regions near transcription start sites (TSSs); and, surprisingly, this pattern of binding is primarily associated with actively transcribed genes. Computational analysis found potential binding sites for Elk-1, a member of the ETS oncogene family, to be highly enriched around PADI4 binding sites; and coimmunoprecipitation analysis then confirmed that Elk-1 physically associates with PADI4. To better understand how PADI4 may facilitate gene transactivation, we then show that PADI4 interacts with Elk-1 at the c-Fos promoter and that, following Epidermal Growth Factor (EGF) stimulation, PADI4 catalytic activity facilitates Elk-1 phosphorylation, histone H4 acetylation, and c-Fos transcriptional activation. These results define a novel role for PADI4 as a transcription factor co-activator.

ENCODE Tiling Array Analysis Identifies Differentially Expressed Annotated and Novel 5′ Capped RNAs in Hepatitis C Infected Liver

Microarray studies of chronic hepatitis C infection have provided valuable information regarding the host response to viral infection. However, recent studies of the human transcriptome indicate pervasive transcription in previously unannotated regions of the genome and that many RNA transcripts have short or lack 3′ poly(A) ends. We hypothesized that using ENCODE tiling arrays (1% of the genome) in combination with affinity purifying Pol II RNAs by their unique 5′ m7GpppN cap would identify previously undescribed annotated and unannotated genes that are differentially expressed in liver during hepatitis C virus (HCV) infection. Both 5′-capped and poly(A)+ populations of RNA were analyzed using ENCODE tiling arrays. Sixty-four annotated genes were significantly increased in HCV cirrhotic as compared to control liver; twenty-seven (42%) of these genes were identified only by analyzing 5′ capped RNA. Thirty-one annotated genes were significantly decreased; sixteen (50%) of these were identified only by analyzing 5′ capped RNA. Bioinformatic analysis showed that capped RNA produced more consistent results, provided a more extensive expression profile of intronic regions and identified upregulated Pol II transcriptionally active regions in unannotated areas of the genome in HCV cirrhotic liver. Two of these regions were verified by PCR and RACE analysis. qPCR analysis of liver biopsy specimens demonstrated that these unannotated transcripts, as well as IRF1, TRIM22 and MET, were also upregulated in hepatitis C with mild inflammation and no fibrosis. The analysis of 5′ capped RNA in combination with ENCODE tiling arrays provides additional gene expression information and identifies novel upregulated Pol II transcripts not previously described in HCV infected liver. This approach, particularly when combined with new RNA sequencing technologies, should also be useful in further defining Pol II transcripts differentially regulated in specific disease states and in studying RNAs regulated by changes in pre-mRNA splicing or 3′ polyadenylation status.

Tiling Array Analysis of UV Treated Escherichia coli Predicts Novel Differentially Expressed Small Peptides

BackgroundDespite comprehensive investigation, the Escherichia coli SOS response system is not yet fully understood. We have applied custom designed whole genome tiling arrays to measure UV invoked transcriptional changes in E. coli. This study provides a more complete insight into the transcriptome and the UV irradiation response of this microorganism.ResultsWe detected a number of novel differentially expressed transcripts in addition to the expected SOS response genes (such as sulA, recN, uvrA, lexA, umuC and umuD) in the UV treated cells. Several of the differentially expressed transcripts might play important roles in regulation of the cellular response to UV damage. We have predicted 23 novel small peptides from our set of detected non-gene transcripts. Further, three of the predicted peptides were cloned into protein expression vectors to test the biological activity. All three constructs expressed the predicted peptides, in which two of them were highly toxic to the cell. Additionally, a remarkably high overlap with previously in-silico predicted non-coding RNAs (ncRNAs) was detected. Generally we detected a far higher transcriptional activity than the annotation suggests, and these findings correspond with previous transcription mappings from E. coli and other organisms.ConclusionsHere we demonstrate that the E. coli transcriptome consists of far more transcripts than the present annotation suggests, of which many transcripts seem important to the bacterial stress response. Sequence alignment of promoter regions suggest novel regulatory consensus sequences for some of the upregulated genes. Finally, several of the novel transcripts identified in this study encode putative small peptides, which are biologically active.

Systematic Identification of Placental Epigenetic Signatures for the Noninvasive Prenatal Detection of Edwards Syndrome

BackgroundNoninvasive prenatal diagnosis of fetal aneuploidy by maternal plasma analysis is challenging owing to the low fractional and absolute concentrations of fetal DNA in maternal plasma. Previously, we demonstrated for the first time that fetal DNA in maternal plasma could be specifically targeted by epigenetic (DNA methylation) signatures in the placenta. By comparing one such methylated fetal epigenetic marker located on chromosome 21 with another fetal genetic marker located on a reference chromosome in maternal plasma, we could infer the relative dosage of fetal chromosome 21 and noninvasively detect fetal trisomy 21. Here we apply this epigenetic-genetic (EGG) chromosome dosage approach to detect Edwards syndrome (trisomy 18) in the fetus noninvasively.Principal FindingsWe have systematically identified methylated fetal epigenetic markers on chromosome 18 by methylated DNA immunoprecipitation (MeDIP) and tiling array analysis with confirmation using quantitative DNA methylation assays. Methylated DNA sequences from an intergenic region between the VAPA and APCDD1 genes (the VAPA-APCDD1 DNA) were detected in pre-delivery, but not post-delivery, maternal plasma samples. The concentrations correlated positively with those of an established fetal genetic marker, ZFY, in pre-delivery maternal plasma. The ratios of methylated VAPA-APCDD1(chr18) to ZFY(chrY) were higher in maternal plasma samples of 9 male trisomy 18 fetuses than those of 27 male euploid fetuses (Mann-Whitney test, P = 0.029). We defined the cutoff value for detecting trisomy 18 fetuses as mean+1.96 SD of the EGG ratios of the euploid cases. Eight of 9 trisomy 18 and 1 of 27 euploid cases showed EGG ratios higher than the cutoff value, giving a sensitivity of 88.9% and a specificity of 96.3%.ConclusionsOur data have shown that the methylated VAPA-APCDD1 DNA in maternal plasma is predominantly derived from the fetus. We have demonstrated that this novel fetal epigenetic marker in maternal plasma is useful for the noninvasive detection of fetal trisomy 18.

Global effects of the small RNA biogenesis machinery on the Arabidopsis thaliana transcriptome

In Arabidopsis thaliana, four different dicer-like (DCL) proteins have distinct but partially overlapping functions in the biogenesis of microRNAs (miRNAs) and siRNAs from longer, noncoding precursor RNAs. To analyze the impact of different components of the small RNA biogenesis machinery on the transcriptome, we subjected dcl and other mutants impaired in small RNA biogenesis to whole-genome tiling array analysis. We compared both protein-coding genes and noncoding transcripts, including most pri-miRNAs, in two tissues and several stress conditions. Our analysis revealed a surprising number of common targets in dcl1 and dcl2 dcl3 dcl4 triple mutants. Furthermore, our results suggest that the DCL1 is not only involved in miRNA action but also contributes to silencing of a subset of transposons, apparently through an effect on DNA methylation.

ChIPpeakAnno: a Bioconductor package to annotate ChIP-seq and ChIP-chip data

BackgroundChromatin immunoprecipitation (ChIP) followed by high-throughput sequencing (ChIP-seq) or ChIP followed by genome tiling array analysis (ChIP-chip) have become standard technologies for genome-wide identification of DNA-binding protein target sites. A number of algorithms have been developed in parallel that allow identification of binding sites from ChIP-seq or ChIP-chip datasets and subsequent visualization in the University of California Santa Cruz (UCSC) Genome Browser as custom annotation tracks. However, summarizing these tracks can be a daunting task, particularly if there are a large number of binding sites or the binding sites are distributed widely across the genome.ResultsWe have developed ChIPpeakAnno as a Bioconductor package within the statistical programming environment R to facilitate batch annotation of enriched peaks identified from ChIP-seq, ChIP-chip, cap analysis of gene expression (CAGE) or any experiments resulting in a large number of enriched genomic regions. The binding sites annotated with ChIPpeakAnno can be viewed easily as a table, a pie chart or plotted in histogram form, i.e., the distribution of distances to the nearest genes for each set of peaks. In addition, we have implemented functionalities for determining the significance of overlap between replicates or binding sites among transcription factors within a complex, and for drawing Venn diagrams to visualize the extent of the overlap between replicates. Furthermore, the package includes functionalities to retrieve sequences flanking putative binding sites for PCR amplification, cloning, or motif discovery, and to identify Gene Ontology (GO) terms associated with adjacent genes.ConclusionsChIPpeakAnno enables batch annotation of the binding sites identified from ChIP-seq, ChIP-chip, CAGE or any technology that results in a large number of enriched genomic regions within the statistical programming environment R. Allowing users to pass their own annotation data such as a different Chromatin immunoprecipitation (ChIP) preparation and a dataset from literature, or existing annotation packages, such as GenomicFeatures and BSgenome, provides flexibility. Tight integration to the biomaRt package enables up-to-date annotation retrieval from the BioMart database.

T1999 Gene Expression Changes in Hepatitis C Infected Liver as Measured by Encode Tiling Array Analysis of Both 5′ CAP and 3′ poly(A) RNA

T1999 Gene Expression Changes in Hepatitis C Infected Liver as Measured by Encode Tiling Array Analysis of Both 5′ CAP and 3′ poly(A) RNA

Oxytocin stimulates expression of a noncoding RNA tumor marker in a human neuroblastoma cell line

Aims A noncoding RNA, metastasis-associated lung adenocarcinoma transcript 1 (MALAT1), is upregulated in several malignant tumors. Its expression in neuroblastoma, however, is not known, and the regulatory mechanisms of MALAT1 gene expression have not been elucidated. The aim of this study is to clarify how MALAT1 gene expression is altered by extracellular signals in the SK–N–SH neuroblastoma cell line and to define its proximal promoter in order to study the mechanism of MALAT1 gene expression. Methods Transcript amounts were analyzed by real-time semiquantitative polymerase chain reaction (qPCR). Genes coregulated with MALAT1 were identified by DNA microarray analysis. The structure of the MALAT1 transcript was delineated using a tiling microarray, and the 5′-end was determined using the rapid amplification of cDNA ends (RACE) method. We investigated binding of the cyclic AMP-responsive element binding (CREB) transcription factor to the MALAT1 promoter by using chromatin immunoprecipitation (ChIP) followed by tiling array analysis, and the results were confirmed using ChIP-qPCR. Key findings The posterior pituitary hormone oxytocin increased the levels of MALAT1 and immediate early gene transcripts as early as 15 min after stimulation. Although the expression of immediate early genes returned to basal levels after 3 h, MALAT1 transcript levels peaked 6–24 h after stimulation. We identified a shorter transcriptional initiation site and found that CREB binds to the defined proximal promoter of the MALAT1 gene. Significance The expression of the tumor marker MALAT1 ncRNA is sensitive to cell surface receptor activation by oxytocin in a neuroblastoma cell line.

Inhibition of heat shock protein 90 attenuates adenylate cyclase sensitization after chronic morphine treatment

Cellular adaptations to chronic opioid treatment result in enhanced responsiveness of adenylate cyclase and an increase in forskolin- or agonist-stimulated cAMP production. It is, however, not known whether chaperone molecules such as heat shock proteins contribute to this adenylate cyclase sensitization. Here, we report that treatment of cells with geldanamycin, an inhibitor of heat shock protein 90 (Hsp90), led to effective attenuation of morphine-induced adenylate cyclase sensitization. In SK-N-SH human neuroblastoma cells, morphine significantly increased RNA transcript and protein levels of type I adenylate cyclase, leading to sensitization. Whole-genome tiling array analysis revealed that cAMP response element-binding protein, an important mediator for cellular adaptation to morphine, associated with the proximal promoter of Hsp90AB1 not only in SK-N-SH cells but also in rat PC12 and human embryonic kidney cells. Hsp90AB1 transcript and protein levels increased significantly during morphine treatment, and co-application of geldanamycin (0.1–10 nM) effectively suppressed the increase in forskolin-activated adenylate cyclase activation by 56%. Type I adenylate cyclase, but not Hsp90AB1, underwent significant degradation during geldanamycin treatment. These results indicate that Hsp90 is a new pharmacological target for the suppression of adenylate cyclase sensitization induced by chronic morphine treatment.

Decoding pooled RNAi screens by means of barcode tiling arrays

BackgroundRNAi screens via pooled short hairpin RNAs (shRNAs) have recently become a powerful tool for the identification of essential genes in mammalian cells. In the past years, several pooled large-scale shRNA screens have identified a variety of genes involved in cancer cell proliferation. All of those studies employed microarray analysis, utilizing either the shRNA's half hairpin sequence or an additional shRNA-associated 60 nt barcode sequence as a molecular tag. Here we describe a novel method to decode pooled RNAi screens, namely barcode tiling array analysis, and demonstrate how this approach can be used to precisely quantify the abundance of individual shRNAs from a pool.ResultsWe synthesized DNA microarrays with six overlapping 25 nt long tiling probes complementary to each unique 60 nt molecular barcode sequence associated with every shRNA expression construct. By analyzing dilution series of expression constructs we show how our approach allows quantification of shRNA abundance from a pool and how it clearly outperforms the commonly used analysis via the shRNA's half hairpin sequences. We further demonstrate how barcode tiling arrays can be used to predict anti-proliferative effects of individual shRNAs from pooled negative selection screens. Out of a pool of 305 shRNAs, we identified 28 candidate shRNAs to fully or partially impair the viability of the breast carcinoma cell line MDA-MB-231. Individual validation of a subset of eleven shRNA expression constructs with potential inhibitory, as well as non-inhibitory, effects on the cell line proliferation provides further evidence for the accuracy of the barcode tiling approach.ConclusionsIn summary, we present an improved method for the rapid, quantitative and statistically robust analysis of pooled RNAi screens. Our experimental approach, coupled with commercially available lentiviral vector shRNA libraries, has the potential to greatly facilitate the discovery of putative targets for cancer therapy as well as sensitizers of drug toxicity.

A role for microRNAs in the Drosophila circadian clock.

Little is known about the contribution of translational control to circadian rhythms. To address this issue and in particular translational control by microRNAs (miRNAs), we knocked down the miRNA biogenesis pathway in Drosophila circadian tissues. In combination with an increase in circadian-mediated transcription, this severely affected Drosophila behavioral rhythms, indicating that miRNAs function in circadian timekeeping. To identify miRNA-mRNA pairs important for this regulation, immunoprecipitation of AGO1 followed by microarray analysis identified mRNAs under miRNA-mediated control. They included three core clock mRNAs-clock (clk), vrille (vri), and clockworkorange (cwo). To identify miRNAs involved in circadian timekeeping, we exploited circadian cell-specific inhibition of the miRNA biogenesis pathway followed by tiling array analysis. This approach identified miRNAs expressed in fly head circadian tissue. Behavioral and molecular experiments show that one of these miRNAs, the developmental regulator bantam, has a role in the core circadian pacemaker. S2 cell biochemical experiments indicate that bantam regulates the translation of clk through an association with three target sites located within the clk 3' untranslated region (UTR). Moreover, clk transgenes harboring mutated bantam sites in their 3' UTRs rescue rhythms of clk mutant flies much less well than wild-type CLK transgenes.

Genome wide ChIP-chip analyses reveal important roles for CTCF in Drosophila genome organization

Insulators or chromatin boundary elements are defined by their ability to block transcriptional activation by an enhancer and to prevent the spread of active or silenced chromatin. Recent studies have increasingly suggested that insulator proteins play a role in large-scale genome organization. To better understand insulator function on the global scale, we conducted a genome-wide analysis of the binding sites for the insulator protein CTCF in Drosophila by Chromatin Immunoprecipitation (ChIP) followed by a tiling-array analysis. The analysis revealed CTCF binding to many known domain boundaries within the Abd-B gene of the BX-C including previously characterized Fab-8 and MCP insulators, and the Fab-6 region. Based on this finding, we characterized the Fab-6 insulator element. In genome-wide analysis, we found that dCTCF-binding sites are often situated between closely positioned gene promoters, consistent with the role of CTCF as an insulator protein. Importantly, CTCF tends to bind gene promoters just upstream of transcription start sites, in contrast to the predicted binding sites of the insulator protein Su(Hw). These findings suggest that CTCF plays more active roles in regulating gene activity and it functions differently from other insulator proteins in organizing the Drosophila genome.

PeakSeq enables systematic scoring of ChIP-seq experiments relative to controls

Chromatin immunoprecipitation (ChIP) followed by tag sequencing (ChIP-seq) using high-throughput next-generation instrumentation is fast, replacing chromatin immunoprecipitation followed by genome tiling array analysis (ChIP-chip) as the preferred approach for mapping of sites of transcription-factor binding and chromatin modification. Using two deeply sequenced data sets for human RNA polymerase II and STAT1, each with matching input-DNA controls, we describe a general scoring approach to address unique challenges in ChIP-seq data analysis. Our approach is based on the observation that sites of potential binding are strongly correlated with signal peaks in the control, likely revealing features of open chromatin. We develop a two-pass strategy called PeakSeq to compensate for this. A two-pass strategy compensates for signal caused by open chromatin, as revealed by inclusion of the controls. The first pass identifies putative binding sites and compensates for genomic variation in the 'mappability' of sequences. The second pass filters out sites not significantly enriched compared to the normalized control, computing precise enrichments and significances. Our scoring procedure enables us to optimize experimental design by estimating the depth of sequencing required for a desired level of coverage and demonstrating that more than two replicates provides only a marginal gain in information.

Amyloid Precursor Protein Is a Primary Androgen Target Gene That Promotes Prostate Cancer Growth

Androgen receptor (AR) is a critical transcription factor that regulates various target genes and contributes to the pathophysiology of prostate cancer hormone dependently. Here, we identify amyloid precursor protein (APP) as a primary androgen target through chromatin immunoprecipitation (ChIP) combined with genome tiling array analysis (ChIP-chip). ChIP-treated DNA were obtained from prostate cancer LNCaP cells with R1881 or vehicle treatment using AR or acetylated histone H3 antibodies. Ligand-dependent AR binding was further enriched by PCR subtraction. Using chromosome 21/22 arrays, we identified APP as one of the androgen-regulated genes with adjacent functional AR binding sites. APP expression is androgen-inducible in LNCaP cells and APP immunoreactivity was correlated with poor prognosis in patients with prostate cancer. Gain-of-function and loss-of-function studies revealed that APP promotes the tumor growth of prostate cancer. The present study reveals a novel APP-mediated pathway responsible for the androgen-dependent growth of prostate cancer. Our findings will indicate that APP could be a potential molecular target for the diagnosis and treatment of prostate cancer.

Gene Transcriptional Activity Analysis in the 6q Region in Waldenstrom Macroglobulinemia (WM) Using ChIP-on-Chip Technology

Background. Multi-level genetic characterization of Waldenstrom's Macroglobulinemia (WM) is required to improve our understanding of the underlying molecular changes that lead to the initiation and progression of this rare disease. Cytogenetic and molecular studies using gene expression analysis at the mRNA level have demonstrated minimal changes in WM cells. The most commonly identified abnormality in WM is the deletion of the 6q arm. Karyotype studies showed this abnormality in 16% of WM patients, but recent studies have shown that this deletion is present in >50% of patients using cIgMFISH. We therefore sought to investigate molecular changes that occur at the level of the “transcriptome”, indicating genes that have been actively transcribed by PolIymerase II (Poll II) enzyme using ChIP on chip technology. This technique offers the advantage of discovery of genes deregulated in WM with greater confidence than when cDNA (from RNA) is hybridized to tiling arrays.Methods. Experiments were performed using genomewide location analysis with ChIP-on-Chip TranscriptionPath technology by Genpathway (San Diego, CA). TranscriptionPath is a novel assay that identifies and quantifies genomic DNA sequences undergoing active transcription within cells. We studied CD19+ selected WM samples (N=7) and CD19+ control cells obtained from normal bone marrows (N=2). The chromatin was immunoprecipitated using an antibody specific for pol II. DNA was labeled using the Affymetrix GeneChip® WT Double-Stranded DNA Terminal Labeling Kit and hybridized to the human promoter tiling arrays. Raw data from the scans were analyzed using Tiling array analysis algorithm MAT (Johnson et al., 2006) to determine the difference in Pol II binding sites/activities using a bandwidth of 300bp and p value cut off of . RefSeq mapping were performed using Galaxy (Giardine et al., 2005). Gene ontology analysis was conducted using DAVID at http://david.adcc.ncifcrf.gov/. Results: Analysis showed that most samples with WM demonstrated a higher level of transcribed genes compared to normal control. Overall, there were 9304 regions corresponding to 4975 known genes (7405 unique Genebank accession) in the WM samples that had enriched Pol II binding activities compared to normal control. We then focused on the long arm of chromosome 6 to determine changes that occur at the “transcriptome”. We did not analyze the samples for 6q deletion to avoid bias in the analysis. At p-value of level, 193 regions were identified with higher Pol II binding activities in WM patients compared with normal donor. With at least 100 bp overlap, the regions hit 144 genes (144 unique Gene Bank accession), which contained 106 unique genes. Of these, proteasome subunit beta type 1, ubiqutin-conjugating enzyme E2, mitogen activated protein kinase, cyclin C, transcription factor B1, BCL2-associated transcription factor 1, ribosomal protein 12, wilms tumor 1 associated protein (WT1), Forkhead Box O3A (FOXO), CBP/P300-interacting transactivator, Interferon gamma receptor and chemokine receptor CCR6 were among the most overtranscripted genes observed in WM compared to control. In addition, poly(a)- specific ribonuclease (parn)-like domain containing 1 (PNLDC1) showed decreased Pol II binding activity in WM compared to control.Conclusion: TranscriptionPath ChIP on Chip/ Tiling provides the first assay for comprehensive identification of sequences undergoing active transcription across the genome in WM with high sensitivity and specificity. These included both known genes such as regulators of the proteosome/ubiqutin pathways, MAPK pathways, and Bcl-2, and previously unknown genes such as FOXO and WT1. Translocation of FOXO with the MLL gene is associated with secondary acute leukemia. Similarly, WT1 has been associated with leukemias. Interestingly, we did not identify a significant number of genes underexpressed in WM samples compared to control in the 6q. Further confirmation of the role of these genes in the regulation of tumor progression in WM is underway.