In vivo CRISPR screens identify the E3 ligase Cop1 as a modulator of macrophage infiltration and cancer immunotherapy target.

Despite remarkable clinical efficacy of immune checkpoint blockade (ICB) in cancer treatment, ICB benefits for triple-negative breast cancer (TNBC) remain limited. Through pooled in vivo CRISPR knockout (KO) screens in syngeneic TNBC mouse models, we found that deletion of the E3 ubiquitin ligase Cop1 in cancer cells decreases secretion of macrophage-associated chemokines, reduces tumor macrophage infiltration, enhances anti-tumor immunity, and strengthens ICB response. Transcriptomics, epigenomics, and proteomics analyses revealed that Cop1 functions through proteasomal degradation of the C/ebpδ protein. The Cop1 substrate Trib2 functions as a scaffold linking Cop1 and C/ebpδ, which leads to polyubiquitination of C/ebpδ. In addition, deletion of the E3 ubiquitin ligase Cop1 in cancer cells stabilizes C/ebpδ to suppress expression of macrophage chemoattractant genes. Our integrated approach implicates Cop1 as a target for improving cancer immunotherapy efficacy in TNBC by regulating chemokine secretion and macrophage infiltration in the tumor microenvironment.

A Somatically Acquired Enhancer of the Androgen Receptor Is a Noncoding Driver in Advanced Prostate Cancer.

Increased androgen receptor (AR) activity drives therapeutic resistance in advanced prostate cancer. The most common resistance mechanism is amplification of this locus presumably targeting the AR gene. Here, we identify and characterize a somatically acquired AR enhancer located 650 kb centromeric to the AR. Systematic perturbation of this enhancer using genome editing decreased proliferation by suppressing AR levels. Insertion of an additional copy of this region sufficed to increase proliferation under low androgen conditions and to decrease sensitivity to enzalutamide. Epigenetic data generated in localized prostate tumors and benign specimens support the notion that this region is a developmental enhancer. Collectively, these observations underscore the importance of epigenomic profiling in primary specimens and the value of deploying genome editing to functionally characterize noncoding elements. More broadly, this work identifies a therapeutic vulnerability for targeting the AR and emphasizes the importance of regulatory elements as highly recurrent oncogenic drivers.

Targetable leukaemia dependency on noncanonical PI3Kγ signalling.

Phosphoinositide-3-kinase-γ (PI3Kγ) is implicated as a target to repolarize tumour-associated macrophages and promote antitumour immune responses in solid cancers1-4. However, cancer cell-intrinsic roles of PI3Kγ are unclear. Here, by integrating unbiased genome-wide CRISPR interference screening with functional analyses across acute leukaemias, we define a selective dependency on the PI3Kγ complex in a high-risk subset that includes myeloid, lymphoid and dendritic lineages. This dependency is characterized by innate inflammatory signalling and activation of phosphoinositide 3-kinase regulatory subunit 5 (PIK3R5), which encodes a regulatory subunit of PI3Kγ5 and stabilizes the active enzymatic complex. We identify p21 (RAC1)-activated kinase 1 (PAK1) as a noncanonical substrate of PI3Kγ that mediates this cell-intrinsic dependency and find that dephosphorylation of PAK1 by PI3Kγ inhibition impairs mitochondrial oxidative phosphorylation. Treatment with the selective PI3Kγ inhibitor eganelisib is effective in leukaemias with activated PIK3R5. In addition, the combination of eganelisib and cytarabine prolongs survival over either agent alone, even in patient-derived leukaemia xenografts with low baseline PIK3R5 expression, as residual leukaemia cells after cytarabine treatment have elevated G protein-coupled purinergic receptor activity and PAK1 phosphorylation. Together, our study reveals a targetable dependency on PI3Kγ-PAK1 signalling that is amenable to near-term evaluation in patients with acute leukaemia.

HIF1A employs CDK8-mediator to stimulate RNAPII elongation in response to hypoxia.

