Strain-specific modifier genes governing craniofacial phenotypes

The Bone Sialoprotein RGD Domain Modulates and Maintains Periodontal Development

Bone sialoprotein (gene: Ibsp; protein: BSP) is a multifunctional extracellular matrix protein present in bone, cementum, and dentin. Accumulating evidence supports BSP as a key regulator of mineralized tissue formation via evolutionarily conserved functional domains, including a C-terminal integrin-binding Arg-Gly-Asp (RGD) domain implicated in extracellular matrix-cell signaling. Ablation of Ibsp in mice (Ibsp-/-) results in impaired bone growth and mineralization and defective osteoclastogenesis, with effects in the craniofacial region including reduced acellular cementum formation, detachment of the periodontal ligament (PDL), alveolar bone hypomineralization, and severe periodontal breakdown. We hypothesized that BSP-RGD plays an important role in cementum and alveolar bone formation and mineralization, as well as periodontal function. This hypothesis was tested by replacing the RGD motif with a nonfunctional Lys-Ala-Glu (KAE) sequence in (IbspKAE/KAE) mice and OCCM.30 murine (IbspKAE) cementoblasts. The RGD domain was not critical for acellular or cellular cementum formation in IbspKAE/KAE mice. However, PDL volume and thickness were increased, and significantly more tartrate-resistant acid phosphatase-positive osteoclasts were found on alveolar bone surfaces of IbspKAE/KAE mice versus wild type mice. PDL organization was disrupted as indicated by picrosirius red stain, second harmonic generation imaging, dynamic mechanical analysis, and decreased asporin proteoglycan localization. In vitro studies implicated RGD functions in cell migration, adhesion, and mineralization, and this was confirmed by an ossicle implant model where cells lacking BSP-RGD showed substantial defects as compared with controls. In total, the BSP-RGD domain is implicated in periodontal development, though the scale and scope of changes indicated by in vitro studies indicate that other factors may partially compensate for and reduce the phenotypic severity of mice lacking BSP-RGD in vivo.

Gpr63 is a modifier of microcephaly in Ttc21b mouse mutants

The primary cilium is a signaling center critical for proper embryonic development. Previous studies have demonstrated that mice lacking Ttc21b have impaired retrograde trafficking within the cilium and multiple organogenesis phenotypes, including microcephaly. Interestingly, the severity of the microcephaly in Ttc21b^aln/aln homozygous null mutants is considerably affected by the genetic background and mutants on an FVB/NJ (FVB) background develop a forebrain significantly smaller than mutants on a C57BL/6J (B6) background. We performed a Quantitative Trait Locus (QTL) analysis to identify potential genetic modifiers and identified two regions linked to differential forebrain size: modifier of alien QTL1 (Moaq1) on chromosome 4 at 27.8 Mb and Moaq2 on chromosome 6 at 93.6 Mb. These QTLs were validated by constructing congenic strains. Further analysis of Moaq1 identified an orphan G-protein coupled receptor (GPCR), Gpr63, as a candidate gene. We identified a SNP that is polymorphic between the FVB and B6 strains in Gpr63 and creates a missense mutation predicted to be deleterious in the FVB protein. We used CRISPR-Cas9 genome editing to create two lines of FVB congenic mice: one with the B6 sequence of Gpr63 and the other with a deletion allele leading to a truncation of the GPR63 C-terminal tail. We then demonstrated that Gpr63 can localize to the cilium in vitro. These alleles affect ciliary localization of GPR63 in vitro and genetically interact with Ttc21b^aln/aln as Gpr63;Ttc21b double mutants show unique phenotypes including spina bifida aperta and earlier embryonic lethality. This validated Gpr63 as a modifier of multiple Ttc21b neural phenotypes and strongly supports Gpr63 as a causal gene (i.e., a quantitative trait gene, QTG) within the Moaq1 QTL.

TTC21B in humans is a known ciliopathy gene and contributes to the pathophysiology of a number of ciliopathies. Mice homozygous for a null allele of Ttc21b also have a spectrum of ciliopathy phenotypes, including microcephaly (small brain). Further work has shown that the severity of the microcephaly significantly depends on the genetic background of the mouse model. The genetic mechanisms of the Ttc21b pathophysiology and the interacting gene network remain far from understood. As an initial attempt to understand the underlying mechanism(s) underlying the variable effects on brain size, we performed a quantitative trait locus (QTL) analysis and found two regions of genomic significance that correlated with smaller brain size. We confirmed both QTLs with congenic lines. One of the two regions was small enough that we considered candidate genes and hypothesized Gpr63 might be a contributing locus for a number of reasons. We evaluated this hypothesis directly with precise variant creation using genome editing and provide evidence that Ttc21b and Gpr63 do indeed genetically interact. Thus, we have been able to combine classical QTL analysis and genome editing to directly test the resulting hypothesis.

Strain-specific differences in brain gene expression in a hydrocephalic mouse model with motile cilia dysfunction

Congenital hydrocephalus results from cerebrospinal fluid accumulation in the ventricles of the brain and causes severe neurological damage, but the underlying causes are not well understood. It is associated with several syndromes, including primary ciliary dyskinesia (PCD), which is caused by dysfunction of motile cilia. We previously demonstrated that mouse models of PCD lacking ciliary proteins CFAP221, CFAP54 and SPEF2 all have hydrocephalus with a strain-dependent severity. While morphological defects are more severe on the C57BL/6J (B6) background than 129S6/SvEvTac (129), cerebrospinal fluid flow is perturbed on both backgrounds, suggesting that abnormal cilia-driven flow is not the only factor underlying the hydrocephalus phenotype. Here, we performed a microarray analysis on brains from wild type and nm1054 mice lacking CFAP221 on the B6 and 129 backgrounds. Expression differences were observed for a number of genes that cluster into distinct groups based on expression pattern and biological function, many of them implicated in cellular and biochemical processes essential for proper brain development. These include genes known to be functionally relevant to congenital hydrocephalus, as well as formation and function of both motile and sensory cilia. Identification of these genes provides important clues to mechanisms underlying congenital hydrocephalus severity.

Determinants of orofacial clefting I: Effects of 5-Aza-2'-deoxycytidine on cellular processes and gene expression during development of the first branchial arch

In this study, we identify gene targets and cellular events mediating the teratogenic action(s) of 5-Aza-2'-deoxycytidine (AzaD), an inhibitor of DNA methylation, on secondary palate development. Exposure of pregnant mice (on gestation day (GD) 9.5) to AzaD for 12h resulted in the complete penetrance of cleft palate (CP) in fetuses. Analysis of cells of the embryonic first branchial arch (1-BA), in fetuses exposed to AzaD, revealed: 1) significant alteration in expression of genes encoding several morphogenetic factors, cell cycle inhibitors and regulators of apoptosis; 2) a decrease in cell proliferation; and, 3) an increase in apoptosis. Pyrosequencing of selected genes, displaying pronounced differential expression in AzaD-exposed 1-BAs, failed to reveal significant alterations in CpG methylation levels in their putative promoters or gene bodies. CpG methylation analysis suggested that the effects of AzaD on gene expression were likely indirect.