跳到正文
cccskills
免费在 GitHub 公开

bioconductor-extrachips

This package builds on existing tools and adds some simple but extremely useful capabilities for working wth ChIP-Seq data. The focus is on detecting differential binding windows/regions. One set of functions focusses on set-operations retaining mcols for GRanges objects, whilst another group of functions are to aid visualisation of results. Coercion to tibble objects is also implemented.

中文简介正在准备中。当前显示原文说明。

实测评分 未实测

查看安装方法

包含的文件(1)

  • SKILL.md4.7 KB

SKILL.md(原文)

安装前可以查看交给 AI 代理的指令内容。

extraChIPs

Workflows

Standard Workflow

This package builds on existing tools and adds some simple but extremely useful capabilities for working wth ChIP-Seq data. The focus is on detecting differential binding windows/regions. One set of functions focusses on set-operations retaining mcols for GRanges objects, whilst another group of functions are to aid visualisation of results. Coercion to tibble objects is also implemented.

library(extraChIPs)
library(GenomicRanges)

# Define chromosome information
sq <- defineSeqinfo("GRCh37")

# Create mock GRangesList for consensus peak calling
gr1 <- GRanges("chr10:1000-2000", seqinfo = sq)
gr2 <- GRanges("chr10:1200-2200", seqinfo = sq)
peaks <- GRangesList(sample1 = gr1, sample2 = gr2)

# Generate consensus peaks
consensus <- makeConsensus(peaks, p = 0.5)

# Plot overlaps
plotOverlaps(peaks)

Input: A GRangesList of peak calls; Output: A consensus GRanges object and overlap plots.

When to Use

  • Consensus Peak Definition: To define consensus peaks across replicates using makeConsensus.
  • Peak Importing and Filtering: To import and filter peaks with blacklists/greylists using importPeaks.
  • Differential Signal Analysis: To perform differential signal analysis on window/region counts using fitAssayDiff.
  • Gene Mapping: To map peaks/regions to genes and promoters using mapByFeature.
  • Visualization: To visualize peak overlaps with plotOverlaps and signal profiles with plotProfileHeatmap or getProfileData.

When NOT to Use

  • De Novo Motif Discovery: For de novo motif discovery, use memes or universalmotif instead.
  • De Novo Peak Calling: For peak calling itself, use MACS2 or epigraHMM as extraChIPs is designed for downstream analysis of existing peak calls.

Data Requirements

  • Peak Files: Peak files in narrowPeak format or GRanges objects.
  • Read Counts: BAM files for read counting, or a pre-computed RangedSummarizedExperiment object.
  • Annotations: Annotation files (GTF/GFF) for mapping peaks to genes.

Key Parameters

  • p (0.5): Minimum proportion of replicates in which a peak must be present to be included in the consensus.
  • norm ("TMM"): Normalization method in fitAssayDiff (e.g., "TMM" or library-size).
  • fc (1.2): Fold-change threshold incorporated into testing in fitAssayDiff.
  • asRanges (TRUE): Logical indicating whether to return results as a GRanges object.
  • upstream (2500): Upstream distance for defining promoters.
  • downstream (500): Downstream distance for defining promoters.

Best Practices

  • Define Seqinfo: Define a consistent Seqinfo object at the start of the workflow using defineSeqinfo.
  • Exclude Artifacts: Exclude blacklisted and grey-listed regions using importPeaks to avoid false positives.
  • Normalization Check: Test whether group-specific count distributions are similar using quantro before applying TMM normalization.
  • Classify Status: Use addDiffStatus to classify regions as "Increased", "Decreased", or "Unchanged" for downstream visualization.

Common Pitfalls

  • Inappropriate Normalization: Applying TMM normalization when group-specific distributions differ significantly; use quantro to test this assumption first.
  • Metadata Loss: Losing metadata columns during GRanges set operations; use reduceMC or makeConsensus to retain metadata.
  • Mapping Coordinates: Mapping peaks to genes without resetting the core ranges; use colToRanges to restore the original peak boundaries before mapping.

Alternatives

  • DiffBind: For standard affinity-based differential analysis.
  • csaw: For sliding window-based differential binding analysis.
  • ChIPseeker: For peak annotation and visualization.

Citations

  • Ross-Innes et al. 2012, Nature (for DiffBind-style approaches)
  • Hicks and Irizarry 2015, Genome Biol. (for quantro)

References

评价

还没有评价。欢迎分享你的使用感受。

评论

无需登录即可发表。欢迎分享使用感受、提问,或说明在哪些环境下无法运行。

还没有评论。来写第一条吧。

提交后会立即公开。请勿填写邮箱、电话等个人信息。URL 最多 2 个。

0 / 1000 字

此表单受 reCAPTCHA 保护,适用 Google 的隐私权政策和服务条款。

同一仓库的技能

The package contains functions for exploratory oligonucleotide array analysis. The dependence on tkWidgets only concerns few convenience functions. 'affy' is fully functional without it.

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

This tool uses the label-tree function from HyPhy to annotate a phylogenetic tree. It allows users to select a subset of leaves using either a regular expression or a list of sequence names, and then apply a specified label to these selecte

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

This package provides a client for the Bioconductor AnnotationHub web resource. The AnnotationHub web resource provides a central location where genomic files (e.g., VCF, bed, wig) and other resources from standard locations (e.g., UCSC, En

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

apeglm provides Bayesian shrinkage estimators for effect sizes for a variety of GLM models, using approximation of the posterior for individual coefficients.

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

ASURAT is a software for single-cell data analysis. Using ASURAT, one can simultaneously perform unsupervised clustering and biological interpretation in terms of cell type, disease, biological process, and signaling pathway activity. Inputting a single-cell RNA-seq data and knowledge-based databases, such as Cell Ontology, Gene Ontology, KEGG, etc., ASURAT transforms gene expression tables into original multivariate tables, termed sign-by-sample matrices (SSMs).

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

AUCell allows to identify cells with active gene sets (e.g. signatures, gene modules...) in single-cell RNA-seq data. AUCell uses the "Area Under the Curve" (AUC) to calculate whether a critical subset of the input gene set is enriched with

中文简介正在准备中。当前显示原文说明。

bioMate-AI/biomate-bioconductor-kb7932026年6月21日 更新

查看 bioMate-AI 的全部技能

报告此技能的问题