1 项目概述

项目编号 GDE21070121-1_std_1
项目名称 Benincasa hispida ------ 2 DAP-seq
参考基因组 ncbi_GCF_009727055.1
样品名称 BhYABBY4-IP   ¦   B214-input

2 项目介绍

DNA亲和纯化测序（DAP-seq，DNA Affinity Purification sequencing），通过体外构建表达转录因子（TF，Transcription Factor）蛋白，与目标基因组片段结合，由此可以把目的蛋白所结合的基因组DNA片段富集下来。通过与高通量测序技术的结合，对DAP后的DNA产物进行测序分析，从全基因组范围内寻找目的蛋白的DNA结合位点，以高效率的测序手段得到高通量的数据结果。

2.1 建库流程

得到 DNA 样品后，首先对样品进行质量检测。样品质量检测合格后，进行文库构建，具体步骤如下：

1. DNA超声打断至约200bp大小，末端修复，3’ 端加A碱基，连接测序接头。

2. 在DAP实验完成后，进行DNA片段的PCR扩增，产物纯化。

3. 文库构建完成后，对文库质量进行检测，合格的文库用于上机测序。


Fig 2-1-1 DAP-seq实验建库流程图

2.2 信息分析流程

对下机数据进行数据过滤，去掉低质量数据。将过滤后数据与参考基因组进行比对，确认比对质量合格后，提取比对上唯一位置的序列，进行信息分析处理，得到标准信息分析结果和个性化分析结果。具体信息分析流程见图所示。

Fig 2-2-1 DAP-seq分析流程图

广州基迪奥生物科技有限公司



3 数据分析与质控

3.1 过滤信息统计
Tab 3-1-1 Reads 过滤信息统计表
SampleClean_Reads_NumHQ_Clean_Reads_Num(%)Read_lengthadapter(%)low_quality(%)polyA(%)N(%)
BhYABBY4-IP7922165079046930(99.78%)150 + 15037776(0.05%)273888(0.17%)0(0.0%)0(0.0%)
B214-input9405458093842944(99.77%)150 + 15045790(0.05%)331692(0.18%)0(0.0%)0(0.0%)

Fig 3-1-1 样本过滤分析频数图 Fig 3-1-2 样本过滤分析比例图

3.2 碱基组成与质量分析
Tab 3-2-1 过滤前后碱基信息统计表
SampleClean_Data(bp)Q20(%)Q30(%)N(%)GC(%)HQ_Clean_Data(bp)HQ_Q20(%)HQ_Q30(%)HQ_N(%)HQ_GC(%)
BhYABBY4-IP1188324750011584420869(97.49%)11052177423(93.01%)56616(0.0%)3860484867(32.48%)1126343156211008808639(97.74%)10523404343(93.43%)53670(0.0%)3539014128(31.42%)
B214-input1410818700013783848554(97.7%)13180853453(93.43%)66923(0.0%)5113842838(36.25%)1353494891713248166188(97.88%)12686790548(93.73%)64164(0.0%)4817304043(35.59%)

  • BhYABBY4-IP
  • B214-input

Fig 3-2-1 过滤前后碱基分布


广州基迪奥生物科技有限公司



4 比对分析

4.1 比对基因组统计

利用比对软件 Bowtie2[1] 将数据比对到参考基因组，过滤掉低比对质量值的 reads ，最终将比对到基因组上唯一位置的 reads (唯一比对 reads )用于后续的信息分析。

Tab 4-1-1 各样品比对参考基因组统计结果
Sample_IDTotal_ReadsUnMapped_ReadsMapped_ReadsMulti_Mapped_ReadsUnique_Mapped_ReadsDuplicate_Reads
BhYABBY4-IP790469302509031-3.17%76537899-96.83%34464079-43.60%42073820-53.23%24156432-31.56%
B214-input938429443032645-3.23%90810299-96.77%47836147-50.97%42974152-45.79%19225638-21.17%

4.2 基因组测序深度累积分布

