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.

Mechanisms Contributing to Tyrosin Kinase Inhibitor Resistance in GISTs: Toward a Personalized Therapy

Gastrointestinal stromal tumors (GISTs) are the most common mesenchymal tumors in the gastrointestinal tract. This work considers the pharmacological response in GIST patients treated with imatinib by two different angles: the genetic and somatic point of view. We analyzed polymorphisms influence on treatment outcome, keeping in consideration SNPs in genes involved in drug transport and folate pathway. Naturally, all these intriguing results cannot be considered as the only main mechanism in imatinib response. GIST mainly depends by oncogenic gain of function mutations in tyrosin kinase receptor genes, KIT or PDGFRA, and the mutational status of these two genes or acquisition of secondary mutation is considered the main player in GIST development and progression. To this purpose we analyzed the secondary mutations to better understand how these are involved in imatinib resistance. In our analysis we considered both imatinib and the second line treatment, sunitinib, in a subset of progressive patients. KIT/PDGFRA mutation analysis is an important tool for physicians, as specific mutations may guide therapeutic choices. Currently, the only adaptations in treatment strategy include imatinib starting dose of 800 mg/daily in KIT exon-9-mutated GISTs. In the attempt to individualize treatment, genetic polymorphisms represent a novelty in the definition of biomarkers of imatinib response in addition to the use of tumor genotype. Accumulating data indicate a contributing role of pharmacokinetics in imatinib efficacy, as well as initial response, time to progression and acquired resistance. At the same time it is becoming evident that genetic host factors may contribute to the observed pharmacokinetic inter-patient variability. Genetic polymorphisms in transporters and metabolism may affect the activity or stability of the encoded enzymes. Thus, integrating pharmacogenetic data of imatinib transporters and metabolizing genes, whose interplay has yet to be fully unraveled, has the potential to provide further insight into imatinib response/resistance mechanisms.

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.

Computational analyses on the structure-function relation in ion channels

