Abstract
Recent advances in genome-wide association study (GWAS) and sequencing studies have shown that the genetic architecture of complex diseases and traits involves a combination of rare and common genetic variants, distributed throughout the genome. One way to better understand this architecture is to visualize genetic associations across a wide range of allele frequencies. However, there is currently no standardized or consistent graphical representation for effectively illustrating these results.
Here we propose a standardized approach for visualizing the effect size of risk variants across the allele frequency spectrum. The proposed plots have a distinctive trumpet shape, with the majority of variants having low frequency and small effects, while a small number of variants have higher frequency and larger effects. These plots, which we call ‘trumpet plots’, can help to provide new and valuable insights into the genetic basis of traits and diseases, and can help prioritize efforts to discover new risk variants. To demonstrate the utility of trumpet plots in illustrating the relationship between the number of variants, their frequency, and the magnitude of their effects in shaping the genetic architecture of complex diseases and traits, we generated trumpet plots for more than one hundred traits in the UK Biobank. To facilitate their broader use, we have developed an R package ‘TrumpetPlots’ and R Shiny application, available at https://juditgg.shinyapps.io/shinytrumpets/, that allows users to explore these results and submit their own data.
STATEMENT OF NEED Visualizations are powerful tools that have helped the field of genetics to better understand and communicate complex findings. By using visual aids like Manhattan and Volcano plots, genetic variants identified through genome-wide association studies can be more easily pinpointed. With the advancement of genome-wide association and sequencing studies, a mounting number of significant genetic variants, both common and rare, are being discovered. To better understand the relationship between these variants, combining these findings into single visualizations help to observe the relationship between effect size and allele frequency, providing a clearer picture of the genetic architecture of different traits and diseases. However, there is currently no consistent method for illustrating these results. In this paper, we propose a standardized approach for visualizing the effect size of risk variants across the allele frequency spectrum, generate plots for over a hundred traits in the UK Biobank, and provide to the field a R package and R Shiny application to explore their own results.
Availability of supporting source code and requirements
Project name:
o. R package available in project ‘TrumpetPlots’ https://gitlab.com/JuditGG/trumpetplots
o. R shiny app and analyses in the UK Biobank available in project ‘freq_or_plots’ https://gitlab.com/JuditGG/freq_or_plots
Project home page: https://juditgg.shinyapps.io/shinytrumpets/
Operating system(s): Platform independent.
Programming language: R
RRID: Not applicable
License: MIT
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
This work was supported by a grant from the National Institute of Health (R01MH122866) to PFO, by a 2022 NARSAD Young Investigator Grant (Number 30749) by the Brain & Behavior Research Foundation to JGG, and through the computational resources and staff expertise provided by Scientific Computing and the Data Ark (Data Commons) teams at the Icahn School of Medicine at Mount Sinai.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
All data used in this manuscript is publicly available. Rare variant associations are available in supplementary data table 2 of the original publication (DOI: https://doi.org/10.1038/s41586-021-04103-z), GWAS summary statistics are available in the website https://www.nealelab.is/uk-biobank/.
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines, such as any relevant EQUATOR Network research reporting checklist(s) and other pertinent material, if applicable.
Yes
Footnotes
Email: lc105{at}wellesley.edu
Email: lathan.liou{at}icahn.mssm.edu
email: paul.oreilly{at}mssm.edu
email: judit.garciagonzalez{at}mssm.edu
Data Availability
All data used in this manuscript is publicly available. Rare variant associations are available in supplementary data table 2 of the original publication (DOI: https://doi.org/10.1038/s41586-021-04103-z ), GWAS summary statistics are available in the website https://www.nealelab.is/uk-biobank/. The code is freely available at https://gitlab.com/JuditGG/freq_or_plots (UK Biobank analyses), https://gitlab.com/JuditGG/trumpetplots (R package with test data) and https://juditgg.shinyapps.io/shinytrumpets/ (R Shiny application).
List of abbreviations
- GWAS
- genome-wide association study
- LD
- linkage disequilibrium
- Log
- logarithm or logarithmic
