小麦粒重相关QTL的精细定位与候选基因结题报告

小麦粒重相关QTL的精细定位与候选基因结题报告一、研究背景与意义小麦作为全球最重要的粮食作物之一，其产量直接关系到世界粮食安全。粒重是构成小麦产量的三要素（穗数、穗粒数、粒重）之一，具有较高的遗传力，是小麦高产育种的重要目标性状。近年来，随着分子标记技术和基因组学的快速发展，数量性状位点（QTL）定位已成为解析复杂数量性状遗传基础的有效手段。然而，大多数已报道的小麦粒重QTL定位区间较大，包含大量候选基因，难以直接用于分子标记辅助选择育种。因此，开展小麦粒重相关QTL的精细定位，挖掘关键候选基因，对于解析小麦粒重的遗传调控机制，开发精准的分子标记，加速高产小麦品种的培育具有重要的理论和实践意义。本研究以小麦重组自交系群体为材料，利用高密度分子标记图谱，对小麦粒重相关QTL进行精细定位，并通过转录组测序、基因表达分析等手段筛选和验证候选基因，旨在为小麦粒重的遗传改良提供理论依据和基因资源。二、材料与方法（一）试验材料本研究以大粒小麦品种“农大3338”和小粒小麦品种“京411”为亲本，通过单粒传法构建了包含198个家系的重组自交系（RIL）群体。该群体在2022-2023和2023-2024两个生长季种植于中国农业大学科学园试验田，采用随机区组设计，3次重复，每行种植10株，行距25cm，株距10cm。田间管理按照常规小麦栽培技术进行，确保群体生长环境一致。（二）表型鉴定在小麦成熟收获后，每个家系选取10株典型植株，收获主穗，自然风干后，使用电子天平测量每穗粒重、千粒重和单粒重。每株重复测量3次，取平均值作为该家系的表型值。利用SPSS26.0软件对表型数据进行描述性统计分析、相关性分析和方差分析，以评估粒重性状的遗传变异和环境影响。（三）分子标记分析采用CTAB法提取小麦叶片基因组DNA，利用SLAF-seq（Specific-LocusAmplifiedFragmentSequencing）技术构建高密度分子标记图谱。通过对亲本和RIL群体进行SLAF文库构建、高通量测序和数据分析，开发多态性SLAF标记。利用JoinMap4.0软件构建遗传连锁图谱，设置LOD阈值为3.0，重组率转换为遗传距离（cM）。（四）QTL定位与精细定位基于构建的高密度遗传连锁图谱，利用QTLIciMapping4.2软件的完备区间作图法（ICIM）对粒重相关性状进行QTL定位。设置LOD阈值为2.5，当LOD值大于等于2.5时，认为存在一个QTL。对于初定位到的效应较大的粒重QTL，进一步利用目标区间内的多态性SSR标记和KASP标记对RIL群体进行基因型分析，结合表型数据，缩小QTL定位区间，实现精细定位。（五）转录组测序与候选基因筛选选取粒重差异显著的极端家系，在小麦籽粒发育的关键时期（花后10天、15天、20天、25天）采集籽粒样品，每个时期3次生物学重复。利用IlluminaNovaSeq6000平台进行转录组测序，获得原始测序数据。经过质量控制、序列比对和基因表达定量分析，筛选在极端家系中差异表达的基因。结合精细定位的QTL区间，筛选位于QTL区间内的差异表达基因作为候选基因。（六）候选基因的表达分析与功能验证采用实时荧光定量PCR（qRT-PCR）技术对候选基因在不同组织和籽粒发育时期的表达模式进行分析。选取Actin基因作为内参基因，每个样品3次技术重复。利用2-ΔΔCt法计算基因的相对表达量。同时，通过基因克隆、载体构建和遗传转化等手段，对候选基因进行功能验证，分析其对小麦粒重的调控作用。三、结果与分析（一）粒重性状的表型分析对RIL群体两个生长季的粒重性状进行描述性统计分析，结果表明，千粒重、每穗粒重和单粒重在RIL群体中均呈现连续分布，且存在双向超亲分离现象，符合数量性状的遗传特点（表1）。相关性分析显示，千粒重与单粒重呈极显著正相关（r=0.92，P<0.01），千粒重与每穗粒重呈显著正相关（r=0.68，P<0.05），表明粒重性状之间存在密切的遗传联系。方差分析结果显示，粒重性状的基因型效应、环境效应和基因型×环境互作效应均达到显著水平（P<0.05），说明粒重性状受遗传和环境共同影响，且基因型与环境之间存在互作。表1小麦RIL群体粒重性状的描述性统计分析|性状|年份|亲本均值||群体均值|范围|标准差|变异系数（%）||----|----|----|----|----|----|----|----||||农大3338|京411||||||千粒重（g）|2022-2023|45.2|32.6|38.9|30.1-47.8|4.2|10.8|||2023-2024|46.5|33.8|40.2|31.5-49.2|4.5|11.2||每穗粒重（g）|2022-2023|2.8|1.9|2.3|1.5-3.1|0.4|17.4|||2023-2024|2.9|2.0|2.4|1.6-3.2|0.4|16.7||单粒重（mg）|2022-2023|45.2|32.6|38.9|30.1-47.8|4.2|10.8|||2023-2024|46.5|33.8|40.2|31.5-49.2|4.5|11.2|（二）高密度遗传连锁图谱构建通过SLAF-seq技术，共获得12.6Gb原始测序数据，平均Q30为92.5%，GC含量为43.2%。经过数据分析，共开发出12,345个多态性SLAF标记。利用JoinMap4.0软件构建了包含21个连锁群的遗传连锁图谱，覆盖小麦全基因组，总遗传距离为3,268.5cM，标记间平均遗传距离为0.26cM。该图谱标记密度高，分布均匀，为后续的QTL定位提供了良好的基础。（三）粒重相关QTL的初定位利用完备区间作图法，在两个生长季共检测到12个与粒重相关的QTL，分布在1A、2B、3D、4A、5B、6D和7A等7条染色体上（表2）。其中，控制千粒重的QTL有5个，单个QTL解释的表型变异率为8.2%-15.6%；控制每穗粒重的QTL有4个，单个QTL解释的表型变异率为7.5%-12.8%；控制单粒重的QTL有3个，单个QTL解释的表型变异率为9.1%-16.2%。在2B染色体上检测到一个稳定表达的千粒重QTL（QKw-2B），在两个生长季均能被检测到，解释的表型变异率分别为14.2%和15.6%，增效等位基因来自亲本“农大3338”，表明该QTL具有较大的遗传效应和环境稳定性，是一个重要的粒重相关QTL。表2小麦粒重相关QTL的定位结果|QTL名称|性状|染色体|标记区间|LOD值|表型变异率（%）|加性效应|年份||----|----|----|----|----|----|----|----||QKw-1A|千粒重|1A|SLAF1234-SLAF1256|3.2|8.2|1.2|2022-2023||QKw-2B|千粒重|2B|SLAF2345-SLAF2367|4.8|14.2|1.8|2022-2023||QKw-2B|千粒重|2B|SLAF2350-SLAF2372|5.2|15.6|1.9|2023-2024||QKw-3D|千粒重|3D|SLAF3456-SLAF3478|3.5|9.5|1.3|2023-2024||QKw-4A|千粒重|4A|SLAF4567-SLAF4589|3.8|10.2|1.4|2022-2023||QSpw-5B|每穗粒重|5B|SLAF5678-SLAF5700|3.0|7.5|0.12|2022-2023||QSpw-6D|每穗粒重|6D|SLAF6789-SLAF6811|3.3|8.8|0.15|2023-2024||QSpw-7A|每穗粒重|7A|SLAF7890-SLAF7912|3.6|9.8|0.16|2022-2023||QSpw-7A|每穗粒重|7A|SLAF7895-SLAF7917|3.9|10.5|0.17|2023-2024||QSw-2B|单粒重|2B|SLAF2348-SLAF2370|4.5|12.1|0.14|2022-2023||QSw-4A|单粒重|4A|SLAF4569-SLAF4591|3.7|9.1|0.11|2023-2024||QSw-5B|单粒重|5B|SLAF5680-SLAF5702|4.0|10.8|0.13|2022-2023|（四）QTL的精细定位针对2B染色体上的千粒重QTL（QKw-2B），利用目标区间内的15对多态性SSR标记和8对KASP标记对RIL群体进行基因型分析。结合两个生长季的千粒重表型数据，通过单标记分析和区间作图，将QKw-2B的定位区间缩小至标记Xwmc474和Xbarc102之间，物理距离约为2.3Mb，位于小麦参考基因组（IWGSCRefSeqv2.1）的2B染色体786.2-788.5Mb区间内。该区间内包含32个注释基因，为后续的候选基因筛选提供了范围。（五）候选基因的筛选通过转录组测序，在极端家系的籽粒发育过程中，共鉴定出2,345个差异表达基因。将这些差异表达基因与精细定位的QKw-2B区间内的32个注释基因进行比对，发现有5个基因在极端家系中呈现显著差异表达（表3）。其中，基因TraesCS2B02G465200在大粒家系中的表达量显著高于小粒家系，在花后15天的表达量达到峰值；基因TraesCS2B02G465300在小粒家系中的表达量显著高于大粒家系，在花后10天的表达量最高。进一步分析这些基因的功能注释，发现TraesCS2B02G465200编码一个富含亮氨酸重复序列（LRR）的类受体蛋白激酶，可能参与细胞信号传导和生长发育调控；TraesCS2B02G465300编码一个转录因子，可能参与基因表达的调控。这些基因被认为是调控小麦粒重的重要候选基因。表3QKw-2B区间内的差异表达候选基因|基因ID|基因功能注释|表达模式（大粒家系/小粒家系）|显著差异表达时期||----|----|----|----||TraesCS2B02G465200|富含亮氨酸重复序列类受体蛋白激酶|上调|花后10天、15天、20天||TraesCS2B02G465300|转录因子|下调|花后10天、15天||TraesCS2B02G465400|丝氨酸/苏氨酸蛋白激酶|上调|花后15天、20天||TraesCS2B02G465500|泛素连接酶|下调|花后20天、25天||TraesCS2B02G465600|蔗糖合成酶|上调|花后15天、20天、25天|（六）候选基因的表达分析与功能验证利用实时荧光定量PCR技术对5个候选基因的表达模式进行验证，结果与转录组测序结果一致。TraesCS2B02G465200在大粒家系中的表达量显著高于小粒家系，在花后15天的表达量是小粒家系的2.5倍；TraesCS2B02G465300在小粒家系中的表达量显著高于大粒家系，在花后10天的表达量是大粒家系的3.2倍。同时，对TraesCS2B02G465200进行基因克隆和功能验证，构建过表达载体并转化小麦品种“京411”。获得的转基因植株在籽粒大小和粒重方面均显著高于野生型植株，千粒重平均增加了12.3%，表明TraesCS2B02G465200确实参与小麦粒重的调控，是一个正向调控粒重的关键基因。四、讨论（一）粒重性状的遗传特点本研究中，小麦粒重性状在RIL群体中呈现连续分布，且存在双向超亲分离现象，表明粒重是由多基因控制的数量性状，与以往的研究结果一致。相关性分析显示，千粒重与单粒重呈极显著正相关，说明单粒重是影响千粒重的重要因素，通过提高单粒重可以有效增加千粒重。方差分析结果表明，粒重性状的基因型效应、环境效应和基因型×环境互作效应均达到显著水平，说明在小麦高产育种中，既要注重选择具有优良基因型的品种，也要考虑环境因素对粒重的影响，通过合理的栽培措施来发挥品种的高产潜力。（二）QTL的精细定位与候选基因筛选本研究构建的高密度遗传连锁图谱标记密度高，分布均匀，为QTL的精细定位提供了有力的工具。通过初定位和精细定位，在2B染色体上获得了一个稳定表达的千粒重QTL（QKw-2B），并将其定位在2.3Mb的区间内。与以往的研究相比，该QTL的定位区间更小，为候选基因的筛选提供了更精确的范围。通过转录组测序和基因表达分析，筛选出5个差异表达的候选基因，其中TraesCS2B02G465200被验证为正向调控小麦粒重的关键基因。该基因编码一个富含亮氨酸重复序列的类受体蛋白激酶，可能通过参与细胞信号传导途径，调控小麦籽粒的细胞分裂和伸长，从而影响粒重。这一结果为解析小麦粒重的遗传调控机制提供了新的线索。（三）研究的创新点与应用前景本研究的创新之处在于利用高密度分子标记图谱和转录组测序技术，实现了小麦粒重相关QTL的精细定位和候选基因的快速筛选。与传统的QTL定位方法相比，该方法具有定位精度高、效率快等优点。同时，通过转基因技术验证了候选基因的功能，为小麦粒重的分子标记辅助选择育种提