Exploring Subclinical Phenotypic Features in Twin Pairs Discordant for Cleft Lip and Palate

Symmetry and fluctuation of cell movements in neural crest-derived facial mesenchyme

In the face, symmetry is established when bilateral streams of neural crest cells leave the neural tube at the same time, follow identical migration routes and then give rise to the facial prominences. However, developmental instability exists, particularly surrounding the steps of lip fusion. The causes of instability are unknown but inability to cope with developmental fluctuations are a likely cause of congenital malformations, such as non-syndromic orofacial clefts. Here, we tracked cell movements over time in the frontonasal mass, which forms the facial midline and participates in lip fusion, using live-cell imaging of chick embryos. Our mathematical examination of cell velocity vectors uncovered temporal fluctuations in several parameters, including order/disorder, symmetry/asymmetry and divergence/convergence. We found that treatment with a Rho GTPase inhibitor completely disrupted the temporal fluctuations in all measures and blocked morphogenesis. Thus, we discovered that genetic control of symmetry extends to mesenchymal cell movements and that these movements are of the type that could be perturbed in asymmetrical malformations, such as non-syndromic cleft lip.

This article has an associated ‘The people behind the papers’ interview.

Highlighted Article: Live imaging of the chick embryo face followed by mathematical analysis of mesenchymal cell tracks captures novel fluctuations between states of order/disorder as well as symmetry/asymmetry, revealing developmental instabilities that are part of normal morphogenesis.

Six2 regulates Pax9 expression, palatogenesis and craniofacial bone formation

In this study, we investigated the role of the transcription factor Six2 in palate development. Six2 was selected using the SysFACE tool to predict genes from the 2p21 locus, a region associated with clefting in humans by GWAS, that are likely to be involved in palatogenesis. We functionally validated the predicted role of Six2 in palatogenesis by showing that 22% of Six2 null embryos develop cleft palate. Six2 contributes to palatogenesis by promoting mesenchymal cell proliferation and regulating bone formation. The clefting phenotype in Six2^-/- embryos is similar to Pax9 null embryos, so we examined the functional relationship of these two genes. Mechanistically, SIX2 binds to a PAX9 5' upstream regulatory element and activates PAX9 expression. In addition, we identified a human SIX2 coding variant (p.Gly264Glu) in a proband with cleft palate. We show this missense mutation affects the stability of the SIX2 protein and leads to decreased PAX9 expression. The low penetrance of clefting in the Six2 null mouse combined with the mutation in one patient with cleft palate underscores the potential combinatorial interactions of other genes in clefting. Our study demonstrates that Six2 interacts with the developmental gene regulatory network in the developing palate.

Nonlinear gene expression-phenotype relationships contribute to variation and clefting in the A/WySn mouse

BACKGROUND

Cleft lip and palate is one of the most common human birth defects, but the underlying etiology is poorly understood. The A/WySn mouse is a spontaneously occurring model of multigenic clefting in which 20% to 30% of individuals develop an orofacial cleft. Recent work has shown altered methylation at a specific retrotransposon insertion downstream of the Wnt9b locus in clefting animals, which results in decreased Wnt9b expression.

RESULTS

Using a newly developed protocol that allows us to measure morphology, gene expression, and DNA methylation in the same embryo, we relate gene expression in an individual embryo directly to its three-dimensional morphology for the first time. We find that methylation at the retrotransposon relates to Wnt9b expression and morphology. IAP methylation relates to shape of the nasal process in a manner consistent with clefting. Embryos with low IAP methylation exhibit increased among-individual variance in facial shape.

CONCLUSIONS

Methylation and gene expression relate nonlinearly to nasal process morphology. Individuals at one end of a continuum of phenotypic states display a clinical phenotype and increased phenotypic variation. Variable penetrance and expressivity in this model is likely determined both by among-individual variation in methylation and changes in phenotypic robustness along the underlying liability distribution for orofacial clefting.

