Automated quantification of DNA demethylation effects in cells via 3D mapping of nuclear signatures and population homogeneity assessment

DNA methylation topology differentiates between normal and malignant in cell models, resected human tissues, and exfoliated sputum cells of lung epithelium

Background

Global DNA hypomethylation is a prominent feature of cancer cells including lung cancer, that has not been widely explored towards cancer diagnosis. In this study we assess the comparative distribution of global DNA methylation in normal cells versus cancer cells in various specimen models.

Methods

We used in situ immunofluorescence labeling of overall 5-methylcytosine (5mC) and covisualization of global DNA (gDNA) by 4’,6-diamidino-2-phenylindole (DAPI), confocal microscopy and 3D image analysis to derive 5mC/DAPI colocalization patterns in human cell lines (BEAS-2B, A549, H157) and upper respiratory epithelial cells derived from various sources (i.e., sputum from healthy and cancer patients, and resected tissues from normal parenchyma and lung tumors).

Results

By introducing 5mC/DAPI colocalization index as a metric we could distinguish between normal epithelial cells and aberrantly hypomethylated cancer cells. Cultured lung cancer cells (H157 and A549) had significantly lower indices compared to normal cells (BEAS-2B). Furthermore, we were able to identify such extensively hypomethylated low-index cells in tumor tissues and the matching sputum from cancer patients. In contrast, the indices of cells derived from sputum of healthy individuals had more similarity to epithelial cells of normal parenchyma and the phenotypically normal BEAS-2B cells.

Conclusions

The results suggest that 5mC topology using high-resolution image cytometry shows potential for identifying hypomethylated cancerous cells in human tissues and amongst normal cells in matching sputum, which may render a valuable surrogate for biopsied tissues. This promising feature deserves further validation in more comprehensive studies.

Biomedical Applications of Translational Optical Imaging: From Molecules to Humans

Light is a powerful investigational tool in biomedicine, at all levels of structural organization. Its multitude of features (intensity, wavelength, polarization, interference, coherence, timing, non-linear absorption, and even interactions with itself) able to create contrast, and thus images that detail the makeup and functioning of the living state can and should be combined for maximum effect, especially if one seeks simultaneously high spatiotemporal resolution and discrimination ability within a living organism. The resulting high relevance should be directed towards a better understanding, detection of abnormalities, and ultimately cogent, precise, and effective intervention. The new optical methods and their combinations needed to address modern surgery in the operating room of the future, and major diseases such as cancer and neurodegeneration are reviewed here, with emphasis on our own work and highlighting selected applications focusing on quantitation, early detection, treatment assessment, and clinical relevance, and more generally matching the quality of the optical detection approach to the complexity of the disease. This should provide guidance for future advanced theranostics, emphasizing a tighter coupling-spatially and temporally-between detection, diagnosis, and treatment, in the hope that technologic sophistication such as that of a Mars rover can be translationally deployed in the clinic, for saving and improving lives.

Covisualization of Global DNA Methylation/Hydroxymethylation and Protein Biomarkers for Ultrahigh-Definition Epigenetic Phenotyping of Stem Cells

DNA methylation and DNA hydroxymethylation are genomic-scale key regulatory modifications in cellular differentiation and are skewed in complex diseases. Therefore, analyzing the nuclear distribution of globally methylated and hydroxymethylated DNA in conjunction with relevant cellular components, such as protein biomarkers, may well add cell-by-cell-specific spatial and temporal information to quantitative molecular data for the discovery of signaling networks in stem cell differentiation and their exploitation in the therapeutic reprogramming of cells. Fluorescence imaging provides an optical approach that has become an essential tool in this context. The in situ fluorescent covisualization of globally methylated and hydroxymethylated DNA (5-methylcytosine = 5mC, 5-hydroxymethylcytosine = 5hmC), global DNA (gDNA), and proteins can be challenging, as the immunofluorescence detection of 5mC and 5hmC sites requires thorough denaturing of double-stranded DNA for antigen retrieval. The protocol we present overcomes this obstacle through optimization of the necessary cell processing to delineate cytosine variants and gDNA while preserving the three-dimensional (3-D) structure of the cells and in connection the immunostaining of protein biomarkers and DNA counterstaining, making it suitable for ultrahigh definition (UHD) imaging of single cells by confocal and super-resolution microscopy, 3-D visualization, and high-content cytometry.

