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.

Studies of OPA1 pathogenic mechanisms in Dominant Optic Atrophy and novel protein function in mitochondrial DNA stability

The mitochondrion is an essential cytoplasmic organelle that provides most of the energy necessary for eukaryotic cell physiology. Mitochondrial structure and functions are maintained by proteins of both mitochondrial and nuclear origin. These organelles are organized in an extended network that dynamically fuses and divides. Mitochondrial morphology results from the equilibrium between fusion and fission processes, controlled by a family of “mitochondria-shaping” proteins. It is becoming clear that defects in mitochondrial dynamics can impair mitochondrial respiration, morphology and motility, leading to apoptotic cell death in vitro and more or less severe neurodegenerative disorders in vivo in humans. Mutations in OPA1, a nuclear encoded mitochondrial protein, cause autosomal Dominant Optic Atrophy (DOA), a heterogeneous blinding disease characterized by retinal ganglion cell degeneration leading to optic neuropathy (Delettre et al., 2000; Alexander et al., 2000). OPA1 is a mitochondrial dynamin-related guanosine triphosphatase (GTPase) protein involved in mitochondrial network dynamics, cytochrome c storage and apoptosis. This protein is anchored or associated on the inner mitochondrial membrane facing the intermembrane space. Eight OPA1 isoforms resulting from alternative splicing combinations of exon 4, 4b and 5b have been described (Delettre et al., 2001). These variants greatly vary among diverse organs and the presence of specific isoforms has been associated with various mitochondrial functions. The different spliced exons encode domains included in the amino-terminal region and contribute to determine OPA1 functions (Olichon et al., 2006). It has been shown that exon 4, that is conserved throughout evolution, confers functions to OPA1 involved in maintenance of the mitochondrial membrane potential and in the fusion of the network. Conversely, exon 4b and exon 5b, which are vertebrate specific, are involved in regulation of cytochrome c release from mitochondria, and activation of apoptosis, a process restricted to vertebrates (Olichon et al., 2007). While Mgm1p has been identified thanks to its role in mtDNA maintenance, it is only recently that OPA1 has been linked to mtDNA stability. Missense mutations in OPA1 cause accumulation of multiple deletions in skeletal muscle. The syndrome associated to these mutations (DOA-1 plus) is complex, consisting of a combination of dominant optic atrophy, progressive external ophtalmoplegia, peripheral neuropathy, ataxia and deafness (Amati- Bonneau et al., 2008; Hudson et al., 2008). OPA1 is the fifth gene associated with mtDNA “breakage syndrome” together with ANT1, PolG1-2 and TYMP (Spinazzola et al., 2009). In this thesis we show for the first time that specific OPA1 isoforms associated to exon 4b are important for mtDNA stability, by anchoring the nucleoids to the inner mitochondrial membrane. Our results clearly demonstrate that OPA1 isoforms including exon 4b are intimately associated to the maintenance of the mitochondrial genome, as their silencing leads to mtDNA depletion. The mechanism leading to mtDNA loss is associated with replication inhibition in cells where exon 4b containing isoforms were down-regulated. Furthermore silencing of exon 4b associated isoforms is responsible for alteration in mtDNA-nucleoids distribution in the mitochondrial network. In this study it was evidenced that OPA1 exon 4b isoform is cleaved to provide a 10kd peptide embedded in the inner membrane by a second transmembrane domain, that seems to be crucial for mitochondrial genome maintenance and does correspond to the second transmembrane domain of the yeasts orthologue encoded by MGM1 or Msp1, which is also mandatory for this process (Diot et al., 2009; Herlan et al., 2003). Furthermore in this thesis we show that the NT-OPA1-exon 4b peptide co-immuno-precipitates with mtDNA and specifically interacts with two major components of the mitochondrial nucleoids: the polymerase gamma and Tfam. Thus, from these experiments the conclusion is that NT-OPA1- exon 4b peptide contributes to the nucleoid anchoring in the inner mitochondrial membrane, a process that is required for the initiation of mtDNA replication and for the distribution of nucleoids along the network. These data provide new crucial insights in understanding the mechanism involved in maintenance of mtDNA integrity, because they clearly demonstrate that, besides genes implicated in mtDNA replications (i.e. polymerase gamma, Tfam, twinkle and genes involved in the nucleotide pool metabolism), OPA1 and mitochondrial membrane dynamics play also an important role. Noticeably, the effect on mtDNA is different depending on the specific OPA1 isoforms down-regulated, suggesting the involvement of two different combined mechanisms. Over two hundred OPA1 mutations, spread throughout the coding region of the gene, have been described to date, including substitutions, deletions or insertions. Some mutations are predicted to generate a truncated protein inducing haploinsufficiency, whereas the missense nucleotide substitutions result in aminoacidic changes which affect conserved positions of the OPA1 protein. So far, the functional consequences of OPA1 mutations in cells from DOA patients are poorly understood. Phosphorus MR spectroscopy in patients with the c.2708delTTAG deletion revealed a defect in oxidative phosphorylation in muscles (Lodi et al., 2004). An energetic impairment has been also show in fibroblasts with the severe OPA1 R445H mutation (Amati-Bonneau et al., 2005). It has been previously reported by our group that OPA1 mutations leading to haploinsufficiency are associated in fibroblasts to an oxidative phosphorylation dysfunction, mainly involving the respiratory complex I (Zanna et al., 2008). In this study we have evaluated the energetic efficiency of a panel of skin fibroblasts derived from DOA patients, five fibroblast cell lines with OPA1 mutations causing haploinsufficiency (DOA-H) and two cell lines bearing mis-sense aminoacidic substitutions (DOA-AA), and compared with control fibroblasts. Although both types of DOA fibroblasts maintained a similar ATP content when incubated in a glucose-free medium, i.e. when forced to utilize the oxidative phosphorylation only to produce ATP, the mitochondrial ATP synthesis through complex I, measured in digitonin-permeabilized cells, was significantly reduced in cells with OPA1 haploinsufficiency only, whereas it was similar to controls in cells with the missense substitutions. Furthermore, evaluation of the mitochondrial membrane potential (DYm) in the two fibroblast lines DOA-AA and in two DOA-H fibroblasts, namely those bearing the c.2819-2A>C mutation and the c.2708delTTAG microdeletion, revealed an anomalous depolarizing response to oligomycin in DOA-H cell lines only. This finding clearly supports the hypothesis that these mutations cause a significant alteration in the respiratory chain function, which can be unmasked only when the operation of the ATP synthase is prevented. Noticeably, oligomycin-induced depolarization in these cells was almost completely prevented by preincubation with cyclosporin A, a well known inhibitor of the permeability transition pore (PTP). This results is very important because it suggests for the first time that the voltage threshold for PTP opening is altered in DOA-H fibroblasts. Although this issue has not yet been addressed in the present study, several are the mechanisms that have been proposed to lead to PTP deregulation, including in particular increased reactive oxygen species production and alteration of Ca2+ homeostasis, whose role in DOA fibroblasts PTP opening is currently under investigation. Identification of the mechanisms leading to altered threshold for PTP regulation will help our understanding of the pathophysiology of DOA, but also provide a strategy for therapeutic intervention.

