Accelerating PROTACs Discovery Through a Direct-to-Biology Platform Enabled by Modular Photoclick Chemistry.

Proteolysis targeting chimeras (PROTACs) have emerged as a promising strategy for drug discovery and exploring protein functions, offering a revolutionary therapeutic modality. Currently, the predominant approach to PROTACs discovery mainly relies on an empirical design-synthesis-evaluation process involving numerous cycles of labor-intensive synthesis-purification and bioassay data collection. Therefore, the development of innovative methods to expedite PROTAC synthesis and exploration of chemical space remains highly desired. Here, a direct-to-biology strategy is reported to streamline the synthesis of PROTAC libraries on plates, enabling the seamless transfer of reaction products to cell-based bioassays without the need for additional purification. By integrating amide coupling and light-induced primary amines and o-nitrobenzyl alcohols cyclization (PANAC) photoclick chemistry into a plate-based synthetic process, this strategy produces PROTAC libraries with high efficiency and structural diversity. Moreover, by employing this platform for PROTACs screening, we smoothly found potent PROTACs effectively inhibit triple-negative breast cancer (TNBC) cell growth and induce rapid, selective targeted degradation of cyclin-dependent kinase 9 (CDK9). The study introduces a versatile platform for assembling PROTACs on plates, followed by direct biological evaluation. This approach provides a promising opportunity for high-throughput synthesis of PROTAC libraries, thereby enhancing the efficiency of exploring chemical space and accelerating the discovery of PROTACs.

Obesity is associated with depot-specific alterations in adipocyte DNA methylation and gene expression.

The present study aimed to identify genes exhibiting concomitant obesity-dependent changes in DNA methylation and gene expression in adipose tissues in the mouse using diet-induced obese (DIO) C57BL/6J and genetically obese ob/ob mice as models. Mature adipocytes were isolated from epididymal and inguinal adipose tissues of ob/ob and DIO C57BL/6J mice. DNA methylation was analyzed by MeDIP-sequencing and gene expression by microarray analysis. The majority of differentially methylated regions (DMRs) were hypomethylated in obese mice. Global methylation of long interspersed elements indicated that hypomethylation did not reflect methyl donor deficiency. In both DIO and ob/ob mice, we observed more obesity-associated methylation changes in epididymal than in inguinal adipocytes. Assignment of DMRs to promoter, exon, intron and intergenic regions demonstrated that DIO-induced changes in DNA methylation in C57BL/6J mice occurred primarily in exons, whereas inguinal adipocytes of ob/ob mice exhibited a higher enrichment of DMRs in promoter regions than in other regions of the genome, suggesting an influence of leptin on DNA methylation in inguinal adipocytes. We observed altered methylation and expression of 9 genes in epididymal adipocytes, including the known obesity-associated genes, Ehd2 and Kctd15, and a novel candidate gene, Irf8, possibly involved in immune type 1/type2 balance. The use of 2 obesity models enabled us to dissociate changes associated with high fat feeding from those associated with obesity per se. This information will be of value in future studies on the mechanisms governing the development of obesity and changes in adipocyte function associated with obesity.

C14orf93 (RTFC) is identified as a novel susceptibility gene for familial nonmedullary thyroid cancer.

The genetic causes for familial nonmedullary thyroid cancer (FNMTC) remain largely unknown. Through genetic linkage analysis and exome sequencing, C14orf93 (RTFC), PYGL, and BMP4 were identified as susceptibility gene candidates in a FNMTC family. By examining the expression and the oncogenic functions of these candidate genes, PYGL and BMP4 were excluded. We further characterized the functions of the uncharacterized gene RTFC in thyroid cancer. RTFC promotes thyroid cancer cell survival under starving conditions, and thyroid cancer cell migration. The R115Q, V205M and G209D RTFC mutants enhance the colony forming capacity of thyroid cancer cells, and are able to transform normal thyroid cells. In summary, our data suggest the roles of RTFC in thyroid carcinogenesis.

