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).

Natural compounds Camptothecin and Triptolide: highly specific enzyme inhibitors and tools to dissect transcriptional functions

With this work I elucidated new and unexpected mechanisms of two strong and highly specific transcription inhibitors: Triptolide and Campthotecin. Triptolide (TPL) is a diterpene epoxide derived from the Chinese plant Trypterigium Wilfoordii Hook F. TPL inhibits the ATPase activity of XPB, a subunit of the general transcription factor TFIIH. In this thesis I found that degradation of Rbp1 (the largest subunit of RNA Polymerase II) caused by TPL treatments, is preceded by an hyperphosphorylation event at serine 5 of the carboxy-terminal domain (CTD) of Rbp1. This event is concomitant with a block of RNA Polymerase II at promoters of active genes. The enzyme responsible for Ser5 hyperphosphorylation event is CDK7. Notably, CDK7 downregulation rescued both Ser5 hyperphosphorylation and Rbp1 degradation triggered by TPL. Camptothecin (CPT), derived from the plant Camptotheca acuminata, specifically inhibits topoisomerase 1 (Top1). We first found that CPT induced antisense transcription at divergent CpG islands promoter. Interestingly, by immunofluorescence experiments, CPT was found to induce a burst of R loop structures (DNA/RNA hybrids) at nucleoli and mitochondria. We then decided to investigate the role of Top1 in R loop homeostasis through a short interfering RNA approach (RNAi). Using DNA/RNA immunoprecipitation techniques coupled to NGS I found that Top1 depletion induces an increase of R loops at a genome-wide level. We found that such increase occurs on the entire gene body. At a subset of loci R loops resulted particularly stressed after Top1 depletion: some of these genes showed the formation of new R loops structures, whereas other loci showed a reduction of R loops. Interestingly we found that new peaks usually appear at tandem or divergent genes in the entire gene body, while losses of R loop peaks seems to be a feature specific of 3’ end regions of convergent genes.

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

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

Statistical Physics and Modeling of Human Mobility

In this thesis, we extend some ideas of statistical physics to describe the properties of human mobility. By using a database containing GPS measures of individual paths (position, velocity and covered space at a spatial scale of 2 Km or a time scale of 30 sec), which includes the 2% of the private vehicles in Italy, we succeed in determining some statistical empirical laws pointing out "universal" characteristics of human mobility. Developing simple stochastic models suggesting possible explanations of the empirical observations, we are able to indicate what are the key quantities and cognitive features that are ruling individuals' mobility. To understand the features of individual dynamics, we have studied different aspects of urban mobility from a physical point of view. We discuss the implications of the Benford's law emerging from the distribution of times elapsed between successive trips. We observe how the daily travel-time budget is related with many aspects of the urban environment, and describe how the daily mobility budget is then spent. We link the scaling properties of individual mobility networks to the inhomogeneous average durations of the activities that are performed, and those of the networks describing people's common use of space with the fractional dimension of the urban territory. We study entropy measures of individual mobility patterns, showing that they carry almost the same information of the related mobility networks, but are also influenced by a hierarchy among the activities performed. We discover that Wardrop's principles are violated as drivers have only incomplete information on traffic state and therefore rely on knowledge on the average travel-times. We propose an assimilation model to solve the intrinsic scattering of GPS data on the street network, permitting the real-time reconstruction of traffic state at a urban scale.

Design and Synthesis of Naphthalene Diimides Based Small Molecules as Anticancer Agents: Targeting the Polyamine Transporter, G-quadruplex Structures and HDAC

Cancer is a multifactorial disease characterized by a very complex etiology. Basing on its complex nature, a promising therapeutic strategy could be based by the “Multi-Target-Directed Ligand” (MTDL) approach, based on the assumption that a single molecule could hit several targets responsible for the pathology. Several agents acting on DNA are clinically used, but the severe deriving side effects limit their therapeutic application. G-quadruplex structures are DNA secondary structures located in key zones of human genome; targeting quadruplex structures could allow obtaining an anticancer therapy more free from side effects. In the last years it has been proved that epigenetic modulation can control the expression of human genes, playing a crucial role in carcinogenesis and, in particular, an abnormal expression of histone deacetylase enzymes are related to tumor onset and progression. This thesis deals with the design and synthesis of new naphthalene diimide (NDI) derivatives endowed with anticancer activity, interacting with DNA together with other targets implicated in cancer development, such as HDACs. NDI-polyamine and NDI-polyamine-hydroxamic acid conjugates have been designed with the aim to provide potential MTDLs, in order to create molecules able simultaneously to interact with different targets involved in this pathology, specifically the G-quadruplex structures and HDAC, and to exploit the polyamine transport system to get selectively into cancer cells. Macrocyclic NDIs have been designed with the aim to improve the quadruplex targeting profile of the disubstituted NDIs. These compounds proved the ability to induce a high and selective stabilization of the quadruplex structures, together with cytotoxic activities in the micromolar range. Finally, trisubstituted NDIs have been developed as G-quadruplex-binders, potentially effective against pancreatic adenocarcinoma. In conclusion, all these studies may represent a promising starting point for the development of new interesting molecules useful for the treatment of cancer, underlining the versatility of the NDI scaffold.

