Abstract Summary Locuszoomr is an R package for visualising and creating publication-ready regional gene locus plots similar to those produced by the original web interface ‘LocusZoom’, but running locally in R. Genetic or genomic data with gene annotation tracks are plotted via R base graphics or ‘ggplot2’, allowing flexibility and easy customisation. Modular plotting functions enable layering of multiple GWAS plots such as for PheWAS, comparing GWAS with eQTL signals, and aligning multiple locus plots on the same page. Interactive plots can be visualised using ‘plotly’ enabling quick identification of data points for labelling. The ‘LDlink’ API (https://ldlink.nih.gov/) can be programmatically queried to obtain linkage disequilibrium (LD) data from the 1000 Genomes Project which is overlaid on plots. The package allows users to query Ensembl databases from Bioconductor or AnnotationHub for genomic information, so it can be used to show gene track information for any species available in Ensembl. Helper functions enable recombination rate data from UCSC to be overlaid on plots. The package has a detailed vignette and includes embedded example GWAS data to allow users to try out different features. Locuszoomr is also fast: identifying 50 peaks in a GWAS of around 8 million SNPs takes around 0.6 seconds. Exporting 50 locus plots with 1Mbp window with recombination rate takes around 30 seconds on Intel 9th gen CPU or 14 seconds on Apple ARM M3. Availability and implementation The R package locuszoomr is open-source and available from CRAN: https://cran.r-project.org/package=locuszoomr and GitHub: https://github.com/myles-lewis/locuszoomr
Read full abstract