Ion channels are protein molecules, embedded in the lipid bilayer of the cell membranes. They act as powerful sensing elements switching chemicalphysical stimuli into ion-fluxes. At a glance, ion channels are water-filled pores, which can open and close in response to different stimuli (gating), and one once open select the permeating ion species (selectivity). They play a crucial role in several physiological functions, like nerve transmission, muscular contraction, and secretion. Besides, ion channels can be used in technological applications for different purpose (sensing of organic molecules, DNA sequencing). As a result, there is remarkable interest in understanding the molecular determinants of the channel functioning. Nowadays, both the functional and the structural characteristics of ion channels can be experimentally solved. The purpose of this thesis was to investigate the structure-function relation in ion channels, by computational techniques. Most of the analyses focused on the mechanisms of ion conduction, and the numerical methodologies to compute the channel conductance. The standard techniques for atomistic simulation of complex molecular systems (Molecular Dynamics) cannot be routinely used to calculate ion fluxes in membrane channels, because of the high computational resources needed. The main step forward of the PhD research activity was the development of a computational algorithm for the calculation of ion fluxes in protein channels. The algorithm - based on the electrodiffusion theory - is computational inexpensive, and was used for an extensive analysis on the molecular determinants of the channel conductance. The first record of ion-fluxes through a single protein channel dates back to 1976, and since then measuring the single channel conductance has become a standard experimental procedure. Chapter 1 introduces ion channels, and the experimental techniques used to measure the channel currents. The abundance of functional data (channel currents) does not match with an equal abundance of structural data. The bacterial potassium channel KcsA was the first selective ion channels to be experimentally solved (1998), and after KcsA the structures of four different potassium channels were revealed. These experimental data inspired a new era in ion channel modeling. Once the atomic structures of channels are known, it is possible to define mathematical models based on physical descriptions of the molecular systems. These physically based models can provide an atomic description of ion channel functioning, and predict the effect of structural changes. Chapter 2 introduces the computation methods used throughout the thesis to model ion channels functioning at the atomic level. In Chapter 3 and Chapter 4 the ion conduction through potassium channels is analyzed, by an approach based on the Poisson-Nernst-Planck electrodiffusion theory. In the electrodiffusion theory ion conduction is modeled by the drift-diffusion equations, thus describing the ion distributions by continuum functions. The numerical solver of the Poisson- Nernst-Planck equations was tested in the KcsA potassium channel (Chapter 3), and then used to analyze how the atomic structure of the intracellular vestibule of potassium channels affects the conductance (Chapter 4). As a major result, a correlation between the channel conductance and the potassium concentration in the intracellular vestibule emerged. The atomic structure of the channel modulates the potassium concentration in the vestibule, thus its conductance. This mechanism explains the phenotype of the BK potassium channels, a sub-family of potassium channels with high single channel conductance. The functional role of the intracellular vestibule is also the subject of Chapter 5, where the affinity of the potassium channels hEag1 (involved in tumour-cell proliferation) and hErg (important in the cardiac cycle) for several pharmaceutical drugs was compared. Both experimental measurements and molecular modeling were used in order to identify differences in the blocking mechanism of the two channels, which could be exploited in the synthesis of selective blockers. The experimental data pointed out the different role of residue mutations in the blockage of hEag1 and hErg, and the molecular modeling provided a possible explanation based on different binding sites in the intracellular vestibule. Modeling ion channels at the molecular levels relates the functioning of a channel to its atomic structure (Chapters 3-5), and can also be useful to predict the structure of ion channels (Chapter 6-7). In Chapter 6 the structure of the KcsA potassium channel depleted from potassium ions is analyzed by molecular dynamics simulations. Recently, a surprisingly high osmotic permeability of the KcsA channel was experimentally measured. All the available crystallographic structure of KcsA refers to a channel occupied by potassium ions. To conduct water molecules potassium ions must be expelled from KcsA. The structure of the potassium-depleted KcsA channel and the mechanism of water permeation are still unknown, and have been investigated by numerical simulations. Molecular dynamics of KcsA identified a possible atomic structure of the potassium-depleted KcsA channel, and a mechanism for water permeation. The depletion from potassium ions is an extreme situation for potassium channels, unlikely in physiological conditions. However, the simulation of such an extreme condition could help to identify the structural conformations, so the functional states, accessible to potassium ion channels. The last chapter of the thesis deals with the atomic structure of the !- Hemolysin channel. !-Hemolysin is the major determinant of the Staphylococcus Aureus toxicity, and is also the prototype channel for a possible usage in technological applications. The atomic structure of !- Hemolysin was revealed by X-Ray crystallography, but several experimental evidences suggest the presence of an alternative atomic structure. This alternative structure was predicted, combining experimental measurements of single channel currents and numerical simulations. This thesis is organized in two parts, in the first part an overview on ion channels and on the numerical methods adopted throughout the thesis is provided, while the second part describes the research projects tackled in the course of the PhD programme. The aim of the research activity was to relate the functional characteristics of ion channels to their atomic structure. In presenting the different research projects, the role of numerical simulations to analyze the structure-function relation in ion channels is highlighted.

Stem cell pathways in the basal-like breast carcinoma: the role of beta-catenin and HIF-1alpha

Basal-like tumor is an aggressive breast carcinoma subtype that displays an expression signature similar to that of the basal/myoepithelial cells of the breast tissue. Basal-like carcinoma are characterized by over-expression of the Epidermal Growth Factor receptor (EGFR), high frequency of p53 mutations, cytoplasmic/nuclear localization of beta-catenin, overexpression of the Hypoxia inducible factor (HIF)-1alpha target Carbonic Anhydrase isoenzime 9 (CA9) and a gene expression pattern similar to that of normal and cancer stem cells, including the over-expression of the mammary stem cell markers CD44. In this study we investigated the role of p53, EGFR, beta-catenin and HIF-1alpha in the regulation of stem cell features and genes associated with the basal-like gene expression profile. The findings reported in this investigation indicate that p53 inactivation in ductal breast carcinoma cells leads to increased EGFR mRNA and protein levels. In our experimental model, EGFR overexpression induces beta-catenin cytoplasmatic stabilization and transcriptional activity and, by that, leads to increased aggressive features including mammosphere (MS) forming and growth capacity, invasive potential and overexpression of the mammary stem cell gene CD44. Moreover we found that EGFR/beta-catenin axis promotes hypoxia survival in breast carcinoma cells via increased CA9 expression. Indeed beta-catenin positively regulates CA9 expression upon hypoxia exposure. Interestingly we found that beta-catenin inhibits HIF-1alpha transcriptional activity. Looking for the mechanism, we found that CA9 expression is promoted by HIF-1alpha and cytoplasmatic beta-catenin further increased it post-transcriptionally, via direct mRNA binding and stabilization. These data reveal a functional beta-catenin/HIF-1alpha interplay among hallmarks of basal-like tumors and unveil a new functional role for cytoplasmic beta-catenin in the phenotype of such tumors. Therefore it can be proposed that the interplay here described among EGFR/beta-catenin and HIF-1alpha may play a role in breast cancer stem cell survival and function.

