CN114350772A - 一种基于染色质开放性测序技术研究植原体病害的方法 - Google Patents
一种基于染色质开放性测序技术研究植原体病害的方法 Download PDFInfo
- Publication number
- CN114350772A CN114350772A CN202210043846.8A CN202210043846A CN114350772A CN 114350772 A CN114350772 A CN 114350772A CN 202210043846 A CN202210043846 A CN 202210043846A CN 114350772 A CN114350772 A CN 114350772A
- Authority
- CN
- China
- Prior art keywords
- phytoplasma
- sequencing
- analysis
- peak
- plant sample
- 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
- 241001468265 Candidatus Phytoplasma Species 0.000 title claims abstract description 46
- 201000010099 disease Diseases 0.000 title claims abstract description 42
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 title claims abstract description 42
- 238000012163 sequencing technique Methods 0.000 title claims abstract description 35
- 108010077544 Chromatin Proteins 0.000 title claims abstract description 29
- 210000003483 chromatin Anatomy 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 title claims abstract description 20
- 108090000623 proteins and genes Proteins 0.000 claims abstract description 23
- 108091023040 Transcription factor Proteins 0.000 claims abstract description 19
- 102000040945 Transcription factor Human genes 0.000 claims abstract description 19
- 238000012216 screening Methods 0.000 claims abstract description 15
- 238000004458 analytical method Methods 0.000 claims abstract description 13
- 210000003855 cell nucleus Anatomy 0.000 claims abstract description 9
- 239000012634 fragment Substances 0.000 claims abstract description 6
- 238000012165 high-throughput sequencing Methods 0.000 claims abstract description 6
- 108020004414 DNA Proteins 0.000 claims abstract description 4
- 108010012306 Tn5 transposase Proteins 0.000 claims abstract description 4
- 238000003766 bioinformatics method Methods 0.000 claims abstract description 4
- 210000004027 cell Anatomy 0.000 claims abstract description 4
- 238000001976 enzyme digestion Methods 0.000 claims abstract description 4
- 238000002360 preparation method Methods 0.000 claims abstract description 4
- 238000000746 purification Methods 0.000 claims abstract description 4
- 238000005336 cracking Methods 0.000 claims abstract description 3
- 241000196324 Embryophyta Species 0.000 claims description 19
- 238000010201 enrichment analysis Methods 0.000 claims description 10
- 241001629585 Paulownia fortunei Species 0.000 claims description 9
- 108700009124 Transcription Initiation Site Proteins 0.000 claims description 6
- MBABOKRGFJTBAE-UHFFFAOYSA-N methyl methanesulfonate Chemical compound COS(C)(=O)=O MBABOKRGFJTBAE-UHFFFAOYSA-N 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 238000011144 upstream manufacturing Methods 0.000 claims description 4
- 238000007405 data analysis Methods 0.000 claims description 2
- 238000000556 factor analysis Methods 0.000 claims description 2
- 230000008506 pathogenesis Effects 0.000 abstract description 5
- 238000012795 verification Methods 0.000 abstract description 2
- 244000055346 Paulownia Species 0.000 description 8
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 4
- 241001481828 Glyptocephalus cynoglossus Species 0.000 description 4
- 244000007853 Sarothamnus scoparius Species 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 239000006228 supernatant Substances 0.000 description 4
- 208000024891 symptom Diseases 0.000 description 4
- 239000003153 chemical reaction reagent Substances 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 239000008188 pellet Substances 0.000 description 3
- 239000000725 suspension Substances 0.000 description 3
- 241001430197 Mollicutes Species 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000000338 in vitro Methods 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 230000037361 pathway Effects 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 239000002023 wood Substances 0.000 description 2
- 108700039691 Genetic Promoter Regions Proteins 0.000 description 1
- GLNADSQYFUSGOU-GPTZEZBUSA-J Trypan blue Chemical compound [Na+].[Na+].[Na+].[Na+].C1=C(S([O-])(=O)=O)C=C2C=C(S([O-])(=O)=O)C(/N=N/C3=CC=C(C=C3C)C=3C=C(C(=CC=3)\N=N\C=3C(=CC4=CC(=CC(N)=C4C=3O)S([O-])(=O)=O)S([O-])(=O)=O)C)=C(O)C2=C1N GLNADSQYFUSGOU-GPTZEZBUSA-J 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000023852 carbohydrate metabolic process Effects 0.000 description 1
- 235000021256 carbohydrate metabolism Nutrition 0.000 description 1
- 210000002421 cell wall Anatomy 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000037149 energy metabolism Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 244000005700 microbiome Species 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 210000004940 nucleus Anatomy 0.000 description 1
- 230000003950 pathogenic mechanism Effects 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000019491 signal transduction Effects 0.000 description 1
Images
Landscapes
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本发明公开了一种基于染色质开放性测序技术研究植原体病害的方法，属于分子生物学领域；包括下述步骤：(1)样本处理：分别取健康植物样本与植原体侵染植物样本，裂解植物细胞获得细胞核；(2)文库制备及测序：利用步骤(1)获得的细胞核，使用Tn5转座酶酶切并纯化，回收DNA片段，构建高通量测序文库后进行测序；(3)生物信息学分析：对健康植物样本与植原体侵染植物样本的测序结果进行Peak分析，筛选植原体病害相关基因。本发明方法可为研究植原体病害发病相关转录因子功能的预测和验证奠定基础，提供借鉴。
Description
技术领域
本发明属于分子生物学技术领域，具体涉及一种基于染色质开放性测序技术研究植原体病害的方法。
背景技术
植原体(Phytoplasma)是一类没有细胞壁，不能离体培养的柔膜菌纲(Mollicutes)微生物，可诱发1,000多种植物病害，并造成严重的的生态破坏和经济损失，但目前尚未阐明植原体的具体致病机理。由于植原体不能体外培养，严重限制了植原体病害发病机理的阐明。
泡桐，树皮灰色、灰褐色或灰黑色，幼时平滑，老时纵裂；假二杈分枝，单叶，对生，叶大，卵形，全缘或有浅裂，具长柄，柄上有绒毛；喜光，较耐阴，喜温暖气候，耐寒性不强；其生长速度快，耐瘠薄，材性优良，在防风固沙、保障粮食安全和木材供应等方面起着重要的作用。
泡桐丛枝病是植原体病害的一种，制约泡桐产业发展的重要因素，丛枝病的防治的有效方法仍未建立，还需要开辟新的途径进行研究。因此，如何提供一种泡桐丛枝病的分析方法，为泡桐丛枝病及其他植原体病害的防治提供依据是本领域亟待解决的问题。
发明内容
本发明公开了一种基于染色质开放性测序技术研究植原体病害的方法，有助于了解植原体病害发生过程中染色质开放性的变化，为进一步探究植原体病害的发病机理打下基础和提供研究思路，有助于泡植原体病害的防治研究。
为了实现上述目的，本发明采用如下技术方案：
一种基于染色质开放性测序技术研究植原体病害的方法，包括下述步骤：
(1)样本处理
分别取健康植物样本与植原体侵染植物样本，裂解植物细胞获得细胞核；
(2)文库制备及测序
利用步骤(1)获得的细胞核，使用Tn5转座酶酶切并纯化，回收DNA片段，构建高通量测序文库后进行测序；
(3)生物信息学分析
对健康植物样本与植原体侵染植物样本的测序结果进行Peak分析，筛选植原体病害相关基因。
进一步地，步骤(1)中，
白花泡桐丛枝病苗经60mg·L-1甲基磺酸甲酯处理后恢复健康状态后取顶芽作为健康植物样本，并以长势一致的白花泡桐丛枝病苗顶芽作为植原体侵染植物样本。
进一步地，步骤(2)中，
构建好的高通量测序文库用Illumina HiSeq4000平台进行PE150测序。
进一步地，步骤(3)中，
数据分析包括Peak的基本分析、Peak关联基因分析及转录因子分析、GO富集分析、KEGG富集分析。
进一步地，步骤(3)中，
Peak关联基因分析：从Peak上下游5kb中，找到最近的转录起始位点对应的基因。
进一步地，步骤(3)中，
根据健康植物样本与植原体侵染植物样本Peak关联基因筛选差异Peak；
筛选标准条件为：FDR值小于0.05，并且Fold change大于0；
从差异Peak区域鉴定出植原体病害相关转录因子家族。
进一步地，步骤(3)中，
对差异Peak关联基因区域进行GO富集分析、KEGG富集分析，筛选染色质开放性增强的区域和染色质开放性减弱的区域。
进一步地，步骤(3)中，
对染色质开放性增强的区域和染色质开放性减弱的区域内peak关联基因进行转录因子的鉴定与分类，从中筛选植原体病害相关基因。
综上所述，本发明通过染色质开放性测序技术，揭示了植原体病害发病过程染色质可及性发生显著改变，并预测出与染色质开放区域互作的转录因子。为之后研究植原体病害发病相关转录因子功能的预测和验证奠定了基础，为其他植原体病害研究提供借鉴。
附图说明
附图1所示为样品对比图片；其中，
A.未使用MMS试剂处理的白花泡桐丛枝病苗，具有明显的丛枝症状；
B.使用60mg·L-1MMS试剂处理后的白花泡桐丛枝病苗，丛枝症状恢复为健康状态。
附图2所示为开放性增强区域检测到的转录因子家族；
附图3所示为开放性减弱区域检测到的转录因子家族。
具体实施方式
下面对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1
一种基于染色质开放性测序技术研究植原体病害的方法，包括下述步骤：
1.样本处理
取白花泡桐丛枝病苗，分别设置使用60mg·L-1甲基磺酸甲酯(MMS)处理组(PFIM60)和未使用60mg·L-1甲基磺酸甲酯(MMS)处理组(PFI)，培养30d；如图1所示，未使用60mg·L-1甲基磺酸甲酯处理组白花泡桐丛枝病苗具有明显的丛枝症状，使用60mg·L-1MMS试剂处理后的白花泡桐丛枝病苗丛枝症状恢复为健康状态。
取每组长势均匀的泡桐顶芽作为样本，每组设置3个平行试验，液氮速冻保存于-80℃冰箱，为后续试验做准备。
取0.1g干净的样本，用液氮研磨，转移材料到1.5mL预冷的NPB中(Sigma-Aldrich)，裂解15min；用双层Miracloth滤布过滤，并收集液体于新的离心管中，之后，6000g，4℃，离心10min；小心移去上清，并用1mL冰冷的NEB2(aladdin)重悬沉淀，转移悬液到1.5mL离心管中，12000g，4℃，离心10min；移去上清，并用300μL的NEB3重悬沉淀；小心地将悬液加入装有300μL预冷NEB3(aladdin)的1.5mL离心管中，16000g，4℃，离心10min；小心除去上清，并用1mL预冷的NPB重悬沉淀，使用终浓度0.04％的台盼蓝染色，在显微镜下进行细胞核计数；取50000细胞核对应体积的NPB悬液于新的EP管中；500g，4℃，离心10min，弃上清，获得细胞核。
2.文库制备及测序
利用获得的细胞核，使用Tn5转座酶酶切并纯化，最后回收DNA片段，对回收片段进行末端修复、3'端加A、连接测序接头，片段大小选择以及PCA扩增等步骤获得标准的高通量测序文库，文库质检合格后采用IlluminaHiSeq4000进行测序，测序读长为2×150bp。
3.生物信息学分析
Illumina HiSeq 4000测序结果显示染色质开放性测序文库有约为8千万的reads数(表1)。使用deeptools软件分析基因转录起始位点(TSS)上下游2kb的reads分布情况，发现2个样本均在TSS附近均有很明显的富集。
表1测序统计结果
对测序数据进行Peak分析，Peak在基因组范围内的分布区域显示：PFI中有26.52％，PFIM60有27.9％，分布在启动子区。
Peak关联基因分析：从Peak上下游5kb中，找到最近的转录起始位点对应的基因，结果显示，PFI和PFIM60分别鉴定到2,192和2,021个基因。
通过PFI vs.PFIM60 Peak关联基因对比，筛选差异Peak，筛选标准条件为：FDR值小于0.05，并且Fold change大于0。
对差异Peak关联基因区域进行GO富集分析、KEGG富集分析，结果显示，染色质开放性增强的区域主要涉及环境信号传导、遗传信息处理以及细胞加工等通路；开放性减弱的区域主要涉及信号传导、碳水化合物代谢和能量代谢等通路。
使用iTAK软件(version 1.7a)对peak关联的基因进行转录因子(Transcriptionfactors)的鉴定与分类。开放性增强的区域鉴定到28个转录因子家族(图2)，开放性减弱的区域鉴定到88个转录因子家族(图3)。
可根据上述结果进行泡桐丛枝病相关转录因子的筛选，使得泡桐丛枝病相关转录因子的筛选更为快速、精准。
对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对上述实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。
Claims (8)
1.一种基于染色质开放性测序技术研究植原体病害的方法，其特征在于，包括下述步骤：
(1)样本处理
分别取健康植物样本与植原体侵染植物样本，裂解植物细胞获得细胞核；
(2)文库制备及测序
利用步骤(1)获得的细胞核，使用Tn5转座酶酶切并纯化，回收DNA片段，构建高通量测序文库后进行测序；
(3)生物信息学分析
对健康植物样本与植原体侵染植物样本的测序结果进行Peak分析，筛选植原体病害相关基因。
2.根据权利要求1所述的一种基于染色质开放性测序技术研究植原体病害的方法，其特征在于，
步骤(1)中，
白花泡桐丛枝病苗经60mg·L-1甲基磺酸甲酯处理后恢复健康状态后取顶芽作为健康植物样本，并以长势一致的白花泡桐丛枝病苗顶芽作为植原体侵染植物样本。
3.根据权利要求1所述的一种基于染色质开放性测序技术研究植原体病害的方法，其特征在于，
步骤(2)中，
构建好的高通量测序文库用Illumina HiSeq 4000平台进行PE150测序。
4.根据权利要求1所述的一种基于染色质开放性测序技术研究植原体病害的方法，其特征在于，
步骤(3)中，数据分析包括Peak的基本分析、Peak关联基因分析及转录因子分析、GO富集分析、KEGG富集分析。
5.根据权利要求4所述的一种基于染色质开放性测序技术研究植原体病害的方法，其特征在于，
步骤(3)中，
Peak关联基因分析：从Peak上下游5kb中，找到最近的转录起始位点对应的基因。
6.根据权利要求5所述的一种基于染色质开放性测序技术研究植原体病害的方法，其特征在于，
步骤(3)中，
根据健康植物样本与植原体侵染植物样本Peak关联基因筛选差异Peak；
筛选标准条件为：FDR值小于0.05，并且Fold change大于0；
从差异Peak区域鉴定出植原体病害相关转录因子家族。
7.根据权利要求6所述的一种基于染色质开放性测序技术研究植原体病害的方法，其特征在于，
步骤(3)中，
对差异Peak关联基因区域进行GO富集分析、KEGG富集分析，筛选染色质开放性增强的区域和染色质开放性减弱的区域。
8.根据权利要求6所述的一种基于染色质开放性测序技术研究植原体病害的方法，其特征在于，
步骤(3)中，
对染色质开放性增强的区域和染色质开放性减弱的区域内peak关联基因进行转录因子的鉴定与分类，从中筛选植原体病害相关基因。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210006757 | 2022-01-05 | ||
| CN2022100067576 | 2022-01-05 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114350772A true CN114350772A (zh) | 2022-04-15 |
Family
ID=81090876
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210043846.8A Pending CN114350772A (zh) | 2022-01-05 | 2022-01-14 | 一种基于染色质开放性测序技术研究植原体病害的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114350772A (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111763721A (zh) * | 2020-07-14 | 2020-10-13 | 武汉爱基百客生物科技有限公司 | 一种适用于植物的ATAC-seq方法 |
| CN113667769A (zh) * | 2021-08-13 | 2021-11-19 | 河南农业大学 | 一种泡桐丛枝病相关基因网络分析方法 |
-
2022
- 2022-01-14 CN CN202210043846.8A patent/CN114350772A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111763721A (zh) * | 2020-07-14 | 2020-10-13 | 武汉爱基百客生物科技有限公司 | 一种适用于植物的ATAC-seq方法 |
| CN113667769A (zh) * | 2021-08-13 | 2021-11-19 | 河南农业大学 | 一种泡桐丛枝病相关基因网络分析方法 |
Non-Patent Citations (4)
| Title |
|---|
| BAJIC 等: "Identification of Open Chromatin Regions in Plant Genomes Using ATAC-Seq" * |
| WANG 等: "Protocol for assaying chromatin accessibility using ATAC-seq in plants." * |
| YAN 等: "Genome-wide analysis of three histone marks and gene expression in Paulownia fortunei with phytoplasma infection" * |
| 谢博洋: "甲基磺酸甲酯对‘豫杂一号’泡桐丛枝病发生基因表达变化的影响" * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Moghe et al. | Evolutionary routes to biochemical innovation revealed by integrative analysis of a plant-defense related specialized metabolic pathway | |
| Zeng et al. | Identification and analysis of ribosome-associated lncRNAs using ribosome profiling data | |
| Kiss et al. | Biogenesis and intranuclear trafficking of human box C/D and H/ACA RNPs | |
| CN108410907A (zh) | 一种基于CRISPR/Cas9技术实现HMGCR基因敲除的方法 | |
| CN105297142B (zh) | 同时对单细胞基因组和转录组构库及测序的方法基于单细胞整合基因组学的测序方法及应用 | |
| WO2017083766A1 (en) | High-throughput crispr-based library screening | |
| Qian et al. | The Dynamics of DNA methylation in the maize (Zea mays L.) inbred line B73 response to heat stress at the seedling stage | |
| Kelliher et al. | Maize germinal cell initials accommodate hypoxia and precociously express meiotic genes | |
| Pei et al. | Plant 3D genomics: the exploration and application of chromatin organization | |
| Hu et al. | Multiplex CRISPR-Cas9 editing of DNA methyltransferases in rice uncovers a class of non-CG methylation specific for GC-rich regions | |
| Dluhošová et al. | Red clover (Trifolium pratense) and zigzag clover (T. medium)–a picture of genomic similarities and differences | |
| CN111763721A (zh) | 一种适用于植物的ATAC-seq方法 | |
| Qiao et al. | Phenotypic and comparative transcriptome analysis of different ploidy plants in Dendrocalamus latiflorus Munro | |
| Shieh et al. | Computational assessment of miRNA binding to low and high expression HLA-DPB1 allelic sequences | |
| Gueno et al. | Chromatin landscape associated with sexual differentiation in a UV sex determination system | |
| Stark et al. | Extraction of high‐quality, high‐molecular‐weight DNA depends heavily on cell homogenization methods in green microalgae | |
| CN113667769A (zh) | 一种泡桐丛枝病相关基因网络分析方法 | |
| CN114107444A (zh) | 一种发掘调控植物发育关键调控因子的方法及其应用 | |
| CN114350772A (zh) | 一种基于染色质开放性测序技术研究植原体病害的方法 | |
| Orantes-Bonilla et al. | Transgressive and parental dominant gene expression and cytosine methylation during seed development in Brassica napus hybrids | |
| Song et al. | Transcriptome analysis of differential gene expression in the longissimus dorsi muscle from Debao and landrace pigs based on RNA-sequencing | |
| Tian et al. | Temporal-spatial transcriptome analyses provide insights into the development of petaloid androecium in Canna indica | |
| Wang et al. | Non-coding RNA expression analysis revealed the molecular mechanism of flag leaf heterosis in inter-subspecific hybrid rice | |
| Chen et al. | Genome-wide high-throughput chromosome conformation capture analysis reveals hierarchical chromatin interactions during early somatic embryogenesis | |
| Ariizumi et al. | Systematic development of tomato bioresources in Japan |
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 |