Nucleotide distance influences co-methylation between nearby CpG sites.

The tendency of individual CpG sites to be methylated is distinctive, non-random and well-regulated throughout the genome. We investigated the structural and spatial factors influencing CpGs methylation by performing an ultra-deep targeted methylation analysis on human, mouse and zebrafish genes. We found that methylation is not a random process and that closer neighboring CpG sites are more likely to share the same methylation status. Moreover, if the distance between CpGs increases, the degree of co-methylation decreases. We set up a simulation model to analyze the contribution of both the intrinsic susceptibility and the distance effect on the probability of a CpG to be methylated. Our finding suggests that the establishment of a specific methylation pattern follows a universal rule that must take into account of the synergistic and dynamic interplay of these two main factors: the intrinsic methylation susceptibility of specific CpG and the nucleotide distance between two CpG sites.

Whole transcriptome targeted gene quantification provides new insights on pulmonary sarcomatoid carcinomas.

Pulmonary sarcomatoid carcinomas (PSC) are a rare group of lung cancer with a median overall survival of 9-12 months. PSC are divided into five histotypes, challenging to diagnose and treat. The identification of PSC biomarkers is warranted, but PSC molecular profile remains to be defined. Herein, a targeted whole transcriptome analysis was performed on 14 PSC samples, evaluated also for the presence of the main oncogene mutations and rearrangements. PSC expression data were compared with transcriptome data of lung adenocarcinomas (LUAD) and squamous cell carcinomas (LUSC) from The Cancer Genome Atlas. Deregulated genes were used for pathway enrichment analysis; the most representative genes were tested by immunohistochemistry (IHC) in an independent cohort (30 PSC, 31 LUAD, 31 LUSC). All PSC cases were investigated for PD-L1 expression. Thirty-eight genes deregulated in PSC were identified, among these IGJ and SLMAP were confirmed by IHC. Moreover, Forkhead box signaling and Fanconi anemia pathways were specifically enriched in PSC. Finally, some PSC harboured alterations in genes targetable by tyrosine kinase inhibitors, as EGFR and MET. We provide a deep molecular characterization of PSC; the identification of specific molecular profiles, besides increasing our knowledge on PSC biology, might suggest new strategies to improve patients management.

DNA Methylation variability among individuals is related to CpGs cluster density and evolutionary signatures.

BACKGROUND: In recent years, epigenetics has gained a central role in the understanding of the process of natural selection. It is now clear how environmental impacts on the methylome could promote methylation variability with direct effects on disease etiology as well as phenotypic and genotypic variations in evolutionary processes. To identify possible factors influencing inter-individual methylation variability, we studied methylation values standard deviation of 166 healthy individuals searching for possible associations with genomic features and evolutionary signatures. RESULTS: We analyzed methylation variability values in relation to CpG cluster density and we found a strong association between them (p-value < 2.2 × 10- 16). Furthermore, we found that genes related to CpGs with high methylation variability values were enriched for immunological pathways; instead, those associated with low ones were enriched for pathways related to basic cellular functions. Finally, we found an association between methylation variability values and signals of both ancient (p-value < 2.2 × 10- 16) and recent selective pressure (p-value < 1 × 10- 4). CONCLUSION: Our results indicate the presence of an intricate interplay between genetics, epigenetic code and evolutionary constraints in humans.

Tracking the evolution of epialleles during neural differentiation and brain development: D-Aspartate oxidase as a model gene.

