Integrated requirement of non-specific and sequence-specific DNA binding in Myc-driven transcription.

Eukaryotic transcription factors recognize specific DNA sequence motifs, but are also endowed with generic, non-specific DNA-binding activity. How these binding modes are integrated to determine select transcriptional outputs remains unresolved. We addressed this question by site-directed mutagenesis of the Myc transcription factor. Impairment of non-specific DNA backbone contacts caused pervasive loss of genome interactions and gene regulation, associated with increased intra-nuclear mobility of the Myc protein in murine cells. In contrast, a mutant lacking base-specific contacts retained DNA-binding and mobility profiles comparable to those of the wild-type protein, but failed to recognize its consensus binding motif (E-box) and could not activate Myc-target genes. Incidentally, this mutant gained weak affinity for an alternative motif, driving aberrant activation of different genes. Altogether, our data show that non-specific DNA binding is required to engage onto genomic regulatory regions; sequence recognition in turn contributes to transcriptional activation, acting at distinct levels: stabilization and positioning of Myc onto DNA, and-unexpectedly-promotion of its transcriptional activity. Hence, seemingly pervasive genome interaction profiles, as detected by ChIP-seq, actually encompass diverse DNA-binding modalities, driving defined, sequence-dependent transcriptional responses.

PML modulates epigenetic composition of chromatin to regulate expression of pro-metastatic genes in triple-negative breast cancer.

The promyelocytic leukemia (PML) protein organizes nuclear aggregates known as PML nuclear bodies (PML-NBs), where many transcription factors localize to be regulated. In addition, associations of PML and PML-NBs with chromatin are described in various cell types, further implicating PML in transcriptional regulation. However, a complete understanding of the functional consequences of PML association to DNA in cellular contexts where it promotes relevant phenotypes is still lacking. We examined PML chromatin association in triple-negative breast cancer (TNBC) cell lines, where it exerts important oncogenic functions. We find that PML associates discontinuously with large heterochromatic PML-associated domains (PADs) that contain discrete gene-rich euchromatic sub-domains locally depleted of PML. PML promotes heterochromatic organization in PADs and expression of pro-metastatic genes embedded in these sub-domains. Importantly, this occurs outside PML-NBs, suggesting that nucleoplasmic PML exerts a relevant gene regulatory function. We also find that PML plays indirect regulatory roles in TNBC cells by promoting the expression of pro-metastatic genes outside PADs. Our findings suggest that PML is an important transcriptional regulator of pro-oncogenic metagenes in TNBC cells, via transcriptional regulation and epigenetic organization of heterochromatin domains that embed regions of local transcriptional activity.

Joubert syndrome-derived induced pluripotent stem cells show altered neuronal differentiation in vitro.

Joubert syndrome (JS) is a recessively inherited congenital ataxia characterized by hypotonia, psychomotor delay, abnormal ocular movements, intellectual disability, and a peculiar cerebellar and brainstem malformation, the "molar tooth sign." Over 40 causative genes have been reported, all encoding for proteins implicated in the structure or functioning of the primary cilium, a subcellular organelle widely present in embryonic and adult tissues. In this paper, we developed an in vitro neuronal differentiation model using patient-derived induced pluripotent stem cells (iPSCs), to evaluate possible neurodevelopmental defects in JS. To this end, iPSCs from four JS patients harboring mutations in distinct JS genes (AHI1, CPLANE1, TMEM67, and CC2D2A) were differentiated alongside healthy control cells to obtain mid-hindbrain precursors and cerebellar granule cells. Differentiation was monitored over 31 days through the detection of lineage-specific marker expression by qRT-PCR, immunofluorescence, and transcriptomics analysis. All JS patient-derived iPSCs, regardless of the mutant gene, showed a similar impairment to differentiate into mid-hindbrain and cerebellar granule cells when compared to healthy controls. In addition, analysis of primary cilium count and morphology showed notable ciliary defects in all differentiating JS patient-derived iPSCs compared to controls. These results confirm that patient-derived iPSCs are an accessible and relevant in vitro model to analyze cellular phenotypes connected to the presence of JS gene mutations in a neuronal context.

Accurate and highly interpretable prediction of gene expression from histone modifications.

BACKGROUND: Histone Mark Modifications (HMs) are crucial actors in gene regulation, as they actively remodel chromatin to modulate transcriptional activity: aberrant combinatorial patterns of HMs have been connected with several diseases, including cancer. HMs are, however, reversible modifications: understanding their role in disease would allow the design of 'epigenetic drugs' for specific, non-invasive treatments. Standard statistical techniques were not entirely successful in extracting representative features from raw HM signals over gene locations. On the other hand, deep learning approaches allow for effective automatic feature extraction, but at the expense of model interpretation. RESULTS: Here, we propose ShallowChrome, a novel computational pipeline to model transcriptional regulation via HMs in both an accurate and interpretable way. We attain state-of-the-art results on the binary classification of gene transcriptional states over 56 cell-types from the REMC database, largely outperforming recent deep learning approaches. We interpret our models by extracting insightful gene-specific regulative patterns, and we analyse them for the specific case of the PAX5 gene over three differentiated blood cell lines. Finally, we compare the patterns we obtained with the characteristic emission patterns of ChromHMM, and show that ShallowChrome is able to coherently rank groups of chromatin states w.r.t. their transcriptional activity. CONCLUSIONS: In this work we demonstrate that it is possible to model HM-modulated gene expression regulation in a highly accurate, yet interpretable way. Our feature extraction algorithm leverages on data downstream the identification of enriched regions to retrieve gene-wise, statistically significant and dynamically located features for each HM. These features are highly predictive of gene transcriptional state, and allow for accurate modeling by computationally efficient logistic regression models. These models allow a direct inspection and a rigorous interpretation, helping to formulate quantifiable hypotheses.

Innovative digital solution supporting sun protection and vitamin D synthesis by using satellite-based monitoring of solar radiation.

Public health campaigns advise minimising UV radiation (UVR) exposure to prevent skin cancer and precancer, e.g. actinic keratosis (AK). A 3-day clinical field study, in Brazil, was performed to evaluate the mobile app Sun4Health® by siHealth Ltd. The app performs real-time monitoring of both erythemal and vitamin D-effective solar radiation doses using satellite data, enabling personalised recommendations on optimal sun exposure time and sunscreen use. When coupled to a wearable device, the app also provides body-site specific recommendations ("3D" version). 59 healthy volunteers were randomised into 3 groups, each given a different app providing: (1) ultraviolet index only (control app), (2) personalised recommendations and sun overexposure alerts (Sun4Health® app), (3) as (2) but connected via Bluetooth to a wearable device to monitor sun exposure in 3D (Sun4Health®-3D app). Participants were offered sunscreens (SPF 30 and 50) to use at their discretion. Erythema, quantified by reflectance spectroscopy, was assessed daily in the mornings and evenings on six body sites. Serum vitamin D (25(OH)D3) was measured before and after the study. Mean increase of erythema (Mexameter® units ± SD) of all exposed body sites combined over 3 days showed 55.76 ± 47.47 for group 1, 40.27 ± 37.91 for group 2 and 37.12 ± 30.69 for group 3 (p < 0.05 for all groups). Mean increase of serum 25(OH)D3 (nmol/l ± SD) showed 1.32 ± 36.49 for group 1, 6.38 ± 21.19 for group 2 and 18.68 ± 35.45 for group 3 (p > 0.05 for all groups). The results show that the Sun4Health® app is safe to use and can modify behaviour to reduce skin erythema (sunburn) yet not decreasing vitamin D status.

Cooperation Between MYC and ß-Catenin in Liver Tumorigenesis Requires Yap/Taz.