Investigating the role of single point mutations in the human proteome: a computational study

In the post genomic era with the massive production of biological data the understanding of factors affecting protein stability is one of the most important and challenging tasks for highlighting the role of mutations in relation to human maladies. The problem is at the basis of what is referred to as molecular medicine with the underlying idea that pathologies can be detailed at a molecular level. To this purpose scientific efforts focus on characterising mutations that hamper protein functions and by these affect biological processes at the basis of cell physiology. New techniques have been developed with the aim of detailing single nucleotide polymorphisms (SNPs) at large in all the human chromosomes and by this information in specific databases are exponentially increasing. Eventually mutations that can be found at the DNA level, when occurring in transcribed regions may then lead to mutated proteins and this can be a serious medical problem, largely affecting the phenotype. Bioinformatics tools are urgently needed to cope with the flood of genomic data stored in database and in order to analyse the role of SNPs at the protein level. In principle several experimental and theoretical observations are suggesting that protein stability in the solvent-protein space is responsible of the correct protein functioning. Then mutations that are found disease related during DNA analysis are often assumed to perturb protein stability as well. However so far no extensive analysis at the proteome level has investigated whether this is the case. Also computationally methods have been developed to infer whether a mutation is disease related and independently whether it affects protein stability. Therefore whether the perturbation of protein stability is related to what it is routinely referred to as a disease is still a big question mark. In this work we have tried for the first time to explore the relation among mutations at the protein level and their relevance to diseases with a large-scale computational study of the data from different databases. To this aim in the first part of the thesis for each mutation type we have derived two probabilistic indices (for 141 out of 150 possible SNPs): the perturbing index (Pp), which indicates the probability that a given mutation effects protein stability considering all the “in vitro” thermodynamic data available and the disease index (Pd), which indicates the probability of a mutation to be disease related, given all the mutations that have been clinically associated so far. We find with a robust statistics that the two indexes correlate with the exception of all the mutations that are somatic cancer related. By this each mutation of the 150 can be coded by two values that allow a direct comparison with data base information. Furthermore we also implement computational methods that starting from the protein structure is suited to predict the effect of a mutation on protein stability and find that overpasses a set of other predictors performing the same task. The predictor is based on support vector machines and takes as input protein tertiary structures. We show that the predicted data well correlate with the data from the databases. All our efforts therefore add to the SNP annotation process and more importantly found the relationship among protein stability perturbation and the human variome leading to the diseasome.

