Targeting ATP2B1 impairs PI3K/Akt/FOXO signaling and reduces SARS-COV-2 infection and replication.

ATP2B1 is a known regulator of calcium (Ca2+) cellular export and homeostasis. Diminished levels of intracellular Ca2+ content have been suggested to impair SARS-CoV-2 replication. Here, we demonstrate that a nontoxic caloxin-derivative compound (PI-7) reduces intracellular Ca2+ levels and impairs SARS-CoV-2 infection. Furthermore, a rare homozygous intronic variant of ATP2B1 is shown to be associated with the severity of COVID-19. The mechanism of action during SARS-CoV-2 infection involves the PI3K/Akt signaling pathway activation, inactivation of FOXO3 transcription factor function, and subsequent transcriptional inhibition of the membrane and reticulum Ca2+ pumps ATP2B1 and ATP2A1, respectively. The pharmacological action of compound PI-7 on sustaining both ATP2B1 and ATP2A1 expression reduces the intracellular cytoplasmic Ca2+ pool and thus negatively influences SARS-CoV-2 replication and propagation. As compound PI-7 lacks toxicity in vitro, its prophylactic use as a therapeutic agent against COVID-19 is envisioned here.

8-oxodG accumulation within super-enhancers marks fragile CTCF-mediated chromatin loops.

8-Oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG), a major product of the DNA oxidization process, has been proposed to have an epigenetic function in gene regulation and has been associated with genome instability. NGS-based methodologies are contributing to the characterization of the 8-oxodG function in the genome. However, the 8-oxodG epigenetic role at a genomic level and the mechanisms controlling the genomic 8-oxodG accumulation/maintenance have not yet been fully characterized. In this study, we report the identification and characterization of a set of enhancer regions accumulating 8-oxodG in human epithelial cells. We found that these oxidized enhancers are mainly super-enhancers and are associated with bidirectional-transcribed enhancer RNAs and DNA Damage Response activation. Moreover, using ChIA-PET and HiC data, we identified specific CTCF-mediated chromatin loops in which the oxidized enhancer and promoter regions physically associate. Oxidized enhancers and their associated chromatin loops accumulate endogenous double-strand breaks which are in turn repaired by NHEJ pathway through a transcription-dependent mechanism. Our work suggests that 8-oxodG accumulation in enhancers-promoters pairs occurs in a transcription-dependent manner and provides novel mechanistic insights on the intrinsic fragility of chromatin loops containing oxidized enhancers-promoters interactions.

Development of a multi-targeted chemotherapeutic approach based on G-quadruplex stabilisation and carbonic anhydrase inhibition.

A novel class of compounds designed to hit two anti-tumour targets, G-quadruplex structures and human carbonic anhydrases (hCAs) IX and XII is proposed. The induction/stabilisation of G-quadruplex structures by small molecules has emerged as an anticancer strategy, disrupting telomere maintenance and reducing oncogene expression. hCAs IX and XII are well-established anti-tumour targets, upregulated in many hypoxic tumours and contributing to metastasis. The ligands reported feature a berberine G-quadruplex stabiliser scaffold connected to a moiety inhibiting hCAs IX and XII. In vitro experiments showed that our compounds selectively stabilise G-quadruplex structures and inhibit hCAs IX and XII. The crystal structure of a telomeric G-quadruplex in complex with one of these ligands was obtained, shedding light on the ligand/target interaction mode. The most promising ligands showed significant cytotoxicity against CA IX-positive HeLa cancer cells in hypoxia, and the ability to stabilise G-quadruplexes within tumour cells.

The Intertwined Role of 8-oxodG and G4 in Transcription Regulation.

The guanine base in nucleic acids is, among the other bases, the most susceptible to being converted into 8-Oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG) when exposed to reactive oxygen species. In double-helix DNA, 8-oxodG can pair with adenine; hence, it may cause a G > T (C > A) mutation; it is frequently referred to as a form of DNA damage and promptly corrected by DNA repair mechanisms. Moreover, 8-oxodG has recently been redefined as an epigenetic factor that impacts transcriptional regulatory elements and other epigenetic modifications. It has been proposed that 8-oxodG exerts epigenetic control through interplay with the G-quadruplex (G4), a non-canonical DNA structure, in transcription regulatory regions. In this review, we focused on the epigenetic roles of 8-oxodG and the G4 and explored their interplay at the genomic level.

Genome-wide mapping of genomic DNA damage: methods and implications.

Exposures from the external and internal environments lead to the modification of genomic DNA, which is implicated in the cause of numerous diseases, including cancer, cardiovascular, pulmonary and neurodegenerative diseases, together with ageing. However, the precise mechanism(s) linking the presence of damage, to impact upon cellular function and pathogenesis, is far from clear. Genomic location of specific forms of damage is likely to be highly informative in understanding this process, as the impact of downstream events (e.g. mutation, microsatellite instability, altered methylation and gene expression) on cellular function will be positional-events at key locations will have the greatest impact. However, until recently, methods for assessing DNA damage determined the totality of damage in the genomic location, with no positional information. The technique of "mapping DNA adductomics" describes the molecular approaches that map a variety of forms of DNA damage, to specific locations across the nuclear and mitochondrial genomes. We propose that integrated comparison of this information with other genome-wide data, such as mutational hotspots for specific genotoxins, tumour-specific mutation patterns and chromatin organisation and transcriptional activity in non-cancerous lesions (such as nevi), pre-cancerous conditions (such as polyps) and tumours, will improve our understanding of how environmental toxins lead to cancer. Adopting an analogous approach for non-cancer diseases, including the development of genome-wide assays for other cellular outcomes of DNA damage, will improve our understanding of the role of DNA damage in pathogenesis more generally.

The genomic landscape of 8-oxodG reveals enrichment at specific inherently fragile promoters.

8-Oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG) is the most common marker of oxidative stress and its accumulation within the genome has been associated with major human health issues such as cancer, aging, cardiovascular and neurodegenerative diseases. The characterization of the different genomic sites where 8-oxodG accumulates and the mechanisms underlying its formation are still poorly understood. Using OxiDIP-seq, we recently derived the genome-wide distribution of 8-oxodG in human non-tumorigenic epithelial breast cells (MCF10A). Here, we identify a subset of human promoters that accumulate 8-oxodG under steady-state condition. 8-oxodG nucleotides co-localize with double strand breaks (DSBs) at bidirectional and CG skewed promoters and their density correlate with RNA Polymerase II co-occupancy and transcription. Furthermore, by performing OxiDIP-seq in quiescent (G0) cells, we found a strong reduction of oxidatively-generated damage in the majority of 8-oxodG-positive promoters in the absence of DNA replication. Overall, our results suggest that the accumulation of 8-oxodG at gene promoters occurs through DNA replication-dependent or -independent mechanisms, with a possible contribution to the formation of cancer-associated translocation events.

Genome-wide mapping of 8-oxo-7,8-dihydro-2'-deoxyguanosine reveals accumulation of oxidatively-generated damage at DNA replication origins within transcribed long genes of mammalian cells.

8-Oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG) is one of the major DNA modifications and a potent pre-mutagenic lesion prone to mispair with 2'-deoxyadenosine (dA). Several thousand residues of 8-oxodG are constitutively generated in the genome of mammalian cells, but their genomic distribution has not yet been fully characterized. Here, by using OxiDIP-Seq, a highly sensitive methodology that uses immuno-precipitation with efficient anti-8-oxodG antibodies combined with high-throughput sequencing, we report the genome-wide distribution of 8-oxodG in human non-tumorigenic epithelial breast cells (MCF10A), and mouse embryonic fibroblasts (MEFs). OxiDIP-Seq revealed sites of 8-oxodG accumulation overlapping with γH2AX ChIP-Seq signals within the gene body of transcribed long genes, particularly at the DNA replication origins contained therein. We propose that the presence of persistent single-stranded DNA, as a consequence of transcription-replication clashes at these sites, determines local vulnerability to DNA oxidation and/or its slow repair. This oxidatively-generated damage, likely in combination with other kinds of lesion, might contribute to the formation of DNA double strand breaks and activation of DNA damage response.

ELIXIR-IT HPC@CINECA: high performance computing resources for the bioinformatics community.

BACKGROUND: The advent of Next Generation Sequencing (NGS) technologies and the concomitant reduction in sequencing costs allows unprecedented high throughput profiling of biological systems in a cost-efficient manner. Modern biological experiments are increasingly becoming both data and computationally intensive and the wealth of publicly available biological data is introducing bioinformatics into the "Big Data" era. For these reasons, the effective application of High Performance Computing (HPC) architectures is becoming progressively more recognized also by bioinformaticians. Here we describe HPC resources provisioning pilot programs dedicated to bioinformaticians, run by the Italian Node of ELIXIR (ELIXIR-IT) in collaboration with CINECA, the main Italian supercomputing center. RESULTS: Starting from April 2016, CINECA and ELIXIR-IT launched the pilot Call "ELIXIR-IT HPC@CINECA", offering streamlined access to HPC resources for bioinformatics. Resources are made available either through web front-ends to dedicated workflows developed at CINECA or by providing direct access to the High Performance Computing systems through a standard command-line interface tailored for bioinformatics data analysis. This allows to offer to the biomedical research community a production scale environment, continuously updated with the latest available versions of publicly available reference datasets and bioinformatic tools. Currently, 63 research projects have gained access to the HPC@CINECA program, for a total handout of ~ 8 Millions of CPU/hours and, for data storage, ~ 100 TB of permanent and ~ 300 TB of temporary space. CONCLUSIONS: Three years after the beginning of the ELIXIR-IT HPC@CINECA program, we can appreciate its impact over the Italian bioinformatics community and draw some considerations. Several Italian researchers who applied to the program have gained access to one of the top-ranking public scientific supercomputing facilities in Europe. Those investigators had the opportunity to sensibly reduce computational turnaround times in their research projects and to process massive amounts of data, pursuing research approaches that would have been otherwise difficult or impossible to undertake. Moreover, by taking advantage of the wealth of documentation and training material provided by CINECA, participants had the opportunity to improve their skills in the usage of HPC systems and be better positioned to apply to similar EU programs of greater scale, such as PRACE. To illustrate the effective usage and impact of the resources awarded by the program - in different research applications - we report five successful use cases, which have already published their findings in peer-reviewed journals.

Blastic plasmacytoid dendritic cell neoplasm: genomics mark epigenetic dysregulation as a primary therapeutic target.

Blastic plasmacytoid dendritic cell neoplasm (BPDCN) is a rare and aggressive hematologic malignancy for which there is still no effective therapy. In order to identify genetic alterations useful for a new treatment design, we used whole-exome sequencing to analyze 14 BPDCN patients and the patient-derived CAL-1 cell line. The functional enrichment analysis of mutational data reported the epigenetic regulatory program to be the most significantly undermined (P<0.0001). In particular, twenty-five epigenetic modifiers were found mutated (e.g. ASXL1, TET2, SUZ12, ARID1A, PHF2, CHD8); ASXL1 was the most frequently affected (28.6% of cases). To evaluate the impact of the identified epigenetic mutations at the gene-expression and Histone H3 lysine 27 trimethylation/acetylation levels, we performed additional RNA and pathology tissue-chromatin immunoprecipitation sequencing experiments. The patients displayed enrichment in gene signatures regulated by methylation and modifiable by decitabine administration, shared common H3K27-acetylated regions, and had a set of cell-cycle genes aberrantly up-regulated and marked by promoter acetylation. Collectively, the integration of sequencing data showed the potential of a therapy based on epigenetic agents. Through the adoption of a preclinical BPDCN mouse model, established by CAL-1 cell line xenografting, we demonstrated the efficacy of the combination of the epigenetic drugs 5'-azacytidine and decitabine in controlling disease progression in vivo.

The histone LSD1 demethylase in stemness and cancer transcription programs.

DNA and histone chromatin modifying enzymes play a crucial role in chromatin remodeling in several biological processes. Lysine-specific demethylase 1 (LSD1), the first identified histone demethylase, is a relevant player in the regulation of a broad spectrum of biological processes including development, cellular differentiation, embryonic pluripotency and cancer. Here, we review recent insights on the role of LSD1 activity in chromatin regulatory complexes, its functional role in the epigenetic changes during embryonic development, in the establishment and maintenance of stemness and during cancer progression.

