Y-chromosome and Surname Analyses for Reconstructing Past Population Structures: The Sardinian Population as a Test Case.

Many anthropological, linguistic, genetic and genomic analyses have been carried out to evaluate the potential impact that evolutionary forces had in shaping the present-day Sardinian gene pool, the main outlier in the genetic landscape of Europe. However, due to the homogenizing effect of internal movements, which have intensified over the past fifty years, only partial information has been obtained about the main demographic events. To overcome this limitation, we analyzed the male-specific region of the Y chromosome in three population samples obtained by reallocating a large number of Sardinian subjects to the place of origin of their monophyletic surnames, which are paternally transmitted through generations in most of the populations, much like the Y chromosome. Three Y-chromosome founding lineages, G2-L91, I2-M26 and R1b-V88, were identified as strongly contributing to the definition of the outlying position of Sardinians in the European genetic context and marking a significant differentiation within the island. The present distribution of these lineages does not always mirror that detected in ancient DNAs. Our results show that the analysis of the Y-chromosome gene pool coupled with a sampling method based on the origin of the family name, is an efficient approach to unravelling past heterogeneity, often hidden by recent movements, in the gene pool of modern populations. Furthermore, the reconstruction and comparison of past genetic isolates represent a starting point to better assess the genetic information deriving from the increasing number of available ancient DNA samples.

Ripples on the surface. Surnames and genes in Sicily and Southern Italy.

BACKGROUND: Southern Italy and Sicily played a key role in the peopling history of the Mediterranean. While genetic research showed the remarkable homogeneity of these regions, surname-based studies instead suggested low population mobility, hence potential structuring. AIM: In order to better understand these different patterns, this study (1) thoroughly analysed the surname structure of Sicily and Southern Italy and (2) tested its relationships with a wide set of molecular markers. SUBJECTS AND METHODS: Surname data were collected from 1213 municipalities and compared to uniparental and autosomal genetic markers typed in ∼300 individuals from 8-10 populations. Surname analyses were performed using different multivariate methods, while comparisons with genetic data relied on correlation tests. RESULTS: Surnames were clearly structured according to regional geographic patterns, which likely emerged because of recent isolation-by-distance-like population dynamics. In general, genetic markers, hinting at a pervasive homogeneity, did not correlate with surname distribution. However, long autosomal haplotypes (>5 cM) that compared to genotypic (SNPs) data identify more "recent" relatedness, showing a clear association with surname patterns. CONCLUSION: The apparent contradiction between surname structure and genetic homogeneity was resolved by figuring surnames as recent "ripples" deposited on a vast and ancient homogeneous genetic "surface".

Consanguinity and late fertility: spatial analysis reveals positive association patterns.

The role of consanguinity on human complex traits is an important and controversial issue. In this work we focused on the Sardinian population and examined the effect of consanguineous unions on late female fertility. During the last century the island has been characterized by a high incidence of marriages between relatives, favoured by socio economic conditions and geographical isolation, and by high fertility despite a widespread tendency to delay reproduction. Through spatial analysis techniques, we explored the geographical heterogeneity of consanguinity and late fertility, and identified in Central-Eastern Sardinia a common area with an excess of both traits, where the traits are positively associated. We found that their association did not significantly affect women's fertility in the area, despite the expected negative role of both traits. Intriguingly, this critical zone corresponds well to areas reported by previous studies as being peculiar for a high frequency of centenarians and for lower risk in pregnancy outcome. The proposed approach can be generally exploited to identify target populations on which socioeconomic, biodemographic and genetic data can be collected at the individual level, and deeper analyses carried out to disentangle the determinants of complex biological traits and to investigate their association.

Characterization of the biological processes shaping the genetic structure of the Italian population.

BACKGROUND: The genetic structure of human populations is the outcome of the combined action of different processes such as demographic dynamics and natural selection. Several efforts toward the characterization of population genetic architectures and the identification of adaptation signatures were recently made. In this study, we provide a genome-wide depiction of the Italian population structure and the analysis of the major determinants of the current existing genetic variation. RESULTS: We defined and characterized 210 genomic loci associated with the first Principal Component calculated on the Italian genotypic data and correlated to the North-south genetic gradient. Using a gene-enrichment approach we identified the immune function as primarily involved in the Italian population differentiation and we described a locus on chromosome 13 showing combined evidence of North-south diversification in allele frequencies and signs of recent positive selection. In this region our bioinformatics analysis pinpointed an uncharacterized long intergenic non-coding (lincRNA), whose expression appeared specific for immune-related tissues suggesting its relevance for the immune function. CONCLUSIONS: Our study, combining population genetic analyses with biological insights provides a description of the Italian genetic structure that in future could contribute to the evaluation of complex diseases risk in the population context.