Distinct epigenomes in CD4+ T cells of newborns, middle-ages and centenarians.

Age-related variations in genes and microRNAs expression and DNA methylation have been reported respectively; however, their interactions during aging are unclear. We therefore investigated alterations in the transcriptomes, miRNAomes and DNA methylomes in the same CD4+T cells from newborn (NB), middle-aged (MA) and long-lived (LL) individuals to elucidate the molecular changes and their interactions. A total 659 genes showed significantly expression changes across NB, MA and LL individuals, in which we identified four age-related co-expression modules with three hub networks of co-expressed genes and non-coding RNAs. Moreover, we identified 9835 differentially methylated regions (DMRs) including 7015 hypermethylated and 2820 hypomethylated DMRs in the NB compared with the MA, and 12,362 DMRs including 4809 hypermethylated and 7553 hypomethylated DMRs in the MA compared with the LL. The integrated analysis revealed a potential relationship between genes transcription and DNA methylation for many age- or immune-related genes, suggesting that DNA methylation-dependent transcription regulation is involved in development and functions of T cells during aging. Our results reveals age-related transcription and methylation changes and their interactions in human T cells from the cradle to the grave. Longitudinal work is required to establish the relationship between identified age-associated genes/DNA methylation and T cells aging phenotypes.

Genomic heterogeneity of multiple synchronous lung cancer.

Multiple synchronous lung cancers (MSLCs) present a clinical dilemma as to whether individual tumours represent intrapulmonary metastases or independent tumours. In this study we analyse genomic profiles of 15 lung adenocarcinomas and one regional lymph node metastasis from 6 patients with MSLC. All 15 lung tumours demonstrate distinct genomic profiles, suggesting all are independent primary tumours, which are consistent with comprehensive histopathological assessment in 5 of the 6 patients. Lung tumours of the same individuals are no more similar to each other than are lung adenocarcinomas of different patients from TCGA cohort matched for tumour size and smoking status. Several known cancer-associated genes have different mutations in different tumours from the same patients. These findings suggest that in the context of identical constitutional genetic background and environmental exposure, different lung cancers in the same individual may have distinct genomic profiles and can be driven by distinct molecular events.

Integrative genomic and functional analysis of human oral squamous cell carcinoma cell lines reveals synergistic effects of FAT1 and CASP8 inactivation.

Oral squamous cell carcinoma (OSCC) is genetically highly heterogeneous, which contributes to the challenges of treatment. To create an in vitro model that accurately reflects this heterogeneity, we generated a panel of HPV-negative OSCC cell lines. By whole exome sequencing of the lines and matched patient blood samples, we demonstrate that the mutational spectrum of the lines is representative of primary OSCC in The Cancer Genome Atlas. We show that loss of function mutations in FAT1 (an atypical cadherin) and CASP8 (Caspase 8) frequently occur in the same tumour. OSCC cells with inactivating FAT1 mutations exhibited reduced intercellular adhesion. Knockdown of FAT1 and CASP8 individually or in combination in OSCC cells led to increased cell migration and clonal growth, resistance to Staurosporine-induced apoptosis and, in some cases, increased terminal differentiation. The OSCC lines thus represent a valuable resource for elucidating the impact of different mutations on tumour behaviour.

Vibrio anguillarum Is Genetically and Phenotypically Unaffected by Long-Term Continuous Exposure to the Antibacterial Compound Tropodithietic Acid.

