In Vivo Chemical Probing for G-Quadruplex Formation.

While DNA inside the cells is predominantly canonical right-handed double helix, guanine-rich DNAs have potential to fold into four-stranded structures that contain stacks of G-quartets (G4 DNA quadruplex). Genome sequencing has revealed G4 sequences tend to localize at the gene control regions, especially in the promoters of oncogenes. A growing body of evidence indicates that G4 DNA quadruplexes might have important regulatory roles in genome function, highlighting the need for techniques to detect genome-wide folding of DNA into this structure. Potassium permanganate in vivo treatment of cells results in oxidizing of nucleotides in single-stranded DNA regions that accompany G4 DNA quadruplexes formation, providing an excellent probe for the conformational state of DNA inside the living cells. Here, we describe a permanganate-based methodology to detect G4 DNA quadruplex, genome-wide. This methodology combined with high-throughput sequencing provides a snapshot of the DNA conformation over the whole genome in vivo.

Carbonic anhydrase 9 expression is associated with poor prognosis, tumor proliferation, and radiosensitivity of thymic carcinomas.

INTRODUCTION: Thymic epithelial tumors (TETs) comprise several histologies of thymoma and thymic carcinomas (TCs), and TC frequently metastasizes and causes death. We therefore aimed here to identify key molecules closely related to prognosis and their biological roles in high-risk TETs, particularly TCs. RESULTS: RNA sequence analysis demonstrated that hypoxia-related genes were highly expressed in TETs. The expression of the hypoxia-related gene CA9 was noteworthy, particularly in TCs. Immunohistochemical analysis revealed that CA9 was expressed in 81.0% of TCs and 20.7% of all TET samples. CA9 expression was significantly associated with Masaoka stage, WHO classification, and recurrence-free survival after tumor resection (P = 0.005). The down-regulation of CA9 transcription in TC cell lines by small interfering RNAs significantly inhibited CA9 expression, which inhibited proliferation and increased sensitivity to irradiation. CONCLUSIONS: CA9 expression may serve as a significant prognostic marker of TETs and therefore represents a potential target for the development of novel drugs and radiation-sensitizing therapy designed to improve the outcomes of patients with TCs. MATERIALS AND METHODS: We performed comprehensive transcriptome sequencing of 23 TETs and physiologic thymic specimens to identify genes highly and specifically expressed in high-risk TETs, particulary TCs. We performed immunohistochemical analysis of 179 consecutive surgically resected TETs to evaluate the significance of the association of protein expression with clinicopathological features and prognosis. The biological significance of the most promising prognostic marker was further studied using the TC cell lines, Ty-82 and MP57.

Carboxypeptidase A4 accumulation is associated with an aggressive phenotype and poor prognosis in triple-negative breast cancer.

Using whole transcriptome analysis and a lentiviral short hairpin RNA screening library, carboxypeptidase A4 (CPA4) was identified as a novel marker in breast cancer and a therapeutic target in triple­negative breast cancer (TNBC) in the present study. Immunohistochemistry was used to evaluate the presence of CPA4, estrogen receptor, progesterone receptor, human epidermal growth factor receptor 2, Ki67, epidermal growth factor receptor, cytokeratin 5/6, aldehyde dehydrogenase 1, cluster of differentiation (CD)44, CD24, claudins, E­cadherin, vimentin and androgen receptor in 221 cases of breast cancer, including 68 TNBC cases. The effects of CPA4 on the viability and migration ability of TNBC cells were analyzed using RNA interference methods. Increased CPA4 expression, specifically in the cytoplasm of cancer tissue cells, was detected. Furthermore, high CPA4 expression in TNBC cases was associated with low expression of E­cadherin and with the expression of cancer stem cell markers (high CD44/low CD24). Patients with TNBC and high levels of CPA4 expression had a significantly poorer prognosis compared with those with low CPA4 expression. Notably, viability and migration were reduced, but E­cadherin expression was upregulated in CPA4­suppressed TNBC cells. The present data suggested that CPA4 may be a novel inducer for epithelial­mesenchymal transition, which is characterized by the downregulation of E­cadherin and mesenchymal phenotypes. To conclude, CPA4 may be a marker for poor prognosis and a promising therapeutic target in TNBC with aggressive phenotypes.