Cu-Zn superoxide dismutase activates muscarinic acetylcholine M1 receptor pathway in neuroblastoma cells.

Muscarinic receptors (mAChRs) control several neuronal functions and are widely expressed in the central nervous system (CNS): M1 subtype represents the predominant mAChR in the CNS. Previously, we showed that antioxidant enzyme Cu-Zn superoxide dismutase (SOD1) is secreted by many cellular lines and specifically interacts with cell surface membrane of human neuroblastoma SK-N-BE cells thus activating phospholipase C (PLC) transduction pathway and increasing intracellular calcium concentration ([Ca(2+)](i)). In addition, we demonstrated that a small amount of SOD1 is contained in large core dense vesicles and that it is secreted in response to depolarization induced by elevated extracellular K(+) concentration. In the present study, we investigated the involvement of muscarinic M1 receptors in SOD1-induced activation of PLC transduction pathway. We showed that, in SK-N-BE cells, SOD1 was able to activate muscarinic M1 receptor producing a phosphorylation of ERK 1/2 and Akt in dose- and time-dependent manner. Interestingly, in the presence of the M1 antagonist pirenzepine, ERK 1/2 and Akt phosphorylation induced by SOD1 was remarkably prevented. This effect was mimicked by knocking-down M1 receptor using two sequences of RNA silencing (siRNA). At functional level, siRNAs against M1 receptor were able to prevent the increase in [Ca(2+)](i) induced by SOD1. The same inhibitory effect on [Ca(2+)](i) changes was produced by the M1 antagonist pirenzepine. Collectively, the results of this study demonstrated that SOD1 could activate a transductional pathway through the involvement of M1 muscarinic receptor.

Downregulation of praja2 restrains endocytosis and boosts tyrosine kinase receptors in kidney cancer.

Clear cell renal cell carcinoma (ccRCC) is the most common kidney cancer in the adult population. Late diagnosis, resistance to therapeutics and recurrence of metastatic lesions account for the highest mortality rate among kidney cancer patients. Identifying novel biomarkers for early cancer detection and elucidating the mechanisms underlying ccRCC will provide clues to treat this aggressive malignant tumor. Here, we report that the ubiquitin ligase praja2 forms a complex with-and ubiquitylates the AP2 adapter complex, contributing to receptor endocytosis and clearance. In human RCC tissues and cells, downregulation of praja2 by oncogenic miRNAs (oncomiRs) and the proteasome markedly impairs endocytosis and clearance of the epidermal growth factor receptor (EGFR), and amplifies downstream mitogenic and proliferative signaling. Restoring praja2 levels in RCC cells downregulates EGFR, rewires cancer cell metabolism and ultimately inhibits tumor cell growth and metastasis. Accordingly, genetic ablation of praja2 in mice upregulates RTKs (i.e. EGFR and VEGFR) and induces epithelial and vascular alterations in the kidney tissue.In summary, our findings identify a regulatory loop between oncomiRs and the ubiquitin proteasome system that finely controls RTKs endocytosis and clearance, positively impacting mitogenic signaling and kidney cancer growth.

Sequence-specific double strand breaks trigger P-TEFb-dependent Rpb1-CTD hyperphosphorylation.

Double strand DNA breaks (DSBs) are one of the most challenging forms of DNA damage which, if left unrepaired, can trigger cellular death and can contribute to cancer. A number of studies have been focused on DNA-damage response (DDR) mechanisms, and most of them rely on the induction of DSBs triggered by chemical compounds or radiations. However, genotoxic drugs and radiation treatments of cultured cell lines induce random DSBs throughout the genome, thus heterogeneously across the cell population, leading to variability of the cellular response. To overcome this aspect, we used here a recently described cell-based DSBs system whereby, upon induction of an inducible restriction enzyme, hundreds of site-specific DSBs are generated across the genome. We show here that sequence-specific DSBs are sufficient to activate the positive transcription elongation factor b (P-TEFb), to trigger hyperphosphorylation of the largest RNA polymerase II carboxyl-terminal-domain (Rpb1-CTD) and to induce activation of p53-transcriptional axis resulting in cell cycle arrest.