UNLABELLED: Minimizing the use of antibiotics in the food production chain is essential for limiting the development and spread of antibiotic-resistant bacteria. One alternative intervention strategy is the use of probiotic bacteria, and bacteria of the marine Roseobacter clade are capable of antagonizing fish-pathogenic vibrios in fish larvae and live feed cultures for fish larvae. The antibacterial compound tropodithietic acid (TDA), an antiporter that disrupts the proton motive force, is key in the antibacterial activity of several roseobacters. Introducing probiotics on a larger scale requires understanding of any potential side effects of long-term exposure of the pathogen to the probionts or any compounds they produce. Here we exposed the fish pathogen Vibrio anguillarum to TDA for several hundred generations in an adaptive evolution experiment. No tolerance or resistance arose during the 90 days of exposure, and whole-genome sequencing of TDA-exposed lineages and clones revealed few mutational changes, compared to lineages grown without TDA. Amino acid-changing mutations were found in two to six different genes per clone; however, no mutations appeared unique to the TDA-exposed lineages or clones. None of the virulence genes of V. anguillarum was affected, and infectivity assays using fish cell lines indicated that the TDA-exposed lineages and clones were less invasive than the wild-type strain. Thus, long-term TDA exposure does not appear to result in TDA resistance and the physiology of V. anguillarum appears unaffected, supporting the application of TDA-producing roseobacters as probiotics in aquaculture. IMPORTANCE: It is important to limit the use of antibiotics in our food production, to reduce the risk of bacteria developing antibiotic resistance. We showed previously that marine bacteria of the Roseobacter clade can prevent or reduce bacterial diseases in fish larvae, acting as probiotics. Roseobacters produce the antimicrobial compound tropodithietic acid (TDA), and we were concerned regarding whether long-term exposure to this compound could induce resistance or affect the disease-causing ability of the fish pathogen. Therefore, we exposed the fish pathogen Vibrio anguillarum to increasing TDA concentrations over 3 months. We did not see the development of any resistance to TDA, and subsequent infection assays revealed that none of the TDA-exposed clones had increased virulence toward fish cells. Hence, this study supports the use of roseobacters as a non-risk-based disease control measure in aquaculture.

Whole exome sequencing of rare aggressive breast cancer histologies.

Little is known about mutational landscape of rare breast cancer (BC) subtypes. The aim of the study was to apply next generation sequencing to three different subtypes of rare BCs in order to identify new genes related to cancer progression. We performed whole exome and targeted sequencing of 29 micropapillary, 23 metaplastic, and 27 pleomorphic lobular BCs. Micropapillary BCs exhibit a profile comparable to common BCs: PIK3CA, TP53, GATA3, and MAP2K4 were the most frequently mutated genes. Metaplastic BCs presented a high frequency of TP53 (78 %) and PIK3CA (48 %) mutations and were recurrently mutated on KDM6A (13 %), a gene involved in histone demethylation. Pleomorphic lobular carcinoma exhibited high mutation rate of PIK3CA (30 %), TP53 (22 %), and CDH1 (41 %) and also presented mutations in PYGM, a gene involved in glycogen metabolism, in 8 out of 27 samples (30 %). Further analyses of publicly available datasets showed that PYGM is dramatically underexpressed in common cancers as compared to normal tissues and that low expression in tumors is correlated with poor relapse-free survival. Immunohistochemical staining on formalin-fixed paraffin-embedded tissues available in our cohort of patients confirmed higher PYGM expression in normal breast tissue compared to equivalent tumoral zone. Next generation sequencing methods applied on rare cancer subtypes can serve as a useful tool in order to uncover new potential therapeutic targets. Sequencing of pleomorphic lobular carcinoma identified a high rate of alterations in PYGM. These findings emphasize the role of glycogen metabolism in cancer progression.

Analysis of a four generation family reveals the widespread sequence-dependent maintenance of allelic DNA methylation in somatic and germ cells.

Differential methylation of the homologous chromosomes, a well-known mechanism leading to genomic imprinting and X-chromosome inactivation, is widely reported at the non-imprinted regions on autosomes. To evaluate the transgenerational DNA methylation patterns in human, we analyzed the DNA methylomes of somatic and germ cells in a four-generation family. We found that allelic asymmetry of DNA methylation was pervasive at the non-imprinted loci and was likely regulated by cis-acting genetic variants. We also observed that the allelic methylation patterns for the vast majority of the cis-regulated loci were shared between the somatic and germ cells from the same individual. These results demonstrated the interaction between genetic and epigenetic variations and suggested the possibility of widespread sequence-dependent transmission of DNA methylation during spermatogenesis.

