欢迎访问作物学报,今天是

作物学报 ›› 2021, Vol. 47 ›› Issue (12): 2423-2439.doi: 10.3724/SP.J.1006.2021.02084

• 作物遗传育种·种质资源·分子遗传学 • 上一篇    下一篇

基于多平台基因表达数据的水稻干旱和盐胁迫相关基因预测

刘亚文1(), 张红燕1,2,*(), 曹丹2, 李兰芝2   

  1. 1湖南农业大学信息与智能科学技术学院, 湖南长沙 410128
    2湖南农业大学 / 湖南省农业大数据分析与决策工程技术研究中心, 湖南长沙 410128
  • 收稿日期:2020-11-30 接受日期:2021-04-26 出版日期:2021-12-12 网络出版日期:2021-06-02
  • 通讯作者: 张红燕
  • 作者简介:E-mail: lyw20201022@163.com
  • 基金资助:
    湖南省教育厅科学研究重点项目(18A105);长沙市工业科技特派员项目(201845);湖南农业大学“双一流”建设项目(SYL2019075)

Prediction of drought and salt stress-related genes in rice based on multi-platform gene expression data

LIU Ya-Wen1(), ZHANG Hong-Yan1,2,*(), CAO Dan2, LI Lan-Zhi2   

  1. 1College of Information and Intelligence, Hunan Agricultural University, Changsha 410128, Hunan, China
    2Hunan Engineering and Technology Research Centre for Agricultural Big Data Analysis and Decision-making, Hunan Agricultural University, Changsha 410128, Hunan, China
  • Received:2020-11-30 Accepted:2021-04-26 Published:2021-12-12 Published online:2021-06-02
  • Contact: ZHANG Hong-Yan
  • Supported by:
    Key Scientific Research Project of Hunan Education Department(18A105);Special Commissioner Project of Changsha City for Industrial Science and Technology(201845);“Double First-class” Construction Project of Hunan Agricultural University(SYL2019075)

摘要:

基于多平台基因表达数据挖掘水稻胁迫相关基因, 可增加关键基因预测的可靠性, 获得更具普适意义的结果。本研究从NCBI数据库中收集了与水稻非生物胁迫相关的94份affymetrix基因芯片数据和42份RNA-seq转录组数据。首先对同一类型同一胁迫相关的多个数据集以数据转换法融合, 得到干旱胁迫相关的affymetrix数据集D_affy和RNA-seq数据集D_rnaseq, 盐胁迫相关的affymetrix数据集S_affy和RNA-seq数据集S_rnaseq; 接着对4个数据集分别基于Pearson线性相关系数的经典WGCNA法和基于MIC非线性相关系数的改进WGCNA法进行基因共表达网络分析, 共获取胁迫相关的8个Hub基因集; 进一步, 对同一胁迫相关的Hub基因进行整合分析, 得到最终的水稻干旱胁迫相关Hub基因1936个、盐胁迫相关的Hub基因1504个。最后, 从预测性能、富集分析、文献报道、STRING在线互作分析和Cytoscape可视化分析等多角度解析Hub基因的生物学意义。结果显示: Hub基因整体预测性能较优, 且大多富集到了与干旱/盐胁迫相关的通路上, 其中有文献已报道的干旱胁迫响应基因31个和盐胁迫响应基因22个。此外, 通过对Hub基因的互作分析, 预测得到11个干旱胁迫候选基因和5个盐胁迫候选基因。本研究为“高维度、小样本”的农作物基因测序数据的有效分析提供了新思路, 实验结果为抗逆水稻品种研究提供了参考。

关键词: 水稻, 多平台, 干旱胁迫, 盐胁迫, WGCNA-MIC

Abstract:

Mining stress-related genes based on multi-platform gene expression data in rice can increase the reliability of key genes prediction and obtain more universally meaningful results. In this study, 94 affymetrix microarray data and 42 RNA-seq transcriptome data related to rice abiotic stress were collected from NCBI databases. First, multiple datasets related to the same stress on the same type were fused by data conversion method to obtain the affymetrix data set D_affy and RNA-seq data set D_rnaseq related to drought stress, and the affymetrix data set S_affy and the RNA-seq data set S_rnaseq related to salt stress. Then, the four datasets were analyzed by the classical WGCNA method based on Pearson's linear correlation coefficient and the improved WGCNA method based on the MIC nonlinear correlation coefficient respectively, and the eight Hub gene sets related to stress were obtained. Further, the integration analysis of stress-related Hub genes yielded the final 1936 drought stress-related Hub genes and 1504 salt stress-related Hub genes. Finally, the biological significance of Hub gene was analyzed from multiple perspectives, including prediction performance, enrichment analysis, literature report, STRING online interaction analysis, and Cytoscape visualization analysis. The results revealed that the overall prediction performance of Hub genes was better, and most of them were enriched in the pathways related to drought/salt stress. Among them, there were 31 drought stress response genes and 22 salt stress response genes reported in the literatures. In addition, 11 drought stress candidate genes and 5 salt stress candidate genes were predicted using the interaction analysis of Hub genes. In conclusion, This study provides a new idea for the effective analysis of “high-dimensional, small-sample” crop gene sequencing data, and the experimental results provide a reference for the study of stress-resistant rice varieties.

