利用高世代转录组测序挖掘控制南方大豆皱叶症候选基因

doi:10.3864/j.issn.0578-1752.2024.15.002

摘要/Abstract

摘要：

【目的】南方大豆皱叶症（southern soybean crinkle leaf disease，SSCLD）严重时可导致大豆减产40%左右，利用高世代分离群体转录组测序技术（high-generation segregating populations RNA-seq，HGRNA-seq）挖掘控制SSCLD的候选基因，为揭示SSCLD发生的分子机制提供数据支撑。【方法】以2个南方皱叶症症级差异材料及其衍生的F_2:7株系为研究材料，对双亲进行重测序，对12个F_2:7单株（7个皱叶，5个正常叶）进行单个样本转录组测序，利用转录组及亲本的SNP/InDel数据进行混池法定位分析，利用转录组表达量差异数据进行GO和KEGG功能注释和富集分析，并在定位区间附近开发7个KASP分子标记，利用230个F₂群体构建局部连锁图谱，对转录组数据定位结果进行验证。联合基因定位和转录组分析结果，筛选控制SSCLD的候选基因。【结果】利用混池法把控制皱叶症的基因位点CL12定位在大豆第12染色体末端39 231 651-40 705 115 bp的1 473 464 bp区间内，利用F₂群体将控制皱叶症的基因定位于39 743 275-40 948 295 bp的1 205 020 bp区间内，与混池法定位结果基本一致。GO注释结果显示，代谢过程包括免疫系统过程，以及对刺激的反应等，细胞组分主要与膜等相关，KEGG注释结果显示，生物系统通路中主要包括植物-病原菌互作和环境适应等通路。GO富集的表达差异基因主要同跨膜受体蛋白活性、蛋白磷酸化及信号受体活性等方面相关，KEGG富集到的最多差异表达基因（differently expressed genes，DEGs）主要在植物-病原菌互作和植物MAPK信号通路上。结合南方大豆皱叶症诱因特点，选择定位候选区间内与抗病等相关且在外显子上存在非同义突变或表达量有差异的基因作为候选基因，结合qRT-PCR验证，最终确定GLYMA_12G223100、GLYMA_12G223900、GLYMA_12G224100、GLYMA_12G231800和GLYMA_12G233000等5个基因为控制南方大豆皱叶症的候选基因。【结论】HGRNA-seq实现了RNA-seq和BSA-seq相结合，成功挖掘了控制SSCLD的候选基因。

关键词: 大豆, 转录组, 皱叶症, 定位, 分子标记, 二代测序

Abstract:

【Objective】 Severe SSCLD can lead to about 40% yield reduction of soybean. In this study, candidate genes associated with controlling SSCLD were identified by HGRNA-seq technology to provide data support for revealing the molecular mechanisms of SSCLD.【Method】Two parents with SSCLD resistance difference and their derived F_2:7 lines were used as the experimental materials. Two parents and 12 F_2:7 lines (7 crinkled leaves and 5 normal leaves) were separately subjected to re-sequenced and transcriptome sequencing. Mixed pool method was carried out for gene localization by using SNP/InDel data from the parents and their progenies. GO and KEGG functional annotation and enrichment analysis was carried out by using differentially expressed genes (DEGs) from RNA-seq. 7 KASP molecular markers were developed in the vicinity of localization interval, and 230 F₂ populations were used to construct a localized linkage map to verify the localization results of the transcriptome data. Candidate genes controlling SSCLD were screened by combining gene mapping and transcriptome analysis.【Result】The gene controlling SSCLD named CL12 was located on chromosome 12 within a 1 473 464 bp interval by using mixed pool method, ranging from 39 231 651 bp to 40 705 115 bp. At the same time, the gene was localized in the 1 205 020 bp interval from 39 743 275 to 40 948 295 bp by using the F₂ population, which was basically consistent with the results of the mixed-pool method. The GO annotation results showed that the metabolic processes included immune system processes and responses to stimuli, and cellular components were mainly related to membranes, etc. The KEGG annotation results showed that the biosystem pathways mainly included plant-pathogen interaction and environmental adaptation, etc. GO-enriched DEGs were mainly related to the activity of transmembrane receptor proteins, protein phosphorylation, and signaling receptors, etc., KEGG enriched DEGs were mainly related to plant-pathogen interaction and MAPK signaling pathway. Combined with the characteristics of the causal factors of SSCLD, genes within the candidate interval associated with disease resistance which had non-synonymous mutations in the coding exons, or showed differential expression were selected as candidate genes for SSCLD resistance. GLYMA_12G223100, GLYMA_12G223900, GLYMA_12G224100, GLYMA_12G231800 and GLYMA_12G233000 were finally identified as the candidate genes controlling SSCLD by qRT-PCR.【Conclusion】HGRNA-seq realized the combination of RNA-seq and BSA-seq, and successfully mined the candidate genes controlling SSCLD.

