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.

Analysis of Copy Number Variants identifies new candidate genes for Autism Spectrum Disorder and Intellectual Disability

Autism spectrum disorder (ASD) and Intellectual Disability (ID) are complex neuropsychiatric disorders characterized by extensive clinical and genetic heterogeneity and with overlapping risk factors. The aim of my project was to further investigate the role of Copy Numbers Variants (CNVs), identified through genome-wide studies performed by the Autism Geome Project (AGP) and the CHERISH consortium in large cohorts of ASD and ID cases, respectively. Specifically, I focused on four rare genic CNVs, selected on the basis of their impact on interesting ASD/ID candidate genes: a) a compound heterozygous deletion involving CTNNA3, predicted to cause the lack of functional protein; b) a 15q13.3 duplication containing CHRNA7; c) a 2q31.1 microdeletion encompassing KLHL23, SSB and METTL5; d) Lastly, I investigated the putative imprinting regulation of the CADPS2 gene, disrupted by a maternal deletion in two siblings with ASD and ID. This study provides further evidence for the role of CTNNA3, CHRNA7, KLHL23 and CADPS2 as ASD and/or ID susceptibility genes, and highlights that rare genetic variation contributes to disease risk in different ways: some rare mutations, such as those impacting CTNNA3, act in a recessive mode of inheritance, while other CNVs, such as those occurring in the 15q13.3 region, are implicated in multiple developmental and/or neurological disorders possibly interacting with other susceptibility variants elsewhere in the genome. On the other hand, the discovery of a tissue-specific monoallelic expression for the CADPS2 gene, implicates the involvement of epigenetic regulatory mechanisms as risk factors conferring susceptibility to ASD/ID.

Studio di Mutazioni Geniche e Variazioni Copy-Number implicate nel deficit combinato degli Ormoni Ipofisari(CPHD)

La crescita normale di un individuo è il risultato dell’azione coordinata di molteplici ormoni e recettori codificati da geni e a tal proposito, discreto interesse è stato dato ai geni tipici dell’asse del GH. Tuttavia altri geni, più a monte di questi e responsabili dello sviluppo dell’ipofisi contribuiscono alla crescita normale o patologica. Alcuni geni studiati sono POU1F1, PROP1, LHX3, LHX4, HESX1, SOX3 e svariate loro mutazioni sono state identificate come causa di panipopituarismo (CPHD=Combined Pituitary Hormone Deficiency). In realtà la ricerca genetica non spiega ancora molte anomalie ipofisarie e molte mutazioni devono ancora essere identificate. Uno degli scopi del dottorato, svoltosi nel laboratorio di Genetica molecolare di Pediatria, è stata l’identificazione di mutazioni geniche da un gruppo di pazienti CPHD considerando in particolare i geni POU1F1, LHX3, SOX3, non ancora messi a punto presso il laboratorio. L’approccio sperimentale si è basato sulle seguenti fasi: prelievo delle informazioni di sequenza da GeneBank, progettazione di primers per amplificare le porzioni esoniche, messa a punto delle fasi della PCR e del sequenziamento, analisi della sequenza e confronto con le informazioni di sequenza depositate allo scopo di rintracciare eventuali mutazioni o varianti. La bassa percentuale di mutazioni in questi geni non ha permesso finora di rintracciare mutazioni nelle porzioni esoniche salvo che in un soggetto, nell’esone 6 di LHX3b (nuova mutazione, recessiva eterozigote, c.1248A>G implicata nella mutazione p.T377A della sequenza proteica). Un metodo di screening di questa mutazione impiegando l’enzima di restrizione SacII è stato usato, senza rilevare nessun altra occorrenza dell’allele mutato in 53 soggetti di controllo. Oltre alla messa a punto del sequenziamento e di alcune tecniche di analisi di singoli SNP o piccoli INDELs per i 3 geni, la ricerca svolta è stata orientata all’impiego di metodi di rilevamento di riarrangiamenti genetici comportanti ampie delezioni e/o variazioni del copy-number di esoni/interi geni detto MLPA (Multiplex Ligation-dependent Probe Amplification) e progettato da MRC-Holland. Il sequenziamento infatti non permette di rilevare tali alterazioni quando sono ampie ed in eterozigosi. Per esempio, in un’ampia delezione in eterozigosi, l’intervallo delimitato dai primers usati per la PCR può non includere totalmente la porzione interessata da delezione su un cromosoma cosicché la PCR ed il sequnziamento si basano solo sulle informazioni dell’altro cromosoma non deleto. Un vantaggio della tecnica MLPA, è l’analisi contemporanea di una quarantina di siti posti su svariati geni. Questa metodo tuttavia può essere affetto da un certo margine di errore spesso dipendente dalla qualità del DNA e dovrebbe essere affiancato e validato da altre tecniche più impegnativa dal punto di vista sperimentale ma più solide, per esempio la Real Time PCR detta anche PCR quantitativa (qPCR). In laboratorio, grazie all’MLPA si è verificata la condizione di delezione eterozigote di un paziente “storico” per il gene GH1 e la stessa mutazione è stata rilevata anche con la qPCR usando lo strumento Corbett Rotor Gene 6000 (Explera). Invece un’analisi solo con la qPCR di variazioni del copy-number (CNV) per SOX3 in pazienti maschili non ha ancora evidenziato anomalie. Entrambe le tecniche hanno aspetti interessanti, il miglior approccio al momento sembra un’analisi iniziale di pazienti con l’MLPA, seguita dalla verifica di un eventuale esito anomalo impiegando la real-time PCR.

Analisi di “Copy Number Variants” ed identificazione di nuovi geni candidati per l’Autismo e Ritardo Mentale

I disturbi dello spettro autistico (DSA) ed il ritardo mentale (RM) sono caratterizzati da un’eziologia genetica complessa ed eterogenea. Grazie ai recenti sviluppi nella ricerca genomica, è stato possibile dimostrare il ruolo di numerose copy number variants (CNVs) nella patogenesi di questi disturbi, anche se nella maggior parte dei casi l’eziologia rimane ancora sconosciuta. Questo lavoro riguarda l’identificazione e la caratterizzazione dei CNVs in famiglie con DSA e RM. E’ stata studiata una microdelezione in 7q31 che coinvolge i geni IMMP2L e DOCK4, trasmessa dalla madre con dislessia a due figli con autismo ed una figlia con dislessia. Nella stessa famiglia segrega una seconda microdelezione in 2q14 che inattiva il gene CNTNAP5 ed è trasmessa dal padre (con tratti autistici) ai due figli con autismo. Abbiamo quindi ipotizzato che i geni DOCK4 e CNTNAP5 potessero essere implicati, rispettivamente, nella suscettibilità a dislessia e DSA. Lo screening di numerosi individui affetti ha supportato la nostra ipotesi, con l’identificazione di una nuova microdelezione di DOCK4 che segrega con la dislessia, e 3 nuove varianti missenso in CNTNAP5 in individui con autismo. Dall’analisi genomica comparativa su array (aCGH) di individui con RM, è stata identificata una delezione nella regione 7q31.32, che coinvolge il gene CADPS2, in due fratelli con RM e tratti autistici, probabilmente ereditata dalla madre. Lo screening di mutazione di questo gene in individui con autismo o RM, ha portato all’identificazione di 3 varianti non sinonime, assenti nei controlli, ed ereditate per via materna. Poiché CADPS2 risiede in una regione genomica che contiene loci soggetti ad imprinting, abbiamo ipotizzato che il gene CADPS2 possa essere anch’esso caratterizzato da imprinting, con espressione monoallelica materna. Lo studio di espressione di CADPS2 in cellule del sangue ha avvalorato questa ipotesi, implicando perciò CADPS2 come un nuovo gene di suscettibilità per il RM e DSA.

Genetica molecolare dell'autismo: studi di associazione ed analisi di geni candidati