Key words: rice, multi-platform, drought stress, salt stress, WGCNA-MIC

图1

水稻affy数据处理流程"

附表1

来源于NCBI的Affymetrix基因芯片数据数据"

胁迫
Stress
GEO关联 GEO accession 测序平台
Platform
样本数(对照组/胁迫组)
Samples (control/stress)
干旱
Drought
GSE6901 Affymetrix Rice Genome Array (GPL2025) N = 6 (3/3 drought)
GSE21651 Affymetrix Rice Genome Array (GPL2025) N = 8 (4/4 drought)
GSE23211 Affymetrix Rice Genome Array (GPL2025) N = 12 (6/6 drought)
GSE26280 Affymetrix Rice Genome Array (GPL2025) N = 36 (18/18 drought)

Salt
GSE6901 Affymetrix Rice Genome Array (GPL2025) N = 6 (3/3 salt)
GSE14403 Affymetrix Rice Genome Array (GPL2025) N = 18 (9/9 salt)
GSE16108 Affymetrix Rice Genome Array (GPL2025) N = 8 (4/4 salt)

附表2

来源于NCBI的RNA-seq数据(RNA-seq)"

样品编号
Library name
样品描述
Library description
样品类型
Library type
单/双端
Library layout
水稻品种
Rice genotype
SRR3647331 1-leaf_34_day_Salt RNAseq Single Nipponbare
SRR3647329 1-leaf_34_day_Salt RNAseq Single Nipponbare
SRR3647327 1-leaf_34_day_Salt RNAseq Single Nipponbare
SRR3647330 1-leaf_34_day_Control RNAseq Single Nipponbare
SRR3647328 1-leaf_34_day_Control RNAseq Single Nipponbare
SRR3647326 1-leaf_34_day_Control RNAseq Single Nipponbare
ERR266228 1-Seedling_shoots_2_weeks_1h_Control RNAseq Single Nipponbare
ERR266233 1-Seedling_shoots_2_weeks_1h_Control RNAseq Single Nipponbare
ERR266230 1-Seedling_shoots_2_weeks_1h_Control RNAseq Single Nipponbare
ERR266225 1-Seedling_shoots_2_weeks_24h_Control RNAseq Single Nipponbare
ERR266234 1-Seedling_shoots_2_weeks_24h_Control RNAseq Single Nipponbare
ERR266232 1-Seedling_shoots_2_weeks_24h_Control RNAseq Single Nipponbare
ERR266229 1-Seedling_shoots_2_weeks_5h_Control RNAseq Single Nipponbare
ERR266223 1-Seedling_shoots_2_weeks_5h_Control RNAseq Single Nipponbare
ERR266222 1-Seedling_shoots_2_weeks_5h_Control RNAseq Single Nipponbare
ERR266237 1-Seedling_shoots_2_weeks_1h_Salt RNAseq Single Nipponbare
ERR266236 1-Seedling_shoots_2_weeks_1h_Salt RNAseq Single Nipponbare
ERR266235 1-Seedling_shoots_2_weeks_1h_Salt RNAseq Single Nipponbare
ERR266238 1-Seedling_shoots_2_weeks_24h_Salt RNAseq Single Nipponbare
ERR266226 1-Seedling_shoots_2_weeks_24h_Salt RNAseq Single Nipponbare
ERR266231 1-Seedling_shoots_2_weeks_24h_Salt RNAseq Single Nipponbare
ERR266227 1-Seedling_shoots_2_weeks_5h_Salt RNAseq Single Nipponbare
ERR266224 1-Seedling_shoots_2_weeks_5h_Salt RNAseq Single Nipponbare
ERR266221 1-Seedling_shoots_2_weeks_5h_Salt RNAseq Single Nipponbare
SRR7054183 2-Inflorescence_Control RNAseq Paired Nipponbare
SRR7054182 2-Inflorescence_Control RNAseq Paired Nipponbare
SRR7054181 2-Inflorescence_Control RNAseq Paired Nipponbare
SRR7054180 2-Inflorescence_Control RNAseq Paired Nipponbare
SRR7054179 2-Inflorescence_Drought RNAseq Paired Nipponbare
SRR7054178 2-Inflorescence_Drought RNAseq Paired Nipponbare
SRR7054177 2-Inflorescence_Drought RNAseq Paired Nipponbare
SRR7054176 2-Inflorescence_Drought RNAseq Paired Nipponbare
SRR3051752 Drought stress_rep1 RNAseq Paired Nipponbare
SRR3051753 Drought stress_rep2 RNAseq Paired Nipponbare
SRR3051754 Drought stress_rep3 RNAseq Paired Nipponbare
SRR3051755 Well-water_rep1 RNAseq Paired Nipponbare
SRR3051756 Well-water_rep2 RNAseq Paired Nipponbare
SRR3051757 Well-water_rep3 RNAseq Paired Nipponbare
SRR3051740 Well-water_rep1 RNAseq Paired Nipponbare
SRR3051741 Well-water_rep2 RNAseq Paired Nipponbare
SRR3051742 Well-water_rep3 RNAseq Paired Nipponbare
SRR3051743 Drought stress_rep1 RNAseq Paired Nipponbare
SRR3051744 Drought stress_rep2 RNAseq Paired Nipponbare
SRR3051745 Drought stress_rep3 RNAseq Paired Nipponbare