Key words: soybean, transcriptome, crinkle leaf disease, localization, molecular markers, second generation sequencing

陈文杰, 陈渊, 韦清源, 汤复跃, 郭小红, 陈淑芳, 覃夏燕, 韦荣昌, 梁江. 利用高世代转录组测序挖掘控制南方大豆皱叶症候选基因[J]. 中国农业科学, 2024, 57(15): 2914-2930.

CHEN WenJie, CHEN Yuan, WEI QingYuan, TANG FuYue, GUO XiaoHong, CHEN ShuFang, QIN XiaYan, WEI RongChang, LIANG Jiang. Identification of Candidate Genes Controlling SSCLD by Utilizing High-Generation Segregating Populations RNA-seq[J]. Scientia Agricultura Sinica, 2024, 57(15): 2914-2930.

图/表 13

图1

表1

表2

表3

图2

表4

图3

图4

图5

图6

图7

表5

图8

参考文献 53

[1]	陈文杰, 陈渊, 韦清源, 郭小红, 汤复跃, 赵团结, 梁江. 南方大豆皱叶发生时叶片形态变化及其对产量性状的影响. 南方农业学报, 2022, 53(2): 460-468.
	CHEN W J, CHEN Y, WEI Q Y, GUO X H, TANG F Y, ZHAO T J, LIANG J. Changes of leaf morphology during the occurrence of wrinkled leaves in southern soybean and its effects on yield traits. Journal of Southern Agriculture, 2022, 53(2): 460-468. (in Chinese)
[2]	REDDY M R, RONAGHI A, BRYANT J A. Differential responses of soybean genotypes to excess manganese in an acid soil. Plant and Soil, 1991, 134(2): 221-226.
[3]	SANTOS E F, KONDO SANTINI J M, PAIXÃO A P, JÚNIOR E F, LAVRES J, CAMPOS M, REIS A R. Physiological highlights of manganese toxicity symptoms in soybean plants: Mn toxicity responses. Plant Physiology and Biochemistry, 2017, 113: 6-19. doi: S0981-9428(17)30034-7 pmid: 28157580
[4]	KELLEY K B, WAX L M, HAGER A G, RIECHERS D E. Soybean response to plant growth regulator herbicides is affected by other postemergence herbicides. Weed Science, 2005, 53(1): 101-112.
[5]	HOPKINS E W. Leaf-wrinkle, a nutritional disorder of soybean. Plant Physiology, 1933, 8(2): 333-336.
[6]	SAMRETWANICH K, KITTIPAKORN K, CHIEMSOMBAT P, IKEGAMI M. Complete nucleotide sequence and genome organization of soybean crinkle leaf virus. Journal of Phytopathology, 2001, 149(6): 333-336.
[7]	CAMPOS R E, BEJERMAN N, NOME C, LAGUNA I G, RODRÍGUEZ PARDINA P. Bean yellow mosaic virus in soybean from Argentina. Journal of Phytopathology, 2014, 162(5): 322-325.
[8]	OCHAR K, SU B H, ZHOU M M, LIU Z X, GAO H W, LAMLOM S F, QIU L J. Identification of the genetic locus associated with the crinkled leaf phenotype in a soybean (Glycine max L.) mutant by BSA-Seq technology. Journal of Integrative Agriculture, 2022, 21(12): 3524-3539.
[9]	SONG X F, WEI H C, CHENG W, YANG S X, ZHAO Y X, LI X, LUO D, ZHANG H, FENG X Z. Development of INDEL markers for genetic mapping based on whole genome resequencing in soybean. G3 (Bethesda, Md), 2015, 5(12): 2793-2799.
[10]	WANG Y Q, CHEN W, ZHANG Y, LIU M F, KONG J J, YU Z P, JAFFER A M, GAI J Y, ZHAO T J. Identification of two duplicated loci controlling a disease-like rugose leaf phenotype in soybean. Crop Science, 2016, 56(4): 1611-1618.
[11]	陈文杰, 陈渊, 韦清源, 郭小红, 汤复跃, 杨萌, 叶万典, 梁江. 一种大豆皱叶症发生特性及材料症级鉴定. 大豆科学, 2020, 39(3): 431-441.
	CHEN W J, CHEN Y, WEI Q Y, GUO X H, TANG F Y, YANG M, YE W D, LIANG J. Occurrence characteristics and identification of resistant materials for soybean crinkle leaf disease. Soybean Science, 2020, 39(3): 431-441. (in Chinese)
[12]	陈文杰, 梁江, 宁德娇, 韦清源, 汤复跃, 郭小红, 梁俊, 陈渊. 南方大豆皱叶症诱因分析研究. 大豆科学, 2022, 41(3): 300-307.
	CHEN W J, LIANG J, NING D J, WEI Q Y, TANG F Y, GUO X H, LIANG J, CHEN Y. Study on inducement of soybean crinkle leaf in southern China. Soybean Science, 2022, 41(3): 300-307. (in Chinese)
[13]	陈文杰, 梁江, 韦清源, 汤复跃, 郭小红, 陈渊. 基于转录组测序的辅助鉴定南方大豆皱叶症分子标记开发. 南方农业学报, 2023, 54(6): 1587-1597.
	CHEN W J, LIANG J, WEI Q Y, TANG F Y, GUO X H, CHEN Y. Development of molecular markers based on transcriptome sequencing for assisting the identification of southern soybean crinkle leaf disease. Journal of Southern Agriculture, 2023, 54(6): 1587-1597. (in Chinese)
