|
|
Genes with expression level regulated by genotype (N2 vs CB4856) and age at L3 larva and Late reproduction stage (96 hours at 24 centigrade). |
For model 2, authors used 100 permutations to estimate the FDR threshold. Per permutation, genotypes and ages were independently randomly distributed, keeping the among-gene structure intact. Then for each spot (23,232) on the array, model 2 was tested. The obtained P-values were used to estimate a threshold for each of the explanatory factors. Authors also used a genome-wide threshold of -log10 P-value = 2, which resembles an FDR of 0.072 and 0.060 for marker and the interaction age-marker for the developing worms and FDR of 0.050 and 0.065 for marker and age-marker for the aging worms. For the physiological age effect, authors used a log10 P-value = 8 in developing worms (0.012 FDR) and -log10 P-value = 6 (0.032 FDR). |
WBPaper00040858:eQTL_age_regulated_developing
|
|
|
Genes with expression level regulated by genotype (N2 vs CB4856) at old adults stage (214 hours at 24 centigrade). |
For model 2, authors used 100 permutations to estimate the FDR threshold. Per permutation, genotypes and ages were independently randomly distributed, keeping the among-gene structure intact. Then for each spot (23,232) on the array, model 2 was tested. The obtained P-values were used to estimate a threshold for each of the explanatory factors. Authors also used a genome-wide threshold of -log10 P-value = 2, which resembles an FDR of 0.072 and 0.060 for marker and the interaction age-marker for the developing worms and FDR of 0.050 and 0.065 for marker and age-marker for the aging worms. For the physiological age effect, authors used a log10 P-value = 8 in developing worms (0.012 FDR) and -log10 P-value = 6 (0.032 FDR). |
WBPaper00040858:eQTL_regulated_aging
|
|
|
Genes with expression level regulated by genotype (N2 vs CB4856) at Late reproduction stage (96 hours at 24 centigrade). |
Authors permuted transcript values and used a genome-wide threshold of log10 P-value = 2, which resembles a false discovery rate (FDR) of 0.0118. |
WBPaper00040858:eQTL_regulated_reproductive
|
|
|
Genes found to be regulated by low-copy overexpression of sir-2.1 with p < 0.014. |
N.A. |
WBPaper00026929:sir-2.1_overexpression_regulated
|
|
|
Embryonic class (E): genes that significantly increase in abundance at some point during embryogenesis. |
A modified Welch F statistic was used for ANOVA. For each gene, regressed error estimates were substituted for observed error estimates. The substitution is justified by the lack of consistency among the most and least variable genes at each time point. Regressed error estimates were abundance-dependent pooled error estimates that represented a median error estimate from a window of genes of similar abundance to the gene of interest. A randomization test was used to compute the probability Pg of the observed F statistic for gene g under the null hypothesis that developmental time had no effect on expression. P-values were not corrected for multiple testing. |
[cgc5767]:expression_class_E
|
|
|
Strictly embryonic class (SE): genes that are the subset of embryonic genes that are not also classified as maternal. |
A modified Welch F statistic was used for ANOVA. For each gene, regressed error estimates were substituted for observed error estimates. The substitution is justified by the lack of consistency among the most and least variable genes at each time point. Regressed error estimates were abundance-dependent pooled error estimates that represented a median error estimate from a window of genes of similar abundance to the gene of interest. A randomization test was used to compute the probability Pg of the observed F statistic for gene g under the null hypothesis that developmental time had no effect on expression. P-values were not corrected for multiple testing. |
[cgc5767]:expression_class_SE
|
|
|
Transcripts that showed significantly increased expression in him-17(me24) comparing to in N2 animals. |
DESeq2, fold change > 2, FDR < 0.05. |
WBPaper00064769:him-17(me24)_upregulated
|
|
|
Genes regulated by spr-5 (greater than 2-fold change between spr-5(by101) generations 1, 13, and 26). |
Significantly differentially regulated genes were selected by using a 2-fold ifference along with intensity values > 1000. |
WBPaper00033101:spr-5_regulated
|
|
|
Genes with differential expression under 0.5mg/l Chlorpyrifos (CPF) treatment at 16 centigrade. |
To identify the differentially expressed genes in each treatment authors used linear models per toxicant and temperature (gene expression = Toxicant (effect) + error). The lm function in R stats package was used to implement the linear models analysis with recommended default options. For threshold determination authors used a permutation approach. For each of the 23,232 permutations used authors randomly picked a transcript (array spot), which could only be picked once. Authors combined all the expression values of this transcript and randomly distributed them over the replicates and used them in the linear model. In this way authors obtained a threshold for each of the toxicants. Authors used a -log10 p-value 2 as common threshold for the analysis, which resembles to the following FDR per toxicant: 0.0155 for CPF at 24 centigrade, 0.0148 for DZN at 24 centigrade, 0.0168 for CPF+DZN at 24 centigrade, 0.0142 for CPF at 16 centigrade, 0.0151 for DZN at 16 centigrade, and 0.0148 for CPF+DZN, at 16 centigrade. |
WBPaper00040210:Chlorpyrifos_16C_regulated
|
|
|
Genes with differential expression under 0.5mg/l Chlorpyrifos (CPF) and 0.5mg/l Diazinon (DZN) treatment at 24 centigrade. |
To identify the differentially expressed genes in each treatment authors used linear models per toxicant and temperature (gene expression = Toxicant (effect) + error). The lm function in R stats package was used to implement the linear models analysis with recommended default options. For threshold determination authors used a permutation approach. For each of the 23,232 permutations used authors randomly picked a transcript (array spot), which could only be picked once. Authors combined all the expression values of this transcript and randomly distributed them over the replicates and used them in the linear model. In this way authors obtained a threshold for each of the toxicants. Authors used a -log10 p-value 2 as common threshold for the analysis, which resembles to the following FDR per toxicant: 0.0155 for CPF at 24 centigrade, 0.0148 for DZN at 24 centigrade, 0.0168 for CPF+DZN at 24 centigrade, 0.0142 for CPF at 16 centigrade, 0.0151 for DZN at 16 centigrade, and 0.0148 for CPF+DZN, at 16 centigrade. |
WBPaper00040210:Chlorpyrifos_Diazinon_24C_regulated
|
|
|
Genes with expression level regulated by genotype (N2 vs CB4856) at Old adults stage (214 hours at 24 centigrade). |
Authors permuted transcript values and used a genome-wide threshold of log10 P-value = 2, which resembles a false discovery rate (FDR) of 0.0136. |
WBPaper00040858:eQTL_regulated_old
|
|
|
Embryonic (E) subclasses are based on the earliest significant increase(abbreviated pi for primary increase). |
A modified Welch F statistic was used for ANOVA. For each gene, regressed error estimates were substituted for observed error estimates. The substitution is justified by the lack of consistency among the most and least variable genes at each time point. Regressed error estimates were abundance-dependent pooled error estimates that represented a median error estimate from a window of genes of similar abundance to the gene of interest. A randomization test was used to compute the probability Pg of the observed F statistic for gene g under the null hypothesis that developmental time had no effect on expression. P-values were not corrected for multiple testing. |
[cgc5767]:expression_class_E_pi(101_min)
|
|
|
Expression Pattern Group B, enriched for genes involved in embryonic development. These patterns have in common that they all have genes of which the expression goes up after the juvenile stage. The expression of the genes in these patterns remains high or even goes up after reproduction. |
The significance (P 0.0001) of the relative age (time) was used to determine if a gene was differentially expressed between the three age (time) groups. The effect of this factor explaining gene expression differences was used to determine if the expression went up or down during the two age/time periods (t1 - t2 and t2 -t3). Authors used a permutation approach to determine the thresholds for the different mapping strategies. For each of the used models for eQTL mapping, authors used 23,000 permutations. For each permutation, authors randomly picked a spot; each spot could only be picked once. The gene expression and relative lifespan values were than randomly distributed over the RILs (and time points) and used for mapping. In this way, authors obtained a threshold for each of the explaining factors. For the single time points, authors used a FDR of 0.01 to adjust for multiple testing. The genome-wide threshold for this FDR is -log10 P = 3.8 for each of the three time points. For the combined models (t1 to t2 and t2 to t3), authors used a genome-wide threshold of -log10 P = 4, which resembles an FDR of 0.006, 0.001, and 0.006 for marker, age, and the interaction between marker and age, respectively. To determine the threshold for the single gene examples, authors used 1000 permutations as in the genome-wide threshold. The difference is that they use the gene under study in all of the permutations. The P-values for the gene specific thresholds were determined at FDR = 0.05. |
WBPaper00036286:Pattern_B
|
|
|
mRNAs that were significantly enriched in the AIN-1 immunoprecipitation samples, compared to the control total mRNAs in the input extracts (p < 0.01). |
Signals from replicates of total worm lysates from wt and strains containing the ain-2::gfp or the ain-2 promoter::gfp transgene were mean normalized and averaged respectively to generate standard profiles of gene expression in these worm strains. Authors then calculated the ratio of signal of each gene from each IP sample to the standard gene expression profile of the corresponding worm strain. Based on this ratio, a percentile rank of each gene relative to all genes in each IP replicate was calculated. The percentile ranks in the three replicates of each IP were averaged. Student t test was utilized to determine if the average percentile ranks of enrichment of individual genes were significantly higher (p value) than the mean enrichment of all genes in the IP samples. To determine the AIN-1 or AIN-2 associated genes, we used the following criteria: (1) average percentile ranks of enrichment is greater than the mean enrichment of all genes in AIN-1 or AIN-2 IP with p < 0.01; (2) average signal in AIN-1 or AIN-2 IP replicates is greater than the background signal (including 2X standard deviation (SD)) (Background signal and SD were calculated based on signals from empty spots on each microarray); (3) criteria 1 is not be satisfied for the same gene in the corresponding control IP. |
WBPaper00031252:AIN-1_IP_enriched
|
|
|
Expression Pattern Group E, enriched for genes involved in dephosphorylation. |
The significance (P 0.0001) of the relative age (time) was used to determine if a gene was differentially expressed between the three age (time) groups. The effect of this factor explaining gene expression differences was used to determine if the expression went up or down during the two age/time periods (t1 - t2 and t2 -t3). Authors used a permutation approach to determine the thresholds for the different mapping strategies. For each of the used models for eQTL mapping, authors used 23,000 permutations. For each permutation, authors randomly picked a spot; each spot could only be picked once. The gene expression and relative lifespan values were than randomly distributed over the RILs (and time points) and used for mapping. In this way, authors obtained a threshold for each of the explaining factors. For the single time points, authors used a FDR of 0.01 to adjust for multiple testing. The genome-wide threshold for this FDR is -log10 P = 3.8 for each of the three time points. For the combined models (t1 to t2 and t2 to t3), authors used a genome-wide threshold of -log10 P = 4, which resembles an FDR of 0.006, 0.001, and 0.006 for marker, age, and the interaction between marker and age, respectively. To determine the threshold for the single gene examples, authors used 1000 permutations as in the genome-wide threshold. The difference is that they use the gene under study in all of the permutations. The P-values for the gene specific thresholds were determined at FDR = 0.05. |
WBPaper00036286:Pattern_E
|
|
Bacteria infection: Erwinia carotovora |
Genes down-regulated in animals infected with Erwinia carotovora compared to the E. coli OP50 control after 24h of infection. |
MAANOVA and BRB-Array-Tools. |
WBPaper00030985:Erwinia_carotovora_downregulated
|
|
|
mRNAs that were significantly enriched in the AIN-2 immunoprecipitation samples, compared to the control total mRNAs in the input extracts (p < 0.01). |
Signals from replicates of total worm lysates from wt and strains containing the ain-2::gfp or the ain-2 promoter::gfp transgene were mean normalized and averaged respectively to generate standard profiles of gene expression in these worm strains. Authors then calculated the ratio of signal of each gene from each IP sample to the standard gene expression profile of the corresponding worm strain. Based on this ratio, a percentile rank of each gene relative to all genes in each IP replicate was calculated. The percentile ranks in the three replicates of each IP were averaged. Student t test was utilized to determine if the average percentile ranks of enrichment of individual genes were significantly higher (p value) than the mean enrichment of all genes in the IP samples. To determine the AIN-1 or AIN-2 associated genes, we used the following criteria: (1) average percentile ranks of enrichment is greater than the mean enrichment of all genes in AIN-1 or AIN-2 IP with p < 0.01; (2) average signal in AIN-1 or AIN-2 IP replicates is greater than the background signal (including 2X standard deviation (SD)) (Background signal and SD were calculated based on signals from empty spots on each microarray); (3) criteria 1 is not be satisfied for the same gene in the corresponding control IP. |
WBPaper00031252:AIN-2_IP_enriched
|
|
|
Down-regulated genes under 0.5mg/l CPF treatment at 16 centigrade. |
The Rank Product package was used to identify the differentially expressed genes between controls and treatment in each experiment. Briefly, genes were ranked based on up- or downregulation by the treatment in each experiment. Then, for each gene a combined probability was calculated as a rank product (RP). The RP values were used to rank the genes based on how likely it was to observe them by chance at that particular position on the list of differentially expressed genes. The RP can be interpreted as a p-value. To determine significance levels, the RP method uses a permutation-based estimation procedure to transform the p-value into an e-value that addresses the multiple testing problem derived from testing many genes simultaneously. Genes with a percentage of false-positives (PFP) < 0.05 were considered differentially expressed between treatments and control in each experiment. This method has the advantage to identify genes with a response to the toxicants even when the absolute effect of the response was low. Because authors used sub-lethal concentrations of the toxicants, methods that use thresholds based on absolute fold change would not identify small changes in gene expression. Moreover, RP has proved to be a robust method for comparing microarray data from different sources and experiments. |
WBPaper00037113:CPF_16C_down-regulated
|
|
|
Proteins interacting with TEBP-1-FLAG-GFP (referred as DTN-2 by authors) during immunoprecipitation. |
N.A. |
WBPaper00060908:TEBP-1_interacting
|
|
|
Proteins interacting with TEBP-2-FLAG-GFP (referred as DTN-2 by authors) during immunoprecipitation. |
N.A. |
WBPaper00060908:TEBP-2_interacting
|
|
|
Genome-wide analysis of developmental and sex-regulated gene expression profile. |
self-organizing map |
cgc4489_group_19
|
|
|
Transcripts that showed significantly increased expression in scc-1(ubs19) comparing to in control animals. |
DESeq2 v.1.34, fold change > 2, FDR < 0.05. |
WBPaper00067078:scc-1(ubs19)_upregulated
|
|
|
Strictly embryonic (SE) subclasses are based on the earliest significant increase(abbreviated pi for primary increase). |
A modified Welch F statistic was used for ANOVA. For each gene, regressed error estimates were substituted for observed error estimates. The substitution is justified by the lack of consistency among the most and least variable genes at each time point. Regressed error estimates were abundance-dependent pooled error estimates that represented a median error estimate from a window of genes of similar abundance to the gene of interest. A randomization test was used to compute the probability Pg of the observed F statistic for gene g under the null hypothesis that developmental time had no effect on expression. P-values were not corrected for multiple testing. |
[cgc5767]:expression_class_SE_pi(101_min)
|
|
|
Developmentally modulated gene cluster. |
self-organizing map |
cgc4386_cluster_5_3
|
|
|
Early embryonic development gene expression profile. |
QT clustering |
[cgc5767]:cluster_81
|
