Koji Tsumagari

Early de novo DNA methylation and prolonged demethylation in the muscle lineage.

Myogenic cell cultures derived from muscle biopsies are excellent models for human cell differentiation. We report the first comprehensive analysis of myogenesis-specific DNA hyper- and hypo-methylation throughout the genome for human muscle progenitor cells (both myoblasts and myotubes) and skeletal muscle tissue vs. 30 non-muscle samples using reduced representation bisulfite sequencing. We also focused on four genes with extensive hyper- or hypo-methylation in the muscle lineage (PAX3, TBX1, MYH7B/MIR499 and OBSCN) to compare DNA methylation, DNaseI hypersensitivity, histone modification, and CTCF binding profiles. We found that myogenic hypermethylation was strongly associated with homeobox or T-box genes and muscle hypomethylation with contractile fiber genes. Nonetheless, there was no simple relationship between differential gene expression and myogenic differential methylation, rather only for subsets of these genes, such as contractile fiber genes. Skeletal muscle retained ~30% of the hypomethylated sites but only ~3% of hypermethylated sites seen in myogenic progenitor cells. By enzymatic assays, skeletal muscle was 2-fold enriched globally in genomic 5-hydroxymethylcytosine (5-hmC) vs. myoblasts or myotubes and was the only sample type enriched in 5-hmC at tested myogenic hypermethylated sites in PAX3/CCDC140 andTBX1. TET1 and TET2 RNAs, which are involved in generation of 5-hmC and DNA demethylation, were strongly upregulated in myoblasts and myotubes. Our findings implicate de novo methylation predominantly before the myoblast stage and demethylation before and after the myotube stage in control of transcription and co-transcriptional RNA processing. They also suggest that, in muscle, TET1 or TET2 are involved in active demethylation and in formation of stable 5-hmC residues.

Epigenetics of a tandem DNA repeat: chromatin DNaseI sensitivity and opposite methylation changes in cancers

DNA methylation and chromatin DNaseI sensitivity were analyzed in and adjacent to D4Z4 repeat arrays, which consist of 1 to ∼100 tandem 3.3-kb units at subtelomeric 4q and 10q. D4Z4 displayed hypomethylation in some cancers and hypermethylation in others relative to normal tissues. Surprisingly, in cancers with extensive D4Z4 methylation there was a barrier to hypermethylation spreading to the beginning of this disease-associated array (facioscapulohumeral muscular dystrophy, FSHD) despite sequence conservation in repeat units throughout the array. We infer a different chromatin structure at the proximal end of the array than at interior repeats, consistent with results from chromatin DNaseI sensitivity assays indicating a boundary element near the beginning of the array. The relative chromatin DNaseI sensitivity in FSHD and control myoblasts and lymphoblasts was as follows: a non-genic D4Z4-adjacent sequence (p13E-11, array-proximal)> untranscribed gene standards > D4Z4 arrays> constitutive heterochromatin (satellite 2; P < 10−4 for all comparisons). Cancers displaying D4Z4 hypermethylation also exhibited a hypermethylation-resistant subregion within the 3.3-kb D4Z4 repeat units. This subregion contains runs of G that form G-quadruplexes in vitro. Unusual DNA structures might contribute to topological constraints that link short 4q D4Z4 arrays to FSHD and make long ones phenotypically neutral.

Expression profiling to predict postoperative prognosis for estrogen receptor-negative breast cancers by analysis of 25,344 genes on a cDNA microarray

