SPECC1L regulates palate development downstream of IRF6.

SPECC1L mutations have been identified in patients with rare atypical orofacial clefts and with syndromic cleft lip and/or palate (CL/P). These mutations cluster in the second coiled-coil and calponin homology domains of SPECC1L and severely affect the ability of SPECC1L to associate with microtubules. We previously showed that gene-trap knockout of Specc1l in mouse results in early embryonic lethality. We now present a truncation mutant mouse allele, Specc1lΔC510, that results in perinatal lethality. Specc1lΔC510/ΔC510 homozygotes showed abnormal palate rugae but did not show cleft palate. However, when crossed with a gene-trap allele, Specc1lcGT/ΔC510 compound heterozygotes showed a palate elevation delay with incompletely penetrant cleft palate. Specc1lcGT/ΔC510 embryos exhibit transient oral epithelial adhesions at E13.5, which may delay shelf elevation. Consistent with oral adhesions, we show periderm layer abnormalities, including ectopic apical expression of adherens junction markers, similar to Irf6 hypomorphic mutants and Arhgap29 heterozygotes. Indeed, SPECC1L expression is drastically reduced in Irf6 mutant palatal shelves. Finally, we wanted to determine if SPECC1L deficiency also contributed to non-syndromic (ns) CL/P. We sequenced 62 Caucasian, 89 Filipino, 90 Ethiopian, 90 Nigerian and 95 Japanese patients with nsCL/P and identified three rare coding variants (p.Ala86Thr, p.Met91Iso and p.Arg546Gln) in six individuals. These variants reside outside of SPECC1L coiled-coil domains and result in milder functional defects than variants associated with syndromic clefting. Together, our data indicate that palate elevation is sensitive to deficiency of SPECC1L dosage and function and that SPECC1L cytoskeletal protein functions downstream of IRF6 in palatogenesis.

The TFAP2A-IRF6-GRHL3 genetic pathway is conserved in neurulation.

Mutations in IRF6, TFAP2A and GRHL3 cause orofacial clefting syndromes in humans. However, Tfap2a and Grhl3 are also required for neurulation in mice. Here, we found that homeostasis of Irf6 is also required for development of the neural tube and associated structures. Over-expression of Irf6 caused exencephaly, a rostral neural tube defect, through suppression of Tfap2a and Grhl3 expression. Conversely, loss of Irf6 function caused a curly tail and coincided with a reduction of Tfap2a and Grhl3 expression in tail tissues. To test whether Irf6 function in neurulation was conserved, we sequenced samples obtained from human cases of spina bifida and anencephaly. We found two likely disease-causing variants in two samples from patients with spina bifida. Overall, these data suggest that the Tfap2a-Irf6-Grhl3 genetic pathway is shared by two embryologically distinct morphogenetic events that previously were considered independent during mammalian development. In addition, these data suggest new candidates to delineate the genetic architecture of neural tube defects and new therapeutic targets to prevent this common birth defect.

Secondary Palate Development in the Dog (Canis lupus familiaris).

OBJECTIVE: To investigate the gestational timing of morphologic events in normal canine secondary palate development as a baseline for studies in dog models of isolated cleft palate (CP). METHODS: Beagle and beagle/cocker spaniel-hybrid fetal dogs were obtained by cesarean-section on various days of gestation, timed from the initial rise of serum progesterone concentration. Morphology of fetal heads was determined by examining serial coronal sections. RESULTS: On gestational day 35 (d35), the palatal shelves pointed ventrally alongside the tongue. On d36, palatal shelves were elongated and elevated to a horizontal position above the tongue but were not touching. On d37, palatine shelves and vomer were touching, but the medial epithelial seam (MES) between the apposed shelves remained. Immunostaining with epithelial protein markers showed that the MES gradually dissolved and was replaced by mesenchyme during d37-d44, and palate fusion was complete by d44. Examination of remnant MES suggested that fusion of palatal shelves began in mid-palate and moved rostrally and caudally. CONCLUSION: Palate development occurs in dogs in the steps described in humans and mice, but palate closure occurs at an intermediate time in gestation. These normative data will form the basis of future studies to determine pathophysiologic mechanisms in dog models of CP. Added clinical significance is the enhancement of dogs as a large animal model to test new approaches for palate repair, with the obvious advantage of achieving full maturity within 2 years rather than 2 decades.

