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

作物学报 ›› 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] 胡赵, 钱润, 谢丰璞, 应素平. 水稻SPX基因家族鉴定及响应磷处理的表达分析[J]. 作物学报, 2026, 52(6): 1902-1912.
[2] 邹仪妹, 徐敏, 汪海洋, 姚辉, 王加峰, 刘浩, 任代胜. 两系不育系水稻幼苗根系响应盐胁迫的转录因子调控网络鉴定[J]. 作物学报, 2026, 52(6): 1728-1742.
[3] 闫安, 蒋昆炜, 王蓉圆, 田林, 张璐, 王韵, 徐建龙. 水稻剑叶小维管束数基因SVN7的鉴定与克隆[J]. 作物学报, 2026, 52(5): 1364-1372.
[4] 宋裕祯, Bheel Chander Kumar, 王跃, 张颖星, 郭娟, Khound Rituraj, Santra Dipak Kumar, 曹晓宁, 王瑞云. 糜子AP2亚家族全基因组鉴定及PmAP2-1和PmAP2-9耐盐功能分析[J]. 作物学报, 2026, 52(4): 1127-1139.
[5] 陈伟, 卫万娟, 赵其兵, 常东伟, 余凌波, 翟鹏飞, 冯志明, 陈宗祥, 任仰涛, 杨鹏, 刘海浪, 李珍富, 杨永乐, 金彦刚, 左示敏. 利用CRISPR/Cas9编辑Hd6基因创制优质早熟水稻新种质[J]. 作物学报, 2026, 52(4): 1046-1056.
[6] 石少阶, 刘凯, 陈姿夷, 王卉颖, 李三和, 周雷, 游艾青. 水稻矮化多分蘖基因DMT1的克隆与功能分析[J]. 作物学报, 2026, 52(4): 1022-1034.
[7] 覃奕琰, 付瑶, 苏畅, 李娜, 徐静茹, 程笑然, 张琪, 赵明辉. OsST41调控水稻苗期耐盐性的功能分析[J]. 作物学报, 2026, 52(3): 802-812.
[8] 叶凡, 李帅, 李思宇, 陈云, 窦超银, 刘立军. 不同节水灌溉方式对东北稻区水稻产量和群体质量的影响[J]. 作物学报, 2026, 52(3): 895-907.
[9] 张晴, 杨昱, 郭茜, 岳霈尧, 殷丛丛, 牛景萍, 赵晋忠, 杜维俊, 岳爱琴. 大豆GmARA6a的克隆及响应盐胁迫的功能分析[J]. 作物学报, 2026, 52(2): 480-493.
[10] 杨飚, 杜帅康, 张继旺, 石瑛, 张丽莉. 马铃薯III类POD基因家族的全基因组鉴定及其表达谱分析[J]. 作物学报, 2026, 52(2): 405-420.
[11] 刘吉昌, 李思烨, 李雪婷, 王洪章, 刘鹏, 张吉旺, 赵斌, 任佰朝, 任昊. 盐胁迫对不同耐盐型夏玉米品种根系生长及养分吸收效率的影响[J]. 作物学报, 2026, 52(2): 565-577.
[12] 孔娜, 刘涛, 刘文婷, 陈刚, 文利超, 邓智超, 郭梅, 李伟, 郭永峰. 烟草NtCEP7基因克隆及其编码小肽在苗期抗旱中的作用分析[J]. 作物学报, 2026, 52(1): 249-261.
[13] 刘海波, 张蕾, 王立琦, 石晓丽, 周文莹, 崔国贤, 佘玮. 苎麻BnGCL1基因响应干旱胁迫的功能研究[J]. 作物学报, 2026, 52(1): 14-27.
[14] 王雅致, 杨飚, 季香林, 石瑛, 张丽莉. 二倍体马铃薯抗旱资源鉴定及抗旱基因初步筛选[J]. 作物学报, 2026, 52(1): 72-84.
[15] 胡城祯, 高维东, 孔斌雪, 王建飞, 车卓, 杨德龙, 陈涛. 小麦TaAPC11基因家族鉴定及TaAPC11-5B参与干旱胁迫的生物学功能研究[J]. 作物学报, 2026, 52(1): 148-164.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!