Whole Exome Sequencing Identifies Frequent Somatic Mutations in Cell-Cell Adhesion Genes in Chinese Patients with Lung Squamous Cell Carcinoma.

Lung squamous cell carcinoma (SQCC) accounts for about 30% of all lung cancer cases. Understanding of mutational landscape for this subtype of lung cancer in Chinese patients is currently limited. We performed whole exome sequencing in samples from 100 patients with lung SQCCs to search for somatic mutations and the subsequent target capture sequencing in another 98 samples for validation. We identified 20 significantly mutated genes, including TP53, CDH10, NFE2L2 and PTEN. Pathways with frequently mutated genes included those of cell-cell adhesion/Wnt/Hippo in 76%, oxidative stress response in 21%, and phosphatidylinositol-3-OH kinase in 36% of the tested tumor samples. Mutations of Chromatin regulatory factor genes were identified at a lower frequency. In functional assays, we observed that knockdown of CDH10 promoted cell proliferation, soft-agar colony formation, cell migration and cell invasion, and overexpression of CDH10 inhibited cell proliferation. This mutational landscape of lung SQCC in Chinese patients improves our current understanding of lung carcinogenesis, early diagnosis and personalized therapy.

Immature truncated O-glycophenotype of cancer directly induces oncogenic features.

Aberrant expression of immature truncated O-glycans is a characteristic feature observed on virtually all epithelial cancer cells, and a very high frequency is observed in early epithelial premalignant lesions that precede the development of adenocarcinomas. Expression of the truncated O-glycan structures Tn and sialyl-Tn is strongly associated with poor prognosis and overall low survival. The genetic and biosynthetic mechanisms leading to accumulation of truncated O-glycans are not fully understood and include mutation or dysregulation of glycosyltransferases involved in elongation of O-glycans, as well as relocation of glycosyltransferases controlling initiation of O-glycosylation from Golgi to endoplasmic reticulum. Truncated O-glycans have been proposed to play functional roles for cancer-cell invasiveness, but our understanding of the biological functions of aberrant glycosylation in cancer is still highly limited. Here, we used exome sequencing of most glycosyltransferases in a large series of primary and metastatic pancreatic cancers to rule out somatic mutations as a cause of expression of truncated O-glycans. Instead, we found hypermethylation of core 1 ß3-Gal-T-specific molecular chaperone, a key chaperone for O-glycan elongation, as the most prevalent cause. We next used gene editing to produce isogenic cell systems with and without homogenous truncated O-glycans that enabled, to our knowledge, the first polyomic and side-by-side evaluation of the cancer O-glycophenotype in an organotypic tissue model and in xenografts. The results strongly suggest that truncation of O-glycans directly induces oncogenic features of cell growth and invasion. The study provides support for targeting cancer-specific truncated O-glycans with immunotherapeutic measures.

Identification of genomic alterations in oesophageal squamous cell cancer.

