Evolutionary epigenetics of modern human populations

Epigenetic variability is a new mechanism for the study of human microevolution, because it creates both phenotypic diversity within an individual and within population. This mechanism constitutes an important reservoir for adaptation in response to new stimuli and recent studies have demonstrated that selective pressures shape not only the genetic code but also DNA methylation profiles. The aim of this thesis is the study of the role of DNA methylation changes in human adaptive processes, considering the Italian peninsula and macro-geographical areas. A whole-genome analysis of DNA methylation profile across the Italian penisula identified some genes whose methylation levels differ between individuals of different Italian districts (South, Centre and North of Italy). These genes are involved in nitrogen compound metabolism and genes involved in pathogens response. Considering individuals with different macro-geographical origins (individuals of Asians, European and African ancestry) more significant DMRs (differentially methylated regions) were identified and are located in genes involved in glucoronidation, in immune response as well as in cell comunication processes. A "profile" of each ancestry (African, Asian and European) was described. Moreover a deepen analysis of three candidate genes (KRTCAP3, MAD1L and BRSK2) in a cohort of individuals of different countries (Morocco, Nigeria, China and Philippines) living in Bologna, was performed in order to explore genetic and epigenetic diversity. Moreover this thesis have paved the way for the application of DNA methylation for the study of hystorical remains and in particular for the age-estimation of individuals starting from biological samples (such as teeth or blood). Noteworthy, a mathematical model that considered methylation values of DNA extracted from cementum and pulp of living individuals can estimate chronological age with high accuracy (median absolute difference between age estimated from DNA methylation and chronological age was 1.2 years).

Mitochondrial Genomics: structure, inheritance and phylogenetic utility of the Mitochondrial genome in Bivalvia and Insecta

This PhD Thesis is the result of my research activity in the last three years. My main research interest was centered on the evolution of mitochondrial genome (mtDNA), and on its usefulness as a phylogeographic and phylogenetic marker at different taxonomic levels in different taxa of Metazoa. From a methodological standpoint, my main effort was dedicated to the sequencing of complete mitochondrial genomes, and the approach to whole-genome sequencing was based on the application of Long-PCR and shotgun sequences. Moreover, this research project is a part of a bigger sequencing project of mtDNAs in many different Metazoans’ taxa, and I mostly dedicated myself to sequence and analyze mtDNAs in selected taxa of bivalves and hexapods (Insecta). Sequences of bivalve mtDNAs are particularly limited, and my study contributed to extend the sampling. Moreover, I used the bivalve Musculista senhousia as model taxon to investigate the molecular mechanisms and the evolutionary significance of their aberrant mode of mitochondrial inheritance (Doubly Uniparental Inheritance, see below). In Insects, I focused my attention on the Genus Bacillus (Insecta Phasmida). A detailed phylogenetic analysis was performed in order to assess phylogenetic relationships within the genus, and to investigate the placement of Phasmida in the phylogenetic tree of Insecta. The main goal of this part of my study was to add to the taxonomic coverage of sequenced mtDNAs in basal insects, which were only partially analyzed.

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.

Interactions of G-quadruplex binders and Topoisomerase I inhibitors in cancer cells

Top1-DNA cleavage complexes (Top1ccs) trigger an accumulation of antisense RNAPII transcripts specifically at active divergent CpG-island promoters in a replication independent and Top1 dependent manner, leading to transcription-dependent genome instability and altered transcription regulation. Using different cancer cell lines of colon and osteo origins, we show that they display different sensitivity to CPT and G4 binder that is independent from Top1 level. To look at the interactions between Top1 and G4, we show that co-treatment with G4 binders potentiate the cell cytotoxicity of CPT regardless of the treatment sequences. Potentiation is indicated by a reduced inhibition concentration (IC50) with a more profound cytotoxicity in CPT-resistant cell lines, HCT15 and U2OS, hence, indicating an interaction between Top1inhibitor and G4 binders. Moreover, computational analysis confirmed the present of G4 motifs in genes with CPT-induced antisense transcription. G4 motifs are present mostly 5000 bp upstream from transcription start site and notably lower in genes. Comparisons between genes with no antisense transcription and genes with antisense transcription show that G4 motifs in this region are notably lower in the genes with antisense transcripts. Since CPT increases negative supercoils at promoters of intermediate activity, the formation of G4 is also increased in CPT-treated cells. Suprisingly, formation of G4 is regulated in parallel to the transient stabilization of R-loops, indicating a role in response to CPT-induced stress. G4 formation is highly elevated in Pyridostatin treated cells, which previous study shows increased formation of γH2Ax foci. This effect is also seen in the CPT-resistant cell lines, HCT15, indicating that the formation is a general event in response to CPT. We also show that R-loop formation is greatly increased in Pyridostatin treated cells. In order to study the role of R-loops and G4 structures in Top1cc-dependant repair pathway, we inhibited tyrosyl-phosphodiestrase 1 (TDP-1) using a TDP-1 inhibitor.

