東鄉(xiāng)野生稻與日本晴多態(tài)性標(biāo)記的開發(fā)

馬小定 唐江紅 張佳妮 崔 迪 李 慧 黎毛毛,* 韓龍植,*

?

東鄉(xiāng)野生稻與日本晴多態(tài)性標(biāo)記的開發(fā)

馬小定1唐江紅2張佳妮2崔 迪1李 慧3黎毛毛3,*韓龍植1,*

1中國農(nóng)業(yè)科學(xué)院作物科學(xué)研究所 / 農(nóng)作物基因資源與基因改良國家重大科學(xué)工程, 北京 100081;2重慶師范大學(xué), 重慶 401331;3江西省農(nóng)業(yè)科學(xué)院水稻研究所, 江西南昌 330200

東鄉(xiāng)野生稻是分布于全球最北端的一種普通野生稻, 與亞洲栽培稻基因組差異較大, 目前缺少覆蓋其全基因組的分子標(biāo)記。本文以東鄉(xiāng)野生稻和日本晴為材料, 通過篩選已有的1017個標(biāo)記, 并利用東鄉(xiāng)野生稻基因組重測序信息設(shè)計的217個InDel標(biāo)記, 共檢測出203個標(biāo)記在東鄉(xiāng)野生稻與日本晴間呈現(xiàn)多態(tài)性。這些標(biāo)記均勻分布于12條染色體, 平均間隔1.9 Mb, 基本覆蓋東鄉(xiāng)野生稻全基因組區(qū)域。通過對秈粳亞種的檢驗分析, 發(fā)現(xiàn)該套多態(tài)性分子標(biāo)記在東鄉(xiāng)野生稻與粳稻雜交后代群體基因型分析上具有較高的應(yīng)用價值。本研究結(jié)果為發(fā)掘東鄉(xiāng)野生稻的有利基因以及分子標(biāo)記輔助選擇育種提供了有力的工具。

InDel標(biāo)記; 基因組重測序; 粳稻

東鄉(xiāng)野生稻于1978年在我國江西省東鄉(xiāng)縣東源鄉(xiāng)被發(fā)現(xiàn), 當(dāng)?shù)胤Q之為“野禾”, 后經(jīng)水稻和植物學(xué)相關(guān)專家鑒定, 確定為普通野生稻[1-2], 而且是迄今為止發(fā)現(xiàn)的分布于全球最北端(28°14￠N)的普通野生稻。東鄉(xiāng)野生稻具有豐富的遺傳多樣性, 而且攜帶耐寒、耐旱、耐貧瘠、廣親和性、野敗育性恢復(fù)性、胞質(zhì)雄性不育、抗病蟲和高產(chǎn)等相關(guān)基因, 對水稻基礎(chǔ)研究和產(chǎn)業(yè)發(fā)展有著舉足輕重的作用, 被譽(yù)為“野生植物類的大熊貓”[3]。

分子標(biāo)記是繼形態(tài)標(biāo)記、細(xì)胞標(biāo)記和生化標(biāo)記后發(fā)展起來的一種簡單易用的遺傳標(biāo)記?，F(xiàn)在水稻研究中被廣泛使用的分子標(biāo)記有簡單序列重復(fù)(Simple Sequence Repeat, SSR)標(biāo)記、插入/缺失(Insertion-Deletion, InDel)標(biāo)記和單核苷酸多態(tài)性(Single Nucleotide Polymorphism, SNP)標(biāo)記。SSR標(biāo)記是根據(jù)廣泛存在于基因組序列中的、由1~6個核苷酸為重復(fù)單位組成的串聯(lián)重復(fù)序列而設(shè)計的分子標(biāo)記。由于其分布廣泛、多態(tài)性豐富, 在不同作物功能基因組研究中發(fā)揮了重要的作用。InDel標(biāo)記是指根據(jù)相同位點(diǎn)序列由不同數(shù)目的核苷酸插入或缺失形成的差異而設(shè)計的多態(tài)性標(biāo)記[4]。InDel標(biāo)記具有基因組差異位點(diǎn)豐富(例如, 日本晴和93-11基因組序列中平均每953 bp就存在1個InDel[5])、染色體位置明確、條帶清晰、結(jié)果可靠等優(yōu)點(diǎn)[6-7]。SNP標(biāo)記是在基因組水平上單個核苷酸的變異引起的DNA序列多態(tài)性, 形式包括單堿基的缺失、插入、轉(zhuǎn)換及顛換等。SNP標(biāo)記遺傳穩(wěn)定性高, 位點(diǎn)豐富且分布廣泛, 檢測快速并且容易實(shí)現(xiàn)自動化分析[8-10]。但是目前由于檢測方法及成本的限制, 還未被廣泛采用。除了上述分子標(biāo)記外, 近年來由于測序技術(shù)的發(fā)展和成本的降低, 一些研究機(jī)構(gòu)開始使用全基因組測序技術(shù)進(jìn)行起源進(jìn)化和功能基因組學(xué)的相關(guān)研究。全基因組重測序是對已知基因組序列的物種進(jìn)行不同個體的基因測序, 并在此基礎(chǔ)上對個體或者群體進(jìn)行差異性分析。通過全基因組重測序分析, 可以找到大量的SNP、InDel、結(jié)構(gòu)變異(Structure Variation, SV)和拷貝數(shù)變異(Copy Number Variation, CNV)等位點(diǎn)[11]。

前期, 雖然已有一些研究者利用東鄉(xiāng)野生稻為供體構(gòu)建了染色體片段置換系, 但是使用的主要是SSR標(biāo)記, 而且基因組覆蓋不完全[12]。主要是由于東鄉(xiāng)野生稻和亞洲栽培稻基因組存在較大的差異, 多態(tài)性分子標(biāo)記數(shù)量和質(zhì)量不佳。本研究以東鄉(xiāng)野生稻和日本晴為材料, 通過篩選現(xiàn)有分子標(biāo)記, 以及利用東鄉(xiāng)野生稻基因組重測序數(shù)據(jù)設(shè)計的InDel標(biāo)記, 形成一套覆蓋全基因組的、條帶清晰、結(jié)果可靠的分子標(biāo)記, 為東鄉(xiāng)野生稻的基礎(chǔ)研究和育種利用提供有力的遺傳工具。

1 材料與方法

1.1 研究材料

東鄉(xiāng)野生稻取自江西省農(nóng)業(yè)科學(xué)院水稻研究所保存的東鄉(xiāng)野生稻庵家山居群中的C35單株, 日本晴來源于中國農(nóng)業(yè)科學(xué)院作物科學(xué)研究所國家種質(zhì)庫。按照常規(guī)大田種植方法進(jìn)行東鄉(xiāng)野生稻的種植和栽培稻日本晴育苗移栽。

1.2 基因組重測序

分別從3個東鄉(xiāng)野生稻植株上取幼嫩的葉片, 然后立即在液氮中冷凍并儲存在–80℃冰箱備用。用DNeasy plant mini kit (QIAGEN)試劑盒提取基因組總DNA, 經(jīng)NanoDrop、Qubit和瓊脂糖凝膠檢測DNA提取的質(zhì)量和濃度。5 μg DNA用于北京諾禾致源科技股份有限公司(北京)基因組重測序(Illumina HiSeq2500)。測序下機(jī)的原始數(shù)據(jù)經(jīng)過濾, 去除接頭序列和低質(zhì)量測序數(shù)據(jù)獲得有效數(shù)據(jù)。有效測序數(shù)據(jù)通過BWA軟件[13]比對到參考基因組(ftp://ftp.ensemblgenomes.org/pub/plants/release-28/fasta/ oryza_sativa/dna/), SAMTOOLs軟件[14]檢測長度小于50 bp的小片段InDel, ANNOVAR軟件[15]對檢測出的InDel注釋。

1.3 InDel標(biāo)記的設(shè)計

根據(jù)東鄉(xiāng)野生稻基因組重測序結(jié)果, 對比日本晴基因組序列信息(IRGSP-1.0.28版), 以每隔1.5~3.0 Mb的距離在各染色體上查找對應(yīng)的插入缺失位點(diǎn), 在位點(diǎn)前后各選擇120 bp序列, 利用Primer Premier 5.0設(shè)計引物。引物長度一般在19~30 bp之間,m值控制在(60±1)℃, 產(chǎn)物設(shè)計長度110~180 bp。

1.4 DNA提取、PCR擴(kuò)增及產(chǎn)物檢測

分別從東鄉(xiāng)野生稻和日本晴植株上取新生的嫩葉, 采用常規(guī)CTAB法提取基因組DNA[16]。

PCR體系10 μL包含模板DNA 0.5 μL、10×PCR buffer 1.0 μL、10 mmol L–1dNTPs 0.3 μL、10 μmol L–1正反引物各0.5 μL、DNA聚合酶0.5 U, 用ddH2O補(bǔ)齊至10 μL。PCR程序為94℃ 預(yù)變性5 min; 94℃ 30 s, 59℃ 30 s, 72℃ 30 s, 35個循環(huán); 72℃延伸5 min; 16℃保存。

PCR產(chǎn)物加2 μL 6×Loading buffer, 混勻后取3.5 μL于8%聚丙烯酰胺凝膠電泳槽中, 180 V穩(wěn)壓電壓10 min后調(diào)電壓至250 V, 再電泳大約70 min, 電泳結(jié)束后用銀染法染膠顯色[17]。

2 結(jié)果與分析

2.1 多態(tài)性分子標(biāo)記開發(fā)

針對東鄉(xiāng)野生稻、日本晴和東鄉(xiāng)野生稻/日本晴雜合體(將兩者DNA等量混合, 模擬雜合體基因組DNA)樣本, 利用實(shí)驗室現(xiàn)有的661對SSR標(biāo)記引物和356對InDel標(biāo)記引物進(jìn)行多態(tài)性檢測, 獲得有多態(tài)性的標(biāo)記81 個, 包括42個SSR標(biāo)記和39個InDel標(biāo)記, 多態(tài)率為8%。

為了獲得更多的多態(tài)性標(biāo)記, 對東鄉(xiāng)野生稻進(jìn)行了基因組重測序, 通過與對照基因組比對, 分別在染色體基因上游、外顯子、內(nèi)含子和基因下游區(qū)共檢測到373,501個InDel。根據(jù)此結(jié)果, 設(shè)計了217對InDel標(biāo)記引物, 經(jīng)過檢測, 其中122個標(biāo)記顯示出多態(tài)性, 多態(tài)率為56% (圖1)。綜合以前和本次設(shè)計的引物, 共獲得203個東鄉(xiāng)野生稻與日本晴間有多態(tài)的分子標(biāo)記。每條染色體上分子標(biāo)記數(shù)目為9~30個, 平均間隔為1.2~2.6 Mb (表1和附表1)。

圖1 多態(tài)性分子標(biāo)記篩選

1: 東鄉(xiāng)野生稻; 2: 日本晴; 3: 東鄉(xiāng)野生稻/日本晴雜合體; M: 標(biāo)準(zhǔn)分子量。

1: Dongxiang wild rice; 2: Nipponbare; 3: artificial F1hybrid between Dongxiang wild rice and Nipponbare; M: DNA size marker.

表1 各染色體多態(tài)性標(biāo)記分布

為了方便研究者直觀地查看標(biāo)記在染色體上的分布, 利用MapDraw (2.1)軟件, 以標(biāo)記的物理位置為參數(shù)制作了203個多態(tài)性標(biāo)記在染色體上的分布圖譜(圖2)。由圖譜可知, 標(biāo)記在染色體上的分布比較均勻, 基本完全覆蓋了12條水稻染色體。

2.2 多態(tài)性分子標(biāo)記的應(yīng)用

為檢驗所開發(fā)的多態(tài)性分子標(biāo)記在其他亞洲栽培稻中的應(yīng)用價值,隨機(jī)選取了59個多態(tài)性分子標(biāo)記, 以5個秈稻和5個粳稻品種為樣本, 以東鄉(xiāng)野生稻和日本晴為對照檢測(圖 3)發(fā)現(xiàn), 在粳稻品種中有84.7%的基因型與日本晴一致, 11.6%的基因型與東鄉(xiāng)野生稻一致; 在秈稻品種中有73.9%的基因型與東鄉(xiāng)野生稻一致, 13.5%的基因型與日本晴一致。而在檢測到的異于對照基因型數(shù)量中, 秈稻(10.9%)要遠(yuǎn)遠(yuǎn)高于粳稻品種(1.7%)(表2和附表2)。因此, 如果利用本套多態(tài)性分子標(biāo)記進(jìn)行基因型分析, 其在分析東鄉(xiāng)野生稻/粳稻雜交后代群體基因型上的應(yīng)用價值高于分析東鄉(xiāng)野生稻/秈稻雜交后代群體。

3 討論

分子標(biāo)記是現(xiàn)代遺傳學(xué)和分子生物學(xué)研究最常用的技術(shù)之一, 隨著技術(shù)的發(fā)展, 分子標(biāo)記的種類也在不斷發(fā)展, 例如SNP、CNV標(biāo)記等。東鄉(xiāng)野生稻做為一種普通野生稻, 與亞洲栽培稻基因組具有較大的差異。直接利用常規(guī)的SSR標(biāo)記以及基于秈稻與粳稻設(shè)計的InDel標(biāo)記來研究東鄉(xiāng)野生稻, 或者開展分子標(biāo)記輔助選擇育種具有一定的局限性。例如, 本研究直接篩選常用的SSR和InDel標(biāo)記, 僅僅只有8%的引物存在多態(tài)性。因此, 開發(fā)一套覆蓋東鄉(xiāng)野生稻全基因組的, 并且與亞洲栽培稻間存在多態(tài)性的分子標(biāo)記至關(guān)重要。本研究在篩選現(xiàn)有分子標(biāo)記的基礎(chǔ)上, 利用東鄉(xiāng)野生稻基因組重測序信息設(shè)計的InDel標(biāo)記, 共獲得了203個多態(tài)性分子標(biāo)記, 均勻分布于12條染色體上, 標(biāo)記平均間隔1.9 Mb, 可以滿足基因初定位和標(biāo)記輔助選擇育種的需要。

表2 多態(tài)性標(biāo)記秈梗亞種基因型鑒定一致性分析

該分析選取了不同省份的5個粳稻(遼鹽241、松粳8、吉粳88、淮稻9號、合系6號)和5個秈稻(黃花占、白玉絲苗、揚(yáng)稻6號、湘晚秈12號、鄂中5號)品種, 利用59個多態(tài)性標(biāo)記進(jìn)行檢測分析。括號中的數(shù)字代表本基因型占總檢測基因型數(shù)目的比例。

Fivecultivars (Liaoyan 241, Songgeng 8, Jigeng 88, Huaidao 9, and Hexi 6) and fivecultivars (Huanghuazhan, Baiyusimiao, Yangdao 6, Xiangwanxian 12, Ezhong 5) were selected for genotype identification by 59 polymorphic molecular markers. The number in parenthesis represents the proportion of current genotype to total detected genotypes.

1: 東鄉(xiāng)野生稻; 2: 日本晴; 3: 遼鹽241; 4: 松粳8; 5: 吉粳88; 6: 淮稻9號; 7: 合系6號; 8: 黃花占; 9: 白玉絲苗; 10: 揚(yáng)稻6號; 11: 湘晚秈12號; 12: 鄂中5號。

1: Dongxiang wild rice; 2: Nipponbare; 3: Liaoyan 241; 4: Songgeng 8; 5: Jigeng 88; 6: Huaidao 9; 7: Hexi 6; 8: Huanghuazhan; 9: Baiyusimiao; 10: Yangdao 6; 11: Xiangwanxian 12; 12: Ezhong 5.

分子標(biāo)記的開發(fā)基本上都是基于2個或者多個特定樣本, 因此獲得的分子標(biāo)記的應(yīng)用范圍具有一定的局限性。本研究獲得的203個多態(tài)性分子標(biāo)記是基于東鄉(xiāng)野生稻和日本晴開發(fā)的, 通過對5個秈稻和5個粳稻品種的檢驗, 發(fā)現(xiàn)其在粳稻品種中檢測的基因型84.7%與日本晴一致, 而在秈稻品種中的基因型73.9%與東鄉(xiāng)野生稻一致。如果利用該套分子標(biāo)記檢測東鄉(xiāng)野生稻與秈稻雜交的后代, 可能會有大約74%的分子標(biāo)記不能表現(xiàn)出多態(tài)性; 而在檢測東鄉(xiāng)野生稻與粳稻雜交后代中, 大約會有85%的標(biāo)記表現(xiàn)出多態(tài)性。因此, 該套分子標(biāo)記在分析東鄉(xiāng)野生稻與粳稻雜交后代群體的基因型上具有較大的應(yīng)用價值。

附表 請見網(wǎng)絡(luò)版: 1) 本刊網(wǎng)站http://zwxb.chinacrops. org/; 2) 中國知網(wǎng)http://www.cnki.net/; 3) 萬方數(shù)據(jù)http:// c.wanfangdata.com.cn/Periodical-zuowxb.aspx。

[1] 鄔柏梁, 何國成, 白國章, 吳藹如. 我省東鄉(xiāng)縣一帶發(fā)現(xiàn)野生稻. 江西農(nóng)業(yè)科技, 1979, (2): 6–7. Wu B L, He G C, Bai G Z, Wu A R. Wild rice found in Dongxiang County in our province., 1979, (2): 6–7 (in Chinese with English abstract).

[2] 潘熙淦, 饒憲章. 東鄉(xiāng)野生稻觀察及特征鑒定報告. 江西農(nóng)業(yè)科技, 1982, (7): 5–9. Pan X G, Rao X Z. Report of Dongxiang wild rice observation and characterization., 1982, (7): 5–9 (in Chinese with English abstract).

[3] 黃依南, 黃國勤. 東鄉(xiāng)野生稻的發(fā)現(xiàn)、價值與保護(hù). 農(nóng)業(yè)資源與環(huán)境學(xué)報, 2012, 29: 13–15. Huang Y N, Huang G Q. The discovery, value and protection of Dongxiang wild rice., 2012, 29: 13–15 (in Chinese with English abstract).

[4] 劉丹, 孫玉友, 魏才強(qiáng), 解忠, 李洪亮, 張巍巍, 程杜娟, 孫國宏, 徐德海. InDel分子標(biāo)記及其在水稻研究中的應(yīng)用. 種子, 2017, (9): 47–52. Liu D, Sun Y Y, Wei C Q, Xie Z, Li H L, Zhang W W, Cheng D J, Sun G H, Xu D H. InDel molecular marker and its application in rice () research., 2017, (9): 47–52 (in Chinese with English abstract).

