How can I add Tukey HSD p-value brackets from a data frame to ggplots from a list of lists?

Question

I have a loop that produces a list of ggplots for any given number of analytes (in this example, 3). An ANOVA is performed and p-value generated and mapped for each ggplot. Howe