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

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.

Ecological genetics and conservation genomics of wolf (Canis lupus)

In the present work, we apply both traditional and Next Generation Sequencing (NGS) tools to investigate some of the most important adaptive traits of wolves (Canis lupus). In the first part, we analyze the variability of three Major Histocompatibility Complex (MHC) class II genes in the Italian wolf population, also studying their possible role in mating choice and their influence on fitness traits. In the second section, as part of a larger canid genome project, we will exploit NGS data to investigate the transcript-level differences between the wolf and the dog genome that can be correlated to domestication.

Molecular and functional characterization of the chemotactic genes in the PCBs-degrader Pseudomonas pseudoalcaligenes KF707

Although bacteria represent the simplest form of life on Earth, they have a great impact on all living beings. For example the degrader bacterium Pseudomonas pseudoalcaligenes KF707 is used in bioremediation procedures for the recovery of polluted sites. Indeed, KF707 strain is know for its ability to degrade biphenyl and polychlorinated biphenyls - to which is chemotactically attracted - and to tolerate the oxydative stress due to toxic metal oxyanions such as tellurite and selenite. Moreover, in bioremediation processes, target compounds can be easily accessible to KF707 through biofilm formation. All these considerations suggest that KF707 is such a unique microorganism and this Thesis work has been focused on determining the molecular nature of some of the peculiar physiological traits of this strain. The genome project provided a large set of informations: putative genes involved in the degradation of aromatic and toxic compounds and associated to stress response were identified. Notably, multiple chemotactic operons and cheA genes were also found. Deleted mutants in the cheA genes were constructed and their role in motility, chemotaxis and biofilm formation were assessed and compared to those previously attributed to a cheA1 gene in a KF707 mutant constructed by a mini-Tn5 transposon insertion and which was impaired in motility and biofilm development. The results of this present Thesis work, taken together, were interpreted to suggest that in Pseudomonas pseudoalcaligenes KF707 strain, multiple factors are involved in these networks and they might play different roles depending on the environmental conditions. The ability of KF707 strain to produce signal molecules possibly involved in cell-to-cell communication, was also investigated: lack of a lux-like QS system - which is conversely widely present in Gram negative bacteria – keeps open the question about the actual molecular nature of KF707 quorum sensing mechanism.

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.

Advanced neuroimaging methodologies to improve connectivity detection in normal and abnormal language brain networks

The language connectome was in-vivo investigated using multimodal non-invasive quantitative MRI. In PPA patients (n=18) recruited by the IRCCS ISNB, Bologna, cortical thickness measures showed a predominant reduction on the left hemisphere (p<0.005) with respect to matched healthy controls (HC) (n=18), and an accuracy of 86.1% in discrimination from Alzheimer’s disease patients (n=18). The left temporal and para-hippocampal gyri significantly correlated (p<0.01) with language fluency. In PPA patients (n=31) recruited by the Northwestern University Chicago, DTI measures were longitudinally evaluated (2-years follow-up) under the supervision of Prof. M. Catani, King’s College London. Significant differences with matched HC (n=27) were found, tract-localized at baseline and widespread in the follow-up. Language assessment scores correlated with arcuate (AF) and uncinate (UF) fasciculi DTI measures. In left-ischemic stroke patients (n=16) recruited by the NatBrainLab, King’s College London, language recovery was longitudinally evaluated (6-months follow-up). Using arterial spin labelling imaging a significant correlation (p<0.01) between language recovery and cerebral blood flow asymmetry, was found in the middle cerebral artery perfusion, towards the right. In HC (n=29) recruited by the DIBINEM Functional MR Unit, University of Bologna, an along-tract algorithm was developed suitable for different tractography methods, using the Laplacian operator. A higher left superior temporal gyrus and precentral operculum AF connectivity was found (Talozzi L et al., 2018), and lateralized UF projections towards the left dorsal orbital cortex. In HC (n=50) recruited in the Human Connectome Project, a new tractography-driven approach was developed for left association fibres, using a principal component analysis. The first component discriminated cortical areas typically connected by the AF, suggesting a good discrimination of cortical areas sharing a similar connectivity pattern. The evaluation of morphological, microstructural and metabolic measures could be used as in-vivo biomarkers to monitor language impairment related to neurodegeneration or as surrogate of cognitive rehabilitation/interventional treatment efficacy.

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.

Biological effects of bioactive components and extracts derived from edible plants commonly used in human nutrition

The main aim of this PhD research project was the evaluation of the biological effects of bioactive compounds derived from edible plants, with particular attention on their possibility to counteract oxidative damage and inflammation. After a preliminary study of in vitro antioxidant activity, regarding the modification eventually occurring after home freezing and cooking of edible vegetables, cultured mammalian cells were used as experimental model systems. Soluble extract and essential oils derived from different cultivars of Brassicaceae and Lamiaceae were tested as possible tools for the counteraction of the oxidative damage due to reactive oxygen species (ROS), underlining differences related to cultivar and agronomic techniques. Since accumulating evidence indicates that phytochemicals exhibit several additional properties in complex biological systems, a nutrigenomic approach was used to further explain the biological activity of a green tea extract, and to evidence the anti-inflammatory role of bioactive compounds derived from different foods. Overall, results obtained could contribute to a better understanding of the potential health benefit of plant foods.

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

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