The transcription factor HIF1A is a key mediator of the cellular response to hypoxia. Despite the importance of HIF1A in homeostasis and various pathologies, little is known about how it regulates RNA polymerase II (RNAPII). We report here that HIF1A employs a specific variant of the Mediator complex to stimulate RNAPII elongation. The Mediator-associated kinase CDK8, but not the paralog CDK19, is required for induction of many HIF1A target genes. HIF1A induces binding of CDK8-Mediator and the super elongation complex (SEC), containing AFF4 and CDK9, to alleviate RNAPII pausing. CDK8 is dispensable for HIF1A chromatin binding and histone acetylation, but it is essential for binding of SEC and RNAPII elongation. Global analysis of active RNAPII reveals that hypoxia-inducible genes are paused and active prior to their induction. Our results provide a mechanistic link between HIF1A and CDK8, two potent oncogenes, in the cellular response to hypoxia.

Synthetic Lethal and Resistance Interactions with BET Bromodomain Inhibitors in Triple-Negative Breast Cancer.

BET bromodomain inhibitors (BBDIs) are candidate therapeutic agents for triple-negative breast cancer (TNBC) and other cancer types, but inherent and acquired resistance to BBDIs limits their potential clinical use. Using CRISPR and small-molecule inhibitor screens combined with comprehensive molecular profiling of BBDI response and resistance, we identified synthetic lethal interactions with BBDIs and genes that, when deleted, confer resistance. We observed synergy with regulators of cell cycle progression, YAP, AXL, and SRC signaling, and chemotherapeutic agents. We also uncovered functional similarities and differences among BRD2, BRD4, and BRD7. Although deletion of BRD2 enhances sensitivity to BBDIs, BRD7 loss leads to gain of TEAD-YAP chromatin binding and luminal features associated with BBDI resistance. Single-cell RNA-seq, ATAC-seq, and cellular barcoding analysis of BBDI responses in sensitive and resistant cell lines highlight significant heterogeneity among samples and demonstrate that BBDI resistance can be pre-existing or acquired.

Cancer-stromal cell interactions in breast cancer brain metastases induce glycocalyx-mediated resistance to HER2-targeting therapies.

Brain metastatic breast cancer is particularly lethal largely due to therapeutic resistance. Almost half of the patients with metastatic HER2-positive breast cancer develop brain metastases, representing a major clinical challenge. We previously described that cancer-associated fibroblasts are an important source of resistance in primary tumors. Here, we report that breast cancer brain metastasis stromal cell interactions in 3D cocultures induce therapeutic resistance to HER2-targeting agents, particularly to the small molecule inhibitor of HER2/EGFR neratinib. We investigated the underlying mechanisms using a synthetic Notch reporter system enabling the sorting of cancer cells that directly interact with stromal cells. We identified mucins and bulky glycoprotein synthesis as top-up-regulated genes and pathways by comparing the gene expression and chromatin profiles of stroma-contact and no-contact cancer cells before and after neratinib treatment. Glycoprotein gene signatures were also enriched in human brain metastases compared to primary tumors. We confirmed increased glycocalyx surrounding cocultures by immunofluorescence and showed that mucinase treatment increased sensitivity to neratinib by enabling a more efficient inhibition of EGFR/HER2 signaling in cancer cells. Overexpression of truncated MUC1 lacking the intracellular domain as a model of increased glycocalyx-induced resistance to neratinib both in cell culture and in experimental brain metastases in immunodeficient mice. Our results highlight the importance of glycoproteins as a resistance mechanism to HER2-targeting therapies in breast cancer brain metastases.

MIRA: joint regulatory modeling of multimodal expression and chromatin accessibility in single cells.

Rigorously comparing gene expression and chromatin accessibility in the same single cells could illuminate the logic of how coupling or decoupling of these mechanisms regulates fate commitment. Here we present MIRA, probabilistic multimodal models for integrated regulatory analysis, a comprehensive methodology that systematically contrasts transcription and accessibility to infer the regulatory circuitry driving cells along cell state trajectories. MIRA leverages topic modeling of cell states and regulatory potential modeling of individual gene loci. MIRA thereby represents cell states in an efficient and interpretable latent space, infers high-fidelity cell state trees, determines key regulators of fate decisions at branch points and exposes the variable influence of local accessibility on transcription at distinct loci. Applied to epidermal differentiation and embryonic brain development from two different multimodal platforms, MIRA revealed that early developmental genes were tightly regulated by local chromatin landscape whereas terminal fate genes were titrated without requiring extensive chromatin remodeling.

Inhibition of EZH2 transactivation function sensitizes solid tumors to genotoxic stress.