Autism is a neurodevelpmental disorder characterized by impaired verbal communication, limited reciprocal social interaction, restricted interests and repetitive behaviours. Twin and family studies indicate a large genetic contribution to ASDs (Autism Spectrum Disorders). During my Ph.D. I have been involved in several projects in which I used different genetic approaches in order to identify susceptibility genes in autism on chromosomes 2, 7 and X: 1)High-density SNP association and CNV analysis of two Autism Susceptibility Loci. The International Molecular Genetic Study of Autism Consortium (IMGSAC) previously identified linkage loci on chromosomes 7 and 2, termed AUTS1 and AUTS5, respectively. In this study, we evaluated the patterns of linkage disequilibrium (LD) and the distribution of haplotype blocks, utilising data from the HapMap project, across the two strongest peaks of linkage on chromosome 2 and 7. More than 3000 SNPs have been selected in each locus in all known genes, as well as SNPs in non-genic highly conserved sequences. All markers have been genotyped to perform a high-density association analysis and to explore copy number variation within these regions. The study sample consisted of 127 and 126 multiplex families, showing linkage to the AUTS1 and AUTS5 regions, respectively, and 188 gender-matched controls. Association and CNV analysis implicated several new genes, including IMMP2L and DOCK4 on chromosome 7 and ZNF533 and NOSTRIN on the chromosome 2. Particularly, my contribution to this project focused on the characterization of the best candidate gene in each locus: On the AUTS5 locus I carried out a transcript study of ZNF533 in different human tissues to verify which isoforms and start exons were expressed. High transcript variability and a new exon, never described before, has been identified in this analysis. Furthermore, I selected 31 probands for the risk haplotype and performed a mutation screen of all known exons in order to identify novel coding variants associated to autism. On the AUTS1 locus a duplication was detected in one multiplex family that was transmitted from father to an affected son. This duplication interrupts two genes: IMMP2L and DOCK4 and warranted further analysis. Thus, I performed a screening of the cohort of IMGSAC collection (285 multiplex families), using a QMPSF assay (Quantitative Multiplex PCR of Short fluorescent Fragments) to analyse if CNVs in this genic region segregate with autism phenotype and compare their frequency with a sample of 475 UK controls. Evidence for a role of DOCK4 in autism susceptibility was supported by independent replication of association at rs2217262 and the finding of a deletion segregating in a sib-pair family. 2)Analysis of X chromosome inactivation. Skewed X chromosome inactivation (XCI) is observed in females carrying gene mutations involved in several X-linked syndromes. We aimed to estimate the role of X-linked genes in ASD susceptibility by ascertaining the XCI pattern in a sample of 543 informative mothers of children with ASD and in a sample of 164 affected girls. The study sample included families from different european consortia. I analysed the XCI inactivation pattern in a sample of italian mothers from singletons families with ASD and also a control groups (144 adult females and 40 young females). We observed no significant excess of skewed XCI in families with ASD. Interestingly, two mothers and one girl carrying known mutations in X-linked genes (NLGN3, ATRX, MECP2) showed highly skewed XCI, suggesting that ascertainment of XCI could reveal families with X-linked mutations. Linkage analysis was carried out in the subgroup of multiplex families with skewed XCI (≥80:20) and a modest increased allele sharing was obtained in the Xq27-Xq28 region, with a peak Z score of 1.75 close to rs719489. In this region FMR1 and MECP2 have been associated in some cases with austim and therefore represent candidates for the disorder. I performed a mutation screen of MECP2 in 33 unrelated probands from IMGSAC and italian families, showing XCI skewness. Recently, Xq28 duplications including MECP2, have been identified in families with MR, with asymptomatic carrier females showing extreme (>85%) skewing of XCI. For these reason I used the sample of probands from X-skewed families to perform CNV analysis by Real-time quantitative PCR. No duplications have been found in our sample. I have also confirmed all data using as alternative method the MLPA assay (Multiplex Ligation dependent Probe Amplification). 3)ASMT as functional candidate gene for autism. Recently, a possible involvement of the acetylserotonin O-methyltransferase (ASMT) gene in susceptibility to ASDs has been reported: mutation screening of the ASMT gene in 250 individuals from the PARIS collection revealed several rare variants with a likely functional role; Moreover, significant association was reported for two SNPs (rs4446909 and rs5989681) located in one of the two alternative promoters of the gene. To further investigate these findings, I carried out a replication study using a sample of 263 affected individuals from the IMGSAC collection and 390 control individuals. Several rare mutations were identified, including the splice site mutation IVS5+2T>C and the L326F substitution previously reported by Melke et al (2007), but the same rare variants have been found also in control individuals in our study. Interestingly, a new R319X stop mutation was found in a single autism proband of Italian origin and is absent from the entire control sample. Furthermore, no replication has been found in our case-control study typing the SNPs on the ASMT promoter B.

Molecular bases, pathogenic mechanisms and possible therapeutic approach in Leber's Hereditary Optic Neuropathy

Leber’s hereditary optic neuropathy (LHON) is a mitochondrial disease characterized by a rapid loss of central vision and optic atrophy, due to the selective degeneration of retinal ganglion cells. The age of onset is around 20, and the degenerative process is fast and usually the second eye becomes affected in weeks or months. Even if this pathology is well known and has been well characterized, there are still open questions on its pathophysiology, such as the male prevalence, the incomplete penetrance and the tissue selectivity. This maternally inherited disease is caused by mutations in mitochondrial encoded genes of NADH ubiquinone oxidoreductase (complex I) of the respiratory chain. The 90% of LHON cases are caused by one of the three common mitochondrial DNA mutations (11778/ND4, 14484/ND6 and 3460/ND1) and the remaining 10% is caused by rare pathogenic mutations, reported in literature in one or few families. Moreover, there is also a small subset of patients reported with new putative pathogenic nucleotide changes, which awaits to be confirmed. We here clarify some molecular aspects of LHON, mainly the incomplete penetrance and the role of rare mtDNA mutations or variants on LHON expression, and attempt a possible therapeutic approach using the cybrids cell model. We generated novel structural models for mitochondrial encoded complex I subunits and a conservation analysis and pathogenicity prediction have been carried out for LHON reported mutations. This in-silico approach allowed us to locate LHON pathogenic mutations in defined and conserved protein domains and can be a useful tool in the analysis of novel mtDNA variants with unclear pathogenic/functional role. Four rare LHON pathogenic mutations have been identified, confirming that the ND1 and ND6 genes are mutational hot spots for LHON. All mutations were previously described at least once and we validated their pathogenic role, suggesting the need for their screening in LHON diagnostic protocols. Two novel mtDNA variants with a possible pathogenic role have been also identified in two independent branches of a large pedigree. Functional studies are necessary to define their contribution to LHON in this family. It also been demonstrated that the combination of mtDNA rare polymorphic variants is relevant in determining the maternal recurrence of myoclonus in unrelated LHON pedigrees. Thus, we suggest that particular mtDNA backgrounds and /or the presence of specific rare mutations may increase the pathogenic potential of the primary LHON mutations, thereby giving rise to the extraocular clinical features characteristic of the LHON “plus” phenotype. We identified the first molecular parameter that clearly discriminates LHON affected individuals from asymptomatic carriers, the mtDNA copy number. This provides a valuable mechanism for future investigations on variable penetrance in LHON. However, the increased mtDNA content in LHON individuals was not correlated to the functional polymorphism G1444A of PGC-1 alpha, the master regulator of mitochondrial biogenesis, but may be due to gene expression of genes involved in this signaling pathway, such as PGC-1 alpha/beta and Tfam. Future studies will be necessary to identify the biochemical effects of rare pathogenic mutations and to validate the novel candidate mutations here described, in terms of cellular bioenergetic characterization of these variants. Moreover, we were not able to induce mitochondrial biogenesis in cybrids cell lines using bezafibrate. However, other cell line models are available, such as fibroblasts harboring LHON mutations, or other approaches can be used to trigger the mitochondrial biogenesis.

