|
miRNA with decreased expression in N2 L3 larva comparing to in N2 L2 larva. |
DEGseq and GFOLD were used to analyze miRNAs expression. Authors chose miRNAs which had more than two-fold difference in expression (P < 0.001, q < 0.01 of Storey) from DEGseq, and miRNAs which had more than two-fold difference in expression (GFOLD score > 0 for up-regulation and GFOLD score < 0 for down-regulation) from GFOLD outcomes. Then authors obtained the intersection of up-regulated miRNAs and down-regulated miRNAs for each stage from the chosen miRNAs, respectively. |
WBPaper00045316:miRNA_N2_L3_vs_L2_downregulated_adult
|
|
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 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_regulated_developing
|
|
Transcripts that showed significantly increased expression in rrf-3(pk1426) comparing to in N2 at L3 larva stage. |
DESeq2v 1.18.1, fold change > 1.5, adjusted p-value < 0.01. |
WBPaper00056169:rrf-3(pk1426)_upregulated_L3
|
|
Expression Pattern Group C, enriched for genes involved in metabolic processes. |
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_C
|
|
Genes that showed oscillating mRNA expression level throughout the 16 hour time courses from L3 larva to young adult. |
The following three lines of R code were used to perform the classification: increasing <-2*amplitude-PC1 < -1.7; oscillating <-!increasing & (amplitude > 0.55); flat <-!increasing & !oscillating; Note that the amplitude of a sinusoidal wave corresponds to only half the fold change between trough and peak. |
WBPaper00044736:oscillating_dev_expression
|
|
Genes that showed expression levels higher than the corresponding reference sample (L3/L4 all cell reference). |
A Mann-Whitney U test with an empirical background model and FDR correction for multiple testing was used to detect expressed transcripts (Benjamini and Hochberg 1995). Genes and TARs with an FDR <= 0.05 were reported as expressed above background. Authors detected differentially expressed transcripts using a method based on linear models. Genes and TARs were called differentially expressed if the FDR was <= 0.05 and the fold change (FC) >= 2.0. To more strictly correct for potential false-positives resulting from multiple sample comparisons, authors divided individual FDR estimates by the number of samplesor sample comparisons, respectively. This resulted in an adjusted FDR of 1.3 * 0.0001 for expression above background and of 7.4 * 0.0001 for differential expression. Authors called genes selectively enriched in a given tissue if they met the following requirements: (1) enriched expression in a given tissue (FDR <= 0.05 and FC >= 2.0), (2) fold change versus reference among the upper 40% of the positive FC range observed for this gene across all tissues, and (3) fold-change entropy among the lower 40% of the distribution observed for all genes. |
WBPaper00037950:dopaminergic-neurons_L3-L4-larva_expressed
|
|
Genes that showed expression levels higher than the corresponding reference sample (L3/L4 all cell reference). |
A Mann-Whitney U test with an empirical background model and FDR correction for multiple testing was used to detect expressed transcripts (Benjamini and Hochberg 1995). Genes and TARs with an FDR <= 0.05 were reported as expressed above background. Authors detected differentially expressed transcripts using a method based on linear models. Genes and TARs were called differentially expressed if the FDR was <= 0.05 and the fold change (FC) >= 2.0. To more strictly correct for potential false-positives resulting from multiple sample comparisons, authors divided individual FDR estimates by the number of samplesor sample comparisons, respectively. This resulted in an adjusted FDR of 1.3 * 0.0001 for expression above background and of 7.4 * 0.0001 for differential expression. Authors called genes selectively enriched in a given tissue if they met the following requirements: (1) enriched expression in a given tissue (FDR <= 0.05 and FC >= 2.0), (2) fold change versus reference among the upper 40% of the positive FC range observed for this gene across all tissues, and (3) fold-change entropy among the lower 40% of the distribution observed for all genes. |
WBPaper00037950:hypodermis_L3-L4-larva_expressed
|
|
Genes that showed expression levels higher than the corresponding reference sample (L3/L4 all cell reference). |
A Mann-Whitney U test with an empirical background model and FDR correction for multiple testing was used to detect expressed transcripts (Benjamini and Hochberg 1995). Genes and TARs with an FDR <= 0.05 were reported as expressed above background. Authors detected differentially expressed transcripts using a method based on linear models. Genes and TARs were called differentially expressed if the FDR was <= 0.05 and the fold change (FC) >= 2.0. To more strictly correct for potential false-positives resulting from multiple sample comparisons, authors divided individual FDR estimates by the number of samplesor sample comparisons, respectively. This resulted in an adjusted FDR of 1.3 * 0.0001 for expression above background and of 7.4 * 0.0001 for differential expression. Authors called genes selectively enriched in a given tissue if they met the following requirements: (1) enriched expression in a given tissue (FDR <= 0.05 and FC >= 2.0), (2) fold change versus reference among the upper 40% of the positive FC range observed for this gene across all tissues, and (3) fold-change entropy among the lower 40% of the distribution observed for all genes. |
WBPaper00037950:PVD-OLL-neurons_L3-L4-larva_expressed
|
|
Transcripts that showed significantly increased expression in dpy-21(e428) comparing to in N2 during L3 stage. |
DESeq v1.6.3. Fold change > 1.5. |
WBPaper00050370:dpy-21(e428)_L3_upregulated
|
|
Transcripts that showed significantly increased expression in blmp-1(tm548) comparing to in N2. |
DESeq2, fold change > 2, p-value < 0.05. |
WBPaper00062347:blmp-1(tm548)_upregulated
|
|
Genes that showed decreased expression in nhr-8(hd117) comparing to N2. |
Expression data were read in the statistical programming environment R and processed for background correction and quantile normalization and summarized for genes that differed in expression >= 1.5-fold (N2 versus nhr-8) and a p value <= 0.01 using the Limma package. |
WBPaper00044030:nhr-8(hd117)_downregulated
|
|
Transcripts that showed significantly decreased expression in animals carrying hsp-16.2p-LIN-29 (induced at early L3 larva stage), comparing to animals carrying empty vectors. |
DESeq2 v.1.12.3, FDR < 0.05, fold change > 1.7. |
WBPaper00059191:hsp-16.2p-LIN-29_downregulated
|
|
Transcripts that showed significantly decreased expression in rrf-3(pk1426) comparing to in N2 at L3 larva stage. |
DESeq2v 1.18.1, fold change > 1.5, adjusted p-value < 0.01. |
WBPaper00056169:rrf-3(pk1426)_downregulated_L3
|
|
Transcripts that showed significantly increased expression in ints-5(RNAi) comparing to in N2. |
DESeq2 and edgeR. |
WBPaper00056284:ints-5(RNAi)_upregulated
|
|
Transcripts that showed significantly increased expression in lin-22(ot269) comparing to in N2 at L3 larva. |
Differences in gene expression were then calculated using the negative binomial test in the DESeq package (FDR = 0.1). |
WBPaper00053295:lin-22(ot269)_upregulated
|
|
Expression Pattern Group D, enriched for genes involved in catabolic processes. |
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_D
|
Exposure to Au-NP. |
Genes with differeiential expression after exposed to Au-NP. |
Hierarchical clustering was performed in Partek to confirm that the samples match to the treatment groups. Analysis of variance (ANOVA) was used to partition the variance due to treatment from technical and biological noise. The list of differentially expressed genes was generated by identifying the genes showing fold change of more than 1.5 and less than -1.5 at p < 0.05 with and without multiple sample correction, False Discovery Rate (FDR). False discovery correction according to Benjamini and Hochberg produced a list of 37 significant transcripts. FDR was not applied when selecting differential expressed genes, because this approach can increase the type II error and result in elimination of the genes responsive to the treatment. |
WBPaper00040821:Au-NP_regulated
|
|
Genes with expression level regulated by genotype (N2 vs CB4856) at L3 larva stage |
Authors permuted transcript values and used a genome-wide threshold of log10 P-value = 2, which resembles a false discovery rate (FDR) of 0.0129. |
WBPaper00040858:eQTL_regulated_juvenile
|
|
Genes that showed significantly decreased experssion after 22.5 hours of treatment in 200nM delta7-dafachronic acid comparing with in ethanol vehicle control. |
To identify the differentially expressed genes, we applied Significance Analysis of Microarrays (SAM) analysis using the R package samr [46]. Genes with median false discovery <5% and fold changes >2.0 were considered differentially expressed. |
WBPaper00046548:dafachronic-acid_downregulated
|
|
Transcripts that showed significantly increased expression in ints-1(RNAi) animals comparing to control RNAi animals, at late L3 larva stage. |
DESeq2, v1.12.0, fold change > 2, FDR < 0.05 |
WBPaper00057068:ints-4(RNAi)_upregulated
|
|
Transcripts that showed significantly increased expression in lin-22(icb38) comparing to in N2 at L3 larva. |
Differences in gene expression were then calculated using the negative binomial test in the DESeq package (FDR = 0.1). |
WBPaper00053295:lin-22(icb38)_upregulated
|
|
Transcripts that showed significantly increased expression in lin-22(ot267) comparing to in N2 at L3 larva. |
Differences in gene expression were then calculated using the negative binomial test in the DESeq package (FDR = 0.1). |
WBPaper00053295:lin-22(ot267)_upregulated
|
|
Genes with expression 1.5X depleted in PVD and OLL neurons. |
Data sets were normalized by RMA and transcripts showing relative PVD enrichment (>=1.5X) vs. the reference sample were identified by SAM analysis (False Discovery Rate, FDR < 1%) |
WBPaper00036375:depleted_in_PVD_OLL
|
|
Genes enriched in AMsh glia. |
Two different statistical methods were used for differential gene expression analysis, DESeq and voom. For DEseq analysis, DESeq2 was applied to normalize count matrix and to perform differential gene expression on the counts using negative binomial distribution; for voom analysis, edgeR was applied to normalize count matrix, and voom was applied for gene differentiation analysis. Significant genes from both analyses were combined. To identify transcripts enriched in AMsh glia compared to other cells (control AMsh versus control non-AMsh), authors used a fold change of 3.5 and an adjusted p value threshold of <0.05. |
WBPaper00049489:AMsh-glia_enriched
|
|
Genes eniched in hermaphrodites by comparing tra-2(ar221ts);xol-1(y9) XX phenotype males (at 25C) with either N2 or xol-1(y9) XX hermaphrodites during L3 larva stage |
Student t-test. An mRNA was considered to have sex-enriched expression if it showed at least a two-fold enrichment within a 99% confidence interval (P < 0.01). |
WBPaper00026629:L3_male_enriched
|
|
Significantly up-regulated genes (>= 1.5-fold change, q value < 0.05) in blmp-1(tm548) mutant versus N2 L3 larvae as obtained from RNA sequencing analysis (n = 3). |
Differential gene expression was performed with cuffdiff - which is part of the cufflinks package (v2.0.2) Genes with a FPKM value >0.5 in at least one of the N2 and blmp1 samples were considered expressed. Those with a q-value < 0.05 and fold change > +/- 1.5 were classified as differentially expressed. |
WBPaper00045017:blmp-1(tm548)_upregulated
|
|
Transcripts that showed significantly increased expression in sir-2.1(ok434) comparing to in N2 at L3 larva stage. |
DESeq2, fold change > 2, FDR < 0.05. |
WBPaper00062335:sir-2.1(ok434)_upregulated_L3
|
|
Genes that showed flat mRNA expression level throughout the 16 hour time courses from L3 larva to young adult. |
The following three lines of R code were used to perform the classification: increasing <-2*amplitude-PC1 < -1.7; oscillating <-!increasing & (amplitude > 0.55); flat <-!increasing & !oscillating; Note that the amplitude of a sinusoidal wave corresponds to only half the fold change between trough and peak. |
WBPaper00044736:flat_dev_expression
|
|
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
|