Mechanism of resistance to tyrosine kinase inhibitors in philadelphia-positive acute lymphblastic leukaemia (all): from genetic alterations to impaired RNA editing

The Ph chromosome is the most frequent cytogenetic aberration associated with adult ALL and it represents the single most significant adverse prognostic marker. Despite imatinib has led to significant improvements in the treatment of patients with Ph+ ALL, in the majority of cases resistance developed quickly and disease progressed. Some mechanisms of resistance have been widely described but the full knowledge of contributing factors, driving both the disease and resistance, remains to be defined. The observation of rapid development of lymphoblastic leukemia in mice expressing altered Ikaros (Ik) isoforms represented the background of this study. Ikaros is a zinc finger transcription factor required for normal hemopoietic differentiation and proliferation, particularly in the lymphoid lineages. By means of alternative splicing, Ikaros encodes several proteins that differ in their abilities to bind to a consensus DNA-binding site. Shorter, DNA nonbinding isoforms exert a dominant negative effect, inhibiting the ability of longer heterodimer partners to bind DNA. The differential expression pattern of Ik isoforms in Ph+ ALL patients was analyzed in order to determine if molecular abnormalities involving the Ik gene could associate with resistance to imatinib and dasatinib. Bone marrow and peripheral blood samples from 46 adult patients (median age 55 yrs, 18-76) with Ph+ ALL at diagnosis and during treatment with imatinib (16 pts) or dasatinib (30 pts) were collected. We set up a fast, high-throughput method based on capillary electrophoresis technology to detect and quantify splice variants. 41% Ph+ ALL patients expressed high levels of the non DNA-binding dominant negative Ik6 isoform lacking critical N-terminal zinc-fingers which display abnormal subcellular compartmentalization pattern. Nuclear extracts from patients expressed Ik6 failed to bind DNA in mobility shift assay using a DNA probe containing an Ikaros-specific DNA binding sequence. In 59% Ph+ ALL patients there was the coexistence in the same PCR sample and at the same time of many splice variants corresponded to Ik1, Ik2, Ik4, Ik4A, Ik5A, Ik6, Ik6 and Ik8 isoforms. In these patients aberrant full-length Ikaros isoforms in Ph+ ALL characterized by a 60-bp insertion immediately downstream of exon 3 and a recurring 30-bp in-frame deletion at the end of exon 7 involving most frequently the Ik2, Ik4 isoforms were also identified. Both the insertion and deletion were due to the selection of alternative splice donor and acceptor sites. The molecular monitoring of minimal residual disease showed for the first time in vivo that the Ik6 expression strongly correlated with the BCR-ABL transcript levels suggesting that this alteration could depend on the Bcr-Abl activity. Patient-derived leukaemia cells expressed dominant-negative Ik6 at diagnosis and at the time of relapse, but never during remission. In order to mechanistically demonstrated whether in vitro the overexpression of Ik6 impairs the response to tyrosine kinase inhibitors (TKIs) and contributes to resistance, an imatinib-sensitive Ik6-negative Ph+ ALL cell line (SUP-B15) was transfected with the complete Ik6 DNA coding sequence. The expression of Ik6 strongly increased proliferation and inhibited apoptosis in TKI sensitive cells establishing a previously unknown link between specific molecular defects that involve the Ikaros gene and the resistance to TKIs in Ph+ ALL patients. Amplification and genomic sequence analysis of the exon splice junction regions showed the presence of 2 single nucleotide polymorphisms (SNPs): rs10251980 [A/G] in the exon2/3 splice junction and of rs10262731 [A/G] in the exon 7/8 splice junction in 50% and 36% of patients, respectively. A variant of the rs11329346 [-/C], in 16% of patients was also found. Other two different single nucleotide substitutions not recognized as SNP were observed. Some mutations were predicted by computational analyses (RESCUE approach) to alter cis-splicing elements. In conclusion, these findings demonstrated that the post-transcriptional regulation of alternative splicing of Ikaros gene is defective in the majority of Ph+ ALL patients treated with TKIs. The overexpression of Ik6 blocking B-cell differentiation could contribute to resistance opening a time frame, during which leukaemia cells acquire secondary transforming events that confer definitive resistance to imatinib and dasatinib.

Depicting the role of CDKN2A/ARF alterations in adult BCR-ABL1-positive acute lymphoblastic leukemia patients: from genomic deletions to prognostic impact