BACKGROUND AND AIMS: Activation of MYC and catenin beta-1 (CTNNB1, encoding ß-catenin) can co-occur in liver cancer, but how these oncogenes cooperate in tumorigenesis remains unclear. APPROACH AND RESULTS: We generated a mouse model allowing conditional activation of MYC and WNT/ß-catenin signaling (through either ß-catenin activation or loss of APC - adenomatous polyposis coli) upon expression of CRE recombinase in the liver and monitored their effects on hepatocyte proliferation, apoptosis, gene expression profiles, and tumorigenesis. Activation of WNT/ß-catenin signaling strongly accelerated MYC-driven carcinogenesis in the liver. Both pathways also cooperated in promoting cellular transformation in vitro, demonstrating their cell-autonomous action. Short-term induction of MYC and ß-catenin in hepatocytes, followed by RNA-sequencing profiling, allowed the identification of a "Myc/ß-catenin signature," composed of a discrete set of Myc-activated genes whose expression increased in the presence of active ß-catenin. Notably, this signature enriched for targets of Yes-associated protein (Yap) and transcriptional coactivator with PDZ-binding motif (Taz), two transcriptional coactivators known to be activated by WNT/ß-catenin signaling and to cooperate with MYC in mitogenic activation and liver transformation. Consistent with these regulatory connections, Yap/Taz accumulated upon Myc/ß-catenin activation and were required not only for the ensuing proliferative response, but also for tumor cell growth and survival. Finally, the Myc/ß-catenin signature was enriched in a subset of human hepatocellular carcinomas characterized by comparatively poor prognosis. CONCLUSIONS: Myc and ß-catenin show a strong cooperative action in liver carcinogenesis, with Yap and Taz serving as mediators of this effect. These findings warrant efforts toward therapeutic targeting of Yap/Taz in aggressive liver tumors marked by elevated Myc/ß-catenin activity.

An early Myc-dependent transcriptional program orchestrates cell growth during B-cell activation.

Upon activation, lymphocytes exit quiescence and undergo substantial increases in cell size, accompanied by activation of energy-producing and anabolic pathways, widespread chromatin decompaction, and elevated transcriptional activity. These changes depend upon prior induction of the Myc transcription factor, but how Myc controls them remains unclear. We addressed this issue by profiling the response to LPS stimulation in wild-type and c-myc-deleted primary mouse B-cells. Myc is rapidly induced, becomes detectable on virtually all active promoters and enhancers, but has no direct impact on global transcriptional activity. Instead, Myc contributes to the swift up- and down-regulation of several hundred genes, including many known regulators of the aforementioned cellular processes. Myc-activated promoters are enriched for E-box consensus motifs, bind Myc at the highest levels, and show enhanced RNA Polymerase II recruitment, the opposite being true at down-regulated loci. Remarkably, the Myc-dependent signature identified in activated B-cells is also enriched in Myc-driven B-cell lymphomas: hence, besides modulation of new cancer-specific programs, the oncogenic action of Myc may largely rely on sustained deregulation of its normal physiological targets.

Chest CT in the emergency department for suspected COVID-19 pneumonia.