[14]	TAKAGI H, ABE A, YOSHIDA K, KOSUGI S, NATSUME S, MITSUOKA C, UEMURA A, UTSUSHI H, TAMIRU M, TAKUNO S, INNAN H, CANO L M, KAMOUN S, TERAUCHI R. QTL-seq: Rapid mapping of quantitative trait loci in rice by whole genome resequencing of DNA from two bulked populations. The Plant Journal, 2013, 74(1): 174-183. doi: 10.1111/tpj.12105 pmid: 23289725
[15]	YAN M, LI F, SUN Q P, ZHAO J R, MA Y. Identification of chilling-tolerant genes in maize via bulked segregant analysis sequencing. Environmental and Experimental Botany, 2023, 208: 105234.
[16]	KLEIN H, XIAO Y G, CONKLIN P A, GOVINDARAJULU R, KELLY J A, SCANLON M J, WHIPPLE C J, BARTLETT M. Bulked-segregant analysis coupled to whole genome sequencing (BSA-seq) for rapid gene cloning in maize. G3, 2018, 8(11): 3583-3592.
[17]	PANDEY M K, KHAN A W, SINGH V K, VISHWAKARMA M K, SHASIDHAR Y, KUMAR V, GARG V, BHAT R S, CHITIKINENI A, JANILA P, GUO B Z, VARSHNEY R K. QTL-seq approach identified genomic regions and diagnostic markers for rust and late leaf spot resistance in groundnut (Arachis hypogaea L.). Plant Biotechnology Journal, 2017, 15(8): 927-941.
[18]	LUO H Y, PANDEY M K, KHAN A W, GUO J B, WU B, CAI Y, HUANG L, ZHOU X J, CHEN Y N, CHEN W G, LIU N, LEI Y, LIAO B S, VARSHNEY R K, JIANG H F. Discovery of genomic regions and candidate genes controlling shelling percentage using QTL-seq approach in cultivated peanut (Arachis hypogaea L.). Plant Biotechnology Journal, 2019, 17(7): 1248-1260.
[19]	SONG J, LI Z, LIU Z X, GUO Y, QIU L J. Next-generation sequencing from bulked-segregant analysis accelerates the simultaneous identification of two qualitative genes in soybean. Frontiers in Plant Science, 2017, 8: 919. doi: 10.3389/fpls.2017.00919 pmid: 28620406
[20]	KONG L P, LU S J, WANG Y P, FANG C, WANG F F, NAN H Y, SU T, LI S C, ZHANG F G, LI X M, ZHAO X H, YUAN X H, LIU B H, KONG F J. Quantitative trait locus mapping of flowering time and maturity in soybean using next-generation sequencing-based analysis. Frontiers in Plant Science, 2018, 9: 995. doi: 10.3389/fpls.2018.00995 pmid: 30050550
[21]	ZHANG X L, WANG W B, GUO N, ZHANG Y Y, BU Y P, ZHAO J M, XING H. Combining QTL-seq and linkage mapping to fine map a wild soybean allele characteristic of greater plant height. BMC Genomics, 2018, 19(1): 226.
[22]	DA SILVA M P, ZACCARON A Z, BLUHM B H, RUPE J C, WOOD L, MOZZONI L A, MASON R E, YINGLING S, PEREIRA A. Bulked segregant analysis using next-generation sequencing for identification of genetic loci for charcoal rot resistance in soybean. Physiological and Molecular Plant Pathology, 2020, 109: 101440.
[23]	陈浣, 夏菲, 吴新儒, 孙玉合. 集群分离分析法在植物基因定位上的应用. 基因组学与应用生物学, 2016, 35(6): 1546-1551.
	CHEN H, XIA F, WU X R, SUN Y H. The application of bulked segregant analysis in plant gene mapping. Genomics and Applied Biology, 2016, 35(6): 1546-1551. (in Chinese)
[24]	AMINFAR Z, RABIEI B, TOHIDFAR M, MIRJALILI M H. Identification of key genes involved in the biosynthesis of triterpenic acids in the mint family. Scientific Reports, 2019, 9(1): 15826.
[25]	MA Y F, MACKON E, JEAZET DONGHO EPSE MACKON G C, ZHAO Y T, LI Q F, DAI X G, YAO Y H, XIA X Z, NONG B X, LIU P Q. Combined analysis of BSA-Seq based mapping, RNA-Seq, and metabolomic unraveled candidate genes associated with panicle grain number in rice (Oryza sativa L.). Biomolecules, 2022, 12(7): 918.
[26]	GUO Z H, CAI L J, CHEN Z Q, WANG R Y, ZHANG L M, GUAN S W, ZHANG S H, MA W D, LIU C X, PAN G J. Identification of candidate genes controlling chilling tolerance of rice in the cold region at the booting stage by BSA-Seq and RNA-Seq. Royal Society Open Science, 2020, 7(11): 201081.
[27]	ZHENG R, DENG M, LV D, TONG B, LIU Y Q, LUO H B. Combined BSA-seq and RNA-seq reveal genes associated with the visual stay-green of maize (Zea mays L.). International Journal of Molecular Sciences, 2023, 24(24): 17617.
[28]	YAN P S, LI W H, ZHOU E X, XING Y, LI B, LIU J, ZHANG Z H, DING D, FU Z Y, XIE H L, TANG J H. Integrating BSA-seq with RNA-seq reveals a novel fasciated Ear5 mutant in maize. International Journal of Molecular Sciences, 2023, 24(2): 1182.
[29]	李世宽, 洪慧龙, 付佳祺, 谷勇哲, 孙如建, 邱丽娟. BSA-Seq结合RNA-Seq技术挖掘大豆叶片提前黄化衰老基因. 作物学报, 2024, 50(2): 294-309. doi: 10.3724/SP.J.1006.2023.34062
	LI S K, HONG H L, FU J Q, GU Y Z, SU R J, QIU L J. Mine the genes of premature yellowing and aging in soybean leaves by BSA-seq combined with RNA-seq technology. Acta Agronomica Sinica, 2024, 50(2): 294-309. (in Chinese)