Estrogen receptor (ER) status is an essential determinant of clinical and biological behavior of human breast cancers. While ER‐positive breast cancers respond well to adjuvant hormone therapy, ER‐negative tumors are generally resistant. To date, no attempts have succeeded in finding molecular markers for classifying ER‐negative breast cancers with respect to postoperative prognosis. To identify a set of prognostic markers for this type of cancer, we used a cDNA microarray consisting of 25,344 human genes to investigate expression profiles of ten primary breast cancers from patients who had died of breast cancer within 5 years after surgery (5y‐D) and 10 from patients who had survived disease‐free for more than 5 years (5y‐S). Sets of genes characterizing each group were identified by Mann‐Whitney and random‐permutation tests. We documented 71 genes with higher expression in the 5y‐D group than in the 5y‐S group, and 15 with higher expression in the 5y‐S group than in the 5y‐D group. Semi‐quantitative RT‐PCR experiments were carried out to confirm the results of the microarray analysis. We established a scoring system for predicting postoperative prognosis of ER‐negative breast cancers on the basis of aberrant gene expression. The list of genes reported here provides valuable information with regard to progression of breast cancer and is a source of possible target molecules for development of novel drugs to treat patients with ER‐negative breast cancers.

A DNA repeat, NBL2, is hypermethylated in some cancers but hypomethylated in others.

Hypermethylation at certain CpG-rich promoters and hypomethylation at repeated DNA sequences are very frequently found in cancers. We provide the first report that a DNA sequence (NBL2) can be either extensively hypermethylated or hypomethylated in cancer. Previously, it was shown that NBL2, a complex tandem DNA repeat in the acrocentric chromosomes, is hypomethylated at NotI sites in >70% of neuroblastomas and hepatocellular carcinomas and in cells from ICF syndrome (DNMT3B-deficiency) patients. Unexpectedly, by Southern blot analysis of 18 ovarian carcinomas, 51 Wilms tumors, and various somatic control tissues, we found that >70% of the cancers exhibited large increases in methylation at HhaI sites in NBL2 compared with all the controls. In contrast, 17% of the carcinomas showed major decreases in methylation at HhaI and NotI sites. The intermediate levels of methylation at HhaI sites in somatic controls enabled this discovery of cancer-linked hypermethylation and hypomethylation in NBL2. In a comparison of ovarian epithelial carcinomas, low malignant potential tumors, and cystadenomas, NBL2 hypermethylation at HhaI sites was significantly related to the degree of malignancy, and hypomethylation was seen only in the carcinomas. By RT-PCR, we found NBL2 transcripts at low levels in a few cancers and undetectable in various normal tissues. In the tumors there was no association of NBL2 hypomethylation and transcription, but this may reflect NBL2’s lack of identifiable promoter elements and our evidence for run-through transcription from adjacent sequences into NBL2. The propensity of NBL2 sequences to become either hypermethylated or hypomethylated in cancer suggests that these opposite epigenetic changes share an early step during carcinogenesis and that cancer-linked hypermethylation might be spontaneously reversible.

DNA methylation and differentiation: HOX genes in muscle cells

BackgroundTight regulation of homeobox genes is essential for vertebrate development. In a study of genome-wide differential methylation, we recently found that homeobox genes, including those in the HOX gene clusters, were highly overrepresented among the genes with hypermethylation in the skeletal muscle lineage. Methylation was analyzed by reduced representation bisulfite sequencing (RRBS) of postnatal myoblasts, myotubes and adult skeletal muscle tissue and 30 types of non-muscle-cell cultures or tissues.ResultsIn this study, we found that myogenic hypermethylation was present in specific subregions of all four HOX gene clusters and was associated with various chromatin epigenetic features. Although the 3′ half of the HOXD cluster was silenced and enriched in polycomb repression-associated H3 lysine 27 trimethylation in most examined cell types, including myoblasts and myotubes, myogenic samples were unusual in also displaying much DNA methylation in this region. In contrast, both HOXA and HOXC clusters displayed myogenic hypermethylation bordering a central region containing many genes preferentially expressed in myogenic progenitor cells and consisting largely of chromatin with modifications typical of promoters and enhancers in these cells. A particularly interesting example of myogenic hypermethylation was HOTAIR, a HOXC noncoding RNA gene, which can silence HOXD genes in trans via recruitment of polycomb proteins. In myogenic progenitor cells, the preferential expression of HOTAIR was associated with hypermethylation immediately downstream of the gene. Other HOX gene regions also displayed myogenic DNA hypermethylation despite being moderately expressed in myogenic cells. Analysis of representative myogenic hypermethylated sites for 5-hydroxymethylcytosine revealed little or none of this base, except for an intragenic site in HOXB5 which was specifically enriched in this base in skeletal muscle tissue, whereas myoblasts had predominantly 5-methylcytosine at the same CpG site.ConclusionsOur results suggest that myogenic hypermethylation of HOX genes helps fine-tune HOX sense and antisense gene expression through effects on 5′ promoters, intragenic and intergenic enhancers and internal promoters. Myogenic hypermethylation might also affect the relative abundance of different RNA isoforms, facilitate transcription termination, help stop the spread of activation-associated chromatin domains and stabilize repressive chromatin structures.