Intrinsic uncoupling in the ATP synthase of Escherichia coli. Studies on WT and ε-truncated mutants

The H+/ATP ratio in the catalysis of ATP synthase has generally been considered a fixed parameter. However, Melandri and coworkers have recently shown that, in the ATP synthase of the photosynthetic bacterium Rb.capsulatus, this ratio can significantly decrease during ATP hydrolysis when the concentration of either ADP or Pi is maintained at a low level (Turina et al., 2004). The present work has dealt with the ATP synthase of E.coli, looking for evidence of this phenomenon of intrinsic uncoupling in this organism as well. First of all, we have shown that the DCCD-sensitive ATP hydrolysis activity of E.coli internal membranes was strongly inhibited by ADP and Pi, with a half-maximal effect in the submicromolar range for ADP and at 140 µM for Pi. In contrast to this monotonic inhibition, however, the proton pumping activity of the enzyme, as estimated under the same conditions by the fluorescence quenching of the ΔpH-sensitive probe ACMA, showed a clearly biphasic progression, both for Pi, increasing from 0 up to approximately 200 µM, and for ADP, increasing from 0 up to a few µM. We have interpreted these results as indicating that the occupancy of ADP and Pi binding sites shifts the enzyme from a partially uncoupled state to a fully coupled state, and we expect that the ADP- and Pi-modulated intrinsic uncoupling is likely to be a general feature of prokaryotic ATP synthases. Moreover, the biphasicity of the proton pumping data suggested that two Pi binding sites are involved. In order to verify whether the same behaviour could be observed in the isolated enzyme, we have purified the ATP synthase of E.coli and reconstituted it into liposomes. Similarly as observed in the internal membrane preparation, in the isolated and reconstituted enzyme it was possible to observe inhibition of the hydrolytic activity by ADP and Pi (with half-maximal effects at few µM for ADP and at 400 µM for Pi) with a concomitant stimulation of proton pumping. Both the inhibition of ATP hydrolysis and the stimulation of proton pumping as a function of Pi were lost upon ADP removal by an ADP trap. These data have made it possible to conclude that the results obtained in E.coli internal membranes are not due to the artefactual interference of enzymatic activities other than the ones of the ATP synthase. In addition, data obtained with liposomes have allowed a calibration of the ACMA signal by ΔpH transitions of known extent, leading to a quantitative evaluation of the proton pumping data. Finally, we have focused our efforts on searching for a possible structural candidate involved in the phenomenon of intrinsic uncoupling. The ε-subunit of the ATP-synthase is known as an endogenous inhibitor of the hydrolysis activity of the complex and appears to undergo drastic conformational changes between a non-inhibitory form (down-state) and an inhibitory form (up-state)(Rodgers & Wilce, 2000; Gibbons et al., 2000). In addition, the results of Cipriano & Dunn (2006) indicated that the C-terminal domain of this subunit played an important role in the coupling mechanism of the pump, and those of Capaldi et al. (2001), Suzuki et al. (2003) were consistent with the down-state showing a higher hydrolysis-to-synthesis ratio than the up-state. Therefore, we decided to search for modulation of pumping efficiency in a C-terminally truncated ε mutant. A low copy number expression vector has been built, carrying an extra copy of uncC, with the aim of generating an ε-overexpressing E.coli strain in which normal levels of assembly of the mutated ATP-synthase complex would be promoted. We have then compared the ATP hydrolysis and the proton pumping activity in membranes prepared from these ε-overexpressing E.coli strains, which carried either the WT ε subunit or the ε88-stop truncated form. Both strains yielded well energized membranes. Noticeably, they showed a marked difference in the inhibition of hydrolysis by Pi, this effect being largely lost in the truncated mutant. However, pre-incubation of the mutated enzyme with ADP at low nanomolar concentrations (apparent Kd = 0.7nM) restored the hydrolysis inhibition, together with the modulation of intrinsic uncoupling by Pi, indicating that, contrary to wild-type, during membrane preparation the truncated mutant had lost the ADP bound at this high-affinity site, evidently due to a lower affinity (and/or higher release) for ADP of the mutant relative to wild type. Therefore, one of the effects of the C-terminal domain of ε appears to be to modulate the affinity of at least one of the binding sites for ADP. The lack of this domain does not appear so much to influence the modulability of coupling efficiency, but instead the extent of this modulation. At higher preincubated ADP concentrations (apparent Kd = 117nM), the only observed effects were inhibition of both hydrolysis and synthesis, providing a direct proof that two ADP-binding sites on the enzyme are involved in the inhibition of hydrolysis, of which only the one at higher affinity also modulates the coupling efficiency.

