Name Mode Size
.github 040000
R 040000
inst 040000
man 040000
src 040000
tests 040000
vignettes 040000
.gitignore 100644 0 kb
DESCRIPTION 100644 1 kb
NAMESPACE 100644 1 kb
NEWS 100644 2 kb 100644 7 kb
# Fast, epiallele-aware methylation<br> caller and reporter <a href=""><img src="vignettes/epialleleR_logo.svg" alt="logo" align="right" height="50%" /></a> [![](]( [![](]( [![](]( [![install with bioconda](]( ## Introduction *`epialleleR`* is an R package for calling and reporting cytosine methylation and hypermethylated variant epiallele frequencies (VEF) at the level of genomic regions or individual cytosines in next-generation sequencing data using binary alignment map (BAM) files as an input. See below for additional functionality. ### Current Features * calling cytosine methylation and saving calls in BAM file (*`callMethylation`*) * creating sample BAM files given mandatory and optional BAM fields (*`simulateBam`*) * conventional reporting of cytosine methylation (*`generateCytosineReport`*) * calling the hypermethylated variant epiallele frequency (VEF) at the level of genomic regions (*`generate[Bed|Amplicon|Capture]Report`*) or individual cytosines (*`generateCytosineReport`*) * calculating linearised Methylated Haplotype Load (lMHL, *`generateMhlReport`*) * extracting methylation patterns for genomic region of interest (*`extractPatterns`*) * testing for the association between epiallele methylation status and sequence variations (*`generateVcfReport`*) * assessing the distribution of per-read beta values for genomic regions of interest (*`generateBedEcdf`*) ### Recent improvements ##### v1.10 [BioC 3.18] * inputs both single-end and paired-end sequencing alignments * makes and stores methylation calls * creates sample BAM files * reports linearised MHL ##### v1.4 [BioC 3.15] * significant speed-up * method to visualize methylation patterns ##### v1.2 [BioC 3.14] * even faster and more memory-efficient BAM loading (by means of HTSlib) * min.baseq parameter to reduce the effect of low quality bases on methylation or SNV calling (in v1.0 the output of *`generateVcfReport`* was equivalent to the one of `samtools mpileup -Q 0 ...`) check out NEWS for more! ------- ## Installation ### install via Bioconductor ```r if (!requireNamespace("BiocManager", quietly = TRUE)) install.packages("BiocManager") BiocManager::install("epialleleR") ``` ### Install the latest version via install_github ```r library(devtools) install_github("BBCG/epialleleR", build_vignettes=FALSE, repos=BiocManager::repositories(), dependencies=TRUE, type="source") ``` ------- ## Using the package Please read *`epialleleR`* vignette at [GitHub pages]( or within the R environment: `vignette("epialleleR", package="epialleleR")`, or consult the function's help pages for the extensive information on usage, parameters and output values. Comparison of beta, VEF and lMHL values for various use cases is given by the [values]( vignette (`vignette("values", package="epialleleR")`) Very brief synopsis: ```r library(epialleleR) # make methylation calls if necessary callMethylation( input.bam.file=system.file("extdata", "test", "dragen-se-unsort-xg.bam", package="epialleleR"), output.bam.file=tempfile(pattern="output-", fileext=".bam"), genome=system.file("extdata", "test", "reference.fasta.gz", package="epialleleR") ) # make a sample BAM file from scratch simulateBam(output.bam.file=tempfile(pattern="simulated-", fileext=".bam"), pos=c(1, 2), XM=c("ZZZzzZZZ", "ZZzzzzZZ"), XG=c("CT", "AG")) # or use external files amplicon.bam <- system.file("extdata", "amplicon010meth.bam", package="epialleleR") amplicon.bed <- system.file("extdata", "amplicon.bed", package="epialleleR") amplicon.vcf <- system.file("extdata", "amplicon.vcf.gz", package="epialleleR") # preload the data <- preprocessBam(amplicon.bam) # methylation patterns, check vignettes or method description for plotting them patterns <- extractPatterns(bam=amplicon.bam, bed=amplicon.bed, bed.row=3) # CpG VEF report for individual bases <- generateCytosineReport( # BED-guided VEF report for genomic ranges <- generateBedReport(bam=amplicon.bam, bed=amplicon.bed, bed.type="capture") # VCF report <- generateVcfReport(bam=amplicon.bam, bed=amplicon.bed, vcf=amplicon.vcf,"NCBI") # lMHL report <- generateMhlReport(bam=amplicon.bam) ``` ------- ### Citing the *`epialleleR`* package Oleksii Nikolaienko, Per Eystein Lønning, Stian Knappskog, *epialleleR*: an R/Bioconductor package for sensitive allele-specific methylation analysis in NGS data. *GigaScience*, Volume 12, 2023, giad087, []( Data: [GSE201690]( ### Our experimental studies that use the package Per Eystein Lonning, Oleksii Nikolaienko, Kathy Pan, Allison W. Kurian, Hans Petter Petter Eikesdal, Mary Pettinger, Garnet L Anderson, Ross L Prentice, Rowan T. Chlebowski, and Stian Knappskog. Constitutional *BRCA1* methylation and risk of incident triple-negative breast cancer and high-grade serous ovarian cancer. *JAMA Oncology* 2022. []( Oleksii Nikolaienko, Hans P. Eikesdal, Bjørnar Gilje, Steinar Lundgren, Egil S. Blix, Helge Espelid, Jürgen Geisler, Stephanie Geisler, Emiel A.M. Janssen, Synnøve Yndestad, Laura Minsaas, Beryl Leirvaag, Reidun Lillestøl, Stian Knappskog, Per E. Lønning. Prenatal *BRCA1* epimutations contribute significantly to triple-negative breast cancer development. *medRxiv* 2023.05.14.23289949. []( Data: [GSE243966]( ### *`epialleleR`* at Bioconductor [release](, [development version]( ------- License --------- Artistic License 2.0