BSA联合转录组分析发掘西瓜叶片黄化候选基因

叶片是植物重要的功能器官之一,不仅是植株进行光合作用的主要场所,也可作为重要的形态标记,应用于育种中。叶片颜色作为形态标记,不仅可用于苗期杂种的清除,亦可用于种子纯度的测定。以西瓜全生育期叶片黄化突变体纯合自交系ly104为母本(P_(1))、绿叶自交系w3为父本(P_(2)),通过杂交创制F_(1)代、F_(2)代、BC_(1)代群体。遗传分析结果表明,该突变体的叶片黄化由单隐性基因控制。采用混合分组分析(BSA)进行初定位,通过简化基因组测序(RAD)开发全基因组单核苷酸多态性(SNP)标记构建西瓜高密度遗传图谱,将西瓜叶片黄化基因定位于2号染色体13950306~15517591 bp(大小约为1.57 Mb)。以西瓜97103v2为参考基因组,该区间包含24个注释基因。对P_(1)(P1Y)、P_(2)(P2G)和F 2代群体中黄叶(F2Y)、绿叶(F2G)株系进行转录组水平分析,结果表明,目标区间内基因Cla97C02G035950、Cla97C02G036010、Cla97C02G036020、Cla97C02G036060在黄化叶片与正常绿叶材料中的表达量差异显著,可能是西瓜叶片的黄化候选基因。研究结果可为进一步解析西瓜叶片黄化基因功能和生物学特性奠定重要基础。

基于文献计量学BSA在作物育种领域的应用现状与展望

为了解国内外集群分离分析法(bulked segregant analysis,BSA)在作物育种领域的研究现状与前沿动态,客观反映各国家、机构及研究人员在该领域的影响力,本文基于文献计量学分析方法使用文献计量分析软件CiteSpace对WOS(Web of Science)数据库2000~2023年2111项和CNKI(China National Knowledge Infrastructure)数据库2003~2023年446项研究成果进行关键词共现分析、突显词分析、关键词聚类分析、聚类时间线分析及作者共被引分析。结果表明:BSA在作物育种领域应用的研究成果在国内外的发文趋势相同,国内外期刊发文量均逐年上升;在发文量国家排序中,中国排名第一、美国第二、印度第三。在CNKI数据库中华中农业大学的发文量最多,而在WOS数据库中中国农业科学院的发文量最多。国外BSA在作物育种领域应用的文章主要集中在植物科学、农学、园艺学和遗传学等学科,而国内主要集中在作物学、植物保护学、园艺学、生物学等学科;Michlmore RW、Kosambi DD和Li H这3位作者在该领域的影响力最高,而Michlmore RW、Lander ES、Li H这3位作者与其他作者有更密切的联系。国内研究的热点作物和性状分别是水稻(Oryza sativa)、大豆(Glycine max)、玉米(Zea mays L.)和抗病性、株高;国外研究的热点作物和性状分别是水稻、拟南芥(Arabidopsis thaliana)、小麦(Triticum aestivum L.)和抗病性。目前,BSA在国内集中应用于作物目标性状候选基因及作物突变体突变基因的定位和功能验证,而国外则集中应用于作物目标性状基因的精细定位和功能验证及遗传机制的解析。此外,BSA在作物育种领域应用的研究前沿分析表明,未来在该领域热点研究的对象为水稻、花生(Arachis hypogaea L.)、陆地棉(Gossypium hirsutum L.)、作物突变体和作物代谢产物。

BSA-Seq结合RNA-Seq技术挖掘大豆叶片提前黄化衰老基因

大豆产量与其生殖生长的持续时间呈正相关,延缓其开花后叶片的衰老,提高其生理功能,有利于增加植物的粒重。叶片黄化是植物衰老的显著特征之一。关于在大豆鼓粒后期叶片黄化的研究鲜见报道。本研究对鼓粒后期叶片提前黄化突变体ly和野生型ofc杂交组合进行遗传分析,结果表明,大豆鼓粒后期叶片提前黄化性状受单隐性核基因控制。通过BSA-Seq在19号染色体得到一个2.23Mb的初定位区间,经开发标记图位克隆后将区间缩短至1.75 Mb,区间内有219个基因,再结合RNA-Seq分析,得到了区间内12个候选基因,其中有4个SNP变异基因和8个差异表达基因。本研究结果为大豆鼓粒后期叶片黄化衰老基因的克隆奠定了基础。

