## 3 Results
### Differential Gene Expression Analysis  
To investigate the P impact on the gene expression across the different seed maturity stages, the following comparisons were conducted; S1-S2, S1-S3 and S2-S3. Regarding the P treatments, when comparisons were made, the P0 was used as the control treatment. Analysis of differential gene expression unveiled a grand total of 83,840 DEGs (DEGs) across all three P treatments, with 19,894 of these being unique. Within each P treatment group, we identified the following numbers of total DEGs and unique DEGs: the P0 treatment had 29,388 total DEGs and 15,680 unique DEGs; the P1 treatment had 28,313 total DEGs, with 15,218 unique DEGs; and the P2 treatment had 26,139 total DEGs, of which 14,385 were unique (Table 2).  
\begin{table}
\begin{tabular}{c c c c c c c c c c} \hline  & & **Total** & **Timmed** & **Total** & **Uniquely** & **Total** & **Mapped** & **Unassigned** & **Total** & **GC** \\  & & **Pair-End** & **Reads** & **Mapped** & **Mapped** & **Read** & **Read** & **Read** & **Number of** & **Content** \\  & & **Reads** & **Reads** & **Reads** & **Counts** & **Counts** & **Counts** & **Counts** & **Sequences** & **(\%)** \\ \hline  & S1 & 50,086,008 & 118,230 & 32,426,409 & 23,685,889 & 54,353,049 & 40,701,860 & 13,651,189 & 165,024,834 & 45.42 \\ P0 & S2 & 65,599,460 & 18,419 & 49,903,703 & 37,588,444 & 105,556,965 & 68,655,129 & 36,901,836 & 231,006,326 & 41.5 \\  & S3 & 71,596,343 & 20,118 & 53,199,567 & 42,520,287 & 92,483,302 & 75,966,636 & 16,516,666 & 249,591,820 & 42.83 \\  & S1 & 65,682,850 & 66,212 & 50,554,374 & 36,592,221 & 87,648,388 & 63,465,738 & 24,191,650 & 232,474,448 & 44.25 \\ P1 & S2 & 71,084,675 & 26,382 & 56,404,380 & 41,925,419 & 115,486,415 & 77,279,411 & 38,190,734 & 254,978,110 & 42.08 \\  & S3 & 71,439,019 & 26,239 & 53,275,790 & 41,108,021 & 92,959,843 & 73,545,382 & 19,414,461 & 249,429,618 & 43.25 \\  & S1 & 62,031,641 & 78,202 & 45,407,960 & 35,883,609 & 80,993,238 & 63,143,521 & 17,849,717 & 214,879,202 & 44 \\ P2 & S2 & 71,636,281 & 17,285 & 55,215,203 & 43,125,086 & 119,371,267 & 79,339,202 & 40,032,065 & 253,702,968 & 41.08 \\  & S3 & 69,727,261 & 25,783 & 53,976,207 & 43,758,093 & 96,997,450 & 78,404,807 & 18,592,643 & 247,406,936 & 42.92 \\ \hline \end{tabular}
\end{table}
Table 1: Summarized report for all the samples per seed maturity stage (S1, S2 and S3) in each P treatment (P0, P1 and P2) from data obtained from the MultiQC report following RNA sequencing procedure.  
Figure 1: The bar chart indicates the differences among the pod maturity stages (S1, S2 and S3) in each P treatment (P0, P1 and P2) from the summarized data obtained from the MultiQC report following RNA sequencing procedure.  
Of the total 19,894 unique genes, 2432 were only identified in the P0 treatment, 1375 were found in P1, and 1099 in P2, while all three treatments shared a common 10,401 DEGs (Figure 2). The comparison between the S1 and S3 in the P0 treatment (S1-S3) retained the highest numbers of upregulated and downregulated DEGs, having 8044 downregulated and 4680 upregulated, while the lowest in-rank DEGs were revealed in P2 treatment between the first seed maturity stage and the second (S1-S2), with 3074 down-regulated and 1943 up-regulated DEGs (Figure 3 and Table 3).