小麦远缘杂交后代小片段易位系中发掘外源功能候选基因的方法技术

本发明专利技术提供小麦远缘杂交后代小片段易位系中发掘外源功能候选基因的方法，该方法仅需对小麦远缘杂交衍生的小片段易位系及其小麦受体和野生近缘种供体三个样品进行转录组测序，结合生物信息学分析、比较基因组学和功能标记开发验证，可实现从小麦远缘杂交衍生系中快速挖掘外源功能候选基因。本发明专利技术不依赖小麦近缘植物的参考基因组序列，实现了低成本且高效地筛选功能候选基因，为功能基因的克隆奠定基础。

A Method for Discovering Candidate Genes for Foreign Functions in Small Fragment Translocation Lines of Distant Hybrid Progenies of Wheat

【技术实现步骤摘要】
小麦远缘杂交后代小片段易位系中发掘外源功能候选基因的方法
本专利技术涉及基因组学、转录组学及生物信息学领域，具体地说，涉及小麦远缘杂交后代小片段易位系中发掘外源功能候选基因的方法。
技术介绍
从小麦进化历史和育种进程可以看出，远缘杂交是拓宽小麦遗传基础，创制新种质，实现小麦高产、稳产育种的有效途径之一。国内外学者已经实现了小麦族的所有属与小麦间的杂交，获得了大量的且具有可利用价值的小麦新种质，从这些新种质中发掘新的功能基因并加以利用，将在小麦抗性育种中发挥着举足轻重的作用。由于多数小麦近缘植物没有参考基因组序列，且在小麦背景下，对外源功能基因进行图位克隆的研究相对困难。有研究人员通过突变体筛选、抗病基因富集测序和易位染色体分拣测序组装等技术，成功从簇毛麦的6V染色体上克隆了抗白粉病基因Pm21。类似研究工作耗时长、工作量大、成本高，且相关的研究主要集中在抗病相关的基因，但与小麦复杂的数量性状相关的基因克隆还没有相关的研究报道。因此，如何有效地在小麦背景下快速筛选外源功能基因，特别是与产量等复杂性状相关的基因，是亟需解决的问题。
技术实现思路
本专利技术的目的是提供小麦远缘杂交后代小片段易位系中发掘外源功能候选基因的方法。为了实现本专利技术目的，本专利技术提供小麦远缘杂交后代小片段易位系中发掘外源功能候选基因的方法，包括以下步骤：A、小麦近缘植物的全长转录组测序：获取小麦近缘植物多个组织且包含所研究的目的性状发生前和发生时的组织，分别提取各个组织mRNA且等量混合；为了获得小麦近缘植物的全长转录组参考序列，利用第三代测序技术和平台对混合样品进行全长转录组的建库和测序。B、小麦远缘杂交后代小片段易位系和小麦受体的转录组测序：分别获取小麦远缘杂交后代小片段易位系和小麦受体植株多个组织且包含所研究的目的性状发生前和发生时的组织，分别提取各个组织mRNA且等量混合；为了对转录本的表达进行定量，利用第二代测序技术和平台对各自的混合样品进行高通量转录组的建库和测序，分别得到小麦受体的转录本序列和易位系的转录本测序数据。C、易位系中外源染色体片段中特异表达转录本分析：对小麦近缘植物的全长转录组测序数据进行生物信息学分析(生物信息学分析包括：使用IsoSeq3软件处理原始测序数据。具体步骤如下：首先，使用ccs算法进行分类，产生每个零模式波导孔中的共有序列，保留具有至少一个完整通道的序列用于后续分析；其次，使用Lima算法获得全长序列；最后，进行聚类和矫正)，获得非嵌合且全长的转录本序列，去除冗余后，与小麦参考基因组序列和注释文件进行整合作为总的参考序列；将小麦受体的转录本序列以及易位系的转录本序列分别与所述总的参考序列进行比对，并过滤，去重复，仅保留唯一比对的序列比对结果，然后统计比对到每个转录本上的序列数目，最后进行小麦受体和易位系之间转录本的序列差异表达显著性检验，根据易位系和小麦受体之间转录本的序列比对数目，使用参数log2(差异倍数)≤-4和假阳性率<0.05进行过滤(为了获得可信的以及统计上显著差异表达的转录本，必须控制表达的差异倍数以及多重检验的P-value)，获得易位系中外源染色体片段中特异表达的转录本，作为候选转录本序列。D、分子标记开发：将获得的候选转录本序列与小麦参考基因序列进行比对，根据比对结果找出候选转录本与小麦基因组上同源基因之间的差异，根据差异序列设计并开发分子标记及其检测引物；利用设计的引物在小麦近缘植物、小麦受体和易位系中进行PCR验证，确定所获得的外源染色体片段中特异表达的转录本的真实性。E、候选基因预测：对于步骤D中获得的真实候选转录本序列，与小麦参考基因序列进行比对，进行比较基因组学分析，构建易位系中外源染色体片段与小麦基因组间的比较基因组学图谱；对共线性候选区间内的基因进行功能注释，获得外源功能候选基因。可选地，步骤C中，对小麦近缘植物的全长转录组测序数据进行生物信息学分析，获得非嵌合且全长的转录本序列，利用CD-HIT-EST软件去冗余，然后与小麦参考基因组序列和注释文件进行整合作为总的参考序列。利用STAR软件将小麦受体的转录本序列以及易位系的转录本序列分别与所述总的参考进行序列比对，并过滤，去重复，仅保留唯一比对的序列比对结果。然后，用featureCounts软件统计比对到每个转录本上的序列数目。最后，用DESeq2软件进行小麦受体和易位系之间转录本的差异表达显著性检验。步骤D中，将获得的候选转录本序列通过BLAST与小麦基因序列进行比对。步骤D中，所开发的分子标记包括但不限于EST、KASP标记。在本专利技术的一个具体实施方式中，提供小麦-冰草易位系中发掘外源功能候选基因的方法，包括以下步骤：S1、冰草的全长转录组测序：获取冰草Z559根、茎、叶、幼穗、籽粒多个组织，分别提取各个组织mRNA且等量混合，构建PacBio测序文库，利用PacBiosequel平台的2个cell对混合样品进行全长转录组测序。共获得了约15G大小的测序数据量。S2、小麦受体Fukuhokumugi和易位系的转录组测序：分别获取Fukuhokumugi和易位系植株的根、茎、叶、籽粒(包括灌浆前、中、后期)多个组织，然后提取各个组织mRNA且等量混合，对Fukuhokumugi和易位系各自的混合样品构建高通量测序文库，并使用IlluminaHiSeq2500高通量测序平台进行转录组测序；S3、易位系中冰草外源染色体易位片段中特异表达的转录本分析：整合冰草全长转录组测序数据和小麦参考基因组序列作为总的参考序列，用STAR软件将Fukuhokumugi和易位系的转录本序列与所述总的参考序列进行比对，仅保留唯一比对的序列比对结果，然后用featureCounts软件分别计算Fukuhokumugi和易位系中比对到每个转录本上的序列数目，最后使用DESeq2软件进行Fukuhokumugi和易位系之间转录本的序列差异表达显著性检验，用log2(差异倍数)≤-4且假阳性率<0.05进行过滤，获得易位系中外源染色体片段中特异表达的转录本，作为候选转录本序列。共获得普冰3035中冰草外源染色体易位片段中特异表达的12个转录本(表1)。S4、分子标记开发：将获得的候选转录本序列与小麦参考基因序列进行比对，并根据比对结果找出候选转录本与小麦A/B/D基因组上同源基因之间的差异，根据差异序列设计EST标记及其检测引物；利用设计的EST标记在冰草Z559、Fukuhokumugi和易位系中进行PCR验证，分别对每个候选转录本开发了相对应的多态性分子标记(图1)，从而验证外源染色体片段中特异表达的转录本的真实性。S5、候选基因预测：对于步骤S4中获得的真实候选转录本序列，与小麦参考基因序列进行比对，进行比较基因组学分析，构建易位系中外源染色体片段与小麦A/B/D基因组间的比较基因组学图谱(图2)；对共线性候选区间内的基因进行功能注释，获得外源功能候选基因。本专利技术中，所述小麦-冰草易位系为普冰3035。步骤S1还包括对冰草全长转录组测序数据进行过滤、质控；具体方法为主要使用IsoSeq3软件处理原始测序数据，首先，使用ccs算法进行分类，产生每个零模式波导孔中的共有序列，保留具有至少一个完整通道的序列用于后续分析；其本文档来自技高网...