Involvement of microRNAs in Androgen Receptor-dependent Breast Cancers

Triple negative breast cancer (TNBC) is a very aggressive tumor subtype characterized by the lack of expression of estrogen receptor 1 (ESR1), due in the most of cases to an increased expression of DNA methyltransferases (DNMTs) and hypermethylation in CpG islands, resulting in gene silencing. Furthermore, in ESR1- negative breast cancers, androgen receptor (AR) is highly expressed and some studies suggest that it can drive tumor progression and might represent a therapeutic target. A correlation between microRNAs, small non-coding RNAs that regulate gene expression, and DNMTs was investigated in a TNBC cell line to restore a normal methylation pattern of ESR1, leading to its re-expression and conferring again sensitivity to selective estrogen receptor modulators (SERMs). miR-148A and miR-29B were found to be involved in the reduction of the expression of DNMT1 and DNMT3A and in a slight increase of ESR1 expression, but not at protein level. Then, we found a down-regulation of AR by miRs-7, -9, -27a, -27b, -29a, -29b, -29c, -127-3p, -127-5p and -376 at 48h post transfection and an up-regulation by miR-15a and miR-16 at every time considered. We concomitantly investigated a possible increase of Tamoxifen, Herceptin and Metformin sensitivity after AR silencing in MDA-MB 453 and T-47D cell lines. Cells seemed more sensitive when silenced for AR only in MDA-MB-453 at 24h post Tamoxifen treatment. Studies on Metformin have basically confirmed an increase of drug sensitivity due to AR silencing in both cell lines. Analysis of Herceptin showed how MDA-MB 453 samples silenced for AR have a slight decrease in the percentage of proliferating cells, demonstrating a possible increase in the response to treatment. These preliminary data provide the basis for further study of the modulation of the expression of AR by microRNAs and it will be interesting to understand the molecular mechanisms underlying these interactions.

Glycomics as an innovative technology to identify biomarkers of aging

Introduction. Glycomic analysis allows investigating on the global glycome within body fluids (as serum/plasma), this could eventually lead to identify new types of disease biomarkers, or as in this study, biomarkers of human aging studying specific aging models. Recent studies demonstrated that the plasma N-glycome is modified during human aging, suggesting that measurements of log-ratio of two serum/plasma N-glycans (NGA2F and NA2F), named GlycoAge test could provide a non-invasive biomarker of aging. Down syndrome (DS) is a genetic disorder in which multiple major aspects of senescent phenotype occur much earlier than in healthy age-matched subjects and has been often defined as an accelerated aging syndrome. The aim of this study was to compare plasma N-glycome of patients affected by DS with age- and sex matched non-affected controls, represented by their siblings (DSS), in order to assess if DS is characterized by a specific N-glycomic pattern. Therefore, in order to investigate if N-glycans changes that occur in DS were able to reveal an accelerated aging in DS patients, we enrolled the mothers (DSM) of the DS and DSS, representing the non-affected control group with a different chronological age respect to DS. We applied two different N-glycomics approaches on the same samples: first, in order to study the complete plasma N-glycome we applied a new high-sensitive protocol based on a MALDI-TOF-MS approach, second, we used DSA-FACE technology. Results: MALDI-TOF/MS analysis detected a specific N-glycomics signature for DS, characterized by an increase of fucosylated and bisecting species. Moreover, in DS the abundance of agalactosylated (as NA2F) species was similar or higher than their mothers. The measurement of GlycoAge test with DSA-FACE, validated also by MALDI-TOF, demonstrated a strongly association with age, moreover in DS, it’s value was similar to their mothers, and significantly higher than their age- and sex matched not-affected siblings