以比对后得到的唯一比对序列为分析对象，分析其在参考基因组上的覆盖分布，统计基因组位点的深度信息，得到基因组上测序深度统计结果。

  • BhYABBY4-IP
  • B214-input

Fig 4-2-1 测序深度累积分布统计图


4.3 基因组测序深度分布

利用 deepTools[2] 软件对转录起始位点（TSS）到转录终止位点（TES）区间，及基因上下游 2k 区间的所有 reads 进行统计，设置大小为 50bp 的窗口，计算每个窗口内的平均 reads 深度，并绘制出 reads 相对于基因特定位置的分布折线图。

  • BhYABBY4-IP
  • B214-input

Fig 4-3-1 测序深度分布图


4.4 测序饱和度分析

测序饱和度分析是用来衡量一个样品的测序量是否达到饱和的标准。随着测序量（reads数量）的增多，检测到的基因数也随之上升。当测序量达到某个值时，其检测到的基因数增长速度趋于平缓，说明检测到的基因数趋于饱和，如下图所示：

  • BhYABBY4-IP
  • B214-input

Fig 4-4-1 测序饱和度分布


4.5 Reads在染色体上的分布

将唯一比对、去重复后的 reads 比对到基因组上各个染色体（分正负链）的密度进行统计。

  • BhYABBY4-IP
  • B214-input

Fig 4-5-1 reads在染色体上的分布图


广州基迪奥生物科技有限公司



5 Peak分析

测序分析 DAP 后的 DNA 产物，从全基因组范围内寻找目的蛋白与 DNA 的结合位点，利用 MACS2[3] 分析软件在全基因组范围进行 peak（DAP Sequencing富集区域）扫描（peak calling），阈值为 q-value<0.05 ，并对 peak 在基因组上的位置信息， peak 区域序列信息等进行分析，筛选出 peak 相关基因。我们将对各处理组分别检测peak。如果一个处理组有多个生物学重复，MACS软件将合并多个生物学重复的数据，输出统一的peak数据。

5.1 PeakCalling
Tab 5-1-1 各样本 peak 数目统计表
SampleIdPeakNumberTotalLengthAverageLengthTotalPileupAveragePileupGenomeRatio(%)FRiP(%)
BhYABBY4-IP4098812068686294814496191.34%10.32%

Tab 5-1-2 样本 BhYABBY4-IP peaks结果信息统计表（前10行）
chrstartendlengthabs_summitpileup-log10(pvalue)fold_enrichment-log10(qvalue)name
NC_052349.119556537134714.003.947452.640131.93038BhYABBY4-IP_peak_1
NC_052349.113111524214143419.006.321493.274863.70891BhYABBY4-IP_peak_2
NC_052349.177157910196780717.005.371573.045562.99084BhYABBY4-IP_peak_3
NC_052349.112664712681817212676718.006.149263.293263.59354BhYABBY4-IP_peak_4
NC_052349.112710912729718912724521.008.220873.865175.22626BhYABBY4-IP_peak_5
NC_052349.114935214959824714938617.004.282472.586832.18159BhYABBY4-IP_peak_6
NC_052349.117460617478217717475919.006.903533.520184.16497BhYABBY4-IP_peak_7
NC_052349.119489419550160819507228.0014.035795.3326010.10866BhYABBY4-IP_peak_8
NC_052349.119631919661329519639228.007.529623.048824.67374BhYABBY4-IP_peak_9
NC_052349.120372220394222120383619.006.919183.526884.17848BhYABBY4-IP_peak_10

5.2 Peak长度分布
Fig 5-2-1 BhYABBY4-IP

5.3 Peak深度分布
Fig 5-3-1 BhYABBY4-IP

5.4 Peak富集倍数分布
Fig 5-4-1 BhYABBY4-IP

5.5 Peak置信程度分布
Fig 5-5-1 BhYABBY4-IP

广州基迪奥生物科技有限公司



6 Peak注释

