Comparative genomic sequencing analysis of a region in human chromosome 11p15.3, mouse chromosome 7 and analysis of the novel STK33, Stk33 gene

The comparative genomic sequence analysis of a region in human chromosome 11p15.3 and its homologous segment in mouse chromosome 7 between ST5 and LMO1 genes has been performed. 158,201 bases were sequenced in the mouse and compared with the syntenic region in human, partially available in the public databases. The analysed region exhibits the typical eukaryotic genomic structure and compared with the close neighbouring regions, strikingly reflexes the mosaic pattern distribution of (G+C) and repeats content despites its relative short size. Within this region the novel gene STK33 was discovered (Stk33 in the mouse), that codes for a serine/threonine kinase. The finding of this gene constitutes an excellent example of the strength of the comparative sequencing approach. Poor gene-predictions in the mouse genomic sequence were corrected and improved by the comparison with the unordered data from the human genomic sequence publicly available. Phylogenetical analysis suggests that STK33 belongs to the calcium/calmodulin-dependent protein kinases group and seems to be a novelty in the chordate lineage. The gene, as a whole, seems to evolve under purifying selection whereas some regions appear to be under strong positive selection. Both human and mouse versions of serine/threonine kinase 33, consists of seventeen exons highly conserved in the coding regions, particularly in those coding for the core protein kinase domain. Also the exon/intron structure in the coding regions of the gene is conserved between human and mouse. The existence and functionality of the gene is supported by the presence of entries in the EST databases and was in vivo fully confirmed by isolating specific transcripts from human uterus total RNA and from several mouse tissues. Strong evidence for alternative splicing was found, which may result in tissue-specific starting points of transcription and in some extent, different protein N-termini. RT-PCR and hybridisation experiments suggest that STK33/Stk33 is differentially expressed in a few tissues and in relative low levels. STK33 has been shown to be reproducibly down-regulated in tumor tissues, particularly in ovarian tumors. RNA in-situ hybridisation experiments using mouse Stk33-specific probes showed expression in dividing cells from lung and germinal epithelium and possibly also in macrophages from kidney and lungs. Preliminary experimentation with antibodies designed in this work, performed in parallel to the preparation of this manuscript, seems to confirm this expression pattern. The fact that the chromosomal region 11p15 in which STK33 is located may be associated with several human diseases including tumor development, suggest further investigation is necessary to establish the role of STK33 in human health.

Bioinformatic methods in applied genomic research