[30]	曾维英, 赖振光, 孙祖东, 杨守臻, 陈怀珠, 唐向民. 基于BSA-Seq和RNA-Seq方法鉴定大豆抗豆卷叶螟候选基因. 作物学报, 2021, 47(8): 1460-1471. doi: 10.3724/SP.J.1006.2021.04195
	ZENG W Y, LAI Z G, SUN Z D, YANG S Z, CHEN H Z, TANG X M. Identification of the candidate genes of soybean resistance to bean pyralid (Lamprosema indicata Fabricius) by BSA-Seq and RNA-Seq. Acta Agronomica Sinica, 2021, 47(8): 1460-1471. (in Chinese)
[31]	NAN Y Y, XIE Y Y, HE H Y, WU H, GAO L X, ATIF A, ZHANG Y F, TIAN H, HUI J, GAO Y J. Integrated BSA-seq and RNA-seq analysis to identify candidate genes associated with nitrogen utilization efficiency (NUtE) in rapeseed (Brassica napus L.). International Journal of Biological Macromolecules, 2024, 254(Pt1): 127771.
[32]	YE S H, YAN L, MA X W, CHEN Y P, WU L M, MA T T, ZHAO L, YI B, MA C Z, TU J X, SHEN J X, FU T D, WEN J. Combined BSA-seq based mapping and RNA-seq profiling reveal candidate genes associated with plant architecture in Brassica napus. International Journal of Molecular Sciences, 2022, 23(5): 2472.
[33]	GAO Y B, YUAN Y H, ZHANG X Y, SONG H, YANG Q H, YANG P, GAO X L, GAO J F, FENG B L. Conuping BSA-seq and RNA-seq reveal the molecular pathway and genes associated with the plant height of foxtail millet (Setaria italica). International Journal of Molecular Sciences, 2022, 23(19): 11824.
[34]	LIU S Z, YEH C T, TANG H M, NETTLETON D, SCHNABLE P S. Gene mapping via bulked segregant RNA-seq (BSR-seq). PLoS ONE, 2012, 7(5): e36406.
[35]	KIM D, PERTEA G, TRAPNELL C, PIMENTEL H, KELLEY R, SALZBERG S L. TopHat2: Accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biology, 2013, 14(4): R36.
[36]	TRAPNELL C, WILLIAMS B A, PERTEA G, MORTAZAVI A, KWAN G, VAN BAREN M J, SALZBERG S L, WOLD B J, PACHTER L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechnology, 2010, 28(5): 511-515. doi: 10.1038/nbt.1621 pmid: 20436464
[37]	LI B, DEWEY C N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics, 2011, 12: 323. doi: 10.1186/1471-2105-12-323 pmid: 21816040
[38]	TANG H B, SATO P F K, RAMIREZ F, KLOPFENSTEIN D, MUNGALL C, YUNES J M, PEDERSEN B S. GOATOOLS: Tools for gene ontology. Zenodo, 2015.
[39]	LI H, DURBIN R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics (Oxford, England), 2010, 26(5): 589-595.
[40]	MCKENNA A, HANNA M, BANKS E, SIVACHENKO A, CIBULSKIS K, KERNYTSKY A, GARIMELLA K, ALTSHULER D, GABRIEL S, DALY M, DEPRISTO M A. The genome analysis toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Research, 2010, 20(9): 1297-1303. doi: 10.1101/gr.107524.110 pmid: 20644199
[41]	YANG L M, WANG J G, HAN Z H, LEI L, LIU H L, ZHENG H L, XIN W, ZOU D T. Combining QTL-seq and linkage mapping to fine map a candidate gene in qCTS6 for cold tolerance at the seedling stage in rice. BMC Plant Biology, 2021, 21(1): 278.
[42]	LEI L, ZHENG H L, BI Y L, YANG L M, LIU H L, WANG J G, SUN J, ZHAO H W, LI X W, LI J M, LAI Y C, ZOU D T. Identification of a major QTL and candidate gene analysis of salt tolerance at the bud burst stage in rice (Oryza sativa L.) using QTL-seq and RNA-seq. Rice, 2020, 13(1): 55.
[43]	LI Y H, LIN R M, HU J H, SHI X H, QIU D, WU P P, GOITOM G H, WANG S Q, ZHANG H J, YANG L, LIU H W, WU Q H, XIE J Z, ZHOU Y, LIU Z Y, LI H J. Mapping of wheat stripe rust resistance gene Yr041133 by BSR-seq analysis. The Crop Journal, 2022, 10(2): 447-455.
[44]	WANG Y, ZHANG H Z, XIE J Z, GUO B M, CHEN Y X, ZHANG H Y, LU P, WU Q H, LI M M, ZHANG D Y, GUO G H, YANG J, ZHANG P P, ZHANG Y, WANG X C, ZHAO H, CAO T J, LIU Z Y. Mapping stripe rust resistance genes by BSR-Seq: YrMM58 and YrHY1 on chromosome 2AS in Chinese wheat lines Mengmai 58 and Huaiyang 1 are Yr17. The Crop Journal, 2018, 6(1): 91-98.
[45]	LI M M, LI B B, GUO G H, CHEN Y X, XIE J Z, LU P, WU Q H, ZHANG D Y, ZHANG H Z, YANG J, ZHANG P P, ZHANG Y, LIU Z Y. Mapping a leaf senescence gene els1 by BSR-Seq in common wheat. The Crop Journal, 2018, 6(3): 236-243.
[46]	WU L R, ZHU T, HE H G, CAO X Y, LI H S, XU H X, JIA M S, ZHANG L P, SONG J C, MIRZAGHADERI G, LIU C, MA P T. Genetic dissection of the powdery mildew resistance in wheat breeding line LS5082 using BSR-Seq. The Crop Journal, 2022, 10(4): 1120-1130.
[47]	陈正杰, 宛永璐, 钟文娟, 陈四维, 周永航, 石盛佳, 蒋理, 戢沛城, 杨泽湖, 毛正轩, 牟方生. 基于大豆基因组重测序的InDel标记开发及应用. 植物遗传资源学报, 2021, 22(3): 815-833. doi: 10.13430/j.cnki.jpgr.20201028003
	CHEN Z J, WAN Y L, ZHONG W J, CHEN S W, ZHOU Y H, SHI S J, JIANG L, JI P C, YANG Z H, MAO Z X, MOU F S. Development and application of soybean InDel markers based on whole-genome resequencing datasets. Journal of Plant Genetic Resources, 2021, 22(3): 815-833. (in Chinese)