CorrelaGenes: a new tool for the interpretation of the human transcriptome.

BACKGROUND: The amount of gene expression data available in public repositories has grown exponentially in the last years, now requiring new data mining tools to transform them in information easily accessible to biologists. RESULTS: By exploiting expression data publicly available in the Gene Expression Omnibus (GEO) database, we developed a new bioinformatics tool aimed at the identification of genes whose expression appeared simultaneously altered in different experimental conditions, thus suggesting co-regulation or coordinated action in the same biological process. To accomplish this task, we used the 978 human GEO Curated DataSets and we manually performed the selection of 2,109 pair-wise comparisons based on their biological rationale. The lists of differentially expressed genes, obtained from the selected comparisons, were stored in a PostgreSQL database and used as data source for the CorrelaGenes tool. Our application uses a customized Association Rule Mining (ARM) algorithm to identify sets of genes showing expression profiles correlated with a gene of interest. The significance of the correlation is measured coupling the Lift, a well-known standard ARM index, and the χ(2) p value. The manually curated selection of the comparisons and the developed algorithm constitute a new approach in the field of gene expression profiling studies. Simulation performed on 100 randomly selected target genes allowed us to evaluate the efficiency of the procedure and to obtain preliminary data demonstrating the consistency of the results. CONCLUSIONS: The preliminary results of the simulation showed how CorrelaGenes could contribute to the characterization of molecular pathways and biological processes integrating data obtained from other applications and available in public repositories.

Phosphorylation of SRSF1 is modulated by replicational stress.

DNA ligase I-deficient 46BR.1G1 cells show a delay in the maturation of replicative intermediates resulting in the accumulation of single- and double-stranded DNA breaks. As a consequence the ataxia telangiectasia mutated protein kinase (ATM) is constitutively phosphorylated at a basal level. Here, we use 46BR.1G1 cells as a model system to study the cell response to chronic replication-dependent DNA damage. Starting from a proteomic approach, we demonstrate that the phosphorylation level of factors controlling constitutive and alternative splicing is affected by the damage elicited by DNA ligase I deficiency. In particular, we show that SRSF1 is hyperphosphorylated in 46BR.1G1 cells compared to control fibroblasts. This hyperphosphorylation can be partially prevented by inhibiting ATM activity with caffeine. Notably, hyperphosphorylation of SRSF1 affects the subnuclear distribution of the protein and the alternative splicing pattern of target genes. We also unveil a modulation of SRSF1 phosphorylation after exposure of MRC-5V1 control fibroblasts to different exogenous sources of DNA damage. Altogether, our observations indicate that a relevant aspect of the cell response to DNA damage involves the post-translational regulation of splicing factor SRSF1 which is associated with a shift in the alternative splicing program of target genes to control cell survival or cell death.

Long-term outcome and lineage-specific chimerism in 194 patients with Wiskott-Aldrich syndrome treated by hematopoietic cell transplantation in the period 1980-2009: an international collaborative study.

In this retrospective collaborative study, we have analyzed long-term outcome and donor cell engraftment in 194 patients with Wiskott-Aldrich syndrome (WAS) who have been treated by hematopoietic cell transplantation (HCT) in the period 1980- 2009. Overall survival was 84.0% and was even higher (89.1% 5-year survival) for those who received HCT since the year 2000, reflecting recent improvement of outcome after transplantation from mismatched family donors and for patients who received HCT from an unrelated donor at older than 5 years. Patients who went to transplantation in better clinical conditions had a lower rate of post-HCT complications. Retrospective analysis of lineage-specific donor cell engraftment showed that stable full donor chimerism was attained by 72.3% of the patients who survived for at least 1 year after HCT. Mixed chimerism was associated with an increased risk of incomplete reconstitution of lymphocyte count and post-HCT autoimmunity, and myeloid donor cell chimerism < 50% was associated with persistent thrombocytopenia. These observations indicate continuous improvement of outcome after HCT for WAS and may have important implications for the development of novel protocols aiming to obtain full correction of the disease and reduce post-HCT complications.

Genetic variability at COMT but not at OPRM1 and UGT2B7 loci modulates morphine analgesic response in acute postoperative pain.

