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.

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.

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

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

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

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

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.

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.

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.

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.

Identification of molecular biomarkers in esophageal adenocarcinoma

Esophageal adenocarcinoma (EAC) is a severe cancer that has been on the rise in Western nations over the past few decades. It has a high mortality rate and the 5-year survival rate is only 35%–45%. EAC has been included in a group of tumors with one of the highest rates of copy number alterations (CNAs), somatic structural rearrangements, high mutation frequency, with different mutational signatures, and with epigenetic mechanisms. The vast heterogeneity of EAC mutations makes it challenging to comprehend the biology that underlies tumor onset and development, identify prognostic biomarkers, and define a molecular classification to stratify patients. The only way to resolve the current disagreements is through an exhaustive molecular analysis of EAC. We examined the genetic profile of 164 patients' esophageal adenocarcinoma samples (without chemo-radiotherapy). The included patients did not receive neoadjuvant therapies, which can change the genetic and molecular composition of the tumor. Using next-generation sequencing technologies (NGS) at high coverage, we examined a custom panel of 26 cancer-related genes. Over the entire cohort, 337 variants were found, with the TP53 gene showing the most frequent alteration (67.27%). Poorer cancer-specific survival was associated with missense mutations in the TP53 gene (Log Rank P=0.0197). We discovered HNF1alpha gene disruptive mutations in 7 cases that were also affected by other gene changes. We started to investigate its role in EAC cell lines by silencing HNF1alpha to mimic our EAC cohort and we use Seahorse technique to analyze its role in the metabolism in esophageal cell. No significant changes were found in transfected cell lines. We conclude by finding that a particular class of TP53 mutations (missense changes) adversely impacted cancer-specific survival in EAC. HNF1alpha, a new EAC-mutated gene, was found, but more research is required to fully understand its function as a tumor suppressor gene.

Phenotypic and molecular investigation of circulating tumor cells (CTCs) isolated from breast cancer patients at single cell resolution

Despite the paramount advances in cancer research, breast cancer (BC) still ranks one of the leading causes of cancer-related death worldwide. Thanks to the screening campaign started in developed countries, BC is often diagnosed at early stages (non-metastatic BC, nmBC), but disease relapse occurrence even after decades and at distant sites is not an uncommon phenomenon. Conversely, metastatic BC (mBC) is considered an incurable disease. The major perpetrators of tumor spread to secondary organs are circulating tumor cells (CTCs), a rare population of cells detectable in the peripheral blood of oncologic patients. In this study, CTCs from patients diagnosed with luminal nmBC and mBC (hormone receptor positive, Human Epidermal Growth Factor Receptor 2 (HER2) negative) were characterized at both phenotypic and molecular levels. To better understand the molecular mechanisms underlying their biology and their metastatic potential, next-generation sequencing (NGS) analyses were performed at single-cell resolution to assess copy number aberrations (CNAs), single nucleotide variants (SNVs) and gene expression profiling. The findings of this study arise hints in CTC detection, and pave the way to new application in CTC research.

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.

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.