[48]	FANG Y D, WANG L W, SAPEY E, FU S, WU T T, ZENG H Y, SUN X G, QIAN S Q, KHAN M A A, YUAN S, WU C X, HOU W S, SUN S, HAN T F. Speed-breeding system in soybean: integrating off-site generation advancement, fresh seeding, and marker-assisted selection. Frontiers in Plant Science, 2021, 12: 717077.
[49]	郭小红, 房裕东, 韦清源, 汤复跃, 陈文杰, 梁江, 韩天富, 陈渊. 黄淮海夏大豆品种中黄39在南宁的三季繁育可行性研究. 大豆科学, 2020, 39(6): 856-861.
	GUO X H, FANG Y D, WEI Q Y, TANG F Y, CHEN W J, LIANG J, HAN T F, CHEN Y. Three successive generation feasibility of Huang-Huai-Hai summer soybean cultivar Zhonghuang 39 in Nanning of South China. Soybean Science, 2020, 39(6): 856-861. (in Chinese)
[50]	ROBLES J A, QURESHI S E, STEPHEN S J, WILSON S R, BURDEN C J, TAYLOR J M. Efficient experimental design and analysis strategies for the detection of differential expression using RNA-Sequencing. BMC Genomics, 2012, 13: 484. doi: 10.1186/1471-2164-13-484 pmid: 22985019
[51]	ASANO T, YOSHIOKA Y, KUREI S, SAKAMOTO W, MACHIDA Y, SODMERGEN. A mutation of the CRUMPLED LEAF gene that encodes a protein localized in the outer envelope membrane of plastids affects the pattern of cell division, cell differentiation, and plastid division in Arabidopsis. The Plant Journal, 2004, 38(3): 448-459.
[52]	CHEN Y L, ASANO T, FUJIWARA M T, YOSHIDA S, MACHIDA Y, YOSHIOKA Y. Plant cells without detectable plastids are generated in the crumpled leaf mutant of Arabidopsis thaliana. Plant & Cell Physiology, 2009, 50(5): 956-969.
[53]	ŠIMKOVÁ K, KIM C H, GACEK K, BARUAH A, LALOI C, APEL K. The chloroplast division mutant caa33 of Arabidopsis thaliana reveals the crucial impact of chloroplast homeostasis on stress acclimation and retrograde plastid-to-nucleus signaling. The Plant Journal, 2012, 69(4): 701-712.

