A statistical model for unwarping of 1-D electrophoresis gels

A statistical model is proposed which relates density profiles in 1-D electrophoresis gels, such as those produced by pulsed-field gel electrophoresis (PFGE), to databases of profiles of known genotypes. The warp in each gel lane is described by a trend that is linear in its parameters plus a first-order autoregressive process, and density differences are modelled by a mixture of two normal distributions. Maximum likelihood estimates are computed efficiently by a recursive algorithm that alternates between dynamic time warping to align individual lanes and generalised-least-squares regression to ensure that the warp is smooth between lanes. The method, illustrated using PFGE of Escherichia coli O157 strains, automatically unwarps and classifies gel lanes, and facilitates manual identification of new genotypes.