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

作物学报 ›› 2021, Vol. 47 ›› Issue (11): 2121-2133.doi: 10.3724/SP.J.1006.2021.04249

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

利用动态转录组学挖掘大豆百粒重候选基因

曾健(), 徐先超, 徐昱斐, 王秀成, 于海燕, 冯贝贝, 邢光南*()   

  1. 南京农业大学大豆研究所 / 国家大豆改良中心 / 农业农村部大豆生物学与遗传育种重点实验室(综合) / 作物遗传与种质创新国家重点实验室/江苏省现代作物生产协同创新中心, 江苏南京 210095
  • 收稿日期:2020-11-20 接受日期:2021-03-19 出版日期:2021-11-12 网络出版日期:2021-04-01
  • 通讯作者: 邢光南
  • 作者简介:E-mail: 2019101129@njau.edu.cn
  • 基金资助:
    国家重点研发计划项目(2016YFD0100201);国家自然科学基金项目(31571694);中央高校基本科研业务费专项资金(KYT201801);长江学者和创新团队发展计划(PCSIRT_17R55);高等学校学科创新引智基地项目111(B08025);国家现代农业产业技术体系(大豆)建设专项(CARS-04);江苏省优势学科建设工程专项和江苏省JCIC-MCP项目

Utilization of dynamic transcriptomics analysis for candidate gene mining of 100-seed weight in soybean

ZENG Jian(), XU Xian-Chao, XU Yu-Fei, WANG Xiu-Cheng, YU Hai-Yan, FENG Bei-Bei, XING Guang-Nan*()   

  1. Soybean Research Institute, Nanjing Agricultural University / National Center for Soybean Improvement / Key Laboratory for Biology and Genetic Improvement of Soybean (General), Ministry of Agriculture and Rural Affairs / National Key Laboratory for Crop Genetic and Germplasm Enhancement / Jiangsu Collaborative Innovation Center for Modern Crop Production, Nanjing 210095, Jiangsu, China
  • Received:2020-11-20 Accepted:2021-03-19 Published:2021-11-12 Published online:2021-04-01
  • Contact: XING Guang-Nan
  • Supported by:
    National Key Research and Development Program of China(2016YFD0100201);National Natural Science Foundation of China(31571694);Fundamental Research Funds for Central Universities(KYT201801);MOE Program for Changjiang Scholars and Innovative Research Team in University(PCSIRT_17R55);Project of Intellectual Base for Discipline Innovation in Colleges and Universities(B08025);China Agriculture Research System(CARS-04);iangsu Higher Education PAPD Program, and the Jiangsu JCIC-MCP Program

RichHTML

36

PDF

1405

摘要：

百粒重是影响大豆产量的重要农艺性状, 揭示其分子基础发掘关键候选基因对大豆改良具有重要意义。本研究通过对12个大豆品种籽粒发育3个时期共36个样本的转录组数据进行加权基因共表达网络分析(weighted gene co-expression network analysis, WGCNA), 得到20个基因共表达模块, 与百粒重及4个粒形性状关联后发现green模块与表型最为相关, 之后根据Gene Significance (GS)值和Eigengene Connectivity (kME)值筛选出13个green模块内的核心基因(hub gene); 然后对2组百粒重存在极显著差异的大豆品种的籽粒发育3个时期分别进行基因差异表达分析发现大豆在籽粒发育前中期可能通过MAPK信号通路调节百粒重大小; 之后对其进行SNPs/InDels挖掘并根据Gene Ontology (GO)注释发现green模块内的Glyma.14G043900和Glyma.15G217400由于SNP变异造成同义以及非同义突变, 且存在调控基因表达相关的GO Terms以及锌指结构域, 表明它们可能通过调控hub基因和差异表达基因调控大豆百粒重及粒形性状。Glyma.15G217400位于已报道的4个百粒重QTL中, 而Glyma.14G043900位于已报道的一个籽粒蛋白含量及一个油脂含量QTL中。通过比对大豆公共数据库发现这2个基因的百粒重增效等位变异受到人工选择, 其频率从野生大豆到地方品种再到育成品种的过程中逐渐升高。这些结果为进一步发掘大豆百粒重候选基因及其表达调控机制提供了新思路。