Endogenous Avian Leukosis Virus in Combination with Serotype 2 Marek's Disease Virus Significantly Boosted the Incidence of Lymphoid Leukosis-Like Bursal Lymphomas in Susceptible Chickens.

In 2010, sporadic cases of avian leukosis virus (ALV)-like bursal lymphoma, also known as spontaneous lymphoid leukosis (LL)-like tumors, were identified in two commercial broiler breeder flocks in the absence of exogenous ALV infection. Two individual ALV subgroup E (ALV-E) field strains, designated AF227 and AF229, were isolated from two different breeder farms. The role of these ALV-E field isolates in development of and the potential joint impact in conjunction with a Marek's disease virus (MDV) vaccine (SB-1) were further characterized in chickens of an experimental line and commercial broiler breeders. The experimental line 0.TVB*S1, commonly known as the rapid feathering-susceptible (RFS) line, of chickens lacks all endogenous ALV and is fully susceptible to all subgroups of ALV, including ALV-E. Spontaneous LL-like tumors occurred following infection with AF227, AF229, and a reference ALV-E strain, RAV60, in RFS chickens. Vaccination with serotype 2 MDV, SB-1, in addition to AF227 or AF229 inoculation, significantly enhanced the spontaneous LL-like tumor incidence in the RFS chickens. The spontaneous LL-like tumor incidence jumped from 14% by AF227 alone to 42 to 43% by AF227 in combination with SB-1 in the RFS chickens under controlled conditions. RNA-sequencing analysis of the LL-like lymphomas and nonmalignant bursa tissues of the RFS line of birds identified hundreds of differentially expressed genes that are reportedly involved in key biological processes and pathways, including signaling and signal transduction pathways. The data from this study suggested that both ALV-E and MDV-2 play an important role in enhancement of the spontaneous LL-like tumors in susceptible chickens. The underlying mechanism may be complex and involved in many chicken genes and pathways, including signal transduction pathways and immune system processes, in addition to reported viral genes.IMPORTANCE Lymphoid leukosis (LL)-like lymphoma is a low-incidence yet costly and poorly understood disease of domestic chickens. The observed unique characteristics of LL-like lymphomas are that the incidence of the disease is chicken line dependent; pathologically, it appeared to mimic avian leukosis but is free of exogenous ALV infection; inoculation of the nonpathogenic ALV-E or MDV-2 (SB-1) boosts the incidence of the disease; and inoculation of both the nonpathogenic ALV-E and SB-1 escalates it to much higher levels. This study was designed to test the impact of two new ALV-E isolates, recently derived from commercial broiler breeder flocks, in combination with the nonpathogenic SB-1 on LL-like lymphoma incidences in both an experimental egg layer line of chickens and a commercial broiler breeder line of chickens under a controlled condition. Data from this study provided an additional piece of experimental evidence on the potency of nonpathogenic ALV-E, MDV-2, and ALV-E plus MDV-2 in boosting the incidence of LL-like lymphomas in susceptible chickens. This study also generated the first piece of genomic evidence that suggests host transcriptomic variation plays an important role in modulating LL-like lymphoma formation.

Rapid functional analysis of computationally complex rare human IRF6 gene variants using a novel zebrafish model.

Large-scale sequencing efforts have captured a rapidly growing catalogue of genetic variations. However, the accurate establishment of gene variant pathogenicity remains a central challenge in translating personal genomics information to clinical decisions. Interferon Regulatory Factor 6 (IRF6) gene variants are significant genetic contributors to orofacial clefts. Although approximately three hundred IRF6 gene variants have been documented, their effects on protein functions remain difficult to interpret. Here, we demonstrate the protein functions of human IRF6 missense gene variants could be rapidly assessed in detail by their abilities to rescue the irf6 -/- phenotype in zebrafish through variant mRNA microinjections at the one-cell stage. The results revealed many missense variants previously predicted by traditional statistical and computational tools to be loss-of-function and pathogenic retained partial or full protein function and rescued the zebrafish irf6 -/- periderm rupture phenotype. Through mRNA dosage titration and analysis of the Exome Aggregation Consortium (ExAC) database, IRF6 missense variants were grouped by their abilities to rescue at various dosages into three functional categories: wild type function, reduced function, and complete loss-of-function. This sensitive and specific biological assay was able to address the nuanced functional significances of IRF6 missense gene variants and overcome many limitations faced by current statistical and computational tools in assigning variant protein function and pathogenicity. Furthermore, it unlocked the possibility for characterizing yet undiscovered human IRF6 missense gene variants from orofacial cleft patients, and illustrated a generalizable functional genomics paradigm in personalized medicine.

Identification of functional variants for cleft lip with or without cleft palate in or near PAX7, FGFR2, and NOG by targeted sequencing of GWAS loci.

Although genome-wide association studies (GWASs) for nonsyndromic orofacial clefts have identified multiple strongly associated regions, the causal variants are unknown. To address this, we selected 13 regions from GWASs and other studies, performed targeted sequencing in 1,409 Asian and European trios, and carried out a series of statistical and functional analyses. Within a cluster of strongly associated common variants near NOG, we found that one, rs227727, disrupts enhancer activity. We furthermore identified significant clusters of non-coding rare variants near NTN1 and NOG and found several rare coding variants likely to affect protein function, including four nonsense variants in ARHGAP29. We confirmed 48 de novo mutations and, based on best biological evidence available, chose two of these for functional assays. One mutation in PAX7 disrupted the DNA binding of the encoded transcription factor in an in vitro assay. The second, a non-coding mutation, disrupted the activity of a neural crest enhancer downstream of FGFR2 both in vitro and in vivo. This targeted sequencing study provides strong functional evidence implicating several specific variants as primary contributory risk alleles for nonsyndromic clefting in humans.

Altered thymocyte and T cell development in neonatal mice with hyperoxia-induced lung injury.

BACKGROUND: The adaptive immune system of neonates is relatively underdeveloped. The thymus is an essential organ for adaptive T cell development and might be affected during the natural course of oxygen induced lung injury. The effect of prolonged hyperoxia on the thymus, thymocyte and T cell development, and its proliferation has not been studied extensively. METHODS: Neonatal mice were exposed to 85% oxygen (hyperoxia) or room air (normoxia) up to 28 days. Flow cytometry using surface markers were used to assay for thymocyte development and proliferation. RESULTS: Mice exposed to prolonged hyperoxia had evidence of lung injury associated alveolar simplification, a significantly lower mean weight, smaller thymic size, lower mean thymocyte count and higher percentage of apoptotic thymocytes. T cells subpopulation in the thymus showed a significant reduction in the count and proliferation of double positive and double negative T cells. There was a significant reduction in the count and proliferation of single positive CD4+ and CD8+ T cells. CONCLUSIONS: Prolonged hyperoxia in neonatal mice adversely affected thymic size, thymocyte count and altered the distribution of T cells sub-populations. These results are consistent with the hypothesis that prolonged hyperoxia causes defective development of T cells in the thymus.

IRF6 expression in basal epithelium partially rescues Irf6 knockout mice.