Synergistic inhibition of thyroid cancer by suppressing MAPK/PI3K/AKT pathways.

BACKGROUND Although a wide spectrum of inhibitors of the MEK/ERK and PI3K/AKT pathways have been discovered and entered clinical trials, the effects of their individual use in thyroid cancer were often disappointing. We hypothesized that dual targeting of these two pathways would be a safe and effective strategy against aggressive thyroid cancers. METHODS We examined the antiproliferative effects of the MEK/ERK inhibitor AZD6244 and the PI3K/AKT inhibitor GDC0941, individually or in combination, on thyroid cancer cells harboring both the BRAF(V600E) and PIK3CA mutations. The effects of drug exposure on both total and phosphorylated (p-) forms of AKT and ERK were monitored by Western blotting analysis. Effects of these inhibitors on cell-cycle progression and apoptosis were measured by flow cytometry and DNA-fragmentation analyses, respectively. RESULTS We observed significant toxicities to viability of cells with low concentrations of AZD6244 or GDC0941, which were synergistic when the two inhibitors were used in combination (P < 0.01). AZD6244 abrogated p-ERK and GDC0941 abrogated p-AKT levels, confirming their expected target effects. Unexpectedly, monotherapy with AZD6244 resulted in activation of the PI3K signaling pathway in some cancer cell lines and co-exposure to AZD6244 and GDC0941 was necessary to suppress both pathways. Flow cytometry showed G1 arrest. DNA fragmentation analysis showed an increased apoptosis of cells dually treated with the two inhibitors. CONCLUSION Concomitant suppression of MEK/ERK and PI3K/AKT pathways by AZD6244 and GDC0941 abrogates compensatory mechanisms of tumor survival and causes synergistic cytotoxicity in thyroid cancer cells.

DNA Replication Timing Is Maintained Genome-Wide in Primary Human Myoblasts Independent of D4Z4 Contraction in FSH Muscular Dystrophy

Facioscapulohumeral muscular dystrophy (FSHD) is linked to contraction of an array of tandem 3.3-kb repeats (D4Z4) at 4q35.2 from 11-100 copies to 1-10 copies. The extent to which D4Z4 contraction at 4q35.2 affects overall 4q35.2 chromatin organization remains unclear. Because DNA replication timing is highly predictive of long-range chromatin interactions, we generated genome-wide replication-timing profiles for FSHD and control myogenic precursor cells. We compared non-immortalized myoblasts from four FSHD patients and three control individuals to each other and to a variety of other human cell types. This study also represents the first genome-wide comparison of replication timing profiles in non-immortalized human cell cultures. Myoblasts from both control and FSHD individuals all shared a myoblast-specific replication profile. In contrast, male and female individuals were readily distinguished by monoallelic differences in replication timing at DXZ4 and other regions across the X chromosome affected by X inactivation. We conclude that replication timing is a robust cell-type specific feature that is unaffected by FSHD-related D4Z4 contraction.

Hybridization analysis of D4Z4 repeat arrays linked to FSHD.

