|
|
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
|
|
|
Maternal class (M): genes that are called present in at least one of the three PC6 replicates. |
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_M
|
|
|
Transcripts that showed significantly decreased expression in daf-16(mgDf50) comparing to in N2 at L1 larva stage. |
DESeq v1.20.0 was used to analyze differential gene expression. Transcripts with adjusted p-value < 0.05 were considered differentialled expressed. |
WBPaper00048971:daf-16(mgDf50)_downregulated_L1
|
|
|
Genes down regulated by mir-243(n4759). |
RNAs that changed at least 2-fold with a probability of p > 0.05 in three biological replicates were considered differentially regulated between wild-type and mir-243. |
WBPaper00036130:mir-243_down_regulated
|
|
|
Genes expressed in N2. |
Expressed transcripts were identified on the basis of a Present call in 3 out of 4 N2 experiments as determined by Affymetrix MAS 5.0. |
WBPaper00025141:N2_Expressed_Genes
|
|
|
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
|
|
Bacteria diet: Lactobacillus rhamnosus CNCM I-3690 |
Differentially expressed genes between C.elegans fed with L.rhamnosus strain CNCM I-3690 and control strain E.coli OP50 after 3 days of feeding. |
Differential expression was assayed via limma. Genes were considered differentially expressed when the multiple testing adjusted P-value < 0.01. |
WBPaper00041906:I-3690_vs_OP50_3d
|
|
|
Genes that showed significant differential expressed between control and 20 mg\/L Cadmium treatment. |
t-test, p < 0.05. |
WBPaper00036123:Cadmium_regulated
|
|
|
Genes expressed in embryonic motor neurons (identified by unc-4::GFP expressing cells). |
Genes called Present by MAS 5.0 in 2 out of 3 unc-4::GFP hybridizations. |
WBPaper00025141:unc-4::GFP_Expressed_Genes
|
|
|
Transcripts with significantly increased expression after treatment with 0.1mM paraquat vs. control |
Comparisons of each genotype were compared to the wild-type using the Empirical Base (Wright & Simon) algorithm and fold changes were represented on a log2 scale. A threshold of p < 0.05 and a fold change of 1.3 (log2) was set to determine differentially expressed targets. |
WBPaper00045263:0.1mM-paraquat_upregulated
|
|
|
Embryonic transient class (ET): genes that are the subset of embryonic genes in which the latest significant increase is earlier than their latest significant decrease. |
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_ET
|
|
|
Maternal-embryonic class (ME): genes that are in the intersection of the maternal and embryonic classes. |
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_ME
|
|
|
Maternal-embryonic transient class (MET): genes that are in the intersection of the maternal and embryonic transient classes. |
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_MET
|
|
|
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
|
|
|
Genes from eat-2(ad465) animals with significantly increased expression after 72 hours of treatment on growth media with 10uM rapamycin in 2% DMSO. |
Analysis of gene expression data was carried out with the Affymetrix Transcriptome Analysis Console. Data preprocessing (using RMA normalization) and QC metrics were performed using Affymetrix Expression Console TM and manually inspected afterwards. Expression analysis was carried out for each two pairwise conditions. FDR statistical correction for multiple testing resulted in a slightly lower number of DEGs in most cases. P-value < 0.05 and fold change > 2.0 were used to determine differentially expressed genes. |
WBPaper00048989:eat-2(ad465)_rapamycin_upregulated
|
|
|
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
|
|
|
Transcripts that showed significantly decreased expression in ogt-1(ok1474) neuronal cells isolated by FACs comparing to in FACs isolated neuronal cells from wild type. |
DESeq2, fold change > 2, FDR < 0.05. |
WBPaper00066485:ogt-1(ok1474)_downregulated_neuron
|
|
|
Transcripts with significantly decreased expression in nuo-6(qm200) vs. N2, and in nuo-6(qm200);ced-4(n1162) vs. ced-4(n1162). |
Comparisons of each genotype were compared to the wild-type using the Empirical Base (Wright & Simon) algorithm and fold changes were represented on a log2 scale. A threshold of p < 0.05 and a fold change of 1.3 (log2) was set to determine differentially expressed targets. |
WBPaper00045263:nuo-6(qm200)_downregulated
|
|
|
Genes with changed expression in lin-54(n2290) embryo. |
To determine differentially expressed genes between wild-type and mutants, moderated T Statistics in limma was used with p-value <= 0.01, fold change >= 1.5. |
WBPaper00038427:lin-54_embryo_changed
|
|
|
Expression Pattern Group H, 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_H
|
|
|
Embryonic transient (ET) subclasses are based on time of max abundance. |
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_ET_max(66_min)
|
|
|
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(23_min)
|
|
Treatment with 2.0mM of HuminFeed until young adult stage (3 days). |
Gene significantly down-regulated by treatment with 2.0mM of HuminFeed until young adult stage (3 days), with a minimum fold change in gene expression of 0.8. |
For selection of DEGs, an unpaired t -test was performed followed by a significance analysis of microarray (SAM) test including a calculation that estimates the false discovery rate (FDR). FDR, reducing on the one hand type I errors for null associations, was set to a non-stringent level of <12.5%, mainly to guard from an increase of type II error and also based on findings by Levine et al. (2011), which described 12.5% as most acceptable optimum level of FDR, representing the 90th percentile of the normal distribution curve. DEGs exceeding a fold change of 1.25 were further analyzed with respect to their functional clustering. This fold-cut-off was chosen to allow an interpretation that is biologically meaningful, akin to the notion that data of sound technical and experimental quality which returns strong, statistically significant, absolute signal intensities is sufficiently robust to justify a fold-cut-off of >1.2. This analysis was conducted using the functional annotation clustering tool of the Database for Annotation, Visualization, and Integrated Discovery (DAVID; Huang et al., 2007). |
WBPaper00041002:HF_3d_2.0mM_Down
|
|
|
Genes that showed significantly increased expression after exposure to 1mg/L MWCNTs from L1 larva to young adult. |
Transcripts with false discovery rate-corrected p-values < 0.05 and fold change > 2 were defined as differentially expressed. |
WBPaper00049377:MWCNT_upregulated
|
|
|
FBF-associated probe sets (FDR <2.25%) |
SAM assigns a score to each probe set and estimates their false discovery rates (FDRs). SAM deemed 4722 probe sets as significantly enriched at an FDR of 2.25% or lower. |
WBPaper00035905:FBF-1_Associated
|
|
Bacteria infection: Serratia marcescens |
Caenorhabditis elegans Genes with expression levels changed significantly after treatment of Serratia marcescens. |
Differential expression were calculated by empirical eBayes method using eBayes function. P_value <= 0.01 and log2 fold change > 1 were used to call differentially expressed genes in all datasets. |
WBPaper00041606:CE_S.marcescens_regulated
|
|
|
Early embryonic development gene expression profile. |
QT clustering |
[cgc5767]:cluster_10
|
|
|
C-lineage related expression profile. |
QT clustering |
WBPaper00025032:cluster_14
|
|
|
Genes with differential expression in prg-1(tm872) male gonads comparing to N2. |
The fold differences from the three experiments were averaged, and a Z test [Z = (observed - expected) / SE] was performed. Genes with up- or downregulation greater than 1.5-fold, p < 0.05, in any given mutant were selected for inclusion. |
WBPaper00031903:prg-1_regulated
|