BACKGROUND: Mutations in IRF6, CHUK (IKKA), and RIPK4 can lead to a disease spectrum that includes cutaneous, limb, and craniofacial malformations. Loss of these alleles in the mouse leads to perinatal lethality and severe cutaneous, limb, and craniofacial defects also. Genetic rescue in the mouse has been shown for Ikka and Ripk4. RESULTS: Here, we show partial genetic rescue of Irf6 knockout embryos using the KRT14 promoter to drive Irf6 expression in the basal epithelium. In contrast to Irf6 knockout embryos, rescue embryos survive the immediate perinatal period. Macroscopic examination reveals rescue of skin adhesions between the axial and appendicular skeleton. Unexpectedly, KRT14-driven Irf6 expression does not completely rescue orofacial clefting and adhesions between the palate and tongue, suggesting the importance of cell-autonomous IRF6 expression in periderm. Like knockout embryos, Irf6 rescue embryos also have persistent esophageal adhesions, which likely contribute to postnatal demise. CONCLUSIONS: Together, these data suggest that targeted expression of IRF6 can significantly reduce disease severity, but that a minimum level of Irf6 in both periderm and basal epithelial cells is necessary for orofacial development. Therefore, homologous human and mouse phenotypes are observed for IRF6, IKKA, and RIPK4. In this work, we show that altering the expression level of IRF6 dramatically modified this phenotype in utero. Developmental Dynamics 246:670-681, 2017. © 2017 Wiley Periodicals, Inc.

Generation and characterization of a conditional allele of Interferon Regulatory Factor 6.

Interferon Regulatory Factor 6 (IRF6) is a critical regulator of differentiation, proliferation, and migration of keratinocytes. Mutations in IRF6 cause two autosomal dominant disorders characterized by cleft lip with or without cleft palate. In addition, DNA variation in IRF6 confers significant risk for non-syndromic cleft lip and palate. IRF6 is also implicated in adult onset development and disease processes, including mammary gland development and squamous cell carcinoma. Mice homozygous for a null allele of Irf6 die shortly after birth due to severe skin, limb, and craniofacial defects, thus impeding the study of gene function after birth. To circumvent this, a conditional allele of Irf6 was generated. To validate the functionality of the conditional allele, we used three "deleter" Cre strains: Gdf9-Cre, CAG-Cre, and Ella-Cre. When Cre expression was driven by the Gdf9-Cre or CAG-Cre transgenes, 100% recombination was observed as indicated by DNA genotyping and phenotyping. In contrast, use of the Ella-Cre transgenic line resulted in incomplete recombination, despite expression at the one-cell stage. In sum, we generated a novel tool to delete Irf6 in a tissue specific fashion, allowing for study of gene function past perinatal stages. However, recombination efficiency of this allele was dictated by the Cre-driver used.

Dominant mutations in GRHL3 cause Van der Woude Syndrome and disrupt oral periderm development.

Mutations in interferon regulatory factor 6 (IRF6) account for ∼70% of cases of Van der Woude syndrome (VWS), the most common syndromic form of cleft lip and palate. In 8 of 45 VWS-affected families lacking a mutation in IRF6, we found coding mutations in grainyhead-like 3 (GRHL3). According to a zebrafish-based assay, the disease-associated GRHL3 mutations abrogated periderm development and were consistent with a dominant-negative effect, in contrast to haploinsufficiency seen in most VWS cases caused by IRF6 mutations. In mouse, all embryos lacking Grhl3 exhibited abnormal oral periderm and 17% developed a cleft palate. Analysis of the oral phenotype of double heterozygote (Irf6(+/-);Grhl3(+/-)) murine embryos failed to detect epistasis between the two genes, suggesting that they function in separate but convergent pathways during palatogenesis. Taken together, our data demonstrated that mutations in two genes, IRF6 and GRHL3, can lead to nearly identical phenotypes of orofacial cleft. They supported the hypotheses that both genes are essential for the presence of a functional oral periderm and that failure of this process contributes to VWS.

Toward an orofacial gene regulatory network.

Orofacial clefting is a common birth defect with significant morbidity. A panoply of candidate genes have been discovered through synergy of animal models and human genetics. Among these, variants in interferon regulatory factor 6 (IRF6) cause syndromic orofacial clefting and contribute risk toward isolated cleft lip and palate (1/700 live births). Rare variants in IRF6 can lead to Van der Woude syndrome (1/35,000 live births) and popliteal pterygium syndrome (1/300,000 live births). Furthermore, IRF6 regulates GRHL3 and rare variants in this downstream target can also lead to Van der Woude syndrome. In addition, a common variant (rs642961) in the IRF6 locus is found in 30% of the world's population and contributes risk for isolated orofacial clefting. Biochemical studies revealed that rs642961 abrogates one of four AP-2alpha binding sites. Like IRF6 and GRHL3, rare variants in TFAP2A can also lead to syndromic orofacial clefting with lip pits (branchio-oculo-facial syndrome). The literature suggests that AP-2alpha, IRF6 and GRHL3 are part of a pathway that is essential for lip and palate development. In addition to updating the pathways, players and pursuits, this review will highlight some of the current questions in the study of orofacial clefting.

