Facial development requires a complex and coordinated series of cellular events that, when perturbed, can lead to structural birth defects. A quantitative approach to quickly assess morphological changes could address how genetic or environmental inputs lead to differences in facial shape and promote malformations. Here, we report on a method to rapidly analyze craniofacial development in zebrafish embryos using facial analytics based on a coordinate extrapolation system, termed zFACE. Confocal images capture facial structures and morphometric data are quantified based on anatomical landmarks present during development. The quantitative morphometric data can detect phenotypic variation and inform on changes in facial morphology. We applied this approach to show that loss of smarca4a in developing zebrafish leads to craniofacial anomalies, microcephaly and alterations in brain morphology. These changes are characteristic of Coffin–Siris syndrome, a rare human genetic disorder associated with mutations in SMARCA4. Multivariate analysis of zFACE data facilitated the classification of smarca4a mutants based on changes in specific phenotypic characteristics. Together, zFACE provides a way to rapidly and quantitatively assess the impact of genetic alterations on craniofacial development in zebrafish.

Vertebrate craniofacial development requires the complex orchestration of cellular processes, molecular signals and interactions between different cell types (Chai and Maxson, 2006; Ji et al., 2020). During the course of organ development, interactions between different tissues are critical in the regulation of cellular proliferation, migration and apoptosis to develop the intricate features that comprise the face (Chai and Maxson, 2006; Compagnucci et al., 2021). Failure to efficiently coordinate the specification, migration, proliferation and apoptosis of cells can lead to craniofacial malformations and structural birth defects, including cleft lip and/or palate, craniosynostosis and facial dysostosis (Hennekam et al., 2010; Van Otterloo et al., 2016). Many developmental syndromes, such as Stickler syndrome, Van der Woude syndrome and Coffin–Siris syndrome (CSS), include craniofacial anomalies in addition to other congenital malformations (Gorlin et al., 2011; Twigg and Wilkie, 2015). One impediment to understanding how genetic variants promote craniofacial anomalies has been our ability to visualize the complex and coordinated cellular interactions during craniofacial development in vivo.

Animal models, such as mouse, chicken, frog and zebrafish, provide an important avenue to gain better mechanistic and genetic insights into human craniofacial development (Van Otterloo et al., 2016). Zebrafish offer many advantages as a model system that make them ideal for detailed craniofacial studies – they develop externally and generate large numbers of transparent embryos, which permits unparalleled high-resolution imaging of specific cell types and structures in living vertebrates (Eames et al., 2013; Nichols et al., 2013; Mork and Crump, 2015; Eisenhoffer et al., 2017; Choe et al., 2021; Shull et al., 2022). Zebrafish craniofacial development, particularly craniofacial bone and cartilage formation, has been well characterized and is comparable to amniote craniofacial development. For example, the development of the anterior neurocranium/ethmoid plate is thought to be functionally analogous to palate development in mammals (Swartz et al., 2011; Dougherty et al., 2013). Zebrafish have been used to systematically test the function of genes associated with birth defects in humans (Eberhart et al., 2008; Yuan et al., 2012; Bhatia et al., 2015; Deml et al., 2015; Swindell et al., 2015; Küry et al., 2017; Machado and Eames, 2017; Shaw et al., 2017; Teng et al., 2018; Fernandes and Lovely, 2021; Truong and Artinger, 2021), and an array of conserved genetic variants/mutations already exist that display altered craniofacial development (Amsterdam et al., 2004; McCarthy et al., 2013). Yet, a standardized, unbiased method that can be used to quantitatively assess phenotypic changes in craniofacial structures resulting from these mutations is not currently available.

Geometric morphometrics (GMM) is a quantitative method used to measure and statistically test for variation(s) in shape (Hallgrimsson et al., 2015). GMM methods have been applied to skeletal and soft tissue data to evaluate craniofacial morphology in the clinical setting and to better understand human craniofacial disorders, such as cleft lip and palate, craniosynostosis, ectodermal dysplasia and neurodevelopmental disorders (Weinberg et al., 2009; Goodwin et al., 2014; Cocos and Halazonetis, 2017; Katsadouris and Halazonetis, 2017; Kesterke et al., 2018; Mercan et al., 2020; Rutland et al., 2021). The presentation of craniofacial abnormalities and syndromes often display phenotypic heterogeneity, and the use of GMM provides an important approach that captures critical information that would otherwise be missed using only a discrete categorical variable such as the presence or absence of a normal or abnormal phenotype (Hallgrimsson et al., 2015). For example, GMM has been used to identify facial differences in unaffected relatives of individuals with non-syndromic cleft lip and palate and strain-specific differences in embryonic facial shape underlying susceptibility to developing orofacial clefts in mice (Young et al., 2007; Weinberg et al., 2009).

To leverage the power of the GMM approach for zebrafish embryos, we developed a novel and easily implemented method, facial analytics based on a coordinate extrapolation system (zFACE), to visualize the developing rostral/frontal craniofacial region and analyze quantitative facial morphometric data in a semi-automated way. We applied zFACE to show that loss of smarca4a modifies craniofacial morphology in zebrafish embryos and identified regional differences that contribute to the observed altered facial dimensions. These results support smarca4a mutant zebrafish as an animal model to provide insight into the craniofacial abnormalities associated with CSS, a rare genetic disorder associated with mutations in SMARCA4 (BRG1). Together, our GMM-based approach and its application demonstrate the power of zFACE to extend our understanding of phenotypic changes in craniofacial development associated with genetic alterations in a vertebrate embryo.

zFACE

Zebrafish embryos develop externally and provide a unique opportunity to visualize development (Mork and Crump, 2015; Eisenhoffer et al., 2017; Machado and Eames, 2017). Although zebrafish craniofacial development has been well characterized, most studies have used lateral and ventral views of the developing face. To visualize additional craniofacial anatomical structures, we first established a simple imaging paradigm using 4′,6-diamidino-2-phenylindole (DAPI), a stain that labels the nucleus of every cell, along with a face-on rostral mounting method using low-melt agarose (Fig. 1A-E). This method can be easily adapted to accommodate different developmental stages, when the orientation of the face with respect to the rest of the body varies, by adjusting the mounting angle of the specimen. However, it is important to mount specimens as consistently as possible in the same orientation for a given developmental stage or mutant. The resulting images capture structures such as the neuromasts, olfactory placodes, eyes, oral cavity, and frontonasal, maxillary and mandibular regions (Fig. 1D,E). Thus, this rostral mounting paradigm provides high-resolution images that reveal important facial information compared to the commonly used lateral and ventral views (Fig. 1B,C).

Fig. 1.

Craniofacial features visualized in frontal/rostral mounts by DAPI staining and confocal microscopy. (A-D) Rostral (A,D), lateral (B) and ventral (C) views of developing zebrafish larvae at 5 days post-fertilization (dpf). (E) Schematic representing the facial anatomical structures visualized in the rostral view, including neuromasts, olfactory placodes, eyes, mouth, frontonasal region (FNP), maxillary region (MX) and mandibular region (MD). Scale bars: 50 µm (A-C), 20 µm (D).

Fig. 1.

Craniofacial features visualized in frontal/rostral mounts by DAPI staining and confocal microscopy. (A-D) Rostral (A,D), lateral (B) and ventral (C) views of developing zebrafish larvae at 5 days post-fertilization (dpf). (E) Schematic representing the facial anatomical structures visualized in the rostral view, including neuromasts, olfactory placodes, eyes, mouth, frontonasal region (FNP), maxillary region (MX) and mandibular region (MD). Scale bars: 50 µm (A-C), 20 µm (D).

We applied GMM to quantitatively analyze facial form based on this newfound ability to identify facial features and landmarks using zFACE. Based on the information captured by the confocal images, 26 easily identifiable landmarks during the embryonic and early larval stages of development were established (Fig. 2; Table S1). Given the complex shape of the oral cavity, and that many craniofacial abnormalities affect the mouth, approximately one-third of the landmarks were assigned to this area to provide better resolution of changes in morphology for this region. Additionally, landmarks were chosen and named to parallel those used in human GMM studies (Weinberg et al., 2009; Liang et al., 2013). Using these established landmarks, an automated calculation of 39 different linear distances, angles and areas was generated to determine the localization of landmarks relative to each other. As repeatability of landmark placement is critical for downstream analysis and comparison, we provide 28 images of 5 days post-fertilization (dpf) wild-type (WT) AB larvae for users to assess consistency among each user, as well as variation between different users (see https://doi.org/10.6084/m9.figshare.22732166.v1) (Klingenberg and McIntyre, 1998; Fruciano, 2016; Špelić et al., 2021; Vrdoljak et al., 2021). Additional statistical approaches were then applied to evaluate overall facial shape (Fig. 2; Table S1). Together, these analyses can be used to quantify phenotypic information resulting from genetic perturbation or environmental exposures and provide valuable information about craniofacial form.

Fig. 2.

Facial analytics based on a coordinate extrapolation system (zFACE) workflow for morphometric analysis of the developing zebrafish craniofacial region. Schematic representing the steps in zFACE analysis, in which DAPI-stained larvae are mounted, images are acquired, points are placed on the 26 anatomical landmarks, and either feature-based zFACE calculations or overall shape analyses can be performed.

Fig. 2.

Facial analytics based on a coordinate extrapolation system (zFACE) workflow for morphometric analysis of the developing zebrafish craniofacial region. Schematic representing the steps in zFACE analysis, in which DAPI-stained larvae are mounted, images are acquired, points are placed on the 26 anatomical landmarks, and either feature-based zFACE calculations or overall shape analyses can be performed.

Morphometric analysis of facial development in zebrafish

To assess the ability to detect changes in specific anatomical structures and locations, we used zFACE to investigate facial development over time. Images of rostrally mounted zebrafish embryos from 48 h post-fertilization (hpf) to 6 dpf were acquired by confocal microscopy and compared to scanning electron microscope (SEM) images. The morphology of the soft tissue and anatomical structures was similar in both conditions, indicating that fixation and the DAPI stain, mounting technique or confocal capture does not induce artifacts in facial form (Fig. 3A). The resulting confocal and SEM images revealed the emergence of specific anatomical structures and regional changes that occurred at defined times throughout development. For example, between 2 and 3 dpf, the oral cavity expands into a semi-circle morphology, while the olfactory placodes shift ventrally. Between 3 and 4 dpf, the biggest changes include a narrowing of the face and an enlargement of the mouth, while between 4 and 5 dpf, the midface widens and the oral cavity becomes crescent shaped. Comparison of 5 and 6 dpf confocal images did not discern any major changes in morphology. All of the 26 anatomical landmarks used by zFACE were not present until 3 dpf, and, therefore, quantitative analysis included the course of development from 3 to 6 dpf.

Fig. 3.

Time series of zebrafish craniofacial development. (A) Zebrafish larvae were analyzed from 2 to 6 dpf. Scanning electron microscope (SEM) images confirmed the soft tissue morphology captured by DAPI staining and confocal microscopy used in zFACE. (B) Changes in zFACE features were followed across developmental time points, and alterations in 28 measurements were determined (representative changes are shown in the graphs). The numbers of individual larvae for each group were as follows: n=30 at 3 dpf, n=31 at 4 dpf, n=28 at 5 dpf and n=28 at 6 dpf. Specifically, 20 zFACE measurements changed between 3 and 4 dpf, 12 changed between 4 and 5 dpf, and 6 nominally changed between 5 and 6 dpf. *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001; n.s., not significant (one-way ANOVA). (C) Untransformed zFACE features were subjected to multivariate statistical testing using principal component analysis (PCA). The first two components captured 66% of the variance in the dataset. There were differences in developmental times that were captured across PC1, with no strain-specific differences observed. Scale bars: 50 µm.

Fig. 3.

Time series of zebrafish craniofacial development. (A) Zebrafish larvae were analyzed from 2 to 6 dpf. Scanning electron microscope (SEM) images confirmed the soft tissue morphology captured by DAPI staining and confocal microscopy used in zFACE. (B) Changes in zFACE features were followed across developmental time points, and alterations in 28 measurements were determined (representative changes are shown in the graphs). The numbers of individual larvae for each group were as follows: n=30 at 3 dpf, n=31 at 4 dpf, n=28 at 5 dpf and n=28 at 6 dpf. Specifically, 20 zFACE measurements changed between 3 and 4 dpf, 12 changed between 4 and 5 dpf, and 6 nominally changed between 5 and 6 dpf. *P<0.05, **P<0.01, ***P<0.001, ****P<0.0001; n.s., not significant (one-way ANOVA). (C) Untransformed zFACE features were subjected to multivariate statistical testing using principal component analysis (PCA). The first two components captured 66% of the variance in the dataset. There were differences in developmental times that were captured across PC1, with no strain-specific differences observed. Scale bars: 50 µm.

The resulting zFACE measurements were compared across developmental time points and strains. In total, 23 measurements showed a significant change between consecutive developmental day comparisons (3 and 4 dpf, 4 and 5 dpf, or 5 and 6 dpf). The most changes were observed between 3 and 4 dpf, when 20 measurements, including facial width, mouth width and height, olfactory to mouth angles, as well as multiple neuromast measurements, were observed to be significantly different (Table 1, Fig. 3B). Comparison of 4 and 5 dpf measurements revealed 12 significantly different measurements, which included a further increase in mouth width and a decrease in mouth height, as well as more pronounced changes in labiale inferius, crista philtri and labiale superius angles (Table 1, Fig. 3B). Additionally, when later-stage (5 and 6 dpf) larvae were compared, no measurements met the Bonferroni-corrected P-value cutoff (six were nominally significant), suggesting that facial morphology was very similar between these two time points (Table 1, Fig. 3B). Interestingly, eight measurements showed a daily increase, ten showed a daily decrease, and five showed a non-linear change (Table 1). Comparison of facial development between the two most commonly used laboratory zebrafish strains, AB and TU, showed relatively few significant differences (Table S2) among the 39 different measurements, suggesting that the strains are very similar in both the timing of anatomical changes and overall facial morphology (Suurväli et al., 2020). This analysis identified time points in developing zebrafish when facial dimensions are most dynamic and ones when face shape remains more stable.

Table 1.

Facial analytics based on a coordinate extrapolation system (zFACE) results, showing phenotypic differences in craniofacial morphology between developmental days

Facial analytics based on a coordinate extrapolation system (zFACE) results, showing phenotypic differences in craniofacial morphology between developmental days
Facial analytics based on a coordinate extrapolation system (zFACE) results, showing phenotypic differences in craniofacial morphology between developmental days

To understand which zFACE measurements contribute the most variance throughout development, we performed multivariate principal component analysis (PCA) of the combined AB and TU groups for all time points. Six principal components met the Kaiser cutoff and collectively accounted for 86% of the variation in the dataset (Fig. 3C). Examination of the principal component (PC) loadings was then used to identify the measurements driving variance along each component. Results showed that neuromast width significantly loaded into PC1, whereas mouth height, mouth area, labiale inferius angle, and left and right crista philtri angles significantly loaded into PC2. The labiale superius and left and right chelion angles significantly loaded into PC3 (Table S2). These data suggested that measurements associated with the morphology of the mouth significantly contributed to the differences between developmental stages. Further, examination of PC plots for the first two components revealed clustering of groups by developmental day along PC1 (Fig. 3C), while no strain-specific clusters were observed. These results indicate that unbiased analyses can also be used to understand how facial morphology changes during early craniofacial development and support the results obtained from the automated zFACE calculations.

We next applied shape analysis in MorphoJ (Klingenberg, 2011) to account for the size differences that are expected to vary between developmental time points. A Procrustes superimposition, which transforms shapes so that they are in maximal superimposition, was applied to the zFACE data in an attempt to remove time-dependent variation due to size, position or rotation, and the PCA was repeated; the first four PCs now cumulatively explained 92% of the variance (the first two PCs explained 81%) (Fig. S1A). Examination of the PC plot for the first two components showed very similar results to the PCA using untransformed data, with developmental days varying across PC1 and no strain-based clusters in the data (Fig. S1A). To focus on how facial shape changed as development progressed, discriminant function analysis (DFA) was utilized, and results were summarized by overlaid wireframe representations (Fig. S1B-D). This analysis revealed significant shape changes between 3 and 4 dpf (Procrustes distance=0.14, P<0.0001) and 4 and 5 dpf (Procrustes distance=0.12, P=0.005), but no changes between 5 and 6 dpf (Procrustes distance=0.05, P=0.04 before Bonferroni correction) (Fig. S1B-D). No strain-specific shape differences were observed at any of the developmental time points (Fig. S1E-G). Together, results from these analyses suggest that zFACE represents a robust and sensitive approach for morphometric analysis of facial development in zebrafish embryos.

Application of zFACE to smarca4a mutant zebrafish reveals similarities with CSS

To test the application of zFACE and assess its utility for detection of variation in morphology that occurs after genetic perturbation, we analyzed zebrafish with loss of smarca4aa8, which have been described to have craniofacial anomalies (Link et al., 2000; Gregg et al., 2003; Eroglu et al., 2006). Our rostral confocal images revealed a severely affected facial phenotype in smarca4a homozygous mutant larvae at 5 dpf, whereas heterozygotes and WT larvae showed normal developmental hallmarks (Fig. 4A-C; Fig. S2). smarca4a mutants had a narrow head and face, smaller brain, reduced olfactory pits, open and elongated oral cavities, and small mandibles (Fig. 4C). Automated calculation and comparison of zFACE measurements showed that homozygous mutants significantly differed from heterozygotes and WT controls in 33 out of the 39 zFACE measurements, including reduced facial width, increased facial height, decreased olfactory distance, reduced mouth width and altered oral cavity angles (Fig. 4D; one-way ANOVA, post hoc Tukey's test, P<0.0013 for all; Table 2). All three groups were equal in upper lip width, chin width, mid olfactory to chin height, mouth area, and the difference between chelion and labiale inferius left and right angles, suggesting that smarca4a does not affect these structures.

Fig. 4.

zFACE analysis of zebrafish with loss of smarca4a. (A-C) Confocal images of wild-type (WT; A), heterozygous (Het; B) and homozygous mutant (Mut; C) smarca4a larvae. (D) zFACE analysis pinpointed 33 significantly altered facial measurements between the groups, with examples of dimensions that were unaltered, decreased and increased shown. The numbers of larvae analyzed for each group were n=28 for WT, n=34 for Het and n=32 for Mut. ****P<0.0001; n.s., not significant (one-way ANOVA). (E) Untransformed zFACE features were subjected to multivariate statistical testing using PCA. The first two components explained 68% of the variance in the dataset, and the principal component (PC) plot showed clear separation of the homozygous mutant group from the WT and heterozygote groups. (F) Discriminant function analysis after landmark data were Procrustes transformed further supported that smarca4a homozygous mutants had different average facial shapes compared to either WT or heterozygous larvae, whereas WT and heterozygous groups had the same facial shape. Wireframe representations of facial shapes highlight the alterations to landmarks in the upper face, eyes, oral cavity, midface and lower jaw in the homozygous smarca4a mutants [discriminant function analysis (DFA)]. Scale bar: 50 µm.

Fig. 4.

zFACE analysis of zebrafish with loss of smarca4a. (A-C) Confocal images of wild-type (WT; A), heterozygous (Het; B) and homozygous mutant (Mut; C) smarca4a larvae. (D) zFACE analysis pinpointed 33 significantly altered facial measurements between the groups, with examples of dimensions that were unaltered, decreased and increased shown. The numbers of larvae analyzed for each group were n=28 for WT, n=34 for Het and n=32 for Mut. ****P<0.0001; n.s., not significant (one-way ANOVA). (E) Untransformed zFACE features were subjected to multivariate statistical testing using PCA. The first two components explained 68% of the variance in the dataset, and the principal component (PC) plot showed clear separation of the homozygous mutant group from the WT and heterozygote groups. (F) Discriminant function analysis after landmark data were Procrustes transformed further supported that smarca4a homozygous mutants had different average facial shapes compared to either WT or heterozygous larvae, whereas WT and heterozygous groups had the same facial shape. Wireframe representations of facial shapes highlight the alterations to landmarks in the upper face, eyes, oral cavity, midface and lower jaw in the homozygous smarca4a mutants [discriminant function analysis (DFA)]. Scale bar: 50 µm.

Table 2.

Results of zFACE feature comparison between wild-type, heterozygous and homozygous mutant smarca4a zebrafish embryos

Results of zFACE feature comparison between wild-type, heterozygous and homozygous mutant smarca4a zebrafish embryos
Results of zFACE feature comparison between wild-type, heterozygous and homozygous mutant smarca4a zebrafish embryos

Dimensionality reduction via multivariate PCA of zFACE measurements identified six principal components that met the Kaiser cutoff and cumulatively explained 84% of the variance (Fig. 4E). The first component (PC1) explained 58% of the variance, PC2 explained an additional 11%, and PC3 explained 5%. After varimax rotation, component loadings showed that the measurements neuromast height and area top significantly loaded into PC1, whereas width, neuromast width, average length of olfactory to mouth, and facial area combined loaded into PC2. The score plot (PC plot) for these first two components showed smarca4a homozygous mutants as having higher PC1 scores and separate from the WT and heterozygous groups (Fig. 4E). Using the PCA model, PC scores were predicted for each embryo, and logistic regression was performed with these predicted scores using group as the dependent variable and PC score as the independent variable. Importantly, PC1 score alone could predict whether an embryo was a homozygous smarca4a mutant (when PC1 score was greater than −0.98) but could not distinguish between WT or heterozygote embryos (P=0.86) (Fig. 4E). Because the zFACE measurements are composed of different types of units (angles, areas, distances and differences) and the standard deviation (s.d.) between the measurements is not equal, PCA was also performed after standardizing the data and scaling it to have a mean of 0 and s.d. of 1. This led to very similar PCA results, with six PCs cumulatively explaining 88% of the variability. Additionally, we compared 4 and 6 dpf WT embryos (1 day earlier and 1 day later) to the smarca4a embryos to examine whether the morphological differences in the mutants could be due to delayed facial development. In the PCA plot, all WT embryos plotted separately from the smarca4a homozygous mutants, suggesting that development is not simply delayed and that the abnormal facial phenotype is specific to smarca4a disruption (Fig. S3).

Lastly, we performed Procrustes superimposition of the landmark coordinates in MorphoJ (Klingenberg, 2011). Data were analyzed based on the assumption of object symmetry in the head, and the symmetric portion was evaluated in all analyses (Klingenberg et al., 2002). PCA resulted in PC1 accounting for 85% of the variance, PC2 4% and PC3 another 3%; cumulatively these first three PCs explained 92% of the variance (Fig. S4A). Similar to results from the untransformed PCA score plot, the genotype groups could be clearly distinguished by graphing PC1 versus PC2 (Fig. S4B). Shape changes across PC1 affected mouth sphericity, with landmarks around the oral cavity showing the biggest changes (highest eigenvectors); changes across PC2 are suggestive of a narrower midface, and elongated head and mouth (Fig. S4C).

To compare the three genotype groups in an unbiased manner, we used canonical variate analysis (CVA) as an exploratory method. The results showed significant Mahalanobis and Procrustes distance differences among groups, with the smarca4a homozygotes significantly differing from the heterozygote and WT groups (Mahalanobis distance=19.85 and Procrustes distance=0.33; P<0.0001 for both when homozygotes were compared to WT; and Mahalanobis distance=19.22 and Procrustes distance=0.33; P<0.0001 for both when homozygotes were compared to heterozygotes). Changes in canonical variate 1 (CV1)-involved landmarks around the mouth led to a more round and open mouth, whereas those along CV2 were associated with a narrower and elongated face, similar to the results obtained from the PCA analysis (Fig. S4). Because the WT and heterozygous groups were so similar, DFA, which is similar to CVA but compares only two groups at a time, was performed to focus on the shape changes specific to the homozygous mutants. Results showed a Procrustes distance of 0.33 (P<0.0001) in both the WT and heterozygote to homozygous mutant comparisons (Fig. 4F), and wireframe representations depicted shape change involving a vertically elongated mouth (Fig. 4F), reflecting the phenotype seen in the confocal images (Fig. 4C; Fig. S2).

To better understand the observed phenotypic abnormalities, we performed fluorescent in situ hybridization to examine smarca4a mRNA expression in the craniofacial region. Expression of smarca4a mRNA was observed at 5 dpf in several facial tissues, including the perioral and oral tissues, olfactory placodes and the telencephalon in both WT and heterozygous larvae (Fig. 5A,B). smarca4a−/− mutants, however, showed a dramatically reduced signal in the entire craniofacial region, suggesting that the point mutation (which creates a premature stop codon) leads to no detectable mRNA expression in developing facial structures (Fig. 5C,D) (Gregg et al., 2003). This expression pattern of smarca4a mRNA was consistent with the facial structures that zFACE identified as the most altered in the smarca4a mutants. Based on the observed morphological changes in brain and neural tissues, the width and length of the brain were measured, and mutants showed an increased length-to-width ratio (one-way ANOVA, P<0.0001) (Fig. S5). Together, these results suggest similarities between smarca4a mutants and individuals with CSS, who often present with craniofacial abnormalities such as a smaller mouth, thicker lips and reduced philtrum, as well as intellectual disability and microcephaly (Coffin and Siris, 1985; Tsurusaki et al., 2012).

Fig. 5.

smarca4a mRNA expression at 5 dpf. (A-D) Fluorescent in situ hybridization showed that smarca4a mRNA is expressed in multiple facial tissues and underlying brain structures in WT and heterozygous zebrafish (A,B), whereas smarca4a/− mutants showed dramatically reduced smarca4a mRNA expression (C,D). Scale bars: 50 µm (A,C), 20 µm (B,D).

Fig. 5.

smarca4a mRNA expression at 5 dpf. (A-D) Fluorescent in situ hybridization showed that smarca4a mRNA is expressed in multiple facial tissues and underlying brain structures in WT and heterozygous zebrafish (A,B), whereas smarca4a/− mutants showed dramatically reduced smarca4a mRNA expression (C,D). Scale bars: 50 µm (A,C), 20 µm (B,D).

GMM is a powerful quantitative approach for assessing phenotypic differences arising from alterations in shape and size. Here, we developed a streamlined method called zFACE to visualize craniofacial structures and applied GMM to evaluate the developing zebrafish orofacial region. Although similar methods exist for other model organisms and human studies, there have been few applications for zebrafish, and zFACE fills this gap in a way that builds on existing approaches and facilitates cross-species comparisons. We first used zFACE to characterize and understand changes in facial morphology from days 3 to 6 of zebrafish development and established standards at various time points. Additionally, we applied it to analyze smarca4aa8 (yng) mutants and uncovered morphometric alterations that coincided with where smarca4a was expressed, supporting the use of these mutants to inform on mechanisms driving CSS. Collectively, the development and testing of zFACE presented here supports its use as a powerful quantitative tool to uncover previously unappreciated craniofacial alterations in zebrafish genetic models.

To facilitate widespread use for future studies, we added several key features to make zFACE informative, reliable and easy to implement. A detailed protocol is included at https://doi.org/10.6084/m9.figshare.22732166.v1, along with recommendations for ensuring reproducible capture of structures when imaging and placing landmarks. Template files that automatically calculate measurements, perform basic analyses and plot results from user data are available at https://doi.org/10.6084/m9.figshare.22732166.v1. This feature-focused exploration is useful for quickly identifying regions or structures affected by the experimental variable. For example, reduced size can reflect reduced growth during development and point to possible cell population deficiencies driving effects in a particular region, or globally (Juriloff et al., 2004; Young et al., 2007). The ability to perform shape analysis from the same landmarks in MorphoJ (Klingenberg, 2011), a widely used and open-source morphometric program, offers additional advantages: it is reliable, user friendly, and has been implemented in craniofacial studies in both zebrafish and humans.

Previous studies have relied on lateral and ventral views of zebrafish embryos and larvae, preventing acquisition of information on important aspects of facial morphology. The staining and mounting technique presented here is able to facilitate acquisition of rostral view images. Using this view, standard anatomical locations and nomenclature facilitate comparison with human, mouse and other model organism studies (Young et al., 2007; Weinberg et al., 2008; Dickinson, 2016). However, zFACE is also amendable for expansion to analyze the lateral and ventral views with the selection of custom landmarks that capture relevant structures in these views (Fig. S6). Some challenges for implementing custom landmarks include the variability of new anatomical structures as landmarks, as well as the consideration for mounting angle variation and developmental differences that might be present for the structures of interest. Embryonic/larval zebrafish offer a unique opportunity for detection of small effects due to high-quality imaging capabilities and large numbers of individuals that can be assessed for each condition, providing high-throughput, adequately powered studies. Additionally, we selected a moderate number of landmarks to allow for exploratory analyses such as CVA, for which the number of individuals (n) in each group needs to be higher than the number of landmarks (Hallgrimsson et al., 2015). We also designated a sufficient proportion of the landmarks around the mouth to capture the complex and dynamic shape this structure has during the course of development. Lastly, it is important to point out that more advanced analyses might be warranted for some datasets, and the resulting zFACE data can be examined in software packages such as Geomorph or Morpho in R for more complex statistics (Schlager, 2017). Although zFACE handles multiple testing using the Procrustes-based approaches, Euclidean distance matrix analysis, which represents each specimen as a matrix of linear distances between all possible pairs of landmarks and has the ability to detect variation in specific structures or around specific landmarks, can also be utilized (Hallgrimsson et al., 2015).

Here, we applied zFACE to test whether we could detect phenotypic changes after genetic manipulation. Using a mutant with previously reported craniofacial abnormalities but analyzed using only lateral and ventral views, our analysis identified several altered facial features in smarca4a mutant larvae that contribute to a different overall average face shape compared to that of WT and heterozygous animals. The abnormal facial features, together with the alterations in brain morphology, show parallels to those observed in patients with CSS. CSS is a rare congenital disorder that presents with distinct facial features due to craniofacial abnormalities, fifth digit anomalies, microcephaly and intellectual disability (Kosho et al., 2014; Mardinian et al., 2021). The clinical features are heterogeneous, and variants in several genes encoding proteins in the SWI/SNF complex, including SMARCA4A, have been identified as causal, with variants spanning the gene and disrupting several protein domains [Kosho et al., 2014; Online Mendelian Inheritance in Man (OMIM) #614609]. To date, there are no established animal models for the subtype of CSS with SMARCA4 mutations (Celen et al., 2017; Filatova et al., 2019). Thus, the smarca4a zebrafish larvae provide an opportunity for future studies to better understand disruptions in development and identify potential strategies for therapeutic intervention.

The zebrafish smarca4aa8 mutant, has a point mutation (C-to-A transversion) leading to a premature stop codon early in the gene, preceding all functional domains (Gregg et al., 2003). This mutant has been described as having craniofacial abnormalities, alterations in brain size and patterning and other anomalies in neural crest cell-derived tissues (Eroglu et al., 2006). The current study provides high-resolution phenotypic information to further support these abnormal craniofacial characteristics of smarca4a mutants compared to WT and heterozygous controls. Our results identified important similarities in craniofacial features of the mutant with abnormalities reported for CSS. Our results also found severely abnormal oral cavities, suggesting compromised oropharynx function and ability to feed, which is possibly another reason why these mutants only survive to 7 dpf (Link et al., 2000). The frontal view additionally allowed brain measurements of the telencephalon, which showed an altered length-to-width axis, supporting that mutant smarca4a leads to abnormal brain morphology (Eroglu et al., 2006). Future studies can use a similar approach to integrate phenotypic data with genetic and molecular information for greater mechanistic insights on gene functions during development.

The complexity of craniofacial development and the structures that make up the craniofacial regions make the detection of subtle phenotypic variation challenging, requiring advanced analytical tools. zFACE can be implemented in a way that is sufficiently powered to detect subtle morphological changes and better allow for genotype–phenotype correlations, which can serve an important role in the investigation of multifactorial disorders, gene–environment interactions, pharmacological interventions and other studies. Direct comparisons on how morphology changes based on genetic perturbations can also bridge the gap between results from genome-wide association studies and the biological effects of genetic variation on craniofacial structures. In conclusion, zFACE can facilitate important studies to examine and integrate the morphological effects of genetic, environmental or developmental perturbations in zebrafish studies of craniofacial development.

Experimental model and subject details

Experiments were conducted on larval zebrafish (Danio rerio) maintained under standard laboratory conditions with a cycle of 14 h of light and 10 h of darkness. Larvae were collected and kept in E3 larva medium at 28.5°C and staged as previously described (Westerfield, 1993). The zebrafish used in this study were handled in accordance with the guidelines of the University of Texas MD Anderson Cancer Center Institutional Animal Care and Use Committee and UTHealth Animal Welfare Committee. smarca4aa8 (yng) mutant zebrafish (Gregg et al., 2003) were obtained from the Zebrafish International Resource Center.

Specimen collection, staining, mounting and image acquisition

Specimens were collected and fixed overnight on a shaker in 4% fixation solution prepared from dilution of a 36.5% formaldehyde stock solution (formalin; Sigma-Aldrich) into PBS with a small amount of detergent (0.05% Triton-X 100) (0.05% PBST).

Fixation solution was then removed, the specimens were briefly washed with PBST, and DAPI (Thermo Fisher Scientific) was added to PBST at a 1:1000 dilution. The samples were incubated in the DAPI solution for 30 min and rinsed with 0.05% PBST in a series of three 20 min washes. Samples were then rinsed and stored in PBS at 4°C. Each embryo was mounted in a 35 mm glass-bottom culture dish with a 10 mm microwell (MatTek Corporation) filled with 1% low-melt agarose solution. For the frontal/rostral orientation, each embryo was manipulated by moving the tail to suspend the sample upside down, ensuring that the eyes were on the same plane to reduce variability in mounting angles.

To ensure proper orientation of each sample, the clear visualization of the midline between the brain ventricles and the ability to see folds in lower jaw tissues was used as a guide. If the size of the embryo head was too big, the stage was rotated 45° during image acquisition so both eye lenses could be captured, and the last z-stack ensured inclusion of these structures. Typically, the 10× (magnification) objective was used to capture 40 slices of 3.08 μm thickness at 2% laser power using 600 V, digital offset of 2 and digital gain of 1. The images were 1024×1024 at 8 bit/pixel, and the maximum scan speed and averaging of 2 was used. The same magnification was used for all images within a given dataset, and landmarking was performed by the same user.

Coordinate point system

A 26-point system was utilized for landmarks in the craniofacial region of zebrafish larvae. Each point was assigned a number 1 through 26. Confocal images were opened in ImageJ, points were selected in numerical order using the Point Picker tool, and xy coordinates for each point were extracted.

Facial dimension measurements

Using these xy coordinates and the distance formula
(1)
various distance, angular and area measurements were calculated between landmarks. Angular measurements utilized the formula
(2)
where a is the angle of the BC connecting point. Triangular areas were found using the formula
(3)
All calculations were performed using Microsoft Excel.

Statistical analysis of zFACE measurements

A minimum of 28 images were analyzed per stage/condition to provide adequate power for statistical comparison. For statistical analysis of each zFACE measurement, the GraphPad Prism analysis function was used to run one-way ANOVA with Tukey's post hoc tests for multiple comparisons. Bonferroni correction for 39 tests was applied, and P<0.00128 was considered significant.

Multivariate analysis of the zFACE measurements was performed in Stata (StataCorp) and GraphPad Prism. Data were standardized and PCA was performed. Multiple models were run to thoroughly identify/evaluate the most robust model in Stata, in which data were rotated and PC loadings could be calculated. Once modeling was determined in Stata, a rapid and streamlined analysis was rerun in GraphPad Prism, producing similar results.

Shape analysis of zFACE landmark coordinates

Landmark coordinates were imported into MorphoJ, and standard protocols were utilized to perform Procrustes superimposition on principal axes, as well as PCA, CVA and DFA to compare shape changes across the dataset and between groups (Klingenberg, 2011).

Genotyping of smarca4a mutants

smarca4aa8 (yng) mutant zebrafish (Gregg et al., 2003) were identified by a PCR reaction using the following primer sequences: Forward 5′-CCTGTCATGCCCCCTCAGAC-3′; and Reverse 5′-CCGACCCCCACTTTGAGAAC-3′. The resulting 190 bp band was excised, and a restriction digest using RsaI was performed at 37°C for 4 h. RsaI only cuts the WT band, resulting in the WT larvae having two bands (50+140 bp), heterozygous larvae having three bands (50+140+190 bp) and homozygous mutants having one 190 bp band.

Scanning electron microscopy was performed by the CCSG-funded High Resolution Electron Microscopy Facility (National Cancer Institute CA0166772). We thank Jordan Pietz at the Creative Communications group at MD Anderson Cancer Center for the biomedical visualizations.

Author contributions

Conceptualization: L.M., O.E.R., S.T.L., G.T.E.; Methodology: L.M., O.E.R., P.H.K., S.T.L., G.T.E.; Software: L.M., S.S.H.; Validation: L.M., O.E.R., P.H.K., S.T.L., S.S.H.; Formal analysis: L.M., O.E.R., P.H.K., F.C., S.T.L., S.S.H.; Investigation: L.M., O.E.R., P.H.K., F.C., S.T.L., S.S.H., G.T.E.; Resources: L.M., P.H.K., S.S.H., J.T.H.; Data curation: L.M., G.T.E.; Writing - original draft: L.M., G.T.E.; Writing - review & editing: L.M., O.E.R., P.H.K., S.T.L., S.S.H., J.T.H., G.T.E.; Visualization: L.M., O.E.R., P.H.K., S.T.L., G.T.E.; Supervision: G.T.E.; Project administration: G.T.E.; Funding acquisition: G.T.E.

Funding

This work was supported by the Gulf Coast Consortium (John S. Dunn Collaborative Research Award to G.T.E. and J.T.H.); the National Institute of Dental and Craniofacial Research (R01DE011931 to J.T.H.); the Cancer Prevention and Research Institute of Texas (RR14007 to G.T.E.); the National Institute of General Medical Sciences (R01GM124043 to G.T.E.); and the Linda and Mark Quick Award for Basic Science (to G.T.E.). Open Access funding provided by the University of Texas MD Anderson Cancer Center. Deposited in PMC for immediate release.

Data availability

Protocol and template files are available at https://doi.org/10.6084/m9.figshare.22732166.v1. All other relevant data can be found within the article and its supplementary information.

Amsterdam
,
A.
,
Nissen
,
R. M.
,
Sun
,
Z.
,
Swindell
,
E. C.
,
Farrington
,
S.
and
Hopkins
,
N.
(
2004
).
Identification of 315 genes essential for early zebrafish development
.
Proc. Natl. Acad. Sci. USA
101
,
12792
-
12797
.
Bhatia
,
S.
,
Gordon
,
C. T.
,
Foster
,
R. G.
,
Melin
,
L.
,
Abadie
,
V.
,
Baujat
,
G.
,
Vazquez
,
M.-P.
,
Amiel
,
J.
,
Lyonnet
,
S.
,
van Heyningen
,
V.
et al. 
(
2015
).
Functional assessment of disease-associated regulatory variants in vivo using a versatile dual colour transgenesis strategy in zebrafish
.
PLoS Genet.
11
,
e1005193
.
Celen
,
C.
,
Chuang
,
J. C.
,
Luo
,
X.
,
Nijem
,
N.
,
Walker
,
A. K.
,
Chen
,
F.
,
Zhang
,
S.
,
Chung
,
A. S.
,
Nguyen
,
L. H.
,
Nassour
,
I.
et al. 
(
2017
).
Arid1b haploinsufficient mice reveal neuropsychiatric phenotypes and reversible causes of growth impairment
.
eLife
6
,
e25730
.
Chai
,
Y.
and
Maxson
,
R. E.
, Jr
. (
2006
).
Recent advances in craniofacial morphogenesis
.
Dev. Dyn.
235
:
2353
-
2375
.
Choe
,
C. P.
,
Choi
,
S.-Y.
,
Kee
,
Y.
,
Kim
,
M. J.
,
Kim
,
S.-H.
,
Lee
,
Y.
,
Park
,
H.-C.
and
Ro
,
H.
(
2021
).
Transgenic fluorescent zebrafish lines that have revolutionized biomedical research
.
Lab. Anim. Res.
37
,
26
.
Cocos
,
A.
and
Halazonetis
,
D. J.
(
2017
).
Craniofacial shape differs in patients with tooth agenesis: geometric morphometric analysis
.
Eur. J. Orthod.
39
,
345
-
351
.
Coffin
,
G. S.
and
Siris
,
E.
(
1985
).
The Coffin-Siris syndrome
.
Am. J. Dis. Child.
139
,
12
.
Compagnucci
,
C.
,
Martinus
,
K.
,
Griffin
,
J.
and
Depew
,
M. J.
(
2021
).
Programmed cell death not as sledgehammer but as chisel: apoptosis in normal and abnormal craniofacial patterning and development
.
Front. Cell Dev. Biol.
9
,
717404
.
Deml
,
B.
,
Reis
,
L. M.
,
Muheisen
,
S.
,
Bick
,
D.
and
Semina
,
E. V.
(
2015
).
EFTUD2 deficiency in vertebrates: Identification of a novel human mutation and generation of a zebrafish model
.
Birth Defects Res. A Clin. Mol. Teratol
103
,
630
-
640
.
Dickinson
,
A. J.
(
2016
).
Using frogs faces to dissect the mechanisms underlying human orofacial defects
.
Semin. Cell Dev. Biol.
51
,
54
-
63
.
Dougherty
,
M.
,
Kamel
,
G.
,
Grimaldi
,
M.
,
Gfrerer
,
L.
,
Shubinets
,
V.
,
Ethier
,
R.
,
Hickey
,
G.
,
Cornell
,
R. A.
and
Liao
,
E. C.
(
2013
).
Distinct requirements for wnt9a and irf6 in extension and integration mechanisms during zebrafish palate morphogenesis
.
Development
140
,
76
-
81
.
Eames
,
B. F.
,
DeLaurier
,
A.
,
Ullmann
,
B.
,
Huycke
,
T. R.
,
Nichols
,
J. T.
,
Dowd
,
J.
,
McFadden
,
M.
,
Sasaki
,
M. M.
and
Kimmel
,
C. B.
(
2013
).
FishFace: interactive atlas of zebrafish craniofacial development at cellular resolution
.
BMC Dev. Biol.
13
,
23
.
Eberhart
,
J. K.
,
He
,
X.
,
Swartz
,
M. E.
,
Yan
,
Y.-L.
,
Song
,
H.
,
Boling
,
T. C.
,
Kunerth
,
A. K.
,
Walker
,
M. B.
,
Kimmel
,
C. B.
and
Postlethwait
,
J. H.
(
2008
).
MicroRNA Mirn140 modulates Pdgf signaling during palatogenesis
.
Nat. Genet.
40
,
290
-
298
.
Eisenhoffer
,
G. T.
,
Slattum
,
G.
,
Ruiz
,
O. E.
,
Otsuna
,
H.
,
Bryan
,
C. D.
,
Lopez
,
J.
,
Wagner
,
D. S.
,
Bonkowsky
,
J. L.
,
Chien
,
C.-B.
,
Dorsky
,
R. I.
et al. 
(
2017
).
A toolbox to study epidermal cell types in zebrafish
.
J. Cell Sci.
130
,
269
-
277
.
Eroglu
,
B.
,
Wang
,
G.
,
Tu
,
N.
,
Sun
,
X.
and
Mivechi
,
N. F.
(
2006
).
Critical role of Brg1 member of the SWI/SNF chromatin remodeling complex during neurogenesis and neural crest induction in zebrafish
.
Dev. Dyn.
235
,
2722
-
2735
.
Fernandes
,
Y.
and
Lovely
,
C. B.
(
2021
).
Zebrafish models of fetal alcohol spectrum disorders
.
Genesis
59
,
e23460
.
Filatova
,
A.
,
Rey
,
L. K.
,
Lechler
,
M. B.
,
Schaper
,
J.
,
Hempel
,
M.
,
Posmyk
,
R.
,
Szczaluba
,
K.
,
Santen
,
G. W. E.
,
Wieczorek
,
D.
and
Nuber
,
U. A.
(
2019
).
Mutations in SMARCB1 and in other Coffin-Siris syndrome genes lead to various brain midline defects
.
Nat. Commun.
10
,
2966
.
Fruciano
,
C.
(
2016
).
Measurement error in geometric morphometrics
.
Dev. Genes Evol.
226
,
139
-
158
.
Goodwin
,
A. F.
,
Larson
,
J. R.
,
Jones
,
K. B.
,
Liberton
,
D. K.
,
Landan
,
M.
,
Wang
,
Z.
,
Boekelheide
,
A.
,
Langham
,
M.
,
Mushegyan
,
V.
,
Oberoi
,
S.
et al. 
(
2014
).
Craniofacial morphometric analysis of individuals with X-linked hypohidrotic ectodermal dysplasia
.
Mol. Genet. Genomic Med.
2
,
422
-
429
.
Gregg
,
R. G.
,
Willer
,
G. B.
,
Fadool
,
J. M.
,
Dowling
,
J. E.
and
Link
,
B. A.
(
2003
).
Positional cloning of the young mutation identifies an essential role for the Brahma chromatin remodeling complex in mediating retinal cell differentiation
.
Proc. Natl. Acad. Sci. USA
100
,
6535
-
6540
.
Hallgrimsson
,
B.
,
Percival
,
C. J.
,
Green
,
R.
,
Young
,
N. M.
,
Mio
,
W.
and
Marcucio
,
R.
(
2015
).
Morphometrics, 3D imaging, and craniofacial development
.
Curr. Top. Dev. Biol.
115
,
561
-
597
.
Hennekam
,
R. C. M.
,
Allanson
,
J. E.
,
Krantz
,
I.D.
and
Gorlin
,
R. J.
(
2010
).
Gorlin's syndromes of the head and neck
,
5th edn.
Oxford University Press
.
Ji
,
Y.
,
Garland
,
M. A.
,
Sun
,
B.
,
Zhang
,
S.
,
Reynolds
,
K.
,
McMahon
,
M.
,
Rajakumar
,
R.
,
Islam
,
M. S.
,
Liu
,
Y.
,
Chen
,
Y.
et al. 
(
2020
).
Cellular and developmental basis of orofacial clefts
.
Birth Defects Res.
112
:
1558
-
1587
.
Juriloff
,
D. M.
,
Harris
,
M. J.
and
Dewell
,
S. L.
(
2004
).
A digenic cause of cleft lip in A-strain mice and definition of candidate genes for the two loci
.
Birth Defects Res. A Clin. Mol. Teratol
70
,
509
-
518
.
Katsadouris
,
A.
and
Halazonetis
,
D. J.
(
2017
).
Geometric morphometric analysis of craniofacial growth between the ages of 12 and 14 in normal humans
.
Eur. J. Orthod.
39
,
386
-
394
.
Kesterke
,
M. J.
,
Judd
,
M. A.
,
Mooney
,
M. P.
,
Siegel
,
M. I.
,
Elsalanty
,
M.
,
Howie
,
R. N.
,
Weinberg
,
S. M.
and
Cray
,
J. J.
(
2018
).
Maternal environment and craniofacial growth: geometric morphometric analysis of mandibular shape changes with in utero thyroxine overexposure in mice
.
J. Anat.
233
,
46
-
54
.
Klingenberg
,
C. P.
(
2011
).
MorphoJ: an integrated software package for geometric morphometrics
.
Mol. Ecol. Resour.
11
,
353
-
357
.
Klingenberg
,
C. P.
and
McIntyre
,
G. S.
(
1998
).
Geometric morphometrics of developmental instability: analyzing patterns of fluctuating asymmetry with procrustes methods
.
Evolution
52
,
1363
-
1375
.
Klingenberg
,
C. P.
,
Barluenga
,
M.
and
Meyer
,
A.
(
2002
).
Shape analysis of symmetric structures: quantifying variation among individuals and asymmetry
.
Evolution
56
,
1909
-
1920
.
Kosho
,
T.
,
Okamoto
,
N.
and
Coffin-Siris Syndrome International Collaborators
. (
2014
).
Genotype-phenotype correlation of Coffin-Siris syndrome caused by mutations in SMARCB1, SMARCA4, SMARCE1, and ARID1A
.
Am. J. Med. Genet. C Semin. Med. Genet.
166C
,
262
-
275
.
Küry
,
S.
,
Besnard
,
T.
,
Ebstein
,
F.
,
Khan
,
T. N.
,
Gambin
,
T.
,
Douglas
,
J.
,
Bacino
,
C. A.
,
Craigen
,
W. J.
,
Sanders
,
S. J.
,
Lehmann
,
A.
et al. 
(
2017
).
De novo disruption of the proteasome regulatory subunit PSMD12 causes a syndromic neurodevelopmental disorder
.
Am. J. Hum. Genet.
100
,
352
-
363
.
Liang
,
S.
,
Wu
,
J.
,
Weinberg
,
S. M.
and
Shapiro
,
L. G.
(
2013
).
Improved detection of landmarks on 3D human face data
.
Annu. Int. Conf. IEEE Eng. Med. Biol. Soc.
2013
,
6482
-
6485
.
Link
,
B. A.
,
Fadool
,
J. M.
,
Malicki
,
J.
and
Dowling
,
J. E.
(
2000
).
The zebrafish young mutation acts non-cell-autonomously to uncouple differentiation from specification for all retinal cells
.
Development
127
,
2177
-
2188
.
Machado
,
R. G.
and
Eames
,
B. F.
(
2017
).
Using zebrafish to test the genetic basis of human craniofacial diseases
.
J. Dent. Res.
96
,
1192
-
1199
.
Mardinian
,
K.
,
Adashek
,
J. J.
,
Botta
,
G. P.
,
Kato
,
S.
and
Kurzrock
,
R.
(
2021
).
SMARCA4: implications of an altered chromatin-remodeling gene for cancer development and therapy
.
Mol. Cancer Ther.
20
,
2341
-
2351
.
McCarthy
,
N.
,
Wetherill
,
L.
,
Lovely
,
C. B.
,
Swartz
,
M. E.
,
Foroud
,
T. M.
and
Eberhart
,
J. K.
(
2013
).
Pdgfra protects against ethanol-induced craniofacial defects in a zebrafish model of FASD
.
Development
140
,
3254
-
3265
.
Mercan
,
E.
,
Hopper
,
R. A.
and
Maga
,
A. M.
(
2020
).
Cranial growth in isolated sagittal craniosynostosis compared with normal growth in the first 6 months of age
.
J. Anat.
236
,
105
-
116
.
Mork
,
L.
and
Crump
,
G.
(
2015
).
Zebrafish craniofacial development: a window into early patterning
.
Curr. Top. Dev. Biol.
115
,
235
-
269
.
Nichols
,
J. T.
,
Pan
,
L.
,
Moens
,
C. B.
and
Kimmel
,
C. B.
(
2013
).
barx1 represses joints and promotes cartilage in the craniofacial skeleton
.
Development
140
,
2765
-
2775
.
Rutland
,
J. W.
,
Bellaire
,
C. P.
,
Yao
,
A.
,
Arrighi-Allisan
,
A.
,
Napoli
,
J. G.
,
Delman
,
B. N.
and
Taub
,
P. J.
(
2021
).
The expanding role of geometric morphometrics in craniofacial surgery
.
J. Craniofac. Surg.
32
,
1104
-
1109
.
Schlager
,
S.
(
2017
).
Morpho and Rvcg – Shape Analysis in R
. In
Statistical Shape and Deformation Analysis
(ed.
G.
Zheng
,
S.
Li
and
G.
Szekely
). pp.
217
-
256
.
Academic Press
.
Shaw
,
N. D.
,
Brand
,
H.
,
Kupchinsky
,
Z. A.
,
Bengani
,
H.
,
Plummer
,
L.
,
Jones
,
T. I.
,
Erdin
,
S.
,
Williamson
,
K. A.
,
Rainger
,
J.
,
Stortchevoi
,
A.
et al. 
(
2017
).
SMCHD1 mutations associated with a rare muscular dystrophy can also cause isolated arhinia and Bosma arhinia microphthalmia syndrome
.
Nat. Genet.
49
,
238
-
248
.
Shull
,
L. C.
,
Lencer
,
E. S.
,
Kim
,
H. M.
,
Goyama
,
S.
,
Kurokawa
,
M.
,
Costello
,
J. C.
,
Jones
,
K.
and
Artinger
,
K. B.
(
2022
).
PRDM paralogs antagonistically balance Wnt/β-catenin activity during craniofacial chondrocyte differentiation
.
Development
149
,
dev200082
.
Špelić
,
I.
,
Rezić
,
A.
,
Kanjuh
,
T.
,
Marić
,
A.
,
Maguire
,
I.
,
Simonović
,
P.
,
Radočaj
,
T.
and
Piria
,
M.
(
2021
).
Application of the geometric morphometrics approach in the discrimination of morphological traits between brown trout lineages in the Danube Basin of Croatia
.
Knowledge Management Aquat. Ecosyst.
22
,
422
.
Suurväli
,
J.
,
Whiteley
,
A. R.
,
Zheng
,
Y.
,
Gharbi
,
K.
,
Leptin
,
M.
and
Wiehe
,
T.
(
2020
).
The laboratory domestication of zebrafish: from diverse populations to inbred substrains
.
Mol. Biol. Evol.
37
,
1056
-
1069
.
Swartz
,
M. E.
,
Sheehan-Rooney
,
K.
,
Dixon
,
M. J.
and
Eberhart
,
J. K.
(
2011
).
Examination of a palatogenic gene program in zebrafish
.
Dev. Dyn.
240
,
2204
-
2220
.
Swindell
,
E. C.
,
Yuan
,
Q.
,
Maili
,
L. E.
,
Tandon
,
B.
,
Wagner
,
D. S.
and
Hecht
,
J. T.
(
2015
).
Crispld2 is required for neural crest cell migration and cell viability during zebrafish craniofacial development
.
Genesis
53
,
660
-
667
.
Teng
,
C. S.
,
Ting
,
M.-C.
,
Farmer
,
D. T.
,
Brockop
,
M.
,
Maxson
,
R. E.
and
Crump
,
J. G.
(
2018
).
Altered bone growth dynamics prefigure craniosynostosis in a zebrafish model of Saethre-Chotzen syndrome
.
eLife
7
,
e37024
.
Truong
,
B. T.
and
Artinger
,
K. B.
(
2021
).
The power of zebrafish models for understanding the co-occurrence of craniofacial and limb disorders
.
Genesis
59
,
e23407
.
Tsurusaki
,
Y.
,
Okamoto
,
N.
,
Ohashi
,
H.
,
Kosho
,
T.
,
Imai
,
Y.
,
Hibi-Ko
,
Y.
,
Kaname
,
T.
,
Naritomi
,
K.
,
Kawame
,
H.
,
Wakui
,
K.
et al. 
(
2012
).
Mutations affecting components of the SWI/SNF complex cause Coffin-Siris syndrome
.
Nat. Genet.
44
,
376
-
378
.
Twigg
,
S. R. F.
and
Wilkie
,
A. O. M.
(
2015
).
New insights into craniofacial malformations
.
Hum. Mol. Genet.
24
,
R50
-
R59
.
Van Otterloo
,
E.
,
Williams
,
T.
and
Artinger
,
K. B.
(
2016
).
The old and new face of craniofacial research: How animal models inform human craniofacial genetic and clinical data
.
Dev. Biol.
415
,
171
-
187
.
Vrdoljak
,
J.
,
Imanol Sanchez
,
K.
,
Arreola-Ramos
,
R.
,
Diaz Huesa
,
E. G.
,
Villagra
,
A.
,
Javier Avila
,
L.
and
Morando
,
M.
(
2021
).
Testing repeatability, measurement error and species differentiation when using geometric morphometrics on complex shapes: a case study of Patagonian lizards of the genus Liolaemus (Squamata: Liolaemini)
.
Biol. J. Linn. Soc.
130
,
800
-
812
.
Weinberg
,
S. M.
,
Neiswanger
,
K.
,
Richtsmeier
,
J. T.
,
Maher
,
B. S.
,
Mooney
,
M. P.
,
Siegel
,
M. I.
and
Marazita
,
M. L.
(
2008
).
Three-dimensional morphometric analysis of craniofacial shape in the unaffected relatives of individuals with nonsyndromic orofacial clefts: a possible marker for genetic susceptibility
.
Am. J. Med. Genet. A
146A
:
409
-
420
.
Weinberg
,
S. M.
,
Naidoo
,
S. D.
,
Bardi
,
K. M.
,
Brandon
,
C. A.
,
Neiswanger
,
K.
,
Resick
,
J. M.
,
Martin
,
R. A.
and
Marazita
,
M. L.
(
2009
).
Face shape of unaffected parents with cleft affected offspring: combining three-dimensional surface imaging and geometric morphometrics
.
Orthod. Craniofac. Res.
12
,
271
-
281
.
Westerfield
,
M.
(
1993
).
The Zebrafish Book: A Guide for the Laboratory use of Zebrafish (Brachydanio rerio)
.
Eugene
:
M. Westerfield
.
Young
,
N. M.
,
Wat
,
S.
,
Diewert
,
V. M.
,
Browder
,
L. W.
and
Hallgrímsson
,
B.
(
2007
).
Comparative morphometrics of embryonic facial morphogenesis: implications for cleft-lip etiology
.
Anat. Rec.
290
,
123
-
139
.
Yuan
,
Q.
,
Chiquet
,
B. T.
,
Devault
,
L.
,
Warman
,
M. L.
,
Nakamura
,
Y.
,
Swindell
,
E. C.
and
Hecht
,
J. T.
(
2012
).
Craniofacial abnormalities result from knock down of nonsyndromic clefting gene, crispld2, in zebrafish
.
Genesis
50
,
871
-
881
.

Competing interests

The authors declare no competing or financial interests.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

Supplementary information