油菜耐盐相关性状的全基因组关联分析及其候选基因预测

【目的】通过对甘蓝型油菜耐盐相关性状进行全基因组关联分析,寻找可能与油菜耐盐性相关的SNP位点,发掘与油菜耐盐性有关的候选基因。【方法】以1.2%NaCl溶液作为培养液,去离子水为对照,对307个不同品系的甘蓝型油菜进行发芽试验。播种后7 d测定幼苗根长、鲜重及发芽率,计算盐胁迫下各性状相对值,并作为评价耐盐的指标。结合油菜60K SNP芯片,利用SPAGeDi v1.4软件对该群体307份甘蓝型油菜进行亲缘关系分析,并计算亲缘关系值的矩阵。利用软件STRUCTURE v2.3.4对关联群体进行了群体结构分析。为有效排除假关联的影响,将群体结构和材料间的亲缘关系考虑进关联分析中,同时进行了PCA+K模型、Q+K模型以及K模型3种混合线性模型分析和比较,根据所有SNP的–lg(P)观察值和期望值,确定每个性状GWAS分析的最优模型。采用TASSEL 5.0软件,在最优模型下对307份材料耐盐各性状的相对值分别与SNP标记进行全基因组关联分析。利用油菜基因组数据库,在显著SNP位点侧翼序列200 kb范围内提取基因。根据拟南芥中已经明确功能的耐盐相关基因,筛选出目标基因组区段内与耐盐相关的油菜同源基因。【结果】全基因组关联分析共检测到164个与根长显著关联的SNP位点,23个与鲜重显著关联的SNP位点,38个与发芽率显著关联的SNP位点。其中与根长、鲜重、发芽率最显著关联的SNP位点分别位于染色体A08、A02和A06,贡献率分别为23.84%、18.59%和31.81%。在这些显著SNP位点侧翼序列200 kb范围内发掘出可能与油菜耐盐性有关的50个候选基因。这些候选基因主要包括转录因子MYB、WRKY、ABI1、b ZIP、ERF1、CZF1、XERICO等以及一些下游受转录因子调控的不同功能基因NHX1、PTR3、CAT1、HKT、CAX1、ACER、STH、STO等。在根长和发芽率2个不同耐盐性状的分析结果中均筛选出位于A03染色体上的耐盐基因BnaA03g14410D。另外,这些耐盐候选基因中包含两组串联重复基因,分别是位于A03染色体上的BnaA03g18900D和BnaA03g18910D,位于C09染色体上的BnaC09g19080D、BnaC09g19090D和BnaC09g19100D。除此之外,耐盐候选基因中还包含2个距离非常近(中间只间隔2个基因)的重复基因BnaC02g39600D和BnaC02g39630D。【结论】共检测到225个与耐盐性状显著关联的SNP位点,筛选出50个可能与油菜耐盐性有关的候选基因。     

甘蓝型油菜品质相关性状全基因组关联分析

采用近红外方法,分析104份甘蓝型油菜种子中的芥酸、油酸与硫代葡萄糖苷含量;运用简化基因组技术(ddRADseq)对全部材料进行基因分型;结合性状与基因型的数据,利用Tassel混合线性模型进行全基因组关联分析。结果表明:芥酸含量与硫代葡萄糖苷呈极显著正相关,油酸含量与芥酸、硫代葡萄糖苷呈极显著负相关;群体结构分析将104份材料划分为4个亚群;通过全基因组关联分析发现,与芥酸、油酸含量显著关联的SNP标记主要分布在A08与C03染色体上,与硫代葡萄糖苷含量显著关联的标记分布在A05、A02与C09染色体上。     

苯磺隆胁迫下油菜萌发期相关性状的全基因组关联分析

目的 研究苯磺隆残留对油菜种子萌发的影响,运用全基因组关联分析（genome-wide association analysis,GWAS）揭示苯磺隆胁迫下油菜萌发期相关性状的遗传因子和候选基因,探究油菜在苯磺隆逆境胁迫下的生理形态所反映的基因调控机制,为耐苯磺隆油菜品种的研究提供参考。方法 以241份甘蓝型油菜品种（系）为材料、25 mg·L ^-1苯磺隆溶液为处理液、蒸馏水为对照进行发芽试验。发芽7 d测定并计算相对发芽率、相对根长和相对鲜重。结合芸薹属60K SNP芯片分析群体基因型,通过STRUCTURE软件和TASSEL软件分别对该群体进行群体结构分析以及亲缘关系和LD衰减分析。为有效排除假关联的影响,采用一般线性模型（GLM）和混合线性模型（MLM）中的6种模型进行比较,确定每个性状GWAS分析的最优模型。同时,利用TASSEL软件在最优模型下对241份材料的3个性状分别进行全基因组关联分析,根据关联SNP位点的LD区间序列预测候选基因。 结果 241份品种（系）群体可分为P1（94份材料）和P2（147份材料）2个亚群,其中约56.28%的材料之间的亲缘关系值为0。全基因组关联分析（K+PCA模型）共检测到16个与性状显著关联的SNP位点,这些位点可解释9.42%—13.14%的表型变异率。通过分析显著SNP位点的LD区间与甘蓝型油菜对应的区间序列,筛选出25个候选基因可能与油菜耐苯磺隆有关,其中9个为细胞色素P450家族基因,5个参与谷胱甘肽合成或代谢过程,2个为多药耐药相关蛋白基因。同时发现与相对发芽率显著相关的基因ATGSTU19编码谷胱甘肽转移酶,参与毒素分解过程,在各种胁迫反应中起重要作用。在相对根长和相对鲜重共同鉴定到的候选基因BnaC02g27690D功能未知。结论 共检测到16个SNP位点与耐苯磺隆性状显著关联,筛选出25个候选基因可能与油菜耐苯磺隆有关。     

甘蓝型油菜角果长度全基因组关联分析

【目的】挖掘与油菜角果长度性状显著相关的SNP位点及候选基因,为揭示油菜角果长度性状的遗传基础和分子机制提供理论依据,为油菜产量分子标记辅助选择育种奠定基础。【方法】在江西农业大学试验地和江西省红壤研究所试验地2个环境下考察300份甘蓝型油菜自交系的角果长度性状,利用简化基因组测序技术(specific locus amplified fragment sequencing,SLAF-seq)对300份甘蓝型油菜自交系基因组DNA进行测序并分析,利用获得的均匀分布于甘蓝型油菜基因组上的201 817个群体SNP(single nucleotide polymorphism,SNP)对角果长度性状进行全基因组关联分析(genome-wide association study,GWAS),探测与油菜角果长度显著相关的SNP位点,并基于群体连锁不平衡分析结果搜寻显著SNP位点两侧100 kb范围内的基因,通过BLAST获得关联区域内基因的注释信息,根据注释信息找出与性状相关的候选基因。【结果】农大试验地角果长度表型变异幅度为46.35—107.07 mm;红壤所试验地角果长度表型变异幅度为39.41—101.35 mm,两性状在2个环境下均表现出广泛表型变异。通过一般线性模型(general linear model,GLM)关联分析,农大环境下共检测到121个角果长度显著关联的SNP位点,分布在A04、A06、A08、A09、C02、C03、C06和C09等8条染色体上,其中,A09染色体上分布最多(83个SNP),红壤所环境下检测到22个角果长度显著关联的SNP位点,其中,1个在C09染色体上,其余21个均分布于A09染色体,在两地探测到20个一致性SNP位点;通过混合线性模型(mixed linear model,MLM)分析,农大环境下共检测到5个角果长度显著关联的SNP位点,其中,3个SNP位点与红壤所环境下检测到3个SNP位点一致,所有位点均位于A09染色体上。对MLM关联分析得到的显著SNP位点两侧100 kb区域内基因进行搜寻并进行功能注释,发现多个候选基因参与调节碳水化合物的运输与合成、花器官和种子的发育、信号转导等,它们可能通过上述功能影响油菜角果的生长,导致角果长度的差异。【结论】通过GLM和MLM两种分析方法探测到多个与油菜角果长度性状显著关联的基因位点,并在显著性位点附近搜寻到相关候选基因。     

甘蓝型油菜polCMS育性恢复位点的全基因组关联分析