A systematic review and meta-analysis on protective role of forkhead box E1 (FOXE1) polymorphisms in susceptibility to non-syndromic cleft lip/palate

INTRODUCTION

Several environmental and genetic factors have a role in the aetiology of non-syndromic cleft lip/palate (NSCL/P). This meta-analysis evaluated the association of rs3758249 and rs4460498 forkhead box E1 (FOXE1) polymorphisms with the NSCL/P risk.

MATERIALS AND METHODS

The Scopus, Cochrane Library, Web of Science, and PubMed databases were searched for articles published until March 2019. The analyses were performed by Review Manager 5.3 using the crude odds ratio (OR) and 95% confidence interval (CI) for a strong association between FOXE1 polymorphisms and the risk of NSCL/P.

RESULTS

Out of 161 articles retrieved from the databases, four case-control articles were involved in the meta-analysis. The pooled ORs of rs4460498 polymorphism based on allelic, homozygous, heterozygous, dominant, and recessive models were 0.74 (95% CI: 0.69, 0.80; P<0.00001), 0.43 (95% CI: 0.30, 0.61; P<0.00001), 0.66 (95% CI: 0.55, 0.80; P<0.0001), 0.66 (95% CI: 0.59, 0.73; P<0.00001), and 0.70 (95% CI: 0.60, 0.82; P<0.0001), respectively; whereas, the pooled OR of rs3758249 polymorphism were 0.86 (95% CI: 0.71, 1.04; P=0.12), 0.68 (95% CI: 0.57, 0.82; P<0.0001), 0.79 (95% CI: 0.57, 1.09; P=0.15), 0.79 (95% CI: 0.58, 1.08; P=0.14), and 0.80 (95% CI: 0.68, 0.95; P=0.010) for the afore-mentioned models, respectively.

CONCLUSIONS

The results showed that the T allele, TT, and CT genotypes of rs4460498 polymorphism were significantly associated with a decreased risk of NSCL/P; whereas, for rs3758249 polymorphism, only the AA genotype had a significant protective role in NSCL/P. Thus, FOXE1 is strongly associated with NSCL/P in the populations.

Association of low-frequency genetic variants in regulatory regions with nonsyndromic orofacial clefts

Genome-wide scans have shown that common risk alleles for orofacial clefts (OFC) tend to be located in noncoding regulatory elements and cumulatively explain only part of the heritability of OFCs. Low-frequency variants may account for some of the "missing" heritability. Therefore, we scanned low-frequency variants located within putative craniofacial enhancers to identify novel OFC risk variants and implicate new regulatory elements in OFC pathogenesis. Analyses were performed in a multiethnic sample of 1,995 cases of cleft lip with or without cleft palate (CL/P), 221 cases with cleft palate (CP) only, and 1,576 unaffected controls. One hundred and nineteen putative craniofacial enhancers identified from ChIP-Seq studies in craniofacial tissues or cell lines contained multiple low-frequency (0.01-1%) variants, which we genotyped in participants using a custom Illumina panel. Two complementary statistical approaches, sequence kernel association test and combined multivariate and collapsing, were used to test association of the aggregated low-frequency variants across each enhancer region with CL/P and CP. We discovered a significant association between CP and a branchial arch enhancer near FOXP1 (mm60; p-value = .0002). Additionally, we observed a suggestive association between CL/P and a forebrain enhancer near FOXE1 (hs1717; p-value = .001). These findings suggest that low-frequency variants in craniofacial enhancer regions contribute to the complex etiology of nonsyndromic OFCs.

Speech Phenotyping in Unaffected Family Members of Individuals With Nonsyndromic Cleft Lip With or Without Palate

OBJECTIVE

Subclinical phenotypes of nonsyndromic cleft lip with or without cleft palate (CL ± P) may be identified from clinically "unaffected" relatives and could be associated with specific cleft-related gene mutations. It has been hypothesized that velopharyngeal insufficiency (VPI) may be a subclinical phenotype of interest in this population, but this has not been explored quantitatively with appropriate control cohorts. The aim of this case-control study was to compare VPI in at-risk clinically unaffected relatives of individuals with nonsyndromic CL ± P with a low-risk matched normative Australian cohort.

PARTICIPANTS

Clinically unaffected (ie, with no overt cleft) first-degree relatives of a proband with nonsyndromic CL ± P (n = 189) and noncleft controls (n = 207).

MAIN OUTCOME MEASURE(S)

Perceptual measures of VPI encompassing resonance, nasal emission, and articulation were evaluated using the Great Ormond Street Speech Assessment. Quantitative measures of VPI were obtained from the Nasometer II using standardized adult and pediatric speech stimuli.

RESULTS

Both perceptual and instrumental measures showed no significant difference (P > .01) between the VPI in unaffected relatives and the noncleft comparison group. Mean nasalance scores for both groups were calculated and reported according to speech stimuli, age, and sex.

CONCLUSIONS

Results suggest that VPI, measured through speech, is not a significant subclinical phenotype of nonsyndromic CL ± P. Therefore, further familial genetic investigations exploring VPI may not yield meaningful results. Exploration across multiple subclinical phenotypes in larger cohorts may enable researchers to better understand the multifaceted nature of this complex and heterogeneous anomaly.

Ultrastructural analysis of collagen fibril diameter distribution in cleft lip

OBJECTIVE

A preliminary study to determine collagen fibril diameter (CF-ED) distribution on medial and lateral sides of cleft lip (CL).

MATERIAL AND METHODS

Tissue samples from medial and lateral sides of CL were fixed in 2.5% glutaraldehyde and 1% osmium tetroxide and embedded in Araldite CY212 resin for transmission electron microscopy. The analysis of CF-ED was performed using the ImageJ program. To characterize the packaging of collagen fibrils (CFs) in the two tissues, we estimated the collagen number density (CF-ND) and fibril-area-fraction (FAF). Differences in measurements across the two sides were calculated using Wilcoxon signed-rank test.

RESULTS

The CF-ED was statistically significantly (p < 0.001) smaller on the medial side (45.69 ± 7.89 nm) than on the lateral side (54.18 ± 7.62 nm). The medial side had a higher CF-ND and a higher percentage of FAF than the lateral side.

CONCLUSION

Our finding of a smaller CF-ED and higher CF-ND and FAF for the medial side suggests possible differences in size and distribution of CFs between medial and lateral sides of CL. This finding provides knowledge toward underlying tissue biomechanics that may help reconstruction of perioral tissue scaffolds, ultimately resulting in better treatment of patients with oral clefts.

Testing the face shape hypothesis in twins discordant for nonsyndromic orofacial clefting

Nonsyndromic orofacial clefts (OFCs) are complex traits characterized by multifactorial inheritance and wide phenotypic variability. Numerous studies have shown subtle differences in the faces of unaffected relatives from cleft families compared to controls, the implication being that such outward differences are an incomplete expression reflecting an underlying genetic predisposition. Twins discordant for OFCs provide a unique opportunity to further test this idea, as the unaffected co-twin shares on average 50% (for dizygotic twins) and 100% (for monozygotic twins) of the genetic risk factors as the affected twin. We used 3D surface imaging and spatially-dense morphometry to compare facial shape in a sample of 44 unaffected co-twins and age- and sex-matched unaffected controls (n = 241). Unaffected co-twins showed statistically significant differences in the midface, lateral upper face, and forehead regions, compared to controls. Furthermore, co-twins were characterized by a distinct pattern of midfacial retrusion, broader upper faces, and greater protrusion of the mandible and brow ridges. This same general facial pattern was shown in both unaffected monozygotic and dizygotic co-twin subsets. These results provide additional support that altered facial shape is a phenotypic marker for OFC susceptibility.