High sensitivity analysis of BRAF mutations in neoplastic and non-neoplastic thyroid lesions

The clonal distribution of BRAFV600E in papillary thyroid carcinoma (PTC) has been recently debated. No information is currently available about precursor lesions of PTCs. My first aim was to establish whether the BRAFV600E mutation occurs as a subclonal event in PTCs. My second aim was to screen BRAF mutations in histologically benign tissue of cases with BRAFV600E or BRAFwt PTCs in order to identify putative precursor lesions of PTCs. Highly sensitive semi-quantitative methods were used: Allele Specific LNA quantitative PCR (ASLNAqPCR) and 454 Next-Generation Sequencing (NGS). For the first aim 155 consecutive formalin-fixed and paraffin-embedded (FFPE) specimens of PTCs were analyzed. The percentage of mutated cells obtained was normalized to the estimated number of neoplastic cells. Three groups of tumors were identified: a first had a percentage of BRAF mutated neoplastic cells > 80%; a second group showed a number of BRAF mutated neoplastic cells < 30%; a third group had a distribution of BRAFV600E between 30-80%. The large presence of BRAFV600E mutated neoplastic cell sub-populations suggests that BRAFV600E may be acquired early during tumorigenesis: therefore, BRAFV600E can be heterogeneously distributed in PTC. For the second aim, two groups were studied: one consisted of 20 cases with BRAFV600E mutated PTC, the other of 9 BRAFwt PTCs. Seventy-five and 23 histologically benign FFPE thyroid specimens were analyzed from the BRAFV600E mutated and BRAFwt PTC groups, respectively. The screening of BRAF mutations identified BRAFV600E in “atypical” cell foci from both groups of patients. “Unusual” BRAF substitutions were observed in histologically benign thyroid associated with BRAFV600E PTCs. These mutations were very uncommon in the group with BRAFwt PTCs and in BRAFV600E PTCs. Therefore, lesions carrying BRAF mutations may represent “abortive” attempts at cancer development: only BRAFV600E boosts neoplastic transformation to PTC. BRAFV600E mutated “atypical foci” may represent precursor lesions of BRAFV600E mutated PTCs.

