Predicting cell types and genetic variations contributing to disease by combining GWAS and epigenetic data

Genome-wide association studies (GWASs) identify single nucleotide polymorphisms (SNPs) that are enriched in individuals suffering from a given disease. Most disease-associated SNPs fall into non-coding regions, so that it is not straightforward to infer phenotype or function; moreover, many SNPs are in tight genetic linkage, so that a SNP identified as associated with a particular disease may not itself be causal, but rather signify the presence of a linked SNP that is functionally relevant to disease pathogenesis. Here, we present an analysis method that takes advantage of the recent rapid accumulation of epigenomics data to address these problems for some SNPs. Using asthma as a prototypic example; we show that non-coding disease-associated SNPs are enriched in genomic regions that function as regulators of transcription, such as enhancers and promoters. Identifying enhancers based on the presence of the histone modification marks such as H3K4me1 in different cell types, we show that the location of enhancers is highly cell-type specific. We use these findings to predict which SNPs are likely to be directly contributing to disease based on their presence in regulatory regions, and in which cell types their effect is expected to be detectable. Moreover, we can also predict which cell types contribute to a disease based on overlap of the disease-associated SNPs with the locations of enhancers present in a given cell type. Finally, we suggest that it will be possible to re-analyze GWAS studies with much higher power by limiting the SNPs considered to those in coding or regulatory regions of cell types relevant to a given disease

Mean first-passage times of non-Markovian random walkers in confinement

The first-passage time (FPT), defined as the time a random walker takes to reach a target point in a confining domain, is a key quantity in the theory of stochastic processes. Its importance comes from its crucial role to quantify the efficiency of processes as varied as diffusion-limited reactions, target search processes or spreading of diseases. Most methods to determine the FPT properties in confined domains have been limited to Markovian (memoryless) processes. However, as soon as the random walker interacts with its environment, memory effects can not be neglected. Examples of non Markovian dynamics include single-file diffusion in narrow channels or the motion of a tracer particle either attached to a polymeric chain or diffusing in simple or complex fluids such as nematics \cite{turiv2013effect}, dense soft colloids or viscoelastic solution. Here, we introduce an analytical approach to calculate, in the limit of a large confining volume, the mean FPT of a Gaussian non-Markovian random walker to a target point. The non-Markovian features of the dynamics are encompassed by determining the statistical properties of the trajectory of the random walker in the future of the first-passage event, which are shown to govern the FPT kinetics.This analysis is applicable to a broad range of stochastic processes, possibly correlated at long-times. Our theoretical predictions are confirmed by numerical simulations for several examples of non-Markovian processes including the emblematic case of the Fractional Brownian Motion in one or higher dimensions. These results show, on the basis of Gaussian processes, the importance of memory effects in first-passage statistics of non-Markovian random walkers in confinement.Comment: Submitted version. Supplementary Information can be found on the Nature website : http://www.nature.com/nature/journal/v534/n7607/full/nature18272.htm

Internalization of Formyl Peptide Receptor in Leukocytes Subject to Fluid Stresses

Human leukocytes retract pseudopods under normal physiologic levels of fluid shear stress even in the absence of any other mediator. To gain more detailed understanding of the mechanisms that regulate this cell behavior, we exposed leukocytes to a steady state laminar shear field in a flow chamber and computed the fluid stresses distribution on the surface of individual cells with and without pseudopod. The surface fluid stress distribution on such cell is quite inhomogeneous. We hypothesized that the local fluid stresses on the cell surface serve to regulate pseudopod retraction by way of membrane receptors, especially the formyl peptide receptor (FPR). Comparison of the receptor distribution and the stress distribution over the surface of the cells indicates that the membrane fluid stress alone is not directly correlated with the extent of regional pseudopod retraction, giving further support to the hypothesis that membrane receptors are involved in the mechanotransduction of leukocytes. We observed that after exposure to fluid shear the FPR was internalized to a small intracellular compartment. This internalization appears to be independent of the original location of the receptor on the surface of the cell and the FPR appears to be more derived from multiple locations on the cell, with both higher and lower fluid stresses. The evidence suggests that FPR involvement in the pseudopod-retraction process is not limited to cell surface regions with the highest fluid shear stress, but rather a more global occurrence over the majority of the cell membrane