Harnessing the power of Drosophila model to obtain new insights on MYC and ABCC1 cellular functions: from neuroblastoma to autism spectrum disorder

MYCN amplification is a genetic hallmark of the childhood tumour neuroblastoma. MYCN-MAX dimers activate the expression of genes promoting cell proliferation. Moreover, MYCN seems to transcriptionally repress cell differentiation even in absence of MAX. We adopted the Drosophila eye as model to investigate the effect of high MYC to MAX expression ratio on cells. We found that dMyc overexpression in eye cell precursors inhibits cell differentiation and induces the ectopic expression of Antennapedia (the wing Hox gene). The further increase of MYC/MAX ratio results in an eye-to-wing homeotic transformation. Notably, dMyc overexpression phenotype is suppressed by low levels of transcriptional co-repressors and MYCN associates to the promoter of Deformed (the eye Hox gene) in proximity to repressive sites. Hence, we envisage that, in presence of high MYC/MAX ratio, the “free MYC” might inhibit Deformed expression, leading in turn to the ectopic expression of Antennapedia. This suggests that MYCN might reinforce its oncogenic role by affecting the physiological homeotic program. Furthermore, poor neuroblastoma outcome associates with a high level of the MRP1 protein, encoded by the ABCC1 gene and known to promote drug efflux in cancer cells. Intriguingly, this correlation persists regardless of chemotherapy and ABCC1 overexpression enhances neuroblastoma cell motility. We found that Drosophila dMRP contributes to the adhesion between the dorsal and ventral epithelia of the wing by inhibiting the function of integrin receptors, well known regulators of cell adhesion and migration. Besides, integrins play a crucial role during synaptogenesis and ABCC1 locus is included in a copy number variable region of the human genome (16p13.11) involved in neuropsychiatric diseases. Interestingly, we found that the altered dMRP/MRP1 level affects nervous system development in Drosophila embryos. These preliminary findings point out novel ABCC1 functions possibly defining ABCC1 contribution to neuroblastoma and to the pathogenicity of 16p13.11 deletion/duplication

EBV Type II latency relies on PARP1 activity and is essential for gastric epithelial cell transformation in EBV-associated Gastric Cancer (EBVaGC)

Epstein-Barr virus (EBV) establishes a lifelong asymptomatic infection by replicating its chromatinized genome, called episome, together with the host genome. EBV exhibits different latency-associated transcriptional repertoires that mirror its three-dimensional structures of the genome. CTCF, Cohesin and PARP1 are involved in maintaining viral latency and establishing episome architecture. Epstein-Barr virus-associated gastric cancer (EBVaGC) represents almost 10% of all gastric cancers globally. EBVaGC exhibit an intermediate viral transcription profile known as "Latency II", expressing specific viral genes and non-coding RNAs. In this study, we investigated the impact of PARP1 inhibition on CTCF/Cohesin binding in Type II latency. We observed a destabilization of the binding of both factors, leading to a disrupted three-dimensional architecture of the episomes and consequently, an altered viral gene expression. Despite sharing the same CTCF binding profile, Type I, II, and III latencies display different 3D episomal structures that correlate with variations in viral gene expression. Additionally, our analysis of H3K27ac-enriched chromatin interactions revealed differences between Type II latency episomes and a link to cellular transformation through docking of the EBV episomes at specific sites of the Human genome, thus promoting oncogene expression. Overall, this work provides insights into the role of PARP1 in maintaining active latency and novel mechanisms of EBV-induced cellular transformation.

Transcriptional functions of DNA Topoisomerases at a genome-wide scale and a single gene levels