This 9p21 locus, encode for important proteins involved in cell cycle regulation and apoptosis containing the p16/CDKN2A (cyclin-dependent kinase inhibitor 2a) tumor suppressor gene and two other related genes, p14/ARF and p15/CDKN2B. This locus, is a major target of inactivation in the pathogenesis of a number of human tumors, both solid and haematologic, and is a frequent site of loss or deletion also in acute lymphoblastic leukemia (ALL) ranging from 18% to 45% 1. In order to explore, at high resolution, the frequency and size of alterations affecting this locus in adult BCR-ABL1-positive ALL and to investigate their prognostic value, 112 patients (101 de novo and 11 relapse cases) were analyzed by genome-wide single nucleotide polymorphisms arrays and gene candidate deep exon sequencing. Paired diagnosis-relapse samples were further available and analyzed for 19 (19%) cases. CDKN2A/ARF and CDKN2B genomic alterations were identified in 29% and 25% of newly diagnosed patients, respectively. Deletions were monoallelic in 72% of cases and in 43% the minimal overlapping region of the lost area spanned only the CDKN2A/2B gene locus. The analysis at the time of relapse showed an almost significant increase in the detection rate of CDKN2A/ARF loss (47%) compared to diagnosis (p = 0.06). Point mutations within the 9p21 locus were found at very low level with only a non-synonymous substition in the exon 2 of CDKN2A. Finally, correlation with clinical outcome showed that deletions of CDKN2A/B are significantly associated with poor outcome in terms of overall survival (p = 0.0206), disease free-survival (p = 0.0010) and cumulative incidence of relapse (p = 0.0014). The inactivation of 9p21 locus by genomic deletions is a frequent event in BCR-ABL1-positive ALL. Deletions are frequently acquired at the leukemia progression and work as a poor prognostic marker.

Different addictive drug exposure induces selective alterations of the Endogenous Opioid System : modulation of transcription and epigenetic mechanisms

Drug addiction manifests clinically as compulsive drug seeking, and cravings that can persist and recur even after extended periods of abstinence. The fundamental principle that unites addictive drugs is that each one enhances synaptic DA by means that dissociate it from normal behavioral control, so that they act to reinforce their own acquisition. Our attention has focused on the study of phenomena associated with the consumption of alcohol and heroin. Alcohol has long been considered an unspecific pharmacological agent, recent molecular pharmacology studies have shown that acts on different primary targets. Through gene expression studies conducted recently it has been shown that the classical opioid receptors are differently involved in the consumption of ethanol and, furthermore, the system nociceptin / NOP, included in the family of endogenous opioid system, and both appear able to play a key role in the initiation of alcohol use in rodents. What emerges is that manipulation of the opioid system, nociceptin, may be useful in the treatment of addictions and there are several evidences that support the use of this strategy. The linkage between gene expression alterations and epigenetic modulation in PDYN and PNOC promoters following alcohol treatment confirm the possible chromatin remodeling mechanism already proposed for alcoholism. In the second part of present study, we also investigated alterations in signaling molecules directly associated with MAPK pathway in a unique collection of postmortem brains from heroin abusers. The interest was focused on understanding the effects that prolonged exposure of heroin can cause in an individual, over the entire MAPK cascade and consequently on the transcription factor ELK1, which is regulated by this pathway. We have shown that the activation of ERK1/2 resulting in Elk-1 phosphorylation in striatal neurons supporting the hypothesis that prolonged exposure to substance abuse causes a dysregulation of MAPK pathway.

Involvement of the opioid system and BDNF in neuroplastic alterations occurring in neuropathic pain conditions

Chronic pain affects one in five adults, reducing quality of life and increasing risk of developing co-morbidities such as depression. Neuropathic pain results by lesions to the nervous system that alter its structure and function leading to spontaneous pain and amplified responses to noxious and innocuous stimuli. The Opioid System is probably the most important system involved in control of nociceptive transmission. Dynorphin and nociceptin systems have been suggested key mediators of some neuropathic pain aspects. An important role also for BDNF has been recently suggested since its involvement in the peripheral and central sensitization phenomena is known. We studied neuroplastic alterations occurring in chronic pain in mice subjected to the chronic constriction injury (CCI). We investigated gene expression alterations of both BDNF and Opioid System at spinal level at different intervals of time. A transient upregulation of pBDNF and pDYN was observed in spinal cord, while increasing upregulation of ppN/OFQ was found in the DRGs of injured mice. Development of neuropathic behavioral signs has been observed in ICR/CD-1 and BDNF+/+ mice, subjected to CCI. A different development of these signs was observed in BDNF+/-. We also studied gene expression changes of investigated systems in different brain areas fourteen days after surgery. We found pBDNF, pDYN, pKOP, ppN/OFQ and pNOP gene expression alterations in several areas of CCI mice. In the same brain regions we also determined bioactive nociceptin peptide levels, and elevated N/OFQ levels were observed in the amygdala area. Histone modifications studies have been performed in BDNF and DYN gene promoters of CCI animal spinal cord showing selected alterations in pDYN gene promoter. In addition, a preliminary characterization of the innovative NOP-EGFP mice was performed. Overall, our results could be useful to understand which and how neuropeptidergic systems are involved in neuroplastic mechanism occurring in neuropathic pain.

