当前位置:首页 >> 生物学 >>

全基因组关联分析在植物中的应用


亿学与生物互程2013,V01.30
Chemis研&B;oenginee湘9

No.6

doi:10.3969/j.issn.1672—5425.2013.06.002

全基因组关联分析在植物中的应用
涂雨辰1”,田云k2,卢向阳1’2 (1.湖南省农业生物工程研究所,湖南长沙410128;2.湖南农业大学生物科学技术学院,湖南长沙410128)

摘要:全基因组关联分析(GwAS)已被广泛应用到植物遗传学和育种相关的重要性状研究中。关联分析是一种 基于连锁不平衡来识别分子标记之间或候选基因与性状之间关系的方法。简单介绍了GwAS的发展背景、研究原理及 研究策略,对GwAS在重要的植物基因位点方面的应用研究进展进行了综述。 关键词:全基因组关联分析(GwAS);基因位点;重要性状;植物 中图分类号:Q
78

文献标识码:A

文章编号:1672—5425【2013)06一o007一04

供更高的分辨率,在基因水平上利用样本的数量关联


GWAs的发展背景
重要的农艺性状如产量、作物品质以及植物抗病

到表型变化的差异。由于高密度单核苷酸(SNP)多态 性的出现,进行全基因组扫描时,能识别很小范围的与 数量性状变异显著相关的单元型域。这些方法的出现 使得研究性状的可能性不断增大。目前,GWAS已经 确定了许多与性状相关的重要位点,为生产实践提供 理论指导。


性是由多个基因控制的,均受连续变化且脆弱的环境 影响。与单基因控制的性状相比,其遗传基础更为复 杂。Visscher等[13指出,众多科学和生物的发现已通 过全基因组关联分析(GWAS)所验证。GWAS是经 典的定量遗传理论的拓展,对基因的研究具有重大意 义。作为一个整体,定量特征是由许多具有同等作用 的微效基因互相影响,通过建立遗传模型和估计遗传 方差,选择统计参数来迸行研究。许多经典的定量遗 传模型都在育种实践中发挥着重要作用,比如籽粒的 淀粉∞]和维生素A原含量嘲、小麦的籽粒大小和研磨

GWAS的研究原理及研究策略
与日益普及的基因组多态性数据相比,GWAS正

成为研究遗传学定量特征的默认方法。通过GWAS, 已经发现了数以百计的遗传变异与人类疾病相关的复 杂的特征,彻底改变了人类基因图谱,并且现在也广泛 应用于植物研究中。 GwAS的研究原理是:在基因水平上通过分子标 记的手段,对整个基因组内的SNP进行综合分析与分 型,再将不同表现的性状变异统计出来,提出假设,并 且验证其与期望性状间的关联性。 GwAS的研究策略类似于传统的候选基因法。 在设计初始GWAS实验时,一般分成两个阶段:首先 对整个基因组SNP基因分型芯片和SNP多态性进行 统计分析,通常能够筛选少量阳性SNPs,然后在更多 样本中对这些阳性SNPs进行基因分型,最后整合两 个阶段的结果进行分析。GwAS两阶段研究策略减 少了工作量和成本,通过重复实验也减小了误测率;但