An Inhibitor of GSK3B and HDACs Kills Pancreatic Cancer Cells and Slows Pancreatic Tumor Growth and Metastasis in Mice

BACKGROUND & AIMS

Growth, progression, and drug resistance of pancreatic ductal adenocarcinomas (PDACs) have been associated with increased levels and activity of glycogen synthase kinase 3 beta (GSK3B) and histone deacetylases (HDACs). We designed and synthesized molecules that simultaneously inhibit the activities of both enzymes. We tested the effects of one of these molecules, Metavert, in pancreatic cancer cells and mice with pancreatic tumors.

METHODS

We tested the ability of Metavert to bind GSK3B and HDACs using surface plasmon resonance. MIA PaCa-2, Bx-PC3, HPAF-II, and HPDE6 cell lines were incubated with different concentrations of Metavert, with or without paclitaxel or gemcitabine, or with other inhibitors of GSK3B and HDACs; cells were analyzed for apoptosis and migration and by immunoblotting, immunofluorescence, and real-time polymerase chain reaction. Krasþ/LSLG12D;Trp53þ/LSLR172H;Pdx-1-Cre (KPC) mice (2 months old) were given injections of Metavert (5 mg/kg, 3 times/week) or vehicle (control). B6.129J mice with tumors grown from UN-KPC961-Luc cells were given injections of Metavert or vehicle. Tumors and metastases were counted and pancreata were analyzed by immunohistochemistry. Glucose metabolism was measured using 13C-glucose tracer and mass spectroscopy and flow cytometry. Cytokine levels in blood samples were measured using multiplexing enzyme-linked immunosorbent assay.

RESULTS

Metavert significantly reduced survival of PDAC cells but not nontransformed cells; the agent reduced markers of the epithelial-to-mesenchymal transition and stem cells in PDAC cell lines. Cells incubated with Metavert in combination with irradiation and paclitaxel or gemcitabine had reduced survival compared with cells incubated with either agent alone; Metavert increased killing of drug-resistant PDAC cells by paclitaxel and gemcitabine. PDAC cells incubated with Metavert acquired normalized glucose metabolism. Administration of Metavert (alone or in combination with gemcitibine) to KPC mice or mice with syngeneic tumors significantly increased their survival times, slowed tumor growth, prevented tumor metastasis, decreased tumor infiltration by tumor-associated macrophages, and decreased blood levels of cytokines.

CONCLUSIONS

In studies of PDAC cells and 2 mouse models of PDAC, we found a dual inhibitor of GSK3B and HDACs (Metavert) to induce cancer cell apoptosis, reduce migration and expression of stem cell markers, and slow growth of tumors and metastases. Metavert had synergistic effects with gemcitabine.

Prostate cancer diagnosis using epigenetic biomarkers, 3D high-content imaging and probabilistic cell-by-cell classifiers

Background

Prostate cancer (PCa) management can benefit from novel concepts/biomarkers for reducing the current 20-30% chance of false-negative diagnosis with standard histopathology of biopsied tissue.

Method

We explored the potential of selected epigenetic markers in combination with validated histopathological markers, 3D high-content imaging, cell-by-cell analysis, and probabilistic classification in generating novel detailed maps of biomarker heterogeneity in patient tissues, and PCa diagnosis. We used consecutive biopsies/radical prostatectomies from five patients for building a database of ∼140,000 analyzed cells across all tissue compartments and for model development; and from five patients and the two well-characterized HPrEpiC primary and LNCaP cancer cell types for model validation.

Results