Investigating the role of Copy Number Variants in Specific Language Impairment and identification of new candidate genes

Specific language impairment (SLI) is a complex neurodevelopmental disorder defined as an unexpected failure to develop normal language abilities for no obvious reason. Copy number variants (CNVs) are an important source of variation in the susceptibility to neuropsychiatric disorders. Therefore, a CNV study within SLI families was performed to investigate the role of structural variants in SLI. Among the identified CNVs, we focused on CNVs on chromosome 15q11-q13, recurrently observed in neuropsychiatric conditions, and a homozygous exonic microdeletion in ZNF277. Since this microdeletion falls within the AUTS1 locus, a region linked to autism spectrum disorders (ASD), we investigated a potential role of ZNF277 in SLI and ASD. Frequency data and expression analysis of the ZNF277 microdeletion suggested that this variant may contribute to the risk of language impairments in a complex manner, that is independent of the autism risk previously described in this region. Moreover, we identified an affected individual with a dihydropyrimidine dehydrogenase (DPD) deficiency, caused by compound heterozygosity of two deleterious variants in the gene DPYD. Since DPYD represents a good candidate gene for both SLI and ASD, we investigated its involvement in the susceptibility to these two disorders, focusing on the splicing variant rs3918290, the most common mutation in the DPD deficiency. We observed a higher frequency of rs3918290 in SLI cases (1.2%), compared to controls (~0.6%), while no difference was observed in a large ASD cohort. DPYD mutation screening in 4 SLI and 7 ASD families carrying the splicing variant identified six known missense changes and a novel variant in the promoter region. These data suggest that the combined effect of the mutations identified in affected individuals may lead to an altered DPD activity and that rare variants in DPYD might contribute to a minority of cases, in conjunction with other genetic or non-genetic factors.

Next-Generation Sequencing-Based Mutations Scanning Strategy of the BCR-ABL Kinase Domain in Patients with PhiladelPhia-Chromosome Positive Leukemias Treated with Tyrosine Kinase Inhibitors

In chronic myeloid leukemia and Philadelphia-positive acute lymphoblastic leukemia patients resistant to tyrosine kinase inhibitors (TKIs), BCR-ABL kinase domain mutation status is an essential component of the therapeutic decision algorithm. The recent development of Ultra-Deep Sequencing approach (UDS) has opened the way to a more accurate characterization of the mutant clones surviving TKIs conjugating assay sensitivity and throughput. We decided to set-up and validated an UDS-based for BCR-ABL KD mutation screening in order to i) resolve qualitatively and quantitatively the complexity and the clonal structure of mutated populations surviving TKIs, ii) study the dynamic of expansion of mutated clones in relation to TKIs therapy, iii) assess whether UDS may allow more sensitive detection of emerging clones, harboring critical 2GTKIs-resistant mutations predicting for an impending relapse, earlier than SS. UDS was performed on a Roche GS Junior instrument, according to an amplicon sequencing design and protocol set up and validated in the framework of the IRON-II (Interlaboratory Robustness of Next-Generation Sequencing) International consortium.Samples from CML and Ph+ ALL patients who had developed resistance to one or multiple TKIs and collected at regular time-points during treatment were selected for this study. Our results indicate the technical feasibility, accuracy and robustness of our UDS-based BCR-ABL KD mutation screening approach. UDS was found to provide a more accurate picture of BCR-ABL KD mutation status, both in terms of presence/absence of mutations and in terms of clonal complexity and showed that BCR-ABL KD mutations detected by SS are only the “tip of iceberg”. In addition UDS may reliably pick 2GTKIs-resistant mutations earlier than SS in a significantly greater proportion of patients.The enhanced sensitivity as well as the possibility to identify low level mutations point the UDS-based approach as an ideal alternative to conventional sequencing for BCR-ABL KD mutation screening in TKIs-resistant Ph+ leukemia patients