Genetic analysis of tumour initiation and progression in a Drosophila model of epithelial cancer

Neoplastic overgrowth depends on the cooperation of several mutations ultimately leading to major rearrangements in cellular behaviour. The molecular crosstalk occurring between precancerous and normal cells strongly influences the early steps of the tumourigenic process as well as later stages of the disease. Precancerous cells are often removed by cell death from normal tissues but the mechanisms responsible for such fundamental safeguard processes remain in part elusive. To gain insight into these phenomena I took advantage of the clonal analysis methods available in Drosophila for studying the phenotypes due to loss of function of the neoplastic tumour suppressor lethal giant larvae (lgl). I found that lgl mutant cells growing in wild-type imaginal wing discs are subject to the phenomenon of cell competition and are eliminated by JNK-dependent cell death because they express very low levels of dMyc oncoprotein compared to those in the surrounding tissue. Indeed, in non-competitive backgrounds lgl mutant clones are able to overgrow and upregulate dMyc, overwhelming the neighbouring tissue and forming tumourous masses that display several cancer hallmarks. These phenotypes are completely abolished by reducing dMyc abundance within mutant cells while increasing it in lgl clones growing in a competitive context re-establishes their tumourigenic potential. Similarly, the neoplastic growth observed upon the oncogenic cooperation between lgl mutation and activated Ras/Raf/MAPK signalling was found to be characterised by and dependent on the ability of cancerous cells to upregulate dMyc with respect to the adjacent normal tissue, through both transcriptional and post-transcriptional mechanisms, thereby confirming its key role in lgl-induced tumourigenesis. These results provide first evidence that the dMyc oncoprotein is required in lgl mutant tissue to promote invasive overgrowth in developing and adult epithelial tissues and that dMyc abundance inside versus outside lgl mutant clones plays a key role in driving neoplastic overgrowth.

Computational methods for genome screening

Motivation An actual issue of great interest, both under a theoretical and an applicative perspective, is the analysis of biological sequences for disclosing the information that they encode. The development of new technologies for genome sequencing in the last years, opened new fundamental problems since huge amounts of biological data still deserve an interpretation. Indeed, the sequencing is only the first step of the genome annotation process that consists in the assignment of biological information to each sequence. Hence given the large amount of available data, in silico methods became useful and necessary in order to extract relevant information from sequences. The availability of data from Genome Projects gave rise to new strategies for tackling the basic problems of computational biology such as the determination of the tridimensional structures of proteins, their biological function and their reciprocal interactions. Results The aim of this work has been the implementation of predictive methods that allow the extraction of information on the properties of genomes and proteins starting from the nucleotide and aminoacidic sequences, by taking advantage of the information provided by the comparison of the genome sequences from different species. In the first part of the work a comprehensive large scale genome comparison of 599 organisms is described. 2,6 million of sequences coming from 551 prokaryotic and 48 eukaryotic genomes were aligned and clustered on the basis of their sequence identity. This procedure led to the identification of classes of proteins that are peculiar to the different groups of organisms. Moreover the adopted similarity threshold produced clusters that are homogeneous on the structural point of view and that can be used for structural annotation of uncharacterized sequences. The second part of the work focuses on the characterization of thermostable proteins and on the development of tools able to predict the thermostability of a protein starting from its sequence. By means of Principal Component Analysis the codon composition of a non redundant database comprising 116 prokaryotic genomes has been analyzed and it has been showed that a cross genomic approach can allow the extraction of common determinants of thermostability at the genome level, leading to an overall accuracy in discriminating thermophilic coding sequences equal to 95%. This result outperform those obtained in previous studies. Moreover, we investigated the effect of multiple mutations on protein thermostability. This issue is of great importance in the field of protein engineering, since thermostable proteins are generally more suitable than their mesostable counterparts in technological applications. A Support Vector Machine based method has been trained to predict if a set of mutations can enhance the thermostability of a given protein sequence. The developed predictor achieves 88% accuracy.