PURPOSE: In overwhelmed emergency departments (EDs) facing COVID-19 outbreak, a swift diagnosis is imperative. CT role was widely debated for its limited specificity. Here we report the diagnostic role of CT in two EDs in Lombardy, epicenter of Italian outbreak. MATERIAL AND METHODS: Admitting chest CT from 142 consecutive patients with suspected COVID-19 were retrospectively analyzed. CT scans were classified in "highly likely," "likely," and "unlikely" COVID-19 pneumonia according to the presence of typical, indeterminate, and atypical findings, or "negative" in the absence of findings, or "alternative diagnosis" when a different diagnosis was found. Nasopharyngeal swab results, turnaround time, and time to positive results were collected. CT diagnostic performances were assessed considering RT-PCR as reference standard. RESULTS: Most of cases (96/142, 68%) were classified as "highly likely" COVID-19 pneumonia. Ten (7%) and seven (5%) patients were classified as "likely" and "unlikely" COVID-19 pneumonia, respectively. In 21 (15%) patients a differential diagnosis was provided, including typical pneumonia, pulmonary edema, neoplasia, and pulmonary embolism. CT was negative in 8/142 (6%) patients. Mean turnaround time for the first COVID-19 RT-PCR was 30 ± 13 h. CT diagnostic accuracy in respect of the first test swab was 79% and increased to 91.5% after repeated swabs and/or BAL, for 18 false-negative first swab. CT performance was good with 76% specificity, 99% sensitivity, 90% positive predictive value and 97% negative predictive value. CONCLUSION: Chest CT was useful to streamline patients' triage while waiting for RT-PCR in the ED, supporting the clinical suspicion of COVID-19 or providing alternative diagnosis.

Integrative analysis of RNA polymerase II and transcriptional dynamics upon MYC activation.

Overexpression of the MYC transcription factor causes its widespread interaction with regulatory elements in the genome but leads to the up- and down-regulation of discrete sets of genes. The molecular determinants of these selective transcriptional responses remain elusive. Here, we present an integrated time-course analysis of transcription and mRNA dynamics following MYC activation in proliferating mouse fibroblasts, based on chromatin immunoprecipitation, metabolic labeling of newly synthesized RNA, extensive sequencing, and mathematical modeling. Transcriptional activation correlated with the highest increases in MYC binding at promoters. Repression followed a reciprocal scenario, with the lowest gains in MYC binding. Altogether, the relative abundance (henceforth, "share") of MYC at promoters was the strongest predictor of transcriptional responses in diverse cell types, predominating over MYC's association with the corepressor ZBTB17 (also known as MIZ1). MYC activation elicited immediate loading of RNA polymerase II (RNAPII) at activated promoters, followed by increases in pause-release, while repressed promoters showed opposite effects. Gains and losses in RNAPII loading were proportional to the changes in the MYC share, suggesting that repression by MYC may be partly indirect, owing to competition for limiting amounts of RNAPII. Secondary to the changes in RNAPII loading, the dynamics of elongation and pre-mRNA processing were also rapidly altered at MYC regulated genes, leading to the transient accumulation of partially or aberrantly processed mRNAs. Altogether, our results shed light on how overexpressed MYC alters the various phases of the RNAPII cycle and the resulting transcriptional response.

Selective transcriptional regulation by Myc in cellular growth control and lymphomagenesis.

The c-myc proto-oncogene product, Myc, is a transcription factor that binds thousands of genomic loci. Recent work suggested that rather than up- and downregulating selected groups of genes, Myc targets all active promoters and enhancers in the genome (a phenomenon termed 'invasion') and acts as a general amplifier of transcription. However, the available data did not readily discriminate between direct and indirect effects of Myc on RNA biogenesis. We addressed this issue with genome-wide chromatin immunoprecipitation and RNA expression profiles during B-cell lymphomagenesis in mice, in cultured B cells and fibroblasts. Consistent with long-standing observations, we detected general increases in total RNA or messenger RNA copies per cell (hereby termed 'amplification') when comparing actively proliferating cells with control quiescent cells: this was true whether cells were stimulated by mitogens (requiring endogenous Myc for a proliferative response) or by deregulated, oncogenic Myc activity. RNA amplification and promoter/enhancer invasion by Myc were separable phenomena that could occur without one another. Moreover, whether or not associated with RNA amplification, Myc drove the differential expression of distinct subsets of target genes. Hence, although having the potential to interact with all active or poised regulatory elements in the genome, Myc does not directly act as a global transcriptional amplifier. Instead, our results indicate that Myc activates and represses transcription of discrete gene sets, leading to changes in cellular state that can in turn feed back on global RNA production and turnover.

MuSERA: Multiple Sample Enriched Region Assessment.