pGQL: A probabilistic graphical query language for gene expression time courses

<p>Abstract</p> <p>Background</p> <p>Timeboxes are graphical user interface widgets that were proposed to specify queries on time course data. As queries can be very easily defined, an exploratory analysis of time course data is greatly facilitated. While timeboxes are effective, they have no provisions for dealing with noisy data or data with fluctuations along the time axis, which is very common in many applications. In particular, this is true for the analysis of gene expression time courses, which are mostly derived from noisy microarray measurements at few unevenly sampled time points. From a data mining point of view the robust handling of data through a sound statistical model is of great importance.</p> <p>Results</p> <p>We propose probabilistic timeboxes, which correspond to a specific class of Hidden Markov Models, that constitutes an established method in data mining. Since HMMs are a particular class of probabilistic graphical models we call our method Probabilistic Graphical Query Language. Its implementation was realized in the free software package pGQL. We evaluate its effectiveness in exploratory analysis on a yeast sporulation data set.</p> <p>Conclusions</p> <p>We introduce a new approach to define dynamic, statistical queries on time course data. It supports an interactive exploration of reasonably large amounts of data and enables users without expert knowledge to specify fairly complex statistical models with ease. The expressivity of our approach is by its statistical nature greater and more robust with respect to amplitude and frequency fluctuation than the prior, deterministic timeboxes.</p

Reply: An inverse association between tumour size and overdiagnosis may explain the results by Bucchi et al

..

A Platform for Processing Expression of Short Time Series (PESTS)

<p>Abstract</p> <p>Background</p> <p>Time course microarray profiles examine the expression of genes over a time domain. They are necessary in order to determine the complete set of genes that are dynamically expressed under given conditions, and to determine the interaction between these genes. Because of cost and resource issues, most time series datasets contain less than 9 points and there are few tools available geared towards the analysis of this type of data.</p> <p>Results</p> <p>To this end, we introduce a platform for Processing Expression of Short Time Series (PESTS). It was designed with a focus on usability and interpretability of analyses for the researcher. As such, it implements several standard techniques for comparability as well as visualization functions. However, it is designed specifically for the unique methods we have developed for significance analysis, multiple test correction and clustering of short time series data. The central tenet of these methods is the use of biologically relevant features for analysis. Features summarize short gene expression profiles, inherently incorporate dependence across time, and allow for both full description of the examined curve and missing data points.</p> <p>Conclusions</p> <p>PESTS is fully generalizable to other types of time series analyses. PESTS implements novel methods as well as several standard techniques for comparability and visualization functions. These features and functionality make PESTS a valuable resource for a researcher's toolkit. PESTS is available to download for free to academic and non-profit users at <url>http://www.mailman.columbia.edu/academic-departments/biostatistics/research-service/software-development</url>.</p

Interaction between genetic and epigenetic variation defines gene expression patterns at the asthma-associated locus 17q12-q21 in lymphoblastoid cell lines

Phenotypic variation results from variation in gene expression, which is modulated by genetic and/or epigenetic factors. To understand the molecular basis of human disease, interaction between genetic and epigenetic factors needs to be taken into account. The asthma-associated region 17q12-q21 harbors three genes, the zona pellucida binding protein 2 (ZPBP2), gasdermin B (GSDMB) and ORM1-like 3 (ORMDL3), that show allele-specific differences in expression levels in lymphoblastoid cell lines (LCLs) and CD4+ T cells. Here, we report a molecular dissection of allele-specific transcriptional regulation of the genes within the chromosomal region 17q12-q21 combining in vitro transfection, formaldehyde-assisted isolation of regulatory elements, chromatin immunoprecipitation and DNA methylation assays in LCLs. We found that a single nucleotide polymorphism rs4795397 influences the activity of ZPBP2 promoter in vitro in an allele-dependent fashion, and also leads to nucleosome repositioning on the asthma-associated allele. However, variable methylation of exon 1 of ZPBP2 masks the strong genetic effect on ZPBP2 promoter activity in LCLs. In contrast, the ORMDL3 promoter is fully unmethylated, which allows detection of genetic effects on its transcription. We conclude that the cis-regulatory effects on 17q12-q21 gene expression result from interaction between several regulatory polymorphisms and epigenetic factors within the cis-regulatory haplotype region

