David Gomez-Cabrero

A beta-mixture quantile normalization method for correcting probe design bias in Illumina Infinium 450 k DNA methylation data

Motivation: The Illumina Infinium 450 k DNA Methylation Beadchip is a prime candidate technology for Epigenome-Wide Association Studies (EWAS). However, a difficulty associated with these beadarrays is that probes come in two different designs, characterized by widely different DNA methylation distributions and dynamic range, which may bias downstream analyses. A key statistical issue is therefore how best to adjust for the two different probe designs. Results: Here we propose a novel model-based intra-array normalization strategy for 450 k data, called BMIQ (Beta MIxture Quantile dilation), to adjust the beta-values of type2 design probes into a statistical distribution characteristic of type1 probes. The strategy involves application of a three-state beta-mixture model to assign probes to methylation states, subsequent transformation of probabilities into quantiles and finally a methylation-dependent dilation transformation to preserve the monotonicity and continuity of the data. We validate our method on cell-line data, fresh frozen and paraffin-embedded tumour tissue samples and demonstrate that BMIQ compares favourably with two competing methods. Specifically, we show that BMIQ improves the robustness of the normalization procedure, reduces the technical variation and bias of type2 probe values and successfully eliminates the type1 enrichment bias caused by the lower dynamic range of type2 probes. BMIQ will be useful as a preprocessing step for any study using the Illumina Infinium 450 k platform. Availability: BMIQ is freely available from http://code.google.com/p/bmiq/. Contact: a.teschendorff@ucl.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online

A survey of best practices for RNA-seq data analysis

RNA-sequencing (RNA-seq) has a wide variety of applications, but no single analysis pipeline can be used in all cases. We review all of the major steps in RNA-seq data analysis, including experimental design, quality control, read alignment, quantification of gene and transcript levels, visualization, differential gene expression, alternative splicing, functional analysis, gene fusion detection and eQTL mapping. We highlight the challenges associated with each step. We discuss the analysis of small RNAs and the integration of RNA-seq with other functional genomics techniques. Finally, we discuss the outlook for novel technologies that are changing the state of the art in transcriptomics.

Data integration in the era of omics: current and future challenges

To integrate heterogeneous and large omics data constitutes not only a conceptual challenge but a practical hurdle in the daily analysis of omics data. With the rise of novel omics technologies and through large-scale consortia projects, biological systems are being further investigated at an unprecedented scale generating heterogeneous and often large data sets. These data-sets encourage researchers to develop novel data integration methodologies. In this introduction we review the definition and characterize current efforts on data integration in the life sciences. We have used a web-survey to assess current research projects on data-integration to tap into the views, needs and challenges as currently perceived by parts of the research community.

An evaluation of analysis pipelines for DNA methylation profiling using the Illumina HumanMethylation450 BeadChip platform

The proper identification of differentially methylated CpGs is central in most epigenetic studies. The Illumina HumanMethylation450 BeadChip is widely used to quantify DNA methylation; nevertheless, the design of an appropriate analysis pipeline faces severe challenges due to the convolution of biological and technical variability and the presence of a signal bias between Infinium I and II probe design types. Despite recent attempts to investigate how to analyze DNA methylation data with such an array design, it has not been possible to perform a comprehensive comparison between different bioinformatics pipelines due to the lack of appropriate data sets having both large sample size and sufficient number of technical replicates. Here we perform such a comparative analysis, targeting the problems of reducing the technical variability, eliminating the probe design bias and reducing the batch effect by exploiting two unpublished data sets, which included technical replicates and were profiled for DNA methylation either on peripheral blood, monocytes or muscle biopsies. We evaluated the performance of different analysis pipelines and demonstrated that: (1) it is critical to correct for the probe design type, since the amplitude of the measured methylation change depends on the underlying chemistry; (2) the effect of different normalization schemes is mixed, and the most effective method in our hands were quantile normalization and Beta Mixture Quantile dilation (BMIQ); (3) it is beneficial to correct for batch effects. In conclusion, our comparative analysis using a comprehensive data set suggests an efficient pipeline for proper identification of differentially methylated CpGs using the Illumina 450K arrays.