Identification of novel genetic alterations in pediatric cytogenetically normal acute myeloid leukemia by next-generation sequencing

Pediatric acute myeloid leukemia (AML) is a molecularly heterogeneous disease that arises from genetic alterations in pathways that regulate self-renewal and myeloid differentiation. While the majority of patients carry recurrent chromosomal translocations, almost 20% of childhood AML do not show any recognizable cytogenetic alteration and are defined as cytogenetically normal (CN)-AML. CN-AML patients have always showed a great variability in response to therapy and overall outcome, underlining the presence of unknown genetic changes, not detectable by conventional analyses, but relevant for pathogenesis, and outcome of AML. The development of novel genome-wide techniques such as next-generation sequencing, have tremendously improved our ability to interrogate the cancer genome. Based on this background, the aim of this research study was to investigate the mutational landscape of pediatric CN-AML patients negative for all the currently known somatic mutations reported in AML through whole-transcriptome sequencing (RNA-seq). RNA-seq performed on diagnostic leukemic blasts from 19 pediatric CN-AML cases revealed a considerable incidence of cryptic chromosomal rearrangements, with the identification of 21 putative fusion genes. Several of the fusion genes that were identified in this study are recurrent and might have a prognostic and/or therapeutic relevance. A paradigm of that is the CBFA2T3-GLIS2 fusion, which has been demonstrated to be a common alteration in pediatric CN-AML, predicting poor outcome. Important findings have been also obtained in the identification of novel therapeutic targets. On one side, the identification of NUP98-JARID1A fusion suggests the use of disulfiram; on the other, here we describe alteration-activating tyrosine kinases, providing functional data supporting the use of tyrosine kinase inhibitors to specifically inhibit leukemia cells. This study provides new insights in the knowledge of genetic alterations underlying pediatric AML, defines novel prognostic markers and putative therapeutic targets, and prospectively ensures a correct risk stratification and risk-adapted therapy also for the “all-neg” AML subgroup.

The biological basis of autism spectrum disorders: evaluation of oxidative stress and erytrocyte membrane alterations

This case-control study involved a total of 29 autistic children (Au) aged 6 to 12 years, and 28 gender and age-matched typically developing children (TD). We evaluated a high number of peripheral oxidative stress parameters, erythrocyte and lymphocyte membrane functional features and membrane lipid composition of erythrocyte. Erythrocyte TBARS, Peroxiredoxin II, Protein Carbonyl Groups and urinary HEL and isoprostane levels were elevated in AU (confirming an imbalance of the redox status of Au); other oxidative stress markers or associated parameters (urinary 8-oxo-dG, plasma Total antioxidant capacity and plasma carbonyl groups, erythrocyte SOD and catalase activities) were unchanged, whilst peroxiredoxin I showed a trend of elevated levels in red blood cells of Au children. A very significant reduction of both erythrocyte and lymphocyte Na+, K+-ATPase activity (NKA), a reduction of erythrocyte membrane fluidity, a reduction of phospatydyl serine exposition on erythrocyte membranes, an alteration in erythrocyte fatty acid membrane profile (increase in MUFA and in ω6/ω3 ratio due to decrease in EPA and DHA) and a reduction of cholesterol content of erythrocyte membrane were found in Au compared to TD, without change in erythrocyte membrane sialic acid content and in lymphocyte membrane fluidity. Some Au clinical features appear to be correlated with these findings; in particular, hyperactivity score appears to be related with some parameters of the lipidomic profile and membrane fluidity, and ADOS and CARS score are inversely related to peroxiredoxin II levels. Oxidative stress and erythrocyte structural and functional alterations may play a role in the pathogenesis of Autism Spectrum Disorders and could be potentially utilized as peripheral biomarkers.