Enriched region (ER) identification is a fundamental step in several next-generation sequencing (NGS) experiment types. Yet, although NGS experimental protocols recommend producing replicate samples for each evaluated condition and their consistency is usually assessed, typically pipelines for ER identification do not consider available NGS replicates. This may alter genome-wide descriptions of ERs, hinder significance of subsequent analyses on detected ERs and eventually preclude biological discoveries that evidence in replicate could support. MuSERA is a broadly useful stand-alone tool for both interactive and batch analysis of combined evidence from ERs in multiple ChIP-seq or DNase-seq replicates. Besides rigorously combining sample replicates to increase statistical significance of detected ERs, it also provides quantitative evaluations and graphical features to assess the biological relevance of each determined ER set within its genomic context; they include genomic annotation of determined ERs, nearest ER distance distribution, global correlation assessment of ERs and an integrated genome browser. We review MuSERA rationale and implementation, and illustrate how sets of significant ERs are expanded by applying MuSERA on replicates for several types of NGS data, including ChIP-seq of transcription factors or histone marks and DNase-seq hypersensitive sites. We show that MuSERA can determine a new, enhanced set of ERs for each sample by locally combining evidence on replicates, and prove how the easy-to-use interactive graphical displays and quantitative evaluations that MuSERA provides effectively support thorough inspection of obtained results and evaluation of their biological content, facilitating their understanding and biological interpretations. MuSERA is freely available at http://www.bioinformatics.deib.polimi.it/MuSERA/.

Cardiosomal microRNAs Are Essential in Post-Infarction Myofibroblast Phenoconversion.

The inclusion of microRNAs (miRNAs) in extracellular microvesicles/exosomes (named cardiosomes when deriving from cardiomyocytes) allows their active transportation and ensures cell-cell communication. We hypothesize that cardiosomal miRNAs play a pivotal role in the activation of myofibroblasts following ischemic injury. Using a murine model of myocardial infarction (MI), we tested our hypothesis by measuring in isolated fibroblasts and cardiosomes the expression levels of a set of miRNAs, which are upregulated in cardiomyocytes post-MI and involved in myofibroblast phenoconversion. We found that miR-195 was significantly upregulated in cardiosomes and in fibroblasts isolated after MI compared with SHAM conditions. Moreover, primary isolated cardiac fibroblasts were activated both when incubated with cardiosomes isolated from ischemic cardiomyocytes and when cultured in conditioned medium of post-MI cardiomyocytes, whereas no significant effect was observed following incubation with cardiosomes or medium from sham cardiomyocytes. Taken together, our findings indicate for the first time that a cardiomyocyte-specific miRNA, transferred to fibroblasts in form of exosomal cargo, is crucial in the activation of myofibroblasts.

FunChIP: an R/Bioconductor package for functional classification of ChIP-seq shapes.