Principal component analysis presented highest covariability for the four biomarkers 4′,6-diamidino-2-phenylindole, 5-methylcytosine, 5-hydroxymethylcytosine, and alpha-methylacyl-CoA racemase in the epithelial tissue compartment. The panel also showed best performance in discriminating between normal and cancer-like cells in prostate tissues with a sensitivity and specificity of 85%, correctly classified 87% of HPrEpiC as healthy and 99% of LNCaP cells as cancer-like, identified a majority of aberrant cells within histopathologically benign tissues at baseline diagnosis of patients that were later diagnosed with adenocarcinoma. Using k-nearest neighbor classifier with cells from an initial patient biopsy, the biomarkers were able to predict cancer stage and grade of prostatic tissue that occurred at later prostatectomy with 79% accuracy.

Conclusion

Our approach showed favorable diagnostic values to identify the portion and pathological category of aberrant cells in a small subset of sampled tissue cells, correlating with the degree of malignancy beyond baseline.

Regulation of microtubule dynamics by DIAPH3 influences amoeboid tumor cell mechanics and sensitivity to taxanes

Taxanes are widely employed chemotherapies for patients with metastatic prostate and breast cancer. Here, we show that loss of Diaphanous-related formin-3 (DIAPH3), frequently associated with metastatic breast and prostate cancers, correlates with increased sensitivity to taxanes. DIAPH3 interacted with microtubules (MT), and its loss altered several parameters of MT dynamics as well as decreased polarized force generation, contractility, and response to substrate stiffness. Silencing of DIAPH3 increased the cytotoxic response to taxanes in prostate and breast cancer cell lines. Analysis of drug activity for tubulin-targeted agents in the NCI-60 cell line panel revealed a uniform positive correlation between reduced DIAPH3 expression and drug sensitivity. Low DIAPH3 expression correlated with improved relapse-free survival in breast cancer patients treated with chemotherapeutic regimens containing taxanes. Our results suggest that inhibition of MT stability arising from DIAPH3 downregulation enhances susceptibility to MT poisons, and that the DIAPH3 network potentially reports taxane sensitivity in human tumors.

Dynamic heterogeneity of DNA methylation and hydroxymethylation in embryonic stem cell populations captured by single-cell 3D high-content analysis

Cell-surface markers and transcription factors are being used in the assessment of stem cell fate and therapeutic safety, but display significant variability in stem cell cultures. We assessed nuclear patterns of 5-hydroxymethylcytosine (5hmC, associated with pluripotency), a second important epigenetic mark, and its combination with 5-methylcytosine (5mC, associated with differentiation), also in comparison to more established markers of pluripotency (Oct-4) and endodermal differentiation (FoxA2, Sox17) in mouse embryonic stem cells (mESC) over a 10-day differentiation course in vitro: by means of confocal and super-resolution imaging together with 3D high-content analysis, an essential tool in single-cell screening.

IN SUMMARY

1) We did not measure any significant correlation of putative markers with global 5mC or 5hmC. 2) While average Oct-4 levels stagnated on a cell-population base (0.015 lnIU/day), Sox17 and FoxA2 increased 22-fold and 3-fold faster, respectively (Sox17: 0.343 lnIU/day; FoxA2: 0.046 lnIU/day). In comparison, global DNA methylation levels increased 4-fold faster (0.068 lnIU/day), and global hydroxymethylation declined at 0.046 lnIU/day, both with a better explanation of the temporal profile. 3) This progression was concomitant with the occurrence of distinct nuclear codistribution patterns that represented a heterogeneous spectrum of states in differentiation; converging to three major coexisting 5mC/5hmC phenotypes by day 10: 5hmC(+)/5mC(-), 5hmC(+)/5mC(+), and 5hmC(-)/5mC(+) cells. 4) Using optical nanoscopy we could delineate the respective topologies of 5mC/5hmC colocalization in subregions of nuclear DNA: in the majority of 5hmC(+)/5mC(+) cells 5hmC and 5mC predominantly occupied mutually exclusive territories resembling euchromatic and heterochromatic regions, respectively. Simultaneously, in a smaller subset of cells we observed a tighter colocalization of the two cytosine variants, presumably delineating chromatin domains in remodeling. We conclude that 1) 5mC emerges as the most differential marker in our model system. 2) However, the combined enrollment of the two DNA modifications provided higher-definition screening and lead to the identification of cell subpopulations based on differential 5hmC/5mC phenotypes corresponding to different 5hmC/5mC ratios. The results encourage: a) assessing the regenerative potential of early-endodermal cells enriched for the three DNA methylation/hydroxymethylation categories, and b) exploring the universality of this type of epigenetic phenotyping across other lineage-specific differentiations.

