A sequence motif within chromatin entry sites directs MSL establishment on the Drosophila X chromosome.

The Drosophila MSL complex associates with active genes specifically on the male X chromosome to acetylate histone H4 at lysine 16 and increase expression approximately 2-fold. To date, no DNA sequence has been discovered to explain the specificity of MSL binding. We hypothesized that sequence-specific targeting occurs at "chromatin entry sites," but the majority of sites are sequence independent. Here we characterize 150 potential entry sites by ChIP-chip and ChIP-seq and discover a GA-rich MSL recognition element (MRE). The motif is only slightly enriched on the X chromosome ( approximately 2-fold), but this is doubled when considering its preferential location within or 3' to active genes (>4-fold enrichment). When inserted on an autosome, a newly identified site can direct local MSL spreading to flanking active genes. These results provide strong evidence for both sequence-dependent and -independent steps in MSL targeting of dosage compensation to the male X chromosome.

Conservation and de novo acquisition of dosage compensation on newly evolved sex chromosomes in Drosophila.

Dosage compensation has arisen in response to the evolution of distinct male (XY) and female (XX) karyotypes. In Drosophila melanogaster, the MSL complex increases male X transcription approximately twofold. X-specific targeting is thought to occur through sequence-dependent binding to chromatin entry sites (CESs), followed by spreading in cis to active genes. We tested this model by asking how newly evolving sex chromosome arms in Drosophila miranda acquired dosage compensation. We found evidence for the creation of new CESs, with the analogous sequence and spacing as in D. melanogaster, providing strong support for the spreading model in the establishment of dosage compensation.

Integrated cistromic and expression analysis of amplified NKX2-1 in lung adenocarcinoma identifies LMO3 as a functional transcriptional target.

The NKX2-1 transcription factor, a regulator of normal lung development, is the most significantly amplified gene in human lung adenocarcinoma. To study the transcriptional impact of NKX2-1 amplification, we generated an expression signature associated with NKX2-1 amplification in human lung adenocarcinoma and analyzed DNA-binding sites of NKX2-1 by genome-wide chromatin immunoprecipitation. Integration of these expression and cistromic analyses identified LMO3, itself encoding a transcription regulator, as a candidate direct transcriptional target of NKX2-1. Further cistromic and overexpression analyses indicated that NKX2-1 can cooperate with the forkhead box transcription factor FOXA1 to regulate LMO3 gene expression. RNAi analysis of NKX2-1-amplified cells compared with nonamplified cells demonstrated that LMO3 mediates cell survival downstream from NKX2-1. Our findings provide new insight into the transcriptional regulatory network of NKX2-1 and suggest that LMO3 is a transcriptional signal transducer in NKX2-1-amplified lung adenocarcinomas.

Sequence analysis of mutations and translocations across breast cancer subtypes.

Breast carcinoma is the leading cause of cancer-related mortality in women worldwide, with an estimated 1.38 million new cases and 458,000 deaths in 2008 alone. This malignancy represents a heterogeneous group of tumours with characteristic molecular features, prognosis and responses to available therapy. Recurrent somatic alterations in breast cancer have been described, including mutations and copy number alterations, notably ERBB2 amplifications, the first successful therapy target defined by a genomic aberration. Previous DNA sequencing studies of breast cancer genomes have revealed additional candidate mutations and gene rearrangements. Here we report the whole-exome sequences of DNA from 103 human breast cancers of diverse subtypes from patients in Mexico and Vietnam compared to matched-normal DNA, together with whole-genome sequences of 22 breast cancer/normal pairs. Beyond confirming recurrent somatic mutations in PIK3CA, TP53, AKT1, GATA3 and MAP3K1, we discovered recurrent mutations in the CBFB transcription factor gene and deletions of its partner RUNX1. Furthermore, we have identified a recurrent MAGI3-AKT3 fusion enriched in triple-negative breast cancer lacking oestrogen and progesterone receptors and ERBB2 expression. The MAGI3-AKT3 fusion leads to constitutive activation of AKT kinase, which is abolished by treatment with an ATP-competitive AKT small-molecule inhibitor.

p53 prevents neurodegeneration by regulating synaptic genes.