Here I will focus on three main topics that best address and include the projects I have been working in during my three year PhD period that I have spent in different research laboratories addressing both computationally and practically important problems all related to modern molecular genomics. The first topic is the use of livestock species (pigs) as a model of obesity, a complex human dysfunction. My efforts here concern the detection and annotation of Single Nucleotide Polymorphisms. I developed a pipeline for mining human and porcine sequences. Starting from a set of human genes related with obesity the platform returns a list of annotated porcine SNPs extracted from a new set of potential obesity-genes. 565 of these SNPs were analyzed on an Illumina chip to test the involvement in obesity on a population composed by more than 500 pigs. Results will be discussed. All the computational analysis and experiments were done in collaboration with the Biocomputing group and Dr.Luca Fontanesi, respectively, under the direction of prof. Rita Casadio at the Bologna University, Italy. The second topic concerns developing a methodology, based on Factor Analysis, to simultaneously mine information from different levels of biological organization. With specific test cases we develop models of the complexity of the mRNA-miRNA molecular interaction in brain tumors measured indirectly by microarray and quantitative PCR. This work was done under the supervision of Prof. Christine Nardini, at the “CAS-MPG Partner Institute for Computational Biology” of Shangai, China (co-founded by the Max Planck Society and the Chinese Academy of Sciences jointly) The third topic concerns the development of a new method to overcome the variety of PCR technologies routinely adopted to characterize unknown flanking DNA regions of a viral integration locus of the human genome after clinical gene therapy. This new method is entirely based on next generation sequencing and it reduces the time required to detect insertion sites, decreasing the complexity of the procedure. This work was done in collaboration with the group of Dr. Manfred Schmidt at the Nationales Centrum für Tumorerkrankungen (Heidelberg, Germany) supervised by Dr. Annette Deichmann and Dr. Ali Nowrouzi. Furthermore I add as an Appendix the description of a R package for gene network reconstruction that I helped to develop for scientific usage (http://www.bioconductor.org/help/bioc-views/release/bioc/html/BUS.html).

Depicting the role of CDKN2A/ARF alterations in adult BCR-ABL1-positive acute lymphoblastic leukemia patients: from genomic deletions to prognostic impact

This 9p21 locus, encode for important proteins involved in cell cycle regulation and apoptosis containing the p16/CDKN2A (cyclin-dependent kinase inhibitor 2a) tumor suppressor gene and two other related genes, p14/ARF and p15/CDKN2B. This locus, is a major target of inactivation in the pathogenesis of a number of human tumors, both solid and haematologic, and is a frequent site of loss or deletion also in acute lymphoblastic leukemia (ALL) ranging from 18% to 45% 1. In order to explore, at high resolution, the frequency and size of alterations affecting this locus in adult BCR-ABL1-positive ALL and to investigate their prognostic value, 112 patients (101 de novo and 11 relapse cases) were analyzed by genome-wide single nucleotide polymorphisms arrays and gene candidate deep exon sequencing. Paired diagnosis-relapse samples were further available and analyzed for 19 (19%) cases. CDKN2A/ARF and CDKN2B genomic alterations were identified in 29% and 25% of newly diagnosed patients, respectively. Deletions were monoallelic in 72% of cases and in 43% the minimal overlapping region of the lost area spanned only the CDKN2A/2B gene locus. The analysis at the time of relapse showed an almost significant increase in the detection rate of CDKN2A/ARF loss (47%) compared to diagnosis (p = 0.06). Point mutations within the 9p21 locus were found at very low level with only a non-synonymous substition in the exon 2 of CDKN2A. Finally, correlation with clinical outcome showed that deletions of CDKN2A/B are significantly associated with poor outcome in terms of overall survival (p = 0.0206), disease free-survival (p = 0.0010) and cumulative incidence of relapse (p = 0.0014). The inactivation of 9p21 locus by genomic deletions is a frequent event in BCR-ABL1-positive ALL. Deletions are frequently acquired at the leukemia progression and work as a poor prognostic marker.

Functional and Genetic characterization of new genomic islands from an E. coli strain associated with neonatal meningitis.

Enterobacteriaceae genomes evolve through mutations, rearrangements and horizontal gene transfer (HGT). The latter evolutionary pathway works through the acquisition DNA (GEI) modules of foreign origin that enhances fitness of the host to a given environment. The genome of E. coli IHE3034, a strain isolated from a case of neonatal meningitis, has recently been sequenced and its subsequent sequence analysis has predicted 18 possible GEIs, of which: 8 have not been previously described, 5 fully meet the pathogenic island definition and at least 10 that seem to be of prophagic origin. In order to study the GEI distribution of our reference strain, we screened for the presence 18 GEIs a panel of 132 strains, representative of E. coli diversity. Also, using an inverse nested PCR approach we identified 9 GEI that can form an extrachromosomal circular intermediate (CI) and their respective attachment sites (att). Further, we set up a qPCR approach that allowed us to determine the excision rates of 5 genomic islands in different growth conditions. Four islands, specific for strains appertaining to the sequence type complex 95 (STC95), have been deleted in order to assess their function in a Dictyostelium discoideum grazing assays. Overall, the distribution data presented here indicate that 16 IHE3034 GEIs are more associated to the STC95 strains. Also the functional and genetic characterization has uncovered that GEI 13, 17 and 19 are involved in the resistance to phagocitation by Dictyostelium d thus suggesting a possible role in the adaptation of the pathogen during certain stages of infection.

Identification of stock units of Mullus barbatus (Linnaeus, 1758) in the Italian seas: development of genomic 2b-RAD markers and biological variation analysis.

This thesis is settled within the STOCKMAPPING project, which represents one of the studies that were developed in the framework of RITMARE Flagship project. The main goals of STOCKMAPPING were the creation of a genomic mapping for stocks of demersal target species and the assembling of a database of population genomic, in order to identify stocks and stocks boundaries. The thesis focuses on three main objectives representing the core for the initial assessment of the methodologies and structure that would be applied to the entire STOCKMAPPING project: individuation of an analytical design to identify and locate stocks and stocks boundaries of Mullus barbatus, application of a multidisciplinary approach to validate biological methods and an initial assessment and improvement for the genotyping by sequencing technique utilized (2b-RAD). The first step is the individuation of an analytical design that has to take in to account the biological characteristics of red mullet and being representative for STOCKMAPPING commitments. In this framework a reduction and selection steps was needed due to budget reduction. Sampling areas were ranked according the individuation of four priorities. To guarantee a multidisciplinary approach the biological data associated to the collected samples were used to investigate differences between sampling areas and GSAs. Genomic techniques were applied to red mullet for the first time so an initial assessment of molecular protocols for DNA extraction and 2b-RAD processing were needed. At the end 192 good quality DNAs have been extracted and eight samples have been processed with 2b-RAD. Utilizing the software Stacks for sequences analyses a great number of SNPs markers among the eight samples have been identified. Several tests have been performed changing the main parameter of the Stacks pipeline in order to identify the most explicative and functional sets of parameters.

Stock boundaries of Parapenaeus longirostris (Lucas, 1846) populations in the Italian seas: Development of genomic 2b-RAD markers and biological variation analysis.

This thesis is developed in the contest of Ritmare project WP1, which main objective is the development of a sustainable fishery through the identification of populations boundaries in commercially important species in Italian Seas. Three main objectives are discussed in order to help reach the main purpose of identification of stock boundaries in Parapenaeus longirostris: 1 -Development of a representative sampling design for Italian seas; 2 -Evaluation of 2b-RAD protocol; 3 -Investigation of populations through biological data analysis. First of all we defined and accomplished a sampling design which properly represents all Italian seas. Then we used information and data about nursery areas distribution, abundance of populations and importance of P. longirostris in local fishery, to develop an experimental design that prioritize the most important areas to maximize the results with actual project funds. We introduced for the first time the use of 2b-RAD on this species, a genotyping method based on sequencing the uniform fragments produced by type IIB restriction endonucleases. Thanks to this method we were able to move from genetics to the more complex genomics. In order to proceed with 2b-RAD we performed several tests to identify the best DNA extraction kit and protocol and finally we were able to extract 192 high quality DNA extracts ready to be processed. We tested 2b-RAD with five samples and after high-throughput sequencing of libraries we used the software “Stacks” to analyze the sequences. We obtained positive results identifying a great number of SNP markers among the five samples. To guarantee a multidisciplinary approach we used the biological data associated to the collected samples to investigate differences between geographical samples. Such approach assures continuity with other project, for instance STOCKMED, which utilize a combination of molecular and biological analysis as well.

Genomic and proteomic approaches in pig meat quality research field

Pig meat and carcass quality is a complex concept determined by environmental and genetic factors concurring to the phenotypic variation in qualitative characteristics of meat (fat content, tenderness, juiciness, flavor,etc). This thesis shows the results of different investigations to study and to analyze pig meat and carcass quality focusing mainly on genomic; moreover proteomic approach has been also used. The aim was to analyze data from association studies between genes considered as candidate and meat and carcass quality in different pig breeds. The approach was used to detect new SNP in genes functionally associated to the studied traits and to confirm as candidate other genes already known. Five polymorphisms (one new SNP in Calponin 1 gene and four additional polymorphism already known in other genes) were considered on chromosome 2 (SSC2). Calponin 1 (CNN1) was associated to the studied traits and furthermore the results reported confirmed the data already known for Lactate dehydrogenase A (LDHA), Low density lipoprotein receptor (LDLR), Myogenic differentiation 1 (MYOD1) e Ubiquitin-like 5 (UBL5), in Italian Large White pigs. Using an in silico search it was possible to detect on SSC2 a new SNP of Deoxyhypusine synthase (DHPS) gene partially overlapping with WD repeat domain 83 (WDR83) gene and significant for the meat pH variation in Italian Large White (ILW) pigs. Perilipin 1 (PLIN1) mapping on chromosome 7 and Perilipin 2 (PLIN2) mapping on chromosome 1 were studied and the results obtained in Duroc breed have shown significant associations with carcass traits. Moreover a study of protein composition of porcine LD muscle, indicated an effect of temperature treatment of carcass, on proteins of the sarcoplasmic fraction and in particular on PGM1 phosphorylation. Future studies on pig meat quality should be based on the integration of different experimental approaches (genomics, proteomics, transcriptomics, etc).

Identifizierung, Quantifizierung und Hemmung von ausgewählten Hefen im Wein

Hefen stellen einen großen und wichtigen Teil der Mikrobiota während der Weinbereitung dar, da ohne ihre alkoholische Fermentation die Umwandlung von Most und Wein nicht möglich wäre. Ferner ist es ihre Vielzahl an Stoffwechselprodukten, die dem Aroma des fertigen Weines eine zusätzliche Komplexität verleihen. Auf der anderen Seite steht durch den Metabolismus verschiedenster so genannter Wildhefen die Gefahr von Qualitätsabstufungen der Weine, was allgemein als „Weinfehler“ betrachtet wird. Ziel dieser Arbeit war zum einen die taxonomische Einordnung von Saccharomyces-Spezies, sowie die Quantifizierung und Hemmung von ausgewählten Wildhefen während der Weinbereitung.rnEin Teil dieser Arbeit umfasste die Identifizierung der nahverwandten Mitglieder der Saccharomyces sensu stricto-Gruppe. Durch den Einsatz des DNA-Fingerpinting-Systems SAPD-PCR konnten alle die Gruppe umfassenden Spezies anhand spezifischer Bandenmuster nachgewiesen werden, wodurch eine Einordnung dieser schwer zu differenzierenden Arten möglich war. Die Differenzierung zwischen den einzelnen Spezies war in jedem Fall deutlicher als dies die Sequenzierung der 5.8S rDNA und ihre flankierenden ITS-Regionen vermochte. Die SAPD-PCR zeichnete sich zudem durch eine geringe Muster-Varianz bei verschiedenen Stämmen einer Art aus und konnte zuverlässig unbekannte Stämme bestimmen und bereits hinterlegte Stämme neu klassifizieren. Zudem konnte mit Hilfe dieses Systems Hybride aus Saccharomyces cerevisiae und S. bayanus bzw. S. cerevisiae und S. kudriavzevii detektiert werden, wenn diese Hybride aus relativ gleichen genomischen Anteilen der Eltern bestanden. rnZusätzlich wurde ein quantitatives PCR-System entwickelt, um die Gattungen Saccharomyces, Hanseniaspora und Brettanomyces in Most und Wein detektieren und quantifizieren zu können. Die hierfür entwickelten Primer zeigten sich spezifisch für die untersuchten Arten. Durch die serielle Verdünnung definierter DNA-Mengen konnte für alle drei Systeme eine Kalibrierungskurve erstellt werden, mit Hilfe derer die tatsächlichen Quantifizierungen durchgeführt wurden. Die qPCR-Analyse lieferte ähnliche Zellzahlen wie Lebendzellzahl-Bestimmungen und wurde nicht von anderen Spezies und von Traubensaft gestört. Die maximal detektierbare Zellzahl betrug 2 x 107 Zellen/ml, während die minimale Detektionsgrenze je nach Art zwischen 1 x 102 Zellen/ml und 1 x 103 Zellen/ml lag. Allerdings konnte eine effektive DNA-Isolierung dieser geringen Zellzahlen nur erreicht werden, wenn die Zellzahl durch artfremde Hefen künstlich erhöht wurde. Die Analyse einer Most-Vergärung mit den drei Spezies zeigte schlussendlich, dass die quantitative PCR sicher und schnell Veränderungen und Sukzessionen detektiert und so ein geeignetes Mittel darstellt, um Populationsdynamiken während der Weinherstellung zu beobachten. rnDer letzte Teil dieser Arbeit befasste sich mit der Inhibierung von Schadhefen durch zellwand-hydrolysierende Enzyme. Es konnte hierbei eine endoglykosidisch wirkende β-1,3-Glucanase aus dem Bakterium Delftia tsuruhatensis isoliert werden. Diese besaß eine ungefähre Masse von 28 kDa, einen isolektrischen Punkt von ca. 4,3 und wirkte mit einer spezifischen Aktivität von 10 U/mg Protein gegen das Glucan Laminarin. Zudem zeigte das Enzym ein Temperaturoptimum von 50 °C und ein pH-Optimum bei pH 4,0. Weinparameter wie erhöhte Konzentrationen an Ethanol, Phenolen und Sulfit beeinflussten die Wirkung des Enzyms nicht oder nur wenig. Neben der allgemeinen Wirkung gegen β-1,3-Glucane konnte hier auch gezeigt werden, dass ebenso gut die β-1,3-Glucane in der Zellwand verschiedener Hefen hydrolysiert wurden. Fluoreszenz- und rasterelektronen-mikroskopische Aufnahmen von Hefezellen nach Inkubation mit der β-1,3-Glucanase zeigten zusätzlich die Zerstörung der Zelloberfläche der Hefen. Die lytische Wirkung des Enzyms wurde an verschiedenen weintypischen Hefen getestet. Hierbei zeigten sich stammspezifische Unterschiede in der Sensitivität gegenüber dem Enzym. Außerdem konnte festgestellt werden, dass sowohl Wachstumsphase als auch Medium der Hefen Einfluss auf deren Zellwand hat und somit auch auf die Wirkung des Enzyms.rn

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.

Studio e valutazione sperimentale di tecniche di radio fingerprinting per navigazione indoor con dispositivi mobili Android

Lo studio effettuato raccoglie informazioni al fine di svolgere un’analisi globale delle tecniche attualmente in uso, o in fase di studio, nel campo della localizzazione di dispositivi all’interno di un ambiente chiuso, ovvero laddove non è possibile sfruttare la copertura del sistema GPS. La panoramica è frutto dell’analisi e dello studio di paper tecnici pubblicati dai collaboratori dell’IEEE, fruibili all’interno del portale IEEE Xplore. A corredo di questo studio è stata sviluppata una applicazione per dispositivi Android basata sulla tecnica del Wi-Fi fingerprint; l’applicazione, che rappresenta un primo approccio alle tecniche di localizzazione, è a tutti gli effetti un sistema standalone per la localizzazione, consente cioè di costruire sia la mappa per la localizzazione, partendo da un ambiente sconosciuto, sia di ottenere la posizione dell’utente all’interno di una mappa conosciuta. La tesi si conclude con una analisi dei risultati e delle performance ottenute dall’applicazione in un uso comune, al fine di poter valutare l’efficacia della tecnica presa in considerazione. I possibili sviluppi futuri sono analizzati in un capitolo a parte e trovano spazio in ambienti nei quali si vogliono offrire servizi "context-based", ovvero basati sulla posizione esatta dell’utente.

Progettazione e implementazione di una piattaforma per la localizzazione indoor basata su tecniche di radio fingerprinting e data-crowdsourcing

Questa tesi ha come scopo principale l'analisi delle diverse tecnologie di localizzazione in ambito indoor, analizzando in particolare l'utilizzo del Wifi RSS Fingerprinting. La tecnica del Wifi RSS Fingerprinting è una tecnica per la localizzazione all'interno di ambienti chiusi, che consiste nella definizione di un 'impronta'(fingerprint) in un punto preciso dell'ambiente(definito reference point), andando a inserire in un database i valori di potenza del segnale ricevuto(RSS) da ogni access point rilevato all'interno di quel determinato reference point. Per l'implementazione di questa tecnica è stato sviluppato un applicativo con un architettura client-server. Il client è stato sviluppato in ambiente Android, realizzando una applicazione per la gestione della fase di salvataggio di nuovi fingerprint e per la fase di localizzazione della posizione corrente, tramite l'utilizzo dei vari fingerprint precedentemente inseriti all'interno del DB. Il server, sviluppato in Node.js(framework Javascript), gestirà le diverse richieste ricevute dal client tramite delle chiamate AJAX, prelevando le informazioni richieste direttamente dal database. All'interno delle applicativo sono stati implementati diversi algoritmi per la localizzazione indoor, in modo da poter verificare l'applicabilità di questo sistema in un ambito reale. Questi algoritmi sono stati in seguito testati per valutare l'accuratezza e la precisione di ciascuno, andando ad individuare gli algoritmi migliori da utilizzare in base a scenari diversi.

Quantitative 1-Step DNA Methylation Analysis with Native Genomic DNA as Template

DNA methylation analysis currently requires complex multistep procedures based on bisulfite conversion of unmethylated cytosines or on methylation-sensitive endonucleases. To facilitate DNA methylation analysis, we have developed a quantitative 1-step assay for DNA methylation analysis.

Deep-sequencing identification of the genomic targets of the cytidine deaminase AID and its cofactor RPA in B lymphocytes

The cytidine deaminase AID hypermutates immunoglobulin genes but can also target oncogenes, leading to tumorigenesis. The extent of AID's promiscuity and its predilection for immunoglobulin genes are unknown. We report here that AID interacted broadly with promoter-proximal sequences associated with stalled polymerases and chromatin-activating marks. In contrast, genomic occupancy of replication protein A (RPA), an AID cofactor, was restricted to immunoglobulin genes. The recruitment of RPA to the immunoglobulin loci was facilitated by phosphorylation of AID at Ser38 and Thr140. We propose that stalled polymerases recruit AID, thereby resulting in low frequencies of hypermutation across the B cell genome. Efficient hypermutation and switch recombination required AID phosphorylation and correlated with recruitment of RPA. Our findings provide a rationale for the oncogenic role of AID in B cell malignancy.

No evidence for large differences in genomic methylation between wild and hatchery steelhead (Oncorhynchus mykiss)

When salmonid fish that have been raised in hatcheries spawn in the wild, they often produce fewer surviving adult offspring than wild fish. Recent data from steelhead (Oncorhynchus mykiss) in the Hood River (Oregon, USA) show that even one or two generations of hatchery culture can result in dramatic declines in fitness. Although intense domestication selection could cause such declines, it is worth considering alternative explanations. One possibility is heritable epigenetic changes induced by the hatchery environment. Here, we show, using methylation-sensitive amplified fragment length polymorphism, that hatchery and wild adult steelhead from the Hood River do not appear to differ substantially in overall levels of genomic methylation. Thus, although altered methylation of specific DNA sites or other epigenetic processes could still be important, the hatchery environment does not appear to cause a global hypo- or hypermethylation of the genome or create a large number of sites that are differentially methylated.