Research focus: Cell division and elongation are key processes in root development. A growing root tip is assumed to be made of three contiguous zones: the division zone (DZ), the elongation zone (EZ) and the mature zone (MZ). Because of the concentric organization of tissues in roots, we here focus on epidermis cell length profiles for delimiting these three developmental zones. We assume that DZ is characterized by small cells, EZ by cells of gradually increasing length and MZ by rather large cells. These zones were identified on 42 maize lateral roots belonging to three genotypes and, for each genotype, to a rather large diversity in terms of development. Methods: We built two types of Gaussian multiple change-point models for each root: heteroscedastic piecewise Gaussian linear models built on the basis of the measured epidermis cell length profiles and Gaussian change in the variance model built on the basis of the residual series deduced from the optimal piecewise linear function. The number of segments was selected using the slope heuristic, a non-asymptotic penalized likelihood criterion. The outputs of the analysis of each root using multiple change-point models (e.g. lengths of DZ and EZ, cell lengths predicted at the limits between DZ and EZ and between EZ and MZ) combined with some morphological features (first root hair position and root diameter in MZ) were then globally analyzed using principal component analysis (PCA).Results: We identified at most three zones corresponding to DZ, EZ and MZ except in arrested roots where DZ and in some cases EZ were missing. In our modeling framework, the limit between two successive zones corresponds to a marked change in slope and in residual standard deviation. In more than 90% of cases, the limits founds on the residual series on the basis of the sole change in residual standard deviation co-localized with the limit found on the original series. In most cases, the limit between EZ and MZ was consistent with the f