DNA damage has been implicated in neurodegenerative disorders, including Alzheimer's disease and other tauopathies, but the consequences of genotoxic stress to postmitotic neurons are poorly understood. Here we demonstrate that p53, a key mediator of the DNA damage response, plays a neuroprotective role in a Drosophila model of tauopathy. Further, through a whole-genome ChIP-chip analysis, we identify genes controlled by p53 in postmitotic neurons. We genetically validate a specific pathway, synaptic function, in p53-mediated neuroprotection. We then demonstrate that the control of synaptic genes by p53 is conserved in mammals. Collectively, our results implicate synaptic function as a central target in p53-dependent protection from neurodegeneration.

Identification of chromatin-associated regulators of MSL complex targeting in Drosophila dosage compensation.

Sex chromosome dosage compensation in Drosophila provides a model for understanding how chromatin organization can modulate coordinate gene regulation. Male Drosophila increase the transcript levels of genes on the single male X approximately two-fold to equal the gene expression in females, which have two X-chromosomes. Dosage compensation is mediated by the Male-Specific Lethal (MSL) histone acetyltransferase complex. Five core components of the MSL complex were identified by genetic screens for genes that are specifically required for male viability and are dispensable for females. However, because dosage compensation must interface with the general transcriptional machinery, it is likely that identifying additional regulators that are not strictly male-specific will be key to understanding the process at a mechanistic level. Such regulators would not have been recovered from previous male-specific lethal screening strategies. Therefore, we have performed a cell culture-based, genome-wide RNAi screen to search for factors required for MSL targeting or function. Here we focus on the discovery of proteins that function to promote MSL complex recruitment to "chromatin entry sites," which are proposed to be the initial sites of MSL targeting. We find that components of the NSL (Non-specific lethal) complex, and a previously unstudied zinc-finger protein, facilitate MSL targeting and display a striking enrichment at MSL entry sites. Identification of these factors provides new insight into how MSL complex establishes the specialized hyperactive chromatin required for dosage compensation in Drosophila.

Sequence-specific targeting of dosage compensation in Drosophila favors an active chromatin context.

The Drosophila MSL complex mediates dosage compensation by increasing transcription of the single X chromosome in males approximately two-fold. This is accomplished through recognition of the X chromosome and subsequent acetylation of histone H4K16 on X-linked genes. Initial binding to the X is thought to occur at "entry sites" that contain a consensus sequence motif ("MSL recognition element" or MRE). However, this motif is only ∼2 fold enriched on X, and only a fraction of the motifs on X are initially targeted. Here we ask whether chromatin context could distinguish between utilized and non-utilized copies of the motif, by comparing their relative enrichment for histone modifications and chromosomal proteins mapped in the modENCODE project. Through a comparative analysis of the chromatin features in male S2 cells (which contain MSL complex) and female Kc cells (which lack the complex), we find that the presence of active chromatin modifications, together with an elevated local GC content in the surrounding sequences, has strong predictive value for functional MSL entry sites, independent of MSL binding. We tested these sites for function in Kc cells by RNAi knockdown of Sxl, resulting in induction of MSL complex. We show that ectopic MSL expression in Kc cells leads to H4K16 acetylation around these sites and a relative increase in X chromosome transcription. Collectively, our results support a model in which a pre-existing active chromatin environment, coincident with H3K36me3, contributes to MSL entry site selection. The consequences of MSL targeting of the male X chromosome include increase in nucleosome lability, enrichment for H4K16 acetylation and JIL-1 kinase, and depletion of linker histone H1 on active X-linked genes. Our analysis can serve as a model for identifying chromatin and local sequence features that may contribute to selection of functional protein binding sites in the genome.

Analyses of clinicopathological, molecular, and prognostic associations of KRAS codon 61 and codon 146 mutations in colorectal cancer: cohort study and literature review.