Drugs that block the activity of the methyltransferase EZH2 are in clinical development for the treatment of non-Hodgkin lymphomas harboring EZH2 gain-of-function mutations that enhance its polycomb repressive function. We have previously reported that EZH2 can act as a transcriptional activator in castration-resistant prostate cancer (CRPC). Now we show that EZH2 inhibitors can also block the transactivation activity of EZH2 and inhibit the growth of CRPC cells. Gene expression and epigenomics profiling of cells treated with EZH2 inhibitors demonstrated that in addition to derepressing gene expression, these compounds also robustly down-regulate a set of DNA damage repair (DDR) genes, especially those involved in the base excision repair (BER) pathway. Methylation of the pioneer factor FOXA1 by EZH2 contributes to the activation of these genes, and interaction with the transcriptional coactivator P300 via the transactivation domain on EZH2 directly turns on the transcription. In addition, CRISPR-Cas9-mediated knockout screens in the presence of EZH2 inhibitors identified these BER genes as the determinants that underlie the growth-inhibitory effect of EZH2 inhibitors. Interrogation of public data from diverse types of solid tumors expressing wild-type EZH2 demonstrated that expression of DDR genes is significantly correlated with EZH2 dependency and cellular sensitivity to EZH2 inhibitors. Consistent with these findings, treatment of CRPC cells with EZH2 inhibitors dramatically enhances their sensitivity to genotoxic stress. These studies reveal a previously unappreciated mechanism of action of EZH2 inhibitors and provide a mechanistic basis for potential combination cancer therapies.

An Enhancer-Driven Stem Cell-Like Program Mediated by SOX9 Blocks Intestinal Differentiation in Colorectal Cancer.

BACKGROUND AND AIMS: Genomic alterations that encourage stem cell activity and hinder proper maturation are central to the development of colorectal cancer (CRC). Key molecular mediators that promote these malignant properties require further elucidation to galvanize translational advances. We therefore aimed to characterize a key factor that blocks intestinal differentiation, define its transcriptional and epigenetic program, and provide preclinical evidence for therapeutic targeting in CRC. METHODS: Intestinal tissue from transgenic mice and patients were analyzed by means of histopathology and immunostaining. Human CRC cells and neoplastic murine organoids were genetically manipulated for functional studies. Gene expression profiling was obtained through RNA sequencing. Histone modifications and transcription factor binding were determined with the use of chromatin immunoprecipitation sequencing. RESULTS: We demonstrate that SRY-box transcription factor 9 (SOX9) promotes CRC by activating a stem cell-like program that hinders intestinal differentiation. Intestinal adenomas and colorectal adenocarcinomas from mouse models and patients demonstrate ectopic and elevated expression of SOX9. Functional experiments indicate a requirement for SOX9 in human CRC cell lines and engineered neoplastic organoids. Disrupting SOX9 activity impairs primary CRC tumor growth by inducing intestinal differentiation. By binding to genome wide enhancers, SOX9 directly activates genes associated with Paneth and stem cell activity, including prominin 1 (PROM1). SOX9 up-regulates PROM1 via a Wnt-responsive intronic enhancer. A pentaspan transmembrane protein, PROM1 uses its first intracellular domain to support stem cell signaling, at least in part through SOX9, reinforcing a PROM1-SOX9 positive feedback loop. CONCLUSIONS: These studies establish SOX9 as a central regulator of an enhancer-driven stem cell-like program and carry important implications for developing therapeutics directed at overcoming differentiation defects in CRC.

Pan-ERBB kinase inhibition augments CDK4/6 inhibitor efficacy in oesophageal squamous cell carcinoma.