【技术保护点】
1.小麦远缘杂交后代小片段易位系中发掘外源功能候选基因的方法，其特征在于，包括以下步骤：A、小麦近缘植物的全长转录组测序：获取小麦近缘植物多个组织且包含所研究的目的性状发生前和发生时的组织，分别提取各个组织mRNA且等量混合；为了获得小麦近缘植物的全长转录组参考序列，利用第三代测序技术和平台对混合样品进行全长转录组的建库和测序；B、小麦远缘杂交后代小片段易位系和小麦受体的转录组测序：分别获取小麦远缘杂交后代小片段易位系和小麦受体植株多个组织且包含所研究的目的性状发生前和发生时的组织，分别提取各个组织mRNA且等量混合；为了对转录本的表达进行定量，利用第二代测序技术和平台对各自的混合样品进行高通量转录组的建库和测序，分别得到小麦受体和易位系的转录本测序数据；C、易位系中外源染色体片段中特异表达转录本分析：对小麦近缘植物的全长转录组测序数据进行生物信息学分析，获得非嵌合且全长的转录本序列，去除冗余后，与小麦参考基因组序列和注释文件进行整合作为总的参考序列；将小麦受体的转录本序列以及易位系的转录本序列分别与所述总的参考序列进行比对，并过滤，去重复，仅保留唯一比对的序列比对结果，然后统计比对到每个转录本上的序列数目，最后进行小麦受体和易位系之间转录本的序列差异表达显著性检验，根据易位系和小麦受体之间转录本的序列比对数目，使用参数log2(差异倍数)≤‑4和假阳性率<0.05进行过滤，获得易位系中外源染色体片段中特异表达的转录本，作为候选转录本序列；D、分子标记开发：将获得的候选转录本序列与小麦参考基因序列进行比对，根据比对结果找出候选转录本与小麦基因组上同源基因之间的差异，根据差异序列设计并开发分子标记及其检测引物；利用设计的引物在小麦近缘植物、小麦受体和易位系中进行PCR验证，确定所获得的外源染色体片段中特异表达的转录本的真实性；E、候选基因预测：对于步骤D中获得的真实候选转录本序列，与小麦参考基因序列进行比对，进行比较基因组学分析，构建易位系中外源染色体片段与小麦基因组间的比较基因组学图谱；对共线性候选区间内的基因进行功能注释，获得外源功能候选基因。...

