1 材料与方法
1.1 材料
1.2 方法
2 结果与分析
2.1 测序质量分析
表1 转录组测序数据统计Tab. 1 Statistics of transcriptome sequencing data |
| 样品名 Sample | 原始数据 Raw datas | 过滤后数据(比例/%) Clean data (Rate/%) | N(比例/%) N (Rate/%) | Q20/% | Q30/% | GC/% |
|---|---|---|---|---|---|---|
| l-0-1 | 56 551 194 | 56 420 470 (99.77) | 924 (0.00) | 97.87 | 93.87 | 54.88 |
| l-0-2 | 44 989 766 | 44 899 834 (99.80) | 104 (0.00) | 97.54 | 93.05 | 55.06 |
| l-0-3 | 62 097 412 | 61 943 386 (99.75) | 862 (0.00) | 97.75 | 93.53 | 54.99 |
| l-6-1 | 57 008 738 | 56 859 498 (99.70) | 704 (0.00) | 97.5 | 93.81 | 53.87 |
| l-6-2 | 44 767 702 | 44 669 556 (99.78) | 604 (0.00) | 98.09 | 94.40 | 55.53 |
| l-6-3 | 46 571 128 | 46 447 750 (99.74) | 582 (0.00) | 97.88 | 93.94 | 53.94 |
| l-24-1 | 44 599 080 | 44 470 182 (99.71) | 554 (0.00) | 97.81 | 93.79 | 55.51 |
| l-24-2 | 48 762 918 | 48 634 878 (99.74) | 162 (0.00) | 97.34 | 92.78 | 55.55 |
| l-24-3 | 58 393 242 | 58 274 628 (99.80) | 1454 (0.00) | 98.10 | 94.36 | 55.77 |
| m-0-1 | 49 707 062 | 49 606 202 (99.80) | 976 (0.00) | 98.06 | 94.33 | 54.97 |
| m-0-2 | 48 996 974 | 48 879 626 (99.76) | 836 (0.00) | 97.86 | 93.81 | 54.70 |
| m-0-3 | 52 782 284 | 52 660 696 (99.77) | 952 (0.00) | 97.93 | 93.99 | 54.56 |
| m-6-1 | 51 771 600 | 51 633 350 (99.73) | 850 (0.00) | 97.66 | 93.35 | 54.72 |
| m-6-2 | 45 978 404 | 45 868 228 (99.76) | 718 (0.00) | 97.63 | 93.26 | 55.35 |
| m-6-3 | 52 962 004 | 52 838 796 (99.77) | 886 (0.00) | 97.92 | 94.01 | 54.67 |
| m-24-1 | 42 810 864 | 42 727 702 (99.81) | 682 (0.00) | 97.95 | 93.96 | 54.55 |
| m-24-2 | 48 389 156 | 48 277 460 (99.77) | 112 (0.00) | 97.40 | 92.77 | 54.38 |
| m-24-3 | 46 940 194 | 46 820 510 (99.75) | 900 (0.00) | 97.84 | 93.78 | 54.61 |
注:N:不确定碱基数;Q20/%、Q30/%:Phred数值大于20、30碱基的数在碱基总数中的百分含量;GC/%:G和C碱基数量的总和在碱基总数的百分比;样品名称中l、m表示58L、58M品系,0、6、24表示0、6、24 h取样时间点,1、2、3表示3个重复。 | |
Note: N: Uncertain bases; Q20/%, Q30/%: The percentage of phred value is greater than 20, 30 bases in total bases; GC/%: The percentage of sum of G and C bases in total bases; In the sample name, l and m represent the 58L and 58M varieties, 0、6、24 represent the sampling time points of 0, 6, 24 h, and 1, 2, 3 represent the three repetitions. |
2.2 样本关系分析
图1 测序样品关系A:所有样品间PCA分析;B:测序样品相关性热图;C:58M品系测序样品间PCA分析;D:58L品系测序样品间PCA分析;E:样品间重复性散点图(以样品l-0-1 vs l-0-2为例);l、m分别表示58L、58M品系,0、6、24分别表示0、6、24 h取样时间点,1、2、3表示3个重复。 Fig. 1 Sequencing sample relationship A: PCA analysis between all samples; B: Correlation heatmap of sequencing samples; C: PCA analysis between m variety sequencing samples; D: PCA analysis between l variety sequencing samples; E: Repeatability scatter plot between samples (in order to sample l-0-1 vs l-0-2 as an example); l and m represent the 58L and 58M varieties, 0, 6, 24 represent the sampling time points of 0, 6, 24 h, and 1, 2, 3 represent the three repetitions. |
2.3 样品差异基因表达情况
图2 样品差异基因韦恩图A:盐胁迫6 h后差异表达基因韦恩图;B:盐胁迫24 h后差异表达基因韦恩图;l、m表示58L、58M品种,0、6、24表示0、6、24 h取样时间点;up表示基因表达量上调的组分,down表示基因表达量下调的组分。 Fig. 2 Venn diagram of differential expression gene in samples A: Venn diagram of differentially expressed genes after 6 h of salt stress; B: Venn diagram of differentially expressed genes after 24 h of salt stress ; l and m represent the 58L and 58M varieties, 0, 6, 24 represent the sampling time points of 0, 6, 24 h; up indicates the component whose gene expression level is up-regulated, and down indicates the component whose gene expression level is down-regulated. |
2.4 目标基因的GO功能富集分析
2.5 目标基因的KEEG代谢途径富集分析
2.6 目标基因转录因子分析
表2 盐胁迫下目标基因中转录因子的表达量Tab. 2 Expression of transcription factors in target genes under salt stress |
| 转录因子 Transcription factor | 58M | 58L | ||||||||
|---|---|---|---|---|---|---|---|---|---|---|
| 0 h | 6 h | 比例*/% | 24 h | 比例*/% | 0 h | 6 h | 比例*/% | 24 h | 比例*/% | |
| NAC073 | 17.30 | 35.97 | 107.90 | 19.80 | 14.45 | 4.39 | 1.83 | -58.35 | 2.51 | -42.85 |
| NAC045 | 2.17 | 7.44 | 243.33 | 2.06 | -4.94 | 1.42 | 0.70 | -50.49 | 1.07 | -24.86 |
| MYB2 | 6.41 | 13.05 | 103.73 | 7.06 | 10.15 | 2.47 | 0.65 | -73.68 | 3.87 | 56.80 |
| MYB308 | 3.04 | 11.37 | 273.55 | 8.40 | 176.14 | 3.99 | 2.71 | -32.21 | 1.76 | -55.92 |
| bHLH51 | 0.63 | 1.32 | 108.53 | 1.46 | 130.17 | 1.07 | 0.44 | -58.48 | 1.78 | 66.82 |
注:*:与对照(0 h)相比,上调比例。 | |
Note: *: The up-regulation ratio compared with the control (0 h). |