PURPOSE: To investigate interindividual variability in response to pain treatment, we characterized postoperative patients for morphine metabolism and for COMT, OPRM1 and UGT2B7 polymorphisms. METHODS: A total of 109 patients treated with morphine were genotyped by DNA sequencing for 12 DNA polymorphisms of the COMT, OPRM1 and UGT2B7 genes. The plasma concentration of morphine and of M3G/M6G metabolites were evaluated by means of reversed phase high-performance liquid chromatography coupled with mass spectrometry. RESULTS: An association between average morphine consumption during the first 24 postoperative hours by patient-controlled analgesia (PCA) and COMT haplotypes was found. Specifically, patients with the diplotype for average pain intensity (APS/APS) required the lowest morphine doses compared to the other subjects (p = 0.011). The APS haplotype contains an adenine corresponding to methionine, instead of valine, at position 158 of the COMT protein. Met/Met homozygous patients consumed significantly lower morphine doses than other subjects (p = 0.014); accordingly, Val158Met genotyping alone might be used in the clinical setting to predict PCA morphine need. Considering both COMT Val158Met and OPRM1 A118G polymorphisms, carriers of both the Met/Met and AA genotypes required less morphine than other subjects, although the difference was not significant. The analysis of UGT2B7 revealed the occurrence of two common haplotypes (G_C_C_A_C and A_T_T_G_T) that did not prove to be related with plasma morphine and M3G/M6G concentration. CONCLUSIONS: By considering COMT, OPRM1, and UGT2B7 genotypes, as well as pharmacokinetic results, only COMT polymorphisms appear to be predictive of morphine need in postoperative pain therapy.

Are myocardial infarction--associated single-nucleotide polymorphisms associated with ischemic stroke?

BACKGROUND AND PURPOSE: Ischemic stroke (IS) shares many common risk factors with coronary artery disease (CAD). We hypothesized that genetic variants associated with myocardial infarction (MI) or CAD may be similarly involved in the etiology of IS. To test this hypothesis, we evaluated whether single-nucleotide polymorphisms (SNPs) at 11 different loci recently associated with MI or CAD through genome-wide association studies were associated with IS. METHODS: Meta-analyses of the associations between the 11 MI-associated SNPs and IS were performed using 6865 cases and 11 395 control subjects recruited from 9 studies. SNPs were either genotyped directly or imputed; in a few cases a surrogate SNP in high linkage disequilibrium was chosen. Logistic regression was performed within each study to obtain study-specific ßs and standard errors. Meta-analysis was conducted using an inverse variance weighted approach assuming a random effect model. RESULTS: Despite having power to detect odds ratio of 1.09-1.14 for overall IS and 1.20-1.32 for major stroke subtypes, none of the SNPs were significantly associated with overall IS and/or stroke subtypes after adjusting for multiple comparisons. CONCLUSIONS: Our results suggest that the major common loci associated with MI risk do not have effects of similar magnitude on overall IS but do not preclude moderate associations restricted to specific IS subtypes. Disparate mechanisms may be critical in the development of acute ischemic coronary and cerebrovascular events.

General method to unravel ancient population structures through surnames, final validation on Italian data.

We analyze the geographic location of 77,451 different Italian surnames (17,579,891 individuals) obtained from the lists of telephone subscribers of the year 1993. By using a specific neural network analysis (Self-Organizing Maps, SOMs), we automatically identify the geographic origin of 49,117 different surnames. To validate the methodology, we compare the results to a study, previously conducted, on the same database, with accurate supervised methods. By comparing the results, we find an overlap of 97%, meaning that the SOMs methodology is highly reliable and well traces back the geographic origin of surnames at the time of their introduction (Late Middle Ages/Renaissance in Italy). SOMs results enables one to distinguish monophyletic surnames from polyphyletic ones, that is surnames having had a single geographic and historic origin from those that started to be in use, with an identical spelling, in different locations (respectively, 76.06% and 21.05% of the total). As we are interested in geographic origins, polyphyletic surnames are excluded from further analyses. By comparing the present location of each monophyletic surname to its inferred geographic origin in late Middle Ages/Renaissance, we measure the extent of the migrations having occurred in Italy since that time. We find that the percentage of individuals presently living in the very area where their surname started to be in use centuries ago is extremely variable (ranging from 22.77% to 77.86% according to the province), thus meaning that self-assessed regional identities seldom correspond to the "autochthony" they imply. For example the upper part of the Thyrennian coast (Northern Latium, Tuscany) has a strong identity but few "autochthonous" inhabitants (∼28%) having been a passageway from the North to the South of Italy.