Identification of novel markers in rheumatoid arthritis through integrated analysis of DNA methylation and microRNA expression

Autoimmune rheumatic diseases are complex disorders, whose etiopathology is attributed to a crosstalk between genetic predisposition and environmental factors. Both variants of autoimmune susceptibility genes and environment are involved in the generation of aberrant epigenetic profiles in a cell-specific manner, which ultimately result in dysregulation of expression. Furthermore, changes in miRNA expression profiles also cause gene dysregulation associated with aberrant phenotypes. In rheumatoid arthritis, several cell types are involved in the destruction of the joints, synovial fibroblasts being among the most important. In this study we performed DNA methylation and miRNA expression screening of a set of rheumatoid arthritis synovial fibroblasts and compared the results with those obtained from osteoarthritis patients with a normal phenotype. DNA methylation screening allowed us to identify changes in novel key target genes like IL6R, CAPN8 and DPP4, as well as several HOX genes. A significant proportion of genes undergoing DNA methylation changes were inversely correlated with expression. miRNA screening revealed the existence of subsets of miRNAs that underwent changes in expression. Integrated analysis highlighted sets of miRNAs that are controlled by DNA methylation, and genes that are regulated by DNA methylation and are targeted by miRNAs with a potential use as clinical markers. Our study enabled the identification of novel dysregulated targets in rheumatoid arthritis synovial fibroblasts and generated a new workflow for the integrated analysis of miRNA and epigenetic control.

An integrative analysis reveals coordinated reprogramming of the epigenome and the transcriptome in human skeletal muscle after training

Regular endurance exercise training induces beneficial functional and health effects in human skeletal muscle. The putative contribution to the training response of the epigenome as a mediator between genes and environment has not been clarified. Here we investigated the contribution of DNA methylation and associated transcriptomic changes in a well-controlled human intervention study. Training effects were mirrored by significant alterations in DNA methylation and gene expression in regions with a homogeneous muscle energetics and remodeling ontology. Moreover, a signature of DNA methylation and gene expression separated the samples based on training and gender. Differential DNA methylation was predominantly observed in enhancers, gene bodies and intergenic regions and less in CpG islands or promoters. We identified transcriptional regulator binding motifs of MRF, MEF2 and ETS proteins in the proximity of the changing sites. A transcriptional network analysis revealed modules harboring distinct ontologies and, interestingly, the overall direction of the changes of methylation within each module was inversely correlated to expression changes. In conclusion, we show that highly consistent and associated modifications in methylation and expression, concordant with observed health-enhancing phenotypic adaptations, are induced by a physiological stimulus.

Cesarean delivery and hematopoietic stem cell epigenetics in the newborn infant: implications for future health?

OBJECTIVE Cesarean section (CS) has been associated with a greater risk for asthma, diabetes, and cancer later in life. Although elective CS continues to rise, it is unclear whether and how it may contribute to compromised future health. Our aim was to investigate the influence of mode of delivery on the epigenetic state in neonatal hematopoietic stem cells. STUDY DESIGN This was an observational study of 64 healthy, singleton, newborn infants (33 boys) born at term. Cord blood was sampled after elective CS (n = 27) and vaginal delivery. Global deoxyribonucleic acid (DNA) methylation in hematopoietic stem cells (CD34+) was determined by luminometric methylation assay, and genome-wide, locus-specific DNA methylation analysis was performed by Illumina Infinium 450K (Illumina, San Diego, CA), validated by bisulfite-pyrosequencing. RESULTS CD34+ cells from infants delivered by CS were globally more DNA methylated (+2%) than DNA from infants delivered vaginally (P = .02). In relation to mode of delivery, a locus-specific analysis identified 343 loci with a difference in DNA methylation of 10% or greater (P < .01). A majority of the differentially methylated loci in neonatal CD34+ cells (76%) were found to be hypermethylated after vaginal delivery. In these infants, the degree of DNA methylation in 3 loci correlated to the duration of labor. The functional relevance of differentially methylated loci involved processes such as immunoglobulin biosynthetic process, regulation of glycolysis and ketone metabolism, and regulation of the response to food. CONCLUSION A possible interpretation is that mode of delivery affects the epigenetic state of neonatal hematopoietic stem cells. Given the functional relevance indicated, our findings may have important implications for health and disease in later life.