OPA1 isoforms and protein domains in the rescue of mitochondrial dysfunctions

Mutations in OPA1 gene have been identified in the majority of patients with Dominant Optic Atrophy (DOA), a blinding disease, and the syndromic form DOA-plus. OPA1 protein is a mitochondrial GTPase involved in various mitochondrial functions, present in humans in eight isoforms, resulting from alternative splicing and proteolytic processing. In this study we have investigated the specific role of each isoform through expression in OPA-/- MEFs, by evaluating their ability to improve the defective mitochondrial phenotypes. All isoforms were able to rescue the energetic efficiency, mitochondrial DNA (mtDNA) content and cristae integrity, but only the presence of both long and short forms could recover the mitochondrial morphology. In order to identify the OPA1 protein domains crucial for its functions, we selected and modified the isoform 1, shown to be one of the most efficient in preserving mitochondrial phenotype, to express three specific OPA1 variants, namely: one with a different N-terminus portion, one unable to generate short form owing to deletion of S1 cleavage site and one with a defective GTPase domain. We demonstrated that the simultaneous presence of the N- and C-terminus of OPA1 was essential for the mtDNA maintenance; a cleavable isoform generating s-forms was necessary to completely rescue the energetic competence and the presence of the C-terminus was sufficient to partially recover the cristae ultrastructure. Lastly, several pathogenic OPA1 mutations were inserted in MEF clones and the biochemical features investigated, to correlate the defective phenotypes with the clinical severity of patients. Our results clearly indicate that this cell model reflects very well the clinical characteristics of the patients, and therefore can be proposed as an useful tool to shed light on the pathomechanism underlying DOA.

Studio sul ruolo dei fattori genetici coinvolti nella Valvola Aortica Bicuspide e/o nell'Aneurisma dell'Aorta Toracica

La Valvola Aortica Bicuspide (BAV) rappresenta la più comune anomalia cardiaca congenita, con un’incidenza dello 0,5%-2% nella popolazione generale. Si caratterizza per la presenza di due cuspidi valvolari anziché tre e comprende diverse forme. La BAV è frequentemente associata agli aneurismi dell’aorta toracica (TAA). La dilatazione dell’aorta espone al rischio di sviluppare le complicanze aortiche acute. Materiali e metodi Sono stati reclutati 20 probandi consecutivi sottoposti a chirurgia della valvola aortica e dell'aorta ascendente presso l'Unità di Chirurgia Cardiaca di Policlinico S.Orsola-Malpighi di TAA associata a BAV. Sono stati esclusi individui con una condizione sindromica predisponente l’aneurisma aortico. Ciascun familiare maggiorenne di primo grado è stato arruolato nello studio. L’analisi di mutazioni dell’intero gene ACTA2 è stata eseguita con la tecnica del “bidirectional direct sequencing”. Nelle forme familiari, l’intera porzione codificante del genoma è stata eseguita usando l’exome sequencing. Risultati Dopo il sequenziamento di tutti i 20 esoni e giunzioni di splicing di ACTA2 nei 20 probandi, non è stata individuata alcuna mutazione. Settantasette familiari di primo grado sono stati arruolati. Sono state identificate cinque forme familiari. In una famiglia è stata trovata una mutazione del gene MYH11 non ritenuta patogenetica. Conclusioni La mancanza di mutazioni, sia nelle forme sporadiche sia in quelle familiari, ci suggerisce che questo gene non è coinvolto nello sviluppo della BAV e TAA e, l’associazione che è stata riportata deve essere considerata occasionale. L’architettura genetica della BAV verosimilmente dovrebbe consistere in svariate differenti varianti genetiche che interagiscono in maniera additiva nel determinare un aumento del rischio.