Mitochondrial dysfunction in a cell model of thyroid oncocytoma

The role of mitochondrial dysfunction in cancer has long been a subject of great interest. In this study, such dysfunction has been examined with regards to thyroid oncocytoma, a rare form of cancer, accounting for less than 5% of all thyroid cancers. A peculiar characteristic of thyroid oncocytic cells is the presence of an abnormally large number of mitochondria in the cytoplasm. Such mitochondrial hyperplasia has also been observed in cells derived from patients suffering from mitochondrial encephalomyopathies, where mutations in the mitochondrial DNA(mtDNA) encoding the respiratory complexes result in oxidative phosphorylation dysfunction. An increase in the number of mitochondria occurs in the latter in order to compensate for the respiratory deficiency. This fact spurred the investigation into the presence of analogous mutations in thyroid oncocytic cells. In this study, the only available cell model of thyroid oncocytoma was utilised, the XTC-1 cell line, established from an oncocytic thyroid metastasis to the breast. In order to assess the energetic efficiency of these cells, they were incubated in a medium lacking glucose and supplemented instead with galactose. When subjected to such conditions, glycolysis is effectively inhibited and the cells are forced to use the mitochondria for energy production. Cell viability experiments revealed that XTC-1 cells were unable to survive in galactose medium. This was in marked contrast to the TPC-1 control cell line, a thyroid tumour cell line which does not display the oncocytic phenotype. In agreement with these findings, subsequent experiments assessing the levels of cellular ATP over incubation time in galactose medium, showed a drastic and continual decrease in ATP levels only in the XTC-1 cell line. Furthermore, experiments on digitonin-permeabilised cells revealed that the respiratory dysfunction in the latter was due to a defect in complex I of the respiratory chain. Subsequent experiments using cybrids demonstrated that this defect could be attributed to the mitochondrially-encoded subunits of complex I as opposed to the nuclearencoded subunits. Confirmation came with mtDNA sequencing, which detected the presence of a novel mutation in the ND1 subunit of complex I. In addition, a mutation in the cytochrome b subunit of complex III of the respiratory chain was detected. The fact that XTC-1 cells are unable to survive when incubated in galactose medium is consistent with the fact that many cancers are largely dependent on glycolysis for energy production. Indeed, numerous studies have shown that glycolytic inhibitors are able to induce apoptosis in various cancer cell lines. Subsequent experiments were therefore performed in order to identify the mode of XTC-1 cell death when subjected to the metabolic stress imposed by the forced use of the mitochondria for energy production. Cell shrinkage and mitochondrial fragmentation were observed in the dying cells, which would indicate an apoptotic type of cell death. Analysis of additional parameters however revealed a lack of both DNA fragmentation and caspase activation, thus excluding a classical apoptotic type of cell death. Interestingly, cleavage of the actin component of the cytoskeleton was observed, implicating the action of proteases in this mode of cell demise. However, experiments employing protease inhibitors failed to identify the specific protease involved. It has been reported in the literature that overexpression of Bcl-2 is able to rescue cells presenting a respiratory deficiency. As the XTC-1 cell line is not only respiration-deficient but also exhibits a marked decrease in Bcl-2 expression, it is a perfect model with which to study the relationship between Bcl-2 and oxidative phosphorylation in respiratory-deficient cells. Contrary to the reported literature studies on various cell lines harbouring defects in the respiratory chain, Bcl-2 overexpression was not shown to increase cell survival or rescue the energetic dysfunction in XTC-1 cells. Interestingly however, it had a noticeable impact on cell adhesion and morphology. Whereas XTC-1 cells shrank and detached from the growth surface under conditions of metabolic stress, Bcl-2-overexpressing XTC-1 cells appeared much healthier and were up to 45% more adherent. The target of Bcl-2 in this setting appeared to be the actin cytoskeleton, as the cleavage observed in XTC-1 cells expressing only endogenous levels of Bcl-2, was inhibited in Bcl-2-overexpressing cells. Thus, although unable to rescue XTC-1 cells in terms of cell viability, Bcl-2 is somehow able to stabilise the cytoskeleton, resulting in modifications in cell morphology and adhesion. The mitochondrial respiratory deficiency observed in cancer cells is thought not only to cause an increased dependency on glycolysis but it is also thought to blunt cellular responses to anticancer agents. The effects of several therapeutic agents were thus assessed for their death-inducing ability in XTC-1 cells. Cell viability experiments clearly showed that the cells were more resistant to stimuli which generate reactive oxygen species (tert-butylhydroperoxide) and to mitochondrial calcium-mediated apoptotic stimuli (C6-ceramide), as opposed to stimuli inflicting DNA damage (cisplatin) and damage to protein kinases(staurosporine). Various studies in the literature have reported that the peroxisome proliferator-activated receptor-coactivator 1(PGC-1α), which plays a fundamental role in mitochondrial biogenesis, is also involved in protecting cells against apoptosis caused by the former two types of stimuli. In accordance with these observations, real-time PCR experiments showed that XTC-1 cells express higher mRNA levels of this coactivator than do the control cells, implicating its importance in drug resistance. In conclusion, this study has revealed that XTC-1 cells, like many cancer cell lines, are characterised by a reduced energetic efficiency due to mitochondrial dysfunction. Said dysfunction has been attributed to mutations in respiratory genes encoded by the mitochondrial genome. Although the mechanism of cell demise in conditions of metabolic stress is unclear, the potential of targeting thyroid oncocytic cancers using glycolytic inhibitors has been illustrated. In addition, the discovery of mtDNA mutations in XTC-1 cells has enabled the use of this cell line as a model with which to study the relationship between Bcl-2 overexpression and oxidative phosphorylation in cells harbouring mtDNA mutations and also to investigate the significance of such mutations in establishing resistance to apoptotic stimuli.

Mitochondrial dysfunction in hereditary optic neuropathies

MITOCHONDRIAL DYSFUNCTION IN HEREDITARY OPTIC NEUROPATHIES Mitochondrial pathologies are a heterogeneous group of clinical manifestations characterized by oxidative phosphorylation impairment. At the beginning of their recognition mitochondrial pathologies were regarded as rare disorders but indeed they are more frequent than originally thought. Due to the unique mitochondria peculiarities mitochondrial pathologies can be caused by mutations in both mitochondrial and nuclear genomes. The poor knowledge of pathologic mechanism of these disorders has not allowed a real development of the “mitochondrial medicine”, that is currently limited to symptoms mitigation. Leber hereditary optic neuropathy (LHON) was the first pathology to be linked to a point mutation in the mtDNA. The mechanism by which point mutations in mitochondrial gene encoding Complex I subunits leads to optic nerve degeneration is still unknown, although is well accepted that other genetic or environmental factors are involved in the modulation of pathology, where a pivotal role is certainly played by oxidative stress. We studied the relationship between the Ala16Val dimorphism in the mitochondrial targeting sequence of nuclear gene SOD2 and the 3460/ND1 LHON mutation. Our results show that, in control population, the heterozygous SOD2 genotype is associated to a higher activity and quantity of MnSOD, particularly with respect to Val homozygotes. Furthermore, we demonstrated that LHON patients harboring at least one Ala allele are characterized by an increased MnSOD activity with respect to relative control population. Since the ATP synthesis rate – severely reduced in LHON patients lymphocytes - is not affected by the SOD2 genotype, we concluded that SOD2 gene could modulate the pathogenicity of LHON mutations through a mechanism associated to an increase of reactive oxygen species production. Autosomal dominant optic atrophy (ADOA) is a pathology linked to mutations in nuclear gene encoding Opa1, a dynamin-related protein localized in the mitochondrial matrix. Although the clinical course is slightly different, the endpoint of ADOA is exactly the same of LHON: optic nerve degeneration with specific involvement of retinal ganglion cells. Opa1 is a relatively new protein, whose major role is the regulation of mitochondrial fusion. Mitochondrial morphology is the results of the equilibrium between two opposite force: fusion and fission, two processes that have to be finely regulated in order to preserve mitochondrial and cellular physiology. We studied fibroblasts deriving from ADOA patients characterized by a new deletion in the GTPase domain of the OPA1 gene. The biochemical characterization of ADOA and control fibroblasts has concerned the evaluation of ATP synthesis rate, mitochondrial membrane potential in different metabolic conditions and the morphological status of mitochondria. Regarding ATP synthesis rate we did not find significant differences between ADOA and control fibroblasts even though a trend toward increased reduction in ADOA samples is observed when fibroblasts are grown in absence of glucose or in the medium containing gramicidin. Furthermore, we found that also in ADOA fibroblasts membrane potential is actively maintained by proton pumping of fully functional respiratory chain complexes. Our results indicate that the mutation found in the pedigree analyzed acts primary impairing the mitochondrial fusion without affecting the energy production, supporting the notion that cell function is tightly linked to mitochondrial morphology. Mitochondrial dysfunctions are acquiring great attention because of their recognized relevance not only in aging but also in age-related pathologies including cancer, cardiovascular disease, type II diabetes, and neurodegenerative disorders. The involvement of mitochondria in such detrimental pathologies that, currently, have become so common enhances the necessity of standardization of therapeutic strategies capable of rescuing the normal mitochondrial function. In order to propose an alternative treatment for energy deficiency-disorders we tested the effect of substrates capable to stimulate the substrate-level phosphorylation on viability and energy availability in different experimental models grown under different metabolic conditions. In fibroblasts, the energy defect was achieved by culturing cells in presence of oligomycin, an inhibitor of ATP synthase complex. NARP cybrids have been used as model of mitochondrial pathology. Cell viability and ATP content have been considered as parameters to assay the capability of exogenous substrate to rescue energy failure. Our results suggest that patients suffering for some forms of ATP synthase deficiency, or characterized by a deficiency in energy production, might benefit from dietary or pharmacological treatment based on supplementation of α-ketoglutarate and aspartate.

Pathogenic mechanisms in mitochondrial optic neuropathies

Leber’s hereditary optic neuropathy (LHON) and Autosomal Dominant Optic Atrophy (ADOA) are the two most common inherited optic neuropathies and both are the result of mitochondrial dysfunctions. Despite the primary mutations causing these disorders are different, being an mtDNA mutation in subunits of complex I in LHON and defects in the nuclear gene encoding the mitochondrial protein OPA1 in ADOA, both pathologies share some peculiar features, such a variable penetrance and tissue-specificity of the pathological processes. Probably, one of the most interesting and unclear aspect of LHON is the variable penetrance. This phenomenon is common in LHON families, most of them being homoplasmic mutant. Inter-family variability of penetrance may be caused by nuclear or mitochondrial ‘secondary’ genetic determinants or other predisposing triggering factors. We identified a compensatory mechanism in LHON patients, able to distinguish affected individuals from unaffected mutation carriers. In fact, carrier individuals resulted more efficient than affected subjects in increasing the mitochondrial biogenesis to compensate for the energetic defect. Thus, the activation of the mitochondrial biogenesis may be a crucial factor in modulating penetrance, determining the fate of subjects harbouring LHON mutations. Furthermore, mtDNA content can be used as a molecular biomarker which, for the first time, clearly differentiates LHON affected from LHON carrier individuals, providing a valid mechanism that may be exploited for development of therapeutic strategies. Although the mitochondrial biogenesis gained a relevant role in LHON pathogenesis, we failed to identify a genetic modifying factor for the variable penetrance in a set of candidate genes involved in the regulation of this process. A more systematic high-throughput approach will be necessary to select the genetic variants responsible for the different efficiency in activating mitochondrial biogenesis. A genetic modifying factor was instead identified in the MnSOD gene. The SNP Ala16Val in this gene seems to modulate LHON penetrance, since the Ala allele in this position significantly predisposes to be affected. Thus, we propose that high MnSOD activity in mitochondria of LHON subjects may produce an overload of H2O2 for the antioxidant machinery, leading to release from mitochondria of this radical and promoting a severe cell damage and death ADOA is due to mutation in the OPA1 gene in the large majority of cases. The causative nuclear defects in the remaining families with DOA have not been identified yet, but a small number of families have been mapped to other chromosomal loci (OPA3, OPA4, OPA5, OPA7, OPA8). Recently, a form of DOA and premature cataract (ADOAC) has been associated to pathogenic mutations of the OPA3 gene, encoding a mitochondrial protein. In the last year OPA3 has been investigated by two different groups, but a clear function for this protein and the pathogenic mechanism leading to ADOAC are still unclear. Our study on OPA3 provides new information about the pattern of expression of the two isoforms OPA3V1 and OPA3V2, and, moreover, suggests that OPA3 may have a different function in mitochondria from OPA1, the major site for ADOA mutations. In fact, based on our results, we propose that OPA3 is not involved in the mitochondrial fusion process, but, on the contrary, it may regulate mitochondrial fission. Furthermore, at difference from OPA1, we excluded a role for OPA3 in mtDNA maintenance and we failed to identify a direct interaction between OPA3 and OPA1. Considering the results from overexpression and silencing of OPA3, we can conclude that the overexpression has more drastic consequences on the cells than silencing, suggesting that OPA3 may cause optic atrophy via a gain-of-function mechanism. These data provide a new starting point for future investigations aimed at identifying the exact function of OPA3 and the pathogenic mechanism causing ADOAC.

Influence of genetic variants and post transcriptional factors on imatinib transporters as determinants of the pharmacological response

Introduction – Although imatinib (IM) is a recognized gold standard in chronic myeloid leukemia (CML) therapy, resistance has emerged in a significant proportion of patients. Aim – The aim of this study was: (1) to investigate the role of genetic variants in genes encoding for IM transporters, as candidate of IM responsiveness and (2) to test the influence of miRNAs on IM response, focusing on efflux transporters. Methods – As a first step, a panel of polymorphisms (SNPs) was genotyped in a subgroup population of 189 patients enrolled in the Tyrosine Kinase Inhibitor Optimization and Selectivity (TOPS) trial. The association with cytogenetic response and molecular response (MR) was assessed for each SNP. As a second step, an in vitro IM-resistant model (K-562 CML cell line) was established. miRNAs profiles were analyzed using Taqman arrays and in silico search was performed for miRNAs deregulated after IM treatment. mRNA and protein expression were quantified using TaqMan realtime PCR and Western blotting, respectively. Results – (1) Among Caucasian patients, ABCB1 rs60023214 significantly correlated with complete MR (P = 0.005). Concerning SNPs combination in IM uptake transporters, the associations with treatment outcomes were statistically significant for both major and complete MR (P = 0.005 and P = 0.01, respectively). (2) ABCB1 protein was not expressed under any conditions of treatment, differently from ABCG2. Two deregulated miRNAs, namely miR-212 and miR-328, were identified to be inversely correlated with ABCG2 (r2= 0.57; p=0.03 and r2=0.47; p=0.06, respectively). Experiments of loss and gain of function confirmed the functional influence of these miRNAs on ABCG2. Conclusion – The multiple candidate gene approach identified single and combination of SNPs that can be proposed as predictor of IM response. The in vitro study suggested that IM resistance could be mediated by miRNA-dependent mechanism. Further studies are needed to validate these preliminary findings.

The role of mitochondria in the regulation of gamma rays induced mTOR-dependent senescence

The first aims of this study were to demonstrate if mitochondrial biogenesis and senescence can be induced simultaneously in cell lines upon exposure to a genotoxic stress, and if the presence of mtDNA mutations which impair the functionality of respiratory complexes can influence the ability of a cell to activate senescence. The data obtained on the oncocytic model XTC.UC1 demonstrated that the presence of mitochondrial dysfunction is involved in the maintenance of a senescent phenotype induced by γ-rays treatment. The involvement of mTORC1 in the regulation of senescence has been shown in this cell line. On the other hand, in cells which do not present mitochondrial dysfunction it has been verified that genotoxic stress determines the activation of both mitochondrial biogenesis and senescence. Further studies are necessary in order to verify if mitochondrial biogenesis sustains the activation of senescence. The second aim of this thesis was to determine the involvement of mTORC1 in the regulation of PGC-1α expression, in order to verify what is the cause of the development of oncocytoma in patients affected by two hereditary cancer syndromes; Cowden and Birt-hogg-Dubé . The study of oncocytic tumors developed by patients affected by these syndromes suggested that the double heterozigosity of the two causative genes, PTEN and FLCN respectively, induce the activation of mTORC1 and therefore the activation of PGC-1α expression. On XTC.UC1 cell line, the most suitable in vitro model, experiments of complementation of PTEN and FLCN were conducted. To date, these results demonstrated that mTORC1 is not involved in the regulation of PGC-1α expression, and PTEN and FLCN seem to have opposite effect on PGC-1α expression.

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.