We performed ultra-deep methylation analysis at single molecule level of the promoter region of developmentally regulated D-Aspartate oxidase (Ddo), as a model gene, during brain development and embryonic stem cell neural differentiation. Single molecule methylation analysis enabled us to establish the effective epiallele composition within mixed or pure brain cell populations. In this framework, an epiallele is defined as a specific combination of methylated CpG within Ddo locus and can represent the epigenetic haplotype revealing a cell-to-cell methylation heterogeneity. Using this approach, we found a high degree of polymorphism of methylated alleles (epipolymorphism) evolving in a remarkably conserved fashion during brain development. The different sets of epialleles mark stage, brain areas, and cell type and unravel the possible role of specific CpGs in favoring or inhibiting local methylation. Undifferentiated embryonic stem cells showed non-organized distribution of epialleles that apparently originated by stochastic methylation events on individual CpGs. Upon neural differentiation, despite detecting no changes in average methylation, we observed that the epiallele distribution was profoundly different, gradually shifting toward organized patterns specific to the glial or neuronal cell types. Our findings provide a deep view of gene methylation heterogeneity in brain cell populations promising to furnish innovative ways to unravel mechanisms underlying methylation patterns generation and alteration in brain diseases.

ampliMethProfiler: a pipeline for the analysis of CpG methylation profiles of targeted deep bisulfite sequenced amplicons.

BACKGROUND: CpG sites in an individual molecule may exist in a binary state (methylated or unmethylated) and each individual DNA molecule, containing a certain number of CpGs, is a combination of these states defining an epihaplotype. Classic quantification based approaches to study DNA methylation are intrinsically unable to fully represent the complexity of the underlying methylation substrate. Epihaplotype based approaches, on the other hand, allow methylation profiles of cell populations to be studied at the single molecule level. For such investigations, next-generation sequencing techniques can be used, both for quantitative and for epihaplotype analysis. Currently available tools for methylation analysis lack output formats that explicitly report CpG methylation profiles at the single molecule level and that have suited statistical tools for their interpretation. RESULTS: Here we present ampliMethProfiler, a python-based pipeline for the extraction and statistical epihaplotype analysis of amplicons from targeted deep bisulfite sequencing of multiple DNA regions. CONCLUSIONS: ampliMethProfiler tool provides an easy and user friendly way to extract and analyze the epihaplotype composition of reads from targeted bisulfite sequencing experiments. ampliMethProfiler is written in python language and requires a local installation of BLAST and (optionally) QIIME tools. It can be run on Linux and OS X platforms. The software is open source and freely available at http://amplimethprofiler.sourceforge.net .

Modeling DNA methylation by analyzing the individual configurations of single molecules.

DNA methylation is often analyzed by reporting the average methylation degree of each cytosine. In this study, we used a single molecule methylation analysis in order to look at the methylation conformation of individual molecules. Using D-aspartate oxidase as a model gene, we performed an in-depth methylation analysis through the developmental stages of 3 different mouse tissues (brain, lung, and gut), where this gene undergoes opposite methylation destiny. This approach allowed us to track both methylation and demethylation processes at high resolution. The complexity of these dynamics was markedly simplified by introducing the concept of methylation classes (MCs), defined as the number of methylated cytosines per molecule, irrespective of their position. The MC concept smooths the stochasticity of the system, allowing a more deterministic description. In this framework, we also propose a mathematical model based on the Markov chain. This model aims to identify the transition probability of a molecule from one MC to another during methylation and demethylation processes. The results of our model suggest that: 1) both processes are ruled by a dominant class of phenomena, namely, the gain or loss of one methyl group at a time; and 2) the probability of a single CpG site becoming methylated or demethylated depends on the methylation status of the whole molecule at that time.

Candidate genes and pathways downstream of PAX8 involved in ovarian high-grade serous carcinoma.

