During development, organisms acquire three-dimensional shapes with important physiological consequences. While the basic mechanisms underlying morphogenesis are known in eukaryotes, it is often difficult to manipulate them in vivo. To circumvent this issue, here we present a study of developing Vibrio cholerae biofilms grown on agar substrates in which the spatiotemporal morphological patterns were altered by varying the agar concentration. Expanding biofilms are initially flat, but later experience a mechanical instability and become wrinkled. Whereas the peripheral region develops ordered radial stripes, the central region acquires a zigzag herringbone-like wrinkle pattern. Depending on the agar concentration, the wrinkles initially appear either in the peripheral region and propagate inward (low agar concentration) or in the central region and propagate outward (high agar concentration). To understand these experimental observations, we developed a model that considers diffusion of nutrients and their uptake by bacteria, bacterial growth/biofilm matrix production, mechanical deformation of both the biofilm and the agar, and the friction between them. Our model demonstrates that depletion of nutrients beneath the central region of the biofilm results in radially-dependent growth profiles, which in turn, produce anisotropic stresses that dictate the morphology of wrinkles. Furthermore, we predict that increasing surface friction (agar concentration) reduces stress anisotropy and shifts the location of the maximum compressive stress, where the wrinkling instability first occurs, toward the center of the biofilm, in agreement with our experimental observations. Our results are broadly applicable to bacterial biofilms with similar morphologies and also provide insight into how other bacterial biofilms form distinct wrinkle patterns.The intricate shapes of organisms are determined by the spatiotemporal patterns of growth as well as the mechanical properties of their underlying biological components [1-3]. Three-dimensional (3D) shape transformations in developing organisms often arise via differential growth of connected tissues [1, 4]. Such asymmetric growth patterns generate compressive stresses within the faster growing tissues, which may cause mechanical instabilities [5][6][7]. Growth-induced mechanical instabilities drive the formation of many convoluted morphologies, such as the gyrification of brains [2, 8, 9], the vilification and looping of guts [10, 11], and the branching of lungs [12] as well as 3D structures of synthetic systems with patterned swelling [5,[13][14][15][16].Biofilms, which are surface-associated bacterial communities encapsulated by a self-produced extracellular matrix [17, 18], also display a variety of 3D developmental morphologies ranging from radial stripes, to concentric rings, to disordered labyrinth and herringbone patterns [19][20][21][22][23]. In the case of Vibrio cholerae, a model * wingreen@princeton.edu † andrej@princeton.edu biofilm former, quantitative imaging revealed a 3...