CN110689925A - 一种基于转录组测序的耐高温玉米基因挖掘的方法 - Google Patents
一种基于转录组测序的耐高温玉米基因挖掘的方法 Download PDFInfo
- Publication number
- CN110689925A CN110689925A CN201910959026.1A CN201910959026A CN110689925A CN 110689925 A CN110689925 A CN 110689925A CN 201910959026 A CN201910959026 A CN 201910959026A CN 110689925 A CN110689925 A CN 110689925A
- Authority
- CN
- China
- Prior art keywords
- corn
- genes
- gene
- temperature
- transcriptome sequencing
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B30/00—ICT specially adapted for sequence analysis involving nucleotides or amino acids
- G16B30/10—Sequence alignment; Homology search
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B40/00—ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Biophysics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Biotechnology (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Molecular Biology (AREA)
- Genetics & Genomics (AREA)
- Artificial Intelligence (AREA)
- Bioethics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Epidemiology (AREA)
- Evolutionary Computation (AREA)
- Public Health (AREA)
- Software Systems (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种基于转录组测序的耐高温玉米基因挖掘的方法,包括高温处理下玉米的基因差异表达,设定差异表达基因的筛选条件,获得得到一定数目的上调差异表达基因或者下调差异表达基因;设定特定的筛选标准,通过BlastX比对,筛选出参与玉米高温胁迫的数个unigene,并对基因的基本功能和其参与调节的信号通路进行分类。本发明以玉米内在的遗传基因作为起点,基于转录组对耐高温玉米的基因进行分析,筛选出耐高温胁迫的基因型,来培育耐高温胁迫的新玉米品种或者种质资源,达到玉米高产稳产的目的。
Description
技术领域
本发明涉及耐高温玉米基因挖掘技术领域,具体涉及一种基于转录组测序的耐高温玉米基因挖掘的方法。
背景技术
玉米是我国重要的粮、经、饲作物,在玉米生产过程中,温度是重要的生态因子,适宜的环境温度是高产所必需的先决条件之一,温度过低或过高都会严重影响玉米的生长发育。黄淮海平原是我国重要的夏玉米产区,花期是玉米对高温胁迫最为敏感的时期,主要影响穗的发育。因此,需要选育和推广耐热性较强的玉米品种,从而减少高温热害损失来实现玉米的高产稳产以及增加经济效益。
发明内容
本发明的目的在于提供一种基于转录组测序的耐高温玉米基因挖掘的方法。
为了实现本发明的上述目的,特采用以下技术方案:
本发明包括以下步骤:
步骤一:在高温环境下处理培养玉米,获得实验组和对照组组织样本;
步骤二:对实验组和对照组组织样本进行转录组测序:从玉米的相关组织中提取RNA,然后对这些组织样品单细胞所有RNA进行测序获得结果,最后整理分析转录组数据,制作包括玉米样本名称、序列的绝对路径及序列方向的样品文件,并对照参考基因组进行分析;
步骤三:转录组测序质量评价:质量评价结果可以反映本转录组测序的质量高低,若本转录组测序的质量较高,则可以进行后续的生物信息分析;
步骤四:高温处理下玉米的基因差异表达:采用GO和KEGG数据库对差异表达基因进行分析,筛选出P值小于0.05的GO和基因通路,GO注释说明基因的基本功能,KEGG说明基因所参与的代谢通路,设定差异表达基因的筛选条件,所述差异表达基因的筛选条件包括筛选指标FDR小于0.01,Log2FC大于1,获得得到一定数目的上调差异表达基因或者下调差异表达基因;
步骤五:对总的差异表达基因的KEGG和GO功能分析:对差异表达基因进行注释可进一步了解基因所参与的功能或代谢通路;
步骤六:对特异上调表达的差异基因的KEGG和GO注释分析:可进一步了解特异上调表达的差异基因所参与的功能或代谢通路;
步骤七:结果分析与玉米耐高温相关基因的挖掘:以Log2FC大于2,差异表达基因长度大于600bp,覆盖度大于85%,同源度大于95%作为筛选条件通过BlastX比对,筛选出参与玉米高温胁迫的数个unigene,并对基因的基本功能和其参与调节的信号通路进行分类,以这些差异表达基因的log10RP-KM作热图分析,筛选出RPKM值大于240的玉米基因判断该基因在玉米中都发挥耐高温的作用。
作为优选,所述玉米组织样本优选的组织是玉米花粉/叶片。
作为优选,所述高温环境为至少40-50度,高温处理于玉米吐丝期前约7天开始进行,通过塑料薄膜覆盖升温,大棚两侧各预留出20%缝隙,用于气体交换,吐丝8天后处理结束后对提取玉米相关组织进行测定。
作为优选,步骤三可将总的clean Reads比对到参考基因组上,若比对率Q30满足标准,样品间的GC含量误差小于1%,则说明转录本测序的质量高则进入下一步骤。
作为优选,分析结果产生后还包括对差异表达基因的qRT-PCR验证:选取几个差异表达基因进行qRT-PCR验证,若DEGs的表达模式与转录组测序结果一致,则说明此次实验测序结果正确。
本发明的有益效果在于:
本发明以玉米内在的遗传基因作为起点,基于转录组对耐高温玉米的基因进行分析,筛选出耐高温胁迫的基因型,来培育耐高温胁迫的新玉米品种或者种质资源,达到玉米高产稳产的目的。
附图说明
图1为本发明基于转录组测序的耐高温玉米基因挖掘的方法的流程结构图;
具体实施方式
下面结合附图对本发明作进一步说明:
如图1所示,将本发明所述的方法的一个实施例应用到对一个耐高温玉米品种良硕88花粉/叶片的耐高温基因的挖掘中。主要步骤如下:
步骤一:在高温,环境正常温度处理下培养玉米品种良硕88,吐丝8天后处理结束,拆除薄膜,每种处理设置3组重复,取玉米花粉/叶片经液氮速冻后,于超低温冰箱中-80℃冷冻保存,用于后续RNA提取;
步骤二:进行转录组测序:从玉米花粉/叶片中提取RNA,然后对这些组织样品单细胞所有RNA进行测序获得结果;
步骤三:转录组测序质量评价:将总的clean Reads比对到参考基因组上,若比对率,Q30满足标准,样品间的GC含量接近,则说明转录本测序的质量高。步骤四:高温处理下玉米的基因差异表达:通过筛选指标FDR小于0.01,Log2FC大于1作为差异基因的筛选条件,获得得到一定数目的上调差异表达基因或者下调差异表达基因。
步骤五:对总的差异表达基因的KEGG和GO功能分析;
步骤六:对特异上调表达的差异基因的KEGG和GO注释分析:可进一步了解特异上调表达的差异基因所参与的功能或代谢通路;
步骤四筛选出的差异表达基因后再进行步骤五、六,根据GO和KEGG分析,筛选出P值小于0.05的GO和基因通路。GO注释主要说明基因的基本功能,KEGG主要说明基因所参与的代谢通路。
步骤七:结果分析与玉米耐高温相关基因的挖掘:以Log2FC大于2,差异表达基因长度大于600bp,覆盖度大于85%,同源度大于95%作为筛选条件,通过BlastX比对,筛选出已有文献报道参与玉米高温胁迫的数个unigene,并对基因的基本功能和其参与调节的信号通路进行分类。以这些差异表达基因的log10RP-KM作热图分析,表明不同的差异基因具有不同的表达模式。根据图表发现某些特定的基因在高温处理下的样品中显著上调表达或者下调表达,而且其RPKM很大,大于240,说明这些基因在玉米中都发挥了耐高温的作用;
步骤八:差异表达基因的qRT-PCR验证:随机选取几个差异表达基因进行qRT-PCR验证,若DEGs的表达模式与转录组测序结果一致,说明此次实验测序结果是可靠的。
本发明以玉米内在的遗传基因作为起点,基于转录组对耐高温玉米的基因进行分析,筛选出耐高温胁迫的基因型,来培育耐高温胁迫的新玉米品种或者种质资源,达到玉米高产稳产的目的。国内外研究人员利用转录组测序技术进行作物耐受各种其他胁迫的基因分析,以及利用转录组对玉米响应高温胁迫的分析,但是利用转录组测序进行耐高温玉米的基因挖掘我们是首例运用。
上述说明也并不仅限于上述举例,本发明未经描述的技术特征可以通过或采用现有技术实现,在此不再赘述;以上实施例及附图仅用于说明本发明的技术方案并非是对本发明的限制,参照优选的实施方式对本发明进行了详细说明,本技术领域的普通技术人员在本发明的实质范围内所作的任何修改、等同替换、改进等也应属于本发明的权利要求保护范围。
Claims (5)
1.一种基于转录组测序的耐高温玉米基因挖掘的方法,其特征在于:
包括以下步骤:
步骤一:在高温环境下处理培养玉米,获得实验组和对照组组织样本。
步骤二:对实验组和对照组组织样本进行转录组测序:从玉米的相关组织中提取RNA,然后对这些组织样品单细胞所有RNA进行测序获得结果,最后整理分析转录组数据,制作包括玉米样本名称、序列的绝对路径及序列方向的样品文件,并对照参考基因组进行分析。
步骤三:转录组测序质量评价:质量评价结果可以反映本转录组测序的质量高低,若本转录组测序的质量较高,则可以进行后续的生物信息分析;
步骤四:高温处理下玉米的基因差异表达:采用GO和KEGG数据库对差异表达基因进行分析,筛选出P值小于0.05的GO和基因通路,GO注释说明基因的基本功能,KEGG说明基因所参与的代谢通路,设定差异表达基因的筛选条件,所述差异表达基因的筛选条件包括筛选指标FDR小于0.01,Log2FC大于1,获得得到一定数目的上调差异表达基因或者下调差异表达基因;
步骤五:对总的差异表达基因的KEGG和GO功能分析:对差异表达基因进行注释可进一步了解基因所参与的功能或代谢通路;
步骤六:对特异上调表达的差异基因的KEGG和GO注释分析:可进一步了解特异上调表达的差异基因所参与的功能或代谢通路;
步骤七:结果分析与玉米耐高温相关基因的挖掘:以Log2FC大于2,差异表达基因长度大于600bp,覆盖度大于85%,同源度大于95%作为筛选条件通过BlastX比对,筛选出参与玉米高温胁迫的数个unigene,并对基因的基本功能和其参与调节的信号通路进行分类,以这些差异表达基因的log10RP-KM作热图分析,筛选出RPKM值大于240的玉米基因判断该基因在玉米中都发挥耐高温的作用。
2.根据权利要求1所述的一种基于转录组测序的耐高温玉米基因挖掘的方法,其特征在于:所述玉米组织样本优选的组织是玉米花粉/叶片。
3.根据权利要求1所述的一种基于转录组测序的耐高温玉米基因挖掘的方法,其特征在于:所述高温环境为至少40-50度,高温处理于玉米吐丝期前约7天开始进行,通过塑料薄膜覆盖升温,大棚两侧各预留出20%缝隙,用于气体交换,吐丝8天后处理结束后对提取玉米相关组织进行测定。
4.根据权利要求1所述的一种基于转录组测序的耐高温玉米基因挖掘的方法,其特征在于:步骤三可将总的clean Reads比对到参考基因组上,若比对率Q30满足标准,样品间的GC含量误差小于1%,则说明转录本测序的质量高则进入下一步骤。
5.根据权利要求1所述的一种基于转录组测序的耐高温玉米基因挖掘的方法,其特征在于:分析结果产生后还包括对差异表达基因的qRT-PCR验证:选取几个差异表达基因进行qRT-PCR验证,若DEGs的表达模式与转录组测序结果一致,则说明此次实验测序结果正确。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910959026.1A CN110689925A (zh) | 2019-10-10 | 2019-10-10 | 一种基于转录组测序的耐高温玉米基因挖掘的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910959026.1A CN110689925A (zh) | 2019-10-10 | 2019-10-10 | 一种基于转录组测序的耐高温玉米基因挖掘的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110689925A true CN110689925A (zh) | 2020-01-14 |
Family
ID=69111998
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910959026.1A Pending CN110689925A (zh) | 2019-10-10 | 2019-10-10 | 一种基于转录组测序的耐高温玉米基因挖掘的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110689925A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113337633A (zh) * | 2021-07-07 | 2021-09-03 | 广西壮族自治区农业科学院 | 间作玉米下花生叶片基因差异表达的比较转录组分析方法 |
CN114015756A (zh) * | 2021-11-04 | 2022-02-08 | 浙江省园林植物与花卉研究所(浙江省萧山棉麻研究所) | 凤丹耐涝差异表达基因及基于转录组测序挖掘的方法 |
CN116179568A (zh) * | 2023-04-23 | 2023-05-30 | 河南大学三亚研究院 | 抑制ZmP5H1和ZmP5H2基因的表达在提高玉米耐热性中的应用 |
CN116515862A (zh) * | 2023-03-30 | 2023-08-01 | 中国热带农业科学院热带作物品种资源研究所 | 一种铝胁迫条件下狗牙根差异表达基因及其挖掘和分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104131107A (zh) * | 2014-08-12 | 2014-11-05 | 江苏省农业科学院 | 高通量筛选玉米内参基因的方法 |
CN109979538A (zh) * | 2019-03-28 | 2019-07-05 | 广州基迪奥生物科技有限公司 | 一种基于10x单细胞转录组测序数据的分析方法 |
CN110055306A (zh) * | 2019-05-16 | 2019-07-26 | 河南省农业科学院粮食作物研究所 | 一种基于转录组测序挖掘玉米耐低氮基因的方法 |
-
2019
- 2019-10-10 CN CN201910959026.1A patent/CN110689925A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104131107A (zh) * | 2014-08-12 | 2014-11-05 | 江苏省农业科学院 | 高通量筛选玉米内参基因的方法 |
CN109979538A (zh) * | 2019-03-28 | 2019-07-05 | 广州基迪奥生物科技有限公司 | 一种基于10x单细胞转录组测序数据的分析方法 |
CN110055306A (zh) * | 2019-05-16 | 2019-07-26 | 河南省农业科学院粮食作物研究所 | 一种基于转录组测序挖掘玉米耐低氮基因的方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113337633A (zh) * | 2021-07-07 | 2021-09-03 | 广西壮族自治区农业科学院 | 间作玉米下花生叶片基因差异表达的比较转录组分析方法 |
CN114015756A (zh) * | 2021-11-04 | 2022-02-08 | 浙江省园林植物与花卉研究所(浙江省萧山棉麻研究所) | 凤丹耐涝差异表达基因及基于转录组测序挖掘的方法 |
CN116515862A (zh) * | 2023-03-30 | 2023-08-01 | 中国热带农业科学院热带作物品种资源研究所 | 一种铝胁迫条件下狗牙根差异表达基因及其挖掘和分析方法 |
CN116179568A (zh) * | 2023-04-23 | 2023-05-30 | 河南大学三亚研究院 | 抑制ZmP5H1和ZmP5H2基因的表达在提高玉米耐热性中的应用 |
CN116179568B (zh) * | 2023-04-23 | 2023-06-23 | 河南大学三亚研究院 | 抑制ZmP5H1和ZmP5H2基因的表达在提高玉米耐热性中的应用 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Gossypium genomics: trends, scope, and utilization for cotton improvement | |
CN110689925A (zh) | 一种基于转录组测序的耐高温玉米基因挖掘的方法 | |
Li et al. | Validation and comparison of reference genes for qPCR normalization of celery (Apium graveolens) at different development stages | |
Wang et al. | Identification of circular RNAs in kiwifruit and their species-specific response to bacterial canker pathogen invasion | |
Li et al. | Wheat breeding history reveals synergistic selection of pleiotropic genomic sites for plant architecture and grain yield | |
Chen et al. | Continuous salt stress-induced long non-coding RNAs and DNA methylation patterns in soybean roots | |
Li et al. | Genomic insights into historical improvement of heterotic groups during modern hybrid maize breeding | |
Sun et al. | Transcriptome analysis of stem development in the tumourous stem mustard Brassica juncea var. tumida Tsen et Lee by RNA sequencing | |
Shi et al. | Identification of candidate genes associated with cell wall digestibility and eQTL (expression quantitative trait loci) analysis in a Flint× Flint maize recombinant inbred line population | |
CN111349713B (zh) | 一套用于检测雪茄烟资源的kasp引物组及其应用 | |
CN110055306A (zh) | 一种基于转录组测序挖掘玉米耐低氮基因的方法 | |
Wu et al. | Genomic footprints of sorghum domestication and breeding selection for multiple end uses | |
Miculan et al. | A forward genetics approach integrating genome‐wide association study and expression quantitative trait locus mapping to dissect leaf development in maize (Zea mays) | |
CN108517368B (zh) | 利用上位性解析毛白杨LncRNA Pto-CRTG及其靶基因Pto-CAD5互作关系的方法及系统 | |
Hu et al. | Resequencing of 388 cassava accessions identifies valuable loci and selection for variation in heterozygosity | |
Busconi et al. | Epigenetic variability among saffron crocus (Crocus sativus L.) accessions characterized by different phenotypes | |
Quan et al. | Genome-wide association studies reveal the coordinated regulatory networks underlying photosynthesis and wood formation in Populus | |
Li et al. | Genomic insights into the genetic basis of cotton breeding in China | |
Han et al. | Genomic insights into genetic improvement of upland cotton in the world’s largest growing region | |
Chen et al. | Genetic diversity of Prunus sibirica L. superior accessions based on the SSR markers developed using restriction-site associated DNA sequencing | |
Yang et al. | Comparative transcriptome analysis of fiber and nonfiber tissues to identify the genes preferentially expressed in fiber development in Gossypium hirsutum | |
Wei et al. | Time-sequential detection of quantitative trait loci and candidate genes underlying the dynamic growth of Salix suchowensis | |
Bhandari et al. | Genetic architecture of fresh-market tomato yield | |
Ge et al. | Genome-wide selection and introgression of Chinese rice varieties during breeding | |
Tonnessen et al. | Shared cis-regulatory architecture identified across defense response genes is associated with broad-spectrum quantitative resistance in rice |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200114 |
|
RJ01 | Rejection of invention patent application after publication |