The role of a boundary in pattern formation from a homogenous state in Turing's reaction–diffusion equations is important, particularly when the domain size is comparable to the pattern scale. Such experimental conditions may be achieved for
in vitro
regeneration of ectodermal appendages such as feathers, via reconstruction of embryonic single cells. This procedure can eliminate a predefined genetic map, such as the midline of chick feather bud formation, leaving uniformly distributed identical cells as a bioengineered skin. Here, the self-organizing nature of multiple feather bud formation was examined in bioengineered 1D-skin samples. Primal formation of feather buds occurred at a fixed length from the skin edge. This formation was numerically recapitulated by a standard two-component reaction-diffusion model, suggesting that the boundary effect caused this observation. The proper boundary conditions were nonstandard, either mixed Dirichlet–Neumann or partial-flux. In addition, the model implies imperfect or hindered bud formation as well as nearly equal distances between buds. In contrast, experimental observations indicated that the skin curvature, which was not included in our model, also strongly affected bud formation. Thus, bioengineered skin may provide an ideal template for modeling a self-organized process from a homogenous state. This study will examine the possible diffusion activities of activator or inhibitor molecular candidates and mechanical activities during cell aggregation, which will advance our understanding of skin appendage regeneration from pluripotent or embryonic stem cells.