Oesophageal cancer is one of the most aggressive cancers and is the sixth leading cause of cancer death worldwide. Approximately 70% of global oesophageal cancer cases occur in China, with oesophageal squamous cell carcinoma (ESCC) being the histopathological form in the vast majority of cases (>90%). Currently, there are limited clinical approaches for the early diagnosis and treatment of ESCC, resulting in a 10% five-year survival rate for patients. However, the full repertoire of genomic events leading to the pathogenesis of ESCC remains unclear. Here we describe a comprehensive genomic analysis of 158 ESCC cases, as part of the International Cancer Genome Consortium research project. We conducted whole-genome sequencing in 17 ESCC cases and whole-exome sequencing in 71 cases, of which 53 cases, plus an additional 70 ESCC cases not used in the whole-genome and whole-exome sequencing, were subjected to array comparative genomic hybridization analysis. We identified eight significantly mutated genes, of which six are well known tumour-associated genes (TP53, RB1, CDKN2A, PIK3CA, NOTCH1, NFE2L2), and two have not previously been described in ESCC (ADAM29 and FAM135B). Notably, FAM135B is identified as a novel cancer-implicated gene as assayed for its ability to promote malignancy of ESCC cells. Additionally, MIR548K, a microRNA encoded in the amplified 11q13.3-13.4 region, is characterized as a novel oncogene, and functional assays demonstrate that MIR548K enhances malignant phenotypes of ESCC cells. Moreover, we have found that several important histone regulator genes (MLL2 (also called KMT2D), ASH1L, MLL3 (KMT2C), SETD1B, CREBBP and EP300) are frequently altered in ESCC. Pathway assessment reveals that somatic aberrations are mainly involved in the Wnt, cell cycle and Notch pathways. Genomic analyses suggest that ESCC and head and neck squamous cell carcinoma share some common pathogenic mechanisms, and ESCC development is associated with alcohol drinking. This study has explored novel biological markers and tumorigenic pathways that would greatly improve therapeutic strategies for ESCC.

Integrative analyses of gene expression and DNA methylation profiles in breast cancer cell line models of tamoxifen-resistance indicate a potential role of cells with stem-like properties.

INTRODUCTION: Development of resistance to tamoxifen is an important clinical issue in the treatment of breast cancer. Tamoxifen resistance may be the result of acquisition of epigenetic regulation within breast cancer cells, such as DNA methylation, resulting in changed mRNA expression of genes pivotal for estrogen-dependent growth. Alternatively, tamoxifen resistance may be due to selection of pre-existing resistant cells, or a combination of the two mechanisms. METHODS: To evaluate the contribution of these possible tamoxifen resistance mechanisms, we applied modified DNA methylation-specific digital karyotyping (MMSDK) and digital gene expression (DGE) in combination with massive parallel sequencing to analyze a well-established tamoxifen-resistant cell line model (TAM(R)), consisting of 4 resistant and one parental cell line. Another tamoxifen-resistant cell line model system (LCC1/LCC2) was used to validate the DNA methylation and gene expression results. RESULTS: Significant differences were observed in global gene expression and DNA methylation profiles between the parental tamoxifen-sensitive cell line and the 4 tamoxifen-resistant TAM(R) sublines. The 4 TAM(R) cell lines exhibited higher methylation levels as well as an inverse relationship between gene expression and DNA methylation in the promoter regions. A panel of genes, including NRIP1, HECA and FIS1, exhibited lower gene expression in resistant vs. parental cells and concurrent increased promoter CGI methylation in resistant vs. parental cell lines. A major part of the methylation, gene expression, and pathway alterations observed in the TAM(R) model were also present in the LCC1/LCC2 cell line model. More importantly, high expression of SOX2 and alterations of other SOX and E2F gene family members, as well as RB-related pocket protein genes in TAMR highlighted stem cell-associated pathways as being central in the resistant cells and imply that cancer-initiating cells/cancer stem-like cells may be involved in tamoxifen resistance in this model. CONCLUSION: Our data highlight the likelihood that resistant cells emerge from cancer-initiating cells/cancer stem-like cells and imply that these cells may gain further advantage in growth via epigenetic mechanisms. Illuminating the expression and DNA methylation features of putative cancer-initiating cells/cancer stem cells may suggest novel strategies to overcome tamoxifen resistance.

Whole-genome DNA methylation in skin lesions from patients with psoriasis vulgaris.