Trascriptome analysis and functional genomics of Neisseria meningitidis in human blood

Neisseria meningitidis (Nm) is the major cause of septicemia and meningococcal meningitis. During the course of infection, it must adapt to different host environments as a crucial factor for survival. Despite the severity of meningococcal sepsis, little is known about how Nm adapts to permit survival and growth in human blood. A previous time-course transcriptome analysis, using an ex vivo model of human whole blood infection, showed that Nm alters the expression of nearly 30% of ORFs of the genome: major dynamic changes were observed in the expression of transcriptional regulators, transport and binding proteins, energy metabolism, and surface-exposed virulence factors. Starting from these data, mutagenesis studies of a subset of up-regulated genes were performed and the mutants were tested for the ability to survive in human whole blood; Nm mutant strains lacking the genes encoding NMB1483, NalP, Mip, NspA, Fur, TbpB, and LctP were sensitive to killing by human blood. Then, the analysis was extended to the whole Nm transcriptome in human blood, using a customized 60-mer oligonucleotide tiling microarray. The application of specifically developed software combined with this new tiling array allowed the identification of different types of regulated transcripts: small intergenic RNAs, antisense RNAs, 5’ and 3’ untranslated regions and operons. The expression of these RNA molecules was confirmed by 5’-3’RACE protocol and specific RT-PCR. Here we describe the complete transcriptome of Nm during incubation in human blood; we were able to identify new proteins important for survival in human blood and also to identify additional roles of previously known virulence factors in aiding survival in blood. In addition the tiling array analysis demonstrated that Nm expresses a set of new transcripts, not previously identified, and suggests the presence of a circuit of regulatory RNA elements used by Nm to adapt to proliferate in human blood.

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.

Activity and mechanisms of action of novel organosulfur derivatives of the HDAC inhibitor Valproic Acid in human experimental models of non-small-cell lung cancer

Non-small-cell lung cancer (NSCLC) represents the leading cause of cancer death worldwide, and 5-year survival is about 16% for patients diagnosed with advanced lung cancer and about 70-90% when the disease is diagnosed and treated at earlier stages. Treatment of NSCLC is changed in the last years with the introduction of targeted agents, such as gefitinib and erlotinib, that have dramatically changed the natural history of NSCLC patients carrying specific mutations in the EGFR gene, or crizotinib, for patients with the EML4-ALK translocation. However, such patients represent only about 15-20% of all NSCLC patients, and for the remaining individuals conventional chemotherapy represents the standard choice yet, but response rate to thise type of treatment is only about 20%. Development of new drugs and new therapeutic approaches are so needed to improve patients outcome. In this project we aimed to analyse the antitumoral activity of two compounds with the ability to inhibit histone deacethylases (ACS 2 and ACS 33), derived from Valproic Acid and conjugated with H2S, in human cancer cell lines derived from NSCLC tissues. We showed that ACS 2 represents the more promising agent. It showed strong antitumoral and pro-apoptotic activities, by inducing membrane depolarization, cytocrome-c release and caspase 3 and 9 activation. It was able to reduce the invasive capacity of cells, through inhibition of metalloproteinases expression, and to induce a reduced chromatin condensation. This last characteristic is probably responsible for the observed high synergistic activity in combination with cisplatin. In conclusion our results highlight the potential role of the ACS 2 compound as new therapeutic option for NSCLC patients, especially in combination with cisplatin. If validated in in vivo models, this compound should be worthy for phase I clinical trials.

The elderly health status and its correlation with ageing biomarkers: the European Project Mark-Age

According to the latest statistics projections formulated by Eurostat, the proportion of elderly EU-27’s population aged over 65 years old is predicted to increase from 17.5 % in 2011 to 29.5 % by 2060. This "population explosion" makes extremely important to identify the different genetic and molecular mechanisms which underpin the morbidity and mortality along with new strategies able to counteract or slow down its progress. In this scenario fits the European Project MARK-AGE whose aim was to identify a robust set of biomarkers of human ageing able to discriminate between chronological and biological ageing and to derive a model for healthy ageing through the analysis of three populations from different European countries, supposed to be characterized by different ageing rate: 1. Subjects representing the “Normal” or “Physiological” aging. 2. Subjects representing the “successful” or “decelerate” aging 3. Subjects representing the “accelerated” aging. The aim of this work was to recruit and characterize volunteers, to perform an accurate analysis of the health status of elderly recruited subjects (60-79 years) verifying any possible dissimilarity in their aging trajectories, to identify a set of robust ageing biomarkers and investigate possible correlations between ageing biomarkers and health status of recruited volunteers. The model proposed by MARK-AGE Project regarding different ageing trajectories has been confirmed and several ageing biomarkers have been identified.