Myc and PI3K/AKT signaling cooperatively repress FOXO3a-dependent PUMA and GADD45a gene expression.

Growth factor withdrawal inhibits cell cycle progression by stimulating expression of growth-arresting genes through the activation of Forkhead box O transcription factors such as FOXO3a, which binds to the FHRE-responsive elements of a number of target genes such as PUMA and GADD45a. Following exposure of cells to growth factors FOXO3a-mediated transcription is rapidly repressed. We determined that repression correlates with activation of PI3K/AKT pathway leading to FOXO3a phosphorylation and release of FOXO3a protein from PUMA and GADD45a chromatin. We show here that Myc significantly and selectively contributes to repression of FOXO-mediated expression of PUMA and GADD45a. We found that in Myc deprived cells inhibition of PUMA and GADD45a following serum stimulation is impaired and that Myc does not interfere with p53 induction of PUMA transcription. We observed that following activation, Myc is rapidly recruited to PUMA and GADD45a chromatin, with a concomitant switch in promoter occupancy from FOXO3a to Myc. Myc recruitment stimulates deacetylation of Histone H3 and H4 and methylation of lysine 9 in H3 (H3K9me2) on both PUMA and GADD45 chromatin. These data highlight a Myc role on cell growth by selectively inhibiting FOXO3a induced transcription of PUMA and GADD45.

Epigenetics and immune cells in medulloblastoma.

Medulloblastoma (MB) is a highly malignant childhood tumor of the cerebellum. Transcriptional and epigenetic signatures have classified MB into four molecular subgroups, further stratified into biologically different subtypes with distinct somatic copy-number aberrations, driver genes, epigenetic alterations, activated pathways, and clinical outcomes. The brain tumor microenvironment (BTME) is of importance to regulate a complex network of cells, including immune cells, involved in cancer progression in brain malignancies. MB was considered with a "cold" immunophenotype due to the low influx of immune cells across the blood brain barrier (BBB). Recently, this assumption has been reconsidered because of the identification of infiltrating immune cells showing immunosuppressive phenotypes in the BTME of MB tumors. Here, we are providing a comprehensive overview of the current status of epigenetics alterations occurring during cancer progression with a description of the genomic landscape of MB by focusing on immune cells within the BTME. We further describe how new immunotherapeutic approaches could influence concurring epigenetic mechanisms of the immunosuppressive cells in BTME. In conclusion, the modulation of these molecular genetic complexes in BTME during cancer progression might enhance the therapeutic benefit, thus firing new weapons to fight MB.

Accumulation of 8-oxodG within the human mitochondrial genome positively associates with transcription.

Mitochondrial DNA (mtDNA) can be subject to internal and environmental stressors that lead to oxidatively generated damage and the formation of 8-oxo-7,8-dihydro-2'-deoxyguanine (8-oxodG). The accumulation of 8-oxodG has been linked to degenerative diseases and aging, as well as cancer. Despite the well-described implications of 8-oxodG in mtDNA for mitochondrial function, there have been no reports of mapping of 8-oxodG across the mitochondrial genome. To address this, we used OxiDIP-Seq and mapped 8-oxodG levels in the mitochondrial genome of human MCF10A cells. Our findings indicated that, under steady-state conditions, 8-oxodG is non-uniformly distributed along the mitochondrial genome, and that the longer non-coding region appeared to be more protected from 8-oxodG accumulation compared with the coding region. However, when the cells have been exposed to oxidative stress, 8-oxodG preferentially accumulated in the coding region which is highly transcribed as H1 transcript. Our data suggest that 8-oxodG accumulation in the mitochondrial genome is positively associated with mitochondrial transcription.

OxiDIP-Seq for Genome-wide Mapping of Damaged DNA Containing 8-Oxo-2'-Deoxyguanosine.