关键词: 大豆, 百粒重, 转录组学, WGCNA, 候选基因

Abstract:

100-seed weight of soybean is an important agronomic trait that affects yield, and it is of great significance to reveal its molecular basis and discover key candidate genes for soybean improvement breeding. In this study, weighted gene co-expression network analysis (WGCNA) was performed on the transcriptome data of 36 samples from 12 soybean varieties at three stages of seed development, and 20 gene co-expression modules were obtained. After correlating with 100-seed weight and four-seed shape traits, the green module was found to be most correlated with the phenotypes. Then 13 hub genes of green module were screened based on the Gene Significance (GS) and Eigengene Connectivity (kME) value. Gene differential expression of two groups of soybean varieties with extremely significant differences in 100-seed weight showed that the MAPK signaling pathway in the early and mid-term of seed development might regulate the 100-seed weight in soybean. According to SNPs/InDels calling and Gene Ontology (GO) annotation, Glyma.14G043900 and Glyma.15G217400 in the green module caused synonymous and non-synonymous coding mutations due to SNP mutations, and there were GO Terms and zinc finger domains related to gene expression regulation. These results suggested that they might regulate the 100-seed weight and seed shape of soybeans by regulating the hub gene and differentially expressed genes. Furthermore, Glyma.15G217400 was located in four reported QTLs of 100-seed weight, while Glyma.14G043900 was located in a reported seed protein content QTL and an oil content QTL. Compared with soybean public database, the increasing 100-seed weight alleles of the two genes were artificially selected and their frequency was gradually increased from wild accessions to landraces, resulting in the improved cultivars. These results provide new ideas for further discovering 100-seed weight candidate gene in soybean and its expression regulation mechanism.

Key words: soybean, 100-seed weight, transcriptomics, WGCNA, candidate gene

图1

软阈值β的确定 A图纵坐标是无尺度网络模型指数, 红线代表R2=0.85。B图纵坐标每一个软阈值对应的平均连接度。图A和图B的横坐标均代表软阈值β。"

表1

SNPs/InDels过滤标准"

SNPs过滤标准
SNPs filter standards		 InDels过滤标准
InDels filter standards
QualByDepth (QD) > 2.0 QualByDepth (QD) > 2.0
RMSMappingQuality (MQ) > 40.0 —
FisherStrand (FS) < 60.0 FisherStrand (FS) < 200.0
StrandOddsRatio (SOR) < 3.0 StrandOddsRatio (SOR) < 10.0
MappingQualityRankSumTest (MQRankSum) > -12.5 MappingQualityRankSumTest (MQRankSum) > -12.5
ReadPosRankSumTest (ReadPosRankSum) > -8.0 ReadPosRankSumTest (ReadPosRankSum) > -8.0

图2

12个大豆品种百粒重的箱型图 横坐标代表12个大豆品种, 纵坐标表示百粒重。一个点代表一次观察值。品种间不同大写字母说明百粒重基于Duncan氏新复极差法在P = 0.01水平的差异显著性。"

图3

模块内基因数目柱状图 横轴代表模块, 纵轴代表模块内基因数目。"

图4

模块与性状关联热图(A)及green模块内2904个基因的表达量热图(B) A图中横轴表示百粒重及粒形性状, 纵轴表示每一个模块的特征向量; 红色的格子代表性状与模块具有正相关性, 蓝色的格子代表性状与模块具有负相关性。B: green模块内2904个基因在12个大豆品种3个时期表达量行标准化后绘制。G、H和I分别为籽粒发育前期、中期和后期。"

图5

green模块内2904个基因的GO富集 横坐标(每列)代表富集基因占输入基因比例, 纵坐标(每行)代表GO Terms。点的大小代表富集的基因个数, 点的颜色代表经费舍尔精确测验矫正后的Padjust值。"

图6

百粒重存在极显著差异的2组大豆品种籽粒发育不同时期基因的表达差异(A)及小粒品种下调基因的MAPK信号通路富集图(B) G、H和I分别代表籽粒发育前期、中期和后期。Down-regulated genes代表小粒品种中下调的基因, Up-regulated genes代表小粒品种中上调的基因。"

图7

green模块内hub基因GS值和kME值的柱形图(A)及不同百粒重大豆品种组间差异表达基因的表达量聚类热图(B) 100-SW: 百粒重; SA: 籽粒面积; SL: 籽粒长度; SP: 籽粒周长; SW: 籽粒宽度。G、H和I分别代表籽粒发育前期、中期和后期。Large和Small分别代表P = 0.01差异极显著的大粒大豆和小粒大豆品种组。图B中加粗的基因代表第1种表达模式的基因。"

图8

SNPs/InDels及DEGs在大豆染色体的分布图 A: SNPs密度; B: InDels密度; C: log2 (Fold Change) of DEGs, 内圈为在小粒品种组中显著上调基因, 外圈为下调基因。"

表2

SNPs/InDels分布区域"

变异位置
Variation region		 SNPs InDels
数目
Count		 比例
Percentage (%)		 数目
Count		 比例
Percentage (%)
DOWNSTREAM 41,636 23.67 5788 26.22
EXON 72,950 41.48 3028 13.72
INTERGENIC 2100 1.19 424 1.92
INTRON 8576 4.88 1625 7.36
SPLICE_SITE_ACCEPTOR 31 0.02 11 0.05
SPLICE_SITE_DONOR 31 0.02 12 0.05
SPLICE_SITE_REGION 334 0.19 69 0.31
TRANSCRIPT — — 1 0.01
UPSTREAM 29,780 16.93 4128 18.70
UTR_3_PRIME 14,683 8.35 4400 19.94
UTR_5_PRIME 5749 3.27 2586 11.72

表3

Green模块内在大粒和小粒品种组间存在非同义SNP突变的10个基因"

基因
Gene		 染色体
Chromosome		 位置
Position		 小粒品种组等位变异
Allele of small 100-SW		 大粒品种组等位变异
Allele of large 100-SW
Glyma.02G139500 Gm02 14480344 G A
Glyma.02G140100 Gm02 14531600 T C
Glyma.06G179600 Gm06 15246959 A G
Glyma.07G018200 Gm07 1456246 A G
Glyma.12G057100 Gm12 4156207 T C
Glyma.14G043900 Gm14 3335050 T C
Glyma.15G217400 Gm15 35930038 A T
Glyma.17G071400 Gm17 5580413 A G
Glyma.17G076400 Gm17 5974992 C T
Glyma.18G237000 Gm18 52573168 T C

表4

候选基因Glyma.14G043900和Glyma.15G217400的GO注释"

基因
Gene		 GO编号
GO ID		 GO类型
GO type		 注释描述
Annotation description
Glyma.14G043900 GO:0006355 生物过程
Biological process		 以DNA为模板的转录调控
Regulation of transcription, DNA-templated
GO:0008270 分子功能
Molecular function		 锌离子结合
Zinc ion binding
Glyma.15G217400 GO:0061158 生物过程
Biological process		 3'端非翻译区介导的mRNA失稳
3'-UTR-mediated mRNA destabilization
GO:0005829 细胞组分
Cellular component		 胞质溶胶
Cytosol
GO:0003730 分子功能
Molecular function		 mRNA 3'端非翻译区结合
mRNA 3'-UTR binding
GO:0046872 分子功能
Molecular function		 金属离子结合
Metal ion binding

图9

green模块内候选基因、hub基因以及差异表达基因的共表达网络 红色点代表候选基因, 黄色点代表hub基因, 绿色点代表DEGs, 蓝色点代表既是DEGs同时又是hub基因。线的颜色深浅代表基因间的weight值大小。黑色的线代表基因间的weight大于0.1, 灰色线代表基因间的weight小于0.1但大于0.01。"

图10

不同大豆类型中Glyma.14G043900和Glyma.15G217400非同义SNP突变等位变异的频率变化 圆点代表百粒重的增效等位变异, 而三角形代表百粒重的减效等位变异。"