Parapoxvirus: basi genomiche della patogenesi

Parapoxvirus (PPV) are member of a genus in the family poxviridae which currently encompasses four species: the prototype orf virus (OV), bovine papular stomatitis virus (BPSV), pseudocowpox virus (PCPV) and parapoxvirus of New Zealand red deer (PVNZ). PPVs cause widespread, but localized diseases of small and large ruminants and they can also be transmitted to man. Knowledge of the molecular biology of PPV is still limited as compared to orthopoxviruses, especially vaccinia virus (VACV). The PPV genome displays a high G+C content and relatively small size for poxvirus. Coventional electron microscopy displays PPV virions with ovoid shape and slightly smaller in size than the brickshaped orthopoxviruses. The most striking feature, which readily enables identification of PPV, is a tubule-like structure that surrounds the particle in a spiral fashion. PPV genome organization and content is very similar to that of other poxviruses, the central region contain 88 genes which are present in all poxviruse, in contrast the terminal regions are variable and contain a set of genes unique to the genus PPV. Genes in the near-terminal regions of the genome are frequently not essential for growth in cultured cells encoding factors with important roles in virushost interactions including modulating host immune responses and determining host range. Recently it was suggested that the open reading frames (ORFs) 109 and 110 of the OV genome have a major role in determining species specificity during natural infection in sheep and goats. This hypothesis is based on the analysis of a few number of sequences of different sheep and goats viral isolates. PPV replicate into the cytoplasm of infected cells and produce three structurally different infectious particles: the intracellular mature virions (IMV), intracellular enveloped virions (IEV) and the extracellular enveloped virions (EEV). The vaccinia A33R and A34R hotologue proteins encoded by the ORFS 109 and 110 are expressed in the envelope of the IEV and EEV. The F1L immunodominant protein of orf virus is the major component of the surface tubule structure of the IMV and can post-translationaly insert into membranes via Cterminal, hydrofobic anchor sequence like its orthologue VACV H3L protein. Moreover the F1L protein binds to glycosaminoglycans on the cell surface and has an important role in IMV adsorption to mammalian cells. In this study we investigated the morphogenesis of the PPV through the construction of a mutant virus deleted of the F1L protein. A study of the deleted virus life cycle was conducted in different type of cells and its morphology was observed with electron microscopy. It was demonstared that F1L protein have important role in morphogenesis and infectivity. Moreover it is essential to determine the spiral fashion of the tubule like structure of the virion surface. Some pathogenetic aspects of the PPV infection were studied, in particular the protein implicated in the host range were analysed in detail. An experimental infection with OV and PCPV was conducted in goats and sheep. After infection, the severity of the lesions were comparable in both the animal species. The OV did not result in severe disease neither in sheep nor in goats, suggesting that host factors, rather than virus strain characteristics, may play an important role in the pathogenesis of the Parapoxvirus infections. The PCPV failed to produce any lesion in both sheep and goats, ruling out the possibility of any recombination between PCPV and OV during natural infection in these animal species. The phylogenetic analysis of the ORFs 109 and 110 from several goats and sheep viral isolates showed a clustering based on the antigenic content of the protein that was independent from species and geographic origin.

Farmacogenomica della Clofarabina nel trattamento delle Leucemie Acute pediatriche: identificazione di nuovi bersagli molecolari e del profilo genomico associato all’efficacia terapeutica del farmaco antitumorale