【目的】甘蓝型油菜波里马细胞质雄性不育(pol CMS)在中国已被广泛应用于杂交种育种,其育性恢复程度表现出受1对主效基因的控制,并受微效修饰基因的影响。通过全基因组关联分析方法挖掘育性恢复位点,并对候选基因进行比较分析。【方法】通过芸薹属60K SNP芯片对308份甘蓝型油菜自然群体进行基因型分型,并用pol CMS系301A作母本,与上述材料分别进行杂交得到308份F1,每份F-_1分别于2013年和2014年进行种植,每年2次重复,于始花期根据花粉育性和花蕊发育情况调查F1植株的育性等级,同时对测交父本自然群体进行群体结构分析和亲缘关系评估,并结合测交父本的基因型分型结果和F_1的育性等级进行全基因组关联分析(GWAS)。从GWAS分析中显著的SNP左右100 kb区间或与显著SNP处于同一单体型块(R~20.5)的区间内预测候选基因,并对候选基因进行QTL比较分析和单体型或等位基因的效应分析。【结果】方差分析结果显示,两年F_1的育性等级存在显著差异(P0.01),但相关分析发现,两年的育性等级存在显著的正相关(r=0.52,P0.001)。群体结构分析显示,所有测交父本被分为3个亚群(冬性、春性和半冬性),亲缘关系分析发现,任何2个材料之间平均亲缘关系值为0.072,73%的任意材料间亲缘关系值小于0.1,其中,约53%的材料亲缘关系值为0。GWAS分析共检测到13个与育性恢复程度显著关联的SNP,构成了6个候选区间,分别位于A01、A09、C03、C06和C08 5条染色体上,单个SNP解释的表型变异介于2.53%—9.96%。从中共预测到6个与育性恢复位点相关的候选基因,其中4个编码的蛋白含有恢复基因特有的PPR保守基序。共线性分析发现,4个候选基因中的2个(Bna A09g46700D和Bna C08g40710D)位于A09和C08染色体部分同源区间,且与已克隆的pol CMS育性恢复位点ORF2同源。另外2个新鉴定到的候选基因(Bna C03g45840D和Bna C06g13000D)连锁的SNP等位基因或单体型变化都与育性等级显著相关(P0.001)。【结论】通过GWAS分析鉴定到多个与油菜育性恢复有关的候选基因,开发基于与这些基因连锁位点或SNP的功能标记将有助于对该不育系统进行恢复系和保持系的筛选。     

油菜主花序角果密度及其相关性状的全基因组关联分析

【目的】油菜高产是育种工作的主要研究目标之一。角果密度、主花序有效角果数等性状与产量都有显著或极显著的正相关关系,是油菜高产育种考查的主要性状。为揭示油菜角果密度及其相关性状的遗传机理和分子机制奠定基础。【方法】以不同遗传背景和地理来源的213份甘蓝型油菜品种（系）构成的自然群体为研究对象,利用芸薹属60K Illumina Infinium SNP芯片对该群体进行基因型分型。分别于2015年和2016年在成熟期调查该群体主花序有效长和主花序有效角果数,计算主花序角果密度。利用Structure 2.3.4软件对该群体进行群体结构分析,Tassel 5.1.0软件分析亲缘关系和染色体连锁不平衡的衰减;然后基于最优模型对主花序角果密度及其相关性状进行全基因组关联分析（genome-wide association analysis,GWAS）,依据关联SNP位点的LD区间序列,预测与性状相关的重要关联候选基因。【结果】群体结构分析显示,213份甘蓝型油菜分为P1和P2亚群,P1亚群包含50份材料（23.5%）,P2亚群包含163份材料（76.5%）,基本上和油菜的地理栽培属性一致;亲缘关系发现约89.74%材料之间的亲缘关系值小于0.2,其中约有59.91%材料的亲缘关系值为0。总体来看,整个自然群体材料之间的亲缘关系比较远。对A、C基因组进行连锁不平衡分析发现,A和C基因组的r2随着遗传距离的增加而下降,A基因组的衰减距离整体比C基因组的衰减距离小。GWAS分析两年数据共检测到17个SNP位点与主花序角果密度及其相关性状关联。其中与主花序角果密度和主花序有效长相关的SNP标记分别有7个和9个,并分别解释11.34%—15.96%和9.67%—13.10%的表型变异;与主花序有效角果数相关联的标记有1个,解释11.56%的表型变异。通过分析关联SNP位点的LD区间与甘蓝型油菜对应的区间序列,找到22个与主花序角果密度及其相关性状有关的候选基因,其中BnaA01g16940D、BnaC01g38800D和BnaA04g09170D等主要通过调控赤霉素和生长素等内源激素的合成和信号转导来控制主花序角果密度及其相关性状;BnaA01g16970D、BnaA03g29180D、BnaA03g29810D、BnaC01g39680D和 BnaC03g32770D通过对分生组织的调控来改变表型;BnaC09g18690D和 BnaC09g09210D等主要通过控制细胞分裂生长等过程改变表型。【结论】检测到17个SNP标记与油菜主花序角果密度、主花序有效长和主花序有效角果数关联,筛选出22个与主花序角果密度及其相关性状有关的候选基因。     