Profilo genomico e di espressione dell'Osteosarcoma. Identificazione dei geni con alto grado di correlazione tra numero di copie ed espressione

To identify the regions of recurrent copy number abnormality in osteosarcoma and their effect on gene expression, we performed an integrated genome-wide high-resolution array CGH (aCGH) and gene expression profiling analysis on 40 human OS tissues and 12 OS cell lines. This analysis identified several recurrent chromosome regions that contain genes that show a gene dosage effect on gene expression. A further search, performed on those genes that were over-expressed and localized in the frequently amplified chromosomal regions, greatly reduced the number of candidate genes while their characterization using gene ontology (GO) analysis suggests the importance of the deregulation of the G1-to-S phase in the development of the disease. We also identified frequent deletions on 3q in the vicinity of LSAMP and performed a fine mapping analysis of the breakpoints. We precisely mapped the breakpoints in several instances and demonstrated that the majority do not involve the LSAMP gene itself, and that they appear to form by a process of non-homologous end joining. In addition, aCGH analysis revealed frequent gains of IGF1R that were highly correlated with its protein level. Blockade of IGF1R in OS cell lines with high copy number gain led to growth inhibition suggesting that IGF1R may be a viable drug target in OS, particularly in patients with copy number driven overexpression of this receptor.

Transposable elements: structure and dynamic of the autonomous retrotransposon R2 in Arthropoda with non-canonical reproduction

Eukaryotic ribosomal DNA constitutes a multi gene family organized in a cluster called nucleolar organizer region (NOR); this region is composed usually by hundreds to thousands of tandemly repeated units. Ribosomal genes, being repeated sequences, evolve following the typical pattern of concerted evolution. The autonomous retroelement R2 inserts in the ribosomal gene 28S, leading to defective 28S rDNA genes. R2 element, being a retrotransposon, performs its activity in the genome multiplying its copy number through a “copy and paste” mechanism called target primed reverse transcription. It consists in the retrotranscription of the element’s mRNA into DNA, then the DNA is integrated in the target site. Since the retrotranscription can be interrupted, but the integration will be carried out anyway, truncated copies of the element will also be present in the genome. The study of these truncated variants is a tool to examine the activity of the element. R2 phylogeny appears, in general, not consistent with that of its hosts, except some cases (e.g. Drosophila spp. and Reticulitermes spp.); moreover R2 is absent in some species (Fugu rubripes, human, mouse, etc.), while other species have more R2 lineages in their genome (the turtle Mauremys reevesii, the Japanese beetle Popilia japonica, etc). R2 elements here presented are isolated in 4 species of notostracan branchiopods and in two species of stick insects, whose reproductive strategies range from strict gonochorism to unisexuality. From sequencing data emerges that in Triops cancriformis (Spanish gonochoric population), in Lepidurus arcticus (two putatively unisexual populations from Iceland) and in Bacillus rossius (gonochoric population from Capalbio) the R2 elements are complete and encode functional proteins, reflecting the general features of this family of transposable elements. On the other hand, R2 from Italian and Austrian populations of T. cancriformis (respectively unisexual and hermaphroditic), Lepidurus lubbocki (two elements within the same Italian population, gonochoric but with unfunctional males) and Bacillus grandii grandii (gonochoric population from Ponte Manghisi) have sequences that encode incomplete or non-functional proteins in which it is possible to recognize only part of the characteristic domains. In Lepidurus couesii (Italian gonochoric populations) different elements were found as in L. lubbocki, and the sequencing is still in progress. Two hypothesis are given to explain the inconsistency of R2/host phylogeny: vertical inheritance of the element followed by extinction/diversification or horizontal transmission. My data support previous study that state the vertical transmission as the most likely explanation; nevertheless horizontal transfer events can’t be excluded. I also studied the element’s activity in Spanish populations of T. cancriformis, in L. lubbocki, in L. arcticus and in gonochoric and parthenogenetic populations of B. rossius. In gonochoric populations of T. cancriformis and B. rossius I found that each individual has its own private set of truncated variants. The situation is the opposite for the remaining hermaphroditic/parthenogenetic species and populations, all individuals sharing – in the so far analyzed samples - the majority of variants. This situation is very interesting, because it isn’t concordant with the Muller’s ratchet theory that hypothesizes the parthenogenetic populations being either devoided of transposable elements or TEs overloaded. My data suggest a possible epigenetic mechanism that can block the retrotransposon activity, and in this way deleterious mutations don’t accumulate.