Psoriasis, a chronic inflammatory skin disorder, is characterized by aberrant keratinocyte proliferation and differentiation in the epidermis. Although the pathogenesis of psoriasis is still incompletely understood, both genetic susceptibilities and environmental triggers are known to act as key players in its development. Several studies have suggested that DNA methylation is involved in the pathogenesis of psoriasis. However, the precise mechanisms underlying the regulation and maintenance of the methylome as well as their relationship with this disease remain poorly characterized. Herein, we used methylated DNA immunoprecipitation sequencing (MeDIP-Seq) to characterize whole-genome DNA methylation patterns in involved and uninvolved skin lesions from patients with psoriasis. The results of our MeDIP-Seq analyses identified differentially methylated regions (DMRs) covering almost the entire genome with sufficient depth and high resolution, showing that the number of hypermethylated DMRs was considerably higher than that of hypomethylated DMRs in involved psoriatic skin samples. Moreover, gene ontology analysis of MeDIP-Seq data showed that the aberrantly methylated genes belonged to several different ontological domains, such as the immune system, cell cycle and apoptosis. The results of the bisulfite-sequencing experiments for the genes PDCD5 and TIMP2 confirmed the methylation status identified by MeDIP-Seq, and the mRNA expression levels of these two genes were consistent with their DNA methylation profiles. To our knowledge, the present study constitutes the first report on MeDIP-Seq in psoriasis. The identification of whole-genome DNA methylation patterns associated with psoriasis provides new insight into the pathogenesis of this complex disease and represents a promising avenue through which to investigate novel therapeutic approaches.

An atlas of DNA methylomes in porcine adipose and muscle tissues.

It is evident that epigenetic factors, especially DNA methylation, have essential roles in obesity development. Here, using pig as a model, we investigate the systematic association between DNA methylation and obesity. We sample eight variant adipose and two distinct skeletal muscle tissues from three pig breeds living within comparable environments but displaying distinct fat level. We generate 1,381 Gb of sequence data from 180 methylated DNA immunoprecipitation libraries, and provide a genome-wide DNA methylation map as well as a gene expression map for adipose and muscle studies. The analysis shows global similarity and difference among breeds, sexes and anatomic locations, and identifies the differentially methylated regions. The differentially methylated regions in promoters are highly associated with obesity development via expression repression of both known obesity-related genes and novel genes. This comprehensive map provides a solid basis for exploring epigenetic mechanisms of adipose deposition and muscle growth.

A systematic analysis on DNA methylation and the expression of both mRNA and microRNA in bladder cancer.

BACKGROUND: DNA methylation aberration and microRNA (miRNA) deregulation have been observed in many types of cancers. A systematic study of methylome and transcriptome in bladder urothelial carcinoma has never been reported. METHODOLOGY/PRINCIPAL FINDINGS: The DNA methylation was profiled by modified methylation-specific digital karyotyping (MMSDK) and the expression of mRNAs and miRNAs was analyzed by digital gene expression (DGE) sequencing in tumors and matched normal adjacent tissues obtained from 9 bladder urothelial carcinoma patients. We found that a set of significantly enriched pathways disrupted in bladder urothelial carcinoma primarily related to "neurogenesis" and "cell differentiation" by integrated analysis of -omics data. Furthermore, we identified an intriguing collection of cancer-related genes that were deregulated at the levels of DNA methylation and mRNA expression, and we validated several of these genes (HIC1, SLIT2, RASAL1, and KRT17) by Bisulfite Sequencing PCR and Reverse Transcription qPCR in a panel of 33 bladder cancer samples. CONCLUSIONS/SIGNIFICANCE: We characterized the profiles between methylome and transcriptome in bladder urothelial carcinoma, identified a set of significantly enriched key pathways, and screened four aberrantly methylated and expressed genes. Conclusively, our findings shed light on a new avenue for basic bladder cancer research.

DNA methylation in liver tumorigenesis in fish from the environment.

