Comparison of expression values between year and Year1, leaf_preflowering
This document has last been compiled on 2021-12-14 21:57:34.
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 Year1) and on data when it was mapped to specific genotype, BTx642 or RTx430, (we will call it BT642Year1). In this report, we will analyze the overlap of DE genes identified between Year1 and BT642Year1, focusing on the results from DE spline analysis.
In this file, we compare the raw fpkm counts of samples from leaf and Preflowering between Year1 and RemapYear1.
Below is the analysis of overlap from the raw fpkm counts between Year1 and BT642Year1:
Load data
Load fpkm counts of all samples
Year1:
## Reading in metadata from file: results/Year1/data/rnaMetaData.txt
## Reading in data from file: results/Year1/data/fpkm_counts.txt## Total number of all samples: 396
## Total number of genes in all samples: 34211
## There are no samples to remove beyond those used in the main experiment
## After filtering all samples: 396 samples, 34211 genesBT642Year1:
## Reading in metadata from file: results/BT642Year1/data/rnaMetaData.txt
## Reading in data from file: results/BT642Year1/data/fpkm_counts.txt## Total number of all samples: 199
## Total number of genes in all samples: 35115
## There are no samples to remove beyond those used in the main experiment
## After filtering all samples: 199 samples, 35115 genesFilter samples not in leaf_preflowering
Year1:
## After filtering: 75 samples, 34211 genesBT642Year1:
## After filtering: 39 samples 35115 genesMatch BTx623 to OrthologGroup
To compare between Year1 and BT642Year1, we need to match the gene naming. In Year1, genes are named with BT623; and in BT642Year1, 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,
## Year1 leaf Preflowering : 75 samples, 34056 genes
## BT642Year1 leaf Preflowering : 39 samples, 35115 genes
## 
Compare gene expressions in Year1 and BT642Year1
In this section, we will compare the log fpkm counts of genes in Year1 and BT642Year1. 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 Year1 and BT642Year1 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 Year1 and BT642Year1
Scatter plots
Here are some genes that are only expressed in one year.
10 genes that were expressed in BT642Year1 but not in Year1 (sort by expression vales in BT642Year1):
## GeneID collapsed value_remap value_623 value_diff
## 15990 SbiBTX642.04G233300 not1-1 9.132268 0 9.132268
## 15991 SbiBTX642.04G233300 not1-1 9.132268 0 9.132268
## 13918 SbiBTX642.04G003700 1-1 6.617504 0 6.617504
## 14751 SbiBTX642.04G091700 1-1 6.604368 0 6.604368
## 23732 SbiBTX642.07G071800 1-1 6.417009 0 6.417009
## 28764 SbiBTX642.09G172500 1-1 5.068241 0 5.068241
## 7769 SbiBTX642.02G264200 1-1 4.195348 0 4.195348
## 18732 SbiBTX642.05G111600 1-1 4.189034 0 4.189034
## 16742 SbiBTX642.04G316200 1-1 4.060912 0 4.060912
## 16017 SbiBTX642.04G236000 not1-1 4.059182 0 4.05918210 genes that were expressed in Year1 but not in BT642Year1 (sort by expression vales in Year1):
## 10 genes that were expressed in BT642Year1 but not in Year1 :## GeneID collapsed value_remap value_623 value_diff
## 21025 SbiBTX642.06G103500 1-1 0 10.314810 -10.314810
## 5447 SbiBTX642.02G010300 1-1 0 5.971773 -5.971773
## 13480 SbiBTX642.03G436300 1-1 0 5.910733 -5.910733
## 1943 SbiBTX642.01G204000 1-1 0 5.715344 -5.715344
## 23056 SbiBTX642.07G007500 1-1 0 5.133810 -5.133810
## 11939 SbiBTX642.03G266300 1-1 0 5.092969 -5.092969
## 8174 SbiBTX642.02G304600 not1-1 0 5.045705 -5.045705
## 7179 SbiBTX642.02G199300 1-1 0 4.920293 -4.920293
## 2383 SbiBTX642.01G252900 1-1 0 4.368768 -4.368768
## 26365 SbiBTX642.08G122900 1-1 0 4.275007 -4.275007Density Plots
Histogram of differences
Mean log expression values in Day63
Compare the average of log fpkm counts of genes in DAY63 in Year1 and BT642Year1
Scatter plots
Here are some genes that are only expressed in one year.
10 genes that were expressed in BT642Year1 but not in Year1 (sort by expression vales in BT642Year1):
## GeneID collapsed value_remap value_623 value_diff
## 15990 SbiBTX642.04G233300 not1-1 7.612653 0 7.612653
## 15991 SbiBTX642.04G233300 not1-1 7.612653 0 7.612653
## 14751 SbiBTX642.04G091700 1-1 6.615860 0 6.615860
## 13918 SbiBTX642.04G003700 1-1 6.153750 0 6.153750
## 19957 SbiBTX642.05G227500 not1-1 4.413482 0 4.413482
## 19958 SbiBTX642.05G227500 not1-1 4.413482 0 4.413482
## 23732 SbiBTX642.07G071800 1-1 4.374597 0 4.374597
## 28764 SbiBTX642.09G172500 1-1 4.194962 0 4.194962
## 7769 SbiBTX642.02G264200 1-1 3.622546 0 3.622546
## 19408 SbiBTX642.05G181900 not1-1 3.500707 0 3.50070710 genes that were expressed in Year1 but not in BT642Year1 (sort by expression vales in Year1):
## 10 genes that were expressed in BT642Year1 but not in Year1 :## GeneID collapsed value_remap value_623 value_diff
## 21025 SbiBTX642.06G103500 1-1 0 9.321446 -9.321446
## 5447 SbiBTX642.02G010300 1-1 0 5.435006 -5.435006
## 13480 SbiBTX642.03G436300 1-1 0 5.220216 -5.220216
## 11939 SbiBTX642.03G266300 1-1 0 4.837798 -4.837798
## 1943 SbiBTX642.01G204000 1-1 0 4.474952 -4.474952
## 23056 SbiBTX642.07G007500 1-1 0 4.448085 -4.448085
## 136 SbiBTX642.01G014300 not1-1 0 3.692819 -3.692819
## 8174 SbiBTX642.02G304600 not1-1 0 3.646247 -3.646247
## 10513 SbiBTX642.03G099500 1-1 0 3.613208 -3.613208
## 12251 SbiBTX642.03G300400 not1-1 0 3.605711 -3.605711Density Plots
Histogram of differences
Log expression values of a single sample
Compare the log fpkm counts of genes from a single sample in Year1 and BT642Year1
## The randomly selected sample: 0817163L17 from BT642:Day077:Preflowering:LeavesScatter plots
Here are some genes that are only expressed in one year.
10 genes that were expressed in BT642Year1 but not in Year1 (sort by expression vales in BT642Year1):
## GeneID collapsed value_remap value_623 value_diff
## 15990 SbiBTX642.04G233300 not1-1 7.577957 0 7.577957
## 15991 SbiBTX642.04G233300 not1-1 7.577957 0 7.577957
## 14751 SbiBTX642.04G091700 1-1 6.125155 0 6.125155
## 13918 SbiBTX642.04G003700 1-1 5.701272 0 5.701272
## 19894 SbiBTX642.05G220700 1-1 5.330558 0 5.330558
## 19957 SbiBTX642.05G227500 not1-1 4.572890 0 4.572890
## 19958 SbiBTX642.05G227500 not1-1 4.572890 0 4.572890
## 28764 SbiBTX642.09G172500 1-1 4.430954 0 4.430954
## 23732 SbiBTX642.07G071800 1-1 4.386811 0 4.386811
## 2592 SbiBTX642.01G288400 1-1 3.814550 0 3.81455010 genes that were expressed in Year1 but not in BT642Year1 (sort by expression vales in Year1):
## 10 genes that were expressed in BT642Year1 but not in Year1 :## GeneID collapsed value_remap value_623 value_diff
## 21025 SbiBTX642.06G103500 1-1 0 9.271696 -9.271696
## 1943 SbiBTX642.01G204000 1-1 0 6.140779 -6.140779
## 5447 SbiBTX642.02G010300 1-1 0 5.536985 -5.536985
## 8174 SbiBTX642.02G304600 not1-1 0 5.373648 -5.373648
## 13480 SbiBTX642.03G436300 1-1 0 4.966707 -4.966707
## 11939 SbiBTX642.03G266300 1-1 0 4.925999 -4.925999
## 22121 SbiBTX642.06G215400 1-1 0 4.694880 -4.694880
## 23726 SbiBTX642.07G071400 not1-1 0 3.894962 -3.894962
## 23727 SbiBTX642.07G071400 not1-1 0 3.894962 -3.894962
## 23056 SbiBTX642.07G007500 1-1 0 3.656496 -3.656496Density Plots
Histogram of differences
