Trumpet plots: Visualizing The Relationship Between Allele Frequency And Effect Size In Genetic Association Studies

Using the logarithmic vs linear scale to represent allele frequencies

In the representation of allele frequencies, the range of values can vary greatly between the smallest and largest frequency. When these associations are plotted on a linear scale, the rare variants can be obscured or difficult to distinguish. To address this issue, we recommend using a logarithmic (log) scale (we use log base 10) for allele frequencies in trumpet plots. Compared to a linear scale, the log scale uses increments that represent a relative increase or decrease, rather than a fixed value increase or decrease. The log scale compresses the allele frequencies that are most common, which results in a more even distribution of values across the scale. This scale of visualization facilitates the identification of important patterns and trends.

Identification of independent significant variants to enhance the interpretation of trumpet plots

Genetic association studies involve testing up to millions of genetic variants for their association with a particular trait or disease. However, many of these variants are correlated with each other due to their physical proximity on the genome, which is known as linkage disequilibrium (LD). This means that many variants nearby to causal variants often show significant associations with the trait under study, due only to their correlation rather than any biological involvement with the trait.

Two methods that can be used to identify independent significant variants in a GWAS are clumping and conditional analysis [14], [15]. Clumping involves selecting a subset of independent significant variants by choosing a lead variant for each LD cluster and then discarding all other variants in that cluster. The lead variant is typically the one with the strongest association with the trait of interest. Conditional analysis, on the other hand, involves identifying independent significant variants after performing a joint analysis of multiple variants together. In this approach, the effect of one variant is conditioned on the effect of other variants, meaning that the association between the trait of interest and one variant is evaluated after accounting for the effect of other variants. This can either be performed as a joint analysis (e.g. a regression) with multiple variants in one model or else as an iterative process, where the lead variant and each other variant in the region are tested jointly one-by-one. By considering the effects of multiple variants jointly in single models, conditional analysis helps to identify independent signals more accurately than clumping, which relies only on correlations between variants to indirectly infer independent signals.

Since the correlation between variants can make it challenging to interpret trumpet plots, we recommend plotting only independent significant variants. These are variants that represent distinct genetic associations with the trait of interest.

Variant-level associations for low allele frequencies

While GWAS are a valuable tool for detecting common genetic variants associated with complex traits or diseases, they have limited power to identify associations with rare variants. To address this limitation, sequencing studies [16]–[18], such as whole exome sequencing and whole genome sequencing, have been used to detect variant-level associations with low allele frequencies.

To further enhance the statistical power of rare-variant associations, a commonly used strategy is to aggregate the rare variants detected into functional genetic units, such as genes, and perform collective variation analysis (e.g. gene burden tests [19]–[21]). However, we found that the reporting of rare variant association analyses varies across studies. Some studies report results at the variant level, while others only report results of the functional unit in which rare variants were aggregated. This makes comparing results across independent studies, and determining the functional significance of specific variants, challenging.

We therefore encourage the reporting of results at the variant level, so that they can be included in visualizations of allele frequency in relation to effect size, such as trumpet plots, aiding study comparisons and variant interpretation. Although the analysis of genetic variants at low frequency is expected to improve with the availability of biobank-scale samples and the development of new methods to reduce biases in association tests, caution should still be exercised when inspecting associations with rare variants, as these can suffer from instability and low power, particularly in the relation to binary traits.

Statistical power considerations

GWAS require careful consideration of statistical power [22], which depends on various factors, including the allele frequency and effect size of variants, represented by the x- and y-axis of trumpet plots, respectively. Common variants, usually defined as having allele frequency greater than 1%, tend to have higher power in association studies because common causal variants are more likely to be present in the sample (either genotyped or imputed), and because their relatively balanced number of alleles is akin to having a larger sample size. Variants with larger effect sizes have higher power because their effects are further from the null hypothesis of zero effect.

We therefore recommend incorporating power curves into trumpet plots, since they visually represent the statistical power across the allele frequency spectrum for a given sample size and effect size [23]. Moreover, power curves can aid in identifying parts of the association testing space in which the power to detect significant associations is low.

Two alternative approaches to illustrate the joint distribution of allele frequency and effect size

One approach to illustrate the relationship between allele frequency and effect size is to plot only positive effects i.e the allele effect for each variant that increases the value of the phenotype. In this case, the effect sizes are always positive, and both the allele and sign of the association regression coefficients (betas or Odds Ratios, ORs) need to be flipped (to the other allele) if they are reported as negative to ensure that the effect size is greater than zero. If this ‘flipping’ is required, then the allele frequency of the other allele should be reported, which will be 1 minus the original allele frequency. In this case, the allele frequencies of the plot range from 0 to 1.

The other approach, which we recommend, allows for both positive (risk allele in the context of disease phenotypes) and negative (protective allele in the context of disease phenotypes) effect sizes and always corresponds to the minor allele. In this case, the effect size of the allele can have either a positive or negative value, and the allele frequencies of the plot range from 0 to 0.5.

Practical example: Generating trumpet plots for 129 traits in the UK Biobank

We examined all continuous UK Biobank traits with available GWAS analyses performed by Benjamin Neale’s group (https://www.nealelab.is/uk-biobank/) and searched whether rare variant associations were available for the same trait (by UK Biobank Field ID) in the exome sequencing analysis performed by the Regeneron team [13].

Common variant associations were extracted from the Neale‘s group GWAS summary statistics. For each GWAS, we extracted the independent variants using COJO GCTA (--cojo-slct command), and a random subset of 4,000 unrelated individuals with European ancestry from the UK Biobank as LD reference panel. We selected independent variants with minor allele frequency >0.01 and association P-value < 5×10⁻⁸ within a 100Kb window.

Rare variants association results were extracted from the supplementary data table (SD2) of the Regeneron study [13]. This study reports results for both burden tests (which typically aggregated variants and indels) and individual rare variant level tests. To ensure that effect sizes reported in our analyses corresponded to individual rare variants, we extracted only results for ‘singleton variants’ with predicted loss of function - including stop-gain, frameshift, stop-lost, start-lost and essential splice variants - and deleterious missense variants. Interactive plots illustrating the relationship between allele frequency (x-axis) and odds ratio (y-axis) were plotted using our R package TrumpetPlots (https://gitlab.com/JuditGG/trumpetplots), which uses the R packages data.table, ggplot and ggplotly (Figure 1).

R Shiny application

We developed a user-friendly web application called shinytrumpets to visualize trumpet plots for our UK Biobank results, as well as any other genetic association results that can be uploaded by the user. With shinytrumpets, researchers with no knowledge of R programming can easily upload and visualize their own datasets.

If a user uploads their own results, shinytrumpets prompts them to upload the input data files and specify the sample size used for the study, such as the GWAS sample size. This information is used to perform power calculations for the visualization. Shinytrumpets offers an intuitive interface for users to explore and download trumpet plots.