利用 ChIPseeker[4] R包，对 peak 相关基因进行注释。

6.1 Peak相关基因分析

根据peak在基因组上的区域信息及基因的注释信息，得到关于peak相关基因，从一定程度上表示了目的蛋白或特定组蛋白修饰可能调控的靶基因区域。

Tab 6-1-1 样本BhYABBY4-IP peak相关基因注释表（前10行）
peak_namechrstartendlengthabs_summitpileupfold_enrichment-log10(pvalue)-log10(qvalue)annotationgeneChrgeneStartgeneEndgeneLengthexonNumbergeneStrandgeneIdtranscriptIddistToTSSSymbolDescriptionKEGG_A_classKEGG_B_classPathwayK_IDGO_ComponentGO_FunctionGO_Process
BhYABBY4-IP_peak_9NC_052349.119631919661329519639228.003.048827.529624.67374Downstream (1-2kb)NC_052349.119830320362253203-ncbi_120074624XM_039027805.17230At2g33170XP_008438405.1 PREDICTED: probable leucine-rich repeat receptor-like protein kinase At5g63930 [Cucumis melo]-----GO:0004672//protein kinase activity;GO:0005515//protein binding;GO:0005524//ATP bindingGO:0006468//protein phosphorylation
BhYABBY4-IP_peak_10NC_052349.120372220394222120383619.003.526886.919184.17848Promoter (<=1kb)NC_052349.119830320362253203-ncbi_120074624XM_039027805.1-214At2g33170XP_008438405.1 PREDICTED: probable leucine-rich repeat receptor-like protein kinase At5g63930 [Cucumis melo]-----GO:0004672//protein kinase activity;GO:0005515//protein binding;GO:0005524//ATP bindingGO:0006468//protein phosphorylation
BhYABBY4-IP_peak_14NC_052349.121886221916430321895317.003.348936.057163.51753Promoter (1-2kb)NC_052349.1220913225139422710+ncbi_120087383XM_039044359.1-1960PAM71-HLXP_008438406.1 PREDICTED: GDT1-like protein 2, chloroplastic [Cucumis melo]----GO:0009512//cytochrome b6f complexGO:0009055//electron carrier activity-
BhYABBY4-IP_peak_15NC_052349.126639026662924026649719.002.708594.902952.64288Promoter (<=1kb)NC_052349.126201626626142465-ncbi_120069147XM_039020829.1-236GLCAT14AXP_022974530.1 beta-glucuronosyltransferase GlcAT14A [Cucurbita maxima]----GO:0016020//membraneGO:0016757//transferase activity, transferring glycosyl groups-
BhYABBY4-IP_peak_16NC_052349.129237529277640229250823.003.966199.008035.84187Promoter (<=1kb)NC_052349.1284761291768700816-ncbi_120073368XM_039026194.1-740SERINC3XP_008438413.1 PREDICTED: probable serine incorporator isoform X2 [Cucumis melo]----GO:0016020//membrane;GO:0016021//integral component of membrane-GO:0007155//cell adhesion;GO:0042246//tissue regeneration
BhYABBY4-IP_peak_18NC_052349.136464036484420536476917.002.437733.914421.91067Downstream (1-2kb)NC_052349.136637536763612622-ncbi_120071386XM_039023646.12867PRP1XP_008438425.1 PREDICTED: probable glutathione S-transferase [Cucumis melo]Metabolism;MetabolismGlobal and overview maps;Metabolism of other amino acidsko01100//Metabolic pathways;ko00480//Glutathione metabolismK00799;K00799-GO:0005515//protein binding;GO:0015035//protein disulfide oxidoreductase activity;GO:0016034//maleylacetoacetate isomerase activityGO:0006559//L-phenylalanine catabolic process;GO:0006572//tyrosine catabolic process;GO:0006749//glutathione metabolic process;GO:0006790//sulfur compound metabolic process;GO:0042965//glutaredoxin biosynthetic process;GO:0045454//cell redox homeostasis;GO:0055114//oxidation-reduction process
BhYABBY4-IP_peak_19NC_052349.136816436853737436830818.003.581936.806364.10294Promoter (<=1kb)NC_052349.136637536763612622-ncbi_120071386XM_039023646.1-672PRP1XP_008438425.1 PREDICTED: probable glutathione S-transferase [Cucumis melo]Metabolism;MetabolismGlobal and overview maps;Metabolism of other amino acidsko01100//Metabolic pathways;ko00480//Glutathione metabolismK00799;K00799-GO:0005515//protein binding;GO:0015035//protein disulfide oxidoreductase activity;GO:0016034//maleylacetoacetate isomerase activityGO:0006559//L-phenylalanine catabolic process;GO:0006572//tyrosine catabolic process;GO:0006749//glutathione metabolic process;GO:0006790//sulfur compound metabolic process;GO:0042965//glutaredoxin biosynthetic process;GO:0045454//cell redox homeostasis;GO:0055114//oxidation-reduction process
BhYABBY4-IP_peak_20NC_052349.136935036954319436945915.002.871454.611352.41714Promoter (1-2kb)NC_052349.136637536763612622-ncbi_120071386XM_039023646.1-1823PRP1XP_008438425.1 PREDICTED: probable glutathione S-transferase [Cucumis melo]Metabolism;MetabolismGlobal and overview maps;Metabolism of other amino acidsko01100//Metabolic pathways;ko00480//Glutathione metabolismK00799;K00799-GO:0005515//protein binding;GO:0015035//protein disulfide oxidoreductase activity;GO:0016034//maleylacetoacetate isomerase activityGO:0006559//L-phenylalanine catabolic process;GO:0006572//tyrosine catabolic process;GO:0006749//glutathione metabolic process;GO:0006790//sulfur compound metabolic process;GO:0042965//glutaredoxin biosynthetic process;GO:0045454//cell redox homeostasis;GO:0055114//oxidation-reduction process
BhYABBY4-IP_peak_26NC_052349.153897453923326053908418.003.111125.721953.25237Intron (XM_039046086.1/ncbi_120088661, intron 10 of 35)NC_052349.15195375443952485936-ncbi_120088661XM_039046086.15311XI-BXP_008438435.1 PREDICTED: myosin-6 isoform X1 [Cucumis melo]----GO:0016459//myosin complexGO:0000155//phosphorelay sensor kinase activity;GO:0003774//motor activity;GO:0004672//protein kinase activity;GO:0005515//protein binding;GO:0005524//ATP bindingGO:0000160//phosphorelay signal transduction system;GO:0006109//regulation of carbohydrate metabolic process;GO:0007059//chromosome segregation;GO:0030261//chromosome condensation
BhYABBY4-IP_peak_27NC_052349.155309855360650955317621.003.836688.152565.16841Downstream (<1kb)NC_052349.155364655669930544-ncbi_120084783XM_039040606.13523BAM1XP_008438436.1 PREDICTED: beta-amylase 1, chloroplastic-like [Cucumis melo]Metabolism;Metabolism;MetabolismGlobal and overview maps;Global and overview maps;Carbohydrate metabolismko01100//Metabolic pathways;ko01110//Biosynthesis of secondary metabolites;ko00500//Starch and sucrose metabolismK01177;K01177;K01177GO:0009341//beta-galactosidase complexGO:0004565//beta-galactosidase activity;GO:0016161//beta-amylase activityGO:0000272//polysaccharide catabolic process;GO:0005975//carbohydrate metabolic process

6.2 Peak在基因功能元件上的分布

利用参考基因组gtf数据库将与peak最相邻的基因（peak中心与基因的TTS距离最近）进行注释，统计peak在基因各功能区域的分布情况，按照promoter（转录起始位点上游2K序列范围）、5'UTR、3'UTR、Exon、Intron、Downstream（转录终止位点下游2K序列范围）、Intergenic（基因间区）的优先级顺序，确定peak与各个功能区域的对应关系。

Tab 6-2-1 各样本peak在基因功能元件上的分布统计表
samplepeakCountpromoter5'UTR3'UTR1st_exonother_exon1st_intronohter_introndownstreamdistal_intergenic
BhYABBY4-IP409885836-14.24%106-0.26%260-0.63%41-0.10%106-0.26%1085-2.65%2534-6.18%2975-7.26%28045-68.42%

Fig 6-2-1 BhYABBY4-IP

Fig 6-2-2 BhYABBY4-IP


Fig 6-2-3 各样本peak在基因功能元件上的分布比例图 Fig 6-2-4 各样本peak相对TSS距离分布比例图

6.3 Peak在染色体上的分布
Fig 6-3-1 BhYABBY4-IP

6.4 peak相关基因 GO富集分析

       Gene Ontology（简称GO）是一个国际标准化的基因功能分类体系，提供了一套动态更新的标准词汇表（controlled vocabulary）来全面描述生物体中基因和基因产物的属性。GO总共有三个ontology（本体），分别描述基因的分子功能（molecular function）、细胞组分（cellular component）、参与的生物过程（biological process）。GO的基本单位是term（词条、节点），每个term都对应一个属性。 GO功能分析一方面给出基因的GO功能分类注释；另一方面给出基因的GO功能显著性富集分析。 首先，我们将基因向GO数据库(http://www.geneontology.org/)的各term映射，并计算每个term的基因数，从而得到具有某个GO功能的基因列表及基因数目统计。然后应用超几何检验，找出与整个基因组背景相比，在基因中显著富集的GO条目。

基因集 细胞组分 分子功能 生物学过程 GO 分类表
BhYABBY4-IP BhYABBY4-IP.C.html BhYABBY4-IP.F.html BhYABBY4-IP.P.html BhYABBY4-IP.Level2.xls

GO富集圈图：(第一圈：富集前20的GOterm,圈外为基因数目的坐标尺。不同的颜色代表不同的Ontology; 第二圈：背景基因中该GOterm的数目以及Q值。基因越多条形越长，Q值越小颜色越红； 第三圈：该GOterm差异基因数量 第四圈：各GOterm的RichFactor值(该GOterm中差异数量除以所有数量),背景网格线，每一格代表0.1)

  • BhYABBY4-IP 富集圈图

Fig 6-4-1 GO 富集圈图


GO 富集分类柱状图：（横坐标为二级GOterm，纵坐标为该term里的基因数量,不同颜色表色不同类型GOterm）

  • BhYABBY4-IP

Fig 6-4-2 GO富集分类柱状图


GO富集气泡图：(利用Q值最小的前20个GOterm来作图，纵坐标为GOterm，横坐标为富集因子(该GOterm中差异数量除以所有数量)，大小表示数量多少，颜色越红Q值越小)");

  • BhYABBY4-IP.C 气泡图
  • BhYABBY4-IP.F 气泡图
  • BhYABBY4-IP.P 气泡图

Fig 6-4-3 GO富集气泡图


GO富集条形图：(利用Q值最小的前20个GOterm来作图，纵坐标为GOterm，横坐标为该GOterm数目占所有差异数目的百分比，颜色越深Q值越小，柱子上的数值为该GOterm数量及Q值");

  • BhYABBY4-IP.C 富集柱形图
  • BhYABBY4-IP.F 富集柱形图
  • BhYABBY4-IP.P 富集柱形图

Fig 6-4-4 GO富集条形图


6.5 peak相关基因 KO富集分析

       在生物体内，不同基因相互协调行使其生物学，基于Pathway的分析有助于更进一步了解基因的生物学功能。KEGG是有关Pathway的主要公共数据库。 Pathway显著性富集分析以KEGG Pathway为单位，应用超几何检验，找出与整个基因组背景相比，在基因中显著性富集的Pathway。通过Pathway显著性富集能确定基因参与的最主要生化代谢途径和信号转导途径。

基因集 Pathway 富集结果 Pathway 注释表
BhYABBY4-IP BhYABBY4-IP.htm BhYABBY4-IP.path.xls

KO富集圈图：(第一圈：富集前20的pathway,圈外为基因数目的坐标尺。不同的颜色代表不同的A class; 第二圈：背景基因中该pathway的数目以及Q值。基因越多条形越长，Q值越小颜色越红； 第三圈：该pathway差异基因数量 第四圈：各pathway的RichFactor值(该pathway中差异数量除以所有数量),背景网格线，每一格代表0.1)

  • BhYABBY4-IP 富集圈图

Fig 6-5-1 KO 富集圈图


KO富集气泡图：(利用Q值最小的前20个pathway来作图，纵坐标为pathway，横坐标为富集因子(该pathway中差异数量除以所有数量)，大小表示数量多少，颜色越红Q值越小)");

  • BhYABBY4-IP

Fig 6-5-2 KO富集气泡图


KO富集条形图：(利用Q值最小的前20个pathway来作图，纵坐标为pathway，横坐标为该pathway数目占所有差异数目的百分比，颜色越深Q值越小，柱子上的数值为该pathway数量及Q值");

  • BhYABBY4-IP 富集柱形图

Fig 6-5-3 KO富集条形图


6.6 Peak以及周边基因结构的可视化

Integrative Genomics Viewer(IGV)是一种高性能的可视化工具，用于基因组中各类数据的整合可视化，可选择全基因组、染色体甚至是单个碱基水平进行浏览，并且能显示其他注释信息。我们提供reads在每条染色体比对结果的bigwig文件，可以使用IGV将结果进行可视化。

Fig 6-6-1 基因结构可视化

广州基迪奥生物科技有限公司



7 Motif分析

转录因子(Transcription Factors, TFs)是指能够以序列特异性方式结合 DNA 并且调节转录的蛋白质，TFs 通过识别特定的 DNA 序列来控制染色质和转录，以形成指导基因组表达的复杂系统。利用 MEME Suite（http://meme-suite.org/）进行 Motif 分析。

7.1 motif检测

我们利用 MEME Suite[5] 中的 MEME-chip 检测 peak 序列中显著 motif 序列。MEME-chip 整合了 MEME（检测8-15bp）与 Dreme（检测3-8bp）功能，可同时检测长 motif 和短 motif 。

  • BhYABBY4-IP motifs









Fig 7-1-1 meme 显著motif序列图


不同颜色代表不同碱基类型，字母的高度代表这个碱基的保守性（字母越高，其在整个位点的频率越高，越保守）。

  • BhYABBY4-IP motifs


Fig 7-1-2 dreme 显著motif序列图


不同颜色代表不同碱基类型，字母的高度代表这个碱基的保守性（字母越高，其在整个位点的频率越高，越保守）。

7.2 各motif到peak顶点（submit）的距离统计

DAP-seq通常会检测到多个显著的motif，但理论上其中只有一个为转录因子的结合位点序列。转录因子结合位点序列理论上应该与peak的顶点重叠，即motif与peak顶点的距离应该趋近于0。通过统计各个motif与peak的顶点的距离，将有利于帮助我们判断哪个motif最有可能为转录因子的结合位点。

我们利用 MEME Suite[5] 中的 fimo 确定 MEME-chip 结果中显著性 Top5 的 motif 在各序列中的位置。

Tab 7-2-1 样本BhYABBY4-IP motif位置统计表（前10行）
motif_idmotif_alt_idsequence_namestartstopstrandscorep-valueq-valuematched_sequencedistToSummitpeak_namechrstartendlengthabs_summitpileupfold_enrichment-log10(pvalue)-log10(qvalue)annotationgeneChrgeneStartgeneEndgeneLengthexonNumbergeneStrandgeneIdtranscriptIddistToTSS
TTTHTYTTYTTYTTTMEME-1NC_052359.11521734715217361-20.57581.04e-090.00373TTTCTCTCCTTCTTT60BhYABBY4-IP_peak_32139NC_052359.115217169152174492811521728722.003.924298.631795.54727Distal IntergenicNC_052359.1152327991523679039922+ncbi_120091818XM_039049959.1-15512
TTTHTYTTYTTYTTTMEME-1NC_052352.13431933634319350-20.16671.3e-090.00373TTTCTCTCCTTCTCT157BhYABBY4-IP_peak_11996NC_052352.134319123343196004783431917918.002.573164.398092.26593Distal IntergenicNC_052352.1343394073434558061748-ncbi_120075660XM_039029254.126401
TTTHTYTTYTTYTTTMEME-1NC_052356.17661476628+19.68185.74e-090.00373tttctcttcttcttt37BhYABBY4-IP_peak_21344NC_052356.176462767492887657717.003.553896.505833.85688Distal IntergenicNC_052356.192196927205251-ncbi_120083073XM_039038607.116143
TTTHTYTTYTTYTTTMEME-1NW_024064913.1150136150150-19.68185.74e-090.00373TTTCTCTTCTTCTTT-40BhYABBY4-IP_peak_40408NW_024064913.114982715029546915019021.004.554549.804626.49649Distal IntergenicNW_024064913.11793851797113271-ncbi_120069666XM_039021448.129521
TTTHTYTTYTTYTTTMEME-1NW_024064891.1864986865000+19.68185.74e-090.00373tttctcttcttcttt-3BhYABBY4-IP_peak_40314NW_024064891.186493886510416786500320.003.995488.255445.23524Distal IntergenicNW_024064891.18381148387556421+ncbi_120069592XM_039021382.126889
TTTHTYTTYTTYTTTMEME-1NC_052352.190769079076921-19.68185.74e-090.00373TTTCTCTTCTTCTTT-39BhYABBY4-IP_peak_10881NC_052352.190768849077341458907696022.003.996348.807785.70189Downstream (1-2kb)NC_052352.19070948907416132147-ncbi_120075543XM_039029029.1-2799
TTTHTYTTYTTYTTTMEME-1NC_052360.195909049590918+19.68185.74e-090.00373tttctcttcttcttt-70BhYABBY4-IP_peak_35832NC_052360.195908809591139260959098824.004.023719.430926.20417Distal IntergenicNC_052360.19651548965442928824+ncbi_120068030XM_039019699.1-60560
TTTHTYTTYTTYTTTMEME-1NC_052355.11553136915531383+19.68185.74e-090.00373tttctcttcttcttt209BhYABBY4-IP_peak_19442NC_052355.115531088155314954081553116019.004.150848.320775.28090Distal IntergenicNC_052355.11543488515440791590710-ncbi_120080991XM_039035674.1-90369
TTTHTYTTYTTYTTTMEME-1NC_052360.12671691426716928-19.68185.74e-090.00373TTTCTCTTCTTCTTT-3BhYABBY4-IP_peak_36639NC_052360.126716599267169423442671693118.002.307203.708691.76631Distal IntergenicNC_052360.1267292912673222129311+ncbi_120092882XM_039051125.1-12360
TTTHTYTTYTTYTTTMEME-1NC_052356.12981736529817379-19.68185.74e-090.00373TTTCTCTTCTTCTTT0BhYABBY4-IP_peak_22653NC_052356.129817225298174452212981736527.005.7966914.7557910.73444Distal IntergenicNC_052356.1297914162979334119262-ncbi_120083489XR_005483436.1-24024

  • BhYABBY4-IP motifs





Fig 7-2-1 motif与peak距离的频率直方图


7.3 motif注释

我们利用 MEME-chip (地址：http://meme-suite.org/) 软件中的 Tomtom 模块将得到的 motif 序列与已知 motif 数据库进行比对，利用已知 motif 来注释检测到的 motif 。

广州基迪奥生物科技有限公司



8 目录结构

upload
├── 1.data_access                                    过滤统计目录
│   ├── all.read.stat.xls                               过滤信息统计表
│   ├── all.data.stat.xls                               过滤前后碱基信息统计表
│   ├── all.count.{png,pdf}                             样本过滤分析频数图
│   ├── all.fill.{png,pdf}                              样本过滤分析比例图
│   ├── *.old.png                                       各样品过滤前各项统计图
│   └── *.new.png                                       各样品过滤后各项统计图
├── 2.alignment                                      比对统计结果
│   ├── sample.alignstat.xls                            各样品比对参考基因组统计结果
│   ├── *.cov.{png,pdf}                                 测序深度累积分布统计图
│   ├── *.tss_tes.{png,pdf}                             测序深度分布图
│   ├── *.saturation.{png,pdf}                          测序饱和度分布图
│   └── *.chr.{png,svg}                                 reads在染色体上的分布图
├── 3.peak_analysis                                  Peak分析结果
│   ├── sample.peakstat.xls                             各样本 peak 数目统计表
│   ├── *.peak.len.{png,pdf}                            Peak长度分布
│   ├── *.peak.dep.{png,pdf}                            Peak深度分布
│   ├── *.peak.fe.{png,pdf}                             Peak富集倍数分布
│   └── *.peak.pv.{png,pdf}                             Peak置信程度分布
├── 4.peak_annotation                                Peak注释
│   ├── all.annostat.xls                                peak在基因功能元件上的分布统计表
│   ├── all.peakanno.bar.{png,pdf}                      peak在基因功能元件上的分布饼图
│   ├── all.peakanno.distToTSS.{png,pdf}                peak相对TSS距离分布比例图
│   ├── *.final.anno.xls                                peak相关基因注释表
│   ├── *.peakanno.bar.{png,pdf}                        peak在基因功能元件上的分布比例图
│   ├── *.peakanno.pie.{png,pdf}                        peak在基因功能元件上的分布饼状图
│   ├── *.peakanno.cov.{png.pdf}                        peak在染色体上的分布图
│   ├── *.distToTss.{png,pdf}                           peak_summit相对基因位置分布图
│   └── enrich                                          富集结果目录
├── 5.motif                                          Motif分析结果
│   ├── meme-chip                                       样本meme-chip结果目录
│   │   └── meme-chip.html                                 样本meme-chip网页报告
│   └── fimo                                            样本fimo结果目录
│       ├── fimo.html                                      样本fimo网页报告
│       ├── *.motif.dist.png                               样本motif距离peakSummit分布图
│       └── *.motif_anno.xls                               样本motif基因注释表
├── src                                              结果报告内容                                    
│   ├── content.html                                    结果报告内容
│   ├── css                                             结题报告css脚本
│   ├── js                                              结题报告js脚本
│   ├── doc                                             结题报告说明文档
│   └── image                                           结题报告图片
└── index.html                                       网页版结题报告

广州基迪奥生物科技有限公司



9 参考文献

广州基迪奥生物科技有限公司



10 附录

10.1 分析方法英文文档

分析方法文档（英文）：DAP-seq_method.pdf

10.2 结果文件查看

*.xls,*.txt ：结果数据表格文件，文件以制表符（Tab）分隔。unix/Linux/Mac用户使用 less 或 more 命令查看；windows用户使用高级文本编辑器Notepad++ 等查看，也可以用Microsoft Excel打开。

*.png：结果图像文件，位图，无损压缩。

*.pdf：结果图像文件，矢量图，可以放大和缩小而不失真，方便用户查看和编辑处理，可使用Adobe Illustrator进行图片编辑，用于文章发表等。

10.3 文章引用与致谢

如果您的研究课题使用了基迪奥的测序和分析服务，我们期望您在论文发表时，在Method部分或Acknowledgements部分引用或提及基迪奥公司。以下语句可供参考：

广州基迪奥生物科技有限公司



帮助文档