Vanillin production from ferulic acid with Pseudomonas fluorescens BF13-1p4

Bioconversion of ferulic acid to vanillin represents an attractive opportunity for replacing synthetic vanillin with a bio-based product, that can be label “natural”, according to current food regulations. Ferulic acid is an abundant phenolic compound in cereals processing by-products, such as wheat bran, where it is linked to the cell wall constituents. In this work, the possibility of producing vanillin from ferulic acid released enzymatically from wheat bran was investigated by using resting cells of Pseudomonas fluorescens strain BF13-1p4 carrying an insertional inactivation of vdh gene and ech and fcs BF13 genes on a low copy number plasmid. Process parameters were optimized both for the biomass production phase and the bioconversion phase using food-grade ferulic acid as substrate and the approach of changing one variable while fixing the others at a certain level followed by the response surface methodology (RSM). Under optimized conditions, vanillin up to 8.46 mM (1.4 g/L) was achieved, whereas highest productivity was 0.53 mmoles vanillin L-1 h-1). Cocktails of a number of commercial enzyme (amylases, xylanases, proteases, feruloyl esterases) combined with bran pre-treatment with steam explosion and instant controlled pressure drop technology were then tested for the release of ferulic acid from wheat bran. The highest ferulic acid release was limited to 15-20 % of the ferulic acid occurring in bran, depending on the treatment conditions. Ferulic acid 1 mM in enzymatic hydrolyzates could be bioconverted into vanillin with molar yield (55.1%) and selectivity (68%) comparable to those obtained with food-grade ferulic acid after purification from reducing sugars with a non polar adsorption resin. Further improvement of ferulic acid recovery from wheat bran is however required to make more attractive the production of natural vanillin from this by-product.

HIF1α regulates Mitochondrial biogenesis and Cellular senescence induced by Gamma Radiation

Cellular response to γ-rays is mediated by ATM-p53 axis. When p53 is phosphorylated, it can transactivate several genes to induce permanent cell cycle arrest (senescence) or apoptosis. Epithelial and mesenchymal cells are more resistant to radiation-induced apoptosis and respond mainly by activating senescence. Hence, tumor cells in a senescent state might remain as “dormant” malignant in fact through disruption of p53 function, cells may overcome growth arrest. Oncocytic features were acquired in the recurring neoplasia after radiation therapy in patient with colonrectal cancer. Oncocytic tumors are characterized by aberrant biogenesis and are mainly non-aggressive neoplasms. Their low proliferation degree can be explained by chronic destabilization of HIF1α, which presides to adaptation to hypoxia and also plays a pivotal role in hypoxia-related radio-resistance. The aim of the present thesis was to verify whether mitochondrial biogenesis can be induced following radiation treatment, in relation of HIF1α status and whether is predictive of a senescence response. In this study was demonstrate that mitochondrial biogenesis parameters like mitochondrial DNA copy number could be used for the prediction of hypoxic status of tissue after radiation treatment. γ-rays induce an increase of mitochondrial mass and function, in response to a genotoxic stress that pushes cells into senescence. Mitochondrial biogenesis is only indirectly regulated by p53, whose activation triggers a MDM2-mediated HIF1α degradation, leading to the release of PGC-1β inhibition by HIF1α. On the other hand, this protein blunts the mitochondrial response to γ-rays as well as the induction of p21-mediated cell senescence, indicating prevalence of the hypoxic over the genotoxic response. Finally in vivo, post-radiotherapy mtDNA copy number increase well correlates with lack of HIF1α increase in the tissue, concluding this may be a useful molecular tool to infer the trigger of a hypoxic response during radiotherapy, which may lead to failure of activation of senescence.

Mitochondrial DNA rearrangements during human aging: a study on liver, muscle and adipose tissue

Aging is characterized by a chronic, low-grade inflammatory state called “inflammaging”. Mitochondria are the main source of reactive oxygen species (ROS), which trigger the production of pro-inflammatory molecules. We are interested in studying the age-related modifications of the mitochondrial DNA (mtDNA), which can be affected by the lifelong exposure to ROS and are responsible of mitochondrial dysfunction. Moreover, increasing evidences show that telomere shortening, naturally occurring with aging, is involved in mtDNA damage processes and thus in the pathogenesis of age-related disorders. Thus the primary aim of this thesis was the analysis of mtDNA copy number, deletion level and integrity in different-age human biopsies from liver, vastus lateralis skeletal muscle of healthy subjects and patients with limited mobility of lower limbs (LMLL), as well as adipose tissue. The telomere length and the expression of nuclear genes related to mitobiogenesis, fusion and fission, mitophagy, mitochondrial protein quality control system, hypoxia, production and protection from ROS were also evaluated. In liver the decrease in mtDNA integrity with age is accompanied with an increase in mtDNA copy number, suggesting the existence of a “compensatory mechanism” able to maintain the functionality of this organ. Different is the case of vastus lateralis muscle, where any “compensatory pathway” is activated and mtDNA integrity and copy number decrease with age, both in healthy subjects and in patients. Interestingly, mtDNA rearrangements do not incur in adipose tissue with advancing age. Finally, in all tissues a marked gender difference appears, suggesting that aging and also gender diversely affect mtDNA rearrangements and telomere length in the three human tissues considered, likely depending on their different metabolic needs and inflammatory status.