[5] Shen Y J, Jiang H, Jin J P, Zhang Z B, Xi B, He Y Y, Wang G, Qian L L, Li X, Yu Q B, Liu H J, Chen D H, Gao J H, Huang H, Shi T L, Yang Z N. Development of genome-wide DNA polymorphism database for map-based cloning of rice genes., 2004, 135: 1198–1205.

[6] 馮芳君, 羅利軍, 李熒, 周立國, 徐小艷, 吳金, 陳宏偉, 陳亮, 梅捍衛(wèi). 水稻InDel和SSR標(biāo)記多態(tài)性的比較分析. 分子植物育種, 2005, 3: 725–730. Feng F J, Luo L J, Li Y, Zhou L G, Xu X Y, Wu J, Chen H W, Chen L, Mei H W. Comparative analysis of polymorphism of InDel and SSR markers in rice., 2005, 3: 725–730 (in Chinese with English abstract).

[7] 初志戰(zhàn), 郭海濱, 曾棟昌, 劉耀光. 秈粳稻基因組295個InDel 標(biāo)記的開發(fā). 作物學(xué)報, 2016, 42: 932–941. Chu Z Z, Guo H B, Zeng D C, Liu Y G. Development of 295 InDel markers forandrice.2016, 42: 932–941 (in Chinese with English abstract).

[8] Taillon-Miller P, Gu Z J, Li Q, Hillier L, Kwok P Y. Overlapping genomic sequences: a treasure trove of single nucleotide polymorphisms., 1998, 8: 748–754.

[9] 唐立群, 肖層林, 王偉平. SNP分子標(biāo)記的研究及其應(yīng)用進(jìn)展. 中國農(nóng)學(xué)通報, 2012, 28(12): 154–158. Tang L Q, Xiao C L, Wang W P. Research and application progress of SNP markers., 2012, 28(12): 154–158 (in Chinese with English abstract).

[10] Yu H H, Xie W B, Wang J, Xing Y Z, Xu C G, Li X H, Xiao J H, Zhang Q F. Gains in QTL detection using an ultra-high density SNP map based on population sequencing relative to traditional RFLP/SSR markers., 2011, 6: e17595.

