跳到正文
cccskills
免费在 GitHub 公开

bioconductor-variantannotation

Annotate variants, compute amino acid coding changes, predict coding outcomes.

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

查看安装方法

包含的文件(1)

  • SKILL.md5.8 KB

SKILL.md(原文)

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

VariantAnnotation

Dependencies & Environment

Package-intrinsic requirements from the Bioconductor landing page — reproduce in any R environment.

  • Version: 1.58.0 · Bioconductor: 3.23 · R: ≥ 4.6
  • Depends: BiocGenerics, MatrixGenerics, Seqinfo, GenomicRanges, SummarizedExperiment, Rsamtools
  • Imports: DBI, Biobase, S4Vectors, IRanges, XVector, Biostrings, AnnotationDbi, rtracklayer, BSgenome, GenomicFeatures, curl
  • System requirements: GNU make
  • Install: BiocManager::install("VariantAnnotation")

When to Use

  • VCF File Parsing: Reading, exploring, and subsetting Variant Call Format (VCF) files using readVcf() and ScanVcfParam().
  • Genomic Context Annotation: Locating variants in and around genes (e.g., coding, intron, 5' UTR, splice site) relative to a transcript database using locateVariants().
  • Coding Consequence Prediction: Computing amino acid coding changes (e.g., synonymous, nonsynonymous, frameshift) for variants in coding regions using predictCoding().
  • Genotype Matrix Conversion: Converting VCF genotype calls or posterior probabilities into a SnpMatrix object for downstream GWAS analysis using genotypeToSnpMatrix().

When NOT to Use

  • For primary variant calling from raw BAM/FASTQ files. Use tools like GATK or bcftools instead.
  • For high-throughput annotation of millions of variants where execution speed is the absolute priority. Command-line tools like Ensembl VEP or SnpEff are faster than R-based loops.

Data Requirements

  • Input Format: Variant Call Format (VCF) text files (can be bgzipped and tabix-indexed).
  • Annotation: A transcript database (TxDb object) matching the VCF's genome assembly.
  • Sequence: A reference genome sequence (BSgenome or fasta file) matching the VCF's assembly for computing codon changes.

Key Parameters

  • param (default: none): A GRanges or ScanVcfParam object passed to readVcf() to restrict import to specific genomic coordinates or VCF fields.
  • info (default: all): Character vector in ScanVcfParam() specifying which INFO fields (e.g., "LDAF", "RSQ") to extract.
  • geno (default: all): Character vector in ScanVcfParam() specifying which FORMAT/genotype fields (e.g., "GT", "DS", "GL") to extract.
  • region (default: none): Constructor passed to locateVariants() specifying the region of interest (e.g., CodingVariants(), AllVariants()).
  • seqSource (default: none): A BSgenome object passed to predictCoding() to retrieve reference sequences.
  • uncertain (default: FALSE): Logical in genotypeToSnpMatrix() indicating whether to use posterior probabilities (GL/GP) instead of called genotypes.

Best Practices

  • Memory Management: Use ScanVcfParam() to read only the necessary info and geno fields, or use TabixFile() with a GRanges object to stream specific genomic coordinates.
  • Chromosome Harmonization: Ensure chromosome names (seqlevels) match exactly between the VCF and the TxDb (e.g., using seqlevels(vcf) <- "chr22") before running locateVariants().
  • Index VCFs: Always index your VCF files using bgzip and tabix (or indexTabix()) to enable fast, random-access querying.
  • Database Integration: Use packages like PolyPhen.Hsapiens.dbSNP131 or SIFT.Hsapiens.dbSNP132 to query pre-computed damage predictions for non-synonymous variants identified by predictCoding().

Common Pitfalls

  • Out of Memory Errors: Attempting to read multi-gigabyte VCF files entirely into memory. Fix: Wrap the file path in a TabixFile and use ScanVcfParam to process the file in chunks or subsets.
  • Empty Overlaps in locateVariants(): The function returns empty results due to mismatched chromosome naming conventions (e.g., "22" vs "chr22"). Fix: Harmonize chromosome names using seqlevels().
  • Incorrect Amino Acid Predictions: predictCoding() returns frameshifts for all variants. Fix: Ensure the BSgenome reference matches the exact assembly version (e.g., hg19) used to call the variants in the VCF.

Alternatives

  • ensemblVEP: An R interface specifically for running the Ensembl Variant Effect Predictor.
  • snpStats: For downstream statistical analysis of SNP data once converted via genotypeToSnpMatrix().

Citations

  • Obenchain V, et al. (2014). VariantAnnotation: a Bioconductor package for functional annotation of genomic variants. Bioinformatics, 30(14), 2075-2076.

References

<!-- biomate-cta -->

Run this on BioMate

This skill is the knowledge layer — when, why, and how to use variantannotation. To run this analysis on your own data with managed compute, automated QC, and reproducible outputs, use BioMate — free to start.

▶ Open variantannotation on BioMate →

评价

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

评论

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

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

提交后会立即公开。请勿填写邮箱、电话等个人信息。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-kb7972026年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-kb7972026年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-kb7972026年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-kb7972026年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-kb7972026年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-kb7972026年6月21日 更新

查看 bioMate-AI 的全部技能

报告此技能的问题