VCF data in RVcfR: a package to manipulate and visualize
VcfR is an R package intended to allow easy manipulation and visualization of variant call format (VCF) data. Functions are provided to rapidly read from and write to VCF files. Once VCF data is read into R a parser function extracts matrices from the VCF data for use with typical R functions. This information can then be used for quality control or other purposes. Additional functions provide visualization of genomic data. Once processing is complete data may be written to a VCF file or converted into other popular R objects (e.g., genlight, DNAbin). VcfR provides a link between VCF data and the R environment connecting familiar software with genomic data.
VcfR is built upon two data structures.
vcfR - S4 class to contain data from a VCF file.
chromR - S4 class to contain variant information (VCF) and optional sequence (FASTA) and annotation (GFF) information.
Functions in vcfR provide the ability to subset VCF data as well as to extract and parse the data. For example, individual genotypes, sequence depths or genotype likelihoods (when provided in the VCF file) can easily be accessed. These tools are provided to aid researchers in rapidly surveying the quality and other characteristics of data provided as VCF data. With this information in hand, researchers should be able to determine criteria for hard filtering in order to attempt to maximize biological variation and minimize technical variation.
Documentation for vcfR can now be found here: vcfR_documentation.
We also have Population genetics and genomics in R which is more general and provides examples of analyses.
If you think you've found a bug, please see reporting an issue.
Knaus, Brian J., and Niklaus J. Grunwald. 2017. VCFR: a package to manipulate and visualize variant call format data in R. Molecular Ecology Resources 17(1):44-53. http://dx.doi.org/10.1111/1755-0998.12549.
Knaus, Brian J., and Niklaus J. Grunwald. 2016. VcfR: an R package to manipulate and visualize VCF format data. bioRxiv: 041277. http://dx.doi.org/10.1101/041277.
Copy number variation article
Knaus, Brian, and Niklaus J. Grünwald. 2018. Inferring variation in copy number using high throughput sequencing data in R. Frontiers in Genetics 9: 123. http://dx.doi.org/10.3389/fgene.2018.00123.
vcfR is available at CRAN. To install use:
The development version can be installed through github:
If you would like the vignettes use:
If you've built the vignettes, you can browse them with:
If you've installed this package with devtools you will probably need to run:
devtools::install(build_vignettes = TRUE)
The devel branch (which may not be stable) can also be installed:
And to build the vignettes:
Software that produce VCF files
A fun part of this project has been learning about how people use vcfR. One facet of this is learning about the software that create VCF files. So I've decided to make a list of these software. If you know of a software that I have not included on this list, particularly if you can report that vcfR works with its files, feel free to let me know!
Genomic variant callers:
- GATK haplotype caller
- GATK MuTect2
- GATK GenotypeGVCFs
Restriction site associated DNA markers (e.g., RADseq, GBS):
Manipulation of VCF data:
- Beagle v4.1 Inputs VCF genotypes and outputs phased genotypes to VCF format
- pegas::read.vcf Population and Evolutionary Genetics Analysis System
- SnpEff Genetic variant annotation and effect prediction toolbox
- Picard A set of command line tools (in Java) for manipulating high-throughput sequencing (HTS) data and formats such as SAM/BAM/CRAM and VCF
- VCF-kit VCF-kit is a command-line based collection of utilities for performing analysis on Variant Call Format (VCF) files.
- VCFtools General manipulation and analysis
- VariantAnnotation::readVcf Bioconductor package for annotating variants
R packages that read VCF data