品质【4]、拟南芥(Ar口6idD加fs£^口zi彻口)的开花期[51和
抗病性[阳等。然而,在“Micro.effect gene”理论中,具 体影响数量性状的基因尚未发现,分子生物学机制的 定量特征变化也未能得出确切结论。近年随着生物技 术的发展,尤其是分子标记技术的出现与发展,人们对 量化特征的认识从基因水平发展到数量性状分析水 平,充分说明理解遗传机制的定量特征已经上升到分 子水平的高度。 此外,植物的数量性状往往受到等位基因多样性 的影响,而传统研究仅能得到有限的基因组相关解析。 GWAS克服传统基因映射方法的一些局限性,通过提

基金项目:国家自然科学基金资助项目(31000536),湖南省自然科学基金重点资助项目(13JJ2028)
收稿日期:2013一03—12

作者简介:涂雨展(1988一),女,湖南长沙人,硕士研究生,研究方向:植物生物化学研究,E-mail:cherr)rtt.p“ncess@163.com;通讯 作者;耻I云,副研究员,E-mail:tianyun79616@j63.coml卢向阳.教授,E-maill】【iang”ngcn@163.com.

万方数据

Il——————————————————————————一涂雨辰等:全基因组关联分析在植物中的应用/2013奄饕6期
是在第一阶段,所有可能与疾病相关联的SNP多态性 的样本量不足。为了寻找更多的易感基因位点,目前 常用的方法是扩大GWAS的样本量,即在同一时间适 当放宽第一阶段的选择标准并扩大验证SNP多态性 的范围。


1,W以一2,Ⅵ铆一3,W靠一4,R^产1)也表现出显著的相关
性[15。8]。而在大豆中,关联到重要性状的重要位点分 别为I、T、W1、R、O[19’2引,这些位点的确定对植物生长 和作物生产实践都起到重大的作用。
3.3

在挖掘水稻重要遗传基因位点方面的应用 世界上约一半人口以大米为主食,但其产量常常

GWAs在植物中的应用

受到大量病原微生物或昆虫的严重威胁[21。2引。为了有 效识别遗传多样性的种植水稻种质资源,挖掘水稻高 品质性状相关的基因、培养遗传改良品种的水稻,对水 稻实践生产意义重大。 近年来开展的全基因组测序分析可通过基因组之 间的差异挖掘出与水稻重要性状相关的位点。研究者 分别对517类不同的水稻种质测序,鉴定出360万个 SNP序列,构成一个水稻基因高密度单体型图。通过 GWAS研究,解释了水稻品种中36%的表型方差,将 其分为以下5类农艺性状:形态特征(分蘖数、叶角 度)、收益率(谷粒宽度、粒长、千粒重和小穗数)、粮食 质量(糊化温度和直链淀粉含量)、着色(顶端颜色、果 皮颜色和外壳颜色)和生理特性(抽穗期、抗旱能力和 种子破碎的程度)。通过GWAS研究,确定了稻籼亚 种的14个农艺性状。研究者将高通量基因组测序技 术应用于950份来自不同国家的水稻种质材料中,通 过基因分型鉴定基因的变异。该体系能高效、准确地确 定多个重要性状的关联位点,分别鉴定出1个控制花期 的位点和10个与粮食主产性状相关的32个新位点[24]。 3.4在拟南芥性状研究中的应用 研究复杂性状的遗传基础,如开花时间、伸长速度 和产量,已经成为改善作物和了解植物适应性的重点。 拟南芥分布广泛,一直是研究植物自然变异和适应性 的极具价值的模式植物。通过全基因组基因分析, Aranzana等确定了拟南芥基因中控制自然变异、开花 时间和病原体抵抗的相关位点。通过分析拟南芥的 107种表型,GWAS搜寻到250 000个SNP位点口51。 通过GWAS来识别潜在的反应变异,其中最为明显的 反应为控制下胚轴的伸长。结果表明,这些性状的变

3.1在玉米育种方面的应用 在过去的20年中,伴随着基因组学以及分子标记 技术的飞速发展,作物遗传育种方面的研究得到有效 且深入的开展。在玉米研究方面,2011年《Nature G占netics》上发表了1篇关于GWAS研究玉米相关性 状的文章。文中通过GWAS得到一个玉米嵌套的关 联映射面板,确定了玉米重要叶片结构的基础架构特 征和一些关键基因,并证明了叶片性状、基因结构是受 ‘基因多效性及环境相互作用的影响。GWAS结果表 明,突变体基因的变化导致更多的叶片向上生长。玉 米的多样性对GWAS来说是个挑战。然而,在玉米
Hap

Map中160万个SNPs位点被鉴定、嵌套关联绘

图(NAM)也有所发展之后,在玉米研究中有效开展 GWAS不再是难题。玉米NAM嵌板是基于一个参 考线,由25条不同的交叉线组成,然后产生5000个重 组自交系(RILs)。作者运用这种新型的玉米NAM 设计并分析向上叶片角度的结构、叶片长度以及宽度, 再运用Joint
stepwise

regression鉴定出30~36个

QTI。s,解释了74.8%~80.3%的表型变异,以及大于 83%的遗传方差。为了进一步剖析这些QTLs,通过

输入160万个Hap

Map SNPs位点进行GwAS分

析,发现27个NAM母系映射到RILs上。在160万 个测试的SNPs中,作者探测到与叶片角度、长度、宽 度相关的SNP位点分别是203、287、295。在其它复 杂的玉米性状里也能观察到类似的基因结构[7’8],与 动物一样‘91,对自花授粉植物有更大的影响‘协12]。这 些结果表明复杂性状的基因结构在异形杂交和自交的 物种中的进化程度相同。 3.2在大麦及大豆分子标记研究中的应用 Lorenz等[13]研究发现,大麦中单独关联到早熟性 状的位点对应着33个染色体区域,并确定了15个与 其它性状相关的重要标记和7个区域[1“。对于每一 个性状的分析,都需要经过包括SSR、SNP和DArT 在内的760个标记来确定。研究发现,共有130个显 著特征关联标记以及8~22个重要位点,这也代表了 62个标记存在于除染色体5d之外的所有染色体上。

异受基因座的影响,符合Ⅵ圯cA5、Ⅵ厄CA9和
RGAl

3组基因的控制模式。砌CCA5和H尼CA9

都参与生长素生物合成,而RGAl是DELLA家族的 成员。下胚轴表型包括在这项研究中,它们是高度的 高R:FR条件(模拟太阳)、高度的低R:FR条件(模拟 阴影)和2个不同指数响应低R:FR条件。GWAS结 果表明,这些性状的变异由许多基因座的低到中度影 响控制。 不同作物重要性状的相关位点见表1。

而观察图谱的位置,最明显的候选基因(Pp小l,‰一
万方数据

涂雨辰等:全基因组关联分析在植物中的应用/2013年簧6啊————————二————————————————————_I盈
表l
Tab-l
The reIated

不同作物重要性状的相关位点 I∞i of important
traits

for different speci器

种类

性状

位点效量
1 1 l 1 1 1

位点
qSw5
GS3 05cl Rc

文献 [26] [26] [26] [26] [26] [26]

谷粒宽度 谷粒长度
顶端颜色 水稻

果皮颜色 淀粉酶含量 凝胶化温度

waxy
ALK

GwAS discovery[J].Am J 4

Hum Genet,2012,90(1):7—24.
structure

结语
GWAS为大量的基因组测序提供了一个很好的

[2]

wen w。Mei H,Fellg F。et a1.Population
mapping
on

and association

chromosome 7 using



diverse panel of Chine5e germ— Genet.2009,119

plasm

平台,由于其技术方法能扩展到在生态环境中且能适 应遗传变异的结构化背景,与人类疾病GWAS相比, 在植物研究中获得了更大的成功[3引。为了能更高效 地运用GwAS技术,首先,在选择基因样本时,应尽量 选择优质品种的优良基因,以提高分析的精确度与有 效性;其次,SNP位点对作物性状的影响基于影响基 因的表达程度,所以应该提高对基因编码区以及相关 调控方式的认识,更确切地了解植物复杂性状的成因; 最后,GWAS能够一次性对重要性状进行轮廓性概 览,适用于复杂性状的研究。在全基因组层面上,开展 多中心、大样本、反复验证的基因与性状的关联研究, 是全面揭示疾病发生、控制其发展的重要手段。 由于植物性状与基因组之间的关系受多种因素的 作用[3¨,GWAS也面临着一些问题,如结果不能完全 解释某些复杂性状等。但随着功能标记开发、反向遗 传学研究及生理学等学科的综合推进,作为一种综合 性的分析方法,GWAS必将迎来更为广阔的发展空 间。在提高自动化和高效率的同时,通过降低成本、挖 掘与植物重要经济性状或生理性状相关的位点, GwAs必将在植物遗传育种中发挥更大的作用。
参考文献:
[1]
Visscher P

office(∞嘲犯ff懈L.)[J].Theor Appl

(3):459—470.

[3]

Ha巧es C E,Rocheford T R,Bai L,et a1.Natunl genetic variation
in Lyco户P玎P。户sfzo玎fyc如即tapped for maize

biofortifIcation[J].

Sdence.2008,319(5861):330一333.

[4]

Breseghello F,SorreIls and milIing quaIity in

M E.Assodation

mpp啦of

ken"l size

wheat(丁'f£fc口辨4Ps£i豫肌L.)cl|ltivars

口].Genetics,2006,172(2):1165—1177. [5]
Olsen K M,Halldorsdottir S S.Stinchcombe J R,et a1.Linkage
disequilibrium mapping of

Am6fdops缸CRy2 flowering time aI—

leles[J].Genetics,2004,167(3):1361—1369. [6]
Ara眦ana M J,Kim s,zhao K,et aI.Genome.wide association
mapping in Ar口6fdops如identifies previously known noweriIlg

time and pathogen resistance e60.

genes[刀.PLos Genet,2005,1(5):

[7]

Buckler E s,Holland
tecture

J B,Bradbury P了,et a1.The genetic archi—

of maize flowering time[J].Science.2009.325(5941):714—

718.

[8]

Laurie C C,Chasalow S D,LeDeaux
tectufe

J R,et aI.The genetic archi.
concen’

of

respon∞to lon争term artmdal 5e王ection for oil

t抬tionirIthemaize kernel[J].Genetics,2004,168(4)12141—2155.

[9]

FIint

J。Mackay T F C Genetic a”Mtecture of q∞ntitatjve traits

inmice。nies,and

humans[J].GenomeRe5,2009.19(5):723—733.

[10]

Koornneef

M。Alonso-Blanco c,Vreugdenhil D.Naturally oc印r.

ring genetic variation in

Ar口6ido户f如咖4zf4拙口[J].AfIllu

Rev

Plant B.01.2004,55(4):141-17Z.

M。Brown M A.McCarthy M I.et a1.Five y龃rs of

[11]

TakahasK Y.Teshima K

M.Yokoi s,et a1.variations in Hdl

万方数据

口耍——————————————————————————一涂雨辰等:全基因组关联分析在植物中的应用/2013芷第6期
pfoteins,Hd3口promoters,and既d1
ute to

expressbn levels contrib—

[22]

of the lnaga“Y S,Ethe“ngtoll G,Geisler K,et al,InVestigation in rice potential for triterpene synthesis and metabolic 448.

diversity of nowering time in cultivated

rice[J].Proc

NatI

through genome mining

Acad Sd

USA,2009.106(11):4555-4560。

engineering[J]。New

Phytol,20ll,19l(2):432。

[12]

Tumer A,Beales JtFaure s,et a1.The pseudo.r鹳ponse regulato。 Pp壬Hl provides adaptation
to

photoperiod

in

barley[J].Sci—

[23]

Huang X H,Zhao

Y,Wei X H.et aI.Genome-wjde association


ence.2005,315(5750):1031-1034.

study of nowering time and grain yield t珀it5 iIl of 5ingle lection of rice

worldwide col?

[13]Lore北A J.Hamblin M T?Jannink J L.Perfomance
nucIeotide polymDrphisms、rersus haplotypes as50ciation analy5i5 e14079. in

germplasm[J].Nat

Genett2012,44:32—39.

for genome-wide

[24]

Arite T,Umehara

M。Is}1ikawa S.et aI.d14,A strigolactone.in-
an

bafley[J].PLos One,2010,5(11):

5ensitive mutant of rice,shows

accelerated outgl口wth of tilI—

ers[J].Plant&lI
K。zhang G P.A genome.wide association analysis

Physiol,2009?50(8):1416-1424.

[14]

Jin



L.wei

[25]

Atwell S,Huang Y S,ViJ

hj丘lmsson BJ?et

a1.Genome_wide

ass0一

of quantitative

tmit loci for protein fraction

content

in

Tibetan

ciation study of 107 phenotypes in

Arn6idD户5i5髓‘IZf口,^口inbred

稍ld barley口].B.otechnoIogy [15]
Pasam R

Letters,2012,34(1):159—165.

lines口].Nature,2010,465(7298):627—631. [26]
Huang X

K,Shama

R,MaIosetti M,et a1.Genome-、llride assoda—


H,Wei

X H.Sang T,et

a1.Genome-wide association

tion studies for agronomical traits in

world wide spring barley

studies of 14 agronomic traits in 2010,42(11):961—967.

rice

landraces[J].Nat

Genet,

coIlection[J].BMc [16]

Plant BioI,2012.12(1):16.

Varshney R K,Paulo

M J。Grando s,et a1.Genome wide assoda—
tolerance托lated traits in crops

[27]

Li

Q,Yang X.Xu S,et a1.Genome-wide association studies iden—

tion孤alyses

for

drought

barley(Hor-

tmed three independent polymorphisms assodated with口-to—
copherol e36807.
content

dP“研伽zgn旭L)[J].Field
180.

R∞earch?2012,126(1)117l一

in

maize

kernels[J].PLoS

0ne,2012,7(5):

[17]

wang H Y,Smith K P,combs E,et a1.Effect of population size
and unbalanced data
sets
on

[28]

Raman H,Stodart Bt Ryan P R,et a1.Genome.wide assodation analyses of common wheat(了’,if{c“研口盯fft吡m L.)germplasm

QTL detection using genome-“de

a5sociation mapping in barley breeding cai and Applied

gemplasm[J].Theoreti-
[29]

identmes multiple

lod for aluminum

re5i5tance[J].Genome,

Genetics,2012,124(1):11l一124.

2010,53(11):957?966.
Yu L stem

[18]

Wang M H,Jiang N,Jia T,et a1.Genome-wide association map—
ping of agronomic and

X,Morgounov A,Wanyera R,et a1.Identification of U999
rust

morphologic

trait5

in

highly and

stnlctufed
Ge—

resistance lod in winter wheat germplasm using

ge-

populations of barley

cuItivars[J].Theoretical

Applied

nome-wide association (4):749-758.

analysis[J].Theor

Appl Genet,2012,125

netics,ZOl2,124(2):233-246.

[19]

Yang

K.Jeong N,Moon J K。et a1.Genetic analysis of gene5
natural

con-

[30]

Hindorff L A.Sethupathy

P,Junkins H A,et a1.Potential etio?

tromng

va血tion

of

s∞d

coat

and

nower∞lor8

in

soy—

logic and functional implications of genom昏wide association loci for human diseases and

bean[J].J [20]

Hered,2010,10l(6):757?768.

tmits[J].Proc

Natl

Acad sci

usA,

Palmer R G。Pfeiffer T in Soybeans


w。Bu5s

G R.et

a1.Qualitathe Genetic8

2009,106(23):9362—9367.

lmprovement,Production,and Use5(3rd

ed)[M].

[31]

Ku C S,Loy E Y,Pawitan Y.et a1.The pursuit of genome.wide association

Madison(WI):ASA。CSSA,and SSSA,2004:137—214.

studies.where

are

we

now?口].J

Hum Genet,2010,

[21]

clark R

M.Genome-wide assodation studies coming of age in

55(4):195—206.

rice[J].Nat Genet,2010.42(11):926—927.

Application of

Genome_Wide Association Study for

ImpOrtant Traits in Plant
TU YI卜chenl”,TlAN Yunl”,LU Xian矿yan91’2

(1.H“,l口行Agrif“Zz“,。nZ BioP刀gi靠P已ri竹g RPsP口,℃^J刀sziz扯£P,C丘口疗gs^口410128,C,li咒口;
2.CoZZPgP D,BiDsffP挖fP口nd了■f|Il以DZDgy,H“以口,2 Agric“Z£“r口Z己,咒i口P,百f£y,(强口,29s^口410128,(’^i以拉) Abstract:Genome.wide association

study(GWAS)has
an

been widely employed
to

to

inVestigate important traits

in plant genetics and breeding.Association study is markers
or

approach

recognize the relationship between molecular
on

candidate genes with traits in



given population based

linkage disequilibrium.The development

background,research principIe and research strategy of
progress of

GWAS

were introduced in this paper,and the appIication

GWAS

in the research of important plant genetic loci was reviewed.

Keywords:genome_wide association

study(GWAS);genetic

loci;important trait;plant

万方数据

全基因组关联分析在植物中的应用
作者: 作者单位: 刊名: 英文刊名: 年,卷(期): 涂雨辰, 田云, 卢向阳, TU Yu-chen, TIAN Yun, LU Xiang-yang 湖南省农业生物工程研究所,湖南长沙410128;湖南农业大学生物科学技术学院,湖南长沙410128 化学与生物工程 Chemistry & Bioengineering 2013,30(6)

本文链接:http://d.g.wanfangdata.com.cn/Periodical_hbhg201306002.aspx


相关文章:
全基因组关联分析在植物中的应用.pdf
化学与生物工程 2013,Vol.30No.6Chemistr&Bioenineerin ygg 7 :/doi10.3969.issn.1672-5425.2013.06.002j 全基因组关联分析在植物中的应用2 22,, 涂雨辰...
全基因组关联分析在植物中的应用_图文.pdf
全基因组关联分析在植物中的应用 - 亿学与生物互程2013,V01.30Chem
全基因组关联分析在水稻遗传育种中的应用_唐富福.pdf
全基因组关联分析在水稻遗传育种中的应用_唐富福_农学_农林牧渔_专业资料。598...2 水稻全基因组 SNP 的检测 SNP 标记以其数量丰 在植物育种和 遗传研究中, ...
关联分析及其在不同分子标记中的应用.doc
关联分析及其在不同分子标记中的应用摘要:关联分析技术在植物中的应用较晚,是...但是,全基因组关联分析花费巨大,耗时较多,目前完成起来困难重重。 对野生甜菜的...
全基因组关联分析-基于全基因组重测序.pdf
全基因组关联分析-基于全基因组重测序_生物学_自然科学_专业资料。对已有参考...已有参考基因组序列的动植物自然群体,建议样本数≥200个 2. 样本间无明显的亚...
关联分析及其在植物研究中的应用.pdf
关联分析及其在植物研究中的应用 - 河北农业科学, 2010, 14 ( 11
关联分析及其在植物遗传学研究中的应用.pdf
关联分析及其在植物遗传学研究中的应用 - 植物学报 Chinese Bullet
关联分析及其在不同分子标记中的应用综述.doc
由于其良好的实用性,关联分析技术已经成为国际 基因组学的研究热点。介绍了 LD 的定义和度量方法,综述了基于分子标记的关联分析在植物 遗传中的应用,并对其在遗传...
全基因组关联分析在畜禽重要经济性状中的研究进展_论文.pdf
全基因组关联分析在畜禽重要经济性状中的研究进展 - 第 38卷 第 10期 20
全基因组关联分析.ppt
association study简 称GWAS 全基因组关联分析 是指在人类全基因组范 围内找出存在的序列变 异,即单核苷酸多态 性(SNP),从中筛选 出与疾病/性状相关...
关联分析及其在植物研究中的应用_图文.pdf
关联分析及其在植物研究中的应用 - 河北 农业科学,2100,1()541:17
全基因组关联分析_图文.ppt
围内找出存在的序列变 异,即单核苷酸多态 性(SNP),从中筛选 出与疾病/性状...全基因组关联分析在植物... 5页 免费 抑郁症全基因组关联分析... 暂无...
关联分析在作物种质资源分子评价中的应用_图文.pdf
关联分析在作物种质资源分子评价中的应用 - 植物遗传资源学报 2007, 8 ( 3 ) : 366 ~372 Journal of Plant Genetic Resources 关联分析在...
植物关联分析方法的研究进展_图文.pdf
关联分析在人类和动植物遗传研究中的应用日益广泛 , 新方法及其软件包不断涌现。...学者们应用不同的统 计模型和参数与非参数估计方法在全基因组开展了 单标记扫描...
全基因组关联分析方法的研究进展_图文.pdf
全基因组关联分析方法的研究进展_生物学_自然科学_...06015201--][43- 种情况在动植物研究中极为常见 ...A.L.[2] 将主成分分析应用于 GW 首先通过主 ...
代谢组学在植物研究领域中的应用Applicationof(精).pdf
代谢组学在植物研究领域中的应用Applicationof(精) - Botani
全基因组关联分析方法的研究进展_郝兴杰_图文.pdf
全基因组关联分析目前已经成为研究复杂性状和疾病遗传...06015201--][43- 种情况在动植物研究中极为常见 ...A.L.[2] 将主成分分析应用于 GW 首先通过主 ...
关联分析.pdf
关联分析及其在植物遗传学研究中的应用谭贤杰1, 2, 吴子恺1, 程伟东2, ...其次, 在全基因组关联分析数理统计方面, 基因型标记与表 3 展望 自2001年关联...
关联分析及其在植物中的研究进展_论文.pdf
关联分析及其在植物中的研究进展 - 阐述了关联分析的方法、各种分子标记在关联分析中的应用关联分析在植物中的研究进展以及基于第三代分子标记的关联分析的基础单...
关联分析在核心种质中的应用_图文.pdf
关联分析在核心种质中的应用 - 植物遗传资源学报 2007, 8( 3) : 366~ 372 Journa l o f P lant G enetic R esources 关联分析在作物种...
更多相关标签: