The cell type specific 5hmC landscape and dynamics of healthy human hematopoiesis and TET2-mutant pre-leukemia

The conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) is a key step in DNA demethylation that is mediated by ten-eleven-translocation (TET) enzymes, which require ascorbate/vitamin C. Here, we report the 5hmC landscape of normal hematopoiesis and identify cell type-specific 5hmC profiles associated with active transcription and chromatin accessibility of key hematopoietic regulators. We utilized CRISPR/Cas9 to model TET2 loss-of-function mutations in primary human HSPCs. Disrupted cells exhibited increased colonies in serial replating, defective erythroid/megakaryocytic differentiation, and in vivo competitive advantage and myeloid skewing coupled with reduction of 5hmC at erythroid-associated gene loci. Azacitidine and ascorbate restored 5hmC abundance and slowed or reverted the expansion of TET2-mutant clones in vivo. These results demonstrate the key role of 5hmC in normal hematopoiesis and TET2-mutant phenotypes and raise the possibility of utilizing these agents to further our understanding of pre-leukemia/clonal hematopoiesis.

Detection of early stage pancreatic cancer using 5-hydroxymethylcytosine signatures in circulating cell free DNA

Pancreatic cancer is often detected late, when curative therapies are no longer possible. Here, we present non-invasive detection of pancreatic ductal adenocarcinoma (PDAC) by 5-hydroxymethylcytosine (5hmC) changes in circulating cell free DNA from a PDAC cohort (n = 64) in comparison with a non-cancer cohort (n = 243). Differential hydroxymethylation is found in thousands of genes, most significantly in genes related to pancreas development or function (GATA4, GATA6, PROX1, ONECUT1, MEIS2), and cancer pathogenesis (YAP1, TEAD1, PROX1, IGF1). cfDNA hydroxymethylome in PDAC cohort is differentially enriched for genes that are commonly de-regulated in PDAC tumors upon activation of KRAS and inactivation of TP53. Regularized regression models built using 5hmC densities in genes perform with AUC of 0.92 (discovery dataset, n = 79) and 0.92-0.94 (two independent test sets, n = 228). Furthermore, tissue-derived 5hmC features can be used to classify PDAC cfDNA (AUC = 0.88). These findings suggest that 5hmC changes enable classification of PDAC even during early stage disease.

Abstract 1372: Detection of early stage pancreatic cancer using 5–hydroxymethylcytosine signatures in circulating cell free DNA

Pancreatic cancers are typically diagnosed at late stage where disease prognosis is poor as exemplified by a 5-year survival rate of 8.2%. Earlier diagnosis would be beneficial by enabling surgical resection or earlier application of therapeutic regimens. We investigated the detection of pancreatic ductal adenocarcinoma (PDAC) in a non-invasive manner by interrogating changes in 5-hydroxymethylated cytosines (5hmC) in circulating cell free DNA in the plasma of a PDAC cohort (n=51) in comparison with a non-cancer cohort (n=41). 5hmC profiles from PDAC and non-cancer samples were generated using a previously published modified hMe-Seal protocol that utilizes chemical labeling of 5hmC by β-glucosyltransferase and allows detection of cell free 5hmC from small amounts of cfDNA (1). We found that 5hmC sites are enriched in a disease and stage specific manner in exons, 3’UTRs and transcription termination sites. Our data show that 5hmC density is reduced in promoters and histone H3K4me3 associated sites with progressive disease suggesting increased transcriptional activity. 5hmC density is differentially represented in thousands of genes, and a stringently filtered set of the most significant genes points to biology related to pancreas (GATA4, GATA6, PROX1, ONECUT1) and/or cancer development (YAP1, TEAD1, PROX1, ONECUT1, ONECUT2, IGF1 and IGF2). Regularized regression models were built using 5hmC densities in a comprehensive set of genes with the most variable 5hmC counts and performed with an AUC = 0.94 - 0.96 on training data. We tested the ability to classify PDAC and non-cancer samples with the Elastic net and Lasso models on three independent pancreatic cancer 5hmC data sets (n = 26, 23 and 7) compared with corresponding independent non-cancer cohorts (n =103, 53 and 10), and found validation performance to be AUC = 0.74 - 0.97. The findings suggest that 5hmC changes enable classification of PDAC patients with high fidelity and are worthy of further investigation on larger cohorts of patient samples. Reference: 1. Song, C. - X. et al. 5 - Hydroxymethylcytosine signatures in cell-free DNA provide information about tumor types and stages. Cell Res 27, 1231 (2017).