Nuclear DNA methylation and chromatin condensation phenotypes are distinct between normally proliferating/aging, rapidly growing/immortal, and senescent cells

This study reports on probing the utility of in situ chromatin texture features such as nuclear DNA methylation and chromatin condensation patterns — visualized by fluorescent staining and evaluated by dedicated three-dimensional (3D) quantitative and high-throughput cell-by-cell image analysis — in assessing the proliferative capacity, i.e. growth behavior of cells: to provide a more dynamic picture of a cell population with potential implications in basic science, cancer diagnostics/prognostics and therapeutic drug development. Two types of primary cells and four different cancer cell lines were propagated and subjected to cell-counting, flow cytometry, confocal imaging, and 3D image analysis at various points in culture. Additionally a subset of primary and cancer cells was accelerated into senescence by oxidative stress. DNA methylation and chromatin condensation levels decreased with declining doubling times when primary cells aged in culture with the lowest levels reached at the stage of proliferative senescence. In comparison, immortal cancer cells with constant but higher doubling times mostly displayed lower and constant levels of the two in situ-derived features. However, stress-induced senescent primary and cancer cells showed similar levels of these features compared with primary cells that had reached natural growth arrest. With regards to global DNA methylation and chromatin condensation levels, aggressively growing cancer cells seem to take an intermediate level between normally proliferating and senescent cells. Thus, normal cells apparently reach cancer-cell equivalent stages of the two parameters at some point in aging, which might challenge phenotypic distinction between these two types of cells. Companion high-resolution molecular profiling could provide information on possible underlying differences that would explain benign versus malign cell growth behaviors.

Epigenetically induced changes in nuclear textural patterns and gelatinase expression in human fibrosarcoma cells

OBJECTIVE

Chromatin texture patterns of tumour cell nuclei can serve as cancer biomarkers, either to define diagnostic classifications or to obtain relevant prognostic information, in a large number of human tumours. Epigenetic mechanisms, mainly DNA methylation and histone post-translational modification, have been shown to influence chromatin packing states, and therefore nuclear texture. The aim of this study was to analyse effects of these two mechanisms on chromatin texture, and also on correlation with gelatinase expression, in human fibrosarcoma tumour cells.

MATERIALS AND METHODS

We investigated effects of DNA hypomethylating agent 5-aza-2'-deoxycytidine (5-azadC) and histone deacetylase inhibitor trichostatin A (TSA) on nuclear textural characteristics of human HT1080 fibrosarcoma cells, evaluated by image cytometry, and expression of gelatinases MMP-2 and MMP-9, two metalloproteinases implicated in cancer progression and metastasis.

RESULTS

5-azadC induced significant variation in chromatin higher order organization, particularly chromatin decondensation, associated with reduction in global DNA methylation, concomitantly with increase in MMP-9, and to a lesser extent, MMP-2 expression. TSA alone did not have any effect on HT1080 cells, but exhibited differential activity when added to cells treated with 5-azadC. When treated with both drugs, nuclei had higher texture abnormalities. In this setting, reduction in MMP-9 expression was observed, whereas MMP-2 expression remained unaffected.

CONCLUSIONS

These data show that hypomethylating drug 5-azadC and histone deacetylase inhibitor TSA were able to induce modulation of higher order chromatin organization and gelatinase expression in human HT1080 fibrosarcoma cells.

3-D DNA methylation phenotypes correlate with cytotoxicity levels in prostate and liver cancer cell models

Background

The spatial organization of the genome is being evaluated as a novel indicator of toxicity in conjunction with drug-induced global DNA hypomethylation and concurrent chromatin reorganization. 3D quantitative DNA methylation imaging (3D-qDMI) was applied as a cell-by-cell high-throughput approach to investigate this matter by assessing genome topology through represented immunofluorescent nuclear distribution patterns of 5-methylcytosine (MeC) and global DNA (4,6-diamidino-2-phenylindole = DAPI) in labeled nuclei.

Methods

Differential progression of global DNA hypomethylation was studied by comparatively dosing zebularine (ZEB) and 5-azacytidine (AZA). Treated and untreated (control) human prostate and liver cancer cells were subjected to confocal scanning microscopy and dedicated 3D image analysis for the following features: differential nuclear MeC/DAPI load and codistribution patterns, cell similarity based on these patterns, and corresponding differences in the topology of low-intensity MeC (LIM) and low in intensity DAPI (LID) sites.

Results

Both agents generated a high fraction of similar MeC phenotypes across applied concentrations. ZEB exerted similar effects at 10–100-fold higher drug concentrations than its AZA analogue: concentration-dependent progression of global cytosine demethylation, validated by measuring differential MeC levels in repeat sequences using MethyLight, and the concurrent increase in nuclear LIM densities correlated with cellular growth reduction and cytotoxicity.

Conclusions

3D-qDMI demonstrated the capability of quantitating dose-dependent drug-induced spatial progression of DNA demethylation in cell nuclei, independent from interphase cell-cycle stages and in conjunction with cytotoxicity. The results support the notion of DNA methylation topology being considered as a potential indicator of causal impacts on chromatin distribution with a conceivable application in epigenetic drug toxicology.

Covisualization of methylcytosine, global DNA, and protein biomarkers for In Situ 3D DNA methylation phenotyping of stem cells

DNA methylation and histone modifications are key regulatory mechanisms in cellular differentiation, and are skewed in complex diseases. Therefore, analyzing the higher nuclear organization of methylated DNA in conjunction with relevant cellular components, such as protein biomarkers, may well add cell-by-cell-specific spatial and temporal information to quantitative molecular data for the discovery of stem cell differentiation-related signaling networks and their exploitation in the therapeutic reprogramming of cells. The in situ fluorescent covisualization of methylated DNA (methylated CG dinucleotides = MeC), global DNA (gDNA), and proteins has been challenging, as the immunofluorescence detection of MeC sites requires thorough denaturing of double-stranded DNA for antigen (methylated carbon-5 of cytosine) retrieval. The protocol we present overcomes this obstacle through optimization of cell membrane permeabilization, acid treatment, and intermediate fixation steps to preserve immunostaining of biomarkers and delineate MeC and gDNA, while conserving the captured three-dimensional (3D) structure of the cells; making it suitable for high-resolution confocal microscopy, 3D visualization, and topological analyses of fixed cultured cells as well as fresh and frozen tissue sections.

DNA methylation topology: potential of a chromatin landmark for epigenetic drug toxicology

Targeting chromatin and its basic components through epigenetic drug therapy has become an increased focus in the treatment of complex diseases. This boost calls for the implementation of high-throughput cell-based assays that exploit the increasing knowledge about epigenetic mechanisms and their interventions for genotoxicity testing of epigenetic drugs. 3D quantitative DNA methylation imaging is a novel approach for detecting drug-induced DNA demethylation and concurrent heterochromatin decondensation/reorganization in cells through the analysis of differential nuclear distribution patterns of methylcytosine and gDNA visualized by fluorescence and processed by machine-learning algorithms. Utilizing 3D DNA methylation patterns is a powerful precursor to a series of fully automatable assays that employ chromatin structure and higher organization as novel pharmacodynamic biomarkers for various epigenetic drug actions.

Early in vitro differentiation of mouse definitive endoderm is not correlated with progressive maturation of nuclear DNA methylation patterns