An etiologic regulatory mutation in IRF6 with loss- and gain-of-function effects.

DNA variation in Interferon Regulatory Factor 6 (IRF6) causes Van der Woude syndrome (VWS), the most common syndromic form of cleft lip and palate (CLP). However, an etiologic variant in IRF6 has been found in only 70% of VWS families. To test whether DNA variants in regulatory elements cause VWS, we sequenced three conserved elements near IRF6 in 70 VWS families that lack an etiologic mutation within IRF6 exons. A rare mutation (350dupA) was found in a conserved IRF6 enhancer element (MCS9.7) in a Brazilian family. The 350dupA mutation abrogated the binding of p63 and E47 transcription factors to cis-overlapping motifs, and significantly disrupted enhancer activity in human cell cultures. Moreover, using a transgenic assay in mice, the 350dupA mutation disrupted the activation of MCS9.7 enhancer element and led to failure of lacZ expression in all head and neck pharyngeal arches. Interestingly, disruption of the p63 Motif1 and/or E47 binding sites by nucleotide substitution did not fully recapitulate the effect of the 350dupA mutation. Rather, we recognized that the 350dupA created a CAAAGT motif, a binding site for Lef1 protein. We showed that Lef1 binds to the mutated site and that overexpression of Lef1/ß-Catenin chimeric protein repressed MCS9.7-350dupA enhancer activity. In conclusion, our data strongly suggest that 350dupA variant is an etiologic mutation in VWS patients and disrupts enhancer activity by a loss- and gain-of-function mechanism, and thus support the rationale for additional screening for regulatory mutations in patients with CLP.

Interferon regulatory factor 6 regulates keratinocyte migration.

Interferon regulatory factor 6 (Irf6) regulates keratinocyte proliferation and differentiation. In this study, we tested the hypothesis that Irf6 regulates cellular migration and adhesion. Irf6-deficient embryos at 10.5 days post-conception failed to close their wound compared with wild-type embryos. In vitro, Irf6-deficient murine embryonic keratinocytes were delayed in closing a scratch wound. Live imaging of the scratch showed deficient directional migration and reduced speed in cells lacking Irf6. To understand the underlying molecular mechanisms, cell-cell and cell-matrix adhesions were investigated. We show that wild-type and Irf6-deficient keratinocytes adhere similarly to all matrices after 60 min. However, Irf6-deficient keratinocytes were consistently larger and more spread, a phenotype that persisted during the scratch-healing process. Interestingly, Irf6-deficient keratinocytes exhibited an increased network of stress fibers and active RhoA compared with that observed in wild-type keratinocytes. Blocking ROCK, a downstream effector of RhoA, rescued the delay in closing scratch wounds. The expression of Arhgap29, a Rho GTPase-activating protein, was reduced in Irf6-deficient keratinocytes. Taken together, these data suggest that Irf6 functions through the RhoA pathway to regulate cellular migration.

Exome sequence identifies RIPK4 as the Bartsocas-Papas syndrome locus.

Pterygium syndromes are complex congenital disorders that encompass several distinct clinical conditions characterized by multiple skin webs affecting the flexural surfaces often accompanied by craniofacial anomalies. In severe forms, such as in the autosomal-recessive Bartsocas-Papas syndrome, early lethality is common, complicating the identification of causative mutations. Using exome sequencing in a consanguineous family, we identified the homozygous mutation c.1127C>A in exon 7 of RIPK4 that resulted in the introduction of the nonsense mutation p.Ser376X into the encoded ankyrin repeat-containing kinase, a protein that is essential for keratinocyte differentiation. Subsequently, we identified a second mutation in exon 2 of RIPK4 (c.242T>A) that resulted in the missense variant p.Ile81Asn in the kinase domain of the protein. We have further demonstrated that RIPK4 is a direct transcriptional target of the protein p63, a master regulator of stratified epithelial development, which acts as a nodal point in the cascade of molecular events that prevent pterygium syndromes.