In quest’ultimi decenni si è assistito ad un notevole miglioramento nella terapia delle Leucemie Acute (LA) pediatriche, nonostante tutto si assiste oggi ad una fase di plateau della curva di sopravvivenza e le leucemie continuano a costituire la principale causa di morte pediatrica per malattia. Ulteriori progressi nel trattamento delle LA potrebbero essere ottenuti mediante studi di farmacogenomica che, identificando le componenti genetiche associate alla risposta individuale ai trattamenti farmacologici, consentono il disegno di terapie personalizzate e tumore-specifiche, ad alta efficacia e bassa tossicità per ciascun paziente. Il lavoro svolto è stato, dunque, finalizzato allo studio della farmacogenomica del farmaco antitumorale Clofarabina (CLO) nel trattamento delle LA pediatriche al fine di identificare marcatori genetici predittivi di risposta delle cellule leucemiche al farmaco, delucidare i meccanismi di resistenza cellulare ed individuare nuovi bersagli verso cui indirizzare terapie più mirate ed efficaci. L’analisi in vitro della sensibilità alla CLO di blasti provenienti da pazienti pediatrici affetti da Leucemia Acuta Linfoblastica (LAL) e Mieloide (LAM) ha consentito l’identificazione di due sottopopolazioni di cellule LAL ad immunofenotipo T a diversa sensibilità alla CLO. Mediante DNA-microarrays, si è identificata la “signature” genetica specificamente associata alla diversa risposta delle cellule LAL-T al farmaco. Successivamente, la caratterizzazione funzionale dei geni differenziali e l’analisi dei pathways hanno consentito l’identificazione specifica di potenziali biomarcatori di risposta terapeutica aprendo nuove prospettive per la comprensione dei meccanismi di resistenza cellulare alla CLO e suggerendo un nuovo bersaglio terapeutico per le forme LAL-T a bassa sensibilità al farmaco. In conclusione, nel lavoro svolto si sono identificati set di geni e pathways di rilievo biologico per la risposta delle cellule LAL-T alla CLO suggerendo marcatori genetici in grado di identificare i soggetti eleggibili per il trattamento o verso cui disegnare terapie innovative. Il lavoro è paradigma per l’applicazione della farmacogenomica in altre neoplasie.

Composti perfluorurati: Valutazione degli effetti biologici e molecolari in modelli cellulari

L’acido perfluorottanoico (PFOA) e l’acido perfluoronanoico (PFNA) sono composti perfluorurati (PFCs) comunemente utilizzati nell’industria, negli ultimi 60 anni, per diverse applicazioni. A causa della loro resistenza alla degradazione, questi composti sono in grado di accumularsi nell’ambiente e negli organismi viventi, da cui possono essere assunti in particolare attraverso la dieta. Le esistenti evidenze sugli effetti dell’esposizione negli animali, tra cui la potenziale cancerogenicità, hanno accresciuto l’interesse sui possibili rischi per la salute nell’uomo. Recenti studi sull’uomo indicano che i PFC sono presenti nel siero, con livelli molto alti soprattutto nei lavoratori cronicamente esposti, e sono associati positivamente al cancro al seno e alla prostata. Inoltre, sono state riportate proprietà estrogen-like e variazioni nei livelli di metilazione sui promotori di alcuni geni. L’esposizione in utero è stata associata positivamente a ipometilazione globale del DNA nel siero cordonale. L’obiettivo di questo studio è stato quello di indagare gli effetti dell’esposizione a questi perfluorurati su linee cellulari tumorali e primarie umane (MOLM-13, RPMI, HEPG2, MCF7,WBC, HMEC e MCF12A), appartenenti a diversi tessuti target, utilizzando un ampio range di concentrazioni (3.12 nM - 500 μM). In particolare, si è valutato: la vitalità, il ciclo cellulare, l’espressione genica, la metilazione globale del DNA e la metilazione gene specifica. Dai risultati è emerso come entrambi i perfluorurati abbiano effetti biologici: PFOA presenta un effetto prevalente citostatico, PFNA prevalentemente citotossico. L’effetto è, però, prevalente sulle linee cellulari primarie di epitelio mammario (HMEC, MCF12A), anche a concentrazioni riscontrate in lavoratori cronicamente esposti (≥31,25 µM). Dall’analisi su queste cellule primarie, non risultano variazioni significative della metilazione globale del DNA alle concentrazioni di 15,6 e 31,25 µM. Emergono invece variazioni sui geni marcatori del cancro al seno, del ciclo cellulare, dell’apoptosi, del pathway di PPAR-α e degli estrogeni, ad una concentrazione di 31,25 µM di entrambi i PFCs.