Role of non-coding RNA alterations in melanoma

Cutaneous melanoma (CM) is a potentially lethal form of skin cancer and its most important histopathologic factor for staging is Breslow thickness (BT). Its correct determination is fundamental for pathologists. A deeper understanding of the molecular processes guiding CM pathogenesis could improve diagnosis, treatment and prognosis. MicroRNAs (miRNAs) play a key role in CM biology. The firs aim was to investigate miRNA expression in reference to BT assessment. We found that the combined miRNA expression of miR-21-5p and miR-146a-5p above or below 1.5 was significantly associated with overall survival and successfully identified all superficially spreading melanoma (SSM) patients with relapsing suggesting that the combined assessment of these miRNAs expression could aid in SSM staging. Secondly, we focus on multiple primary melanoma (MPM) patients, which develop multiple primary melanomas in their lifetime, and represent a model of high-risk CM occurrence. We explored the miRNome of single CM and MPM: CM and MPM present several dysregulated miRNAs, including key miRNAs involved in epithelial-mesenchymal transition. A different miRNA profile was observed between 1st and 2nd melanoma from the same patient. MiRNA target analysis revealed a more differentiated and less invasive status of MPMs compared to CMs. This characterization of the miRNA regulatory network of MPMs highlights molecular features differentiating this subtype from CM. Recently, NGS experiments revealed the existence of miRNA variants (isomiRs) with different length and sequence. We identified a shorter 3’isoform as tenfold over-represented compared to the canonical form of miR-125a-5p. Target analysis revealed that miRNA shortening could change the pattern of target gene regulation. Finally, we study miRNA and isomiR dysregulation in benign nevi (BN) and CM and in CM and melanoma metastasis. The reported non-random dysregulation of specific isomiRs contributes to the understanding of the complex melanoma pathogenesis and serves as the basis for further functional studies.

Clinical impact of next-generation sequencing multi-gene panel highlighting the landscape of germline alterations in ovarian cancer patients

BRCA1 and BRCA2 are the most frequently mutated genes in ovarian cancer (OC), crucial both for the identification of cancer predisposition and therapeutic choices. However, germline variants in other genes could be involved in OC susceptibility. We characterized OC patients to detect mutations in genes other than BRCA1/2 that could be associated with a high risk to develop OC, and that could permit patients to enter the most appropriate treatment and surveillance program. Next-Generation Sequencing analysis with a 94-gene panel was performed on germline DNA of 219 OC patients. We identified 34 pathogenic/likely-pathogenic variants in BRCA1/2 and 38 in other 21 genes. Patients with pathogenic/likely-pathogenic variants in non-BRCA1/2 genes developed mainly OC alone compared to the other groups that developed also breast cancer or other tumors (p=0.001). Clinical correlation analysis showed that low-risk patients were significantly associated with platinum sensitivity (p<0.001). Regarding PARP inhibitors (PARPi) response, patients with pathogenic mutations in non-BRCA1/2 genes had significantly worse PFS and OS. Moreover, a statistically significant worse PFS was found for every increase of one thousand platelets before PARPi treatment. To conclude, knowledge about molecular alterations in genes beyond BRCA1/2 in OC could allow for more personalized diagnostic, predictive, prognostic, and therapeutic strategies for OC patients.

TiMMing: developing an innovative suite of bioinformatic tools to harmonize and track the origin of copy number alterations in the evolutive history of multiple myeloma

Multiple Myeloma (MM) is a hematologic cancer with heterogeneous and complex genomic landscape, where Copy Number Alterations (CNAs) play a key role in the disease's pathogenesis and prognosis. It is of biological and clinical interest to study the temporal occurrence of early alterations, as they play a disease "driver" function by deregulating key tumor pathways. This study presents an innovative bioinformatic tools suite created for harmonizing and tracing the origin of CNAs throughout the evolutionary history of MM. To this aim, large cohorts of newly-diagnosed MM (NDMM, N=1582) and Smoldering-MM (SMM, N=282) were aggregated. The tools developed in this study enable the harmonization of CNAs as obtained from different genomic platforms in such a way that a high statistical power can be obtained. By doing so, the high numerosity of those cohorts was harnessed for the identification of novel genes characterized as "driver" (NFKB2, NOTCH2, MAX, EVI5 and MYC-ME2-enhancer), and the generation of an innovative timing model, implemented with a statistical method to introduce confidence intervals in the CNAs-calls. By applying this model on both NDMM and SMM cohorts, it was possible to identify specific CNAs (1q(CKS1B)amp, 13q(RB1)del, 11q(CCND1)amp and 14q(MAX)del) and categorize them as "early"/ "driver" events. A high level of precision was guaranteed by the narrow confidence intervals in the timing estimates. These CNAs were proposed as critical MM alterations, which play a foundational role in the evolutionary history of both SMM and NDMM. Finally, a multivariate survival model was able to identify the independent genomic alterations with the greatest effect on patients’ survival, including RB1-del, CKS1B-amp, MYC-amp, NOTCH2-amp and TRAF3-del/mut. In conclusion, the alterations that were identified as both "early-drivers” and correlated with patients’ survival were proposed as biomarkers that, if included in wider survival models, could provide a better disease stratification and an improved prognosis definition.