Reconciling coherent oscillation with modulation of irregular spiking activity in selective attention: gamma-range synchronization between sensory and executive cortical areas.

In this computational work, we investigated gamma-band synchronization across cortical circuits associated with selective attention. The model explicitly instantiates a reciprocally connected loop of spiking neurons between a sensory-type (area MT) and an executive-type (prefrontal/parietal) cortical circuit (the source area for top-down attentional signaling). Moreover, unlike models in which neurons behave as clock-like oscillators, in our model single-cell firing is highly irregular (close to Poisson), while local field potential exhibits a population rhythm. In this “sparsely synchronized oscillation” regime, the model reproduces and clarifies multiple observations from behaving animals. Top-down attentional inputs have a profound effect on network oscillatory dynamics while only modestly affecting single-neuron spiking statistics. In addition, attentional synchrony modulations are highly selective: interareal neuronal coherence occurs only when there is a close match between the preferred feature of neurons, the attended feature, and the presented stimulus, a prediction that is experimentally testable. When interareal coherence was abolished, attention-induced gain modulations of sensory neurons were slightly reduced. Therefore, our model reconciles the rate and synchronization effects, and suggests that interareal coherence contributes to large-scale neuronal computation in the brain through modest enhancement of rate modulations as well as a pronounced attention-specific enhancement of neural synchrony.

Pre-B cell to macrophage transdifferentiation without significant promoter DNA methylation changes

Transcription factor-induced lineage reprogramming or transdifferentiation experiments are essential for understanding the plasticity of differentiated cells. These experiments helped to define the specific role of transcription factors in conferring cell identity and played a key role in the development of the regenerative medicine field. We here investigated the acquisition of DNA methylation changes during C/EBPα-induced pre-B cell to macrophage transdifferentiation. Unexpectedly, cell lineage conversion occurred without significant changes in DNA methylation not only in key B cell- and macrophage-specific genes but also throughout the entire set of genes differentially methylated between the two parental cell types. In contrast, active and repressive histone modification marks changed according to the expression levels of these genes. We also demonstrated that C/EBPα and RNA Pol II are associated with the methylated promoters of macrophage-specific genes in reprogrammed macrophages without inducing methylation changes. Our findings not only provide insights about the extent and hierarchy of epigenetic events in pre-B cell to macrophage transdifferentiation but also show an important difference to reprogramming towards pluripotency where promoter DNA demethylation plays a pivotal role.

Serotonin Regulates Performance Nonmonotonically in a Spatial Working Memory Network

The prefrontal cortex (PFC) contains a dense network of serotonergic [serotonin, 5-hydroxytryptamine (5-HT)] axons, and endogenous 5-HT markedly modulates PFC neuronal function via several postsynaptic receptors. The therapeutic action of atypical antipsychotic drugs, acting mainly via 5-HT receptors, also suggests a role for serotonergic neurotransmission in cognitive functions. However, psychopharmacological studies have failed to find a consistent relationship between serotonergic transmission and cognitive functions of the PFC, including spatial working memory (SWM). Here, we built a computational network model to investigate 5-HT modulation of SWM in the PFC. We found that 5-HT modulates networks SWM performance nonmonotonically via 5-HT1A and 5-HT2A receptors, following an inverted U-shape. This relationship may contribute to blur the effects of serotonergic agents in previous SWM group-based behavioral studies. Our simulations also showed that errors occurring at low and high 5-HT concentrations are due to different network dynamics instabilities, suggesting that these 2 conditions can be distinguished experimentally based on their distinct dependency on experimental variables. We inferred specific predictions regarding the expected behavioral effects of serotonergic agents in 2 classic working-memory tasks. Our results underscore the relevance of identifying different error types in SWM tasks in order to reveal the association between neuromodulatory systems and SWM.