OBJECTIVE: Oesophageal squamous cell carcinoma (OSCC), like other squamous carcinomas, harbour highly recurrent cell cycle pathway alterations, especially hyperactivation of the CCND1/CDK4/6 axis, raising the potential for use of existing CDK4/6 inhibitors in these cancers. Although CDK4/6 inhibition has shown striking success when combined with endocrine therapy in oestrogen receptor positive breast cancer, CDK4/6 inhibitor palbociclib monotherapy has not revealed evidence of efficacy to date in OSCC clinical studies. Herein, we sought to elucidate the identification of key dependencies in OSCC as a foundation for the selection of targets whose blockade could be combined with CDK4/6 inhibition. DESIGN: We combined large-scale genomic dependency and pharmaceutical screening datasets with preclinical cell line models, to identified potential combination therapies in squamous cell cancer. RESULTS: We identified sensitivity to inhibitors to the ERBB family of receptor kinases, results clearly extending beyond the previously described minority of tumours with EGFR amplification/dependence, specifically finding a subset of OSCCs with dual dependence on ERBB3 and ERBB2. Subsequently. we demonstrated marked efficacy of combined pan-ERBB and CDK4/6 inhibition in vitro and in vivo. Furthermore, we demonstrated that squamous lineage transcription factor KLF5 facilitated activation of ERBBs in OSCC. CONCLUSION: These results provide clear rationale for development of combined ERBB and CDK4/6 inhibition in these cancers and raises the potential for KLF5 expression as a candidate biomarker to guide the use of these agents. These data suggested that by combining existing Food and Drug Administration (FDA)-approved agents, we have the capacity to improve therapy for OSCC and other squamous cancer.

Non-muscle-invasive micropapillary bladder cancer has a distinct lncRNA profile associated with unfavorable prognosis.

BACKGROUND: Molecular subtyping of bladder cancer has revealed luminal tumors generally have a more favourable prognosis. However, some aggressive forms of variant histology, including micropapillary, are often classified luminal. In previous work, we found long non-coding RNA (lncRNA) expression profiles could identify a subgroup of luminal bladder tumors with less aggressive biology and better outcomes. OBJECTIVE: In the present study, we aimed to investigate whether lncRNA expression profiles could identify high-grade T1 micropapillary bladder cancer with differential outcome. DESIGN, SETTING, AND PARTICIPANTS: LncRNAs were quantified from RNA-seq data from a HGT1 bladder cancer cohort that was enriched for primary micropapillary cases (15/84). Unsupervised consensus clustering of variant lncRNAs identified a three-cluster solution, which was further characterised using a panel of micropapillary-associated biomarkers, molecular subtypes, gene signatures, and survival analysis. A single-sample genomic signature was trained using lasso-penalized logistic regression to classify micropapillary-like gene-expression, as characterised by lncRNA clustering. The genomic classifier (GC) was tested on luminal tumors derived from the TCGA cohort (N = 202). OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: Patient and tumor characteristics were compared between subgroups by using X2 tests and two-sided Wilcoxon rank-sum tests. Primary endpoints were overall, progression-free and high-grade recurrence-free survival, calculated as the date of high-grade T1 disease at TURBT till date of death from any cause, progression, or recurrence, respectively. Survival rates were estimated using weighted Kaplan-Meier (KM) curves. RESULTS AND LIMITATIONS: Primary micropapillary HGT1 showed decreased FGFR3, SHH, and p53 pathway activity relative to tumors with conventional urothelial carcinoma. Many bladder cancer-associated lncRNAs were downregulated in micropapillary tumors, including UCA1, LINC00152, and MALAT1. Unsupervised consensus clustering resulted in a lncRNA cluster 1 (LC1) with worse prognosis that was enriched for primary micropapillary histology and the Luminal Unstable (LumU) molecular subtype. Interestingly, LC1 appeared to better identify aggressive HGT1 disease, compared to stratifying outcomes using primary histologic characteristics. A signature trained to identify LC1 cases showed good performance in the testing cohort, identifying seven cases with significantly worse survival (p < 0.001). Limitations include the retrospective nature of the study and the lack of a validation cohort. CONCLUSIONS: Using the lncRNA transcriptome we identified a subgroup of aggressive HGT1 bladder cancer that was enriched with micropapillary histology. These data suggest that lncRNAs can facilitate the identification of aggressive micropapillary-like tumors, potentially improving patient management.

Response and resistance to BET bromodomain inhibitors in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is a heterogeneous and clinically aggressive disease for which there is no targeted therapy. BET bromodomain inhibitors, which have shown efficacy in several models of cancer, have not been evaluated in TNBC. These inhibitors displace BET bromodomain proteins such as BRD4 from chromatin by competing with their acetyl-lysine recognition modules, leading to inhibition of oncogenic transcriptional programs. Here we report the preferential sensitivity of TNBCs to BET bromodomain inhibition in vitro and in vivo, establishing a rationale for clinical investigation and further motivation to understand mechanisms of resistance. In paired cell lines selected for acquired resistance to BET inhibition from previously sensitive TNBCs, we failed to identify gatekeeper mutations, new driver events or drug pump activation. BET-resistant TNBC cells remain dependent on wild-type BRD4, which supports transcription and cell proliferation in a bromodomain-independent manner. Proteomic studies of resistant TNBC identify strong association with MED1 and hyper-phosphorylation of BRD4 attributable to decreased activity of PP2A, identified here as a principal BRD4 serine phosphatase. Together, these studies provide a rationale for BET inhibition in TNBC and present mechanism-based combination strategies to anticipate clinical drug resistance.