Combining molecular alterations and functional imaging in metastatic castration resistant prostate cancer treated with taxanes

The treatment of metastatic castration-resistant prostate cancer (mCRPC) is currently characterized by several drugs with different mechanisms of action, such as new generation hormonal agents (abiraterone, enzalutamide), chemotherapy (docetaxel, cabazitaxel), PARP inhibitors (olaparib) and radiometabolic therapies (radium-223, LuPSMA). There is an urgent need to identify biomarkers to guide personalized therapy in mCRPC. In recent years, the status of androgen receptor (AR) gene detected in liquid biopsy has been associated with outcomes in patients treated with abiraterone or enzalutamide. More recently, plasma tumor DNA (ptDNA) and its changes during treatment have been identified as early indicators of response to anticancer treatments. Recent works also suggested a potential role of tumor-related metabolic parameters of 18Fluoro-Choline Positron Emission Tomography (F18CH-PET)-computed tomography (CT) as a prognostic tool in mCRCP. Other clinical features, such as the presence of visceral metastases, have been correlated with outcome in mCRPC patients. Recent studies conducted by our research group have designed and validated a prognostic model based on the combination of molecular characteristics (ptDNA levels), metabolic features found in basal FCH PET scans (metabolic tumor volume values, MTV), clinical parameters (absence or presence of visceral metastases), and laboratory tests (serum lactate dehydrogenase levels, LDH). Within this PhD project, 30 patients affected by mCRPC, pre-treated with abiraterone or enzalutamide, candidate for taxane-based treatments (docetaxel or cabazitaxel), have been prospectively evaluated. The prognostic model previously described was applied to this population, to interrogate its prognostic power in a more advanced cohort of patients, resulting in a further external validation of the tool.

Identification of outcome-oriented cut-offs for copy number alterations in multiple myeloma: predictive biomarkers with improved prognostic accuracy

Aim of the present study was to develop a statistical approach to define the best cut-off Copy number alterations (CNAs) calling from genomic data provided by high throughput experiments, able to predict a specific clinical end-point (early relapse, 18 months) in the context of Multiple Myeloma (MM). 743 newly diagnosed MM patients with SNPs array-derived genomic and clinical data were included in the study. CNAs were called both by a conventional (classic, CL) and an outcome-oriented (OO) method, and Progression Free Survival (PFS) hazard ratios of CNAs called by the two approaches were compared. The OO approach successfully identified patients at higher risk of relapse and the univariate survival analysis showed stronger prognostic effects for OO-defined high-risk alterations, as compared to that defined by CL approach, statistically significant for 12 CNAs. Overall, 155/743 patients relapsed within 18 months from the therapy start. A small number of OO-defined CNAs were significantly recurrent in early-relapsed patients (ER-CNAs) - amp1q, amp2p, del2p, del12p, del17p, del19p -. Two groups of patients were identified either carrying or not ≥1 ER-CNAs (249 vs. 494, respectively), the first one with significantly shorter PFS and overall survivals (OS) (PFS HR 2.15, p<0001; OS HR 2.37, p<0.0001). The risk of relapse defined by the presence of ≥1 ER-CNAs was independent from those conferred both by R-IIS 3 (HR=1.51; p=0.01) and by low quality (< stable disease) clinical response (HR=2.59 p=0.004). Notably, the type of induction therapy was not descriptive, suggesting that ER is strongly related to patients’ baseline genomic architecture. In conclusion, the OO- approach employed allowed to define CNAs-specific dynamic clonality cut-offs, improving the CNAs calls’ accuracy to identify MM patients with the highest probability to ER. As being outcome-dependent, the OO-approach is dynamic and might be adjusted according to the selected outcome variable of interest.