8-oxo-7,8-dihydro-2'-deoxyguanosine (8-oxodG) is considered to be a premutagenic DNA lesion generated by 2'-deoxyguanosine (dG) oxidation due to reactive oxygen species (ROS). In recent years, the 8-oxodG distribution in human, mouse, and yeast genomes has been underlined using various next-generation sequencing (NGS)-based strategies. The present study reports the OxiDIP-Seq protocol, which combines specific 8-oxodG immuno-precipitation of single-stranded DNA with NGS, and the pipeline analysis that allows the genome-wide 8-oxodG distribution in mammalian cells. The development of this OxiDIP-Seq method increases knowledge on the oxidative DNA damage/repair field, providing a high-resolution map of 8-oxodG in human cells.

A novel workflow for the qualitative analysis of DNA methylation data.

DNA methylation is an epigenetic modification that plays a pivotal role in major biological mechanisms, such as gene regulation, genomic imprinting, and genome stability. Different combinations of methylated cytosines for a given DNA locus generate different epialleles and alterations of these latter have been associated with several pathological conditions. Existing computational methods and statistical tests relevant to DNA methylation analysis are mostly based on the comparison of average CpG sites methylation levels and they often neglect non-CG methylation. Here, we present EpiStatProfiler, an R package that allows the analysis of CpG and non-CpG based epialleles starting from bisulfite sequencing data through a collection of dedicated extraction functions and statistical tests. EpiStatProfiler is provided with a set of useful auxiliary features, such as customizable genomic ranges, strand-specific epialleles analysis, locus annotation and gene set enrichment analysis. We showcase the package functionalities on two public datasets by identifying putative relevant loci in mice harboring the Huntington's disease-causing Htt gene mutation and in Ctcf +/- mice compared to their wild-type counterparts. To our knowledge, EpiStatProfiler is the first package providing functionalities dedicated to the analysis of epialleles composition derived from any kind of bisulfite sequencing experiment.

Identification of Cdk8 and Cdkn2d as New Prame-Target Genes in 2C-like Embryonic Stem Cells.

Embryonic stem cells (ESCs) present a characteristic pluripotency heterogeneity correspondent to specific metastates. We recently demonstrated that retinoic acid (RA) induces an increase in a specific 2C-like metastate marked by target genes specific to the two-cell embryo stage in preimplantation. Prame (Preferentially expressed antigen in melanoma) is one of the principal actors of the pluripotency stage with a specific role in RA responsiveness. Additionally, PRAME is overexpressed in a variety of cancers, but its molecular functions are poorly understood. To further investigate Prame's downstream targets, we used a chromatin immunoprecipitation sequencing (ChIP-seq) assay in RA-enriched 2C-like metastates and identified two specific target genes, Cdk8 and Cdkn2d, bound by Prame. These two targets, involved in cancer dedifferentiation and pluripotency, have been further validated in RA-resistant ESCs. Here, we observed for the first time that Prame controls the Cdk8 and Cdkn2d genes in ESCs after RA treatment, shedding light on the regulatory network behind the establishment of naïve pluripotency.

MC profiling: a novel approach to analyze DNA methylation heterogeneity in genome-wide bisulfite sequencing data.

DNA methylation is an epigenetic mark implicated in crucial biological processes. Most of the knowledge about DNA methylation is based on bulk experiments, in which DNA methylation of genomic regions is reported as average methylation. However, average methylation does not inform on how methylated cytosines are distributed in each single DNA molecule. Here, we propose Methylation Class (MC) profiling as a genome-wide approach to the study of DNA methylation heterogeneity from bulk bisulfite sequencing experiments. The proposed approach is built on the concept of MCs, groups of DNA molecules sharing the same number of methylated cytosines. The relative abundances of MCs from sequencing reads incorporates the information on the average methylation, and directly informs on the methylation level of each molecule. By applying our approach to publicly available bisulfite-sequencing datasets, we individuated cell-to-cell differences as the prevalent contributor to methylation heterogeneity. Moreover, we individuated signatures of loci undergoing imprinting and X-inactivation, and highlighted differences between the two processes. When applying MC profiling to compare different conditions, we identified methylation changes occurring in regions with almost constant average methylation. Altogether, our results indicate that MC profiling can provide useful insights on the epigenetic status and its evolution at multiple genomic regions.