Estrogen receptor signaling is reprogrammed during breast tumorigenesis.

Limited knowledge of the changes in estrogen receptor (ER) signaling during the transformation of the normal mammary gland to breast cancer hinders the development of effective prevention and treatment strategies. Differences in estrogen signaling between normal human primary breast epithelial cells and primary breast tumors obtained immediately following surgical excision were explored. Transcriptional profiling of normal ER+ mature luminal mammary epithelial cells and ER+ breast tumors revealed significant difference in the response to estrogen stimulation. Consistent with these differences in gene expression, the normal and tumor ER cistromes were distinct and sufficient to segregate normal breast tissues from breast tumors. The selective enrichment of the DNA binding motif GRHL2 in the breast cancer-specific ER cistrome suggests that it may play a role in the differential function of ER in breast cancer. Depletion of GRHL2 resulted in altered ER binding and differential transcriptional responses to estrogen stimulation. Furthermore, GRHL2 was demonstrated to be essential for estrogen-stimulated proliferation of ER+ breast cancer cells. DLC1 was also identified as an estrogen-induced tumor suppressor in the normal mammary gland with decreased expression in breast cancer. In clinical cohorts, loss of DLC1 and gain of GRHL2 expression are associated with ER+ breast cancer and are independently predictive for worse survival. This study suggests that normal ER signaling is lost and tumor-specific ER signaling is gained during breast tumorigenesis. Unraveling these changes in ER signaling during breast cancer progression should aid the development of more effective prevention strategies and targeted therapeutics.

Subependymal giant cell astrocytomas are characterized by mTORC1 hyperactivation, a very low somatic mutation rate, and a unique gene expression profile.

Subependymal giant-cell astrocytomas (SEGAs) are slow-growing brain tumors that are a hallmark feature seen in 5-10% of patients with Tuberous Sclerosis Complex (TSC). Though histologically benign, they can cause serious neurologic symptoms, leading to death if untreated. SEGAs consistently show biallelic loss of TSC1 or TSC2. Herein, we aimed to define other somatic events beyond TSC1/TSC2 loss and identify potential transcriptional drivers that contribute to SEGA formation. Paired tumor-normal whole-exome sequencing was performed on 21 resected SEGAs from 20 TSC patients. Pathogenic variants in TSC1/TSC2 were identified in 19/21 (90%) SEGAs. Copy neutral loss of heterozygosity (size range: 2.2-46 Mb) was seen in 76% (16/21) of SEGAs (44% chr9q and 56% chr16p). An average of 1.4 other somatic variants (range 0-7) per tumor were identified, unlikely of pathogenic significance. Whole transcriptome RNA-sequencing analyses revealed 190 common differentially expressed genes in SEGA (n = 16, 13 from a prior study) in pairwise comparison to each of: low grade diffuse gliomas (n = 530) and glioblastoma (n = 171) from The Cancer Genome Atlas (TCGA) consortium, ganglioglioma (n = 10), TSC cortical tubers (n = 15), and multiple normal tissues. Among these, homeobox transcription factors (TFs) HMX3, HMX2, VAX1, SIX3; and TFs IRF6 and EOMES were all expressed >12-fold higher in SEGAs (FDR/q-value < 0.05). Immunohistochemistry supported the specificity of IRF6, VAX1, SIX3 for SEGAs in comparison to other tumor entities and normal brain. We conclude that SEGAs have an extremely low somatic mutation rate, suggesting that TSC1/TSC2 loss is sufficient to drive tumor growth. The unique and highly expressed SEGA-specific TFs likely reflect the neuroepithelial cell of origin, and may also contribute to the transcriptional and epigenetic state that enables SEGA growth following two-hit loss of TSC1 or TSC2 and mTORC1 activation.

Broadly permissive intestinal chromatin underlies lateral inhibition and cell plasticity.