甘蓝型油菜GPDH基因家族的全基因组鉴定

在甘蓝型油菜基因组中鉴定出18个BnGPDH基因,参考拟南芥AtGPDH基因家族分类信息将它们分为BnGPDHp、BnGPDHm和BnGPDHc 3个亚家族。基因结构、保守模块分析发现亚家族内的基因成员具有较高的保守性;染色体位置定位显示BnGPDH分布在10条染色体上;qRT-PCR分析发现BnGPDHc家族成员在所检测的组织部位中具有不同的表达模式。此外,在花后21 d(day after flower, DAF)和30DAF的种子中 BnGPDHc1-2、 BnGPDHc2-1、 BnGPDHc2-2的相对表达量在高油材料中要显著高于低油材料,表明这些基因可能与油菜种子中三脂酰甘油的形成相关。     

人工海水胁迫下不同甘蓝型油菜品种发芽能力的差异

通过人工海水胁迫研究了不同甘蓝型油菜品种发芽能力的差异,结果表明:随海水盐度的增加,油菜种子萌发能力逐渐减弱;杂交油菜与常规油菜耐海水盐度能力相当,但中等盐度海水对杂交油菜的萌发抑制作用大于常规油菜;甘蓝型油菜品种种子萌发时耐海水盐度为1.31%,耐海水半致死盐度为1.70%,耐海水极限盐度为2.16%;甘蓝型油菜耐盐能力在品种间存在显著差异.     

整合关联分析和共表达网络分析挖掘甘蓝型油菜籽粒质量候选基因

甘蓝型油菜是中国重要的油料作物,籽粒质量是油菜产量构成的重要因素。本研究利用A-D Test模型对496份油菜的籽粒质量进行全基因组关联分析,共检测到19个显著位点,联合解释34.1%的表型变异。整合A-D Test及前期混合线性模型和一般线性模型结果后得到71个位点,联合解释50.1%的表型变异。在22个位点置信区间内找到ARF2、NPC6、TTG2和WRI1等已被报道的拟南芥中的籽粒质量基因的同源基因。同时,利用大粒品种中双11和中小粒品种中油821的籽粒和角果皮转录组数据进行加权基因共表达网络分析,构建了13个共表达模块,其中紫色(purple)和洋红色(magenta)模块与籽粒质量表型显著相关。GO富集分析结果表明,2个模块在L-苯丙氨酸氨基转移酶活性、硫双加氧酶活性、焦磷酸酶活性和RNA解旋酶活性等显著富集。2个模块中的枢纽基因BnaA06g00850D、BnaA01g00990D、BnaC06g10000D和BnaC02g44260D等的同源基因为ETHE1、DAR1、GLN1;1和SMG7等拟南芥籽粒质量已知基因。整合全基因组关联分析和加权基因共表达网络分析的分析结果,...     

水稻耐淹成苗率相关性状全基因组的关联分析

背景 耐淹成苗率低是限制直播稻产量的重要因素,挖掘高种子活力、低氧萌发能力强的水稻材料是提高耐淹成苗率的关键,但控制耐淹成苗率的遗传位点的挖掘仍然比较有限。目的 利用来源广泛的自然种质,分析影响耐淹成苗率的关键表型性状,挖掘相关的遗传位点和候选基因,为直播稻耐淹成苗机理研究提供一定的理论和材料基础。方法 以200份来源广泛的水稻种质为材料,在有氧环境下进行发芽试验,测量种子发芽率、发芽指数和活力指数;在低氧条件下测量芽鞘长和芽鞘直径;进行耐淹成苗试验,水深10 cm,20 d后测量耐淹成苗率。分析各性状间的相关性,挖掘影响耐淹成苗率的关键性状;利用简化基因组测序对以上6个表型进行全基因组关联分析,鉴定与性状显著关联的SNP位点,并在关联区间内筛选候选基因;对02428和YZX 2份材料进行有氧、无氧以及氧气含量转换条件下的转录组检测,结合全基因组关联分析结果,分析候选基因的表达模式差异。结果 种子活力、芽鞘表型和耐淹成苗率在200份材料间存在广泛的遗传变异,其中,芽鞘长和活力指数的变异系数最大;相关分析结果表明,芽鞘长、活力指数与耐淹成苗率呈极显著正相关;通过全基因组关联分析,共鉴定出8个与活力指数显著关联的位点,15个与芽鞘长显著关联的位点;结合基因组注释,在关联区间筛选出6个与活力指数相关的候选基因,7个与芽鞘长度相关的候选基因;进一步比较13个基因在有氧、无氧及氧气转换条件下的表达模式以及表达量的变化,发现Os02g0657000、Os03g0592500、Os08g0380100表达量变化显著,表现出对氧气处理的敏感性。结论 种子活力、芽鞘长与耐淹成苗率密切相关,可作为筛选高耐淹成苗水稻材料的重要性状。全基因组关联分析、转录组分析与基因表达模式比较的联合应用可提高候选基因的筛选效率。水稻耐淹成苗过程可能受到与逆境胁迫、光合作用相关基因的调控。     

