CN114410816A - 一种适用木薯抗病研究的内参基因的筛选方法及应用 - Google Patents
一种适用木薯抗病研究的内参基因的筛选方法及应用 Download PDFInfo
- Publication number
- CN114410816A CN114410816A CN202210006997.6A CN202210006997A CN114410816A CN 114410816 A CN114410816 A CN 114410816A CN 202210006997 A CN202210006997 A CN 202210006997A CN 114410816 A CN114410816 A CN 114410816A
- Authority
- CN
- China
- Prior art keywords
- cassava
- reference gene
- manes
- genes
- different
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6888—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms
- C12Q1/6895—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for detection or identification of organisms for plants, fungi or algae
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6844—Nucleic acid amplification reactions
- C12Q1/6851—Quantitative amplification
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/13—Plant traits
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/158—Expression markers
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/166—Oligonucleotides used as internal standards, controls or normalisation probes
Landscapes
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Organic Chemistry (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Engineering & Computer Science (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- Biotechnology (AREA)
- Molecular Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Immunology (AREA)
- Biophysics (AREA)
- Physics & Mathematics (AREA)
- Genetics & Genomics (AREA)
- Biochemistry (AREA)
- Microbiology (AREA)
- General Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Botany (AREA)
- Mycology (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明提供一种适用木薯抗病研究的稳定内参基因的筛选方法及应用,包括:1)选取木薯候选内参基因,进行内参引物设计;2)木薯组培苗进行外源激素诱导处理、鞭毛蛋白处理,及不同品种的不同组织部位取样,RNA提取,合成cDNA样本;3)对候选内参基因引物采用cDNA样本进行qRT‑PCR反应,采用geNorm法、NormFinder法和BestKeeper法和ΔCT法,对候选内参基因在外源激素诱导处理、鞭毛蛋白处理后和不同木薯组织部位的表达稳定性分析,筛选出特定条件下木薯内参基因Manes.05G133600表达最稳定,可用作实时荧光定量PCR的稳定内参,为木薯的抗病分子生物学研究奠定方法学基础。
Description
技术领域
本发明涉及植物基因工程技术领域,特别涉及一种适用木薯抗病研究的稳定内参基因的筛选方法及应用。
背景技术
实时荧光定量PCR(Quantitative Real-time PCR)是一种在DNA扩增反应中,以荧光化学物质检测每次聚合酶链式反应(PCR)循环后产物总量的方法。近年来,由于实时荧光定量PCR具有灵敏性、高效性和准确性的特点,被广泛用于分子生物学研究中分析不同基因的表达模式,成为研究基因表达的通用方法,然而实时荧光定量PCR的准确性会受到样品生理状态、样品RNA、反转录时cDNA的合成效率等因素的影响,因而需要引入内参基因(reference gene) 进行均一化用以消除样本之间的差异,从而获得目的基因的真正特异性表达。
在研究基因表达水平时,常用管家基因(house keeper)作为内参基因。管家基因又称看家基因、持家基因,是指不受外部环境影响,在所有条件下均能稳定、大量表达的一类基因,其产物是维持细胞基本生命活动所需的。然而,大量研究结果表明,并没有表达绝对稳定的基因,任何一种管家基因的所谓恒定表达都只是在一定类型的细胞或实验因素作用下“有范围”的恒定。使用不合适的内参基因对目的基因的表达量分析可能会产生极大的影响,一方面使基因表达的微小差异难以发现,另一方面出现错误甚至相反的结论,因此,为了得到更可靠的实验结果,实验中选择较为合适的内参基因进行校正显得尤为重要,稳定的内参是实时荧光定量PCR分析准确性的先决条件,而不同植物对于内参基因的筛选鉴定存在较大的差异。
木薯(Manihot esculeuta Crantz)是热带地区重要的粮食作物,同时也是一种重要的、很有发展潜力的再生能源作物。木薯在生产上的病害,严重影响了木薯产业的健康发展,不断阐明木薯抗病机理、挖掘抗性基因资源是实现木薯产业可持续发展的有效途径。而现有对于木薯相关的内参基因的筛选鉴定的研究则未见报道,针对木薯各器官的不同处理的表达变化分析,筛选鉴定适用于木薯抗病研究的稳定内参基因,为进一步开展木薯的抗病分子生物学研究奠定方法学基础具有重要的实际意义。
发明内容
鉴以此,本发明提出一种适用木薯抗病研究的稳定内参基因的筛选方法及应用,本发明筛选了13个木薯候选内参基因,并采用运用GeNorm、 Norm-Finder、BestKeeper和ΔCT分析方法,对候选的13个内参基因在不同基因型木薯不同组织、不同激素以及细菌鞭毛蛋白Flg22处理下的表达稳定性分析,发现基因Manes.05G133600在全部处理条件下是最稳定的内参,可用作实时荧光定量PCR的稳定的内参。
本发明的技术方案是这样实现的:
一种适用木薯抗病研究的内参基因的筛选方法,包括如下步骤:
S1、选取不同品种、不同组织中的木薯转录组原始数据集中的基因,计算 RPKM值作为不同基因的表达量;并根据其表达量,对其基因表达量的变异系数CV值进行排序,筛选出CV值低于0.3的基因,作为候选内参基因;并对其候选内参基因进行实时荧光定量PCR内参基因引物设计;
S2、将木薯组培苗分别进行外源激素诱导处理、鞭毛蛋白处理后取样,以及对长势均一的不同品种木薯组培苗的不同组织部位的进行取样,分别作为实验材料,进行RNA提取,合成cDNA样本;
S3、对候选内参基因引物分别采用cDNA样本进行qRT-PCR反应,检测表达量,根据其对应Ct值判断基因的表达丰度;并分别采用geNorm法、NormFinder 法和BestKeeper法和ΔCT法,对候选内参基因在外源激素诱导处理、鞭毛蛋白处理后和不同木薯组织部位的表达稳定性分析,根据分析结果,对候选内参基因进行综合稳定性排序,筛选出稳定表达的木薯内参基因。
进一步说明,步骤S2中,所述外源激素诱导处理包括采用茉莉酸、水杨酸和乙烯利的三种外源激素诱导处理。
进一步说明,所述外源激素诱导处理,具体为:对木薯组培苗分别采用质量浓度为100μmol/L的茉莉酸、水杨酸和乙烯利进行整株喷雾接种,并在不同时间点对其叶片进行取样。
进一步说明,步骤S2中,所述鞭毛蛋白处理是采用质量浓度为0.8~1.2μmol/L 的细菌鞭毛蛋白Flg22进行处理。
进一步说明,所述鞭毛蛋白处理,具体为,取木薯组培苗的幼嫩叶片,制备原生质体,采用1μmol/L的细菌鞭毛蛋白Flg22进行接种处理,并在不同时间点对其叶片进行取样。
进一步说明,所述木薯原生质体的制备方法,包括如下步骤:
(1)将木薯叶片在甘露醇中切成细条,转入酶解液中,26℃40rpm黑暗条件下酶解16h;酶解液为:1.6%纤维素酶,0.8%离析酶,9%甘露醇,0.02 mol/L KCl,0.2mol/L MES,0.1%BSA,0.01mol/L CaCl2;
(2)采用W5溶液润洗过的纱布过滤酶解产物,弃滤渣,将所得液体室温 100g,离心3min,收集离心管底部的原生质体;W5溶液为:0.45%NaCl,0.92% CaCl2,2.5mmol/L KCl,0.1mol/L MES;
(3)再次加入10mL W5溶液悬浮原生质体,离心力100g,离心3min,重复洗涤2次,即得到木薯原生质体。
进一步说明,步骤S2中,所述木薯组培苗的不同组织部位是经过MS培养基60d的木薯根、茎、叶的组织部位。
一种根据上述适用木薯抗病研究的稳定内参基因的筛选方法筛选获得的内参基因,包括Manes.05G133600、Manes.02G183100、Manes.14G008100、 Manes.05G144500、Manes.02G137500和EF1α。
更优选的,所述内参基因为Manes.05G133600。
进一步说明,Manes.05G133600、Manes.02G183100、Manes.14G008100、Manes.05G144500、Manes.02G137500和EF1α作为内参基因,用于木薯组织发育及病害胁迫的抗病性分子生物学研究中。
与现有技术相比,本发明的有益效果是:
本发明运用生物信息学预测,结合组织发育及病害胁迫下表达情况,筛选了候选内参基因的基础上,分别对检测其候选内参基因在不同不品种,不同组织部位(根、茎、叶),以及经水杨酸(SA)、茉莉酸(JA)、乙烯利(ETH) 三种激素处理、细菌鞭毛蛋白Flg22处理下的表达情况,通过采用geNorm法、 NormFinder法和BestKeeper法和ΔCT法,对其候选内参基因在木薯各器官及各处理中的表达变化进行比较分析,同时采用实时定量PCR对不同内参基因下目的基因的表达稳定性进行分析,筛选出特定条件下表达稳定的内参基因,为进一步开展木薯的抗病分子生物学研究奠定方法学基础。
经分析结果表明,基因Manes.05G133600在全部处理特定条件下是最稳定的内参,该基因的鉴定为木薯抗病分子机理研究提供了稳定的内参,其在病原菌下的表达稳定;Manes.02G183100、Manes.14G008100、Manes.05G144500、 Manes.02G137500和EF1α具有良好的表达稳定性,可用作一般实时荧光定量 PCR的内参。
附图说明
图1为本发明实施例的候选内参基因引物扩增检测的电泳图;其中,图中的第1条带是Marker,后面2-14分别对应1-13个内参引物的RT-PCR结果;
图2为本发明实施例的候选内参基因对应的木薯RNA的检测的电泳图;其中,第1条带是DNA Marker,后面2-14分别对应图1中提取的RNA,说明RNA 提取质量好。
图3为本发明实施例的13个候选内参基因在不同样本中的Ct值分布图;
图4为本发明实施例的13个候选内参基因的表达稳定性geNorm分析图;
图5为本发明实施例的13个候选内参基因的表达稳定性NormFinder分析图;
图6为本发明实施例的13个候选内参基因的表达稳定性BestKeeper分析图;
图7为本发明实施例的13个候选内参基因的表达稳定性ΔCT法分析图;
图8为本发明实施例的13个候选内参基因的稳定性综合分析柱状图;
图9为本发明实施例的MeRBOHA/MeRBOHD基因在不同内参基因不同处理下的表达模式。
具体实施方式
为了更好理解本发明技术内容,下面提供具体实施例,对本发明做进一步的说明。
本发明实施例所用的实验方法如无特殊说明,均为常规方法。
本发明实施例所用的材料、试剂等,如无特殊说明,均可从商业途径得到。
实施例-适用木薯抗病研究的内参基因的筛选方法
我们运用生物信息学预测筛选了7个新的木薯内参,结合在组织发育及病害胁迫下表达相对较稳定的6个木薯内参,检测其在阿根廷7号(Arg7)和华南8号(SC8)两个不品种,根、茎、叶组织,水杨酸(SA)、茉莉酸(JA)、乙烯利(ETH)三种激素处理,以及细菌鞭毛蛋白Flg22处理下的表达情况,并利用geNorm法、NormFinder法和BestKeeper法和ΔCT法,对13个内参基因在木薯各器官及各处理中的表达变化进行比较分析,并采用实时定量PCR对不同内参基因下目的基因的表达稳定性进行分析,筛选出特定条件下表达相对最稳定的内参基因。
实施例1-木薯内参基因的筛选和引物设计
木薯转录组原始数据集及基因组数据下载于NCBI高通量测序数据库SRA(Sequence Read Archive),包括转录组原始数据集SSR1298996(W14块根), SSR1298998(W14叶片),SSR1299000(KU50叶片),SSR129001、SSR129002 及SSR129003(KU50块根不同发育阶段),SSR129005、SSR1299006、SSR1299007 (Arg7块根不同发育阶段),SSR1299009(Arg7叶片),SSR1302047(W14 块根)。
为了探究木薯不同品种不同组织中稳定表达的内参基因,采用RPKM值 (ReadsPer Kb per Million reads)衡量11份木薯转录组原始数据集中的所有基因的表达量,其计算公式如下:(C值代表唯一匹配到某基因的reads 数量;N值代表唯一匹配到所有转录本的reads数量;L值代表匹配到的某基因 CDS碱基数量)。
本发明通过对不同实验中获得的RPKM值的计算,根据11份转录组数据所得的RPKM计算,得出木薯基因组中每个基因表达量的平均值和标准差 (Standard deviation,SD);同时,根据基因表达量的标准差与平均数的比值,得到每个基因表达量的变异系数(Coefficient of Variance,CV),采用变异系数评价基因表达稳定性,当某个基因表达量的CV值低于0.3时,可以认为该基因在不同品种不同组织中稳定表达。本发明对其基因表达量的变异系数CV值进行排序,筛选出7个CV值低于0.3的基因;并同时结合现有在组织发育及病害胁迫下表达相对较稳定的6个木薯内参基因:Ubq10(Manes.07G019300)、Tub(Manes.08G061700)、EF1a(Manes.15G054700)等,作为候选内参基因,其信息如表1;并对其候选内参基因进行实时荧光定量PCR内参基因引物设计,其引物序列信息如表2;
对这13个基因在病原菌以及抗病相关激素处理下的表达稳定性进行检测,表1。
表1候选参考基因的信息
表2 13个候选参考基因的引物和扩增
实施例2-木薯材料处理与RNA提取、cDNA合成
(1)实验材料获取
本发明采用木薯阿根廷7号(Arg7)和华南8号(SC8)的组培苗作为原始实验材料,组培苗采用MS培养基于25℃培养室培养2个月,光周期为16h光 /8h暗。
将木薯组培苗分别进行外源激素诱导处理、鞭毛蛋白处理后取样,以及对长势均一的不同品种木薯组培苗的不同组织部位的进行取样,作为实验材料;具体如下:
a.外源激素诱导处理采用喷雾法:将组培60d的Arg7木薯苗,分别采用 100μmol/L的茉莉酸(JA)、水杨酸(SA)及乙烯利(ETH)进行整株喷雾接种,每个处理9株,分别在接种前0min与接种后15min,30min,60min,120min 共5个时间点进行混合采集叶片,每个时间点取3片叶。
b.鞭毛蛋白的处理:取Arg7木薯苗幼嫩叶片,制备原生质体,采用1μmol/L 的细菌鞭毛蛋白Flg22处理,于接种前0min,接种后5min,10min,20min,30min 共5个时间点进行取样。
c.取长势均一的SC8和Arg7木薯苗各3株,对其根、茎、叶进行混合采样, 样品立即放入液氮中冻存,随后置于-80℃冰箱冷冻保存备用。
(2)RNA提取与cDNA合成
木薯原生质体的制备方法:
a.将木薯叶片在甘露醇中切成细条(1mm)后转入酶解液(1.6%纤维素酶, 0.8%离析酶,9%甘露醇,0.02mol/L KCl,0.2mol/L MES,0.1%BSA,0.01mol/L CaCl2)中,26℃,40rpm黑暗条件下酶解16h。
b.用W5溶液(0.45%NaCl,0.92%CaCl2,2.5mmol/L KCl,0.1mol/L MES) 润洗过的纱布过滤酶解产物,弃滤渣,将所得液体室温100g,离心3min(上升和下降速率均为3),轻柔去除废液,收集离心管底部的原生质体。
c.再次加入10mL W5溶液悬浮原生质体,离心力100g,离心3min,重复洗涤2次,即得木薯原生质体,可用于后续RNA的提取。
各样品参照植物总RNA提取试剂盒(TIANGEN)说明书进行总RNA的分离,为除去RNA中残留的基因组DNA,用Dnase I(Ambion)对3mg总RNA进行消化。取消化后的样品进行cDNA第一链的合成,操作参照RevertAid First Stand cDNA synthesis(Thermo)说明书进行。
实施例3-RNA质量检测、引物特异性和扩增效率分析
RNA的质量是保证Q-PCR准确性的关键因素,提取的RNA样品用Nano drop 2000(Thermo Scientific)进行了测定,其RNA的A260/A280比值均在2.0 左右,电泳结果(图2)显示有完整的28S和18S两条带,说明RNA质量和纯度均符合要求。
对设计的内参基因引物(表2)进行RT-PCR反应,测定引物的特异性。电泳结果(图1)显示不同的内参基因的引物只有一条特异性条带,该条带与预期条带大小一致,说明本实验所用引物具有高度的特异性。
引物的扩增效率(E)由公式E=(10-1/slope)计算可得,在进行qRT-PCR之前,对每对引物的qRT-PCR的扩增效率采用了标准曲线法进行了计算。用于计算引物扩增效率的qRT-PCR的cDNA模板来自于各种处理后不同时间点的cDNA等比例混合样本作为样本池,并按照梯度稀释。结果表明:不同引物的标准曲线的决定系数(R2)均在0.99以上,引物的扩增效率在94.32%~111.77%之间(表2)。 13对引物的扩增效率都达到了94%以上,均符合qRT-PCR实验对引物扩增效率 (大于90%)的要求。
实施例4-实时荧光定量PCR分析
25μL RT-PCR反应体系中如下:2.5μL 10×PCR buffer、1.5μL 2.5mmol/L MgCl2、0.5μL 10mmol/LdNTPs、0.5μL 10μmol/Lforward primer、0.5μL 10μmol/Lreverseprimer、0.2μL 5U/μLTaq和1μL模板,加无菌水补至总体积25 μL。反应条件为:94℃预变性3min;94℃变性30s,62℃退火30s,72℃延伸30s,30个循环;72℃延伸5min。RT-PCR反应结束后,进行1.2%琼脂糖凝胶电泳检测扩增产物。
Real-time PCR所用仪器为Rotor-Gene Q(QIAGEN),试剂为SYBR Premix ExTaqTM(TaKaRa)。根据引物Tm值按照TaKaRa试剂说明书,Real-time PCR 扩增程序为:95℃30s;40个循环95℃5s,55℃30s,72℃30s;72℃延伸10min。每个样品进行3次独立的生物学重复。
对13对候选内参基因引物分别采用不同处理后的不同时间点的cDNA样本进行了qRT-PCR反应,检测其表达量,得到对应的Ct值。Ct值与基因的表达丰度呈负相关,即Ct值越低,基因的表达丰度越高。根据其对应Ct值判断基因的表达丰度,13个内参基因在所有的cDNA样本中均有一定丰度的表达,Ct值分布在15-27之间(图3);其中g13基因的表达丰度最高,其Ct值在检测样本中为15.20~20.47;g12、g4和g11的表达丰度次之,Ct值为16.34~24.22;而g9 的表达丰度较低,Ct值为20.52~27.30;其他内参基因的Ct值在检测样本中为18.67~26.97;并进行ΔCT分析内参基因稳定性,其结果如图7所示。
实施例5-木薯内参基因的稳定性分析
分别采用geNorm法、NormFinder法和BestKeeper法和ΔCT法,对候选内参基因在外源激素诱导处理、鞭毛蛋白处理后和不同木薯组织部位的表达稳定性分析。
(1)geNorm法分析内参基因稳定性
geNorm软件研究样品中每个基因相对量比率的变化。基因稳定性的量化值 (M)是由多个内参基因在一组给定样本中的几何平均值和其中任何一个内参基因和所有其他内参基因的平均成对变化值来确定。M值得分越低的内参基因其表达量越稳定,而M值得分越高的内参基因的表达量变化越大,越不稳定。图 4表明,以三种(SA、JA、ETH)激素、Flg22处理、不同组织部位为样品的内参基因表达稳定性最高的是g7(Manes.S058400)和g8(Manes.02G137500),最不稳定的g12(Manes.08G061700),但是全部内参的M值均低于1.5。
(2)NormFinder法分析内参基因稳定性
NormFinder用于基因稳定性分析的算法软件。同geNorm软件相似, NormFinder是通过稳定值(S)的大小对候选内参基因稳定性进行排序,即S值越高,基因表达越不稳定,S值越低,基因表达越稳定。在所有样品中,13对内参基因引物的S值排序为图5。其中,g2(Manes.05G133600)的S值最低为0.304,是表达最稳定的基因;Tubulin(g12,Manes.08G061700)的S值最大(1.022),最不稳定,这与geNorm软件分析结果一致。
(3)BestKeeper法分析内参基因稳定性
BestKeeper软件是通过比较内参基因的Cq值得标准差(SD)与调节系数标准差SD[±x-fold]评判各候选内参基因表达稳定性。SD越小,调节系数标准差越小,基因的稳定性就越好。采用BestKeeper软件对木薯13个内参候选基因的表达稳定性进行分析的结果显示(图6),Ubq10(g11,Manes.07G019300)为最稳定的内参,这与geNorm和NormFinder软件分析结果相反;Manes.S058400为最不稳定的内参。
实施例6-综合分析内参基因稳定性
根据上述四种基因表达稳定性分析方法的分析结果,对13个木薯候选内参基因进行综合稳定性排序,结果如图8,可以看出,g2(Manes.05G133600)基因是最稳定的内参基因,g12(Manes.08G061700)是最不稳定的内参基因。
实施例7-目的基因的表达分析
ROS爆发是植物应对外界逆境环境的典型反应,受到病原菌浸染后 RBOHA、RBOHD会上调表达。为验证本发明所候选的内参基因在均一化其它基因中的作用,选择了ROS爆发中的RBOHA(Manes.14G042600.1)和RBOHD (Manes.09G172500.1)基因分别进行了检测,其结果如图9所示,从目标基因的表达趋势变化来看,采用不同的内参基因均一化,目的基因表达趋势基本一致,其中,g2基因作为内参所获得的qPCR结果更灵敏,在不同处理条件下都能有效的反应目的基因的变化;从不同时间点的稳定性来看g2作为内参目的基因表达量更稳定。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
序列表
<110> 海南大学
<120> 一种适用木薯抗病研究的稳定内参基因的筛选方法及应用
<160> 26
<170> SIPOSequenceListing 1.0
<210> 1
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 1
cttgaatgtg aagcccggtg 20
<210> 2
<211> 22
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 2
caacatcctc tccagtgaca gt 22
<210> 3
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 3
cggtggagtt gaagagcgga 20
<210> 4
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 4
gctgcgacac ctttccatcc 20
<210> 5
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 5
accccactca actcatcgcc 20
<210> 6
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 6
cacctccacg tccaccttca 20
<210> 7
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 7
gtatccctct tcctaatgtc 20
<210> 8
<211> 18
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 8
gaagagagta gcctgatc 18
<210> 9
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 9
ccttctttcc ccaccgagct 20
<210> 10
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 10
cggtcttgtc caatggcagc 20
<210> 11
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 11
ccaggctatg aaccaccaga 20
<210> 12
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 12
ggcggctctt ctccactata 20
<210> 13
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 13
gcaccatgcg agaaacaccc 20
<210> 14
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 14
agataactcc cggtcgctct 20
<210> 15
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 15
aatgcgctcc tacaacaagc 20
<210> 16
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 16
gatcatccgt agcagcctct 20
<210> 17
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 17
attatgcaag cgggacaaac 20
<210> 18
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 18
actccaccgt acatccttgc 20
<210> 19
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 19
tggtcagcac atttgttcgt 20
<210> 20
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 20
agcagacccc gtcattgtag 20
<210> 21
<211> 20
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 21
tgcatctcgt tctccgattg 20
<210> 22
<211> 21
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 22
gcgaagatca gtcgttgttg g 21
<210> 23
<211> 22
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 23
tgccatgttc cgtggaaaga tg 22
<210> 24
<211> 25
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 24
cccctaggtg gaatgtcaca gacac 25
<210> 25
<211> 23
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 25
tccactgggt agattcgctg tta 23
<210> 26
<211> 22
<212> DNA
<213> 人工序列(Artificial Sequence)
<400> 26
tagctccact tgggtccttc tt 22
Claims (10)
1.一种适用木薯抗病研究的内参基因的筛选方法,其特征在于:包括如下步骤:
S1、选取不同品种、不同组织中的木薯转录组原始数据集中的基因,计算RPKM值作为不同基因的表达量;并根据其表达量,对其基因表达量的变异系数CV值进行排序,筛选出CV值低于0.3的基因,作为候选内参基因;并对其候选内参基因进行实时荧光定量PCR内参基因引物设计;
S2、将木薯组培苗分别进行外源激素诱导处理、鞭毛蛋白处理后取样,以及对长势均一的不同品种木薯组培苗的不同组织部位的进行取样,分别作为实验材料,进行RNA提取,合成cDNA样本;
S3、对候选内参基因引物分别采用cDNA样本进行qRT-PCR反应,检测表达量,根据其对应Ct值判断基因的表达丰度;并分别采用geNorm法、NormFinder法和BestKeeper法和ΔCT法,对候选内参基因在外源激素诱导处理、鞭毛蛋白处理后和不同木薯组织部位的表达稳定性分析,根据分析结果,对候选内参基因进行综合稳定性排序,筛选出稳定表达的木薯内参基因。
2.如权利要求1所述的一种适用木薯抗病研究的内参基因的筛选方法,其特征在于:步骤S2中,所述外源激素诱导处理包括采用茉莉酸、水杨酸和乙烯利的三种外源激素诱导处理。
3.如权利要求2所述的一种适用木薯抗病研究的稳定内参基因的筛选方法,其特征在于:所述外源激素诱导处理,具体为:对木薯组培苗分别采用质量浓度为100μmol/L的茉莉酸、水杨酸和乙烯利进行整株喷雾接种,并在不同时间点对其叶片进行取样。
4.如权利要求1所述的一种适用木薯抗病研究的稳定内参基因的筛选方法,其特征在于:步骤S2中,所述鞭毛蛋白处理是采用质量浓度为0.8~1.2μmol/L的细菌鞭毛蛋白Flg22进行处理。
5.如权利要求4所述的一种适用木薯抗病研究的稳定内参基因的筛选方法,其特征在于:所述鞭毛蛋白处理,具体为,取木薯组培苗的幼嫩叶片,制备原生质体,采用1μmol/L的细菌鞭毛蛋白Flg22进行接种处理,并在不同时间点对其叶片进行取样。
6.如权利要求5所述的一种适用木薯抗病研究的稳定内参基因的筛选方法,其特征在于:所述木薯原生质体的制备方法,包括如下步骤:
(1)将木薯叶片在甘露醇中切成细条,转入酶解液中,26℃40rpm黑暗条件下酶解16h;酶解液为:1.6%纤维素酶,0.8%离析酶,9%甘露醇,0.02mol/L KCl,0.2mol/L MES,0.1%BSA,0.01mol/L CaCl2;
(2)采用W5溶液润洗过的纱布过滤酶解产物,弃滤渣,将所得液体室温100g,离心3min,收集离心管底部的原生质体;W5溶液为:0.45%NaCl,0.92%CaCl2,2.5mmol/L KCl,0.1mol/L MES;
(3)再次加入10mL W5溶液悬浮原生质体,离心力100g,离心3min,重复洗涤2次,即得到木薯原生质体。
7.如权利要求1所述的一种适用木薯抗病研究的稳定内参基因的筛选方法,其特征在于:步骤S2中,所述木薯组培苗的不同组织部位是经过MS培养基60d的木薯根、茎、叶的组织部位。
8.一种如权利要求1~7中任意一项所述的适用木薯抗病研究的稳定内参基因的筛选方法筛选获得的内参基因,其特征在于:所述内参基因包括Manes.05G133600、Manes.02G183100、Manes.14G008100、Manes.05G144500、Manes.02G137500和EF1α。
9.如权利要求8中所述的内参基因,其特征在于:所述内参基因为Manes.05G133600。
10.如权利要求8中所述的内参基因的应用,其特征在于:所述Manes.05G133600、Manes.02G183100、Manes.14G008100、Manes.05G144500、Manes.02G137500和EF1α作为内参基因,用于木薯组织发育及病害胁迫的抗病性分子生物学研究中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210006997.6A CN114410816B (zh) | 2022-01-05 | 2022-01-05 | 一种适用木薯抗病研究的内参基因的筛选方法及应用 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210006997.6A CN114410816B (zh) | 2022-01-05 | 2022-01-05 | 一种适用木薯抗病研究的内参基因的筛选方法及应用 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114410816A true CN114410816A (zh) | 2022-04-29 |
CN114410816B CN114410816B (zh) | 2023-09-19 |
Family
ID=81272042
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210006997.6A Active CN114410816B (zh) | 2022-01-05 | 2022-01-05 | 一种适用木薯抗病研究的内参基因的筛选方法及应用 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114410816B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012142299A2 (en) * | 2011-04-13 | 2012-10-18 | Donald Danforth Plant Science Center | Improved cassava |
CN109872777A (zh) * | 2019-03-14 | 2019-06-11 | 江苏省中国科学院植物研究所 | 海滨木槿实时荧光定量pcr内参基因的筛选方法 |
CN112575010A (zh) * | 2020-12-14 | 2021-03-30 | 云南农业大学 | 山药不同组织荧光定量的内参基因及其引物与应用 |
-
2022
- 2022-01-05 CN CN202210006997.6A patent/CN114410816B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2012142299A2 (en) * | 2011-04-13 | 2012-10-18 | Donald Danforth Plant Science Center | Improved cassava |
CN109872777A (zh) * | 2019-03-14 | 2019-06-11 | 江苏省中国科学院植物研究所 | 海滨木槿实时荧光定量pcr内参基因的筛选方法 |
CN112575010A (zh) * | 2020-12-14 | 2021-03-30 | 云南农业大学 | 山药不同组织荧光定量的内参基因及其引物与应用 |
Non-Patent Citations (2)
Title |
---|
ISABEL MORENO ET AL.: ""Reference genes for reliable potyvirus quantitation in cassava and analysis of Cassava brown streak virus load in host varieties"", vol. 177, pages 49 - 54 * |
MEIZHEN HU ET AL.: ""Validation of Reference Genes for Relative Quantitative Gene Expression Studies in Cassava (Manihot esculenta Crantz) by Using Quantitative Real-Time PCR"", vol. 7, pages 1 - 12 * |
Also Published As
Publication number | Publication date |
---|---|
CN114410816B (zh) | 2023-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Validation and comparison of reference genes for qPCR normalization of celery (Apium graveolens) at different development stages | |
Liu et al. | Genome‐wide development of microRNA‐based SSR markers in Medicago truncatula with their transferability analysis and utilization in related legume species | |
Ma et al. | Genome-wide identification and phylogenetic analysis of the ERF gene family in melon | |
CN114717355A (zh) | 一种西瓜全基因组SNP-Panel | |
Guo et al. | Identification and characterization of miRNAome in tobacco (Nicotiana tabacum) by deep sequencing combined with microarray | |
Sreedasyam et al. | JGI Plant Gene Atlas: an updateable transcriptome resource to improve structural annotations and functional gene descriptions across the plant kingdom | |
Lin et al. | Generation and analysis of a large-scale expressed sequence tag database from a full-length enriched cDNA library of developing leaves of Gossypium hirsutum L | |
CN113774065A (zh) | 美国白蛾不同性别成虫荧光定量内参基因及其引物和应用 | |
CN113046467A (zh) | 一组与小麦条锈病抗性显著关联的snp位点及其在遗传育种中的应用 | |
CN116926234A (zh) | 与大豆籽粒油分含量相关的snp分子标记及其应用 | |
CN114736979B (zh) | 与西瓜全缘叶形基因ClLL紧密连锁的分子标记及应用 | |
CN114410816B (zh) | 一种适用木薯抗病研究的内参基因的筛选方法及应用 | |
CN108085409B (zh) | 不同组织中杉木内参基因的筛选方法和筛选基因作为内参基因的应用 | |
CN114752702B (zh) | 一种与油菜钙含量性状QTL紧密连锁的分子标记BnCa-2C2及其应用 | |
CN108424975B (zh) | 基于转录组序列开发的葛根ssr标记引物对及其应用 | |
CN108823207B (zh) | 一种苎麻的Bn-miR43及其应用 | |
CN117230240B (zh) | 与大豆籽粒油分含量相关的InDel位点、分子标记、引物及其应用 | |
Prasad et al. | Advancement in molecular tools of plant population genetics | |
CN116970734B (zh) | 一种棉花多心室控制基因GaMV连锁的SNP位点及其应用 | |
CN114854893B (zh) | 一种与谷子抽穗期性状关联的SNPs标记及其鉴定方法 | |
CN111733276B (zh) | 耐盐基因及其应用 | |
CN115011591B (zh) | 一种用于抑制薇甘菊生长发育的RNAi分子制备方法 | |
CN117265176B (zh) | 与大豆籽粒油分含量相关的snp位点、分子标记及应用 | |
CN108728443B (zh) | 一种苎麻的Bn-miR6及其应用 | |
CN108611334A (zh) | 桑树类糖原合成酶激酶基因及其检测和应用 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |