A thick-plate analysis is carried out to first-order in interface slope for three-dimensional folding of a single linear viscous layer of thickness H and viscosity eta embedded in a uniform viscous medium with viscosity eta-1. Layer and medium are subject to a basic state of homogeneous pure shear with principal directions of strain-rate x and y lying in the plane of the layer. The solution is explicitly carried out for a fold perturbation in interface shape that is symmetric with respect to the principal strain-rate axes, A cos (lx) cos (my), but it is shown to apply to an arbitrary perturbation in the x,y plane. The growth rate of the perturbation is found to be dA/dt = -epsilonBAR(zz)A + (1-R){1 - R2)-[(1 + R2)(e(k) - e(-k) + 2R(e(k) + e(-k)]/2k}-1[(l2/lambda-2)-epsilonBAR(xx) + (m2/lambda-2)-epsilonBAR(yy) - epsilonBAR(zz]A, where R = eta-1/eta, k = lambda-H, lambda-2 = l2 + m2, and EBAR(XX), EBAR(YY), and EBAR(ZZ) are the principal strain-rates of the basic state. The wavenumber of the most rapidly growing perturbation, lambda-d, is the same as that obtained for a cylindrical perturbation (m = 0) in a basic state of plane strain (EBAR(YY = 0) . For maximum rate of shortening parallel to x, the cylindrical fold form with axis normal to x, m/l = 0, grows most rapidly for any ratio of EBAR(YY)/EBAR(XX) < 1. If EBAR(YY) = EBAR (XX), all fold forms grow at the same rate and in particular, an 'egg-carton' fold form is not preferentially amplified. The velocity field for the perturbing flow consists of a poloidal field, which solely determines the growth of the fold form, and a toroidal field, with non-zero component of vorticity about the axis normal to the layer. The latter is required to satisfy both of the two independent shear traction continuity conditions at interfaces. There is no coupling between the two fields. Three-dimensional fold forms in the Appalachian Plateau province in western Pennsylvania are described.