BACKGROUND: KRAS mutations in codons 12 and 13 are established predictive biomarkers for anti-EGFR therapy in colorectal cancer. Previous studies suggest that KRAS codon 61 and 146 mutations may also predict resistance to anti-EGFR therapy in colorectal cancer. However, clinicopathological, molecular, and prognostic features of colorectal carcinoma with KRAS codon 61 or 146 mutation remain unclear. METHODS: We utilized a molecular pathological epidemiology database of 1267 colon and rectal cancers in the Nurse's Health Study and the Health Professionals Follow-up Study. We examined KRAS mutations in codons 12, 13, 61 and 146 (assessed by pyrosequencing), in relation to clinicopathological features, and tumor molecular markers, including BRAF and PIK3CA mutations, CpG island methylator phenotype (CIMP), LINE-1 methylation, and microsatellite instability (MSI). Survival analyses were performed in 1067 BRAF-wild-type cancers to avoid confounding by BRAF mutation. Cox proportional hazards models were used to compute mortality hazard ratio, adjusting for potential confounders, including disease stage, PIK3CA mutation, CIMP, LINE-1 hypomethylation, and MSI. RESULTS: KRAS codon 61 mutations were detected in 19 cases (1.5%), and codon 146 mutations in 40 cases (3.2%). Overall KRAS mutation prevalence in colorectal cancers was 40% (=505/1267). Of interest, compared to KRAS-wild-type, overall, KRAS-mutated cancers more frequently exhibited cecal location (24% vs. 12% in KRAS-wild-type; P < 0.0001), CIMP-low (49% vs. 32% in KRAS-wild-type; P < 0.0001), and PIK3CA mutations (24% vs. 11% in KRAS-wild-type; P < 0.0001). These trends were evident irrespective of mutated codon, though statistical power was limited for codon 61 mutants. Neither KRAS codon 61 nor codon 146 mutation was significantly associated with clinical outcome or prognosis in univariate or multivariate analysis [colorectal cancer-specific mortality hazard ratio (HR) = 0.81, 95% confidence interval (CI) = 0.29-2.26 for codon 61 mutation; colorectal cancer-specific mortality HR = 0.86, 95% CI = 0.42-1.78 for codon 146 mutation]. CONCLUSIONS: Tumors with KRAS mutations in codons 61 and 146 account for an appreciable proportion (approximately 5%) of colorectal cancers, and their clinicopathological and molecular features appear generally similar to KRAS codon 12 or 13 mutated cancers. To further assess clinical utility of KRAS codon 61 and 146 testing, large-scale trials are warranted.

Quantized correlation coefficient for measuring reproducibility of ChIP-chip data.

BACKGROUND: Chromatin immunoprecipitation followed by microarray hybridization (ChIP-chip) is used to study protein-DNA interactions and histone modifications on a genome-scale. To ensure data quality, these experiments are usually performed in replicates, and a correlation coefficient between replicates is used often to assess reproducibility. However, the correlation coefficient can be misleading because it is affected not only by the reproducibility of the signal but also by the amount of binding signal present in the data. RESULTS: We develop the Quantized correlation coefficient (QCC) that is much less dependent on the amount of signal. This involves discretization of data into set of quantiles (quantization), a merging procedure to group the background probes, and recalculation of the Pearson correlation coefficient. This procedure reduces the influence of the background noise on the statistic, which then properly focuses more on the reproducibility of the signal. The performance of this procedure is tested in both simulated and real ChIP-chip data. For replicates with different levels of enrichment over background and coverage, we find that QCC reflects reproducibility more accurately and is more robust than the standard Pearson or Spearman correlation coefficients. The quantization and the merging procedure can also suggest a proper quantile threshold for separating signal from background for further analysis. CONCLUSIONS: To measure reproducibility of ChIP-chip data correctly, a correlation coefficient that is robust to the amount of signal present should be used. QCC is one such measure. The QCC statistic can also be applied in a variety of other contexts for measuring reproducibility, including analysis of array CGH data for DNA copy number and gene expression data.

Small Molecule Inhibition of CPS1 Activity through an Allosteric Pocket.

Carbamoyl phosphate synthetase 1 (CPS1) catalyzes the first step in the ammonia-detoxifying urea cycle, converting ammonia to carbamoyl phosphate under physiologic conditions. In cancer, CPS1 overexpression supports pyrimidine synthesis to promote tumor growth in some cancer types, while in others CPS1 activity prevents the buildup of toxic levels of intratumoral ammonia to allow for sustained tumor growth. Targeted CPS1 inhibitors may, therefore, provide a therapeutic benefit for cancer patients with tumors overexpressing CPS1. Herein, we describe the discovery of small-molecule CPS1 inhibitors that bind to a previously unknown allosteric pocket to block ATP hydrolysis in the first step of carbamoyl phosphate synthesis. CPS1 inhibitors are active in cellular assays, blocking both urea synthesis and CPS1 support of the pyrimidine biosynthetic pathway, while having no activity against CPS2. These newly discovered CPS1 inhibitors are a first step toward providing researchers with valuable tools for probing CPS1 cancer biology.