The DNA topology is an important modifier of DNA functions. Torsional stress is generated when right handed DNA is either over- or underwound, producing structural deformations which drive or are driven by processes such as replication, transcription, recombination and repair. DNA topoisomerases are molecular machines that regulate the topological state of the DNA in the cell. These enzymes accomplish this task by either passing one strand of the DNA through a break in the opposing strand or by passing a region of the duplex from the same or a different molecule through a double-stranded cut generated in the DNA. Because of their ability to cut one or two strands of DNA they are also target for some of the most successful anticancer drugs used in standard combination therapies of human cancers. An effective anticancer drug is Camptothecin (CPT) that specifically targets DNA topoisomerase 1 (TOP 1). The research project of the present thesis has been focused on the role of human TOP 1 during transcription and on the transcriptional consequences associated with TOP 1 inhibition by CPT in human cell lines. Previous findings demonstrate that TOP 1 inhibition by CPT perturbs RNA polymerase (RNAP II) density at promoters and along transcribed genes suggesting an involvement of TOP 1 in RNAP II promoter proximal pausing site. Within the transcription cycle, promoter pausing is a fundamental step the importance of which has been well established as a means of coupling elongation to RNA maturation. By measuring nascent RNA transcripts bound to chromatin, we demonstrated that TOP 1 inhibition by CPT can enhance RNAP II escape from promoter proximal pausing site of the human Hypoxia Inducible Factor 1 (HIF-1) and c-MYC genes in a dose dependent manner. This effect is dependent from Cdk7/Cdk9 activities since it can be reversed by the kinases inhibitor DRB. Since CPT affects RNAP II by promoting the hyperphosphorylation of its Rpb1 subunit the findings suggest that TOP 1inhibition by CPT may increase the activity of Cdks which in turn phosphorylate the Rpb1 subunit of RNAP II enhancing its escape from pausing. Interestingly, the transcriptional consequences of CPT induced topological stress are wider than expected. CPT increased co-transcriptional splicing of exon1 and 2 and markedly affected alternative splicing at exon 11. Surprisingly despite its well-established transcription inhibitory activity, CPT can trigger the production of a novel long RNA (5’aHIF-1) antisense to the human HIF-1 mRNA and a known antisense RNA at the 3’ end of the gene, while decreasing mRNA levels. The effects require TOP 1 and are independent from CPT induced DNA damage. Thus, when the supercoiling imbalance promoted by CPT occurs at promoter, it may trigger deregulation of the RNAP II pausing, increased chromatin accessibility and activation/derepression of antisense transcripts in a Cdks dependent manner. A changed balance of antisense transcripts and mRNAs may regulate the activity of HIF-1 and contribute to the control of tumor progression After focusing our TOP 1 investigations at a single gene level, we have extended the study to the whole genome by developing the “Topo-Seq” approach which generates a map of genome-wide distribution of sites of TOP 1 activity sites in human cells. The preliminary data revealed that TOP 1 preferentially localizes at intragenic regions and in particular at 5’ and 3’ ends of genes. Surprisingly upon TOP 1 downregulation, which impairs protein expression by 80%, TOP 1 molecules are mostly localized around 3’ ends of genes, thus suggesting that its activity is essential at these regions and can be compensate at 5’ ends. The developed procedure is a pioneer tool for the detection of TOP 1 cleavage sites across the genome and can open the way to further investigations of the enzyme roles in different nuclear processes.

Human papillomavirus (HPV) and associated diseases: between applied diagnostic and basic research

Human Papillomavirus (HPV) is the cause of cervical cancers (among these, adenocarcinoma, AdCa) and is associated to a subgroup of oropharyngeal carcinomas (OPSCCs). Even if the risk for cancer development is linked to the infection by some viral genotypes, mainly HPV16 and 18, viral DNA alone seems not to be sufficient for diagnosis. Moreover, the role of the virus in OPSCCs has not been totally clarified yet. In the first part of the thesis, the performances concerning viral genotyping in clinical cervical samples of a new pyrosequencing-based test and a well-known hybridization-based assay have been compared. Similar results between the methods have been obtained. However, the former showed advantages in detecting intratype variants, higher specificity and a broader spectrum of detectable HPV types. The second part deals with the evaluation of virological markers (genotyping, viral oncoproteins expression, viral load, physical state and CpG methylation of HPV16 genome) in the diagnosis/prognosis of cervical AdCa and HPV-associated OPSCCs. HPV16 has been confirmed the most prevalent genotype in both the populations. Interestingly, the mean methylation frequency of viral DNA at the early promoter showed the tendency to be associated to invasion for cervical AdCa and to a worse prognosis for OPSCCs, suggesting a promising role as diagnostic/prognostic biomarker. The experiments of the third part were performed at the DKFZ in Heidelberg (Germany) and dealt with the analysis of the response to IFN-k transfection in HPV16-positive cervical cancer and head&neck carcinoma cell lines to evaluate its potential role as new treatment. After 24h, we observed increased IFN-b expression which lead to the up-regulation of genes involved in the antigens presentation pathway (MHC class I and immunoproteasome) and antiviral response as well, in particular in cervical cancer cell lines. This fact suggested also the presence of different HPV-mediated carcinogenic pathways between the two anatomical districts.