Spatial analysis of the aptitude to late maternity on the island of Sardinia.

This study examines local heterogeneity in the aptitude of Sardinian mothers towards late reproduction, and explores its temporal persistence and association with both post-reproductive longevity and propensity to consanguineous marriage. Data on women's fertility from 1961 and birth records for 1980-1996 from Vital Statistics were analysed by means of the following indicators: the incidence of old mothers at last childbirth, female mortality (1980-2001) at 80 years of age and over and the proportion of consanguineous marriages (1930-1969). A variable kernel-smoothing method was used to create interpretable representations of the true spatial structure of the indicators, and to highlight areas of higher than expected intensity. In particular, an area of reproductive and post-reproductive longevity was identified where the traits combine with a higher tendency to relatedness. Intriguingly, this area corresponds approximately to the geographically and historically well defined central-eastern zone, which was the refuge of Sardinians during past invasions, and overlaps the Ogliastra region, which has been widely studied for its genetic homogeneity.

GC-AG Introns Features in Long Non-coding and Protein-Coding Genes Suggest Their Role in Gene Expression Regulation.

Long non-coding RNAs (lncRNAs) are recognized as an important class of regulatory molecules involved in a variety of biological functions. However, the regulatory mechanisms of long non-coding genes expression are still poorly understood. The characterization of the genomic features of lncRNAs is crucial to get insight into their function. In this study, we exploited recent annotations by GENCODE to characterize the genomic and splicing features of long non-coding genes in comparison with protein-coding ones, both in human and mouse. Our analysis highlighted differences between the two classes of genes in terms of their gene architecture. Significant differences in the splice sites usage were observed between long non-coding and protein-coding genes (PCG). While the frequency of non-canonical GC-AG splice junctions represents about 0.8% of total splice sites in PCGs, we identified a significant enrichment of the GC-AG splice sites in long non-coding genes, both in human (3.0%) and mouse (1.9%). In addition, we found a positional bias of GC-AG splice sites being enriched in the first intron in both classes of genes. Moreover, a significant shorter length and weaker donor and acceptor sites were found comparing GC-AG introns to GT-AG introns. Genes containing at least one GC-AG intron were found conserved in many species, more prone to alternative splicing and a functional analysis pointed toward their enrichment in specific biological processes such as DNA repair. Our study shows for the first time that GC-AG introns are mainly associated with lncRNAs and are preferentially located in the first intron. Additionally, we discovered their regulatory potential indicating the existence of a new mechanism of non-coding and PCGs expression regulation.

Autoantibodies to poly(ADP-ribose) polymerase in centenarians: a reappraisal of Grabar's hypothesis.

BACKGROUND: Centenarians display a characteristic autoantibody profile, this being the absence of organ-specific autoantibodies and an increase in non-organ-specific autoantibodies without any full-blown autoimmune disease. OBJECTIVE: Antibodies directed to the nuclear protein poly(ADP-ribose) polymerase (PARP) were frequently found in the sera of patients affected by autoimmune diseases. This study aims at investigating the presence of circulating autoantibodies directed against PARP-1 in normal subjects, and searching for a possible correlation between level of circulating autoantibodies and age. METHODS: The presence of antibodies to PARP was monitored by ELISA according to a previously developed protocol. Data were analysed by parametric statistics (unpaired t test, chi(2) test). RESULTS: Our study performed on 33 centenarians and 66 subjects of age ranging from 12 to 80 years shows that circulating autoantibodies to the nuclear enzyme PARP, previously described in autoimmune diseases, are present in the sera of normal healthy people and increase with age. CONCLUSIONS: Taking into account the role of PARP-1 in DNA damage and apoptosis, the data are compatible with Grabar's hypothesis, which proposed several decades ago that 'autoreactive antibodies represent a physiological system for disposing the products of metabolism and catabolism', thereby helping to attain longevity.

Identification of a novel antiviral micro-RNA targeting the NS1 protein of the H1N1 pandemic human influenza virus and a corresponding viral escape mutation.

The influenza A virus (IAV) NS1 protein is one of the major regulators of pathogenicity, being able to suppress innate immune response and host protein synthesis. In this study we identified the human micro RNA hsa-miR-1307-3p as a novel potent suppressor of NS1 expression and influenza virus replication. Transcriptomic analysis indicates that hsa-miR-1307-3p also negatively regulates apoptosis and promotes cell proliferation. In addition, we identified a novel mutation in the NS1 gene of A(H1N1)pdm09 strains circulating in Italy in the 2010-11 season, which enabled the virus to escape the hsa-miR-1307-3p inhibition, conferring replicative advantage to the virus in human cells. To the best of our knowledge, this is the first validation of suppression of IAV H1N1 NS1 by a human micro RNA and the first example of an escape mutation from micro RNA-mediated antiviral response for the A(H1N1)pdm09 virus.