Citation Format: Francois Collin, Yuhong Ning, Gulfem D. Guler, Tierney Phillips, Erin McCarthy, Aaron Scott, Chris Ellison, Chin-Jen Ku, Kim Chau, Alan Ashworth, Stephen R. Quake, Samuel Levy. Detection of early stage pancreatic cancer using 5–hydroxymethylcytosine signatures in circulating cell free DNA [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 1372.

Abstract P4-02-08: Global quantitative measures using next-generation sequencing for breast cancer presence outperform individual tumor markers in plasma

Background: Analytically and clinically validated non-invasive blood tests that quantify breast cancer burden and clinical drug response/resistance are greatly needed. Many groups have successfully detected tumor markers in blood using a variety of technologies, including next generation sequencing (NGS). We performed a comprehensive NGS study on a small number of patients to evaluate the value of global versus individual markers for the quantitation of tumor-derived cell free DNA (cfDNA) in plasma.

Methods: DNA isolated from formalin-fixed primary tumor, buffy coat cells, and plasma from 2 patients with metastatic breast cancer were characterized simultaneously for copy number aberrations (CNAs) and differentially methylated regions (DMRs) using whole genome bisulfite sequencing (WBGS), and targeted sequencing-based genotyping of 346 cancer-associated single nucleotide variations (SNVs). CNA and DMR regions were identified from log normalized, GC content corrected counts and DMR data using Poisson and binomial distribution theory and false discovery rate controlling methods. Percent tumor in cfDNA was estimated from the normalized ratio (plasma: primary tumor) of CNA or DMR compared to buffy coat, aggregating over genomic regions. Sample sets from 8 non-metastatic patients were also profiled using the targeted SNV panel in order to compare SNVs between samples and estimate percent tumor cfDNA.

Results: WGBS detected tumor specific alterations in each primary tumor compared to buffy coat. By analyzing the genome using 100 Kb bins, we observed over 1000 bins with detectable CNA signal and, among 56 million CpG sites, over 30,000 DMRs. As expected, 5 or fewer informative somatic SNVs were detected in each patient. Analysis of these somatic changes in plasma revealed that the tumor fraction estimated from SNV detected in cfDNA varied widely between sites originally discovered in the patient’s primary tumor. In contrast, similar estimates of tumor fraction in cfDNA were obtained using CNA and DMR profiles within each patient; both methods yielded similar estimates of over 50% in one patient and less than 10% in the other. For the patient with high tumor fraction, both CNA and DMR profiles contained examples of individual large genomic regions that displayed additional clear aberrations in the plasma compared to the original tumor, such as a striking loss of a >25 Mb region of chromosome 4.

Conclusions: Although individual somatic SNV in cfDNA can be detected in metastatic disease, calculated allelic fraction based on individual SNVs varies greatly within the same patient. Measuring and integrating CNA or DMR across the genome provided more consistent and reliable estimates of tumor DNA fraction in plasma, and also revealed alterations in plasma from patients with metastatic disease that were not prominent in the primary tumor.

Citation Format: Ellen M Beasley, Richard D Abramson, Gregory E Alexander, David Chan, Kristen Bradley, Francois Collin, Michael Crager, Andrew Dei Rossi, Joseph Dorado, Adam Friedman, William J Gibb, Jennie Jeong, Col Jones, C J Ku, Yan Ma, John Morlan, Kunbin Qu, Aibing Rao, Aaron Scott, Haluk Tezcan. Global quantitative measures using next-generation sequencing for breast cancer presence outperform individual tumor markers in plasma [abstract]. In: Proceedings of the Thirty-Seventh Annual CTRC-AACR San Antonio Breast Cancer Symposium: 2014 Dec 9-13; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2015;75(9 Suppl):Abstract nr P4-02-08.

Identifying network motifs that buffer front-to-back signaling in polarized neutrophils

Neutrophil polarity relies on local, mutual inhibition to segregate incompatible signaling circuits to the leading and trailing edges. Mutual inhibition alone should lead to cells having strong fronts and weak backs or vice versa. However, analysis of cell-to-cell variation in human neutrophils revealed that back polarity remains consistent despite changes in front strength. How is this buffering achieved? Pharmacological perturbations and mathematical modeling revealed a functional role for microtubules in buffering back polarity by mediating positive, long-range crosstalk from front to back; loss of microtubules inhibits buffering and results in anticorrelation between front and back signaling. Furthermore, a systematic, computational search of network topologies found that a long-range, positive front-to-back link is necessary for back buffering. Our studies suggest a design principle that can be employed by polarity networks: short-range mutual inhibition establishes distinct signaling regions, after which directed long-range activation insulates one region from variations in the other.

Network crosstalk dynamically changes during neutrophil polarization

How complex signaling networks shape highly coordinated, multistep cellular responses is poorly understood. Here, we made use of a network-perturbation approach to investigate causal influences, or "crosstalk," among signaling modules involved in the cytoskeletal response of neutrophils to chemoattractant. We quantified the intensity and polarity of cytoskeletal marker proteins over time to characterize stereotyped cellular responses. Analyzing the effects of network disruptions revealed that, not only does crosstalk evolve rapidly during polarization, but also that intensity and polarity responses are influenced by different patterns of crosstalk. Interestingly, persistent crosstalk is arranged in a surprisingly simple circuit: a linear cascade from front to back to microtubules influences intensities, and a feed-forward network in the reverse direction influences polarity. Our approach provided a rational strategy for decomposing a complex, dynamically evolving signaling system and revealed evolving paths of causal influence that shape the neutrophil polarization response.

Patterns of basal signaling heterogeneity can distinguish cellular populations with different drug sensitivities

Phenotypic heterogeneity has been widely observed in cellular populations. However, the extent to which heterogeneity contains biologically or clinically important information is not well understood. Here, we investigated whether patterns of basal signaling heterogeneity, in untreated cancer cell populations, could distinguish cellular populations with different drug sensitivities. We modeled cellular heterogeneity as a mixture of stereotyped signaling states, identified based on colocalization patterns of activated signaling molecules from microscopy images. We found that patterns of heterogeneity could be used to separate the most sensitive and resistant populations to paclitaxel within a set of H460 lung cancer clones and within the NCI-60 panel of cancer cell lines, but not for a set of less heterogeneous, immortalized noncancer human bronchial epithelial cell (HBEC) clones. Our results suggest that patterns of signaling heterogeneity, characterized as ensembles of a small number of distinct phenotypic states, can reveal functional differences among cellular populations.

ON IDENTIFYING INFORMATION FROM IMAGE-BASED SPATIAL POLARITY PHENOTYPES IN NEUTROPHILS

Cell polarity is involved in many biological functions such as development, wound healing and immune responses. In human neutrophils, polarization is characterized by the translocation of distinct sets of signaling molecules to opposite ends of the cell and the rapid rearrangement of cytoskeleton to initiate migration. While many image-based studies have described cellular morphology and the intensity level of polarity signaling molecules, systematic characterization of the spatial distribution of polarity signaling molecules has been lacking. Here we designed a collection of analytical features to quantify spatial phenotypes of polarity molecules. We compared our features to commonly used polarity readouts and found that they captured additional aspects of the polarization dynamics that were not contained in the existing features. Our work provides a starting point to identify informative features for the study of neutrophil polarization.

The distance-profile representation and its application to detection of distantly related protein families

Background

Detecting homology between remotely related protein families is an important problem in computational biology since the biological properties of uncharacterized proteins can often be inferred from those of homologous proteins. Many existing approaches address this problem by measuring the similarity between proteins through sequence or structural alignment. However, these methods do not exploit collective aspects of the protein space and the computed scores are often noisy and frequently fail to recognize distantly related protein families.

Results

We describe an algorithm that improves over the state of the art in homology detection by utilizing global information on the proximity of entities in the protein space. Our method relies on a vectorial representation of proteins and protein families and uses structure-specific association measures between proteins and template structures to form a high-dimensional feature vector for each query protein. These vectors are then processed and transformed to sparse feature vectors that are treated as statistical fingerprints of the query proteins. The new representation induces a new metric between proteins measured by the statistical difference between their corresponding probability distributions.

Conclusion

Using several performance measures we show that the new tool considerably improves the performance in recognizing distant homologies compared to existing approaches such as PSIBLAST and FUGUE.

Epidemiology of otitis media in Chinese children

From May to July 2001, 3013 Chinese children, aged 3-6 years old, from 19 kindergartens in Kaohsiung were examined for the prevalence of otitis media. A senior paediatric otolaryngologist performed examinations. Tympanometric and otoscopic records were combined for assessment of middle ear status. The prevalence rate of otitis media was 9.82% for all children, 5.24% for middle ear effusion and 5.34% for retracted tympanic membrane. The prevalence rates of otitis media in different age groups were as follows: 11.3% for 3-year-olds, 12.4% for 4-year-olds, 11.8% for 5-year-olds, and 6.1% for 6-year-olds. Considering the genetic factors and the previous history of otitis media, there was a significant difference between the normal group and the sick group.