Transposon DNA sequences facilitate the tissue-specific gene transfer of circulating tumor DNA between human cells.

The exchange of genes between cells is known to play an important physiological and pathological role in many organisms. We show that circulating tumor DNA (ctDNA) facilitates cell-specific gene transfer between human cancer cells and explain part of the mechanisms behind this phenomenon. As ctDNA migrates into the nucleus, genetic information is transferred. Cell targeting and ctDNA integration require ERVL, SINE or LINE DNA sequences. Chemically manufactured AluSp and MER11C sequences replicated multiple myeloma (MM) ctDNA cell targeting and integration. Additionally, we found that ctDNA may alter the treatment response of MM and pancreatic cancer models. This study shows that retrotransposon DNA sequences promote cancer gene transfer. However, because cell-free DNA has been detected in physiological and other pathological conditions, our findings have a broader impact than just cancer. Furthermore, the discovery that transposon DNA sequences mediate tissue-specific targeting will open up a new avenue for the delivery of genes and therapies.

Natural mutagenesis of human genomes by endogenous retrotransposons.

Two abundant classes of mobile elements, namely Alu and L1 elements, continue to generate new retrotransposon insertions in human genomes. Estimates suggest that these elements have generated millions of new germline insertions in individual human genomes worldwide. Unfortunately, current technologies are not capable of detecting most of these young insertions, and the true extent of germline mutagenesis by endogenous human retrotransposons has been difficult to examine. Here, we describe technologies for detecting these young retrotransposon insertions and demonstrate that such insertions indeed are abundant in human populations. We also found that new somatic L1 insertions occur at high frequencies in human lung cancer genomes. Genome-wide analysis suggests that altered DNA methylation may be responsible for the high levels of L1 mobilization observed in these tumors. Our data indicate that transposon-mediated mutagenesis is extensive in human genomes and is likely to have a major impact on human biology and diseases.

Genetic heterogeneity within collective invasion packs drives leader and follower cell phenotypes.

Collective invasion, the coordinated movement of cohesive packs of cells, has become recognized as a major mode of metastasis for solid tumors. These packs are phenotypically heterogeneous and include specialized cells that lead the invasive pack and others that follow behind. To better understand how these unique cell types cooperate to facilitate collective invasion, we analyzed transcriptomic sequence variation between leader and follower populations isolated from the H1299 non-small cell lung cancer cell line using an image-guided selection technique. We now identify 14 expressed mutations that are selectively enriched in leader or follower cells, suggesting a novel link between genomic and phenotypic heterogeneity within a collectively invading tumor cell population. Functional characterization of two phenotype-specific candidate mutations showed that ARP3 enhances collective invasion by promoting the leader cell phenotype and that wild-type KDM5B suppresses chain-like cooperative behavior. These results demonstrate an important role for distinct genetic variants in establishing leader and follower phenotypes and highlight the necessity of maintaining a capacity for phenotypic plasticity during collective cancer invasion.

Subtype-Independent ANP32E Reduction During Breast Cancer Progression in Accordance with Chromatin Relaxation.

BACKGROUND: Chromatin state provides a clear decipherable blueprint for maintenance of transcriptional patterns, exemplifying a mitotically stable form of cellular programming in dividing cells. In this regard, genomic studies of chromatin states within cancerous tissues have the potential to uncover novel aspects of tumor biology and unique mechanisms associated with disease phenotypes and outcomes. The degree to which chromatin state differences occur in accordance with breast cancer features has not been established. METHODS: We applied a series of unsupervised computational methods to identify chromatin and molecular differences associated with discrete physiologies across human breast cancer tumors. RESULTS: Chromatin patterns alone are capable of stratifying tumors in association with cancer subtype and disease progression. Major differences occur at DNA motifs for the transcription factor FOXA1, in hormone receptor-positive tumors, and motifs for SOX9 in Basal-like tumors. We find that one potential driver of this effect, the histone chaperone ANP32E, is inversely correlated with tumor progression and relaxation of chromatin at FOXA1 binding sites. Tumors with high levels of ANP32E exhibit an immune response and proliferative gene expression signature, whereas tumors with low ANP32E levels appear programmed for differentiation. CONCLUSIONS: Our results indicate that ANP32E may function through chromatin state regulation to control breast cancer differentiation and tumor plasticity. This study sets a precedent for future computational studies of chromatin changes in carcinogenesis.

Prognostic significance of an invasive leader cell-derived mutation cluster on chromosome 16q.

BACKGROUND: Intratumoral heterogeneity is defined by subpopulations with varying genotypes and phenotypes. Specialized, highly invasive leader cells and less invasive follower cells are phenotypically distinct subpopulations that cooperate during collective cancer invasion. Because leader cells are a rare subpopulation that would be missed by bulk sequencing, a novel image-guided genomics platform was used to precisely select this subpopulation. This study identified a novel leader cell mutation signature and tested its ability to predict prognosis in non-small cell lung cancer (NSCLC) patient cohorts. METHODS: Spatiotemporal genomic and cellular analysis was used to isolate and perform RNA sequencing on leader and follower populations from the H1299 NSCLC cell line, and it revealed a leader-specific mutation cluster on chromosome 16q. Genomic data from patients with lung squamous cell carcinoma (LUSC; n = 475) and lung adenocarcinoma (LUAD; n = 501) from The Cancer Genome Atlas were stratified by 16q mutation cluster (16qMC) status (16qMC+ vs 16qMC-) and compared for overall survival (OS), progression-free survival (PFS), and gene set enrichment analysis (GSEA). RESULTS: Poorer OS, poorer PFS, or both were found across all stages and among early-stage patients with 16qMC+ tumors within the LUSC and LUAD cohorts. GSEA revealed 16qMC+ tumors to be enriched for the expression of metastasis- and survival-associated gene sets. CONCLUSIONS: This represents the first leader cell mutation signature identified in patients and has the potential to better stratify high-risk NSCLC and ultimately improve patient outcomes.

R loops: lassoing DNA methylation at CpGi.

In this issue of Molecular Cell, Ginno et al. (2012) describe unusual sequence features at promoter CpG islands that can lead to formation of persistent RNA-DNA hybrids (R loops), which are proposed to prevent genomic DNA methylation.

A hot L1 retrotransposon evades somatic repression and initiates human colorectal cancer.

Although human LINE-1 (L1) elements are actively mobilized in many cancers, a role for somatic L1 retrotransposition in tumor initiation has not been conclusively demonstrated. Here, we identify a novel somatic L1 insertion in the APC tumor suppressor gene that provided us with a unique opportunity to determine whether such insertions can actually initiate colorectal cancer (CRC), and if so, how this might occur. Our data support a model whereby a hot L1 source element on Chromosome 17 of the patient's genome evaded somatic repression in normal colon tissues and thereby initiated CRC by mutating the APC gene. This insertion worked together with a point mutation in the second APC allele to initiate tumorigenesis through the classic two-hit CRC pathway. We also show that L1 source profiles vary considerably depending on the ancestry of an individual, and that population-specific hot L1 elements represent a novel form of cancer risk.

MAX is an epigenetic sensor of 5-carboxylcytosine and is altered in multiple myeloma.

The oncogenic transcription factor MYC and its binding partner MAX regulate gene expression by binding to DNA at enhancer-box (E-box) elements 5΄-CACGTG-3΄. In mammalian genomes, the central E-box CpG has the potential to be methylated at the 5-position of cytosine (5mC), or to undergo further oxidation to the 5-hydroxymethyl (5hmC), 5-formyl (5fC), or 5-carboxyl (5caC) forms. We find that MAX exhibits the greatest affinity for a 5caC or unmodified C-containing E-box, and much reduced affinities for the corresponding 5mC, 5hmC or 5fC forms. Crystallization of MAX with a 5caC modified E-box oligonucleotide revealed that MAX Arg36 recognizes 5caC using a 5caC-Arg-Guanine triad, with the next nearest residue to the carboxylate group being Arg60. In an analysis of >800 primary multiple myelomas, MAX alterations occurred at a frequency of ∼3%, more than half of which were single nucleotide substitutions affecting a basic clamp-like interface important for DNA interaction. Among these, arginines 35, 36 and 60 were the most frequently altered. In vitro binding studies showed that whereas mutation of Arg36 (R36W) or Arg35 (R35H/L) completely abolished DNA binding, mutation of Arg60 (R60Q) significantly reduced DNA binding, but retained a preference for the 5caC modified E-box. Interestingly, MAX alterations define a subset of myeloma patients with lower MYC expression and a better overall prognosis. Together these data indicate that MAX can act as a direct epigenetic sensor of E-box cytosine modification states and that local CpG modification and MAX variants converge to modulate the MAX-MYC transcriptional network.

GC skew defines distinct RNA polymerase pause sites in CpG island promoters.

CpG islands (CGIs) are associated with over half of human gene promoters and are characterized by a unique chromatin environment and high levels of bidirectional transcriptional activity relative to surrounding genomic regions, suggesting that RNA polymerase (Pol II) progression past the CGI boundaries is restricted. Here we describe a novel transcriptional regulatory step wherein Pol II encounters an additional barrier to elongation distinct from the promoter-proximal pause and occurring at the downstream boundary of the CGI domain. For most CGI-associated promoters, Pol II exhibits a dominant pause at either the promoter-proximal or this distal site that correlates, both in position and in intensity, with local regions of high GC skew, a sequence feature known to form unique secondary structures. Upon signal-induced gene activation, long-range enhancer contacts at the dominant pause site are selectively enhanced, suggesting a new role for enhancers at the downstream pause. These data point to an additional level of control over transcriptional output at a subset of CGI-associated genes that is linked to DNA sequence and the integrity of the CGI domain.

Gene integrated set profile analysis: a context-based approach for inferring biological endpoints.

The identification of genes with specific patterns of change (e.g. down-regulated and methylated) as phenotype drivers or samples with similar profiles for a given gene set as drivers of clinical outcome, requires the integration of several genomic data types for which an 'integrate by intersection' (IBI) approach is often applied. In this approach, results from separate analyses of each data type are intersected, which has the limitation of a smaller intersection with more data types. We introduce a new method, GISPA (Gene Integrated Set Profile Analysis) for integrated genomic analysis and its variation, SISPA (Sample Integrated Set Profile Analysis) for defining respective genes and samples with the context of similar, a priori specified molecular profiles. With GISPA, the user defines a molecular profile that is compared among several classes and obtains ranked gene sets that satisfy the profile as drivers of each class. With SISPA, the user defines a gene set that satisfies a profile and obtains sample groups of profile activity. Our results from applying GISPA to human multiple myeloma (MM) cell lines contained genes of known profiles and importance, along with several novel targets, and their further SISPA application to MM coMMpass trial data showed clinical relevance.

Characterization of a Linked Jumonji Domain of the KDM5/JARID1 Family of Histone H3 Lysine 4 Demethylases.

The KDM5/JARID1 family of Fe(II)- and α-ketoglutarate-dependent demethylases remove methyl groups from tri- and dimethylated lysine 4 of histone H3. Accumulating evidence from primary tumors and model systems supports a role for KDM5A (JARID1A/RBP2) and KDM5B (JARID1B/PLU1) as oncogenic drivers. The KDM5 family is unique among the Jumonji domain-containing histone demethylases in that there is an atypical insertion of a DNA-binding ARID domain and a histone-binding PHD domain into the Jumonji domain, which separates the catalytic domain into two fragments (JmjN and JmjC). Here we demonstrate that internal deletion of the ARID and PHD1 domains has a negligible effect on in vitro enzymatic kinetics of the KDM5 family of enzymes. We present a crystal structure of the linked JmjN-JmjC domain from KDM5A, which reveals that the linked domain fully reconstitutes the cofactor (metal ion and α-ketoglutarate) binding characteristics of other structurally characterized Jumonji domain demethylases. Docking studies with GSK-J1, a selective inhibitor of the KDM6/KDM5 subfamilies, identify critical residues for binding of the inhibitor to the reconstituted KDM5 Jumonji domain. Further, we found that GSK-J1 inhibited the demethylase activity of KDM5C with 8.5-fold increased potency compared with that of KDM5B at 1 mm α-ketoglutarate. In contrast, JIB-04 (a pan-inhibitor of the Jumonji demethylase superfamily) had the opposite effect and was ~8-fold more potent against KDM5B than against KDM5C. Interestingly, the relative selectivity of JIB-04 toward KDM5B over KDM5C in vitro translates to a ~10-50-fold greater growth-inhibitory activity against breast cancer cell lines. These data define the minimal requirements for enzymatic activity of the KDM5 family to be the linked JmjN-JmjC domain coupled with the immediate C-terminal helical zinc-binding domain and provide structural characterization of the linked JmjN-JmjC domain for the KDM5 family, which should prove useful in the design of KDM5 demethylase inhibitors with improved potency and selectivity.

The Mechanisms of Generation, Recognition, and Erasure of DNA 5-Methylcytosine and Thymine Oxidations.

One of the most fundamental questions in the control of gene expression in mammals is how the patterns of epigenetic modifications of DNA are generated, recognized, and erased. This includes covalent cytosine methylation of DNA and its associated oxidation states. An array of AdoMet-dependent methyltransferases, Fe(II)- and α-ketoglutarate-dependent dioxygenases, base excision glycosylases, and sequence-specific transcription factors is responsible for changing, maintaining, and interpreting the modification status of specific regions of chromatin. This review focuses on recent developments in characterizing the functional and structural links between the modification status of two DNA bases 5-methylcytosine and thymine (5-methyluracil).

HNF4α regulates claudin-7 protein expression during intestinal epithelial differentiation.

The intestinal epithelium is a dynamic barrier that maintains the distinct environments of intestinal tissue and lumen. Epithelial barrier function is defined principally by tight junctions, which, in turn, depend on the regulated expression of claudin family proteins. Claudins are expressed differentially during intestinal epithelial cell (IEC) differentiation. However, regulatory mechanisms governing claudin expression during epithelial differentiation are incompletely understood. We investigated the molecular mechanisms regulating claudin-7 during IEC differentiation. Claudin-7 expression is increased as epithelial cells differentiate along the intestinal crypt-luminal axis. By using model IECs we observed increased claudin-7 mRNA and nascent heteronuclear RNA levels during differentiation. A screen for potential regulators of the CLDN7 gene during IEC differentiation was performed using a transcription factor/DNA binding array, CLDN7 luciferase reporters, and in silico promoter analysis. We identified hepatocyte nuclear factor 4α as a regulatory factor that bound endogenous CLDN7 promoter in differentiating IECs and stimulated CLDN7 promoter activity. These findings support a role of hepatocyte nuclear factor 4α in controlling claudin-7 expression during IEC differentiation.

Regulation of estrogen receptor alpha by the SET7 lysine methyltransferase.

Estrogen receptor alpha (ER) is a ligand-dependent transcription factor. Upon binding estrogen, ER recruits coactivator complexes with histone acetyltransferase or methyltransferase activities to activate downstream target genes. In addition to histones, coactivators can modify ER itself and other proteins in the transactivation complex. Here, we show that ER is directly methylated at lysine 302 (K302) by the SET7 methyltransferase. SET7-mediated methylation stabilizes ER and is necessary for the efficient recruitment of ER to its target genes and for their transactivation. The SET7-ER complex structure reveals the molecular basis for ER peptide recognition and predicts that modifications or mutations of nearby residues would affect K302 methylation. Indeed, a breast cancer-associated mutation at K303 (K303R) alters methylation at K302 in vitro and in vivo. These findings raise the possibility that generation, recognition, and removal of modifications within the ER hinge region generate "ER modification cassettes" that yield distinct patterns for signaling downstream events.

CRISPR/Cas9 systems have off-target activity with insertions or deletions between target DNA and guide RNA sequences.

CRISPR/Cas9 systems are a versatile tool for genome editing due to the highly efficient targeting of DNA sequences complementary to their RNA guide strands. However, it has been shown that RNA-guided Cas9 nuclease cleaves genomic DNA sequences containing mismatches to the guide strand. A better understanding of the CRISPR/Cas9 specificity is needed to minimize off-target cleavage in large mammalian genomes. Here we show that genomic sites could be cleaved by CRISPR/Cas9 systems when DNA sequences contain insertions ('DNA bulge') or deletions ('RNA bulge') compared to the RNA guide strand, and Cas9 nickases used for paired nicking can also tolerate bulges in one of the guide strands. Variants of single-guide RNAs (sgRNAs) for four endogenous loci were used as model systems, and their cleavage activities were quantified at different positions with 1- to 5-bp bulges. We further investigated 114 putative genomic off-target loci of 27 different sgRNAs and confirmed 15 off-target sites, each harboring a single-base bulge and one to three mismatches to the guide strand. Our results strongly indicate the need to perform comprehensive off-target analysis related to DNA and sgRNA bulges in addition to base mismatches, and suggest specific guidelines for reducing potential off-target cleavage.

A dual role for the histone methyltransferase PR-SET7/SETD8 and histone H4 lysine 20 monomethylation in the local regulation of RNA polymerase II pausing.

RNA polymerase II (Pol II) promoter-proximal pausing plays a critical role in postinitiation transcriptional regulation at many metazoan genes. We showed recently that histone H4 lysine 16 acetylation (H4K16Ac), mediated by the MSL complex, facilitates the release of paused Pol II. In contrast, H4 lysine 20 trimethylation (H4K20me3), mediated by SUV420H2, enforces Pol II pausing by inhibiting MSL recruitment. However, how the balance between H4K16Ac and H4K20me3 is locally regulated remains unclear. Here, we demonstrate that PR-SET7/SETD8, which monomethylates histone H4 lysine 20 (H4K20me1), controls both H4K16Ac and H4K20me3 and in doing so, regulates Pol II pausing dynamics. We find that PR-SET7-mediated H4K20me1 is necessary for the recruitment of the MSL complex, subsequent H4K16Ac, and release of Pol II into active elongation. Although dispensable for SUV420H2 recruitment, PR-SET7-mediated H4K20me1 is required for H4K20me3. Although depletion of SUV420H2 is sufficient to deplete H4K20me3 and relieve an H4K20me3-induced pause, pausing is maintained in the absence of PR-SET7 despite H4K20me3 depletion because of an inability to recruit the MSL complex in the absence of H4K20me1. These findings highlight the requirement for PR-SET7 and H4K20me1 in establishing both the H4K16Ac and H4K20me3 marks and point to a dual role in the local regulation of Pol II pausing.

A heterozygous IDH1R132H/WT mutation induces genome-wide alterations in DNA methylation.

Monoallelic point mutations of the NADP(+)-dependent isocitrate dehydrogenases IDH1 and IDH2 occur frequently in gliomas, acute myeloid leukemias, and chondromas, and display robust association with specific DNA hypermethylation signatures. Here we show that heterozygous expression of the IDH1(R132H) allele is sufficient to induce the genome-wide alterations in DNA methylation characteristic of these tumors. Using a gene-targeting approach, we knocked-in a single copy of the most frequently observed IDH1 mutation, R132H, into a human cancer cell line and profiled changes in DNA methylation at over 27,000 CpG dinucleotides relative to wild-type parental cells. We find that IDH1(R132H/WT) mutation induces widespread alterations in DNA methylation, including hypermethylation of 2010 and hypomethylation of 842 CpG loci. We demonstrate that many of these alterations are consistent with those observed in IDH1-mutant and G-CIMP+ primary gliomas and can segregate IDH wild-type and mutated tumors as well as those exhibiting the G-CIMP phenotype in unsupervised analysis of two primary glioma cohorts. Further, we show that the direction of IDH1(R132H/WT)-mediated DNA methylation change is largely dependent upon preexisting DNA methylation levels, resulting in depletion of moderately methylated loci. Additionally, whereas the levels of multiple histone H3 and H4 methylation modifications were globally increased, consistent with broad inhibition of histone demethylation, hypermethylation at H3K9 in particular accompanied locus-specific DNA hypermethylation at several genes down-regulated in IDH1(R132H/WT) knock-in cells. These data provide insight on epigenetic alterations induced by IDH1 mutations and support a causal role for IDH1(R132H/WT) mutants in driving epigenetic instability in human cancer cells.

Epigenetic changes associated with inflammation in breast cancer patients treated with chemotherapy.

Inflammation has been associated with fatigue during and after various types of breast cancer treatments. We examined whether prior chemotherapy was associated with DNA methylation patterns that could explain persisting inflammation and/or fatigue in women treated for breast cancer. Prior to breast radiation therapy, DNA was extracted from peripheral blood mononuclear cells (PBMCs) of 61 Stage 0-IIIA breast cancer patients who had received partial mastectomy with or without chemotherapy. DNA methylation was assessed at >485,000 CpG sites across the genome along with fatigue and plasma inflammatory markers previously associated with fatigue. Compared to non-chemotherapy-treated, women who had received chemotherapy exhibited significantly decreased methylation at eight CpG sites (p<1.03×10(-7)) including four in exon 11 of transmembrane protein 49 (TMEM49), which demonstrated the largest decreases in methylation. Lower methylation at each identified CpG site was associated with increased plasma soluble tumor necrosis factor receptor 2 (sTNFR2) and interleukin (IL)-6 and mediated the relationship between chemotherapy and increases in these inflammatory biomarkers adjusting for multiple clinical and treatment characteristics. sTNFR2, but not CpG methylation status, was correlated with fatigue. Six months after breast radiation therapy, DNA methylation, inflammatory biomarkers and fatigue assessments were repeated in a subset of subjects (N=39). Reduced methylation in 4 of the 8 identified CpG sites was still observed in chemotherapy versus non-chemotherapy-treated patients, albeit with some decay indicating the dynamic and potentially reversible nature of the changes. Reduced methylation in these 4 CpG sites also continued to correlate with either increased sTNFR2 or IL-6, but not fatigue. In conclusion, prior chemotherapy treatment was associated with decreased methylation of CpG sites in DNA from PBMCs of breast cancer patients, which correlated with increased inflammatory markers prior to and 6months after radiation therapy. Persisting epigenetic changes secondary to chemotherapy may be one factor that contributes to inflammation and its consequences including cancer-related fatigue in vulnerable breast cancer patients.

Recognition and potential mechanisms for replication and erasure of cytosine hydroxymethylation.

Cytosine residues in mammalian DNA occur in at least three forms, cytosine (C), 5-methylcytosine (M; 5mC) and 5-hydroxymethylcytosine (H; 5hmC). During semi-conservative DNA replication, hemi-methylated (M/C) and hemi-hydroxymethylated (H/C) CpG dinucleotides are transiently generated, where only the parental strand is modified and the daughter strand contains native cytosine. Here, we explore the role of DNA methyltransferases (DNMT) and ten eleven translocation (Tet) proteins in perpetuating these states after replication, and the molecular basis of their recognition by methyl-CpG-binding domain (MBD) proteins. Using recombinant proteins and modified double-stranded deoxyoligonucleotides, we show that DNMT1 prefers a hemi-methylated (M/C) substrate (by a factor of >60) over hemi-hydroxymethylated (H/C) and unmodified (C/C) sites, whereas both DNMT3A and DNMT3B have approximately equal activity on all three substrates (C/C, M/C and H/C). Binding of MBD proteins to methylated DNA inhibited Tet1 activity, suggesting that MBD binding may also play a role in regulating the levels of 5hmC. All five MBD proteins generally have reduced binding affinity for 5hmC relative to 5mC in the fully modified context (H/M versus M/M), though their relative abilities to distinguish the two varied considerably. We further show that the deamination product of 5hmC could be excised by thymine DNA glycosylase and MBD4 glycosylases regardless of context.

Epigenetic and inflammatory markers in older adults with cancer: A Young International Society of Geriatric Oncology narrative review.

The number of adults aged ≥ 65 years with cancer is rapidly increasing. Older adults with cancer are susceptible to treatment-related acute and chronic adverse events, resulting in loss of independence, reduction in physical function, and decreased quality of life. Nevertheless, evidence-based interventions to prevent or treat acute and chronic adverse events in older adults with cancer are limited. Several promising blood-based biomarkers related to inflammation and epigenetic modifications are available to identify older adults with cancer who are at increased risk of accelerated aging and physical, functional, and cognitive impairments caused by the cancer and its treatment. Inflammatory changes and epigenetic modifications can be reversible and targeted by lifestyle changes and interventions. Here we discuss ways in which changes in inflammatory and epigenetic pathways influence the aging process and how these pathways can be targeted by interventions aimed at reducing inflammation and aging-associated biological markers. As the number of older adults with cancer entering survivorship continues to increase, it is becoming progressively more important to understand ways in which the benefit from treatment can be enhanced while reducing the effects of accelerated aging.