[1] Lu X, Xiong Q, Cheng T, Li Q T, Liu X L, Bi Y D, Li W, Zhang W K, Ma B, Lai Y C, Du W G, Man W Q, Chen S Y, Zhang J S. A PP2C-1 allele underlying a quantitative trait locus enhances soybean 100-seed weight. Mol Plant, 2017, 10: 670-684.
doi: 10.1016/j.molp.2017.03.006
[2] Li J, Zhao J, Li Y, Gao Y, Hua S, Nadeem M, Sun G, Zhang W, Hou J, Wang X, Qiu L. Identification of a novel seed size associated locus SW9-1 in soybean. Crop J, 2019, 7: 548-559.
doi: 10.1016/j.cj.2018.12.010
[3] 蒋费涛, 王书平, 祁俊生, 熊春霞. 转录组学技术及其在植物系统学上的研究进展. 现代盐化工, 2020, 47(4):14-17.
Jiang F T, Wang S P, Qi J S, Xiong C X. Research progress of transcriptional technology and its advances in plant phylogeny. Mod Salt Chem Ind, 2020, 47(4):14-17 (in Chinese with English abstract).
[4] Rory S, Marta G, James H. RNA sequencing: the teenage years. Nat Rev Genet, 2019, 20: 631-656.
doi: 10.1038/s41576-019-0150-2
[5] 杨宇昕, 桑志勤, 许诚, 代文双, 邹枨. 利用WGCNA进行玉米花期基因共表达模块鉴定. 作物学报, 2019, 45: 161-174.
Yang Y X, Sang Z Q, Xu C, Dai W S, Zou C. Identification of co-expression modules for maize flowering stage genes using WGCNA. Acta Agron Sin, 2019, 45: 161-174 (in Chinese with English abstract).
[6] Hollender C A, Kang C, Darwish O, Geretz A, Matthews B F, Slovin J, Alkharouf N, Liu Z. Floral transcriptomes in woodland strawberry uncover developing receptacle and anther gene networks. Plant Physiol, 2014, 165: 1062-1075.
pmid: 24828307
[7] Greenham K, Guadagno C R, Gehan M A, Mockler T C, Weinig C, Ewers B E, McClung C R. Temporal network analysis identifies early physiological and transcriptomic indicators of mild drought in Brassica rapa. eLife, 2017, 6: e29655.
doi: 10.7554/eLife.29655
[8] Yang S N, Miao L, He J B, Zhang K, Li Y, Gai J Y. Dynamic transcriptome changes related to oil accumulation in developing soybean seeds. Int J Mol Sci, 2019, 20: 2202.
doi: 10.3390/ijms20092202
[9] Lu X, Li Q T, Xiong Q, Li W, Bi Y D, Lai Y C, Liu X L, Man W Q, Zhang W K, Ma B, Chen S Y, Zhang J S. The transcriptomic signature of developing soybean seeds reveals the genetic basis of seed trait adaptation during domestication. Plant J, 2016, 86: 530-544.
doi: 10.1111/tpj.2016.86.issue-6
[10] Lopez-Maestre H, Brinza L, Marchet C, Kielbassa J, Bastien S, Boutigny M, Monnin D, Filali A E, Carareto C M, Vieira C, Picard F, Kremer N, Vavre F, Sagot M F, Lacroix V. SNP calling from RNA-seq data without a reference genome: identification, quantification, differential analysis and impact on the protein sequence. Nucleic Acids Res, 2016, 44: e148.
[11] Rogier O, Chateigner A, Amanzougarene S, Lesage-Descauses M C, Balzergue S, Brunaud V, Caius J, Soubigou-Taconnat L, Jorge V, Segura V. Accuracy of RNAseq based SNP discovery and genotyping in Populus nigra. BMC Genomics, 2018, 19: 909.
doi: 10.1186/s12864-018-5239-z pmid: 30541448
[12] Liu Y C, Du H L, Li P C, Shen Y T, Peng H, Liu S L, Zhou G A, Zhang H K, Liu Z, Shi M, Huang X H, Li Y, Zhang M, Wang Z, Zhu B G, Han B, Liang C Z, Tian Z X. Pan-genome of wild and cultivated soybeans. Cell, 2020, 182: 162-176.
doi: 10.1016/j.cell.2020.05.023
[13] 丁琦, 徐伟, 李蒙, 王秀成, 卢伟, 盖钧镒, 王玲, 邢光南. 基于分水岭和统计矩的大豆籽粒形态参数测量方法. 大豆科学, 2019, 38: 960-967.
Ding Q, Xu W, Li M, Wang X C, Lu W, Gai J Y, Wang L, Xing G N. Measurement method of soybean seed morphological parameters based on watershed and statistical moment. Soybean Sci, 2019, 38: 960-967 (in Chinese with English abstract).
[14] Pertea M, Kim D, Pertea G M, Leek J T, Salzberg S L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat Protoc, 2016, 11: 1650-1667.
doi: 10.1038/nprot.2016.095
[15] Michael I L, Wolfgang H, Simon A. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol, 2014, 15: 550.
pmid: 25516281
[16] Peter L, Steve H. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics, 2008, 9: 1-32.
doi: 10.1186/1471-2105-9-1
[17] Yu G C, Wang L G, Han Y Y, He Q Y. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics, 2012, 16: 284-287.
doi: 10.1089/omi.2011.0118
[18] Paul S, Andrew M, Owen O, Nitin S B, Jonathan T W, Daniel R, Nada A, Benno S, Trey I. Cytoscape: a software environment for integrated. models of biomolecular interaction networks. Genome Res, 2003, 13: 2498-2504.
doi: 10.1101/gr.1239303
[19] Engström P G, Steijger T, Sipos B, Grant G R, Kahles A, Rätsch G, Goldman N, Hubbard T J, Harrow J, Guigó R, Bertone P. Systematic evaluation of spliced alignment programs for RNA-seq data. Nat Methods, 2013, 10: 1185-1191.
doi: 10.1038/nmeth.2722 pmid: 24185836
[20] Alexander D, Carrie A D, Felix S, Jorg D, Chris Z, Sonali J, Philippe B, Mark C, Thomas R G. STAR: ultrafast universal RNA-seq aligner. Bioinformatics, 2013, 29: 15-21.
doi: 10.1093/bioinformatics/bts635 pmid: 23104886
[21] Li H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics, 2011, 27: 2987-2993.
doi: 10.1093/bioinformatics/btr509
[22] Heldenbrand J R, Baheti S, Bockol M A, Drucker T M, Hart S N, Hudson M E, Iyer R K, Kalmbach M T, Kendig K I, Klee E W, Mattson N R, Wieben E D, Wiepert M, Wildman D E, Mainzer L S. Recommendations for performance optimizations when using GATK3.8 and GATK4. BMC Bioinf, 2019, 20: 722.
doi: 10.1186/s12859-019-3277-4
[23] Petr D, Adam A, Goncalo A, Cornelis A A, Eric B, Mark A D, Robert E H, Gerton L, Gabor T M, Stephen T S, Gilean M, Richard D. The variant call format and VCFtools. Bioinformatics, 2011, 27: 2156-2158.
doi: 10.1093/bioinformatics/btr330
[24] Cingolani P, Platts A, Wang L L, Coon M, Nguyen T, Wang L, Land S J, Lu X Y, Ruden D M. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff. Fly, 2012, 6: 80-92.
doi: 10.4161/fly.19695 pmid: 22728672
[25] Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones S J, Marra M A. Circos: an information aesthetic for comparative genomics. Genome Res, 2009, 19: 1639-1645.
doi: 10.1101/gr.092759.109 pmid: 19541911
[26] Torkamaneh D, Laroche J, Valliyodan B, O’Donoughue L, Belzile F. Soybean (Glycine max) Haplotype Map (GmHapMap): a universal resource for soybean translational and functional genomics. Plant Biotechnol J, 2020, 19: 324-334.
doi: 10.1111/pbi.v19.2
[27] Li N, Xu R, Li Y. Molecular networks of seed size control in plants. Annu Rev Plant Biol, 2019, 70: 435-463.
doi: 10.1146/annurev-arplant-050718-095851
[28] Cheng Z J, Zhao X Y, Shao X X, Wang F, Zhou C, Liu Y G, Zhang Y, Zhang X S. Abscisic acid regulates early seed development in Arabidopsis by ABI5-mediated transcription of SHORT HYPOCOTYL UNDER BLUE1. Plant Cell, 2014, 26: 1053-1068.
doi: 10.1105/tpc.113.121566
[29] Koichiro A, Tokunori H, Kanna S I, Miyako U T, Hidemi K, Makoto M. A novel AP2-type transcription factor, SMALL ORGAN SIZE1, controls organ size downstream of an auxin signaling pathway. Plant Cell Physiol, 2014, 55: 897-912.
doi: 10.1093/pcp/pcu023 pmid: 24486766
[30] Gu Y Z, Li W, Jiang H W, Wang Y, Gao H H, Liu M, Chen Q S, Lai Y C, He C Y. Differential expression of a WRKY gene between wild and cultivated soybeans correlates to seed size. J Exp Bot, 2017, 68: 2717-2729.
doi: 10.1093/jxb/erx147
[31] Li N, Li Y H. Signaling pathways of seed size control in plants. Curr Opin Plant Biol, 2016, 33: 23-32.
doi: 10.1016/j.pbi.2016.05.008
[32] Hans W, Ljudmilla B, Ulrich W. Molecular physiology of legume seed development. Annu Rev Plant Biol, 2005, 56: 253-279.
doi: 10.1146/annurev.arplant.56.032604.144201
[33] 倪资园, 蔡丽, 谢皓, 陈学珍. 大豆籽粒发育过程中营养物质的积累动态分析. 北京农学院学报, 2011, 26(2):1-3.
Ni Z Y, Cai L, Xie H, Chen X Z. Analysis of dynamic accumulation on nutrients during maturing process in soybean seed. J Beijing Agric Coll, 2011, 26(2):1-3 (in Chinese with English abstract).
[34] Xu R, Duan P G, Yu H Y, Zhou Z K, Zhang B L, Wang R C, Li J, Zhang G Z, Zhuang S S, Lyu J, Li N, Chai T Y, Tian Z X, Yao S G, Li Y H. Control of grain size and weight by the OsMKKK10-OsMKK4-OsMAPK6 signaling pathway in rice. Mol Plant, 2018, 11: 860-873.
doi: 10.1016/j.molp.2018.04.004
[35] Srivastava A K, Lu Y M, Zinta G, Lang Z B, Zhu J K. UTR-dependent control of gene expression in plants. Trends Plant Sci, 2018, 23: 248-259.
doi: S1360-1385(17)30255-8 pmid: 29223924
[36] Mian M A, Bailey M A, Tamulonis J P, Shipe E R, Carter T E, Parrott W A, Ashley D A, Hussey R S, Boerma H R. Molecular markers associated with seed weight in two soybean populations. Theor Appl Genet, 1996, 93: 1011-1016.
doi: 10.1007/BF00230118 pmid: 24162474
[37] Orf J H, Chase K, Jarvik T, Mansur L M, Cregan P B, Adler F R, Lark K G. Genetics of soybean agronomic traits: I. Comparison of three related recombinant inbred populations. Crop Sci, 1999, 39: 1642-1651.
doi: 10.2135/cropsci1999.3961642x
[38] Li W X, Zheng D H. QTL mapping for major agronomic traits across two years in soybean (Glycine max L. Merr.). J Crop Sci Biotechnol, 2008, 11: 171-176.
[39] Teng W, Han Y, Du Y, Sun D, Zhang Z, Qiu L, Sun G, Li W. QTL analyses of seed weight during the development of soybean (Glycine max L. Merr.). Heredity, 2009, 102: 372-380.
doi: 10.1038/hdy.2008.108 pmid: 18971958
[40] Akond M, Liu S M, Boney M, Kantartzi St K, Meksem K, Bellaloui N, Lightfoot D A, Kassem M A. Identification of quantitative trait loci (QTL) underlying protein, oil, and five major fatty acids’ contents in soybean. Am J Plant Sci, 2014, 5: 158-167.
doi: 10.4236/ajps.2014.51021
[1] 习千辉, 徐梓瑗, 刘梦梦, 王宏艺, 郎凯琳, 井震海, 陈锋, 赵磊. 小麦籽粒铜含量的全基因组关联分析及候选基因预测[J]. 作物学报, 2026, 52(6): 1604-1617.
[2] 王亚, 赵宜婷, 王宙, 杨俊芳, 张宏斌, 曹越. 转录组-代谢组联合分析蓖麻蜡质合成相关基因[J]. 作物学报, 2026, 52(6): 1774-1787.
[3] 金昱何, 王雪菲, 徐张一娃, 缪怡宁, 蒋云杰, 伊莹, 缪德麟, 朱静仪, 钟一帆, 陈铭亨, 方芳, 刘鹏. 外源激素对低温胁迫下大豆叶片叶绿素荧光参数及抗氧化酶系统的影响[J]. 作物学报, 2026, 52(6): 1817-1829.
[4] 唐宽强, 李公允, 宋美毅, 赵雪, 常春玲. 大豆株高性状全基因组关联分析及预测模型构建[J]. 作物学报, 2026, 52(6): 1743-1756.
[5] 姚术, 郭凯悦, 翟慧慧, 姚佳慧, 邓文琪, 闫玲, 黄驰, 高阳, 俞嫣然, 赵振邦, 李英慧, 王晓波, 李佳佳. 大豆苗期耐低铁综合评价及优异种质筛选[J]. 作物学报, 2026, 52(5): 1373-1387.
[6] 杨欣雨, 崔文涛, 迪力尼格尔·阿力木, 汪凯翔, 吴鹏昊, 任姣姣. 玉米穗上叶片数全基因组关联分析和全基因组选择[J]. 作物学报, 2026, 52(5): 1573-1590.
[7] 马亮, 马璐, 张舒钰, 章慧敏, 王仁明, 宋旭东, 张振良, 冒宇翔, 陆虎华, 陈国清, 郝德荣, 周广飞. 玉米苞叶数目转录组分析及候选基因鉴定[J]. 作物学报, 2026, 52(3): 790-801.
[8] 杨宗桃, 杨婷, 王禹童, 艾静, 李燕烨, 刘家勇, 邓军, 赵勇, 张跃彬. 甘蔗CLC基因家族鉴定与表达分析[J]. 作物学报, 2026, 52(3): 722-734.
[9] 张晴, 杨昱, 郭茜, 岳霈尧, 殷丛丛, 牛景萍, 赵晋忠, 杜维俊, 岳爱琴. 大豆GmARA6a的克隆及响应盐胁迫的功能分析[J]. 作物学报, 2026, 52(2): 480-493.
[10] 鲁雅妮, 丁超杰, 张煜, 杜习军, 齐学礼, 胡琳, 许为钢. 河南省200份小麦品种苗期茎基腐病抗性鉴定与全基因组关联分析[J]. 作物学报, 2026, 52(2): 363-375.
[11] 田甲春, 葛霞, 李守强, 李梅, 田世龙, 张亚倩, 程建新, 李玉梅. 低O2高CO2贮藏环境延缓马铃薯块茎衰老的作用机制[J]. 作物学报, 2026, 52(1): 262-278.
[12] 王雅致, 杨飚, 季香林, 石瑛, 张丽莉. 二倍体马铃薯抗旱资源鉴定及抗旱基因初步筛选[J]. 作物学报, 2026, 52(1): 72-84.
[13] 董丽华, 董成艳, 李正楠, 余静, 叶靓, 刘芳, 谭静. 玉米禾谷镰孢穗腐病抗性候选基因的筛选与鉴定[J]. 作物学报, 2026, 52(1): 131-147.
[14] 李璐琪, 程宇坤, 白斌, 雷斌, 耿洪伟. 小麦叶片气孔相关性状全基因组关联分析[J]. 作物学报, 2025, 51(9): 2266-2284.
[15] 吉白璐, 孙艺文, 刘万峰, 钱亚新, 蒋彩虹, 耿锐梅, 刘旦, 程立锐, 杨爱国, 黄立钰, 李晓旭, 蒲文宣, 高军平, 张强, 文柳璎. 烟草脂类合成关键基因NtLPAT的功能验证[J]. 作物学报, 2025, 51(9): 2527-2537.
Viewed
Full text


Abstract

Cited
  Shared   
  Discussed   
No Suggested Reading articles found!
京ICP备15008789号-2