白菜类作物开花时间的全基因组关联分析

【目的】解析白菜类作物开花时间的调控位点,定位白菜类作物开花时间相关的候选基因,为白菜类作物抽薹开花时间的遗传改良提供依据。【方法】以116份白菜类作物组成的自然群体作为研究材料,分别种植在温室与露地2个独立的环境中,进行开花时间调查。同时,提取试验材料的DNA样品进行深度为1.2x的重测序,对测序数据用Pooled Mapping法进行过滤、与参考基因组比对,获得全基因组高密度SNP集合。经过条件过滤后,对高质量的SNP集合进行生物信息分析,包括试验材料的群体结构分析和全基因组连锁不平衡分析。从高质量的SNP集合中,随机挑选出2 000个变异位点,用Phy ML软件以最大似然法对116份试验材料进行系统发育树分析。用全部的高质量SNP集合位点通过软件Haploview进行全基因组连锁不平衡分析。最后,将高质量的SNP集合与开花时间数据结合,通过TASSEL和GAPIT软件包以及R程序语言进行全基因组关联分析。根据强关联峰值信号点位置和连锁不平衡区间定位开花时间候选位点,再通过白菜与同源物种拟南芥的基因共线性关系以及基因功能注释分析来预测白菜类作物开花时间相关的候选基因。【结果】不同种植条件下、不同类型的白菜类作物在开花时间上存在广泛差异。试验材料在露地环境下的开花时间高峰期明显早于温室环境下的材料;试验材料在露地环境下的开花时间总体表现出偏正态分布,而在温室环境下,开花时间各个阶段呈现出较为均衡的分布。温室与露地环境下的开花时间呈显著正相关。通过生物信息学分析最终得到的高质量SNP位点共103万个。试验材料的群体结构分析表明在系统发育树上各亚群内部分布较为集中,不同亚群之间的分布与材料的地理起源密切相关。全基因组衰减平均LD为2.3 kb,表明在116份白菜类作物构建的群体内存在较为频繁的重组和突变。对不同条件下的开花时间进行全基因组关联分析,用复合模型检测到54个(P4)强关联峰值信号点,一般模型检测到87个(P5)。通过进一步分析强关联信号点的连锁不平衡(linkage disequilibrium,LD)区段,得到存在强连锁关系(r20.33)的峰值信号点共33个(温室环境下27个,露地环境下19个)。其中,在温室与露地环境下的共定位位点13个。根据33个关联候选位点,再通过白菜与同源物种拟南芥的基因共线性关系以及基因功能注释分析筛选出白菜类作物开花时间相关的候选基因14个,其中温室与露地环境下共定位候选基因3个(FUL、PHYB和FPF1)。在露地条件下定位到开花关键基因FT1。【结论】不同条件下开花时间的相关性分析表明,遗传效应在开花早晚中起着决定性作用。全基因组关联分析共鉴定出33个与开花时间相关的显著关联信号。通过连锁不平衡分析、白菜与同源物种拟南芥的基因共线性关系以及基因功能注释分析初步鉴定出14个白菜类作物开花时间相关的候选基因。     

不同氮处理下水稻剑叶叶宽的全基因组关联分析