Somatic Mutational Landscape of Splicing Factor Genes and Their Functional Consequences across 33 Cancer Types.

Hotspot mutations in splicing factor genes have been recently reported at high frequency in hematological malignancies, suggesting the importance of RNA splicing in cancer. We analyzed whole-exome sequencing data across 33 tumor types in The Cancer Genome Atlas (TCGA), and we identified 119 splicing factor genes with significant non-silent mutation patterns, including mutation over-representation, recurrent loss of function (tumor suppressor-like), or hotspot mutation profile (oncogene-like). Furthermore, RNA sequencing analysis revealed altered splicing events associated with selected splicing factor mutations. In addition, we were able to identify common gene pathway profiles associated with the presence of these mutations. Our analysis suggests that somatic alteration of genes involved in the RNA-splicing process is common in cancer and may represent an underappreciated hallmark of tumorigenesis.

H3B-8800, an orally available small-molecule splicing modulator, induces lethality in spliceosome-mutant cancers.

Genomic analyses of cancer have identified recurrent point mutations in the RNA splicing factor-encoding genes SF3B1, U2AF1, and SRSF2 that confer an alteration of function. Cancer cells bearing these mutations are preferentially dependent on wild-type (WT) spliceosome function, but clinically relevant means to therapeutically target the spliceosome do not currently exist. Here we describe an orally available modulator of the SF3b complex, H3B-8800, which potently and preferentially kills spliceosome-mutant epithelial and hematologic tumor cells. These killing effects of H3B-8800 are due to its direct interaction with the SF3b complex, as evidenced by loss of H3B-8800 activity in drug-resistant cells bearing mutations in genes encoding SF3b components. Although H3B-8800 modulates WT and mutant spliceosome activity, the preferential killing of spliceosome-mutant cells is due to retention of short, GC-rich introns, which are enriched for genes encoding spliceosome components. These data demonstrate the therapeutic potential of splicing modulation in spliceosome-mutant cancers.

Normalization and experimental design for ChIP-chip data.

BACKGROUND: Chromatin immunoprecipitation on tiling arrays (ChIP-chip) has been widely used to investigate the DNA binding sites for a variety of proteins on a genome-wide scale. However, several issues in the processing and analysis of ChIP-chip data have not been resolved fully, including the effect of background (mock control) subtraction and normalization within and across arrays. RESULTS: The binding profiles of Drosophila male-specific lethal (MSL) complex on a tiling array provide a unique opportunity for investigating these topics, as it is known to bind on the X chromosome but not on the autosomes. These large bound and control regions on the same array allow clear evaluation of analytical methods.We introduce a novel normalization scheme specifically designed for ChIP-chip data from dual-channel arrays and demonstrate that this step is critical for correcting systematic dye-bias that may exist in the data. Subtraction of the mock (non-specific antibody or no antibody) control data is generally needed to eliminate the bias, but appropriate normalization obviates the need for mock experiments and increases the correlation among replicates. The idea underlying the normalization can be used subsequently to estimate the background noise level in each array for normalization across arrays. We demonstrate the effectiveness of the methods with the MSL complex binding data and other publicly available data. CONCLUSION: Proper normalization is essential for ChIP-chip experiments. The proposed normalization technique can correct systematic errors and compensate for the lack of mock control data, thus reducing the experimental cost and producing more accurate results.

Somatic copy number alterations in gastric adenocarcinomas among Asian and Western patients.

Gastric cancer, a leading worldwide cause of cancer mortality, shows high geographic and ethnic variation in incidence rates, which are highest in East Asia. The anatomic locations and clinical behavior also differ by geography, leading to the controversial idea that Eastern and Western forms of the disease are distinct. In view of these differences, we investigated whether gastric cancers from Eastern and Western patients show distinct genomic profiles. We used high-density profiling of somatic copy-number aberrations to analyze the largest collection to date of gastric adenocarcinomas and utilized genotyping data to rigorously annotate ethnic status. The size of this collection allowed us to accurately identify regions of significant copy-number alteration and separately to evaluate tumors arising in Eastern and Western patients. Among molecular subtypes classified by The Cancer Genome Atlas, the frequency of gastric cancers showing chromosomal instability was modestly higher in Western patients. After accounting for this difference, however, gastric cancers arising in Easterners and Westerners have highly similar somatic copy-number patterns. Only one genomic event, focal deletion of the phosphatase gene PTPRD, was significantly enriched in Western cases, though also detected in Eastern cases. Thus, despite the different risk factors and clinical features, gastric cancer appears to be a fundamentally similar disease in both populations and the divergent clinical outcomes cannot be ascribed to different underlying structural somatic genetic aberrations.

Splicing modulators act at the branch point adenosine binding pocket defined by the PHF5A-SF3b complex.

Pladienolide, herboxidiene and spliceostatin have been identified as splicing modulators that target SF3B1 in the SF3b subcomplex. Here we report that PHF5A, another component of this subcomplex, is also targeted by these compounds. Mutations in PHF5A-Y36, SF3B1-K1071, SF3B1-R1074 and SF3B1-V1078 confer resistance to these modulators, suggesting a common interaction site. RNA-seq analysis reveals that PHF5A-Y36C has minimal effect on basal splicing but inhibits the global action of splicing modulators. Moreover, PHF5A-Y36C alters splicing modulator-induced intron-retention/exon-skipping profile, which correlates with the differential GC content between adjacent introns and exons. We determine the crystal structure of human PHF5A demonstrating that Y36 is located on a highly conserved surface. Analysis of the cryo-EM spliceosome Bact complex shows that the resistance mutations cluster in a pocket surrounding the branch point adenosine, suggesting a competitive mode of action. Collectively, we propose that PHF5A-SF3B1 forms a central node for binding to these splicing modulators.

A statistically based density map method for identification and quantification of regional differences in microcolumnarity in the monkey brain.

We present a statistical density map method derived from condensed matter physics to quantify microcolumns, the fundamental computational unit of the cerebral cortex. This method provides measures for microcolumnar strength, width, spacing, length, and periodicity. We applied this method to Nissl-stained 30 microm thick frozen sections from areas 46, TE, and TL of rhesus monkey brains, areas that differ visually in microcolumnarity and are associated with different cognitive functions. Our results indicate that microcolumns in these areas are similar in width, spacing, and periodicity, but are stronger (possess a higher neuronal density) in area TE, as compared to areas TL and 46. We modeled the effect of section orientation on microcolumnar spacing and demonstrated that this method provides an adequate estimate of spacing. We also modeled disruption of microcolumnarity by performing simulations that randomly displace neurons and demonstrated that displacements of only one neuronal diameter effectively eliminate microcolumnar organization. These results indicate that our density map method is sensitive enough to detect and quantify subtle differences in microcolumnar organization that may occur in the context of development, aging, and neuropathology, as well as between areas and species.

SOX15 governs transcription in human stratified epithelia and a subset of esophageal adenocarcinomas.

BACKGROUND & AIMS: Intestinal metaplasia (Barrett's esophagus, BE) is the principal risk factor for esophageal adenocarcinoma (EAC). Study of the basis for BE has centered on intestinal factors, but loss of esophageal identity likely also reflects absence of key squamous-cell factors. As few determinants of stratified epithelial cell-specific gene expression are characterized, it is important to identify the necessary transcription factors. METHODS: We tested regional expression of mRNAs for all putative DNA-binding proteins in the mouse digestive tract and verified esophagus-specific factors in human tissues and cell lines. Integration of diverse data defined a human squamous esophagus-specific transcriptome. We used chromatin immunoprecipitation (ChIP-seq) to locate transcription factor binding sites, computational approaches to profile transcripts in cancer datasets, and immunohistochemistry to reveal protein expression. RESULTS: The transcription factor SOX15 is restricted to esophageal and other murine and human stratified epithelia. SOX15 mRNA levels are attenuated in BE and its depletion in human esophageal cells reduced esophageal transcripts significantly and specifically. SOX15 binding is highly enriched near esophagus-expressed genes, indicating direct transcriptional control. SOX15 and hundreds of genes co-expressed in squamous cells are reactivated in up to 30% of EAC specimens. Genes normally confined to the esophagus or intestine appear in different cells within the same malignant glands. CONCLUSIONS: These data identify a novel transcriptional regulator of stratified epithelial cells and a subtype of EAC with bi-lineage gene expression. Broad activation of squamous-cell genes may shed light on whether EACs arise in the native stratified epithelium or in ectopic columnar cells.

Dynamic epigenetic regulation by menin during pancreatic islet tumor formation.

UNLABELLED: The tumor suppressor gene MEN1 is frequently mutated in sporadic pancreatic neuroendocrine tumors (PanNET) and is responsible for the familial multiple endocrine neoplasia type 1 (MEN-1) cancer syndrome. Menin, the protein product of MEN1, associates with the histone methyltransferases (HMT) MLL1 (KMT2A) and MLL4 (KMT2B) to form menin-HMT complexes in both human and mouse model systems. To elucidate the role of methylation of histone H3 at lysine 4 (H3K4) mediated by menin-HMT complexes during PanNET formation, genome-wide histone H3 lysine 4 trimethylation (H3K4me3) signals were mapped in pancreatic islets using unbiased chromatin immunoprecipitation coupled with next-generation sequencing (ChIP-seq). Integrative analysis of gene expression profiles and histone H3K4me3 levels identified a number of transcripts and target genes dependent on menin. In the absence of Men1, histone H3K27me3 levels are enriched, with a concomitant decrease in H3K4me3 within the promoters of these target genes. In particular, expression of the insulin-like growth factor 2 mRNA binding protein 2 (IGF2BP2) gene is subject to dynamic epigenetic regulation by Men1-dependent histone modification in a time-dependent manner. Decreased expression of IGF2BP2 in Men1-deficient hyperplastic pancreatic islets is partially reversed by ablation of RBP2 (KDM5A), a histone H3K4-specific demethylase of the jumonji, AT-rich interactive domain 1 (JARID1) family. Taken together, these data demonstrate that loss of Men1 in pancreatic islet cells alters the epigenetic landscape of its target genes. IMPLICATIONS: Epigenetic profiling and gene expression analysis in Men1-deficient pancreatic islet cells reveals vital insight into the molecular events that occur during the progression of pancreatic islet tumorigenesis.

Paired exome analysis of Barrett's esophagus and adenocarcinoma.

Barrett's esophagus is thought to progress to esophageal adenocarcinoma (EAC) through a stepwise progression with loss of CDKN2A followed by TP53 inactivation and aneuploidy. Here we present whole-exome sequencing from 25 pairs of EAC and Barrett's esophagus and from 5 patients whose Barrett's esophagus and tumor were extensively sampled. Our analysis showed that oncogene amplification typically occurred as a late event and that TP53 mutations often occurred early in Barrett's esophagus progression, including in non-dysplastic epithelium. Reanalysis of additional EAC exome data showed that the majority (62.5%) of EACs emerged following genome doubling and that tumors with genomic doubling had different patterns of genomic alterations, with more frequent oncogenic amplification and less frequent inactivation of tumor suppressors, including CDKN2A. These data suggest that many EACs emerge not through the gradual accumulation of tumor-suppressor alterations but rather through a more direct path whereby a TP53-mutant cell undergoes genome doubling, followed by the acquisition of oncogenic amplifications.

Cancer-Associated SF3B1 Hotspot Mutations Induce Cryptic 3' Splice Site Selection through Use of a Different Branch Point.

Recurrent mutations in the spliceosome are observed in several human cancers, but their functional and therapeutic significance remains elusive. SF3B1, the most frequently mutated component of the spliceosome in cancer, is involved in the recognition of the branch point sequence (BPS) during selection of the 3' splice site (ss) in RNA splicing. Here, we report that common and tumor-specific splicing aberrations are induced by SF3B1 mutations and establish aberrant 3' ss selection as the most frequent splicing defect. Strikingly, mutant SF3B1 utilizes a BPS that differs from that used by wild-type SF3B1 and requires the canonical 3' ss to enable aberrant splicing during the second step. Approximately 50% of the aberrantly spliced mRNAs are subjected to nonsense-mediated decay resulting in downregulation of gene and protein expression. These findings ascribe functional significance to the consequences of SF3B1 mutations in cancer.