Mutant TP53 modulates metastasis of triple negative breast cancer through adenosine A2b receptor signaling.

PURPOSE: The identification of genes with synthetic lethality in the context of mutant TP53 is a promising strategy for the treatment of basal-like triple negative breast cancer (TNBC). This study investigated regulators of mutant TP53 (R248Q) in basal-like TNBC and their impact on tumorigenesis. EXPERIMENTAL DESIGN: TNBC cells were analyzed by RNA-seq, and synthetic-lethal shRNA knock-down screening, to identify genes related to the expression of mutant TP53. A tissue microarray of 232 breast cancer samples, that included 66 TNBC cases, was used to assess clinicopathological correlates of tumor protein expression. Functional assays were performed in vitro and in vivo to assess the role of ADORA2B in TNBC. RESULTS: Transcriptome profiling identified ADORA2B as up-regulated in basal-like TNBC cell lines with R248Q-mutated TP53, with shRNA-screening suggesting the potential for a synthetic-lethal interaction between these genes. In clinical samples, ADORA2B was highly expressed in 39.4% (26/66) of TNBC patients. ADORA2B-expression was significantly correlated with ER (P < 0.01), PgR (P = 0.027), EGFR (P < 0.01), and tumor size (P = 0.037), and was an independent prognostic factor for outcome (P = 0.036). In line with this, ADORA2B-transduced TNBC cells showed increased tumorigenesis, and ADORA2B knockdown, along with mutant p53 knockdown, decreased metastasis both in vitro and in vivo. Notably, the cytotoxic cyclic peptide SA-I suppressed ADORA2B expression and tumorigenesis in TNBC cell lines. CONCLUSIONS: ADORA2B expression increases the oncogenic potential of basal-like TNBC and is an independent factor for poor outcome. These data suggest that ADORA2B could serve as a prognostic biomarker and a potential therapeutic target for basal-like TNBC.

Caspase14 expression is associated with triple negative phenotypes and cancer stem cell marker expression in breast cancer patients.

BACKGROUND AND OBJECTIVES: The Caspase14 (CASP14) was reported that the low expression of CASP14 in ovarian cancer and colon cancer was associated with cancer progression, on the other hand, that the CASP14 expression in breast cancer was higher than that of non-cancerous tissues. The purpose of this study is to determine the clinical significance of CASP14 in breast cancer. METHODS: We performed immunohistochemistry for CASP14, ER, PgR, HER2, Ki67, EGFR, CK5/6, CD44, CD24, ALDH1, claudins, and androgen receptor in 222 breast cancer patients including 55 TNBC cases, and evaluated the relationship of CASP14, above-mentioned markers, and prognosis. Using public microarray database of breast cancer, the prognostic value of CASP14 was calculated. RESULTS: High CASP14 expression was significantly associated with TNBC subtype (P = 0.015), nuclear grade (P = 0.006), Ki67, EGFR (P < 0.001, P = 0.016), ALDH1, CD44 and CD24 (P < 0.001, P < 0.001, P = 0.001) in 222 breast cancer cases, and the high expression of claudin1 (P = 0.017), and androgen receptor (P = 0.002) in TNBC cases was related to the high CASP14. According to the public database, survival in the high CASP14 breast cancer patients was shorter than low CASP14 patients. CONCLUSIONS: High CASP14 expression is a marker of breast cancer aggressiveness in association with proliferation, TNBC phenotype, and cancer stemness.

Overexpression of B-cell lymphoma 6 alters gene expression profile in a myeloma cell line and is associated with decreased DNA damage response.

B-cell lymphoma 6 (BCL6) attenuates DNA damage response (DDR) through gene repression and facilitates tolerance to genomic instability during immunoglobulin affinity maturation in germinal center (GC) B cells. Although BCL6 expression is repressed through normal differentiation of GC B cells into plasma cells, a recent study showed the ectopic expression of BCL6 in primary multiple myeloma (MM) cells. However, the functional roles of BCL6 in MM cells are largely unknown. Here, we report that overexpression of BCL6 in a MM cell line, KMS12PE, induced transcriptional repression of ataxia telangiectasia mutated (ATM), a DDR signaling kinase, which was associated with a reduction in γH2AX formation after DNA damage. In contrast, transcription of known targets of BCL6 in GC B cells was not affected, suggesting a cell type-specific function of BCL6. To further investigate the effects of BCL6 overexpression on the MM cell line, we undertook mRNA sequence analysis and found an upregulation in the genomic mutator activation-induced cytidine deaminase (AID) with alteration in the gene expression profile, which is suggestive of de-differentiation from plasma cells. Moreover, interleukin-6 exposure to KMS12PE led to upregulation of BCL6 and AID, downregulation of ATM, and attenuation of DDR, which were consistent with the effects of BCL6 overexpression in this MM cell line. Taken together, these results indicated that overexpression of BCL6 alters gene expression profile and confers decreased DDR in MM cells. This phenotypic change could be reproduced by interleukin-6 stimulation, suggesting an important role of external stimuli in inducing genomic instability, which is a hallmark of MM cells.

Permanganate/S1 Nuclease Footprinting Reveals Non-B DNA Structures with Regulatory Potential across a Mammalian Genome.

DNA in cells is predominantly B-form double helix. Though certain DNA sequences in vitro may fold into other structures, such as triplex, left-handed Z form, or quadruplex DNA, the stability and prevalence of these structures in vivo are not known. Here, using computational analysis of sequence motifs, RNA polymerase II binding data, and genome-wide potassium permanganate-dependent nuclease footprinting data, we map thousands of putative non-B DNA sites at high resolution in mouse B cells. Computational analysis associates these non-B DNAs with particular structures and indicates that they form at locations compatible with an involvement in gene regulation. Further analyses support the notion that non-B DNA structure formation influences the occupancy and positioning of nucleosomes in chromatin. These results suggest that non-B DNAs contribute to the control of a variety of critical cellular and organismal processes.

HMGN proteins modulate chromatin regulatory sites and gene expression during activation of naïve B cells.

The activation of naïve B lymphocyte involves rapid and major changes in chromatin organization and gene expression; however, the complete repertoire of nuclear factors affecting these genomic changes is not known. We report that HMGN proteins, which bind to nucleosomes and affect chromatin structure and function, co-localize with, and maintain the intensity of DNase I hypersensitive sites genome wide, in resting but not in activated B cells. Transcription analyses of resting and activated B cells from wild-type and Hmgn(-/-) mice, show that loss of HMGNs dampens the magnitude of the transcriptional response and alters the pattern of gene expression during the course of B-cell activation; defense response genes are most affected at the onset of activation. Our study provides insights into the biological function of the ubiquitous HMGN chromatin binding proteins and into epigenetic processes that affect the fidelity of the transcriptional response during the activation of B cell lymphocytes.

APOBEC3B high expression status is associated with aggressive phenotype in Japanese breast cancers.

BACKGROUND: The members of AID/APOBEC protein family possess cytidine deaminase activity that converts cytidine residue to uridine on DNA and RNA. Recent studies have shown the possible influence of APOBEC3B (A3B) as DNA mutators of breast cancer genome. However, the clinical significance of A3B expression in Japanese breast cancer has not been studied in detail. METHODS: Ninety-three primary breast cancer tissues (74 estrogen-receptor (ER) positive, 3 ER and HER2 positive, 6 HER2 positive, and 10 triple negative) including 37 tumor-normal pairs were assessed for A3B mRNA expression using quantitative real-time RT-PCR. We analyzed the relation between A3B expression, mutation analysis of TP53 and PIK3CA by direct sequencing, polymorphic A3B deletion allele and human papillomavirus (HPV) infection in tumors. RESULTS: A3B mRNA was overexpressed in tumors compared with normal tissue. Patients with high A3B expression were associated with subtype and progression of lymph node metastasis and pathological nuclear grade. However, the expression was not related to any other clinicopathological factors, including mutation of TP53 and PIK3CA, polymorphic A3B deletion allele, HPV infection and survival time. CONCLUSION: The expression of A3B in breast cancer was higher than in non-cancerous tissues and was related to the lymph node metastasis and nuclear grade, which are reliable aggressive phenotype markers in breast cancer. Evaluation of A3B expression in tumor may be a marker for breast cancer with malignant potential.