Genetic evidence for a worldwide chaotic dispersion pattern of the arbovirus vector, Aedes albopictus.

BACKGROUND: Invasive species represent a global concern for their rapid spread and the possibility of infectious disease transmission. This is the case of the global invader Aedes albopictus, the Asian tiger mosquito. This species is a vector of medically important arboviruses, notably chikungunya (CHIKV), dengue (DENV) and Zika (ZIKV). The reconstruction of the complex colonization pattern of this mosquito has great potential for mitigating its spread and, consequently, disease risks. METHODOLOGY/PRINCIPAL FINDINGS: Classical population genetics analyses and Approximate Bayesian Computation (ABC) approaches were combined to disentangle the demographic history of Aedes albopictus populations from representative countries in the Southeast Asian native range and in the recent and more recently colonized areas. In Southeast Asia, the low differentiation and the high co-ancestry values identified among China, Thailand and Japan indicate that, in the native range, these populations maintain high genetic connectivity, revealing their ancestral common origin. China appears to be the oldest population. Outside Southeast Asia, the invasion process in La Réunion, America and the Mediterranean Basin is primarily supported by a chaotic propagule distribution, which cooperates in maintaining a relatively high genetic diversity within the adventive populations. CONCLUSIONS/SIGNIFICANCE: From our data, it appears that independent and also trans-continental introductions of Ae. albopictus may have facilitated the rapid establishment of adventive populations through admixture of unrelated genomes. As a consequence, a great amount of intra-population variability has been detected, and it is likely that this variability may extend to the genetic mechanisms controlling vector competence. Thus, in the context of the invasion process of this mosquito, it is possible that both population ancestry and admixture contribute to create the conditions for the efficient transmission of arboviruses and for outbreak establishment.

An Association Rule Mining Approach to Discover lncRNAs Expression Patterns in Cancer Datasets.

In the past few years, the role of long noncoding RNAs (lncRNAs) in tumor development and progression has been disclosed although their mechanisms of action remain to be elucidated. An important contribution to the comprehension of lncRNAs biology in cancer could be obtained through the integrated analysis of multiple expression datasets. However, the growing availability of public datasets requires new data mining techniques to integrate and describe relationship among data. In this perspective, we explored the powerness of the Association Rule Mining (ARM) approach in gene expression data analysis. By the ARM method, we performed a meta-analysis of cancer-related microarray data which allowed us to identify and characterize a set of ten lncRNAs simultaneously altered in different brain tumor datasets. The expression profiles of the ten lncRNAs appeared to be sufficient to distinguish between cancer and normal tissues. A further characterization of this lncRNAs signature through a comodulation expression analysis suggested that biological processes specific of the nervous system could be compromised.

From surnames to the history of Y chromosomes: the Sardinian population as a paradigm.

A total of 202 Sardinian male subjects were examined for 13 biallelic stable markers, the complex 49a,f/TaqI system and three microsatellites of the Y chromosome in order to investigate, through surname analysis, on a possible territorial heterogeneity inside the island. The study of geographical distribution and linguistic derivation of Sardinian surnames allow us to discover their 'probable place of origin' and reconstruct ancient genetic isolates which borders are, today, no more recognizable. The molecular analysis revealed that about 90% of the Sardinian Y chromosomes fell into haplogroups E-M35, G-M201, I-M26, J-12f2 and R-M269. In contrast with the territorial homogeneity of these haplogroups, when the individuals were distributed according to their birthplace, a significant difference between the three historically and culturally distinct geographical areas into which Sardinia can be subdivided was observed when the individuals were distributed according to the ancestral location of surnames. In particular, the major contribution to this heterogeneity is due to the 'Sardinian-specific' haplogroup I-M26 (almost completely associated with the 49a,f-Ht12/12f2-10Kb/YCAIIa-21/YCAIIb-11 compound haplotype), which shows both a significantly higher incidence in the central-eastern (archaic) area and a significantly lower frequency in the northern area. The results of this study agree with the hypothesis that the ancestral homeland of this specific subset of haplogroup I is the mountainous central-eastern area of Sardinia, where the population underwent a long history of isolation since ancient times, and highlight the informative power of the surname analysis.