Understanding the biology and molecular pathogenesis of ovarian epithelial cancer (EOC) is key to developing improved diagnostic and prognostic indicators and effective therapies. Although research has traditionally focused on the hypothesis that high-grade serous carcinoma (HGSC) arises from the ovarian surface epithelium (OSE), recent studies suggest that additional sites of origin exist and a substantial proportion of cases may arise from precursor lesions located in the Fallopian tubal epithelium (FTE). In FTE cells, the transcription factor PAX8 is a marker of the secretory cell lineage and its expression is retained in 96% of EOC. We have recently reported that PAX8 is involved in the tumorigenic phenotype of ovarian cancer cells. In this study, to uncover genes and pathways downstream of PAX8 involved in ovarian carcinoma we have determined the molecular profiles of ovarian cancer cells and in parallel of Fallopian tube epithelial cells by means of a silencing approach followed by an RNA-seq analysis. Interestingly, we highlighted the involvement of pathways like WNT signaling, epithelial-mesenchymal transition, p53 and apoptosis. We believe that our analysis has led to the identification of candidate genes and pathways regulated by PAX8 that could be additional targets for the therapy of ovarian carcinoma.

The PPARγ2 Pro12Ala variant is protective against progression of nephropathy in people with type 2 diabetes.

OBJECTIVE: Cross-sectional studies suggest the association between diabetic nephropathy and the PPARγ2 Pro12Ala polymorphism of the peroxisome proliferator-activated receptor γ2 (PPARγ2). Prospective data are limited to microalbuminuria and no information on renal function is available to date. The present study evaluates the association between the Pro12Ala polymorphism of PPARγ2 and the progression of albuminuria and decay in glomerular filtration rate (GFR) in type 2 diabetes. PATIENTS AND MEASUREMENTS: We studied 256 patients with an average 5-year follow-up. Among others, urinary albumin excretion rate (UAER) was measured on spot sample, GFR was estimated with the CKD-EPI Equation. RESULTS: Baseline UAER and GFR were similar for carriers or non-carriers of the polymorphism. At follow-up no significant changes from baseline were observed for UAER or eGFR in carriers of the Pro12Ala polymorphism whereas a significant increase in UAER [17 (11.3-37.9) versus 24.5 (13.8-49.9) µg/mg, p < 0.006)] and a significant reduction in the eGFR (82.8 ± 14.5 versus 80.3 ± 17.3 ml/min/1.73, m(2) p = 0.02), were observed in non carriers of the Pro12Ala polymorphism. Progression of nephropathy - defined according to a combined end point of UAER and eGFR- i.e. doubling of baseline UAER to at least 100 µg/mg, or new onset microalbuminuria, or progression from micro to macroalbuminuria, or 25% reduction of eGFR, or annualized eGFR decline >3 ml/min/year - was significantly less frequent in Ala carriers than non carriers (11.4% vs 35.8%; p < 0.01); HR adjusted for baseline age, AER, eGFR, HbA1c, diabetes duration and blood pressure was 0.32 (0.12-0.80). CONCLUSIONS: This study found that among patients with type 2 diabetes, the PPARγ2 Pro12Ala polymorphism is protective against progression of nephropathy and decay of renal function independent of major confounders.

Evidence for evolutionary and nonevolutionary forces shaping the distribution of human genetic variants near transcription start sites.

The regions surrounding transcription start sites (TSSs) of genes play a critical role in the regulation of gene expression. At the same time, current evidence indicates that these regions are particularly stressed by transcription-related mutagenic phenomena. In this work we performed a genome-wide analysis of the distribution of single nucleotide polymorphisms (SNPs) inside the 10 kb region flanking human TSSs by dividing SNPs into four classes according to their frequency (rare, two intermediate classes, and common). We found that, in this 10 kb region, the distribution of variants depends on their frequency and on their localization relative to the TSS. We found that the distribution of variants is generally different for TSSs located inside or outside of CpG islands. We found a significant relationship between the distribution of rare variants and nucleosome occupancy scores. Furthermore, our analysis suggests that evolutionary (purifying selection) and nonevolutionary (biased gene conversion) forces both play a role in determining the relative SNP frequency around TSSs. Finally, we analyzed the potential pathogenicity of each class of variant using the Combined Annotation Dependent Depletion score. In conclusion, this study provides a novel and detailed view of the distribution of genomic variants around TSSs, providing insight into the forces that instigate and maintain variability in such critical regions.

A distinct group of CpG islands shows differential DNA methylation between replicas of the same cell line in vitro.

BACKGROUND: CpG dinucleotide-rich genomic DNA regions, known as CpG islands (CGIs), can be methylated at their cytosine residues as an epigenetic mark that is stably inherited during cell mitosis. Differentially methylated regions (DMRs) are genomic regions showing different degrees of DNA methylation in multiple samples. In this study, we focused our attention on CGIs showing different DNA methylation between two culture replicas of the same cell line. RESULTS: We used methylation data of 35 cell lines from the Encyclopedia of DNA Elements (ENCODE) consortium to identify CpG islands that were differentially methylated between replicas of the same cell line and denoted them Inter Replicas Differentially Methylated CpG islands (IRDM-CGIs). We identified a group of IRDM-CGIs that was consistently shared by different cell lines, and denoted it common IRDM-CGIs. X chromosome CGIs were overrepresented among common IRDM-CGIs. Autosomal IRDM-CGIs were preferentially located in gene bodies and intergenic regions had a lower G + C content, a smaller mean length, and a reduced CpG percentage. Functional analysis of the genes associated with autosomal IRDM-CGIs showed that many of them are involved in DNA binding and development. CONCLUSIONS: Our results show that several specific functional and structural features characterize common IRDM-CGIs. They may represent a specific subset of CGIs that are more prone to being differentially methylated for their intrinsic characteristics.

CpG islands under selective pressure are enriched with H3K4me3, H3K27ac and H3K36me3 histone modifications.

BACKGROUND: Histone modification is an epigenetic mechanism that influences gene regulation in eukaryotes. In particular, histone modifications in CpG islands (CGIs) are associated with different chromatin states and with transcription activity. Changes in gene expression play a crucial role in adaptation and evolution. RESULTS: In this paper, we have studied, using a computational biology approach, the relationship between histone modifications in CGIs and selective pressure in Homo sapiens. We considered three histone modifications: histone H3 lysine 4 trimethylation (H3K4me3), histone H3 lysine 27 acetylation (H3K27ac) and histone H3 lysine 36 trimethylation (H3K36me3), and we used the publicly available genomic-scale histone modification data of thirteen human cell lines. To define regions under selective pressure, we used three distinct signatures that mark selective events from different evolutionary periods. We found that CGIs under selective pressure showed significant enrichments for histone modifications. CONCLUSION: Our result suggests that, CGIs that have undergone selective events are characterized by epigenetic signatures, in particular, histone modifications that are distinct from CGIs with no evidence of selection.

Can telomere shortening in human peripheral blood leukocytes serve as a disease biomarker of Friedreich's ataxia?

Enhanced oxidative stress and inflammation contribute to telomere erosion. Friedreich's ataxia is a neurodegenerative disorder caused by a reduction in frataxin expression that results in mitochondrial dysfunction and oxidative damage. Furthermore, frataxin deficiency induces a strong activation of inflammatory genes and neuronal death. We investigated telomere length (TL) in peripheral blood leukocytes of 37 patients with Friedreich's ataxia and 36 controls. We noted a significant telomere shortening in patients with Friedreich's ataxia compared to healthy controls (p=0.03). We also found a correlation between TL and disease duration (p=0.001). Our observations lead to the hypothesis that the TL of human peripheral blood leukocytes may serve as a biomarker of Friedreich's ataxia that could be used as an outcome measure in clinical trials.

Signs of selective pressure on genetic variants affecting human height.

Many decades of scientific investigation have proved the role of selective pressure in Homo Sapiens at least at the level of individual genes or loci. Nevertheless, there are examples of polygenic traits that are bound to be under selection, but studies devoted to apply population genetics methods to unveil such occurrence are still lacking. Stature provides a relevant example of well-studied polygenic trait for which is now available a genome-wide association study which has identified the genes involved in this trait, and which is known to be under selection. We studied the behavior of F(ST) in a simulated toy model to detect population differentiation on a generic polygenic phenotype under selection. The simulations showed that the set of alleles involved in the trait has a higher mean F(ST) value than those undergoing genetic drift only. In view of this we looked for an increase in the mean F(ST) value of the 180 variants associated to human height. For this set of alleles we found F(ST) to be significantly higher than the genomic background (p = 0.0356). On the basis of a statistical analysis we excluded that the increase was just due to the presence of outliers. We also proved as marginal the role played by local adaptation phenomena, even on different phenotypes in linkage disequilibrium with genetic variants involved in height. The increase of F(ST) for the set of alleles involved in a polygenic trait seems to provide an example of symmetry breaking phenomenon concerning the population differentiation. The splitting in the allele frequencies would be driven by the initial conditions in the population dynamics which are stochastically modified by events like drift, bottlenecks, etc, and other stochastic events like the born of new mutations.

Shorter telomeres in patients with cerebral autosomal dominant arteriopathy and leukoencephalopathy (CADASIL).

CADASIL is a hereditary systemic vasculopathy which affects mainly small cerebral arteries and is caused by mutations in the Notch3 gene. Misfolding of Notch3 is linked to endoplasmic reticulum stress and increased reactive oxygen species, which may result in dysfunction of endothelial cells, inflammation and ischemia. Oxidative stress and inflammation may induce a rapid telomere shortening in peripheral blood leukocytes (PBLs). The aim of this study was to assess the telomere length in PBLs from 29 patients with a genetic diagnosis of CADASIL by using a modified quantitative real-time polymerase chain reaction based assay. PBL telomere length was significantly shorter in CADASIL patients (T/S ratio = 0.17, 95% CI, 0.14-0.20) than in the controls (T/S ratio = 0.31, 95% CI, 0.27-0.35, t-test p < 0.001). Moreover, patients with functional dependence displayed shorter telomeres than those with functional independence (p = 0.039). Our data provide the first evidence that PBL telomere length is shortened in CADASIL disease, and this may be a systemic oxidative stress indicator in CADASIL patients, providing a possible biomarker of disease progression and for future therapeutic strategies.

CpG islands undermethylation in human genomic regions under selective pressure.

DNA methylation at CpG islands (CGIs) is one of the most intensively studied epigenetic mechanisms. It is fundamental for cellular differentiation and control of transcriptional potential. DNA methylation is involved also in several processes that are central to evolutionary biology, including phenotypic plasticity and evolvability. In this study, we explored the relationship between CpG islands methylation and signatures of selective pressure in Homo Sapiens, using a computational biology approach. By analyzing methylation data of 25 cell lines from the Encyclopedia of DNA Elements (ENCODE) Consortium, we compared the DNA methylation of CpG islands in genomic regions under selective pressure with the methylation of CpG islands in the remaining part of the genome. To define genomic regions under selective pressure, we used three different methods, each oriented to provide distinct information about selective events. Independently of the method and of the cell type used, we found evidences of undermethylation of CGIs in human genomic regions under selective pressure. Additionally, by analyzing SNP frequency in CpG islands, we demonstrated that CpG islands in regions under selective pressure show lower genetic variation. Our findings suggest that the CpG islands in regions under selective pressure seem to be somehow more "protected" from methylation when compared with other regions of the genome.

Schizophrenia and vitamin D related genes could have been subject to latitude-driven adaptation.

BACKGROUND: Many natural phenomena are directly or indirectly related to latitude. Living at different latitudes, indeed, has its consequences with being exposed to different climates, diets, light/dark cycles, etc. In humans, one of the best known examples of genetic traits following a latitudinal gradient is skin pigmentation. Nevertheless, also several diseases show latitudinal clinals such as hypertension, cancer, dismetabolic conditions, schizophrenia, Parkinson's disease and many more. RESULTS: We investigated, for the first time on a wide genomic scale, the latitude-driven adaptation phenomena. In particular, we selected a set of genes showing signs of latitude-dependent population differentiation. The biological characterization of these genes showed enrichment for neural-related processes. In light of this, we investigated whether genes associated to neuropsychiatric diseases were enriched by Latitude-Related Genes (LRGs). We found a strong enrichment of LRGs in the set of genes associated to schizophrenia. In an attempt to try to explain this possible link between latitude and schizophrenia, we investigated their associations with vitamin D. We found in a set of vitamin D related genes a significant enrichment of both LRGs and of genes involved in schizophrenia. CONCLUSIONS: Our results suggest a latitude-driven adaptation for both schizophrenia and vitamin D related genes. In addition we confirm, at a molecular level, the link between schizophrenia and vitamin D. Finally, we discuss a model in which schizophrenia is, at least partly, a maladaptive by-product of latitude dependent adaptive changes in vitamin D metabolism.

Uncoupling protein 2 G(-866)A polymorphism: a new gene polymorphism associated with C-reactive protein in type 2 diabetic patients.

BACKGROUND: This study evaluated the relationship between the G(-866)A polymorphism of the uncoupling protein 2 (UCP2) gene and high-sensitivity C reactive protein (hs-CRP) plasma levels in diabetic patients. METHODS: We studied 383 unrelated people with type 2 diabetes aged 40-70 years. Anthropometry, fasting lipids, glucose, HbA1c, and hs-CRP were measured. Participants were genotyped for the G (-866)A polymorphism of the uncoupling protein 2 gene. RESULTS: Hs-CRP (mg/L) increased progressively across the three genotype groups AA, AG, or GG, being respectively 3.0 ± 3.2, 3.6 ± 5.0, and 4.8 ± 5.3 (p for trend = 0.03). Since hs-CRP values were not significantly different between AA and AG genotype, these two groups were pooled for further analyses. Compared to participants with the AA/AG genotypes, homozygotes for the G allele (GG genotype) had significantly higher hs-CRP levels (4.8 ± 5.3 vs 3.5 ± 4.7 mg/L, p = 0.01) and a larger proportion (53.9% vs 46.1%, p = 0.013) of elevated hs-CRP (> 2 mg/L). This was not explained by major confounders such as age, gender, BMI, waist circumference, HbA1c, smoking, or medications use which were comparable in the two genotype groups. CONCLUSIONS: The study shows for the first time, in type 2 diabetic patients, a significant association of hs-CRP levels with the G(-866)A polymorphism of UCP2 beyond the effect of major confounders.

PGC-1alpha down-regulation affects the antioxidant response in Friedreich's ataxia.

BACKGROUND: Cells from individuals with Friedreich's ataxia (FRDA) show reduced activities of antioxidant enzymes and cannot up-regulate their expression when exposed to oxidative stress. This blunted antioxidant response may play a central role in the pathogenesis. We previously reported that Peroxisome Proliferator Activated Receptor Gamma (PPARgamma) Coactivator 1-alpha (PGC-1alpha), a transcriptional master regulator of mitochondrial biogenesis and antioxidant responses, is down-regulated in most cell types from FRDA patients and animal models. METHODOLOGY/PRINCIPAL FINDINGS: We used primary fibroblasts from FRDA patients and the knock in-knock out animal model for the disease (KIKO mouse) to determine basal superoxide dismutase 2 (SOD2) levels and the response to oxidative stress induced by the addition of hydrogen peroxide. We measured the same parameters after pharmacological stimulation of PGC-1alpha. Compared to control cells, PGC-1alpha and SOD2 levels were decreased in FRDA cells and did not change after addition of hydrogen peroxide. PGC-1alpha direct silencing with siRNA in control fibroblasts led to a similar loss of SOD2 response to oxidative stress as observed in FRDA fibroblasts. PGC-1alpha activation with the PPARgamma agonist (Pioglitazone) or with a cAMP-dependent protein kinase (AMPK) agonist (AICAR) restored normal SOD2 induction. Treatment of the KIKO mice with Pioglitazone significantly up-regulates SOD2 in cerebellum and spinal cord. CONCLUSIONS/SIGNIFICANCE: PGC-1alpha down-regulation is likely to contribute to the blunted antioxidant response observed in cells from FRDA patients. This response can be restored by AMPK and PPARgamma agonists, suggesting a potential therapeutic approach for FRDA.