Characterization of mitochondrial genomes in bivalve species with doubly uniparental inheritance of mitochondria

Many bivalve species possess two distinct mtDNA lineages, called F and M, respectively inherited maternally and paternally: this system is called doubly uniparental inheritance (DUI). The main experimental project of my PhD was the quantification of the two mtDNAs during the development of the DUI species Ruditapes philippinarum, from early embryos to sub-adults, using Real-Time qPCR. I identified the time interval in which M mtDNA is lost from female individuals, while it is retained in males (which are heteroplasmic through all of their life cycle). The results also suggested absence of mtDNA replication during early embryogenesis, a process constituting a bottleneck that highly reduces the copy number of mtDNA molecules in cells of developing larvae. In males this bottleneck may produce cells homoplasmic for M mtDNA, and could be considered as a first step of the segregation of M in the male germ line. Another finding was the characterization, in young clams approaching the first reproductive season, of a significant boost in copy number of F mtDNA in females and of M in males. Given the age of animals in which this mtDNA-specific growth was observed, the finding could probably be the outcome of the first round of gonads and gametes production. Other lines of research included the characterization of the unassigned regions in mt genomes of DUI bivalves. These regions can harbor signals involved in the control of replication and/or transcription of the mtDNA molecule, as well as additional open reading frames (ORFs) not related to oxidative phosphorylation. These features in DUI species could be associated to the maintenance of separate inheritance routes for the two mtDNAs. Additional ORFs are also found in other animal mt genomes: I summarized the presence of gene duplications as a co-author in a review focusing on animal mt genomes with unusual gene content.

Analysis of the pig genome for the identification of genomic regions affecting production traits

The aim of this work was to identify markers associated with production traits in the pig genome using different approaches. We focused the attention on Italian Large White pig breed using Genome Wide Association Studies (GWAS) and applying a selective genotyping approach to increase the power of the analyses. Furthermore, we searched the pig genome using Next Generation Sequencing (NSG) Ion Torrent Technology to combine selective genotyping approach and deep sequencing for SNP discovery. Other two studies were carried on with a different approach. Allele frequency changes for SNPs affecting candidate genes and at Genome Wide level were analysed to identify selection signatures driven by selection program during the last 20 years. This approach confirmed that a great number of markers may affect production traits and that they are captured by the classical selection programs. GWAS revealed 123 significant or suggestively significant SNP associated with Back Fat Thickenss and 229 associated with Average Daily Gain. 16 Copy Number Variant Regions resulted more frequent in lean or fat pigs and showed that different copies of those region could have a limited impact on fat. These often appear to be involved in food intake and behavior, beside affecting genes involved in metabolic pathways and their expression. By combining NGS sequencing with selective genotyping approach, new variants where discovered and at least 54 are worth to be analysed in association studies. The study of groups of pigs undergone to stringent selection showed that allele frequency of some loci can drastically change if they are close to traits that are interesting for selection schemes. These approaches could be, in future, integrated in genomic selection plans.

Genotypic variation in Actinidia deliciosa fruit size and carbohydrate content