Chromosomal end-to-end fusions in immortalized mouse embryonic fibroblasts deficient in the DNA-dependent protein kinase catalytic subunit.

Telomeres constitute the ends of the linear eukaryotic chromosomes and are essential for the maintenance of chromosome stability and genome integrity. One of the consequences of an altered telomere structure is the formation of telomeric fusions (TFs), that is aberrant chromosomes in which two elements are fused at their telomeres. Proteins involved in the non-homologous end joining pathway for the repair of the DNA double strand breaks, as the DNA-dependent protein kinase catalytic subunit (DNA-PKcs), contribute to the formation of a functional telomere. To investigate the role of DNA-PKcs in telomere functionality, we studied the frequency of TFs in mouse embryonic fibroblasts obtained from animals in which the DNA-PKcs gene had been inactivated; the analysis was performed prior and after spontaneous immortalization in culture. Our results suggest that DNA-PKcs deficiency has a limited effect, if any, on TF formation in primary cells, while it further increases chromosomal instability in immortalized cells. In fact, the frequency of TFs was significantly higher in immortalized DNA-PKcs mutant cells compared to wild type cells. Together with TFs, we also found metacentric or submetacentric chromosomes in which no telomeric sequences were detected at the joining site. The frequency of this anomaly, that resembles the Robertsonian translocations observed in wild mice, was independent of the DNA-PKcs genotype. This suggests that the formation of these rearranged chromosomes does not rely on a functional DNA-PKcs.

Influence of COMT Val158Met polymorphism on Alzheimer's disease and mild cognitive impairment in Italian patients.

COMT (Catechol-O methyltransferase) gene is one of the key players in synaptic plasticity and in learning and memory mechanisms. A single nucleotide polymorphism (rs4680; G to A) in the COMT coding region causes Val158Met aminoacid substitution in the corresponding protein, with Val allele exhibiting a 3- to 4-fold increase in enzyme activity compared to Met. With the purpose of examining the influence of COMT as a genetic risk factor for cognitive impairment, we analyzed a sample of 248 healthy subjects, 276 patients affected by Alzheimer's disease (AD), and 70 subjects with mild cognitive impairment (MCI), the latter condition possibly representing a prodrome for dementia. All subjects were analyzed for COMT rs4680 polymorphism and APOE genotype. Our study strengthens data showing that APOE ε4 allele is an independent risk factor for AD and also a risk factor for MCI. Neither COMT alleles nor genotypes proved to be independently associated with the risk of AD or MCI in our sample. However, we found an association between COMT GG genotype (Val/Val) and APOE ε4 carrier status and the risk of AD and MCI. In particular, when GG genotype is included into the multinomial analysis, the risk of AD and MCI due to APOE ε4 allele is increased of about 2-3 fold; moreover, the risk conferred by the combination of G and ε4 alleles is more pronounced in male patients. To our knowledge, this synergistic effect is here shown for the first time on a population sample representative of Caucasian patients.

Search for cellular stress biomarkers in lymphocytes from patients with multiple sclerosis: a pilot study.

Multiple Sclerosis (MS) is a chronic disease of the central nervous system, the etiology of which, although not completely known, involves inflammation and autoimmunity. In the present study we aimed at identifying molecular markers of apoptosis, cellular stress and DNA damage in isolated peripheral blood mononuclear cells (PBMCs) of MS patients. The analysis was carried on 19 relapsing-remitting untreated MS patients and 13 healthy individuals. We investigated the emergency-driven synthesis of poly(ADP-ribose) (PAR), the expression level of the constitutive enzyme poly(ADP-ribose) polymerase-1 (PARP-1) and the DNA damage-induced phosphorylation of histone H2AX. PAR accumulation, PARP-1 and phosphorylated H2AX (γH2AX) were detected by immunofluorescence experiments on PBMCs isolated from 19 patients and 13 healthy volunteers. Our results show for the first time a net increased amount in PAR and γH2AX in MS patients compared to healthy individuals. Patients were further subdivided in three groups, according to the neuroimaging (MRI)-based classification of disease phase. Remarkably, we found a positive correlation between the level of γH2AX and MS aggressiveness. In addition, apoptosis in PBMCs was monitored by flow cytometry of both phosphatidylserine exposure (revealed by Annexin V-FITC labeling) and membrane permeability to propidium iodide. Our observations provide the evidence that the number of apoptotic cells was significantly higher in patients compared to healthy individuals, thus suggesting that apoptosis could affect MS lymphocyte function.