Cells differentiate when transcription factors bind accessible cis-regulatory elements to establish specific gene expression programs. In differentiating embryonic stem cells, chromatin at lineage-restricted genes becomes sequentially accessible, probably by means of 'pioneer' transcription factor activity, but tissues may use other strategies in vivo. Lateral inhibition is a pervasive process in which one cell forces a different identity on its neighbours, and it is unclear how chromatin in equipotent progenitors undergoing lateral inhibition quickly enables distinct, transiently reversible cell fates. Here we report the chromatin and transcriptional underpinnings of differentiation in mouse small intestine crypts, where notch signalling mediates lateral inhibition to assign progenitor cells into absorptive or secretory lineages. Transcript profiles in isolated LGR5(+) intestinal stem cells and secretory and absorptive progenitors indicated that each cell population was distinct and the progenitors specified. Nevertheless, secretory and absorptive progenitors showed comparable levels of H3K4me2 and H3K27ac histone marks and DNase I hypersensitivity--signifying accessible, permissive chromatin-at most of the same cis-elements. Enhancers acting uniquely in progenitors were well demarcated in LGR5(+) intestinal stem cells, revealing early priming of chromatin for divergent transcriptional programs, and retained active marks well after lineages were specified. On this chromatin background, ATOH1, a secretory-specific transcription factor, controls lateral inhibition through delta-like notch ligand genes and also drives the expression of numerous secretory lineage genes. Depletion of ATOH1 from specified secretory cells converted them into functional enterocytes, indicating prolonged responsiveness of marked enhancers to the presence or absence of a key transcription factor. Thus, lateral inhibition and intestinal crypt lineage plasticity involve interaction of a lineage-restricted transcription factor with broadly permissive chromatin established in multipotent stem cells.

Embryonic transcription factor SOX9 drives breast cancer endocrine resistance.

The estrogen receptor (ER) drives the growth of most luminal breast cancers and is the primary target of endocrine therapy. Although ER blockade with drugs such as tamoxifen is very effective, a major clinical limitation is the development of endocrine resistance especially in the setting of metastatic disease. Preclinical and clinical observations suggest that even following the development of endocrine resistance, ER signaling continues to exert a pivotal role in tumor progression in the majority of cases. Through the analysis of the ER cistrome in tamoxifen-resistant breast cancer cells, we have uncovered a role for an RUNX2-ER complex that stimulates the transcription of a set of genes, including most notably the stem cell factor SOX9, that promote proliferation and a metastatic phenotype. We show that up-regulation of SOX9 is sufficient to cause relative endocrine resistance. The gain of SOX9 as an ER-regulated gene associated with tamoxifen resistance was validated in a unique set of clinical samples supporting the need for the development of improved ER antagonists.

VIPER: Visualization Pipeline for RNA-seq, a Snakemake workflow for efficient and complete RNA-seq analysis.

BACKGROUND: RNA sequencing has become a ubiquitous technology used throughout life sciences as an effective method of measuring RNA abundance quantitatively in tissues and cells. The increase in use of RNA-seq technology has led to the continuous development of new tools for every step of analysis from alignment to downstream pathway analysis. However, effectively using these analysis tools in a scalable and reproducible way can be challenging, especially for non-experts. RESULTS: Using the workflow management system Snakemake we have developed a user friendly, fast, efficient, and comprehensive pipeline for RNA-seq analysis. VIPER (Visualization Pipeline for RNA-seq analysis) is an analysis workflow that combines some of the most popular tools to take RNA-seq analysis from raw sequencing data, through alignment and quality control, into downstream differential expression and pathway analysis. VIPER has been created in a modular fashion to allow for the rapid incorporation of new tools to expand the capabilities. This capacity has already been exploited to include very recently developed tools that explore immune infiltrate and T-cell CDR (Complementarity-Determining Regions) reconstruction abilities. The pipeline has been conveniently packaged such that minimal computational skills are required to download and install the dozens of software packages that VIPER uses. CONCLUSIONS: VIPER is a comprehensive solution that performs most standard RNA-seq analyses quickly and effectively with a built-in capacity for customization and expansion.

Integration of multiethnic fine-mapping and genomic annotation to prioritize candidate functional SNPs at prostate cancer susceptibility regions.

Interpretation of biological mechanisms underlying genetic risk associations for prostate cancer is complicated by the relatively large number of risk variants (n = 100) and the thousands of surrogate SNPs in linkage disequilibrium. Here, we combined three distinct approaches: multiethnic fine-mapping, putative functional annotation (based upon epigenetic data and genome-encoded features), and expression quantitative trait loci (eQTL) analyses, in an attempt to reduce this complexity. We examined 67 risk regions using genotyping and imputation-based fine-mapping in populations of European (cases/controls: 8600/6946), African (cases/controls: 5327/5136), Japanese (cases/controls: 2563/4391) and Latino (cases/controls: 1034/1046) ancestry. Markers at 55 regions passed a region-specific significance threshold (P-value cutoff range: 3.9 × 10(-4)-5.6 × 10(-3)) and in 30 regions we identified markers that were more significantly associated with risk than the previously reported variants in the multiethnic sample. Novel secondary signals (P < 5.0 × 10(-6)) were also detected in two regions (rs13062436/3q21 and rs17181170/3p12). Among 666 variants in the 55 regions with P-values within one order of magnitude of the most-associated marker, 193 variants (29%) in 48 regions overlapped with epigenetic or other putative functional marks. In 11 of the 55 regions, cis-eQTLs were detected with nearby genes. For 12 of the 55 regions (22%), the most significant region-specific, prostate-cancer associated variant represented the strongest candidate functional variant based on our annotations; the number of regions increased to 20 (36%) and 27 (49%) when examining the 2 and 3 most significantly associated variants in each region, respectively. These results have prioritized subsets of candidate variants for downstream functional evaluation.

Refined DNase-seq protocol and data analysis reveals intrinsic bias in transcription factor footprint identification.

Sequencing of DNase I hypersensitive sites (DNase-seq) is a powerful technique for identifying cis-regulatory elements across the genome. We studied the key experimental parameters to optimize performance of DNase-seq. Sequencing short fragments of 50-100 base pairs (bp) that accumulate in long internucleosome linker regions was more efficient for identifying transcription factor binding sites compared to sequencing longer fragments. We also assessed the potential of DNase-seq to predict transcription factor occupancy via generation of nucleotide-resolution transcription factor footprints. In modeling the sequence-specific DNase I cutting bias, we found a strong effect that varied over more than two orders of magnitude. This indicates that the nucleotide-resolution cleavage patterns at many transcription factor binding sites are derived from intrinsic DNase I cleavage bias rather than from specific protein-DNA interactions. In contrast, quantitative comparison of DNase I hypersensitivity between states can predict transcription factor occupancy associated with particular biological perturbations.

ChiLin: a comprehensive ChIP-seq and DNase-seq quality control and analysis pipeline.

BACKGROUND: Transcription factor binding, histone modification, and chromatin accessibility studies are important approaches to understanding the biology of gene regulation. ChIP-seq and DNase-seq have become the standard techniques for studying protein-DNA interactions and chromatin accessibility respectively, and comprehensive quality control (QC) and analysis tools are critical to extracting the most value from these assay types. Although many analysis and QC tools have been reported, few combine ChIP-seq and DNase-seq data analysis and quality control in a unified framework with a comprehensive and unbiased reference of data quality metrics. RESULTS: ChiLin is a computational pipeline that automates the quality control and data analyses of ChIP-seq and DNase-seq data. It is developed using a flexible and modular software framework that can be easily extended and modified. ChiLin is ideal for batch processing of many datasets and is well suited for large collaborative projects involving ChIP-seq and DNase-seq from different designs. ChiLin generates comprehensive quality control reports that include comparisons with historical data derived from over 23,677 public ChIP-seq and DNase-seq samples (11,265 datasets) from eight literature-based classified categories. To the best of our knowledge, this atlas represents the most comprehensive ChIP-seq and DNase-seq related quality metric resource currently available. These historical metrics provide useful heuristic quality references for experiment across all commonly used assay types. Using representative datasets, we demonstrate the versatility of the pipeline by applying it to different assay types of ChIP-seq data. The pipeline software is available open source at https://github.com/cfce/chilin . CONCLUSION: ChiLin is a scalable and powerful tool to process large batches of ChIP-seq and DNase-seq datasets. The analysis output and quality metrics have been structured into user-friendly directories and reports. We have successfully compiled 23,677 profiles into a comprehensive quality atlas with fine classification for users.