SUBMIT: Systemic therapy with or without up front surgery of the primary tumor in breast cancer patients with distant metastases at initial presentation

<p>Abstract</p> <p>Background</p> <p>Five percent of all patients with breast cancer have distant metastatic disease at initial presentation. Because metastatic breast cancer is considered to be an incurable disease, it is generally treated with a palliative intent. Recent non-randomized studies have demonstrated that (complete) resection of the primary tumor is associated with a significant improvement of the survival of patients with primary metastatic breast cancer. However, other studies have suggested that the claimed survival benefit by surgery may be caused by selection bias. Therefore, a randomized controlled trial will be performed to assess whether breast surgery in patients with primary distant metastatic breast cancer will improve the prognosis.</p> <p>Design</p> <p>Randomization will take place after the diagnosis of primary distant metastatic breast cancer. Patients will either be randomized to up front surgery of the breast tumor followed by systemic therapy or to systemic therapy, followed by delayed local treatment of the breast tumor if clinically indicated.</p> <p>Patients with primary distant metastatic breast cancer, with no prior treatment of the breast cancer, who are 18 years or older and fit enough to undergo surgery and systemic therapy are eligible. Important exclusion criteria are: prior invasive breast cancer, surgical treatment or radiotherapy of this breast tumor before randomization, irresectable T4 tumor and synchronous bilateral breast cancer. The primary endpoint is 2-year survival. Quality of life and local tumor control are among the secondary endpoints.</p> <p>Based on the results of prior research it was calculated that 258 patients are needed in each treatment arm, assuming a power of 80%. Total accrual time is expected to take 60 months. An interim analysis will be performed to assess any clinically significant safety concerns and to determine whether there is evidence that up front surgery is clinically or statistically inferior to systemic therapy with respect to the primary endpoint.</p> <p>Discussion</p> <p>The SUBMIT study is a randomized controlled trial that will provide evidence on whether or not surgery of the primary tumor in breast cancer patients with metastatic disease at initial presentation results in an improved survival.</p> <p>Trial registration</p> <p><a href="http://www.clinicaltrials.gov/ct2/show/NCT01392586">NCT01392586</a>.</p

Mastectomy rates are decreasing in the era of service screening: a population-based study in Italy (1997–2001)

We enrolled all 2162 in situ and 21 148 invasive cases of breast cancer in 17 areas of Italy, diagnosed in 1997–2001. Rates of early cancer increased by 13.7% in the screening age group (50–69 years), and breast conserving surgery by 24.6%. Advanced cancer rates decreased by 19.4%, and mastectomy rates by 24.2%. Service screening did not increase mastectomy rates in the study population

Conditional embryonic lethality to improve the sterile insect technique in Ceratitis capitata (Diptera: Tephritidae)

<p>Abstract</p> <p>Background</p> <p>The sterile insect technique (SIT) is an environment-friendly method used in area-wide pest management of the Mediterranean fruit fly <it>Ceratitis capitata </it>(Wiedemann; Diptera: Tephritidae). Ionizing radiation used to generate reproductive sterility in the mass-reared populations before release leads to reduction of competitiveness.</p> <p>Results</p> <p>Here, we present a first alternative reproductive sterility system for medfly based on transgenic embryonic lethality. This system is dependent on newly isolated medfly promoter/enhancer elements of cellularization-specifically-expressed genes. These elements act differently in expression strength and their ability to drive lethal effector gene activation. Moreover, position effects strongly influence the efficiency of the system. Out of 60 combinations of driver and effector construct integrations, several lines resulted in larval and pupal lethality with one line showing complete embryonic lethality. This line was highly competitive to wildtype medfly in laboratory and field cage tests.</p> <p>Conclusion</p> <p>The high competitiveness of the transgenic lines and the achieved 100% embryonic lethality causing reproductive sterility without the need of irradiation can improve the efficacy of operational medfly SIT programs.</p