本文へ移動
cccskills
無料GitHub で公開

bioconductor-variantannotation

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

インストール方法を見る

含まれるファイル(1)

  • SKILL.md5.8 KB

SKILL.md(原文)

インストールする前に、エージェントに与えられる指示の中身を確認できます。

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 →

レビュー

まだレビューはありません。使ってみた感想をお寄せください。

同じリポジトリのスキル

概要と使いどころ

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-kb8042026年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-kb8042026年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-kb8042026年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-kb8042026年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-kb8042026年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-kb8042026年6月21日 更新

bioMate-AI のスキルをすべて見る

このスキルの問題を報告する