【目的】探索水稻剑叶叶宽调控及其对氮肥响应的遗传机理,为氮高效水稻新品种的培育提供新的种质资源和基因标记。【方法】以134份水稻地方种质资源为关联分析材料,通过基因组重测序发掘获得了3 356 591个分布于全基因组的高密度SNP位点(single nucleotide polymorphism,SNP)。在大田栽培条件下,以施氮量为主区,品种为裂区的设计,设计低氮(不施氮肥,N0)、正常氮(纯氮96 kg·hm~(-2),N1)和高氮(纯氮192 kg·hm~(-2),N2)3种氮肥处理。于水稻成熟期分别调查水稻剑叶叶宽在低、中、高3种氮肥处理下的表现及响应,结合EMMAX软件计算群体亲缘关系矩阵和EIGENSOFT软件分析群体结构,采用纳入亲缘关系矩阵及群体结构的混合线性模型开展全基因组关联分析。【结果】水稻剑叶叶宽在N0、N1、N2 3种氮肥处理下均呈正态分布,并表现丰富的变异。剑叶叶宽受品种差异及氮水平的影响,且与施氮水平呈显著正相关。在N0、N1、N2氮处理下共检测到14个与剑叶叶宽显著相关的SNP位点。其中,低氮处理下检测到的SNP位点的最小等位基因频率均大于0.46,表明此类SNP在关联群体中广泛存在;中氮和高氮水平下检测到的SNP位点的最小等位基因频率均较小,是一类较为稀有的SNP位点。位于第12染色体上的一个SNP(chr12:15 066 507)位点在正常氮及高氮处理下均被检测到,在高氮处理下还检测到的另一显著性位点,其候选区间内包含一候选基因LOC_Os12g25660,该基因与业已报道的叶宽性状相关基因Os BR6ox同属于细胞色素P450家族。根据不同氮处理下剑叶叶宽的响应,鉴定出20与低氮响应有关的SNP位点,8个位点与高氮响应有关。其中与高氮响应的显著性位点中,位于第1染色体显著性峰候选区间包含业已克隆的与氮素利用相关的基因OsATG7。【结论】通过全基因组关联分析共检测到42个与剑叶叶宽及其在不同氮处理下叶宽响应相关的显著性关联位点。     

小麦苗期铅耐受性的全基因组关联分析

[目的]随着工业化的推进,重金属尤其是铅对耕地的污染已成为世界性问题.小麦作为主要粮食作物,其健康生产对保障粮食安全意义重大,筛选铅耐受性强和铅低积累小麦品种、挖掘相关调控基因或QTL区间,为耐铅种质创新和揭示小麦铅耐受性遗传机制奠定基础.[方法]采用140 mg·kg-1的硝酸铅溶液对102份小麦品种(系)进行苗期胁...     

栽培稻芽期耐低温全基因组关联分析

【目的】水稻是重要的粮食作物,芽期是水稻生长发育过程中最脆弱的时期,直播稻遭遇冷害时发芽率大幅降低,减产严重。深入了解耐冷性的遗传机制,为培育芽期强耐受性水稻品种奠定基础。【方法】以世界范围内14个国家代表性的238份水稻种质资源为试验材料,于2021和2022年在沈阳开展表型鉴定试验,统计不同水稻品种在人工气候培养箱15℃低温条件下第1—10天的发芽率和相对发芽率,利用R语言绘制5—10 d的频率直方图,通过表型丰富度Hill值选择宜作关联分析的天数,将发芽率和相对发芽率表型数据与重测序数据相结合,进行基于混合线性模型MLM（QK）的全基因组关联分析,并对所获得的SNP位点进行耐冷候选基因的预测。【结果】发芽率频数分布直方图和表型丰富度计算结果显示第8天发芽率多态性最好,其Hill值为0.84,高于其他几天发芽率（0.48—0.83）,可用于全基因组关联分析;主成分分析结果显示,这些水稻品种可以分为indica、aus、temperate japonica、tropical japonica和aromatic 5个亚群;2个指标进行的GWAS分析检测到3个相同的显著性SNP位点,均位于第4染色体,解释表型的11.9%—25.4%;在上下游各50 kb进行基因搜索,共发现24个相关候选基因,进一步开展LD和单倍型分析,发现LOC_Os04g24840和LOC_Os04g25140的不同单倍型耐冷性之间存在极显著差异。LOC_Os04g24840被编码区SNP分为5个单倍型,且Hap_3的耐冷性显著强于Hap_1;LOC_Os04g25140被编码区SNP分为18个单倍型,且77 bp处的氨基酸变异（S>L）存在籼粳差异。结果表明,编码糖基转移酶的基因LOC_Os04g24840和编码F-box蛋白基因LOC_Os04g25140可能与水稻芽期耐冷性密切相关。【结论】在238份水稻种质资源中共检测到3个与芽期耐冷性显著关联的SNP位点,筛选出2个与水稻芽期耐冷性相关的候选基因。