36 research outputs found
Effect of non-linear transformation on interaction effects.
<p>We defined an outcome <i>Y</i> as a function of a single nucleotide polymorphism <i>G</i> with a minor allele frequency of 0.1, an exposure <i>E</i> normally distributed with mean 5 and variance 1, and a right-skewed normal distributed residual term <i>ε</i>. In the framework of this analysis, TF mRNA level is considered as an exposure <i>E</i>. We generated two datasets of 10,000 individuals for the two scenarios. In a) <i>G</i> and <i>E</i> have only main effects and each explain 20% of the variance of <i>Y</i>. In b) <i>G</i> and <i>E</i> main effects each explain 10% of the outcome variance, but also have an interaction effect explaining 20% of the variance of <i>Y</i>. Upper panels show <i>Y</i> as a function of <i>E</i> by genotypic class and trend slope from a standard linear regression. Lower panels show the same data plotted after a rank-normal transformation (<i>rkt</i>) of <i>Y</i>. Interaction effect (observed as differences in slope by genotypic class) appears or disappears depending on the transformation applied to <i>Y</i>. <i>P</i>-values for interaction are indicated in red.</p
Distribution of interaction test lambdaGC in ECLIPSE.
<p>We derived the genomic inflation factor (<i>λ</i><sub><i>GC</i></sub>) of the standard interaction test using across sub-groups stratified based on P<sub><i>TF</i>.<i>marg</i></sub>, the <i>p</i>-value for association between the target gene and the candidate transcription factors (TFs). Grey bars present the total number of interaction tests falling in each strata. Four approaches were performed: i) no normal rank-transformation of the expression data (<i>std</i>), ii) HC3 correction of the effect estimate variance to account for heteroscedasticity (<i>h</i>3), iii) normal rank-transformation of expression data (<i>rkt</i>), and iv) HC3 correction and normal rank-transformation of expression data (<i>rkt</i>.<i>h</i>3).</p
Gene expression data pre-processing pipeline.
<p>Standard pre-processing methods applied to gene expression data prior to expression quantitative trait locus analysis. Note that alternative strategies are also used. For example step 2 is sometimes skipped and confounding factors (e.g. batch) are included in the model tested as covariates. Others have also applied step 3 before step 2.</p
Robustness comparison.
<p>QQplots over series of 8 million replicates where an outcome <i>Y</i> is simulated as a function of a genetic variant <i>G</i>, an unmeasured exposure <i>E</i>, an interaction between <i>G</i> and <i>E</i>, and in 50% of the replicates a measured exposure <i>Z</i>. In the framework of this analysis, <i>Z</i> and <i>E</i> are considered as measured and unmeasured TF mRNA level, respectively. The validity of five tests is evaluated by comparing the observed -log<sub>10</sub> (p-value) against the expected -log<sub>10</sub> (p-value) when testing for the null interaction between a <i>G</i> and <i>Z</i>. The tests include a standard linear regression using main and interaction terms only (STD), heteroscedasticity consistent-based tests using effect estimates from STD (HC0 and HC3), linear regression using binary-transformed <i>Z</i> (BIN), and a saturated model including a main effect of <i>Z</i><sup>2</sup> and each genotype coded as dummy variable (SAT). We considered coded allele frequency (CAF) of 0.05 (first row), 0.3 (middle row) and 0.5 (bottom row), and sample size <i>N</i> of 100, 500, 1,000 and 5,000. We randomly draw <i>E</i>, <i>Z</i>, and <i>ε</i>, the residual of <i>Y</i> from either a normal or a right-skewed normal distribution. For each scenario we derived the genomic inflation factor <i>λ</i><sub><i>GC</i></sub>.</p
When a true interaction can bias interaction screening.
<p>A quantitative outcome <i>Y</i> is defined as a linear function of a SNP <i>G</i>, an unmeasured exposure <i>E</i>, a measured exposure <i>Z</i>, and an interaction between <i>G</i> by <i>E</i>, with effect <i>γ</i><sub><i>G</i></sub>, <i>γ</i><sub><i>E</i></sub>, <i>γ</i><sub><i>Z</i></sub>, and <i>γ</i><sub><i>GE</i></sub>, respectively (as defined in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0173847#pone.0173847.e001" target="_blank">Eq 1</a>). All predictors were standardized to have mean 0 and variance 1. In the framework of this analysis, TF mRNA level is considered as an exposure <i>E</i>. We vary <i>γ</i><sub><i>GE</i></sub> so that the interaction term explains between 0 and 30% of the variance of <i>Y</i>. For simplicity we assume that, when relevant, the main effect of either <i>G</i>, <i>E</i>, or <i>Z</i> explains the same amount of variance as the interaction effect and set <i>ε</i> so that the variance of <i>Y</i> equals 1. Using this model we simulated series of 10,000 replicates, each including 400 individuals and tested for interaction between <i>G</i> and <i>Z</i> using a model not including the unmeasured exposure <i>E</i> (as defined in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0173847#pone.0173847.e003" target="_blank">Eq 2</a>), in the absence of main effect of the predictors (<i>γ</i><sub><i>G</i></sub> = <i>γ</i><sub><i>E</i></sub> = <i>γ</i><sub><i>Z</i></sub> = 0), panel <i>a</i>) or when including a main effect of <i>G</i> (<i>γ</i><sub><i>G</i></sub> ≠0, panel <i>b</i>), a main effect of <i>E</i> (<i>γ</i><sub><i>E</i></sub> ≠0), panel <i>b</i>), or a main effect of <i>G</i> (<i>γ</i><sub><i>Z</i></sub> ≠0, panel <i>d</i>). Upper panels show the increase in the residual variance of the outcome <i>δ</i> minus <i>ε</i> (so that models are comparable) stratified by genotypic class while increasing the interaction effect <i>γ</i><sub><i>GE</i></sub>. Lower panels show the type I error rate <i>α</i> at a <i>p</i>-value threshold of 0.05 for the interaction tests between <i>G</i> and <i>Z</i> derived for each series of 10,000 replicates.</p
Top 5 interaction signals for four different analytical strategies.
<p>Top 5 interaction signals for four different analytical strategies.</p
Characteristics of common variants identified across 6 blood pressure loci.
<p>Characteristics of common variants identified across 6 blood pressure loci.</p
Forest plot for rs12806040 in relation to DBP.
<p>The beta estimate and standard error are reported from the unconditional weighted analyses. P-values are reported from unconditional unweighted analyses.</p
Characteristics of the study participants.
<p>*PP calculated as SBP-DBP.</p>†<p>MAP calculated as (DBP+(SBP−DBP)/3).</p><p>All BP measures in the table are not adjusted for use of anti-hypertensive medication.</p><p>Characteristics of the study participants.</p
Forest plot for rs11191416 in relation to SBP.
<p>The beta estimate and standard error are reported from the unconditional weighted analyses. P-values are reported from unconditional unweighted analyses.</p