Facioscapulohumeral muscular dystrophy (FSHD) is an autosomal dominant disease involving shortening of D4Z4, an array of tandem 3.3-kb repeat units on chromosome 4. These arrays are in subtelomeric regions of 4q and 10q and have 1–100 units. FSHD is associated with an array of 1–10 units at 4q35. Unambiguous clinical diagnosis of FSHD depends on determining the array length at 4q35, usually with the array-adjacent p13E-11 probe after pulsed-field or linear gel electrophoresis. Complicating factors for molecular diagnosis of FSHD are the phenotypically neutral 10q D4Z4 arrays, cross-hybridizing sequences elsewhere in the genome, deletions including the genomic p13E-11 sequence and part of D4Z4, translocations between 4q and 10q D4Z4 arrays, and the extremely high G + C content of D4Z4 arrays (73%). In this study, we optimized conditions for molecular diagnosis of FSHD with a 1-kb D4Z4 subfragment probe after hybridization with p13E-11. We demonstrate that these hybridization conditions allow the identification of FSHD alleles with deletions of the genomic p13E-11 sequence and aid in determination of the nonpathogenic D4Z4 arrays at 10q. Furthermore, we show that the D4Z4-like sequences present elsewhere in the genome are not tandemly arranged, like those at 4q35 and 10q26.

DNaseI hypersensitivity at gene-poor, FSH dystrophy-linked 4q35.2

A subtelomeric region, 4q35.2, is implicated in facioscapulohumeral muscular dystrophy (FSHD), a dominant disease thought to involve local pathogenic changes in chromatin. FSHD patients have too few copies of a tandem 3.3-kb repeat (D4Z4) at 4q35.2. No phenotype is associated with having few copies of an almost identical repeat at 10q26.3. Standard expression analyses have not given definitive answers as to the genes involved. To investigate the pathogenic effects of short D4Z4 arrays on gene expression in the very gene-poor 4q35.2 and to find chromatin landmarks there for transcription control, unannotated genes and chromatin structure, we mapped DNaseI-hypersensitive (DH) sites in FSHD and control myoblasts. Using custom tiling arrays (DNase-chip), we found unexpectedly many DH sites in the two large gene deserts in this 4-Mb region. One site was seen preferentially in FSHD myoblasts. Several others were mapped >0.7 Mb from genes known to be active in the muscle lineage and were also observed in cultured fibroblasts, but not in lymphoid, myeloid or hepatic cells. Their selective occurrence in cells derived from mesoderm suggests functionality. Our findings indicate that the gene desert regions of 4q35.2 may have functional significance, possibly also to FSHD, despite their paucity of known genes.

Simultaneous suppression of the MAP kinase and NF-κB pathways provides a robust therapeutic potential for thyroid cancer

The MAP kinase and NF-κB signaling pathways play an important role in thyroid cancer tumorigenesis. We aimed to examine the therapeutic potential of dually targeting the two pathways using AZD6244 and Bortezomib in combination. We evaluated their effects on cell proliferation, cell-cycle progression, apoptosis, cell migration assay, and the activation of the MAPK pathway in vitro and the in vivo using tumor size and immunohistochemical changes of Ki67 and ppRB. We found inhibition of cell growth rate by 10%, 20%, and 56% (p <0.05), migration to 55%, 61%, and 29% (p <0.05), and induction of apoptosis to 10%, 15%, and 38% (p <0.05) with AZD6244, Bortezomib, or combination, respectively. Induction of cell cycle arrest occurred only with drug combination. Dual drug treatment in the xenograft model caused a 94% reduction in tumor size (p <0.05) versus 15% with AZD6244 and 34% with Bortezomib (p < 0.05) and also reduced proliferative marker Ki67, and increased pRb dephosphorylation. Our results demonstrate a robust therapeutic potential of combining AZD6244 and Bortezomib as an effective strategy to overcome drug resistance encountered in monotherapy in the treatment of thyroid cancer, strongly supporting clinical trials to further test this strategy.