基于BSA-seq技术挖掘番茄耐裂果候选基因

[目的]明确番茄耐裂基因在染色体上的位置。[方法]以耐裂果番茄14FP5和易裂果番茄14FP26为试材,将二者杂交得到F 1代,F 1自交得到F 2,分别构建2个极端性状的DNA混合池,采用BSA-seq方法进行定位。[结果]利用BSA-seq对F 2群体的2个极端混池进行基因组重测序分析,将番茄耐裂基因定位在5号染色体上,位于34140000~38120000 pb,区间大小为3.98 Mb,候选区间内共注释到38个基因,其中非同义突变基因4个。[结论]该研究结果可为番茄耐裂基因的精细定位提供数据支撑,为番茄耐裂分子育种奠定基础。

基于BSA-seq技术对乌拉尔图小麦条锈病抗病基因的初步定位

【目的】发掘和克隆小麦条锈病抗病基因,为小麦条锈病抗病育种提供遗传基因资源。【方法】以乌拉尔图小麦为材料,通过接种鉴定筛选具有条锈病抗性的材料,利用正向遗传学和混合分组测序分析法(BSA-seq)在乌拉尔图小麦材料中鉴定新的条锈病抗病基因位点。【结果】通过接种鉴定筛选到6个具有小麦条锈病抗性的乌拉尔图小麦材料,通过构建杂交群体结合BSA-seq技术在其中1个抗病材料TuW1中鉴定出1个新的条锈病抗病基因位点,定位于7AS染色体的0~35.2 Mb区间。【结论】通过正向遗传学结合BSA-seq技术在乌拉尔图小麦中鉴定到1个新的条锈病抗病基因位点,为小麦条锈病抗病育种提供了基因资源。

基于BSA全基因组重测序的甘蓝型油菜有限花序突变体候选基因的鉴定

油菜是我国第一大油料作物,机械化生产是我国油菜产业发展的必然趋势。油菜无限花序导致花期长和角果成熟不一致,是影响机械化收获的关键因素。因此,开展甘蓝型油菜有限花序性状的基因挖掘对油菜的遗传改良、培育适合油菜机械化收获的新品种、突破油菜机械化生产的瓶颈具有重大意义。本研究以一个自然变异产生的甘蓝型油菜有限花序突变体Y69为研究对象,通过与中双11杂交构建F2分离群体,在群体中挑选有限花序和无限花序的单株各20株构建混合池,对混合池和亲本开展20×和10×覆盖度的全基因组重测序。根据测序结果,分别筛选出277679个SNP和302625个InDel多态性位点,用于有限花序性状的全基因组定位。利用△SNP-in-dex方法关联分析,共筛选到892个多态性标记位点,涉及683个基因,分布于甘蓝型油菜基因组A09、A10和C093条染色体上的6个关联区间,其中C09染色体上的关联区间峰值最高。通过拟南芥基因组同源序列比对,结合基因功能注释以及序列差异分析,候选基因预测位于A09和C09染色体上的关联区间内包含6个候选基因,分别为Bn⁃aA09g34410D、BnaA09g37880D、BnaA09g38520D、BnaC09g40470D、BnaC09g40480D和BnaC09g49710D,其中,Bn⁃aA09g34410D、BnaA09g37880D和BnaC09g49710D根据基因注释参与花发育和开花时间的控制参与花发育和开花时间调控,且位于C09染色体上的3个基因存在等位变异,推测为控制该有限花序性状变异的关键基因,本研究结果为油菜有限花序后续的基因克隆和功能验证奠定了基础。

基于BSA-ELM模型的建筑项目施工成本预测研究

为降低建筑施工项目的管理成本,提出一种基于BSA-ELM的建筑项目施工成本预测模型。首先介绍ELM神经网络算法和BSA算法;其次利用BSA算法的优势对ELM模型进行优化,从而形成一种基于BSA-ELM的预测模型;最后对基于BSA-ELM建筑项目施工成本预测模型进行测试,将其与传统的预测模型进行对比。结果证明,BSA-ELM预测模型的性能更好,预测精度更高。

利用BSA-Seq技术初步鉴定棉花芽黄候选基因

芽黄性状可作为棉花育种的标记性状,在棉花杂交选育中可减轻棉花杂交制种和鉴定过程,还可以作为叶绿素等生理代谢研究的理想材料。研究棉花芽黄突变体具有重要的理论和实践意义。本研究以芽黄突变体115-23与绿叶正常材料JSH1527为亲本,构建F2分离群体,鉴定后代芽黄性状,选择30个芽黄单株和30个正常绿叶单株构建混合池,对4个样本(2个亲本池和2个混合池)开展全基因组重测序,采用SNP-index进行关联分析,确定芽黄基因的候选区域,筛选芽黄候选基因。结果表明,在2个亲本之间共获得840,612个多态性SNP,利用Δ(SNP-index)方法,在95%的置信区间内将候选区域定位到染色体D02区间,包含93个候选基因,其中4个基因功能注释与叶绿体合成等相关,可能是调控棉花叶色的候选基因。研究结果可为棉花芽黄分子机制及芽黄相关基因的克隆奠定理论基础。

基于BSA-seq技术对豌豆花色基因的精细定位

BSA-seq技术在挖掘农艺性状相关的新基因中已被广泛应用,随着豌豆首个参考基因组问世,将BSA-seq技术结合豌豆基因组的基因定位策略势在必行。本研究利用紫花亲本G0004562、白花亲本G0002930以及F2群体,通过BSA-seq技术对豌豆花色基因进行初步定位,获得31.42Mb定位区间,再通过设计InDel分子标记分析进一步缩小定位区间,最终将目标基因定位在包含19个基因的0.99 Mb区间内,通过基因注释信息推测出Psat6g060480.1为豌豆花色候选基因。本研究结果验证了BSA-seq技术快速高效定位豌豆花色基因的可行性,为利用该技术挖掘豌豆其他重要农艺性状相关基因奠定了基础。

基于BSA分析技术对陆地棉抗旱基因的鉴定

【目的】研究抗旱相关的候选基因,为棉花抗旱性鉴定与评价和抗旱基因的挖掘提供理论依据。【方法】以陆地棉石远321和奎85-174为亲本构建重组自交系中选取的20份极端抗旱和旱敏感材料,将筛选的极端材料分别构建极端抗旱池、极端敏旱池,与亲本进行BSA测序分析,以陆地棉基因组为参考,通过BSA测序分析,挖掘抗旱相关的候选基因。【结果】关联到3个相关的候选区域,候选区域的总长度为6.12 Mb,其中共包含86个基因。结合GO注释、KEGG通路和前期筛选出的5个关键指标,共筛选出5个候选基因。【结论】5个基因在叶片中的表达量显著高于根和茎,并且在抗旱性强材料中的表达量均高于抗旱性弱的材料,且总体呈先上升后下降的趋势。

基于极端混合池(BSA)全基因组重测序的羽衣甘蓝白色叶基因定位

筛选羽衣甘蓝白色叶形成相关的关键基因和位点以及心叶颜色性状的遗传规律,以羽衣甘蓝红叶亲本WR、白叶亲本WB及其构建的F2分离群体为试验材料,从F2群体中挑选红叶和白叶植株各20株,分别构建2个DNA混合池,对子代混合池和亲本分别开展50×和20×覆盖深度的全基因组重测序,定位白叶性状关联区间,并根据基因注释信息预测候选基因。结果表明,重测序获得3 987 718个单核酸多态性(SNP)标记,获得位于第2、3和9染色体的8个显著关联区间,根据基因注释和功能分析,筛选到8个候选基因(Bol030253、Bol029431、Bol012077、Bol007709、Bol030318、Bol030235、Bol030286和Bol005195)。将8个基因确定为羽衣甘蓝白叶的候选基因,可能在羽衣甘蓝叶片颜色形成过程中起着重要作用。

基于BSA-seq方法挖掘大豆再生相关候选基因

转基因生物育种技术可以定向改良大豆品种,为大豆定向育种提供一种新思路。为寻找与大豆再生相关的基因,探索大豆再生规律、提高遗传转化效率,本研究利用再生能力强材料东农50、再生能力弱材料Keburi及其子代RIL群体的200份材料进行大豆器官发生试验,比较不同基因型之间再生能力的差异,筛选出极端材料各20份,后通过BSA-seq(bulked segregant analysis sequencing)技术对大豆再生候选基因进行初步定位,共获得88.04 G的clean data,平均测序深度为20.03×,定位到2 Mb区间,利用GO等数据库对区间内基因进行富集分析,差异表达基因主要富集在纤维素微纤维组织、植物型细胞壁组织或生物发生等20个条目中,其中被显著富集的植物型细胞壁组织或生物发生条目中共有6个基因,对6个基因进行组织表达量分析,在丛生芽伸长期间表达水平较高,说明其在大豆再生过程中发挥作用,可能为影响大豆再生的关键基因。本研究为大豆再生机制研究提供了重要的候选基因信息与必要的材料基础。

利用BSA-Seq筛选甘蓝型油菜菌核病抗性候选基因

油菜菌核病居油菜3大病害之首,它所引起的茎腐烂是导致油菜减产的重要原因之一,因此筛选出调控菌核病的主效基因,对利用分子辅助育种手段提高油菜产量具有重要意义.该研究以甘蓝型油菜高抗材料21Y490为父本与高感材料21Y689杂交,在F2分离群体终花期截取油菜茎秆进行核盘菌室内接种,培养3 d后测量菌斑大小鉴定菌核病抗性.筛选抗性极端表型个体构建DNA混池,与亲本一起重测序提取差异SNP,InDel信息,利用BSA-seq技术、ED算法和Δ-index算法获得两个显著的关联区域,分别位于C09号染色体的17.8~21.2 Mb和31.6~39.1 Mb区间,区间内共有399个基因.最后,利用拟南芥和甘蓝型油菜基因组注释信息对关联候选区间对应的基因功能进行分析,共筛选出与油菜菌核病抗性相关的8个候选基因(BnaC05g50540D,BnaC09g33900D,BnaA04g09130D,BnaA01g29170D,BnaC09g15640D,BnaC09g35310D,BnaA03g09220D,BnaCnng15460D),通过候选基因同源分析、代谢通路分析及共表达基因分析发现,这些基因主要参与植物激素生长调控及免疫系统的信号转导,为下一步抗性基因的克隆及功能研究奠定了基础.

基于BSA-seq结合连锁分析发掘大豆荚粒性状QTLs及候选基因

豆荚和籽粒是大豆重要的产量形成器官,同时,荚粒性状作为大豆重要外观品质,还直接影响其商品性。鉴于此,基于大豆重组自交系群体(recombinant inbred line,RIL)在多种环境条件下的荚粒性状(荚长、荚宽、荚重、粒长、粒宽和粒重),筛选极端家系构建2个混池——大荚大粒池(large pool,LP)和小荚小粒池(small pool,SP);经混合群体分离测序(bulked segregant analysis sequencing,BSA-seq)寻找控制荚粒性状的遗传位点,同时结合前期该群体连锁定位QTLs以及RIL亲本荚粒转录组数据,发掘多环境、多方法、多性状共定位区间以及候选基因。结果表明,构建的LP与SP混池在不同环境下的荚粒性状存在显著差异,BSA-seq结果与Williams82参考基因组的平均比对效率在98%以上,10×基因组覆盖率在94%以上;基于此,在置信度0.95时,经SNP-index、Euclidean distance(ED)和G-statistic算法分别获得27、41和51个关联区域,在置信度0.99时,分别获得5、15和16个关联区域,位于11条染色体(1、3、4、6、7、12、13、15、17、19和20号);经与前期连锁定位QTLs比较发现,在7、19和20号染色体上存在多环境、多方法、多性状共定位遗传位点,且20号染色体SNP标记ss715637629及19号染色体SNP标记ss715635705、ss715635755、ss715635844等在连锁定位与BSA-seq同时被检测到;进一步分析3种BSA-seq算法共关联区域内的候选基因,发现40个基因存在外显子非同义突变,其中15个基因在RIL群体亲本C813和KN7豆荚和籽粒表达量较高,并有3个基因随豆荚发育进程表达量增加。以上结果为开展大豆荚粒性状分子遗传改良以及遗传基础解析奠定了重要基础。

通过二维光谱分析BSA中不同微环境下的色氨酸

在研究小分子与蛋白质相互作用时,常用到模式蛋白BSA。它含有两个不同微环境下的色氨酸残基,这两个色氨酸又分别位于两个不同的结合口袋附近。由于它们的荧光发射谱重叠,所以常规的方法难以区分。为进一步判断小分子究竟结合到了哪个区域,本文尝试采用一种更为简便易行的方法:通过二维相关光谱手段,以BV浓度作为变量,检测了一系列BSA的内源荧光发射谱,并进行了同步、异步分析,获得了一些初步结果。

热处理对等离子喷涂BSAS环境障涂层性能影响

采用自制的固相烧结BSAS粉体,利用大气等离子喷涂工艺在SiC基体表面制备Si/Mullite+BSAS/BSAS三层结构环境障涂层。通过扫描电子显微镜、能谱仪和X射线衍射分析仪等研究环境障涂层在不同热处理温度下涂层的相结构、显微组织演变规律。结果表明:喷涂态BSAS涂层主要由单斜结构BAS相和非晶相组成;经过1100℃热处理后,涂层内部非晶相转变为六方结构BAS相;1200℃热处理后六方结构BAS相逐渐向单斜结构BAS相转变,在1300℃热处理后单斜结构BAS相的含量达到最大;随着热处理温度的进一步上升,出现鳞石英结构的SiO_(2)相和高钡含量的Ba3SiO5相以及副钡长石结构BAS相,1400℃热处理过程中发现硅液滴渗出的现象。

基于BSA-seq和RNA-seq挖掘水稻株高相关QTL

株高是影响水稻产量稳定性的重要因素之一,水稻株高相关QTL的定位及候选基因的挖掘,有助于深入了解株高分子调控机制,进而为培育理想株型的水稻品种奠定基础。以油占8号及广西野生稻为亲本构建的285个CSSL群体为试验对象,结合NGS与BSA-seq等方法,利用SNP和InDel两种分子标记进行关联分析,对可能与株高相关的基因组区段进行定位。结果显示,Δ(SNP-index)分子标记在Chr.7和Chr.10中分别关联到3.205和1.311 Mb大小的候选基因组区域。Δ(InDel-index)标记关联到的基因组区域大小分别为2.848和1.292 Mb,且全部包含于Δ(SNP-index)关联到的区间内。结合GO、KEGG、Uniprot、eggNOG等数据库中的功能注释、高质量多态性位点筛选以及已有的株高相关转录组数据,最终关联到Chr.7中的5个候选基因,包括功能和机理已知的OsTCP21。RT-qPCR结果显示,OsTCP21在高株和矮株中的表达差异与已有研究结果相一致,LOC_Os07g05050和LOC_Os07g02850在高株中表达量较高,LOC_Os07g04220和LOC_Os07g02770在矮株中表达量较高。这5个候选基因在水稻株高性状的调控中起重要的作用,OsTCP21是一个关键调控基因。

基于BSA-seq技术定位绿豆种皮颜色基因