[11] Abyzov A, Urban A E, Snyder M, Gerstein M. CNVnator: an approach to discover, genotype, and characterize typical and atypical CNVs from family and population genome sequencing., 2011, 21: 974–984.

[12] 李德軍. 江西東鄉(xiāng)普通野生稻滲入系的構(gòu)建及高產(chǎn)QTL定位. 中國農(nóng)業(yè)大學(xué)博士學(xué)位論文, 北京, 2003. Li D J. Development of Introgression Lines of Common Wild Rice (Griff) from Dongxiang in Jiangxi Province and Mapping QTLs for Improving Yield. PhD Dissertationof China Agricultural University, Beijing, China, 2003 (in Chinese with English abstract).

[13] Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform., 2009, 25: 1754–1760.

[14] Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R. The Sequence Alignment/ Map format and SAMtools.,2009, 25: 2078–2079.

[15] Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data., 2010, 38: e164.

[16] Murray M G, Thompson W F. Rapid isolation of high molecular weight plant DNA., 1980, 8: 4321–4325.

[17] Bassam B J, Caetano-Anollés G, Gresshoff P M. Fast and sensitive silver staining of DNA in polyacrylamide gels., 1991, 196: 80–83.

Development of molecular markers polymorphic between Dongxiang wild rice andrice cultivar ‘Nipponbare’

MA Xiao-Ding1, TANG Jiang-Hong2, ZHANG Jia-Ni2, CUI Di1, LI Hui3, LI Mao-Mao3,*, and HAN Long-Zhi1,*

1National Key Facility for Crop Gene Resources and Genetic Improvement / Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China;2Chongqing Normal University, Chongqing 401331, China;3Rice Research Institute, Jiangxi Academy of Agricultural Sciences, Nanchang 330200, Jiangxi, China

Dongxiang wild rice is a type of common wild rice () that has the northernmost natural distribution of any wild rice species. The genome of Dongxiang wild rice differs from that of modern cultivated rice varieties (). At present, a set of molecular markers that covers the entire genome of Dongxiang wild rice is lacking. In this study, we used Dongxiang wild rice and therice variety “Nipponbare” as research materials. By screening the existing collection of 1017 SSR and InDel markers and the 217 InDel markers which were designed using Dongxiang wild rice genome resequencing information, we obtained a set of 203 markers polymorphic between Dongxiang wild rice and ‘Nipponbare’, which were relatively uniformly distributed on the 12 rice chromosomes, and basically covered the entire genome, with an average inter-locus distance of 1.9 Mb. Through the genotyping of fiveand fivevarieties, we concluded that those 203 polymorphic molecular markers have a high application value in genotype analysis of Dongxiang wild rice andrice offspring population. These results provide a powerful tool for exploring beneficial genes from Dongxiang wild rice as well as marker-assisted breeding and selection.

InDel marker; genome re-sequencing;rice

2018-07-09;

2018-10-08;

2018-11-06.

10.3724/SP.J.1006.2019.82037

韓龍植, E-mail: hanlongzhi@caas.cn; 黎毛毛, E-mail: lmm3056@163.com

E-mail: maxiaoding@caas.cn

本研究由國家重點(diǎn)研發(fā)計劃項目(2016YFD0100101, 2016YFD0100301), 國家自然科學(xué)基金項目(31501287, 31671664), 國家科技支撐計劃項目(2015BAD01B01-1), 中國農(nóng)業(yè)科學(xué)院科技創(chuàng)新工程項目, 國家農(nóng)作物種質(zhì)資源保護(hù)項目(2017NWB036-01, 2017NWB036-12-2), 國家農(nóng)作物種質(zhì)資源平臺(NICGR2017-001)資助。

This study was supported by the National Key Research and Development Program of China (2016YFD0100101, 2016YFD0100301), the National Natural Science Foundation of China (31501287, 31671664), the National Key Technology Support Program of China (2015BAD01B01-1), the Agricultural Science and Technology Innovation Program of CAAS, the National Crop Germplasm Conservation Project (2017NWB036-01, 2017NWB036-12-2), and the National Crop Germplasm Resources Platform Project (NICGR2017-001).

URL:http://kns.cnki.net/kcms/detail/11.1809.S.20181105.1326.016.html