【技术特征摘要】
1.小麦远缘杂交后代小片段易位系中发掘外源功能候选基因的方法，其特征在于，包括以下步骤：A、小麦近缘植物的全长转录组测序：获取小麦近缘植物多个组织且包含所研究的目的性状发生前和发生时的组织，分别提取各个组织mRNA且等量混合；为了获得小麦近缘植物的全长转录组参考序列，利用第三代测序技术和平台对混合样品进行全长转录组的建库和测序；B、小麦远缘杂交后代小片段易位系和小麦受体的转录组测序：分别获取小麦远缘杂交后代小片段易位系和小麦受体植株多个组织且包含所研究的目的性状发生前和发生时的组织，分别提取各个组织mRNA且等量混合；为了对转录本的表达进行定量，利用第二代测序技术和平台对各自的混合样品进行高通量转录组的建库和测序，分别得到小麦受体和易位系的转录本测序数据；C、易位系中外源染色体片段中特异表达转录本分析：对小麦近缘植物的全长转录组测序数据进行生物信息学分析，获得非嵌合且全长的转录本序列，去除冗余后，与小麦参考基因组序列和注释文件进行整合作为总的参考序列；将小麦受体的转录本序列以及易位系的转录本序列分别与所述总的参考序列进行比对，并过滤，去重复，仅保留唯一比对的序列比对结果，然后统计比对到每个转录本上的序列数目，最后进行小麦受体和易位系之间转录本的序列差异表达显著性检验，根据易位系和小麦受体之间转录本的序列比对数目，使用参数log2(差异倍数)≤-4和假阳性率<0.05进行过滤，获得易位系中外源染色体片段中特异表达的转录本，作为候选转录本序列；D、分子标记开发：将获得的候选转录本序列与小麦参考基因序列进行比对，根据比对结果找出候选转录本与小麦基因组上同源基因之间的差异，根据差异序列设计并开发分子标记及其检测引物；利用设计的引物在小麦近缘植物、小麦受体和易位系中进行PCR验证，确定所获得的外源染色体片段中特异表达的转录本的真实性；E、候选基因预测：对于步骤D中获得的真实候选转录本序列，与小麦参考基因序列进行比对，进行比较基因组学分析，构建易位系中外源染色体片段与小麦基因组间的比较基因组学图谱；对共线性候选区间内的基因进行功能注释，获得外源功能候选基因。2.根据权利要求1所述的方法，其特征在于，步骤C中，对小麦近缘植物的全长转录组测序数据进行生物信息学分析，获得非嵌合且全长的转录本序列，利用CD-HIT-EST软件去冗余，然后与小麦参考基因组序列和注释文件进行整合作为总的参考序列；和/或利用STAR软件将小麦受体的转录本序列以及易位系的转录本序列分别与所述总的参考进行序列比对，并过滤，去重复，仅保留唯一比对的序列比对结果；和/或然后用featureCounts软件统计比对到每个转录本上的序列数目；和/或最后用DESeq2软件进行小麦受体和易位系之间转录本的差异表达显著性检验。3.根据权利要求1所述的方法，其特征在于，步骤D中，将获得的候选转录本序列通过BLAST与小麦基因序列进行比对。4.根据权利要求1-3任一项所述的方法，其特征在于，步骤D中，所开发的分子标记包括EST、KASP标记。5.小麦-冰草易位系中发掘外源功能候选基因的方法，其特征在于，包括以下步骤：S1、冰草的全长转录组测序：获取冰草Z559根、茎、叶、幼穗、籽粒多个组织，分别提取各个组织mRNA且等量混合，构建PacBio测序文库，利用PacBios...

【专利技术属性】
技术研发人员：李立会，周升辉，张锦鹏，韩海明，刘伟华，鲁玉清，杨欣明，李秀全，
申请(专利权)人：中国农业科学院作物科学研究所，
类型：发明
国别省市：北京,11