RAG Represents a Widespread Threat to the Lymphocyte Genome.

The RAG1 endonuclease, together with its cofactor RAG2, is essential for V(D)J recombination but is a potent threat to genome stability. The sources of RAG1 mis-targeting and the mechanisms that have evolved to suppress it are poorly understood. Here, we report that RAG1 associates with chromatin at thousands of active promoters and enhancers in the genome of developing lymphocytes. The mouse and human genomes appear to have responded by reducing the abundance of "cryptic" recombination signals near RAG1 binding sites. This depletion operates specifically on the RSS heptamer, whereas nonamers are enriched at RAG1 binding sites. Reversing this RAG-driven depletion of cleavage sites by insertion of strong recombination signals creates an ectopic hub of RAG-mediated V(D)J recombination and chromosomal translocations. Our findings delineate rules governing RAG binding in the genome, identify areas at risk of RAG-mediated damage, and highlight the evolutionary struggle to accommodate programmed DNA damage in developing lymphocytes.

Interactome maps of mouse gene regulatory domains reveal basic principles of transcriptional regulation.

A key finding of the ENCODE project is that the enhancer landscape of mammalian cells undergoes marked alterations during ontogeny. However, the nature and extent of these changes are unclear. As part of the NIH Mouse Regulome Project, we here combined DNaseI hypersensitivity, ChIP-seq, and ChIA-PET technologies to map the promoter-enhancer interactomes of pluripotent ES cells and differentiated B lymphocytes. We confirm that enhancer usage varies widely across tissues. Unexpectedly, we find that this feature extends to broadly transcribed genes, including Myc and Pim1 cell-cycle regulators, which associate with an entirely different set of enhancers in ES and B cells. By means of high-resolution CpG methylomes, genome editing, and digital footprinting, we show that these enhancers recruit lineage-determining factors. Furthermore, we demonstrate that the turning on and off of enhancers during development correlates with promoter activity. We propose that organisms rely on a dynamic enhancer landscape to control basic cellular functions in a tissue-specific manner.

Global regulation of promoter melting in naive lymphocytes.

Lymphocyte activation is initiated by a global increase in messenger RNA synthesis. However, the mechanisms driving transcriptome amplification during the immune response are unknown. By monitoring single-stranded DNA genome wide, we show that the genome of naive cells is poised for rapid activation. In G0, ∼90% of promoters from genes to be expressed in cycling lymphocytes are polymerase loaded but unmelted and support only basal transcription. Furthermore, the transition from abortive to productive elongation is kinetically limiting, causing polymerases to accumulate nearer to transcription start sites. Resting lymphocytes also limit the expression of the transcription factor IIH complex, including XPB and XPD helicases involved in promoter melting and open complex extension. To date, two rate-limiting steps have been shown to control global gene expression in eukaryotes: preinitiation complex assembly and polymerase pausing. Our studies identify promoter melting as a third key regulatory step and propose that this mechanism ensures a prompt lymphocyte response to invading pathogens.

A genome-wide map of CTCF multivalency redefines the CTCF code.

The "CTCF code" hypothesis posits that CTCF pleiotropic functions are driven by recognition of diverse sequences through combinatorial use of its 11 zinc fingers (ZFs). This model, however, is supported by in vitro binding studies of a limited number of sequences. To study CTCF multivalency in vivo, we define ZF binding requirements at ∼50,000 genomic sites in primary lymphocytes. We find that CTCF reads sequence diversity through ZF clustering. ZFs 4-7 anchor CTCF to ∼80% of targets containing the core motif. Nonconserved flanking sequences are recognized by ZFs 1-2 and ZFs 8-11 clusters, which also stabilize CTCF broadly. Alternatively, ZFs 9-11 associate with a second phylogenetically conserved upstream motif at ∼15% of its sites. Individually, ZFs increase overall binding and chromatin residence time. Unexpectedly, we also uncovered a conserved downstream DNA motif that destabilizes CTCF occupancy. Thus, CTCF associates with a wide array of DNA modules via combinatorial clustering of its 11 ZFs.

Rif1 prevents resection of DNA breaks and promotes immunoglobulin class switching.

DNA double-strand breaks (DSBs) represent a threat to the genome because they can lead to the loss of genetic information and chromosome rearrangements. The DNA repair protein p53 binding protein 1 (53BP1) protects the genome by limiting nucleolytic processing of DSBs by a mechanism that requires its phosphorylation, but whether 53BP1 does so directly is not known. Here, we identify Rap1-interacting factor 1 (Rif1) as an ATM (ataxia-telangiectasia mutated) phosphorylation-dependent interactor of 53BP1 and show that absence of Rif1 results in 5'-3' DNA-end resection in mice. Consistent with enhanced DNA resection, Rif1 deficiency impairs DNA repair in the G(1) and S phases of the cell cycle, interferes with class switch recombination in B lymphocytes, and leads to accumulation of chromosome DSBs.

53BP1 alters the landscape of DNA rearrangements and suppresses AID-induced B cell lymphoma.

Deficiencies in factors that regulate the DNA damage response enhance the incidence of malignancy by destabilizing the genome. However, the precise influence of the DNA damage response on regulation of cancer-associated rearrangements is not well defined. Here we examine the genome-wide impact of tumor protein P53-binding protein 1 (53BP1) deficiency in lymphoma and translocation. While both activation-induced cytidine deaminase (AID) and 53BP1 have been associated with cancer in humans, neither AID overexpression nor loss of 53BP1 is sufficient to produce malignancy. However, the combination of 53BP1 deficiency and AID deregulation results in B cell lymphoma. Deep sequencing of the genome of 53BP1(-/-) cancer cells and translocation capture sequencing (TC-Seq) of primary 53BP1(-/-) B cells revealed that their chromosomal rearrangements differ from those found in wild-type cells in that they show increased DNA end resection. Moreover, loss of 53BP1 alters the translocatome by increasing rearrangements to intergenic regions.

RPA accumulation during class switch recombination represents 5'-3' DNA-end resection during the S-G2/M phase of the cell cycle.

Activation-induced cytidine deaminase (AID) promotes chromosomal translocations by inducing DNA double-strand breaks (DSBs) at immunoglobulin (Ig) genes and oncogenes in the G1 phase. RPA is a single-stranded DNA (ssDNA)-binding protein that associates with resected DSBs in the S phase and facilitates the assembly of factors involved in homologous repair (HR), such as Rad51. Notably, RPA deposition also marks sites of AID-mediated damage, but its role in Ig gene recombination remains unclear. Here, we demonstrate that RPA associates asymmetrically with resected ssDNA in response to lesions created by AID, recombination-activating genes (RAG), or other nucleases. Small amounts of RPA are deposited at AID targets in G1 in an ATM-dependent manner. In contrast, recruitment in the S-G2/M phase is extensive, ATM independent, and associated with Rad51 accumulation. In the S-G2/M phase, RPA increases in nonhomologous-end-joining-deficient lymphocytes, where there is more extensive DNA-end resection. Thus, most RPA recruitment during class switch recombination represents salvage of unrepaired breaks by homology-based pathways during the S-G2/M phase of the cell cycle.

c-Myc is a universal amplifier of expressed genes in lymphocytes and embryonic stem cells.

The c-Myc HLH-bZIP protein has been implicated in physiological or pathological growth, proliferation, apoptosis, metabolism, and differentiation at the cellular, tissue, or organismal levels via regulation of numerous target genes. No principle yet unifies Myc action due partly to an incomplete inventory and functional accounting of Myc's targets. To observe Myc target expression and function in a system where Myc is temporally and physiologically regulated, the transcriptomes and the genome-wide distributions of Myc, RNA polymerase II, and chromatin modifications were compared during lymphocyte activation and in ES cells as well. A remarkably simple rule emerged from this quantitative analysis: Myc is not an on-off specifier of gene activity, but is a nonlinear amplifier of expression, acting universally at active genes, except for immediate early genes that are strongly induced before Myc. This rule of Myc action explains the vast majority of Myc biology observed in literature.