Abnormal skin, limb and craniofacial morphogenesis in mice deficient for interferon regulatory factor 6 (Irf6).

Transcription factor paralogs may share a common role in staged or overlapping expression in specific tissues, as in the Hox family. In other cases, family members have distinct roles in a range of embryologic, differentiation or response pathways (as in the Tbx and Pax families). For the interferon regulatory factor (IRF) family of transcription factors, mice deficient in Irf1, Irf2, Irf3, Irf4, Irf5, Irf7, Irf8 or Irf9 have defects in the immune response but show no embryologic abnormalities. Mice deficient for Irf6 have not been reported, but in humans, mutations in IRF6 cause two mendelian orofacial clefting syndromes, and genetic variation in IRF6 confers risk for isolated cleft lip and palate. Here we report that mice deficient for Irf6 have abnormal skin, limb and craniofacial development. Histological and gene expression analyses indicate that the primary defect is in keratinocyte differentiation and proliferation. This study describes a new role for an IRF family member in epidermal development.

Haploinsufficiency of interferon regulatory factor 6 alters brain morphology in the mouse.

Orofacial clefts are among the commonest birth defects. Among many genetic contributors to orofacial clefting, Interferon Regulatory Factor 6 (IRF6) is unique since mutations in this gene cause Van der Woude (VWS), the most common clefting syndrome. Furthermore, variants in IRF6 contribute to increased risk for non-syndromic cleft lip and/or palate (NSCL/P). Our previous work shows that individuals with either VWS or NSCL/P may have cerebral anomalies (larger anterior, smaller posterior regions), and a smaller cerebellum. The objective of this study was to test the hypothesis that disrupting Irf6 in the mouse will result in quantitative brain changes similar to those reported for humans with VWS and NSCL/P. Male mice heterozygous for Irf6 (Irf6(gt1/+); n = 9) and wild-type (Irf6(+/+) ; n = 6) mice at comparable age underwent a 4.7-T MRI scan to obtain quantitative measures of cortical and subcortical brain structures. There was no difference in total brain volume between groups. However, the frontal cortex was enlarged in the Irf6(gt1/+) mice compared to that of wild types (P = 0.028) while the posterior cortex did not differ. In addition, the volume of the cerebellum of Irf6(gt1/+) mice was decreased (P = 0.004). Mice that were heterozygous for Irf6 showed a similar pattern of brain anomalies previously reported in humans with VWS and NSCL/P. These structural differences were present in the absence of overt oral clefts. These results support a role for IRF6 in brain morphometry and provide evidence for a potential genetic link to abnormal brain development in orofacial clefting.

Pedigree structure and kinship measurements of a mid-Michigan community: a new North American population isolate identified.

Previous studies identified a cluster of individuals with an autosomal recessive form of deafness that reside in a small region of mid-Michigan. We hypothesized that affected members from this community descend from a defined founder population. Using public records and personal interviews, we constructed a genealogical database that includes the affected individuals and their extended families as descendants of 461 settlers who emigrated from the Eifel region of Germany between 1836 and 1875. The genealogical database represents a 13-generation pedigree that includes 27,747 descendants of these settlers. Among these descendants, 13,784 are presumed living. Many of the extant descendants reside in a 90-square-mile area, and 52% were born to parents who share at least one common ancestor. Among those born to related parents, the median kinship coefficient is 3.7 × 10(-3). While the pedigree contains 2,510 founders, 344 of the 461 settlers accounted for 67% of the genome in the extant population. These data suggest that we identified a new population isolate in North America and that, as demonstrated for congenital hearing loss, this rural mid-Michigan community is a new resource to discover heritable factors that contribute to common health-related conditions.

The role of Irf6 in tooth epithelial invagination.

Thickening and the subsequent invagination of the epithelium are an important initial step in ectodermal organ development. Ikkα has been shown to play a critical role in controlling epithelial growth, since Ikkα mutant mice show protrusions (evaginations) of incisor tooth, whisker and hair follicle epithelium rather than invagination. We show here that mutation of the Interferon regulatory factor (Irf) family, Irf6 also results in evagination of incisor epithelium. In common with Ikkα mutants, Irf6 mutant evagination occurs in a NF-κB-independent manner and shows the same molecular changes as those in Ikkα mutants. Irf6 thus also plays a critical role in regulating epithelial invagination. In addition, we also found that canonical Wnt signaling is upregulated in evaginated incisor epithelium of both Ikkα and Irf6 mutant embryos.

Exploring the effects of gene dosage on mandible shape in mice as a model for studying the genetic basis of natural variation.

Mandible shape in the mouse is a complex trait that is influenced by many genetic factors. However, little is known about the action of single genes on adult mandible shape so far, since most developmentally relevant genes are already required during embryogenesis, i.e., knockouts lead to embryonic death or severe deformations, before the mandible is fully formed. We employ here a geometric morphometric approach to identify subtle phenotypic differences caused by dosage effects of candidate genes. We use mouse strains with specific gene modifications (knockouts and knockins) to compare heterozygous animals with controls from the same stock, which is expected to be equivalent to a change of gene expression of the respective locus. Such differences in expression level are also likely to occur as part of the natural variation. We focus on Bmp pathway genes (Bmp4, its antagonist Noggin, and combinations of Bmp5-7 genotypes), but include also two other developmental control genes suspected to affect mandible development in some way (Egfr and Irf6). In addition, we study the effects of Hoxd13, as well as an extracellular matrix constituent (Col2a1). We find that subtle but significant shape differences are caused by differences in gene dosage of several of these genes. The changes seen for Bmp4 and Noggin are partially compatible with the action of these genes known from birds and fish. We find significant shape changes also for Hoxd13, although this gene has so far only been implicated in skeletal patterning processes of the limbs. Comparing the effect sizes of gene dosage changes to the variation found in natural populations of mice as well as quantitative trait loci (QTL) effects on mandible shape, we find that the effect sizes caused by gene dosage changes are at the lower end of the spectrum of natural variation, but larger than the average additive effects found in QTL studies. We conclude that studying gene dosage effects have the potential to provide new insights into aspects of craniofacial development, variation, and evolution.

Comparative analysis of IRF6 variants in families with Van der Woude syndrome and popliteal pterygium syndrome using public whole-exome databases.

PURPOSE: Mutations in the transcription factor IRF6 cause allelic autosomal dominant clefting syndromes, Van der Woude syndrome, and popliteal pterygium syndrome. We compared the distribution of IRF6 coding and splice-site mutations from 549 families with Van der Woude syndrome or popliteal pterygium syndrome with that of variants from the 1000 Genomes and National Heart, Lung, and Blood Institute Exome Sequencing Projects. METHODS: We compiled all published pathogenic IRF6 mutations and performed direct sequencing of IRF6 in families with Van der Woude syndrome or popliteal pterygium syndrome. RESULTS: Although mutations causing Van der Woude syndrome or popliteal pterygium syndrome were nonrandomly distributed with significantly increased frequencies in the DNA-binding domain (P = 0.0001), variants found in controls were rare and evenly distributed in IRF6. Of 194 different missense or nonsense variants described as potentially pathogenic, we identified only two in more than 6,000 controls. PolyPhen and SIFT (sorting intolerant from tolerant) reported 5.9% of missense mutations in patients as benign, suggesting that use of current in silico prediction models to determine function can have significant false negatives. CONCLUSION: Mutation of IRF6 occurs infrequently in controls, suggesting that for IRF6 there is a high probability that disruption of the coding sequence, particularly the DNA-binding domain, will result in syndromic features. Prior associations of coding sequence variants in IRF6 with clefting syndromes have had few false positives.