标记 Markers	Primer_AlleleFAM (5′-3′)	Primer_AlleleHEX (5′-3′)	Primer_Common (5′-3′)	AlleleFAM	AlleleHEX
Chr.12_38720327	gaaggtgaccaagttcatgctaaagcatgctttaatattttggA	gaaggtcggagtcaacggattaagcatgctttaatattttggG	GCAATTGAAAACATAAGTCACGCCTAC	A	G
Chr.12_39553364	gaaggtgaccaagttcatgcttaaatcgtgtttcttttcctcctgC	gaaggtcggagtcaacggatttgtaaatcgtgtttcttttcctcctgT	CCCAACGGCTGCAATGGTC	C	T
Chr.12_39643002	gaaggtgaccaagttcatgctatggttggaccagacaggtaaggG	gaaggtcggagtcaacggattatggttggaccagacaggtaaggA	CTCTTTGCGCTGCTCGTCATTACTC	G	A
Chr.12_39743275	gaaggtgaccaagttcatgctgggagtccttgttatgtattatgatggC	gaaggtcggagtcaacggattgagggagtccttgttatgtattatgatggT	ATCCCCTAAGAAGTTGATACCACAT	C	T
Chr.12_40239112	gaaggtgaccaagttcatgctggaggcgactaaaattgtttccaaT	gaaggtcggagtcaacggattggaggcgactaaaattgtttccaaC	AAATTAGCCAAACACGCATGGTTCA	T	C
Chr.12_40443647	gaaggtgaccaagttcatgctccgaattgtgttttaagtcatatcT	gaaggtcggagtcaacggattccgaattgtgttttaagtcatatcC	GCAATATCAAGAATGAAACTTCT	T	C
Chr.12_40948295	gaaggtgaccaagttcatgctattgaaccatgttagcacagcagtT	gaaggtcggagtcaacggattattgaaccatgttagcacagcagtA	TAACAGAATGCTAACTTGTGATACT	T	A

