|
|
mRNAs that showed decreased expression in 1 cell mebryo comparing to in oocyte, according to RNAseq analysis. |
Gaussian error propagation. As cutoff for the up-regulated genes authors used log2 fold change > 1 and P < 0.05 and as cutoff for the down-regulated genes authors used log2 fold change < -1 and P < 0.05. |
WBPaper00045420:fertilization_downregulated_transcript
|
|
Bacteria infection: Enterococcus faecalis |
Genes with increased expression after 24 hours of infection by E.faecalis Fold changes shown are pathogen vs OP50. |
For RNA-seq and tiling arrays, log2 fold changes between gene expression values of infected versus uninfected nematodes were calculated. For log2 fold changes > 0.00001 the values > 81.25th percentile were defined as up-regulated and for log2 fold changes < -0.00001 the values < 18.75th percentile were defined as down-regulated. |
WBPaper00038438:E.faecalis_24hr_upregulated_TilingArray
|
|
|
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
|
|
|
Genes that showed expression levels higher than the corresponding reference sample (L2 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:A-class-motor-neurons_L2-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 with significantly increased expression in eat-2(ad465) treated with 2% DMSO for 72 hours, comparing to in N2 treated with 2% DMSO for 72 hours. |
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)_upregulated_in-DMSO
|
|
|
Genes significantly enriched in NSM neurons (isolated by FACS) versus the reference, according to tiling array analysis towards total RNA. |
A linear model and moderated t-statistic were used to determine differentially expressed genes as implemented by the limma package (v3.21.4). Enriched list contains only genes significantly enriched in the NSM neurons versus the reference <=1.5X and <= 5% FDR. |
WBPaper00045974:NSM_enriched_totalRNA_tiling
|
|
|
Genes that showed expression levels higher than the corresponding reference sample (Young adult 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:CEP-sheath-cells_Day1-adult_expressed
|
|
|
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 found to be regulated in daf-16(mgDf50) by resveratrol treatment with p < 0.01. |
N.A. |
WBPaper00026929:Resveratrol_regulated_daf-16
|
|
Bacteria infection: Photorhabdus luminescens |
Genes with increased expression after 24 hours of infection by P.lumniescens Fold changes shown are pathogen vs OP50. |
For RNA-seq and tiling arrays, log2 fold changes between gene expression values of infected versus uninfected nematodes were calculated. For log2 fold changes > 0.00001 the values > 81.25th percentile were defined as up-regulated and for log2 fold changes < -0.00001 the values < 18.75th percentile were defined as down-regulated. |
WBPaper00038438:P.lumniescens_24hr_upregulated_TilingArray
|
|
|
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
|
|
|
miRNA targets that are significantly enriched at embryo stage. To generate a global view of the dynamics of miRNA-mediated regulation of gene expression during C. elegans development, authors analyzed the mRNAs in the AIN-2-GFP IP results from five developmental stages. The magnitude of the combined interaction of miRNAs with a given target mRNA was assessed by measuring the fold enrichment of that mRNA in AIN-2 IP samples, relative to the abundance of the mRNA in the corresponding total lysate. Because this enrichment in the IP sample versus total lysate directly reflects the miRISC-associated fraction of a given mRNA, high enrichment indicates the likelihood of strong miRNA-mediated regulation of the mRNA, whereas low or negative enrichment indicates the likelihood of weak or absent miRNA regulation of the mRNA. It is also possible that poor enrichment could reflect interactions that occur only in a rare subset of cells at any given stage of development. Transcripts that were significantly enriched (0
|
For each transcript represented on the microarray, authors first used all 18 data points from all of the 18 microarrays to calculate its average in-stage standard deviation as SD=SQRT((i=1-K (ni-1)SDi2)/(N-K)), degree of freedom (df)=N-K, where K is the number of developmental stages in which the given transcript had at least one non-missing log2(IP/Total) value; ni is the number of non-missing values in stage i among the K stages; SDi is the standard deviation of the log2(IP/Total) values from all the replicates of stage i; N is the total number of non-missing values for this transcript at all stages. A transcript must have at least one stage with at least two non-missing values to be testable. All of the non-missing values of a transcript at each stage were averaged to generate the stage-average values (Mi for stage i). The standard error of Mi was calculated as SEi=SD/SQRT(ni) (df=N-K), where SD is the average in-stage standard deviation calculated above and ni is the number of non-missing values at stage i, as explained above. Based on the SEi, a one-tailed Students t-test was used to calculate the P-value of enrichment in stage i (Test if Mi>0). The T statistic was constructed as T=Mi/SEi (df=N-K). Mi>0 and enrichment P<0.001 were used as the threshold of enrichment for each stage. |
WBPaper00035084:embryo_enriched_AIN-2_IP
|
|
|
|
Expression Pattern Group F, 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_F
|
|
heat shock |
Genes for which heat shock F3 (fraction 3, containing heavy polysomes) versus control F3 is significantly increased. |
Survival during heat stress was analyzed using a non-parametric (Mantel-Haenszel) Log rank test and presented as Kaplan-Meier survival curves (Prism software package). For the TSA analysis all full genome wide expression analyses involved comparisons of two groups and thus on a gene by gene basis, for each set, authors performed two-sample t-tests on the log2 transformed expression measures to derive raw p-values. The Benjamini and Hochberg was then used to report the inference adjusted for multiple comparisons (q-values), in their case controlling false discovery rate (FDR). Statistically significant differential expression was defined as FDR q-value <= 0.05, thereby controlling FDR at <= 0.05 for each set of comparisons. Shown here are heat shock vs control ratio and raw p_value based on two-sample t-test. |
WBPaper00037147:heatshock_upregulated
|
|
|
Genes that showed significantly increased 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_upregulated
|