图2

水稻RNA-seq数据处理流程"

表1

水稻数据集"

数据集
Data set
基因数
No. of genes
总样本数
Total samples
对照组样本数
Control samples
胁迫组样本数
Stress samples
干旱芯片数据
Drought stress-related affymetrix dataset (D_affy)
27,344 62 31 31
盐芯片数据
Salt stress-related affymetrix dataset (S_affy)
27,344 32 16 16
干旱RNA-seq数据
Drought stress-related RNA-seq dataset (D_rnaseq)
29,828 20 10 10
盐RNA-seq数据
Salt stress-related RNA-seq dataset (S_rnaseq)
28,425 22 11 11

图3

D_affy数据WGCNA网络分析结果 A、B: 分别基于WGCNA-P和WGCNA-MIC方法的基因聚类树和模块划分, Dynamic Tree Cut表示由原始计算划分的模块, Merged dynamic表示合并后的结果; C、D: 分别基于WGCNA-P和WGCNA-MIC方法的各模块的ME和MS值, 红色表示模块与干旱胁迫正相关, 蓝色表示模块与干旱胁迫负相关。"

图4

D_rnaseq数据WGCNA网络分析结果 A、B: 分别基于WGCNA-P和WGCNA-MIC方法的基因聚类树和模块划分, Dynamic Tree Cut表示由原始计算划分的模块, Merged dynamic表示合并后的结果; C、D: 分别基于WGCNA-P和WGCNA-MIC方法的各模块的ME和MS值, 红色表示模块与干旱胁迫正相关, 蓝色表示模块与干旱胁迫负相关。"

附图1

S_affy数据WGCNA网络分析结果 A、B: 分别基于WGCNA-P和WGCNA-MIC方法的基因聚类树和模块划分, Dynamic Tree Cut表示由原始计算划分的模块, Merged dynamic表示合并后的结果; C、D: 分别基于WGCNA-P和WGCNA-MIC方法的各模块的ME和MS值, 红色表示模块与盐胁迫正相关, 蓝色表示模块与盐胁迫负相关。"

附图2

S_rnaseq数据WGCNA网络分析结果 A、B: 分别基于WGCNA-P和WGCNA-MIC法的基因聚类树和模块划分, Dynamic Tree Cut表示由原始计算划分的模块, Merged dynamic表示合并后的结果; C、D: 分别基于WGCNA-P和WGCNA-MIC法的各模块的ME和MS值, 红色表示模块与盐胁迫正相关, 蓝色表示模块与盐胁迫负相关。"

表2

Hub基因的分类精度"

胁迫
Stress
Hub基因集
Hub gene set
基因数目
Number of genes
数据集
Data set
平均精度
Average accuracy (%)
干旱Drought D_affy_P 220 D_affy 100
D_affy_MIC 104 D_affy 100
D_rnaseq_P 738 D_rnaseq 90.0
D_rnaseq_MIC 1634 D_rnaseq 91.0
D_meta_hub 1936 D_affy 100
D_rnaseq 96.0
盐Salt S_affy_P 470 S_affy 100
S_affy_MIC 293 S_affy 100
S_rnaseq_P 684 S_rnaseq 81.0
S_rnaseq_MIC 331 S_rnaseq 84.6
S_meta_hub 1504 S_affy 100
S_rnaseq 84.6

表3

Hub基因的GO富集部分分析结果"

胁迫
Stress
GO条目
GO term
基因数目
Number of genes
基因本体
Ontology
描述
Description
P
P-value
干旱
Drought
GO:0010033 12 BP 对有机物的响应
Response to organic substance
9.00E-06
GO:0009719 12 BP 内源性刺激响应
Response to endogenous stimulus
9.00E-06
GO:0009725 12 BP 激素刺激响应
Response to hormone stimulus
9.00E-06
GO:0006721 8 BP 萜类化合物代谢过程
Terpenoid metabolic process
2.80E-05
GO:0007165 37 BP 信号传导Signal transduction 2.20E-05
GO:0009628 15 BP 对非生物刺激的响应
Response to abiotic stimulus
0.00054
GO:0006970 5 BP 渗透胁迫响应
Response to osmotic stress
0.0036
GO:0009415 5 BP 对水的响应Response to water 0.006
GO:0004872 22 MF 受体活性 Receptor activity 9.40E-13
GO:0004713 21 MF 蛋白质酪氨酸激酶活性
Protein tyrosine kinase activity
9.00E-12
GO:0004722 15 MF 蛋白丝氨酸/苏氨酸磷酸酶活性
Protein serine/threonine phosphatase activity
2.70E-05
GO:0008135 12 MF 翻译因子活性, 核酸结合
Translation factor activity, Nucleic acid binding
0.0099
GO:0044424 1120 CC 细胞内成分 Intracellular part 0
GO:0005737 1011 CC 细胞质Cytoplasm 0
GO:0016020 278 CC 薄膜 Membrane 3.10E-23

Salt
GO:0070887 6 BP 细胞对化学刺激的响应
Cellular response to chemical stimulus
2.00E-06
GO:0007275 13 BP 多细胞有机体的发育
Multicellular organismal development
5.30E-08
GO:0043436 52 BP 草酸代谢过程
Oxoacid metabolic process
4.80E-06
GO:0019752 52 BP 羧酸代谢过程
Carboxylic acid metabolic process
4.80E-06
GO:0006082 52 BP 有机酸代谢过程
Organic acid metabolic process
5.10E-06
GO:0010033 11 BP 对有机物的响应
Response to organic substance
4.70E-06
GO:0009719 11 BP 内源性刺激响应
Response to endogenous stimulus
4.70E-06
GO:0009416 6 BP 光刺激响应
Response to light stimulus
0.024
GO:0004872 14 MF 受体活性 Receptor activity 5.70E-08
GO:0008135 18 MF 翻译因子活性, 核酸结合
Translation factor activity, nucleic acid binding
1.00E-06
GO:0003743 11 MF 翻译起始因子活性
Translation initiation factor activity
3.00E-05
GO:0016874 32 MF 连接酶活性 Ligase activity 0.00064
GO:0016020 225 CC 薄膜 Membrane 4.60E-21
GO:0009507 20 CC 叶绿体 Chloroplast 4.30E-15
GO:0005886 18 CC 细胞质膜 Plasma membrane 1.20E-09

表4

已报道与胁迫相关的Hub基因"

胁迫
Stress
基因编号RAP_locus 基因符号
Gene symbol
胁迫
Stress
基因编号
RAP_locus
基因符号
Gene symbol
干旱
Drought
Os05g0455500 OsP5CS; OsP5CS1; OsALDH18B1 干旱
Drought
Os03g0286900 OsRCI2-5
Os02g0766700 OsbZIP23 Os02g0149800 OsPP18
Os08g0112700 OsMADS26 Os03g0267000 OsHSP18.0-CI; OsMSR3; OsSHSP1
Os06g0130100 OsSIK1 Os05g0475400 OsAMTR1
Os09g0552300 OsRPK1 Os08g0408500 OsERF48; OsDRAP1
Os01g0867300 OsABF1; OsbZIP12 Os04g0676700 OsMYB6
Os03g0125100 DSM2 Os12g0597500 OsUAH
Os03g0745000 OsHsfA2a Os06g0316000 Os2H16
Os05g0542500 OsLEA3; OsLEA3-1 Os06g0612800 OsiSAP8
Os02g0671100 MAIF1 Os04g0572400 OsDREB1E
Os06g0211200 OsAREB1; OsbZIP46; OsABF2; ABL1 Os11g0126900 OsNAC10; ONAC122
Os05g0437700 EDT1; OsbZIP40 Os11g0707600 OsGL1-11
Os05g0569300 OsbZIP45 Os03g0805100 SQS
Os08g0196700 OsNF-YA7; OsHAP2A Os05g0213500 OsPYL/RCAR5; OsPYL5; OsPYL11
Os03g0230300 OsSRO1c; BOC1 Os06g0598800 WSL1
Os04g0541700 Oshox22
盐Salt Os03g0348900 OsSRFP1; SDEL2 盐Salt Os02g0121300 OsCYP2; LRT2
Os04g0652400 OsSULTR3; 3;lpa Os03g0272300 OsSDIR1
Os03g0329900 OsPHR1 Os07g0129200 OsPR1a; OsSCP
Os07g0187700 OsPHF1 Os05g0437700 EDT1; OsbZIP40
Os02g0678200 OsSPX-MFS2; OsPSS2 Os08g0557000 OsPIMT1
Os02g0325600 NIGT1 Os06g0693700 OsSIDP366
Os01g0755700 NBIP1 Os04g0676700 OsMYB6
Os10g0545700 OsACR2.1 Os01g0612700 OsLOL2; OsLOL5
Os01g0869900 OsSAPK4; OSPDK Os01g0948400 OsP5CR
Os03g0719900 OsPTR8; OsNPF8.5 Os03g0319300 OsCam1-1
Os09g0434500 OsBIERF1 Os05g0584200 OsLEA5

图5

干旱胁迫相关基因互作网络"

图6

盐胁迫相关基因互作网络"

表5

候选基因在STRING中的注释"

胁迫
Stress
候选基因
Candidate gene
STRING中名称
Name in STRING
注释
Annotation
干旱Drought Os01g0733200 HSF11 热应激转录因子C-1b; 与热休克启动子元件(HSE)的DNA特异性结合的转录调控因子。
Heat stress transcription factor C-1b; transcriptional regulator that specifically binds DNA of heat shock promoter elements (HSE).
Os07g0178600 HSF5 热应激转录因子A-2b; 转录调控因子, 可特异性结合热休克启动子元件(HSE)的DNA; 属于HSF家族, A类亚科。
Heat stress transcription factor A-2b; transcriptional regulator that specifically binds DNA of heat shock promoter elements (HSE); belongs to the HSF family. Class A subfamily.
Os09g0526600 HSFB2C 热应激转录因子B-2c; 转录调控因子, 可特异性结合热休克启动子元件(HSE)的DNA; 属于HSF家族, B类亚科。
Heat stress transcription factor B-2c; transcriptional regulator that specifically binds DNA of heat shock promoter elements (HSE); belongs to the HSF family. Class B subfamily.
Os01g0583100 OS01T0583100-01 可能的蛋白质磷酸酶2C 6; 属于PP2C家族。
Probable protein phosphatase 2C 6; belongs to the PP2C family.
Os03g0231700 OS03T0231700-02 Os03g0231700蛋白; 角鲨烯单加氧酶, 假定表达; cDNA克隆: J033045D18, 完整插入序列。
Os03g0231700 protein; Squalene monooxygenase, putative, expressed; cDNA clone: J033045D18, full insert sequence.
Os03g0376100 OS03T0376100-01 Os03g0376100蛋白。
Os03g0376100 protein.
胁迫
Stress
候选基因
Candidate gene
STRING中名称
Name in STRING
注释
Annotation
Os04g0107900 OS04T0107900-02 Os04g0107900蛋白。
Os04g0107900 protein.
Os01g0840100 OsJ_04024 cDNA克隆: J100050G20, 完整插入序列; 70 kD热激蛋白; Os01g0840100蛋白; 假定的HSP70; 未表征的蛋白质。
cDNA, clone: J100050G20, full insert sequence; 70 kD heat shock protein; Os01g0840100 protein; Putative HSP70; uncharacterized protein.
Os03g0277300 OsJ_10337 热休克同源70 kD蛋白, 假定表达。
Heat shock cognate 70 kD protein, putative, expressed.
Os05g0460000 OsJ_18811 Os05g0460000蛋白; 推定的hsp70; cDNA克隆: J090096I11, 完整插入序列; 属于热休克蛋白70家族。
Os05g0460000 protein; Putative hsp70; cDNA, clone: J090096I11, full insert sequence; Belongs to the heat shock protein 70 family.
Os06g0110200 OsJ_19858 cDNA克隆: 002-135-D09, 完整插入序列; Os06g0110200蛋白; 假定的未表征蛋白OSJNBa0004I20.22; 假定的未表征蛋白P0514G12.46。
cDNA clone: 002-135-D09, full insert sequence; Os06g0110200 protein; putative uncharacterized protein OSJNBa0004I20.22; putative uncharacterized protein P0514G12.46.

Salt
Os06g0727200 CATB 过氧化氢酶同工酶B; 几乎发生在所有有氧呼吸的生物中, 并保护细胞免受过氧化氢的毒性作用。
Catalase isozyme B; occurs in almost all aerobically respiring organisms and serves to protect cells from the toxic effects of hydrogen peroxide.
Os12g0502300 CYCA2-1 Cyclin-A2-1; 属于细胞周期蛋白家族。Cyclin AB亚家族。
Cyclin-A2-1; belongs to the cyclin family. Cyclin AB subfamily.
Os03g0821100 OsJ_13143 热休克同源70 kD蛋白2, 假定表达; 热休克蛋白同源物70; Os03g0821100蛋白; cDNA克隆: J023030D03, 完整插入序列。
Heat shock cognate 70 kD protein 2, putative, expressed; heat shock protein cognate 70; Os03g0821100 protein; cDNA clone:J023030D03, full insert sequence.
Os10g0491801 OsJ_31995 假定泛素/核糖体蛋白S27a融合蛋白;泛素融合蛋白, 假定表达。
Putative ubiquitin/ribosomal protein S27a fusion protein; ubiquitin fusion protein, putative, expressed.
Os02g0775200 RFC3 复制因子C亚基3; 可能参与DNA复制, 从而调节细胞增殖。
Replication factor C subunit 3; may be involved in DNA replication and thus regulate cell proliferation.
[1] Gong Z Z, Xiong L M, Shi H Z, Yang S H, Herrera-Estrella L R, Xu G H, Chao D Y, Li J R, Wang P Y, Qin F, Li J G, Ding Y L, Shi Y T, Wang Y, Yang Y Q, Guo Y, Zhu J K. Plant abiotic stress response and nutrient use efficiency. Sci China (Life Sci Edn), 2020, 63:635-674.
[2] Hossain M R, Bassel G W, Pritchard J, Sharma G P, Ford-Lloyd B V. Trait specific expression profiling of salt stress responsive genes in diverse rice genotypes as determined by modified significance analysis of microarrays. Front Plant Sci, 2016, 7:567.
[3] Zhou Y, Yang P, Cui F L, Zhang F T, Luo X D, Xie J K. Transcriptome analysis of salt stress responsiveness in the seedlings of Dongxiang wild rice (Oryza rufipogon Griff.). PLoS One, 2016, 11(1):e0146242.
doi: 10.1371/journal.pone.0146242
[4] Song T, Das D, Yang F, Chen M X, Tian Y, Cheng C L, Sun C, Xu W F, Zhang J H. Genome-wide transcriptome analysis of roots in two rice varieties in response to alternate wetting and drying irrigation. Crop J, 2020, 8:586-601.
doi: 10.1016/j.cj.2020.01.007
[5] Kaur S, Iquebal M A, Jaiswal S, Tandon G, Sundaram R M, Gautam R K, Suresh K P, Rai A, Kumar D. A meta-analysis of potential candidate genes associated with salinity stress tolerance in rice. Agric Gene, 2016, 1:126-134.
[6] Serin E A R, Nijveen H, Hilhorst H W M, Ligterink W. Learning from co-expression networks: possibilities and challenges. Front Plant Sci, 2016, 7:444.
[7] Zhou G Y, Soufan O, Ewald J, Hancock R E W, Basu N, Xia J G. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res, 2019, 47(W1):W234-W241.
doi: 10.1093/nar/gkz240
[8] 王凯莉, 张礼, 刘学军. 融合多平台表达数据的转录组差异表达分析. 计算机学报, 2018, 41:1195-1210.
Wang K L, Zhang L, Liu X J. Differential expression analysis based on integrating transcriptome expression data from multiple platforms. Chin J Comp, 2018, 41:1195-1210 (in Chinese with English abstract).
[9] Cheng Y F, Li L, Qin Z S, Li X, Qi F. Identification of castration-resistant prostate cancer-related hub genes using weighted gene co-expression network analysis. J Cell Mol Med, 2020, 24:1-12.
doi: 10.1111/jcmm.v24.1
[10] 李旭凯, 李任建, 张宝俊. 利用WGCNA鉴定非生物胁迫相关基因共表达网络. 作物学报, 2019, 45:1349-1364.
Li X K, Li R J, Zhang B J. Identification of rice stress-related gene co-expression modules by WGCNA. Acta Agron Sin, 2019, 45:1349-1364 (in Chinese with English abstract).
[11] Zhu M D, Xie H J, Wei X J, Dossa K, Yu Y Y, Hui S J, Tang G H, Zeng X S, Yu Y H, Hu P S, Wang J L. WGCNA analysis of salt-responsive core transcriptome identifies novel hub genes in rice. Genes, 2019, 10:719.
doi: 10.3390/genes10090719
[12] Lv Y M, Xu L, Dossa K, Zhou K, Zhu M D, Xie H J, Tang S J, Yu Y Y, Guo X Y, Zhou B. Identification of putative drought-responsive genes in rice using gene co-expression analysis. Bioinformation, 2019, 15:480-488.
doi: 10.6026/bioinformation
[13] Hopper D W, Ghan R, Schlauch K A, Cramer G R. Transcriptomic network analyses of leaf dehydration responses identify highly connected ABA and ethylene signaling hubs in three grapevine species differing in drought tolerance. BMC Plant Biol, 2016, 16:118.
doi: 10.1186/s12870-016-0804-6
[14] 秦天元, 孙超, 毕真真, 梁文君, 李鹏程, 张俊莲, 白江平. 基于WGCNA的马铃薯根系抗旱相关共表达模块鉴定和核心基因发掘. 作物学报, 2020, 46:1033-1051.
Qin T Y, Sun C, Bi Z Z, Liang W J, Li P C, Zhang J L, Bai J P. Identification of drought-related co-expression modules and hub genes in potato roots based on WGCNA. Acta Agron Sin, 2020, 46:1033-1051 (in Chinese with English abstract).
[15] Reshef D N, Reshef Y A, Finucane H K, Grossman S R, McVean G, Turnbaugh P J, Lander E S, Mitzenmacher M, Sabeti P C. Detecting novel associations in large data sets. Science, 2011, 334:1518-1524.
doi: 10.1126/science.1205438 pmid: 22174245
[16] Britt C L, Weisburd D. Statistical Power. Handbook of Quantitative Criminology. New York: Springer, 2010. pp 313-32.
[17] Durinck S, Spellman P T, Birney E, Huber W. Mapping identifiers for the integration of genomic datasets with the R/Bioconductor package biomaRt. Nat Protoc, 2009, 4:1184-1191.
doi: 10.1038/nprot.2009.97
[18] Steffen D, Yves M, Arek K, Sean D, Bart D M, Alvis B, Wolfgang H. BioMart and Bioconductor: a powerful link between biological databases and microarray data analysis. Bioinformatics, 2005, 21:3439-3440.
pmid: 16082012
[19] Kroll K W, Mokaram N E, Pelletier A R, Frankhouser D E, Westphal M S, Stump P A, Stump C L, Bundschuh R, Blachly J S, Yan P. Quality control for RNA-Seq (QuaCRS): an integrated quality control pipeline. Cancer Inform, 2014, 13(S3):7-14.
[20] Chen S F, Zhou Y Q, Chen Y R, Gu J. Fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics, 2018, 34:i884-i890.
doi: 10.1093/bioinformatics/bty560
[21] Liao Y, Smyth G K, Shi W. FeatureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics, 2014, 30:923-930.
doi: 10.1093/bioinformatics/btt656
[22] Love M I, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genom Biol, 2014, 15:31-46.
[23] Smita S, Katiyar A, Pandey D M, Chinnusamy V, Archak S, Bansal K C. Identification of conserved drought stress responsive gene-network across tissues and developmental stages in rice. Bioinformation, 2013, 9:72-78.
doi: 10.6026/bioinformation
[24] Shaik R, Ramakrishna W. Machine learning approaches distinguish multiple stress conditions using stress-responsive genes and identify candidate genes for broad resistance in rice. Plant Physiol, 2014, 164:481-495.
doi: 10.1104/pp.113.225862 pmid: 24235132
[25] Tian T, Liu Y, Yan H Y, You Q, Yi X, Du Z, Xu W Y, Su Z. AgriGO v2.0: a GO analysis toolkit for the agricultural community. Nucleic Acids Res, 2017, 45(W1):W122-W129.
doi: 10.1093/nar/gkx382
[26] Szklarczyk D, Gable A L, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva N T, Morris J H, Bork P, Jensen L J, Mering C V. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res, 2019, 47(D1):D607-D613.
doi: 10.1093/nar/gky1131
[27] Shannon P, Markiel A, Ozier O, Baliga N S, Wang J T, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genom Res, 2003, 13:2498-2504.
doi: 10.1101/gr.1239303
[28] 王胜昌, 涂海甫, 胡丹, 吴奈, 岑祥, 熊立仲. 水稻抗非生物逆境功能基因的发掘. 生命科学, 2016, 28:1216-1229.
Wang S C, Tu H F, Hu D, Wu N, Cen X, Xiong L Z. The exploitation of rice functional genes for abiotic stress. Chin Bull Life Sci, 2016, 28:1216-1229 (in Chinese with English abstract).
[29] Zahedi S M, Karimi M, Venditti A. Plants adapted to arid areas: specialized metabolites. Nat Prod Res, 2019: 1-18.
[30] Murad M A, Khan A L, Muneer S. Silicon in horticultural crops: cross-talk, signaling, and tolerance mechanism under salinity stress. Plants, 2020, 9:460.
doi: 10.3390/plants9040460
[31] Sripinyowanich S, Klomsakul P, Boonburapong B, Bangyeekhun T, Asami T, Gu Y H, Buaboocha T, Chadchawan S. Exogenous ABA induces salt tolerance in indica rice (Oryza sativa L.): the role of OsP5CS1 and OsP5CR gene expression during salt stress. Environ Exp Bot, 2013, 86:94-105.
doi: 10.1016/j.envexpbot.2010.01.009
[32] Ngan K G, Pati P K, Richaud F, Parizot B, Bidzinski P, Mai C D, Bès M, Bourrié I, Meynard D, Beeckman T, Selvaraj M G, Manabu I, Genga A M, Brugidou C, Do V N, Guiderdoni E, Morel J B, Gantet P. OsMADS26 negatively regulates resistance to pathogens and drought tolerance in rice. Plant Physiol, 2015, 169:2935-2949.
[33] Zong W, Tang N, Yang J, Peng L, Ma S Q, Xu Y, Li G L, Xiong L Z. Feedback regulation of ABA signaling and biosynthesis by a bZIP transcription factor targets drought-resistance-related genes. Plant Physiol, 2016, 171:2810-2825.
doi: 10.1104/pp.16.00469 pmid: 27325665
[34] Ouyang S Q, Liu Y F, Liu P, Lei G, He S J, Ma B, Zhang W K, Zhang J S, Chen S Y. Receptor-like kinase OsSIK1 improves drought and salt stress tolerance in rice (Oryza sativa) plants. Plant J, 2010, 62:316-329.
doi: 10.1111/j.1365-313X.2010.04146.x
[35] Cheng Y W, Qi Y C, Zhu Q, Chen X, Wang N, Zhao X, Chen H Y, Cui X J, Xu L L, Zhang W. New changes in the plasma-membrane-associated proteome of rice roots under salt stress. Proteomics, 2009, 9:3100-3114.
doi: 10.1002/pmic.200800340
[36] Fang H M, Meng Q L, Xu J W, Tang H J, Tang S Y, Zhang H S, Huang J. Knock-down of stress inducible OsSRFP1 encoding an E3 ubiquitin ligase with transcriptional activation activity confers abiotic stress tolerance through enhancing antioxidant protection in rice. Plant Mol Biol, 2015, 87:441-458.
doi: 10.1007/s11103-015-0294-1
[37] Ouyang J, Cai Z Y, Xia K F, Wang Y Q, Duan J, Zhang M Y. Identification and analysis of eight peptide transporter homologs in rice. Plant Sci, 2010, 179:374-382.
doi: 10.1016/j.plantsci.2010.06.013
[38] Kothari K S, Dansana P K, Jitender G, Tyagi A K. Rice stress associated protein 1 (OsSAP1) interacts with aminotransferase (OsAMTR1) and pathogenesis-related 1a protein (OsSCP) and regulates abiotic stress responses. Front Plant Sci, 2016, 7:1057.
[39] Wei Y D, Xu H B, Diao L R, Zhu Y H, Xie H G, Cai Q H, Wu F X, Wang Z H, Zhang J F, Xie H A. Protein repair L-isoaspartyl methyl transferase 1 ( PIMT1) in rice improves seed longevity by preserving embryo vigor and viability. Plant Mol Biol, 2015, 89:475-492.
doi: 10.1007/s11103-015-0383-1
[40] He S, Tan L L, Hu Z L, Chen G P, Wang G X, Hu T Z. Molecular characterization and functional analysis by heterologous expression in E. coli under diverse abiotic stresses for OsLEA5, the atypical hydrophobic LEA protein from Oryza sativa L. Mol Genet Genom, 2012, 287:39-54.
doi: 10.1007/s00438-011-0660-x
[1] 田甜, 陈丽娟, 何华勤. 基于Meta-QTL和RNA-seq的整合分析挖掘水稻抗稻瘟病候选基因[J]. 作物学报, 2022, 48(6): 1372-1388.
[2] 郑崇珂, 周冠华, 牛淑琳, 和亚男, 孙伟, 谢先芝. 水稻早衰突变体esl-H5的表型鉴定与基因定位[J]. 作物学报, 2022, 48(6): 1389-1400.
[3] 周文期, 强晓霞, 王森, 江静雯, 卫万荣. 水稻OsLPL2/PIR基因抗旱耐盐机制研究[J]. 作物学报, 2022, 48(6): 1401-1415.
[4] 郑小龙, 周菁清, 白杨, 邵雅芳, 章林平, 胡培松, 魏祥进. 粳稻不同穗部籽粒的淀粉与垩白品质差异及分子机制[J]. 作物学报, 2022, 48(6): 1425-1436.
[5] 颜佳倩, 顾逸彪, 薛张逸, 周天阳, 葛芊芊, 张耗, 刘立军, 王志琴, 顾骏飞, 杨建昌, 周振玲, 徐大勇. 耐盐性不同水稻品种对盐胁迫的响应差异及其机制[J]. 作物学报, 2022, 48(6): 1463-1475.
[6] 杨建昌, 李超卿, 江贻. 稻米氨基酸含量和组分及其调控[J]. 作物学报, 2022, 48(5): 1037-1050.
[7] 杨德卫, 王勋, 郑星星, 项信权, 崔海涛, 李生平, 唐定中. OsSAMS1在水稻稻瘟病抗性中的功能研究[J]. 作物学报, 2022, 48(5): 1119-1128.
[8] 朱峥, 王田幸子, 陈悦, 刘玉晴, 燕高伟, 徐珊, 马金姣, 窦世娟, 李莉云, 刘国振. 水稻转录因子WRKY68在Xa21介导的抗白叶枯病反应中发挥正调控作用[J]. 作物学报, 2022, 48(5): 1129-1140.
[9] 王小雷, 李炜星, 欧阳林娟, 徐杰, 陈小荣, 边建民, 胡丽芳, 彭小松, 贺晓鹏, 傅军如, 周大虎, 贺浩华, 孙晓棠, 朱昌兰. 基于染色体片段置换系群体检测水稻株型性状QTL[J]. 作物学报, 2022, 48(5): 1141-1151.
[10] 王霞, 尹晓雨, 于晓明, 刘晓丹. 干旱锻炼对B73自交后代当代干旱胁迫记忆基因表达及其启动子区DNA甲基化的影响[J]. 作物学报, 2022, 48(5): 1191-1198.
[11] 雷新慧, 万晨茜, 陶金才, 冷佳俊, 吴怡欣, 王家乐, 王鹏科, 杨清华, 冯佰利, 高金锋. 褪黑素与2,4-表油菜素内酯浸种对盐胁迫下荞麦发芽与幼苗生长的促进效应[J]. 作物学报, 2022, 48(5): 1210-1221.
[12] 王泽, 周钦阳, 刘聪, 穆悦, 郭威, 丁艳锋, 二宫正士. 基于无人机和地面图像的田间水稻冠层参数估测与评价[J]. 作物学报, 2022, 48(5): 1248-1261.
[13] 陈悦, 孙明哲, 贾博为, 冷月, 孙晓丽. 水稻AP2/ERF转录因子参与逆境胁迫应答的分子机制研究进展[J]. 作物学报, 2022, 48(4): 781-790.
[14] 王吕, 崔月贞, 吴玉红, 郝兴顺, 张春辉, 王俊义, 刘怡欣, 李小刚, 秦宇航. 绿肥稻秆协同还田下氮肥减量的增产和培肥短期效应[J]. 作物学报, 2022, 48(4): 952-961.
[15] 巫燕飞, 胡琴, 周棋, 杜雪竹, 盛锋. 水稻延伸因子复合体家族基因鉴定及非生物胁迫诱导表达模式分析[J]. 作物学报, 2022, 48(3): 644-655.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!