Anova for proteomics

Hi,
I’m new with R and I have a dataset with more than 1000 proteins. Every row contains data for the single protein. In the columns I have 3 values for the control, 3 values for Treatment 1, 3 values for Treatment 2 and 3 values for Treatment 3.
I need to perform ANOVA and TukeyHSD between the 4 treatment groups for every protein.
I think I should use for, but I can’t select row by row the elements that I need every time, considering that aov functions for both anova and Tukey requires values belonging to a single group inserted in the same column. Then, I also need to get a table with p-value of anova and Tukey too
Can anyone explain me in a simple way how to do all of that?

Thanks.