The genome organization in pluripotent cells undergoing the first steps of differentiation is highly relevant to the reprogramming process in differentiation. Considering this fact, chromatin texture patterns that identify cells at the very early stage of lineage commitment could serve as valuable tools in the selection of optimal cell phenotypes for regenerative medicine applications. Here we report on the first-time use of high-resolution three-dimensional fluorescence imaging and comprehensive topological cell-by-cell analyses with a novel image-cytometrical approach towards the identification of in situ global nuclear DNA methylation patterns in early endodermal differentiation of mouse ES cells (up to day 6), and the correlations of these patterns with a set of putative markers for pluripotency and endodermal commitment, and the epithelial and mesenchymal character of cells. Utilizing this in vitro cell system as a model for assessing the relationship between differentiation and nuclear DNA methylation patterns, we found that differentiating cell populations display an increasing number of cells with a gain in DNA methylation load: first within their euchromatin, then extending into heterochromatic areas of the nucleus, which also results in significant changes of methylcytosine/global DNA codistribution patterns. We were also able to co-visualize and quantify the concomitant stochastic marker expression on a per-cell basis, for which we did not measure any correlation to methylcytosine loads or distribution patterns. We observe that the progression of global DNA methylation is not correlated with the standard transcription factors associated with endodermal development. Further studies are needed to determine whether the progression of global methylation could represent a useful signature of cellular differentiation. This concept of tracking epigenetic progression may prove useful in the selection of cell phenotypes for future regenerative medicine applications.

Measuring topology of low-intensity DNA methylation sites for high-throughput assessment of epigenetic drug-induced effects in cancer cells

Epigenetic anti-cancer drugs with demethylating effects have shown to alter genome organization in mammalian cell nuclei. The interest in the development of novel epigenetic drugs has increased the demand for cell-based assays to evaluate drug performance in pre-clinical studies. An imaging-based cytometrical approach that can measure demethylation effects as changes in the spatial nuclear distributions of methylated cytosine and global DNA in cancer cells is introduced in this paper. The cells were studied by immunofluorescence with a specific antibody against 5-methylcytosine (MeC), and 4,6-diamidino-2-phenylindole (DAPI) for delineation of methylated sites and global DNA in nuclei. In the preprocessing step the segmentation of nuclei in three-dimensional images (3-D) is followed by an automated assessment of nuclear DAPI/MeC patterns to exclude dissimilar entities. Next, low-intensity MeC (LIM) and low-intensity DNA (LID) sites of similar nuclei are localized and processed to obtain specific nuclear density profiles. These profiles sampled at half of the total nuclear volume yielded two parameters: LIM(0.5) and LID(0.5). The analysis shows that zebularine and 5-azacytidine-the two tested epigenetic drugs introduce changes in the spatial distribution of low-intensity DNA and MeC signals. LIM(0.5) and LID(0.5) were significantly different (p<0.001) in 5-azacytidine treated (n=660) and zebularine treated (n=496) vs. untreated (n=649) DU145 human prostate cancer cells. In the latter case the LIM sites were predominantly found at the nuclear border, whereas treated populations showed different degrees of increase in LIMs towards the interior nuclear space, in which a large portion of heterochromatin is located. The cell-by-cell evaluation of changes in the spatial reorganization of MeC/DAPI signals revealed that zebularine is a more gentle demethylating agent than 5-azacytidine. Measuring changes in the topology of low-intensity sites can potentially be a valuable component in the high-throughput assessment of demethylation and risk of chromatin reorganization in epigenetic-drug screening tasks.

Evolving generalized Voronoi diagrams for accurate cellular image segmentation

Analyzing cellular morphologies on a cell-by-cell basis is vital for drug discovery, cell biology, and many other biological studies. Interactions between cells in their culture environments cause cells to touch each other in acquired microscopy images. Because of this phenomenon, cell segmentation is a challenging task, especially when the cells are of similar brightness and of highly variable shapes. The concept of topological dependence and the maximum common boundary (MCB) algorithm are presented in our previous work (Yu et al., Cytometry Part A 2009;75A:289-297). However, the MCB algorithm suffers a few shortcomings, such as low computational efficiency and difficulties in generalizing to higher dimensions. To overcome these limitations, we present the evolving generalized Voronoi diagram (EGVD) algorithm. Utilizing image intensity and geometric information, EGVD preserves topological dependence easily in both 2D and 3D images, such that touching cells can be segmented satisfactorily. A systematic comparison with other methods demonstrates that EGVD is accurate and much more efficient.