Increased BDNF promoter methylation in the Wernicke area of suicide subjects.

CONTEXT: Brain-derived neurotrophic factor (BDNF) plays a pivotal role in the pathophysiology of suicidal behavior and BDNF levels are decreased in the brain and plasma of suicide subjects. So far, the mechanisms leading to downregulation of BDNF expression are poorly understood. OBJECTIVES: To test the hypothesis that alterations of DNA methylation could be involved in the dysregulation of BDNF gene expression in the brain of suicide subjects. DESIGN: Three independent quantitative methylation techniques were performed on postmortem samples of brain tissue. BDNF messenger RNA levels were determined by quantitative real-time polymerase chain reaction. SETTING: Academic medical center. PATIENTS OR OTHER PARTICIPANTS: Forty-four suicide completers and 33 nonsuicide control subjects of white ethnicity. MAIN OUTCOME MEASURES: The DNA methylation degree at BDNF promoter IV and the genome-wide DNA methylation levels in the brain's Wernicke area. RESULTS: Postmortem brain samples from suicide subjects showed a statistically significant increase of DNA methylation at specific CpG sites in BDNF promoter/exon IV compared with nonsuicide control subjects (P < .001). Most of the CpG sites lying in the -300/+500 region, on both strands, had low or no methylation, with the exception of a few sites located near the transcriptional start site that had differential methylation, while genome-wide methylation levels were comparable among the subjects. The mean methylation degree at the 4 CpG sites analyzed by pyrosequencing was always less than 12.9% in the 33 nonsuicide control subjects, while in 13 of 44 suicide victims (30%), the mean methylation degree ranged between 13.1% and 34.2%. Higher methylation degree corresponded to lower BDNF messenger RNA levels. CONCLUSIONS: BDNF promoter/exon IV is frequently hypermethylated in the Wernicke area of the postmortem brain of suicide subjects irrespective of genome-wide methylation levels, indicating that a gene-specific increase in DNA methylation could cause or contribute to the downregulation of BDNF expression in suicide subjects. The reported data reveal a novel link between epigenetic alteration in the brain and suicidal behavior.

Genome-wide scan for signatures of human population differentiation and their relationship with natural selection, functional pathways and diseases.

Genetic differences both between individuals and populations are studied for their evolutionary relevance and for their potential medical applications. Most of the genetic differentiation among populations are caused by random drift that should affect all loci across the genome in a similar manner. When a locus shows extraordinary high or low levels of population differentiation, this may be interpreted as evidence for natural selection. The most used measure of population differentiation was devised by Wright and is known as fixation index, or F(ST). We performed a genome-wide estimation of F(ST) on about 4 millions of SNPs from HapMap project data. We demonstrated a heterogeneous distribution of F(ST) values between autosomes and heterochromosomes. When we compared the F(ST) values obtained in this study with another evolutionary measure obtained by comparative interspecific approach, we found that genes under positive selection appeared to show low levels of population differentiation. We applied a gene set approach, widely used for microarray data analysis, to detect functional pathways under selection. We found that one pathway related to antigen processing and presentation showed low levels of F(ST), while several pathways related to cell signalling, growth and morphogenesis showed high F(ST) values. Finally, we detected a signature of selection within genes associated with human complex diseases. These results can help to identify which process occurred during human evolution and adaptation to different environments. They also support the hypothesis that common diseases could have a genetic background shaped by human evolution.