Home

Awesome

:mount_fuji: Fuji plot :mount_fuji:

a circos representation of multiple GWAS results.

Overview

Fuji plot (a circos representation of multiple GWAS results) was developed to efficiently summarize significant SNPs across dozens of traits (Kanai, M. et al., Nat. Genet. 2018). Technically, it is a wrapper script for Circos, generating input data and configuration files from GWAS results.

<p align="center"><img src="https://mkanai.github.io/assets/img/Kanai2018_Fig1.svg" width="640px"></p>

Fuji plot is named after :mount_fuji: Mt. Fuji in Japan :mount_fuji:, whose ridge and skirts greatly resemble our figure!

<p align="center"><img src="https://i.imgur.com/QXlnKZk.jpg" width="480px"></p>

Requirements

Usage

Please add circos to your PATH or set CIRCOS_PATH in fujiplot.R appropriately.

Rscript fujiplot.R input_example/input.txt input_example/traitlist.txt output_example

Inputs

Input files used in Kanai, M. et al. (2018) are provided under ./input_example.

Significant trait-associated loci (input.txt)

This file provides a concatenated list of all significant trait-associated loci (for each trait). It also defines a unique locus across multiple traits to highlight pleiotropic loci. The required fields are as follows:

Trait list (traitlist.txt)

This file provides a list of traits and their categories. It defines an order of traits and a color of each category in a figure. The required fields are as follows:

Outputs

All output files will be created under ./output.

Circos outputs (circos.{png,svg})

This is the raw output from Circos. Unfortunately, we could not figure out the appropriate way to make our highlights (radial lines) over scatter plots. To this end, we edited the svg file using Adobe Illustrator and produced the final figure.

<p align="center"><img src="output_example/circos.svg" width="640px"></p>

Optional: a bar plot (barplot.pdf)

This corresponds to Fig. 1a of Kanai, M. et al. (2018). Again, you need to do some Illustrator edits to make it pretty. You can suppress this output by setting OUTPUT_BARPLOT = FALSE in fujiplot.R.

Citation

When using Fuji plot, please cite the following paper.

Contact

Masahiro Kanai (mkanai@g.harvard.edu)

http://mkanai.github.io/

Links