Comparison of expression values between year and Year2, leaf_preflowering
This document has last been compiled on 2021-12-14 23:44:12.
We ran the same analysis pipeline from preprocessing to differential gene expression (DE) analysis on data when it was mapped to BTx623 (we will call it Year2) and on data when it was mapped to specific genotype, BTx642 or RTx430, (we will call it BT642Year2). In this report, we will analyze the overlap of DE genes identified between Year2 and BT642Year2, focusing on the results from DE spline analysis.
In this file, we compare the raw fpkm counts of samples from leaf and Preflowering between Year2 and RemapYear2.
Below is the analysis of overlap from the raw fpkm counts between Year2 and BT642Year2:
Load data
Load fpkm counts of all samples
Year2:
## Reading in metadata from file: results/Year2/data/rnaMetaData.txt
## Reading in data from file: results/Year2/data/fpkm_counts.txt## Total number of all samples: 578
## Total number of genes in all samples: 34211
## Removing the 50 samples identified as not part of the main experiment, from the all samples
## After filtering all samples: 528 samples, 34211 genesBT642Year2:
## Reading in metadata from file: results/BT642Year2/data/rnaMetaData.txt
## Reading in data from file: results/BT642Year2/data/fpkm_counts.txt## Total number of all samples: 290
## Total number of genes in all samples: 35115
## Removing the 26 samples identified as not part of the main experiment, from the all samples
## After filtering all samples: 264 samples, 35115 genesFilter samples not in leaf_preflowering
Year2:
## After filtering: 98 samples, 34211 genesBT642Year2:
## After filtering: 50 samples 35115 genesMatch BTx623 to OrthologGroup
To compare between Year2 and BT642Year2, we need to match the gene naming. In Year2, genes are named with BT623; and in BT642Year2, genes are names with Ortholog Group. So we implemented mapping the BTx623 gene names to the ortholog groups according to the annotation file.
Mapping BTx623 to OrthologGroup:
## Out of 34211 genes with BT623 names,
## 28538 genes are matched to remap gene names.
## Since 0 out of 34056 remap gene names are mapped to multiple BT623 gene names,
## we collapsed on BT623 with average values of counts so that each remap gene has a unique row.## After mapping the BTx623 gene names to the remap gene names,
## Year2 leaf Preflowering : 98 samples, 34056 genes
## BT642Year2 leaf Preflowering : 50 samples, 35115 genes
## 
Compare gene expressions in Year2 and BT642Year2
In this section, we will compare the log fpkm counts of genes in Year2 and BT642Year2. Particularly, since there are many samples for each gene, we will use the 3rd max log fpkm counts, the average of log fpkm counts in Day063, and the log log fpkm count from a particular sample of each gene to compare.
There are genes that their counts are results from collapsing and we provide comparison of Year2 and BT642Year2 of these genes. Some genes do not have one-to-one relationship of RT430-to-BT642, so their counts were collapsed by summation when we were merging data from two genotypes. They are flaged as Collasped On Ortholog. Some genes do not have one-to-one relationship of BT623-to-OrthologGroup. so their counts were collapsed by averaging when we were mapping BT623 gene names to ortholog groups. They are flaged as Collapsed on BT623. We also provide analysis that focusing on these genes with collapsed values.
3rd max log expression values
Compare the 3rd max log fpkm counts of genes in Year2 and BT642Year2
Scatter plots
Here are some genes that are only expressed in one year.
10 genes that were expressed in BT642Year2 but not in Year2 (sort by expression vales in BT642Year2):
## GeneID collapsed value_remap value_623 value_diff
## 15990 SbiBTX642.04G233300 not1-1 8.764308 0 8.764308
## 15991 SbiBTX642.04G233300 not1-1 8.764308 0 8.764308
## 14751 SbiBTX642.04G091700 1-1 6.987548 0 6.987548
## 13918 SbiBTX642.04G003700 1-1 6.985273 0 6.985273
## 23732 SbiBTX642.07G071800 1-1 5.969703 0 5.969703
## 7769 SbiBTX642.02G264200 1-1 4.623516 0 4.623516
## 3210 SbiBTX642.01G354500 1-1 4.579542 0 4.579542
## 7603 SbiBTX642.02G247100 1-1 4.150560 0 4.150560
## 19408 SbiBTX642.05G181900 not1-1 4.141596 0 4.141596
## 28764 SbiBTX642.09G172500 1-1 3.880686 0 3.88068610 genes that were expressed in Year2 but not in BT642Year2 (sort by expression vales in Year2):
## 10 genes that were expressed in BT642Year2 but not in Year2 :## GeneID collapsed value_remap value_623 value_diff
## 21025 SbiBTX642.06G103500 1-1 0 10.986745 -10.986745
## 1943 SbiBTX642.01G204000 1-1 0 6.240123 -6.240123
## 8174 SbiBTX642.02G304600 not1-1 0 5.598425 -5.598425
## 13480 SbiBTX642.03G436300 1-1 0 5.538538 -5.538538
## 23056 SbiBTX642.07G007500 1-1 0 5.183487 -5.183487
## 10553 SbiBTX642.03G104200 1-1 0 5.076388 -5.076388
## 11939 SbiBTX642.03G266300 1-1 0 5.050502 -5.050502
## 31841 SbiBTX642.10G231500 1-1 0 4.612942 -4.612942
## 7592 SbiBTX642.02G245800 1-1 0 4.539779 -4.539779
## 7179 SbiBTX642.02G199300 1-1 0 4.498889 -4.498889Density Plots
Histogram of differences
Mean log expression values in Day63
Compare the average of log fpkm counts of genes in DAY63 in Year2 and BT642Year2
Scatter plots
Here are some genes that are only expressed in one year.
10 genes that were expressed in BT642Year2 but not in Year2 (sort by expression vales in BT642Year2):
## GeneID collapsed value_remap value_623 value_diff
## 15990 SbiBTX642.04G233300 not1-1 8.058909 0 8.058909
## 15991 SbiBTX642.04G233300 not1-1 8.058909 0 8.058909
## 13918 SbiBTX642.04G003700 1-1 6.606435 0 6.606435
## 14751 SbiBTX642.04G091700 1-1 6.355378 0 6.355378
## 23732 SbiBTX642.07G071800 1-1 5.034489 0 5.034489
## 7769 SbiBTX642.02G264200 1-1 4.166679 0 4.166679
## 7603 SbiBTX642.02G247100 1-1 3.884710 0 3.884710
## 3210 SbiBTX642.01G354500 1-1 3.767260 0 3.767260
## 19957 SbiBTX642.05G227500 not1-1 3.750745 0 3.750745
## 19958 SbiBTX642.05G227500 not1-1 3.750745 0 3.75074510 genes that were expressed in Year2 but not in BT642Year2 (sort by expression vales in Year2):
## 10 genes that were expressed in BT642Year2 but not in Year2 :## GeneID collapsed value_remap value_623 value_diff
## 21025 SbiBTX642.06G103500 1-1 0 10.262791 -10.262791
## 5447 SbiBTX642.02G010300 1-1 0 5.903575 -5.903575
## 11939 SbiBTX642.03G266300 1-1 0 4.855084 -4.855084
## 1943 SbiBTX642.01G204000 1-1 0 4.365491 -4.365491
## 13480 SbiBTX642.03G436300 1-1 0 4.218546 -4.218546
## 23056 SbiBTX642.07G007500 1-1 0 3.970972 -3.970972
## 2383 SbiBTX642.01G252900 1-1 0 3.852990 -3.852990
## 11725 SbiBTX642.03G242600 not1-1 0 3.812745 -3.812745
## 11726 SbiBTX642.03G242600 not1-1 0 3.812745 -3.812745
## 8174 SbiBTX642.02G304600 not1-1 0 3.808631 -3.808631Density Plots
Histogram of differences
Log expression values of a single sample
Compare the log fpkm counts of genes from a single sample in Year2 and BT642Year2
## The randomly selected sample: 0816174L14 from BT642:Day063:Preflowering:LeavesScatter plots