种皮颜色是作物重要的驯化性状和形态标记,绿豆种皮颜色与黄酮类化合物含量有关。挖掘绿豆种皮颜色相关基因有助于开发新品种,提高绿豆种皮的应用价值。本研究以冀绿9号(黑色种皮)和资源330(黄色种皮)为亲本构建F2分离群体,采用BSA-seq方法进行定位。结果表明,SNPs和InDels关联区域交集位于4号染色体共3.26 Mb的区段,包含324个基因,其中非同义突变基因共49个,移码突变基因15个。进一步开发KASP分子标记进行精细定位,筛选出高质量KASP引物11对。最终将绿豆种皮色位点定位于4号染色体上的KA330~KA421之间,物理距离为16、302、330~18、013、421 bp(1.71 Mb)。结合转录组数据分析和qRT-PCR表达分析共筛选出6个候选基因,其中LOC106758748基因注释为MYB90,其功能为参与调控黄酮类化合物生物合成,故可能是调控绿豆种皮颜色的关键基因。本研究结果可为绿豆种皮色相关基因的克隆和在育种中的利用提供一定的理论依据。

基于BSA-seq技术的马铃薯块茎蛋白含量基因定位与分子标记开发

【目的】开发与四倍体马铃薯蛋白含量相关的分子标记,加快选育高蛋白的马铃薯品种,提高马铃薯品种竞争力。【方法】以高蛋白马铃薯品种‘大西洋’为母本、低蛋白品种‘定薯1号’为父本构建分离群体,利用混池和高通量简化基因组测序技术相结合的方法(BSA-seq)对马铃薯块茎蛋白含量进行QTL定位,在定位区间开发与马铃薯块茎蛋白含量紧密连锁的SSR标记,并使用分离群体及四倍体马铃薯品种进行筛选验证。【结果】在2号、4号染色体共定位到3个控制马铃薯块茎蛋白含量的主效位点,分别位于2号染色体18.88~21.59 Mb,区间大小为2.71 Mb;4号染色体8.3~12.84 Mb,区间大小为4.54 Mb;4号染色体65.12~66.39 Mb,区间大小为1.27 Mb。根据定位区域、基因位置及参考基因组信息,使用NR、 Tr EMBL、 KEGG、 GO、 KOG、 swissprot、 PFAM共7个功能数据库对候选基因进行功能注释,共注释到719个候选基因,注释基因显著富集在玉米素合成代谢通路,该代谢通路可阻止蛋白质降解。在2号染色体18.88~21.59 Mb区间内开发分子标记pChr2-4。用分离群体、四倍体马铃薯品种结合田间表型对pChr2-4的准确性进行鉴定,结果显示,pChr2-4在离群体、四倍体马铃薯品种中的准确性分别为73.17%和81.82%,准确度较高。【结论】通过BSA混池测序,检测到3个与马铃薯块茎蛋白含量相关的区间。在2号染色体开发与马铃薯块茎蛋白含量紧密连锁的标记pChr2-4,该标记有助于马铃薯块茎蛋白含量分子标记辅助育种。

基于BSA-seq的拟南芥辐射诱变突变体的基因定位

辐射诱变极易导致染色体结构变异,目前鲜有采用BSA-seq方法定位辐射诱变突变体基因的研究报道。为探索BSA-seq用于辐射诱变突变体的基因定位的可行性,本研究通过辐射诱变筛选得到一个拟南芥突变体,将其与亲本杂交,在F2代分离群体中根据表型筛选单株,构建两个极端表型的子代混池;将两个子代混池与亲本混池进行全基因组测序,使用MutMap、QTL-seq和GPS等多种BSA-seq定位策略对其进行定位分析。结果表明,多种BSA-seq定位策略均取得了一致的结果,在2号染色体上获得7 Mb的初步定位区间;使用IGV软件对该区间内的基因组进行可视化,发现该区间内存在25 189 bp的大片段缺失,缺失区段包含6个基因;通过SnpEff进行突变位点注释,经基因注释信息推测AT2G28610为候选基因;通过遗传学试验验证了该候选基因为突变基因。本研究结果为应用BSA-seq技术定位辐射诱变得到的突变位点提供了参考。