Quantifying the relationship between cell proliferation and morphology during development of the face

Abstract

Morphogenesis requires highly coordinated, complex interactions between cellular processes: proliferation, migration, and apoptosis, along with physical tissue interactions. How these cellular and tissue dynamics drive morphogenesis remains elusive. Three dimensional (3D) microscopic imaging poses great promise, and generates elegant images. However, generating even moderate through-put quantified images is challenging for many reasons. As a result, the association between morphogenesis and cellular processes in 3D developing tissues has not been fully explored. To address this critical gap, we have developed an imaging and image analysis pipeline to enable 3D quantification of cellular dynamics along with 3D morphology for the same individual embryo. Specifically, we focus on how 3D distribution of proliferation relates to morphogenesis during mouse facial development. Our method involves imaging with light-sheet microscopy, automated segmentation of cells and tissues using machine learning-based tools, and quantification of external morphology via geometric morphometrics. Applying this framework, we show that changes in proliferation are tightly correlated to changes in morphology over the course of facial morphogenesis. These analyses illustrate the potential of this pipeline to investigate mechanistic relationships between cellular dynamics and morphogenesis during embryonic development.

Keywords: Convolutional Neural Networks; Developmental biology; Image segmentation; Light-Sheet imaging; Morphometrics; Mouse embryo.

Fig. 1.

Maximum projection intensity images in which Nuclear Green has been used to stain nuclei; phospho-Histone H3 (pHH3) has been stained in magenta to mark cell proliferation. Embryos were harvested at E10–11 from C57Bl/6J mice following timed matings. Embryos were cleared and prepared following the CUBIC protocol, antibody stained, and mounted in 1.5% low-melt agarose. Images were captured on a Zeiss Lightsheet Z1 microscope and analyzed in Arivis. Embryo staging was established by counting tail somites. FNP = Frontonasal Prominence, MxP = Maxillary Prominence, MdP = Mandibular Prominence Representative images of n > 20 samples.

Fig. 2.. Overview of the workflow for creating shape and proliferation atlases from LSFM scans.

Twenty mouse embryos are stained for nuclei (Nuclear Green) and proliferating cells (pHH3). For each specimen, proliferating cell maps and tissues are automatically segmented using CNNs (mesenchyme in teal, neural ectoderm in yellow). For each age group (E10, E10.5, E11, and E11.5), the groupwise affine-registration is performed using the tissue segmentation volume. Then, the resulting transformation for each sample is applied to its proliferating cell map. The shape atlas per age is obtained by majority voting in each voxel. The proliferation atlases are produced by smoothing and normalizing each proliferating cell map, and then averaging the resulting heat maps of the age group.

Fig. 3.

(A) Tissue and proliferation atlases for C57Bl/6J embryos between E10.0 and E11.5. The top row presents the tissue atlases for embryonic ages E10.0, E10.5, E11.0, and E11.5, showing mesenchyme in teal and neural ectoderm in yellow. Voxels are displayed slightly transparent for visualization purposes. The bottom row shows the proliferation map atlases for the same embryonic ages, only in the mesenchyme. (B,D) Principal Component Analysis (PCA) of Procrustes shape coordinates and facial proliferation, respectively, for the whole embryonic dataset. Specimens are coloured by tail somite count, as a proxy for age. Dataset comprises embryos from 8 to 20 somites. (C) Location of anatomical landmarks manually placed on the facial surfaces of embryonic samples in the atlas (average, in red). Dispersion of landmarks for each specimen is also shown and coloured by tail somite count.

Fig. 4.. Segmented cell volume and proliferative fraction along each tissue axis at E10.0 and E10.5 in key tissues.

(A) Segmentation of the maxillary prominence (purple) and frontonasal prominence (orange) on the average tissue atlas. (B) Proliferative heatmaps within segmented tissues in the same view as (A). (C,D) Cell volumes and proliferative cell volumes are compared along each axis and reveal largely homogeneous proliferation in these averaged volumes. The blue line represents the proliferative fraction, or the proportion of proliferating nuclei to total nuclei, and the ribbon represents the standard error between specimens.

Fig. 5.. Voxel-based analysis of relation between proliferation and shape change in the face.

(A) Workflow to relate the shape and proliferation atlases of the previous embryonic age N-1 to the shape of each of the samples belonging to the present embryonic age N. A differential volume (−1,0,+1) for each sample at age N is obtained by comparing to the mesenchyme of the age N-1 atlas (middle row). Then, each voxel of the differential volume at age N was associated with the corresponding voxel of the proliferation atlas of age N-1. Finally, Pearson’s linear correlation coefficient is computed only in a mask corresponding to the mesenchyme of the face. (B) Box-plots showing the proportion of voxels in the fixed mask of the face that presents outward (+1) or inward (−1) remodeling between the consecutive ages studied in this work (E10.0-E11.5). (C) Box-plots presenting the absolute value of the correlation coefficient between the proliferation at age N-1 with the shape change from N-1 to N, in consecutive ages.

Fig. 6.. Voxel-based morphometry for whole embryo heads.

(A) Statistical parametric map showing sagittal slices of craniofacial regions which undergo significant shape change from E10.5 to E11. B. E11 to E11.5. The slice legend is displayed on the bottom right, with each yellow line corresponding to an individual slice in the volume. Hot colours (red-yellow scale) indicate voxel expansion, whereas cold colours (blue-teal scale) represent voxel shrinkage. (C-E) Localizing the relationship between cell proliferation and shape change in anatomical context. (C) Overview of the image analysis workflow. We non-linearly registered each tissue volume to a novel tissue atlas, and use the corresponding transformation to warp the proliferation volume into the atlas space. Next, we calculated the Jacobian determinants of the tissue registration field and visualized significant stage-related shape deviations (t values) with a parametric map. The red-yellow t value scale indicates levels of increasingly large shape change. Next, we exported a volume of t statistics to generate a shape significance mask (i.e., voxels with a t value above the significance threshold), as well as a set of randomly sampled masks (i.e., voxels with randomly sampled t values). The atlas masks are then overlaid (dotted line) onto each warped proliferation volume to relate proliferation to shape change in the anatomical context. (D-E) Density plots comparing the mean proliferation of N=100 randomly sampled masks (bell curves) against the mean proliferation of the significance mask (arrows) for each specimen. E11 and E11.5 are displayed separately.

Fig. 7.. Quantitative relationship between morphology and spatial patterning of proliferation

(A) 2-Block Partial Least Squares (PLS) analysis of Procrustes shape coordinates and proliferation in younger embryos (8–13 somites). Morphs within the scatterplot were constructed from Procrustes coordinates for each specimen. Heatmaps below the scatterplot depict shape changes associated with proliferation, based on the minimum and maximum values along the first latent variable of block 1 in the PLS analysis (shape), which tends toward more medial nasal proliferation. (B) 2-block PLS analysis of Procrustes shape coordinates and proliferation in older embryos (15–20 somites). Morphs were constructed from Procrustes coordinates for each specimen, and heatmaps depict shape changes associated with changes in proliferation along the first axis (latent variable of shape block), which moves more toward the maxilla with increasing age.