flue-cured tobacco lines andhybrids carryingthe introgressed nicotiana rustica region,wz.crop sci.,2015,55:79-86.)；shi 等利用近等基因系进行分子标记分析，发现转育至烤烟的wz渐渗片段位于19号染色体，片段大小约65mb，并开发了可用于wz基因标记辅助选择的kasp和scar标记(shi r,hubert h, dexter-boone a,zeng j,kernodle s p,lewis r s.identification and validation of snp markersassociatedwith wz-mediated phytophthora nicotianae resistancein nicotiana tabacum l.molbreeding,2019,39:102；mccorkle k l,drake-stowe k,lewis r s and shew d.characterizationof phytophthora nicotianae resistance conferred by the introgressed nicotiana rustica region,wz, in flue-cured robacco.plant dis.,2018.102:309-317.)。此外，针对上述属于质量性状的烟草黑胫病抗性单基因(ph和wz)开发的标记，也有相关的授权专利报道(一种用于烟草黑胫病抗性基因ph辅助选择的分子标记及其应用,专利号:zl201410250266.7；一种分子标记辅助选择定向改良烟草黑胫病抗性的育种方法,专利号:zl201610979196.2；一种抗烟草黑胫病ph基因连锁的分子标记及其应用,专利号:zl 2016 1 0054257.4；一种与烟草黑胫病0号、1号生理小种抗性基因bs_t连锁的ssr标记及其应用,专利号:zl 2017 1 1124725.1)。
4.另一类是受微效多基因控制的数量性状抗性，抗源是两个雪茄烟品种florida301和 beinhart 1000-1，且对黑胫病菌0号和1号小种均有抗性，其中beinhart 1000-1的黑胫病抗性强于florida301。vontimitta等进行了beinhart 1000抗黑胫病qtl分析，检测到6个qtl与黑胫病抗性相关，贡献率最大的主效qtl位点位于8号连锁群(vontimitta v,lewis r s.mapping ofquantitative trait loci affecting resistanceto phytophthora nicotianae in tobacco(nicotiana tabacuml.)line beinhart-1000.mol breeding,2012,29:89-98.)。zhang等在beinhart 1000-1中也定位到效应最大qtl(qbs7)，ma等进一步分析了该qtl(定名为phn7.1)的位置和效应，通过开发 snp提高标记密度，利用一系列近等基因系(nil)将定位区间缩小到约3cm的区间，可用于分子标记辅助选择育种和候选基因克隆(zhang y,guo x,yan x,ren m,jiang c,cheng y,wenl,liu d,zhang y,sun m,fengq,yang a,chenga l.identification of stably expressed qtlforresistance to black shank disease in tobacco(nicotiana tabacum l.)line beinhart 1000-1.thecrop j.,2018,6:282-290；ma j m,heim c,humphry m,nifong j m,lewis r s.genetic analysisof phn7.1,a major qtl conferring partialresistance to phytophthora nicotianae in nicotianatabacum.mol breeding,2019,39:11.)。drake-stowe等(2017)对k346抗黑胫病qtl进行了分析，检测到5个qtl与黑胫病抗性相关，其中效应最大的qtl也是位于lg7连锁群，与前述 phn7.1应为同一个qtl(drake-stowe k,bakaher n,goepfert s,philippon b,markr,petersonpandlewisr s.multiple disease resistance loci affect soilborne disease resistance in tobacco (nicotiana tabacum).phytopathology,2017,107:1055-1061.)。
5.针对上述质量性状(单基因控制)和数量性状(微效多基因控制)两大类型抗源，研究较多的则是受单基因控制的质量性状抗性。因受单基因(ph和wz)控制的烟草黑胫病抗性很强、且易于通过育种手段转移到栽培烟草品种中进行应用。但此类黑胫病抗性主要是针对黑胫病0号生理小种的，随着大量具有黑胫病0号生理小种抗性烟草品种的推广应用，导
致目前生产上烟草黑胫病优势生理小种呈现0号到1号的转换趋势。为了培育出同时具有黑胫病0、1 号生理小种抗性且无野生烟草种和黄花烟草种基因组连锁累赘的烟草新品种，受微效多基因控制且呈水平抗性的数量性状抗源受到烟草育种工作者的重视。
6.鉴于此，本发明以高抗黑胫病0、1号生理小种的雪茄烟品种beinhart1000-1和易感黑胫病 0、1号生理小种的烤烟品种红花大金元为亲本，通过杂交、连续套袋自交，构建一个含有341 份株系的烟草重组自交系(rils_f
7:8
)为作图群体，利用数量性状连锁分析(qtl)法和单粒黑胫病菌谷创伤接种法，在烟草全基因组范围内筛选获得与雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁的共显性ssr标记，以加速分子标记辅助选择(mas)在烟草抗黑胫病0、 1号生理小种新品系选育中的精准、高效利用。


技术实现要素：

7.本发明正是为了解决上述问题缺陷，提供一种与雪茄烟beinhart1000-1抗黑胫病基因 qbs01紧密连锁的共显性ssr标记及其应用。
8.本发明采用如下技术方案实现。
9.一种与雪茄烟beinhart1000-1抗黑胫病基因qbs01，本发明所述的抗黑胫病基因qbs01 属于数量性状基因位点且具有烟草黑胫病0号和1号生理小种抗性，其黑胫病抗性的抗源是来自雪茄烟品种beinhart1000-1。
10.一种与雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁的共显性ssr标记，本发明所述的与雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁的共显性ssr标记的编号为 tmc14137和tmc14142，其pcr扩增产物核苷酸序列分别为seq id no.1和seq id no.2、 seq id no.3和seq id no.4所示。
11.本发明所述的分子标记所对应的2个位点的引物序列分别为：
12.tmc14137序列为：
13.tmc14137f：5
’‑
cggaagaaagagagatgaagtgg-3’(seq id no.5)，
14.tmc14137r：5
’‑
atcaacttttgcggtcgaga-3’(seq id no.6)；
15.tmc14142序列为：
16.tmc14142f：5
’‑
gcaataccaatgggcaaaac-3’(seq id no.7)，tmc14142r：5
’‑
tgcaacgcgtgtaccttaaa-3’(seq id no.8)。
17.上述的与雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁的共显性ssr标记的应用在于检测烟草基因组dna中是否存在抗黑胫病基因qbs01。
18.本发明上述的应用，该应用的方法为，以tmc14137序列的引物和tmc14142序列的引物分别扩增待检测烟草基因组dna，检测pcr扩增产物。
19.本发明该应用的方法为，pcr扩增产物中同时含有如seq id no.1和seq id no.3所示序列则表明该待测烟草植株含有源自雪茄烟beinhart1000-1抗黑胫病0、1号生理小种的纯合等位基因，基因型记为bsbs。
20.本发明该应用的方法为，pcr扩增产物中同时含有如seq id no.2和seq id no.4所示序列则为待测烟草植株不含源自雪茄烟beinhart1000-1抗黑胫病0、1号生理小种的纯合等位基因，基因型记为bsbs。
21.本发明该应用的方法为，pcr扩增产物中含有如seq id no.1和seq id no.2所示
序列，或含有如seq id no.3和seq id no.4所示序列，或含有如seq id no.1、seq id no.3和seqid no.4所示序列，或含有如seq id no.3、seq id no.1和seq id no.2所示序列即为含有源自雪茄烟beinhart1000-1抗黑胫病0、1号生理小种的杂合等位基因，基因型记为bsbs。
22.本发明的有益效果为，本发明提供了一种用于检测具有雪茄烟beinhart1000-1抗源的黑胫病抗性基因qbs01的分子标记，与已授权专利和文献报道的黑胫病抗性基因标记相比，本发明所提供的标记更具有显著优势：(1)抗源类型不同，已授权专利的抗源均是受单基因控制的质量性状，来源于非栽培烟草－野生烟草种(n.plumbaginifolia和n.longiflora)和黄花烟草种(n.rustica)；而本发明中所涉及的黑胫病抗源则是受微效多基因控制的数量性状，来源于栽培烟草种－雪茄烟beinhart1000-1。(2)抗性水平不同，已授权专利中所涉及的黑胫病抗病性是垂直抗性，抗性虽强但随着长时间大面积推广应用后，极易失去抗性或促使其他生理小种转变为优势小种，而导致具有含基因ph或wz失去对其他优势生理小种的黑胫病抗性；而本发明中所涉及的黑胫病抗病性则是水平抗性，其对黑胫病0号和1号生理小种均具有抗性且抗性持久。(3)连锁累赘存在与否，已授权专利中的黑胫病抗性基因均与非栽培烟草(野生烟草种和黄花烟草种)基因组片段紧密连锁，导致连锁累赘，且这些连锁累赘通常都是烟叶生产中不利性状，故此，在将已授权专利中的黑胫病抗性基因转移到栽培烟草中时，也因存在无法打破的连锁累赘而使育成的烟草黑胫病抗性品种无法在生产中推广应用；而本发明中所涉及的黑胫病抗性基因qbs01则来源于栽培烟草beinhart1000-1，因此，不存在影响其大面积推广应用的连锁累赘存在。(4)更精细的qtl定位结果，相较于文献报道的雪茄烟 beinhart1000-1黑胫病抗性qtl定位结果，本发明提供的与雪茄烟beinhart1000-1抗黑胫病基因qbs01两侧紧密连锁标记间的遗传距离更小、定位结果更精确，因此，用作烟草黑胫病抗性选育的分子标记进行辅助选择的结果也更精准。
23.下面结合附图和具体实施方式本发明做进一步解释。
附图说明
24.图1是基于烟草重组自交系群体(rils_f
7:8
；红花大金元
×
beinhart1000-1)的黑胫病抗性 qtl分析曲线图。
25.其中，上图为烟草全基因范围内的黑胫病抗性qtl定位图；下图为烟草第7号连锁群上黑胫病抗性主效qtl定位图。qtl定位软件为：winqtlcart v2.5；参数设置：定位方法为cim： composite interval mappingqtl，迭代次数为1000(permutation times＝1000)，显著性为0.01 (significance＝0.01)，步长为0.5cm(walk speed＝0.5cm)。上图中横坐标为烟草基因组中的染色体或连锁群编号；纵坐标为lod值。下图中横坐标为遗传距离(单位：厘摩cm)；纵坐标为lod值。图中的横向虚线是在0.01显著性阈值下的lod值＝7.86；lod曲线最高点为黑胫病抗性主效基因(qbs01_7)。
具体实施方式
26.本发明的第一目的在于提供一种既具有黑胫病0、1号生理小种抗性，又属于数量性状基因位点(qtl)的烟草黑胫病新抗源，其抗源来自与雪茄烟品种beinhart1000-1；第二
目的在于提供一种与雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁的共显性ssr标记；第三目的在于提供所述的与雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁的共显性ssr标记的应用。
27.本发明的第一目的是这样实现的，所述的烟草抗黑胫病基因qbs01是来源于雪茄烟品种 beinhart1000-1，以其为抗病亲本与感黑胫病的其余优良烟草品种进行杂交，使雪茄烟品种 beinhart1000-1中的抗黑胫病0、1号生理小种基因qbs01转移到待改良的烟草品种中，进而培育具有基因qbs01而抗黑胫病0、1号生理小种的烟草新品种(系)。
28.本发明的第二目的是这样实现的，所述的与雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁的共显性ssr标记的编号为tmc14137和tmc14142，其pcr扩增产物核苷酸序列分别为 seq id no.1和seq id no.2、seq id no.3和seq id no.4所示。
29.本发明的第三目的是这样实现的，所述的与雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁的共显性ssr标记在检测烟草基因组dna中是否存在来源于雪茄烟beinhart1000-1抗黑胫病基因qbs01中的应用。
30.为简捷、精准、高效的选择具有同时抗烟草黑胫病0、1号生理小种的品种，有针对性的选择含抗黑胫病基因qbs01的后代材料，本发明提供一种用于检测抗黑胫病基因qbs01的分子标记tmc14137和tmc14142，该分子标记采用数量性状连锁分析(qtl)法和单粒黑胫病菌谷创伤接种法，在烟草全基因组范围内筛选获得与雪茄烟beinhart1000-1抗黑胫病基因qbs01 基因连锁的共显性ssr标记，可用于抗黑胫病基因qbs01的辅助选择，以提高分子标记辅助选择的效率及抗黑胫病0、1号生理小种烟草品种选育的效率。
31.本发明以高抗黑胫病0、1号生理小种的雪茄烟品种beinhart1000-1和易感黑胫病0、1号生理小种的烤烟品种红花大金元为亲本，通过杂交、连续套袋自交，构建一个含有341份株系的烟草重组自交系(rils_f
7:8
)为作图群体，利用数量性状连锁分析(qtl)法和单粒黑胫病菌谷创伤接种法，在烟草全基因组范围内筛选获得与雪茄烟beinhart1000-1抗黑胫病基因 qbs01紧密连锁的共显性ssr标记，以加速分子标记辅助选择(mas)在烟草抗黑胫病0、1 号生理小种新品系选育中的精准、高效利用。
32.本发明所述的与雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁的共显性ssr标记的编号为tmc14137和tmc14142，其pcr扩增产物核苷酸序列分别为seq id no.1和seq id no.2、 seq id no.3和seq id no.4所示。
33.所述的分子标记所对应的2个位点的引物序列分别为：
34.tmc14137序列为
35.tmc14137f：5
’‑
cggaagaaagagagatgaagtgg-3’，
36.tmc14137r：5
’‑
atcaacttttgcggtcgaga-3’；
37.tmc14142序列为
38.tmc14142f：5
’‑
gcaataccaatgggcaaaac-3’，
39.tmc14142r：5
’‑
tgcaacgcgtgtaccttaaa-3’。
40.本发明所述的与雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁的共显性ssr标记的应用为所述的与雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁的共显性ssr标记在检测烟草基因组dna中是否存在抗黑胫病基因qbs01中的应用。
41.所述的与雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁的共显性ssr标记
的应用是分别以tmc14137序列的引物和tmc14142序列的引物扩增待检测烟草基因组dna，检测pcr扩增产物，如果pcr扩增产物中同时含有如seq id no.1和seq id no.3所示序列即为该烟草植株具有黑胫病抗性的纯合等位基因bsbs；如果pcr扩增产物中同时含有如seq id no.2和seqid no.4所示序列则为待测烟草植株不存在黑胫病抗性的纯合等位基因bsbs；如果pcr扩增产物中含有如seq id no.1和seq id no.2所示序列，或含有如seq id no.3和seq id no.4所示序列，或含有如seq id no.1、seq id no.3和seq id no.4所示序列，或含有如seq id no.3、 seq id no.1和seq id no.2所示序列即为待测烟草植株中含有黑胫病抗性的杂合基因bsbs。
42.下面以具体实施案例对本发明做进一步说明：
43.实施例1
44.采用数量性状连锁分析(qtl)法并结合单粒黑胫病菌谷创伤接种法，在烟草全基因组范围内筛选与雪茄烟beinhart1000-1抗黑胫病基因qbs01连锁的共显性ssr标记
45.一、实验材料以综合性状优良但易感黑胫病0、1号生理小种的烤烟品种红花大金元为母本，以高抗黑胫病0、 1号生理小种的雪茄烟品种beinhart1000-1为父本，经杂交、连续自交，获得341份株系的重组自交系(rils_f
7:8
)作为遗传作图群体。
46.二、亲本及rils_f
7:8
群体黑胫病0、1号生理小种抗性表型数据获得
47.采用单粒黑胫病菌谷创伤接种法，具体如下：
48.(a)烟苗的准备：取hd、beinhart1000-1、f1及341份rils_f
7:8
种子，按常规漂浮育苗方式育苗，50～60天后移栽至装有育苗基质的花盆中，1株/盆，8株/处理，4次重复/生理小种，置温室煅苗7天，烟苗成活后，剪叶、整苗，待接种黑胫病菌。
49.(b)接种方法：采用单粒菌谷创伤根茎部接种法。烟草疫霉由本实验室分离鉴定，其中pn 119 为0号生理小种、pn 68为1号生理小种，病菌活化2～3次后，在燕麦培养基(oa)平板上28℃恒温培养箱培养7～10天，待培养基表面长满黑胫病菌菌丝后，接种到燕麦固体培养基上28℃恒温培养箱培养14～18天，待菌丝长满燕麦粒后，掏出、分散成单粒，作为接种体。将烟苗根茎部一侧的基质掏开，用一次性注射器针头创伤，将1粒长满菌丝的燕麦粒贴于烟苗根茎部的创伤口处，再用基质填平，然后浇灌少量的蒸馏水使基质表面完全湿润，小花盆置于盛自来水的托盘中维持湿度。每批次接种黑胫病0、1号生理小种均以抗病亲本beinhart1000-1、感病品种红花大金元及其子一代(f1)同样处理作对照。接种后的烟苗放置在28℃、相对湿度75％以上、12小时光照12小时黑暗的人工气候室培养。
50.(c)病情调查及数据处理：接种后3天观察发病情况，按国家标准gb/t 23222-2008烟草病虫害分级及调查方法分别在接种后第7、14和21天进行3次调查，统计病情指数(di：diseaseindex)。利用每次调查的平均病情指数和调查时间间隔计算每个株系的病害流行曲线下面积 (audpc：area under disease progress curve)。最后，针对每个生理小种，3次调查的每个株系平均di值和audpc值作为4个处理条件(即接种后第7天、14天、21天的平均di值和 audpc值)。上述获得的341份rils_f
7:8
黑胫病0、1号生理小种抗性di值及audpc值作为表型数据，用于下一步的qtl连锁分析。
51.三、ssr标记分析
52.烟草基因组dna提取：采用常规ctab法或植物组织dna提取试剂盒均可，方法可参考已有的文献或试剂盒中的说明书。
53.pcr扩增及电泳检测：pcr扩增体系是常规的体系可参照已发表的文献，其中，本发明所提供的标记退火温度均为60℃；pcr扩增程序信息可参照相关文献；电泳检测也是采用常规的方法，可参照已发表的相关文献。
54.利用本实验室分别基于烤烟红花大金元和雪茄烟beinhart1000-1基因组信息开发的约50000个 ssr标记，对rils_f
7:8
群体的双亲(红花大金元和beinhart1000-1)和子一代(f1)进行多态性筛选，最终，筛选获得2001个多态性ssr标记。再利用筛选获得2001个多态ssr标记，对341 份rils_f
7:8
样品进行基因型分析。其次，利用遗传连锁作图软件joinmap 4.0对341份rils_f
7:8
样品的基因型数据进行连锁分析，绘制一张含有24条连锁群且均匀分布1974个ssr标记，覆盖烟草基因组长度为3213.138cm的高质量雪茄烟遗传连锁图谱，作为rils_f
7:8
群体的基因型值，用于下一步的qtl连锁分析。
55.四、雪茄烟beinhart1000-1抗黑胫病基因(qbs01)的全基因组qtl分析利用qtl定位分析软件winqtlcart v2.5对rils_f
7:8
群体的基因型数据和表型数据，对抗黑胫病基因qbs01进行全基因组qtl扫描。其中，相关参数设置为：定位方法选cim：composite interval mapping qtl，迭代次数为1000(permutation times＝1000)，显著性为0.01 (significance＝0.01)，步长为0.5cm(walk speed＝0.5cm)。最后，在全基因组范围内的 lod＝7.86条件下，位于第7号连锁群的146.50cm处定位获得烟草黑胫病0、1号生理小种抗性的1个主效qtl(暂命名为qbs01_7)。该主效qtl可解释约20.6～43.70％的表型变异率，且此时的lod值约为16.14～42.71，详见图1和表1。
56.表1源于beinhart1000-1抗黑胫病qtl(qbs01_7)信息统计
[0057][0058]
注：bs_di0_7d、bs_di0_14d、bs_di0_21d和bs_audpc_0分别表示烟草rils群体在接种烟草黑胫病0号生理小种第7天、14天、21天后的平均病情指数(di)和基于上述3个病情指数计算获得的病害流行曲线下面积(audpc)；同理，bs_di1_7d、bs_di1_14d、bs_di1_21d和bs_audpc_1分别表示烟草黑胫病1号生理小种的3个时期di值和audpc值。
[0059]
实施例2
[0060]
共显性连锁标记在rils_f
8:9
群体单株中的验证
[0061]
利用获得的与雪茄烟beinhart1000-1抗黑胫病基因qbs01两侧紧密连锁的共显性ssr标记 tmc14137和tmc14142，对苗期的rils_f
8:9
群体(红花大金元
×
beinhart1000-1)单株进行基因型分析，获得rils_f
8:9
群体各单株的基因型数据；另一方面，采用单粒黑胫病菌谷创伤接种法对各株系的黑胫病0、1号生理小种进行接种鉴定并获得各株系的黑胫病0、1号生理小种抗性表型值。最后，分析341份rils_f
8:9
群体的基因型数据与黑胫病0、1号生理小种抗性表型值，发现本发明公布的两个共显性ssr标记tmc14137和tmc14142的基因型值与表型值完全吻合，即，一致率达100％。具体的分析方法为：通过单粒黑胫病菌谷创伤接种法鉴定获得的各株系黑胫病0、1号生理小种di值和audpc值低于或等于抗病亲本beinhart1000-1时，该株系的基因型中也是同时呈现如seq id no.1(339bp)和seq id no.3(277bp)所示序列即为抗病纯合基因型bsbs；鉴定获得的各株系黑胫病0、1号生理小种di值和audpc值等于或高于感病亲本红花大金元时，该株系的基因型中也是同时呈现如id no.2(317bp)和seq idno.4(229bp)所示序列即为感病纯合基因型bsbs；而当鉴定获得的各株系黑胫病0、1号生理小种di值和audpc值介于感病亲本红花大金元与抗病亲本beinhart1000-1之间，也即与子一代(f1)相近时，该株系的基因型中也是同时呈现如seq id no.1和seq id no.2所示序列，或含有seq id no.3和seq id no.4所示序列，或含有seq id no.1、seq id no.3和seq idno.4所示序列，或含有seq id no.3、seq id no.1和seq id no.2所示序列即为抗病杂合基因型bsbs。
[0062]
以上结果表明，共显性标记tmc14137和tmc14142分别与源自雪茄烟beinhart1000-1抗黑胫病基因qbs01紧密连锁，且该两个标记位于目的基因/qtl(qbs01_7)两侧。利用上述两个共显性紧密连锁ssr标记，不仅可以精准、高效、便捷且低成本的实现烟草任何生育期的黑胫病0、 1号生理小种抗性鉴定，而且也可清晰鉴别待测植株中的黑胫病抗性基因型状态，此既提高了具有黑胫病0、1号生理小种抗性烤烟新品种选育的科学性、可预测性，又加速了育种进程。
[0063]
序列表
[0064]
seq id no.1：
[0065]
cggaagaaagagagatgaagtggatttgatgaatagcaataaaaagagtaataatgaca aaagttttttttgacagaaaataagaaaaaagttcaaaaatttgaggaaaaagataaa ctatcatactgaatcatgtcacatcatcttttatattcccaactttatatttatatatagatat atatatatatatatatatagatatatatatatatatatatatatatatatatatatatatatatataga tccttatctagatatatttatcttacatattttgtctatttaaagctccacatttatacattat ttttctctcgaccgcaaaagttgat
[0066]
seq id no.2：
[0067]
cggaagaaagagagatgaagtggatttgatgaatagcaataaaaagagtaataatgaca aaagttttttttgacagaaaataagaaaaaagttcaaaaatttgaggaaaaagataaa ctatcatactgaatcatgtcacatcatcttttatattcccaactttatatttatatatagatat atatatatatatatatatatatatatatatatatatatatatatagatccttatctagatatatttat cttacatattttgtctatttaaagctccacatttatacattatttttctctcgaccgcaaaa gttgat
[0068]
seq id no.3：
[0069]
gcaataccaatgggcaaaacttacatctatacttaaaatatgtatctttacatatattgtca cttgg
atttactctttagttttcttagttggtatacataaattatacacgagttatggatga ttatacattattttatatgaagttatatatatatatatatatatatatatatatatatatatatatat atatatatatatatatatttatatatacatatacataatagttttaatttttttaaatactagtc tttaaggtacacgcgttgca
[0070]
seq id no.4：
[0071]
gcaataccaatgggcaaaacttacatctatacttaaaatatgtatctttacatatattgtca cttggatttactctttagttttcttagttggtatacataaattatacacgagttatggatga ttatacattattttatatgaagttatatatatatatatatatatatacatatacataatagtttta atttttttaaatactagtctttaaggtacacgcgttgca
[0072]
seq id no.5：
[0073]
cggaagaaagagagatgaagtgg
[0074]
seq id no.6：
[0075]
atcaacttttgcggtcgaga
[0076]
seq id no.7：
[0077]
gcaataccaatgggcaaaac
[0078]
seq id no.8：
[0079]
tgcaacgcgtgtaccttaaa
[0080]
以上所述的仅是本发明的部分具体实施例,方案中公知的具体内容或常识在此未作过多描述。应当指出，上述实施例不以任何方式限制本发明，对于本领域的技术人员来说，凡是采用等同替换或等效变换的方式获得的技术方案均落在本发明的保护范围内。本技术要求的保护范围应当以其权利要求的内容为准，说明书中的具体实施方式等记载可以用于解释权利要求的内容。