样本 Sample	过滤得到的reads数 Clean reads	过滤后的总碱基数 Clean bases	Q30 (%)	GC含量 GC content (%)	定位到基因组总Reads占比 Mapped Reads radio (%)
GC8	37329976	11235520244	91.72	35.15	99.54
Y2017-1	36749923	11054768304	92.28	35.69	98.64
C_10	48149604	7164122179	96.23	45.71	94.20
C_11	47618902	7090388779	96.57	45.95	94.19
C_18	45335046	6734925207	96.65	45.6	95.18
C_23	49603440	7381504895	96.56	46.14	95.12
C_24	49929876	7430549886	96.52	45.69	95.35
C_3	41997550	6249241100	96.59	45.94	94.70
C_6	45105078	6703705471	96.55	46.07	94.78
N_1	48155648	7150543340	96.44	45.87	95.00
N_3	44567696	6617278005	96.44	45.86	93.82
N_5	47817526	7091029900	96.66	46.07	94.92
N_6	42285914	6296598949	95.26	45.57	94.52
N_7	40742934	6065420513	96.49	46.24	94.91

染色体 Chromosome	长度 Length (bp)	来自GC8 From GC8	来自Y2017-1 From Y2017-1	相同位点 Same site	杂合位点 Hybridization site	SNP/InDel数 No. of SNP/ InDel site	位点总数 Total numbers of SNP/InDel site	纯合率 Purity ratio (%)	标记密度 Density of SNP/ InDel (SNP/Mb)
Chr.1	57932355	150	285	607	11	435	1053	98.96	7.70
Chr.2	50400358	143	218	616	10	361	987	98.99	7.36
Chr.3	46951866	68	334	1107	159	402	1668	90.47	11.95
Chr.4	51203389	159	156	602	26	315	943	97.24	6.66
Chr.5	42274530	299	1	854	239	300	1393	82.84	12.75
Chr.6	50945864	326	217	1196	54	543	1793	96.99	11.72
Chr.7	44949256	479	125	1091	188	604	1883	90.02	17.62
Chr.8	47227184	119	501	820	476	620	1916	75.16	23.21
Chr.9	50572668	339	201	984	284	540	1808	84.29	16.29
Chr.10	51638687	84	229	616	70	313	999	92.99	7.42
Chr.11	39643745	35	277	584	27	312	923	97.07	8.55
Chr.12	41531199	106	77	690	131	183	1004	86.95	7.56
Chr.13	45225048	159	88	1399	422	247	2068	79.59	14.79
Chr.14	49893278	240	176	679	103	416	1198	91.40	10.40
Chr.15	53754295	84	271	955	185	355	1495	87.63	10.05
Chr.16	38112070	637	191	1193	157	828	2178	92.79	25.84
Chr.17	41740656	24	425	597	23	449	1069	97.85	11.31
Chr.18	58286270	404	357	1975	159	761	2895	94.51	15.78
Chr.19	51272880	94	94	924	25	188	1137	97.80	4.15
Chr.20	47846026	323	14	776	16	337	1129	98.58	7.38
总计/平均 Total/mean	961401624	4272	4237	18265	2765	8509	29539	91.61	11.73

基因编号 Gene ID	标记 Marker	正向引物 Forward primer (5′-3′)	反向引物 Reverse primer (5′-3′)	退火温度 Tm (℃)
GLYMA_12G223900	qPCR-12G223900	TGTGTAGTGGCTTGTCATAT	CTTAGTCTTATGGAGCAATCAC	60
GLYMA_12G224000	qPCR-12G224000	ATTGATGTTGGTGCTTCTTG	TGAATGTGTTGGTGCTTCTA	60
GLYMA_12G224100	qPCR-12G224100	CATCAACCCAACTCACTCAA	GTGGCAACATAAACTCTTGT	60
GLYMA_12G233000	qPCR-12G233000	GCTTCAACCGATCAACTTAT	CCACCATACATATAGCAAGC	60