Bacteria picked up on fingers after multiple surface contacts: non-normal data, repeated measures, crossed participants Intro
I have participants who are repeatedly touching contaminated surfaces with E. coli in two conditions (A=wearing gloves, B=no gloves). I want to know if there's a difference between the amount of bacteria on their fingertips with and without gloves but also between the number of contacts. Both factors are within-participant.
Experimental Method:
Participants (n=35) touch each square once with the same finger for a maximum of 8 contacts (see figure a).

I then swab the finger of the participant and measure the bacteria on the fingertip after each contact. They then use a new finger to touch a different number of surfaces and so on from 1 to 8 contacts (see figure b).
Here is the real data: real data
The data is non-normal so see marginal distribution of bacteria|NumberContacts below. x=bacteria. Each facet is a different number of contacts.

MODEL
Trying  from lme4::glmer based on amoeba's suggestions using Gamma(link="log") and polynomial for NumberContacts:
cfug<-glmer(CFU ~ Gloves + poly(NumberContacts,2) + (-1+NumberContacts|Participant),
            data=(K,CFU<4E5),
           family=Gamma(link="log")
            )
plot(cfug)

NB. Gamma(link="inverse") won't run saying PIRLS step-halving failed to reduce deviance.
Results:
Fitted vs residuals for cfug


qqp(resid(cfug))


Question:
Is my glmer model properly defined to incorporate the random effects of each participant and the fact that everyone does both experiment A followed by experiment B?
Addition:
Autocorrelation seems to exist between participants. This is probably because they weren't tested on the same day and the flask of bacteria grows and declines over time. Does it matter?
acf(CFU,lag=35) shows an significant correlation between one participant and the next.

 A: Firstly, good work on your graph; it gives a clear representation of the data, so you can already see the kind of pattern in the data based on the number of contacts and the use or absence of gloves.$^\dagger$  Looking at this graph, I would think you will get good results with a basic log-polynomial model, with random effects for the participants.  The model you have chosen looks reasonable, but you might also consider adding a quadratic term for the number of contacts.
As to whether to use MASS:glmmPQL or lme4:glmer for your model, my understanding is that both these functions will fit the same model (so long as you set the model equation, distribution and link function the same) but they use different estimation methods to find the fit.  I could be mistaken, but my understanding from the documentation is that glmmPQL uses penalised quasi-likelihood as described in Wolfinger and O'Connell (1993), whereas glmer uses Gauss-Hermite quadrature.  If you are worried about it you can fit your model with both methods and check that they give the same coefficient estimates and that way you will have greater confidence that the fitting algorithm has converged to the true MLEs of the coefficients.


Should NumberContacts be a categorical factor?

This variable has a natural ordering that appears from your plots to have a smooth relationship with the response variable, so you could reasonably treat it as a numeric variable.  If you were to include factor(NumberContacts) then you will not constrain its form and you will not lose many degrees-of-freedom.  You could even use the interaction Gloves*factor(NumberContacts) without losing too many degrees-of-freedom.  However, it is worth considering whether using a factor variable would involve over-fitting the data.  Given that there is a fairly smooth relationship in your plot a simple linear function or quadratic would get good results without over-fitting.


How do you make Participant a random slope but not intercept variable?

You have already put your response variable on a log-scale by using a logarithmic link function, so an intercept effect for Participant is giving a multiplicative effect on the response.  If you were to give this a random slope interacting with NumberContacts then it would have a power-based effect on the response.  If you want this then you can get it with (~ -1 + NumberContacts|Participant) which will remove the intercept but add a slope based on the number of contacts.


Should I use Box-Cox for transforming my data? (e.g. lambda=0.779)

If in doubt, try fitting a model with this transformation and see how it compares to other models using appropriate goodness-of-fit statistics.  If you are going to use this transformation it would be better to leave the parameter $\lambda$ as a free parameter and let it be estimated as part of your model, rather than pre-specifying a value.


Should I include weights for variance?

Start by looking at your residual plot to see if there is evidence of heteroscedasticity.  Based on the plots you have already included it looks to me like this is not an issue, so you don't need to add any weights for the variance.  If in doubt you could add weights using a simple linear function and then perform a statistical test to see if the slope of the weighting is flat.  That would amount to a formal test of heteroscedasticity, which would give you some backup for your choice.


Should I include autocorrelation in NumberContacts?

If you have already included a random effect term for the participant then it would probably be a bad idea to add an auto-correlation term on the number of contacts.  Your experiment uses a different finger for different numbers of contacts so you would not expect autocorrelation for the case where you have already accounted for the participant.  Adding an autocorrelation term in addition to the participant effect would mean that you think there is a conditional dependence between the outcome of different fingers, based on the number of contacts, even for a given participant.

$^\dagger$ Your graph is good in showing the relationships, but you could improve it aesthetically by adding a title and subtitle information and giving it better axis labels.  You could also simplify your legend by removing its title and changing 'Yes' to 'Gloves' and 'No' to 'No Gloves'.
A: Indeed, it is reasonable to argue that measurements taken from one participant are not independent from those taken from another participant. For example, some people might tend to press their finger with more (or less) force, which would affect all their measurements accross each number of contacts.
So the 2-way repeated-measures ANOVA would be an acceptable model to apply in this case.
Alternatively, one could also apply a mixed-effects model with participant as a random factor. This is a more advanced and a more sophisticated solution.