In a global and increasingly competitive fresh produce market, more attention is being given to fruit quality traits and consumer satisfaction. Kiwifruit occupies a niche position in the worldwide market, when compared to apples, oranges or bananas. It is a fruit with extraordinarily good nutritional traits, and its benefits to human health have been widely described. Until recently, international trade in kiwifruit was restricted to a single cultivar, but different types of kiwifruit are now becoming available in the market. Effective programmes of kiwifruit improvement start by considering the requirements of consumers, and recent surveys indicate that sweeter fruit with better flavour are generally preferred. There is a strong correlation between at-harvest dry matter and starch content, and soluble solid concentration and flavour when fruit are eating ripe. This suggests that carbon accumulation strongly influences the development of kiwifruit taste. The overall aim of the present study was to determine what factors affect carbon accumulation during Actinidia deliciosa berry development. One way of doing this is by comparing kiwifruit genotypes that differ greatly in their ability to accumulate dry matter in their fruit. Starch is the major component of dry matter content. It was hypothesized that genotypes were different in sink strength. Sink strength, by definition, is the effect of sink size and sink activity. Chapter 1 reviews fruit growth, kiwifruit growth and development and carbon metabolism. Chapter 2 describes the materials and methods used. Chapter 3, 4, 5 and 6 describes different types of experimental work. Chapter 7 contains the final discussions and the conclusions Three Actinidia deliciosa breeding populations were analysed in detail to confirm that observed differences in dry matter content were genetically determined. Fruit of the different genotypes differed in dry matter content mainly because of differences in starch concentrations and dry weight accumulation rates, irrespective of fruit size. More detailed experiments were therefore carried out on genotypes which varied most in fruit starch concentrations to determine why sink strengths were so different. The kiwifruit berry comprises three tissues which differ in dry matter content. It was initially hypothesised that observed differences in starch content could be due to a larger proportion of one or other of these tissues, for example, of the central core which is highest in dry matter content. The study results showed that this was not the case. Sink size, intended as cell number or cell size, was then investigated. The outer pericarp makes up about 60% of berry weight in ‘Hayward’ kiwifruit. The outer pericarp contains two types of parenchyma cells: large cells with low starch concentration, and small cells with high starch concentration. Large cell, small cell and total cell densities in the outer pericarp were shown to be not correlated with either dry matter content or fruit size but further investigation of volume proportion among cell types seemed justified. It was then shown that genotypes with fruit having higher dry matter contents also had a higher proportion of small cells. However, the higher proportion of small cell volume could only explain half of the observed differences in starch content. So, sink activity, intended as sucrose to starch metabolism, was investigated. In transiently starch storing sinks, such as tomato fruit and potato tubers, a pivotal role in carbon metabolism has been attributed to sucrose cleaving enzymes (mainly sucrose synthase and cell wall invertase) and to ADP-glucose pyrophosphorylase (the committed step in starch synthesis). Studies on tomato and potato genotypes differing in starch content or in final fruit soluble solid concentrations have demonstrated a strong link with either sucrose synthase or ADP-glucose pyrophosphorylase, at both enzyme activity and gene expression levels, depending on the case. Little is known about sucrose cleaving enzyme and ADP-glucose pyrophosphorylase isoforms. The HortResearch Actinidia EST database was then screened to identify sequences putatively encoding for sucrose synthase, invertase and ADP-glucose pyrophosphorylase isoforms and specific primers were designed. Sucrose synthase, invertase and ADP-glucose pyrophosphorylase isoform transcript levels were anlayzed throughout fruit development of a selection of four genotypes (two high dry matter and two low dry matter). High dry matter genotypes showed higher amounts of sucrose synthase transcripts (SUS1, SUS2 or both) and higher ADP-glucose pyrophosphorylase (AGPL4, large subunit 4) gene expression, mainly early in fruit development. SUS1- like gene expression has been linked with starch biosynthesis in several crop (tomato, potato and maize). An enhancement of its transcript level early in fruit development of high dry matter genotypes means that more activated glucose (UDP-glucose) is available for starch synthesis. This can be then correlated to the higher starch observed since soon after the onset of net starch accumulation. The higher expression level of AGPL4 observed in high dry matter genotypes suggests an involvement of this subunit in drive carbon flux into starch. Changes in both enzymes (SUSY and AGPse) are then responsible of higher starch concentrations. Low dry matter genotypes showed generally higher vacuolar invertase gene expression (and also enzyme activity), early in fruit development. This alternative cleavage strategy can possibly contribute to energy loss, in that invertases’ products are not adenylated, and further reactions and transport are needed to convert carbon into starch. Although these elements match well with observed differences in starch contents, other factors could be involved in carbon metabolism control. From the microarray experiment, in fact, several kinases and transcription factors have been found to be differentially expressed. Sink strength is known to be modified by application of regulators. In ‘Hayward’ kiwifruit, the synthetic cytokinin CPPU (N-(2-Chloro-4-Pyridyl)-N-Phenylurea) promotes a dramatic increase in fruit size, whereas dry matter content decreases. The behaviour of CPPU-treated ‘Hayward’ kiwifruit was similar to that of fruit from low dry matter genotypes: dry matter and starch concentrations were lower. However, the CPPU effect was strongly source limited, whereas in genotype variation it was not. Moreover, CPPU-treated fruit gene expression (at sucrose cleavage and AGPase levels) was similar to that in high dry matter genotypes. It was therefore concluded that CPPU promotes both sink size and sink activity, but at different “speeds” and this ends in the observed decrease in dry matter content and starch concentration. The lower “speed” in sink activity is probably due to a differential partitioning of activated glucose between starch storage and cell wall synthesis to sustain cell expansion. Starch is the main carbohydrate accumulated in growing Actinidia deliciosa fruit. Results obtained in the present study suggest that sucrose synthase and AGPase enzymes contribute to sucrose to starch conversion, and differences in their gene expression levels, mainly early in fruit development, strongly affect the rate at which starch is therefore accumulated. This results are interesting in that starch and Actinidia deliciosa fruit quality are tightly connected.