SUMMARY: Chromatin Immunoprecipitation followed by sequencing (ChIP-seq) generates local accumulations of sequencing reads on the genome ("peaks"), which correspond to specific protein-DNA interactions or chromatin modifications. Peaks are detected by considering their total area above a background signal, usually neglecting their shapes, which instead may convey additional biological information. We present FunChIP, an R/Bioconductor package for clustering peaks according to a functional representation of their shapes: after approximating their profiles with cubic B-splines, FunChIP minimizes their functional distance and classifies the peaks applying a k-mean alignment and clustering algorithm. The whole pipeline is user-friendly and provides visualization functions for a quick inspection of the results. An application to the transcription factor Myc in 3T9 murine fibroblasts shows that clusters of peaks with different shapes are associated with different genomic locations and different transcriptional regulatory activity. AVAILABILITY AND IMPLEMENTATION: The package is implemented in R and is available under Artistic Licence 2.0 from the Bioconductor website (http://bioconductor.org/packages/FunChIP). CONTACT: marco.morelli@iit.it. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Mutual epithelium-macrophage dependency in liver carcinogenesis mediated by ST18.

The ST18 gene has been proposed to act either as a tumor suppressor or as an oncogene in different human cancers, but direct evidence for its role in tumorigenesis has been lacking thus far. Here, we demonstrate that ST18 is critical for tumor progression and maintenance in a mouse model of liver cancer, based on oncogenic transformation and adoptive transfer of primary precursor cells (hepatoblasts). ST18 messenger RNA (mRNA) and protein were detectable neither in normal liver nor in cultured hepatoblasts, but were readily expressed after subcutaneous engraftment and tumor growth. ST18 expression in liver cells was induced by inflammatory cues, including acute or chronic inflammation in vivo, as well as coculture with macrophages in vitro. Knocking down the ST18 mRNA in transplanted hepatoblasts delayed tumor progression. Induction of ST18 knockdown in pre-established tumors caused rapid tumor involution associated with pervasive morphological changes, proliferative arrest, and apoptosis in tumor cells, as well as depletion of tumor-associated macrophages, vascular ectasia, and hemorrhage. Reciprocally, systemic depletion of macrophages in recipient animals had very similar phenotypic consequences, impairing either tumor development or maintenance, and suppressing ST18 expression in hepatoblasts. Finally, RNA sequencing of ST18-depleted tumors before involution revealed down-regulation of inflammatory response genes, pointing to the suppression of nuclear factor kappa B-dependent transcription. CONCLUSION: ST18 expression in epithelial cells is induced by tumor-associated macrophages, contributing to the reciprocal feed-forward loop between both cell types in liver tumorigenesis. Our findings warrant the exploration of means to interfere with ST18-dependent epithelium-macrophage interactions in a therapeutic setting. (Hepatology 2017;65:1708-1719).

Distinct gene expression profiles associated with Notch ligands Delta-like 4 and Jagged1 in plaque material from peripheral artery disease patients: a pilot study.

BACKGROUND: The lack of early diagnosis, progression markers and effective pharmacological treatment has dramatic unfavourable effects on clinical outcomes in patients with peripheral artery disease (PAD). Addressing these issues will require dissecting the molecular mechanisms underlying this disease. We sought to characterize the Notch signaling and atherosclerosis relevant markers in lesions from femoral arteries of symptomatic PAD patients. METHODS: Plaque material from the common femoral, superficial femoral or popliteal arteries of 20 patients was removed by directional atherectomy. RNA was obtained from 9 out of 20 samples and analysed by quantitative reverse transcriptase-polymerase chain reaction (qRT-PCR). RESULTS: We detected expression of Notch ligands Delta-like 4 (Dll4) and Jagged1 (Jag1), of Notch target genes Hes1, Hey1, Hey2, HeyL and of markers of plaque inflammation and stability such as vascular cell adhesion molecule 1 (VCAM1), smooth muscle 22 (SM22), cyclooxygenase 2 (COX2), Bcl2, CD68 and miRNAs 21-5p, 125a-5p, 126-5p,146-5p, 155-5p, 424-5p. We found an "inflamed plaque" gene expression profile characterized by high Dll4 associated to medium/high CD68, COX2, VCAM1, Hes1, miR126-5p, miR146a-5p, miR155-5p, miR424-5p and low Jag1, SM22, Bcl2, Hey2, HeyL, miR125a-5p (2/9 patients) and a "stable plaque" profile characterized by high Jag1 associated to medium/high Hey2, HeyL, SM22, Bcl2, miR125a and low Dll4, CD68, COX2, VCAM1, miR126-5p, miR146a-5p, miR155-5p, miR424-5p (3/9 patients). The remaining patients (4/9) showed a plaque profile with intermediate characteristics. CONCLUSIONS: This study reveals the existence of a gene signature associated to Notch activation by specific ligands that could be predictive of PAD progression.

Serum From Advanced Heart Failure Patients Promotes Angiogenic Sprouting and Affects the Notch Pathway in Human Endothelial Cells.

It is unknown whether components present in heart failure (HF) patients' serum provide an angiogenic stimulus. We sought to determine whether serum from HF patients affects angiogenesis and its major modulator, the Notch pathway, in human umbilical vein endothelial cells (HUVECs). In cells treated with serum from healthy subjects or from patients at different HF stage we determined: (1) Sprouting angiogenesis, by measuring cells network (closed tubes) in collagen gel. (2) Protein levels of Notch receptors 1, 2, 4, and ligands Jagged1, Delta-like4. We found a higher number of closed tubes in HUVECs treated with advanced HF patients serum in comparison with cells treated with serum from mild HF patients or controls. Furthermore, as indicated by the reduction of the active form of Notch4 (N4IC) and of Jagged1, advanced HF patients serum inhibited Notch signalling in HUVECs in comparison with mild HF patients' serum and controls. The circulating levels of NT-proBNP (N-terminal of the pro-hormone brain natriuretic peptide), a marker for the detection and evalutation of HF, were positively correlated with the number of closed tubes (r = 0.485) and negatively with Notch4IC and Jagged1 levels in sera-treated cells (r = -0.526 and r = -0.604, respectively). In conclusion, we found that sera from advanced HF patients promote sprouting angiogenesis and dysregulate Notch signaling in HUVECs. Our study provides in vitro evidence of an angiogenic stimulus arising during HF progression and suggests a role for the Notch pathway in it. J. Cell. Physiol. 231: 2700-2710, 2016. © 2016 Wiley Periodicals, Inc.

Using combined evidence from replicates to evaluate ChIP-seq peaks.

MOTIVATION: Chromatin Immunoprecipitation followed by sequencing (ChIP-seq) detects genome-wide DNA-protein interactions and chromatin modifications, returning enriched regions (ERs), usually associated with a significance score. Moderately significant interactions can correspond to true, weak interactions, or to false positives; replicates of a ChIP-seq experiment can provide co-localised evidence to decide between the two cases. We designed a general methodological framework to rigorously combine the evidence of ERs in ChIP-seq replicates, with the option to set a significance threshold on the repeated evidence and a minimum number of samples bearing this evidence. RESULTS: We applied our method to Myc transcription factor ChIP-seq datasets in K562 cells available in the ENCODE project. Using replicates, we could extend up to 3 times the ER number with respect to single-sample analysis with equivalent significance threshold. We validated the 'rescued' ERs by checking for the overlap with open chromatin regions and for the enrichment of the motif that Myc binds with strongest affinity; we compared our results with alternative methods (IDR and jMOSAiCS), obtaining more validated peaks than the former and less peaks than latter, but with a better validation. AVAILABILITY AND IMPLEMENTATION: An implementation of the proposed method and its source code under GPLv3 license are freely available at http://www.bioinformatics.deib.polimi.it/MSPC/ and http://mspc.codeplex.com/, respectively. CONTACT: marco.morelli@iit.it SUPPLEMENTARY INFORMATION: Supplementary Material are available at Bioinformatics online.

INSPEcT: a computational tool to infer mRNA synthesis, processing and degradation dynamics from RNA- and 4sU-seq time course experiments.

MOTIVATION: Cellular mRNA levels originate from the combined action of multiple regulatory processes, which can be recapitulated by the rates of pre-mRNA synthesis, pre-mRNA processing and mRNA degradation. Recent experimental and computational advances set the basis to study these intertwined levels of regulation. Nevertheless, software for the comprehensive quantification of RNA dynamics is still lacking. RESULTS: INSPEcT is an R package for the integrative analysis of RNA- and 4sU-seq data to study the dynamics of transcriptional regulation. INSPEcT provides gene-level quantification of these rates, and a modeling framework to identify which of these regulatory processes are most likely to explain the observed mRNA and pre-mRNA concentrations. Software performance is tested on a synthetic dataset, instrumental to guide the choice of the modeling parameters and the experimental design. AVAILABILITY AND IMPLEMENTATION: INSPEcT is submitted to Bioconductor and is currently available as Supplementary Additional File S1. CONTACT: mattia.pelizzola@iit.it SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Speckle-tracking global longitudinal strain as an early predictor of cardiotoxicity in breast carcinoma.

PURPOSE: Recent development in anticancer therapies for breast carcinoma allowed an improvement in patients' survival, notwithstanding a parallel increase of cardiovascular morbidity. Cardiotoxicity of anticancer therapies represents a relevant problem due to its insidious onset and potentially irreversible cardiac damage. The aim of the present study was to test whether 2D speckle tracking analysis can help in predicting overt systolic dysfunction. METHODS: A "real world" cohort of 69 patients with breast carcinoma undergoing adjuvant and/or neo-adjuvant chemotherapy was tested 2D-speckle tracking analysis before the beginning of chemotherapy and every 3 months for 1 year. Clinical data, 12-lead ECGs, and lab tests were collected according to the same visit protocol. RESULTS: Over 1-year follow-up, 19 patients (27 %) developed cardiac dysfunction according to the CREC criteria, with an average onset time from enrolment of 6.8 months. A global longitudinal strain (GLS) threshold ≥-16 % at 3 months from chemotherapy was able to predict subsequent systolic dysfunction development with high sensitivity (80 %) and specificity (90 %) and a negative predictive value of 92 %. After the introduction of cardioprotective drugs, left ventricular ejection fraction (LVEF) progressively recovered, while alterations of GLS persisted at 1-year follow-up. CONCLUSIONS: Strain imaging with 2D speckle tracking allows the identification of patients at low-risk for chemotherapy-related systolic dysfunction and can help optimizing resources allocations and improving follow-up quality. GLS can also provide a more accurate prognostic index of resolved systolic dysfunction when compared to standard LVEF.

methylPipe and compEpiTools: a suite of R packages for the integrative analysis of epigenomics data.

BACKGROUND: Numerous methods are available to profile several epigenetic marks, providing data with different genome coverage and resolution. Large epigenomic datasets are then generated, and often combined with other high-throughput data, including RNA-seq, ChIP-seq for transcription factors (TFs) binding and DNase-seq experiments. Despite the numerous computational tools covering specific steps in the analysis of large-scale epigenomics data, comprehensive software solutions for their integrative analysis are still missing. Multiple tools must be identified and combined to jointly analyze histone marks, TFs binding and other -omics data together with DNA methylation data, complicating the analysis of these data and their integration with publicly available datasets. RESULTS: To overcome the burden of integrating various data types with multiple tools, we developed two companion R/Bioconductor packages. The former, methylPipe, is tailored to the analysis of high- or low-resolution DNA methylomes in several species, accommodating (hydroxy-)methyl-cytosines in both CpG and non-CpG sequence context. The analysis of multiple whole-genome bisulfite sequencing experiments is supported, while maintaining the ability of integrating targeted genomic data. The latter, compEpiTools, seamlessly incorporates the results obtained with methylPipe and supports their integration with other epigenomics data. It provides a number of methods to score these data in regions of interest, leading to the identification of enhancers, lncRNAs, and RNAPII stalling/elongation dynamics. Moreover, it allows a fast and comprehensive annotation of the resulting genomic regions, and the association of the corresponding genes with non-redundant GeneOntology terms. Finally, the package includes a flexible method based on heatmaps for the integration of various data types, combining annotation tracks with continuous or categorical data tracks. CONCLUSIONS: methylPipe and compEpiTools provide a comprehensive Bioconductor-compliant solution for the integrative analysis of heterogeneous epigenomics data. These packages are instrumental in providing biologists with minimal R skills a complete toolkit facilitating the analysis of their own data, or in accelerating the analyses performed by more experienced bioinformaticians.