Role of αvβ3 – Integrin and TLR2 in the innate response to Herpes Simplex Virus infection and Delivery of retargeted oncolytic Herpes Simplex via carrier cells

In the first part of my thesis I studied the mechanism of initiation of the innate response to HSV-1. Innate immune response is the first line of defense set up by the cell to counteract pathogens infection and it is elicited by the activation of a number of membrane or intracellular receptors and sensors, collectively indicated as PRRs, Patter Recognition Receptors. We reported that the HSV pathogen-associated molecular patterns (PAMP) that activate Toll-like receptor 2 (TLR2) and lead to the initiation of innate response are the virion glycoproteins gH/gL and gB, which constitute the conserved fusion core apparatus across the Herpesvirus. Specifically gH/gL is sufficient to initiate a signaling cascade which leads to NF-κB activation. Then, by gain and loss-of-function approaches, we found that αvβ3-integrin is a sensor of and plays a crucial role in the innate defense against HSV-1. We showed that αvβ3-integrin signals through a pathway that concurs with TLR2, affects activation/induction of interferons type 1, NF-κB, and a polarized set of cytokines and receptors. Thus, we demonstrated that gH/gL is sufficient to induce IFN1 and NF-κB via this pathway. From these data, we proposed that αvβ3-integrin is considered a class of non-TLR pattern recognition receptors. In the second part of my thesis I studied the capacity of human mesenchymal stromal cells isolated by fetal membranes (FM-hMSCs) to be used as carrier cells for the delivery of retargeted R-LM249 virus. The use of systemically administrated carrier cells to deliver oncolytic viruses to tumoral targets is a promising strategy in oncolytic virotherapy. We observed that FM-hMSCs can be infected by R-LM249 and we optimized the infection condition; then we demonstrate that stromal cells sustain the replication of retargeted R-LM249 and spread it to target tumoral cells. From these preliminary data FM-hMSCs resulted suitable to be used as carrier cells

Rationale for an adjunctive therapy with fenofibrate in pharmacoresistant nocturnal frontal lobe epilepsy (NFLE).

Nocturnal Frontal Lobe Epilepsy (NFLE) is characterized by onset during infancy or childhood with persistence in adulthood, family history of similar nocturnal episodes simulating non-REM parasomnias (sleep terrors or sleepwalking), general absence of morphological substrates, often by normal interictal electroencephalographical recordings (EEGs) during wakefulness. A family history of epilepsy may be present with Mendelian autosomal dominant inheritance has been described in some families. Recent studies indicate the involvement of neuronal nicotinic acetylcholine receptors (nAChRs) in the molecular mechanisms of NFLE. Mutations in the genes encoding for the α4 (CHRNA4) and ß2 (CHRNB2) subunits of the nAChR induce changes in the biophysical properties of nAChR, resulting generally in a “gain of function”. Preclinical studies report that activation of a nuclear receptor called type peroxisome proliferator-activated receptor (PPAR-α) by endogenous molecules or by medications (e.g. fenofibrate) reduces the activity of the nAChR and, therefore, may decrease the frequency of seizures. Thus, we hypothesize that negative modulation of nAChRs might represent a therapeutic strategy to be explored for pharmacological treatment of this form of epilepsy, which only partially responds to conventional antiepileptic drugs. In fact, carbamazepine, the current medication for NFLE, abolishes the seizures only in one third of the patients. The aim of the project is: 1)_to verify the clinical efficacy of adjunctive therapy with fenofibrate in pharmacoresistant NFLE and ADNFLE patients; focousing on the analysis of the polysomnographic action of the PPAR- agonist (fenofibrate). 2)_to demonstrate the subtended mechanism of efficacy by means of electrophysiological and behavioral experiments in an animal model of the disease: particularly, transgenic mice carrying the mutation in the nAChR 4 subunit (Chrna4S252F) homologous to that found in the humans. Given that a PPAR-α agonist, FENOFIBRATE, already clinically utilized for lipid metabolism disorders, provides a promising therapeutic avenue in the treatment of NFLE\ADNFLE.