The link between environment, alteration in DNA methylation and cancer has been well established in humans; yet, it is under-studied in unsequenced non-model organisms. The occurrence of liver tumors in the flatfish dab collected at certain UK sampling sites exceeds 20%, yet the causative agents and the molecular mechanisms of tumor formation are not known, especially regarding the balance between epigenetic and genetic factors. Methylated DNA Immunoprecipitation (MeDIP) combined with de novo high-throughput DNA sequencing were used to investigate DNA methylation changes in dab hepatocellular adenoma tumors for the first time in an unsequenced species. Novel custom-made dab gene expression arrays were designed and used to determine the relationship between DNA methylation and gene expression. In addition, the confirmatory techniques of bisulfite sequencing PCR (BSP) and RT-PCR were applied. Genes involved in pathways related to cancer, including apoptosis, wnt/ß-catenin signaling and genomic and non-genomic estrogen responses, were altered both in methylation and transcription. Global methylation was statistically significantly 1.8-fold reduced in hepatocellular adenoma and non-cancerous surrounding tissues compared with liver from non-cancer bearing dab. Based on the identified changes and chemical exposure data, our study supports the epigenetic model of cancer. We hypothesize that chronic exposure to a mixture of environmental contaminants contributes to a global hypomethylation followed by further epigenetic and genomic changes. The findings suggest a link between environment, epigenetics and cancer in fish tumors in the wild and show the utility of this methodology for studies in non-model organisms.

Genome-wide mapping of DNA methylation in chicken.

Cytosine DNA methylation is an important epigenetic modification termed as the fifth base that functions in diverse processes. Till now, the genome-wide DNA methylation maps of many organisms has been reported, such as human, Arabidopsis, rice and silkworm, but the methylation pattern of bird remains rarely studied. Here we show the genome-wide DNA methylation map of bird, using the chicken as a model organism and an immunocapturing approach followed by high-throughput sequencing. In both of the red jungle fowl and the avian broiler, DNA methylation was described separately for the liver and muscle tissue. Generally, chicken displays analogous methylation pattern with that of animals and plants. DNA methylation is enriched in the gene body regions and the repetitive sequences, and depleted in the transcription start site (TSS) and the transcription termination site (TTS). Most of the CpG islands in the chicken genome are kept in unmethylated state. Promoter methylation is negatively correlated with the gene expression level, indicating its suppressive role in regulating gene transcription. This work contributes to our understanding of epigenetics in birds.

An improved method for genome wide DNA methylation profiling correlated to transcription and genomic instability in two breast cancer cell lines.

BACKGROUND: DNA methylation is a widely studied epigenetic mechanism known to correlate with gene repression and genomic stability. Development of sensitive methods for global detection of DNA methylation events is of particular importance. RESULTS: We here describe a technique, called modified methylation-specific digital karyotyping (MMSDK) based on methylation-specific digital karyotyping (MSDK) with a novel sequencing approach. Briefly, after a tandem digestion of genomic DNA with a methylation-sensitive mapping enzyme and a fragmenting enzyme, short sequence tags are obtained. These tags are amplified, followed by direct, massively parallel sequencing (Solexa 1G Genome Analyzer). This method allows high-throughput and low-cost genome-wide DNA methylation mapping. We applied this method to investigate global DNA methylation profiles for widely used breast cancer cell lines, MCF-7 and MDA-MB-231, which are representatives for luminal-like and mesenchymal-like cancer types, respectively. By comparison, a highly similar overall DNA methylation pattern was revealed for the two cell lines. However a cohort of individual genomic loci with significantly different DNA methylation status between two cell lines was identified. Furthermore, we revealed a genome-wide significant correlation between gene expression and the methylation status of gene promoters with CpG islands (CGIs) in the two cancer cell lines, and a correlation of gene expression and the methylation status of promoters without CGIs in MCF-7 cells. CONCLUSION: The MMSDK method will be a valuable tool to increase the current knowledge of genome wide DNA methylation profiles.