De novo identification of differentially methylated regions in the human genome

Background: The identification and characterisation of differentially methylated regions (DMRs) between phenotypes in the human genome is of prime interest in epigenetics. We present a novel method, DMRcate, that fits replicated methylation measurements from the Illumina HM450K BeadChip (or 450K array) spatially across the genome using a Gaussian kernel. DMRcate identifies and ranks the most differentially methylated regions across the genome based on tunable kernel smoothing of the differential methylation (DM) signal. The method is agnostic to both genomic annotation and local change in the direction of the DM signal, removes the bias incurred from irregularly spaced methylation sites, and assigns significance to each DMR called via comparison to a null model. Results: We show that, for both simulated and real data, the predictive performance of DMRcate is superior to those of Bumphunter and Probe Lasso, and commensurate with that of comb-p. For the real data, we validate all array-derived DMRs from the candidate methods on a suite of DMRs derived from whole-genome bisulfite sequencing called from the same DNA samples, using two separate phenotype comparisons. Conclusions: The agglomeration of genomically localised individual methylation sites into discrete DMRs is currently best served by a combination of DM-signal smoothing and subsequent threshold specification. The findings also suggest the design of the 450K array shows preference for CpG sites that are more likely to be differentially methylated, but its overall coverage does not adequately reflect the depth and complexity of methylation signatures afforded by sequencing. For the convenience of the research community we have created a user-friendly R software package called DMRcate, downloadable from Bioconductor and compatible with existing preprocessing packages, which allows others to apply the same DMR-finding method on 450K array data

Spatial Organization and Molecular Correlation of Tumor-Infiltrating Lymphocytes Using Deep Learning on Pathology Images

Beyond sample curation and basic pathologic characterization, the digitized H&E-stained images of TCGA samples remain underutilized. To highlight this resource, we present mappings of tumorinfiltrating lymphocytes (TILs) based on H&E images from 13 TCGA tumor types. These TIL maps are derived through computational staining using a convolutional neural network trained to classify patches of images. Affinity propagation revealed local spatial structure in TIL patterns and correlation with overall survival. TIL map structural patterns were grouped using standard histopathological parameters. These patterns are enriched in particular T cell subpopulations derived from molecular measures. TIL densities and spatial structure were differentially enriched among tumor types, immune subtypes, and tumor molecular subtypes, implying that spatial infiltrate state could reflect particular tumor cell aberration states. Obtaining spatial lymphocytic patterns linked to the rich genomic characterization of TCGA samples demonstrates one use for the TCGA image archives with insights into the tumor-immune microenvironment

Disruption of a GATA4/Ankrd1 Signaling Axis in Cardiomyocytes Leads to Sarcomere Disarray: Implications for Anthracycline Cardiomyopathy

Doxorubicin (Adriamycin) is an effective anti-cancer drug, but its clinical usage is limited by a dose-dependent cardiotoxicity characterized by widespread sarcomere disarray and loss of myofilaments. Cardiac ankyrin repeat protein (CARP, ANKRD1) is a transcriptional regulatory protein that is extremely susceptible to doxorubicin; however, the mechanism(s) of doxorubicin-induced CARP depletion and its specific role in cardiomyocytes have not been completely defined. We report that doxorubicin treatment in cardiomyocytes resulted in inhibition of CARP transcription, depletion of CARP protein levels, inhibition of myofilament gene transcription, and marked sarcomere disarray. Knockdown of CARP with small interfering RNA (siRNA) similarly inhibited myofilament gene transcription and disrupted cardiomyocyte sarcomere structure. Adenoviral overexpression of CARP, however, was unable to rescue the doxorubicin-induced sarcomere disarray phenotype. Doxorubicin also induced depletion of the cardiac transcription factor GATA4 in cardiomyocytes. CARP expression is regulated in part by GATA4, prompting us to examine the relationship between GATA4 and CARP in cardiomyocytes. We show in co-transfection experiments that GATA4 operates upstream of CARP by activating the proximal CARP promoter. GATA4-siRNA knockdown in cardiomyocytes inhibited CARP expression and myofilament gene transcription, and induced extensive sarcomere disarray. Adenoviral overexpression of GATA4 (AdV-GATA4) in cardiomyocytes prior to doxorubicin exposure maintained GATA4 levels, modestly restored CARP levels, and attenuated sarcomere disarray. Interestingly, siRNA-mediated depletion of CARP completely abolished the Adv-GATA4 rescue of the doxorubicin-induced sarcomere phenotype. These data demonstrate co-dependent roles for GATA4 and CARP in regulating sarcomere gene expression and maintaining sarcomeric organization in cardiomyocytes in culture. The data further suggests that concurrent depletion of GATA4 and CARP in cardiomyocytes by doxorubicin contributes in large part to myofibrillar disarray and the overall pathophysiology of anthracycline cardiomyopathy

Sarcomeric and nuclear localization of CARP in ARVMs.

<p><b>A:</b> Low magnification (40×) images of ARVMs co-stained with antibodies to CARP (green) and myomesin (red). <b>B:</b> High magnification (100×) of ARVMs immunostained for CARP (green) and myomesin (red, left image) or α-actinin (red, right image), and shown below the image is the corresponding fluorescence intensity along the white arrow in the red and green channels.</p

Doxorubicin inhibits CARP expression at the transcriptional level.

<p><b>A:</b> ARVMs were pretreated with 10 µg/ml cycloheximide (Cyclo), a protein synthesis inhibitor, in the presence or absence of 1 µM doxorubicin (Doxo) and cell lysates analyzed by immunoblot for CARP and actin and corresponding densitometry analysis is shown below. <b>B:</b> Comparison of CARP mRNA decay (quantified by RT-PCR) in ARVMs pretreated with 5 µg/ml actinomycin D (act D) in the presence or absence of 1 µM doxorubicin. <b>C:</b> NRVMs were transfected with a CARP promoter luciferase reporter (CARP-pGL3) and treated with increasing concentrations of doxorubicin. Cell lysates were assayed for luciferase activity and values were normalized to a promoterless control (pGL3 basic). Shown are mean±SD from 4 independent experiments. The qRT-PCR and luciferase-reporter experiments were performed in triplicate. * <i>P</i><0.05, ANOVA.</p

CARP knockdown in ARVMs induces sarcomere disarray.

<p><b>A:</b> ARVMs were transfected with 50 nM nonsilencing siRNA (top panel) or 50 nM CARP-siRNA (bottom panel and fixed at 48 h for immunofluorescence imaging. Cells were stained for CARP (green) myomesin (blue), and filamentous actin (red).</p