penalizedclr: an R package for penalized conditional logistic regression for integration of multiple omics layers.

BACKGROUND: The matched case-control design, up until recently mostly pertinent to epidemiological studies, is becoming customary in biomedical applications as well. For instance, in omics studies, it is quite common to compare cancer and healthy tissue from the same patient. Furthermore, researchers today routinely collect data from various and variable sources that they wish to relate to the case-control status. This highlights the need to develop and implement statistical methods that can take these tendencies into account. RESULTS: We present an R package penalizedclr, that provides an implementation of the penalized conditional logistic regression model for analyzing matched case-control studies. It allows for different penalties for different blocks of covariates, and it is therefore particularly useful in the presence of multi-source omics data. Both L1 and L2 penalties are implemented. Additionally, the package implements stability selection for variable selection in the considered regression model. CONCLUSIONS: The proposed method fills a gap in the available software for fitting high-dimensional conditional logistic regression models accounting for the matched design and block structure of predictors/features. The output consists of a set of selected variables that are significantly associated with case-control status. These variables can then be investigated in terms of functional interpretation or validation in further, more targeted studies.

ScreenDMT reveals linoleic acid diols replicably associate with BMI and stimulate adipocyte calcium fluxes.

Activating brown adipose tissue (BAT) improves systemic metabolism, making it a promising target for metabolic syndrome. BAT is activated by 12, 13-dihydroxy-9Z-octadecenoic acid (12, 13-diHOME), which we previously identified to be inversely associated with BMI and which directly improves metabolism in multiple tissues. Here we profile plasma lipidomics from a cohort of 83 people and test which lipids' association with BMI replicates in a concordant direction using our novel tool ScreenDMT, whose power and validity we demonstrate via mathematical proofs and simulations. We find that the linoleic acid diols 12, 13-diHOME and 9, 10-diHOME both replicably inversely associate with BMI and mechanistically activate calcium fluxes in mouse brown and white adipocytes in vitro, which implicates this pathway and 9, 10-diHOME as candidate therapeutic targets. ScreenDMT can be applied to test directional mediation, directional replication, and qualitative interactions, such as identifying biomarkers whose association is shared (replication) or opposite (qualitative interaction) across diverse populations.

A Bartlett-type correction for likelihood ratio tests with application to testing equality of Gaussian graphical models.

This work defines a new correction for the likelihood ratio test for a two-sample problem within the multivariate normal context. This correction applies to decomposable graphical models, where testing equality of distributions can be decomposed into lower dimensional problems.

Pre-diagnostic DNA methylation in blood leucocytes in cutaneous melanoma; a nested case-control study within the Norwegian Women and Cancer cohort.

The prognosis of cutaneous melanoma depends on early detection, and good biomarkers for melanoma risk may provide a valuable tool to detect melanoma development at a pre-clinical stage. By studying the epigenetic profile in pre-diagnostic blood samples of melanoma cases and cancer free controls, we aimed to identify DNA methylation sites conferring melanoma risk. DNA methylation was measured at 775,528 CpG sites using the Illumina EPIC array in whole blood in incident melanoma cases (n = 183) and matched cancer-free controls (n = 183) in the Norwegian Women and Cancer cohort. Phenotypic information and ultraviolet radiation exposure were obtained from questionnaires. Epigenome wide association (EWAS) was analyzed in future melanoma cases and controls with conditional logistic regression, with correction for multiple testing using the false discovery rate (FDR). We extended the analysis by including a public data set on melanoma (GSE120878), and combining these different data sets using a version of covariate modulated FDR (AdaPT). The analysis on future melanoma cases and controls did not identify any genome wide significant CpG sites (0.85 ≤ padj ≤ 0.99). In the restricted AdaPT analysis, 7 CpG sites were suggestive at the FDR level of 0.15. These CpG sites may potentially be used as pre-diagnostic biomarkers of melanoma risk.

On optimal two-stage testing of multiple mediators.

Mediation analysis in high-dimensional settings often involves identifying potential mediators among a large number of measured variables. For this purpose, a two-step familywise error rate procedure called ScreenMin has been recently proposed. In ScreenMin, variables are first screened and only those that pass the screening are tested. The proposed data-independent threshold for selection has been shown to guarantee asymptotic familywise error rate. In this work, we investigate the impact of the threshold on the finite-sample familywise error rate. We derive a power maximizing threshold and show that it is well approximated by an adaptive threshold of Wang et al. (2016, arXiv preprint arXiv:1610.03330). We illustrate the investigated procedures on a case-control study examining the effect of fish intake on the risk of colorectal adenoma. We also apply our procedure in the context of replicability analysis to identify single nucleotide polymorphisms (SNP) associated with crop yield in two distinct environments.

High-throughput mediation analysis of human proteome and metabolome identifies mediators of post-bariatric surgical diabetes control.

To improve the power of mediation in high-throughput studies, here we introduce High-throughput mediation analysis (Hitman), which accounts for direction of mediation and applies empirical Bayesian linear modeling. We apply Hitman in a retrospective, exploratory analysis of the SLIMM-T2D clinical trial in which participants with type 2 diabetes were randomized to Roux-en-Y gastric bypass (RYGB) or nonsurgical diabetes/weight management, and fasting plasma proteome and metabolome were assayed up to 3 years. RYGB caused greater improvement in HbA1c, which was mediated by growth hormone receptor (GHR). GHR's mediation is more significant than clinical mediators, including BMI. GHR decreases at 3 months postoperatively alongside increased insulin-like growth factor binding proteins IGFBP1/BP2; plasma GH increased at 1 year. Experimental validation indicates (1) hepatic GHR expression decreases in post-bariatric rats; (2) GHR knockdown in primary hepatocytes decreases gluconeogenic gene expression and glucose production. Thus, RYGB may induce resistance to diabetogenic effects of GH signaling.Trial Registration: Clinicaltrials.gov NCT01073020.

Leveraging three-dimensional chromatin architecture for effective reconstruction of enhancer-target gene regulatory interactions.

A growing amount of evidence in literature suggests that germline sequence variants and somatic mutations in non-coding distal regulatory elements may be crucial for defining disease risk and prognostic stratification of patients, in genetic disorders as well as in cancer. Their functional interpretation is challenging because genome-wide enhancer-target gene (ETG) pairing is an open problem in genomics. The solutions proposed so far do not account for the hierarchy of structural domains which define chromatin three-dimensional (3D) architecture. Here we introduce a change of perspective based on the definition of multi-scale structural chromatin domains, integrated in a statistical framework to define ETG pairs. In this work (i) we develop a computational and statistical framework to reconstruct a comprehensive map of ETG pairs leveraging functional genomics data; (ii) we demonstrate that the incorporation of chromatin 3D architecture information improves ETG pairing accuracy and (iii) we use multiple experimental datasets to extensively benchmark our method against previous solutions for the genome-wide reconstruction of ETG pairs. This solution will facilitate the annotation and interpretation of sequence variants in distal non-coding regulatory elements. We expect this to be especially helpful in clinically oriented applications of whole genome sequencing in cancer and undiagnosed genetic diseases research.

Lifetime Ultraviolet Radiation Exposure and DNA Methylation in Blood Leukocytes: The Norwegian Women and Cancer Study.

Ultraviolet radiation (UVR) exposure is a leading cause of skin cancers and an ubiquitous environmental exposure. However, the molecular mechanisms relating UVR exposure to melanoma is not fully understood. We aimed to investigate if lifetime UVR exposure could be robustly associated to DNA methylation (DNAm). We assessed DNAm in whole blood in three data sets (n = 183, 191, and 125) from the Norwegian Woman and Cancer cohort, using Illumina platforms. We studied genome-wide DNAm, targeted analyses of CpG sites indicated in the literature, global methylation, and accelerated aging. Lifetime history of UVR exposure (residential ambient UVR, sunburns, sunbathing vacations and indoor tanning) was collected by questionnaires. We used one data set for discovery and the other two for replication. One CpG site showed a genome-wide significant association to cumulative UVR exposure (cg01884057) (pnominal = 3.96e-08), but was not replicated in any of the two replication sets (pnominal ≥ 0.42). Two CpG sites (cg05860019, cg00033666) showed suggestive associations with the other UVR exposures. We performed extensive analyses of the association between long-term UVR exposure and DNAm. There was no indication of a robust effect of past UVR exposure on DNAm.

SourceSet: A graphical model approach to identify primary genes in perturbed biological pathways.

Topological gene-set analysis has emerged as a powerful means for omic data interpretation. Although numerous methods for identifying dysregulated genes have been proposed, few of them aim to distinguish genes that are the real source of perturbation from those that merely respond to the signal dysregulation. Here, we propose a new method, called SourceSet, able to distinguish between the primary and the secondary dysregulation within a Gaussian graphical model context. The proposed method compares gene expression profiles in the control and in the perturbed condition and detects the differences in both the mean and the covariance parameters with a series of likelihood ratio tests. The resulting evidence is used to infer the primary and the secondary set, i.e. the genes responsible for the primary dysregulation, and the genes affected by the perturbation through network propagation. The proposed method demonstrates high specificity and sensitivity in different simulated scenarios and on several real biological case studies. In order to fit into the more traditional pathway analysis framework, SourceSet R package also extends the analysis from a single to multiple pathways and provides several graphical outputs, including Cytoscape visualization to browse the results.

Global test for high-dimensional mediation: Testing groups of potential mediators.

We address the problem of testing whether a possibly high-dimensional vector may act as a mediator between some exposure variable and the outcome of interest. We propose a global test for mediation, which combines a global test with the intersection-union principle. We discuss theoretical properties of our approach and conduct simulation studies that demonstrate that it performs equally well or better than its competitor. We also propose a multiple testing procedure, ScreenMin, that provides asymptotic control of either familywise error rate or false discovery rate when multiple groups of potential mediators are tested simultaneously. We apply our approach to data from a large Norwegian cohort study, where we look at the hypothesis that smoking increases the risk of lung cancer by modifying the level of DNA methylation.

simPATHy: a new method for simulating data from perturbed biological PATHways.

Summary: In the omic era, one of the main aims is to discover groups of functionally related genes that drive the difference between different conditions. To this end, a plethora of potentially useful multivariate statistical approaches has been proposed, but their evaluation is hindered by the absence of a gold standard. Here, we propose a method for simulating biological data ­ gene expression, RPKM/FPKM or protein abundances ­ from two conditions, namely, a reference condition and a perturbation of it. Our approach is built upon probabilistic graphical models and is thus especially suited for testing topological approaches. Availability and Implementation: The simPATHy is an R package, it is open source and freely available on CRAN. Contacts: elisa.salviato.2@studenti.unipd.it or chiara.romualdi@unipd.it Supplementary Information: Supplementary data are available at Bioinformatics online.

The role of the frontal aslant tract and premotor connections in visually guided hand movements.

Functional neuroimaging and brain lesion studies demonstrate that secondary motor areas of the frontal lobe play a crucial role in the cortical control of hand movements. However, no study so far has examined frontal white matter connections of the secondary motor network, namely the frontal aslant tract, connecting the supplementary motor complex and the posterior inferior frontal regions, and the U-shaped dorsal and ventral premotor fibers running through the middle frontal gyrus. The aim of the current study is to explore the involvement of the short frontal lobe connections in reaching and reach-to-grasp movements in 32 right-handed healthy subjects by correlating tractography data based on spherical deconvolution approach with kinematical data. We showed that individual differences in the microstructure of the bilateral frontal aslant tract, bilateral ventral and left dorsal premotor tracts were associated with kinematic features of hand actions. Furthermore, bilateral ventral premotor connections were also involved in the closing grip phase necessary for determining efficient and stable grasping of the target object. This work suggests for the first time that hand kinematics and visuomotor processing are associated with the anatomy of the short frontal lobe connections.

Graphical modeling for gene set analysis: A critical appraisal.

Current demand for understanding the behavior of groups of related genes, combined with the greater availability of data, has led to an increased focus on statistical methods in gene set analysis. In this paper, we aim to perform a critical appraisal of the methodology based on graphical models developed in Massa et al. (2010) that uses pathway signaling networks as a starting point to develop statistically sound procedures for gene set analysis. We pay attention to the potential of the methodology with respect to the organizational aspects of dealing with such complex but highly informative starting structures, that is pathways. We focus on three themes: the translation of a biological pathway into a graph suitable for modeling, the role of shrinkage when more genes than samples are obtained, the evaluation of respondence of the statistical models to the biological expectations. To study the impact of shrinkage, two simulation studies will be run. To evaluate the biological expectation we will use data from a network with known behavior that offer the possibility of carrying out a realistic check of respondence of the model to changes in the experimental conditions.