Zinc finger transcription factor zDC is a negative regulator required to prevent activation of classical dendritic cells in the steady state.

Classical dendritic cells (cDCs) process and present antigens to T cells. Under steady-state conditions, antigen presentation by cDCs induces tolerance. In contrast, during infection or inflammation, cDCs become activated, express higher levels of cell surface MHC molecules, and induce strong adaptive immune responses. We recently identified a cDC-restricted zinc finger transcription factor, zDC (also known as Zbtb46 or Btbd4), that is not expressed by other immune cell populations, including plasmacytoid DCs, monocytes, or macrophages. We define the zDC consensus DNA binding motif and the genes regulated by zDC using chromatin immunoprecipitation and deep sequencing. By deleting zDC from the mouse genome, we show that zDC is primarily a negative regulator of cDC gene expression. zDC deficiency alters the cDC subset composition in the spleen in favor of CD8(+) DCs, up-regulates activation pathways in steady-state cDCs, including elevated MHC II expression, and enhances cDC production of vascular endothelial growth factor leading to increased vascularization of skin-draining lymph nodes. Consistent with these observations, zDC protein expression is rapidly down-regulated after TLR stimulation. Thus, zDC is a TLR-responsive, cDC-specific transcriptional repressor that is in part responsible for preventing cDC maturation in the steady state.

Mouse model of endemic Burkitt translocations reveals the long-range boundaries of Ig-mediated oncogene deregulation.

Human Burkitt lymphomas are divided into two main clinical variants: the endemic form, affecting African children infected with malaria and the Epstein-Barr virus, and the sporadic form, distributed across the rest of the world. However, whereas sporadic translocations decapitate Myc from 5' proximal regulatory elements, most endemic events occur hundreds of kilobases away from Myc. The origin of these rearrangements and how they deregulate oncogenes at such distances remain unclear. We here recapitulate endemic Burkitt lymphoma-like translocations in plasmacytomas from uracil N-glycosylase and activation-induced cytidine deaminase-deficient mice. Mapping of translocation breakpoints using an acetylated histone H3 lysine 9 chromatin immunoprecipitation sequencing approach reveals Igh fusions up to ∼350 kb upstream of Myc or the related oncogene Mycn. A comprehensive analysis of epigenetic marks, PolII recruitment, and transcription in tumor cells demonstrates that the 3' Igh enhancer (Eα) vastly remodels ∼450 kb of chromatin into translocated sequences, leading to significant polymerase occupancy and constitutive oncogene expression. We show that this long-range epigenetic reprogramming is directly proportional to the physical interaction of Eα with translocated sites. Our studies thus uncover the extent of epigenetic remodeling by Ig 3' enhancers and provide a rationale for the long-range deregulation of translocated oncogenes in endemic Burkitt lymphomas. The data also shed light on the origin of endemic-like chromosomal rearrangements.

DNA damage defines sites of recurrent chromosomal translocations in B lymphocytes.

Recurrent chromosomal translocations underlie both haematopoietic and solid tumours. Their origin has been ascribed to selection of random rearrangements, targeted DNA damage, or frequent nuclear interactions between translocation partners; however, the relative contribution of each of these elements has not been measured directly or on a large scale. Here we examine the role of nuclear architecture and frequency of DNA damage in the genesis of chromosomal translocations by measuring these parameters simultaneously in cultured mouse B lymphocytes. In the absence of recurrent DNA damage, translocations between Igh or Myc and all other genes are directly related to their contact frequency. Conversely, translocations associated with recurrent site-directed DNA damage are proportional to the rate of DNA break formation, as measured by replication protein A accumulation at the site of damage. Thus, non-targeted rearrangements reflect nuclear organization whereas DNA break formation governs the location and frequency of recurrent translocations, including those driving B-cell malignancies.