Gene Discovery in Nonsyndromic Cleft Lip with or without Cleft Palate

Nonsyndromic cleft lip with or without cleft palate (NSCLP), a common, complex orofacial birth defect that affects approximately 4,000 newborns each year in the United States, is caused by both genetic and environmental factors. Orofacial clefts affect the mouth and nose, causing severe deformity of the face, which require medical, dental and speech therapies. Despite having substantial genetic liability, less than 25% of the genetic contribute to NSCLP has been identified. The studies described in this thesis were performed to identify genes that contribute to NSCLP and to demonstrate the role of these genes in normal craniofacial development. Using genome scan and candidate gene approaches, novel associations with NSCLP were identified. These include MYH9 (7 SNPs, 0.009≤p<0.05), Wnt3A (4 SNPs, 0.001≤p≤0.005), Wnt11 (2 SNPs, 0.001≤p≤0.01) and CRISPLD2 (4 SNPs, 0.001≤p<0.05). The most interesting findings were for CRISPLD2. This gene is expressed in the fused mouse palate at E17.5. In zebrafish, crispld2 localized to the craniofacial region by one day post fertilization. Morpholino knockdown of crispld2 resulted in a lower survival rates and altered neural crest cell (NCC) clustering. Because NCCs form the tissues that populate the craniofacies, this NCC abnormality resulted in cartilage abnormalities of the jaw including fewer ceratobranchial cartilages forming the lower jaw (three pairs compared to five) and broader craniofacies compared to wild-type zebrafish. These findings suggest that the CRISPLD2 gene plays an important role in normal craniofacial development and perturbation of this gene in humans contributes to orofacial clefting. Overall, these results are important because they contribute to our understanding of normal craniofacial development and orofacial clefting etiology, information that can be used to develop better methods to diagnose, counsel and potentially treat NSCLP patients.

Mutation and association analysis of the PVR and PVRL2 genes in patients with non-syndromic cleft lip and palate.

Orofacial clefts (OFC; MIM 119530) are among the most common major birth defects. Here, we carried out mutation screening of the PVR and PVRL2 genes, which are both located at an OFC linkage region at 19q13 (OFC3) and are closely related to PVRL1, which has been associated with both syndromic and non-syndromic cleft lip and palate (nsCLP). We screened a total of 73 nsCLP patients and 105 non-cleft controls from the USA for variants in PVR and PVRL2, including all exons and encompassing all isoforms. We identified four variants in PVR and five in PVRL2. One non-synonymous PVR variant, A67T, was more frequent among nsCLP patients than among normal controls, but this difference did not achieve statistical significance.

Mutation analysis of the PVRL1 gene in caucasians with nonsyndromic cleft lip/palate.

Nonsyndromic cleft lip with or without cleft palate (nsCL/P, MIM 119530) is perhaps the most common major birth defect. Homozygous PVRL1 loss-of-function mutations result in an autosomal recessive CL/P syndrome, CLPED1, and a PVRL1 nonsense mutation is associated with sporadic nsCL/P in Northern Venezuela. To address the more general role of PVRL1 variation in risk of nsCL/P, we carried out mutation analysis of PVRL1 in North American and Australian nsCL/P cases and population-matched controls. We identified a total of 15 variants, 5 of which were seen in both populations and 1 of which, an in-frame insertion at Glu442, was more frequent in patients than in controls in both populations, though the difference was not statistically significant. Another variant, which is specific to the PVRL1 beta (HIgR) isoform, S447L, was marginally associated with nsCL/P in North American Caucasian patients, but not in Australian patients, and overall variants that affect the beta-isoform were significantly more frequent among North American patients. One Australian patient had a splice junction mutation of PVRL1. Our results suggest that PVRL1 may play a minor role in susceptibility to the occurrence of nsCL/P in some Caucasian populations, and that variation involving the beta (HIgR) isoform might have particular importance for risk of orofacial clefts. Nevertheless, these results underscore the need for studies that involve very large numbers when assessing the possible role of rare variants in risk of complex traits such as nsCL/P.

Genomic screening identifies novel linkages and provides further evidence for a role of MYH9 in nonsyndromic cleft lip and palate.

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth anomaly that requires prolonged multidisciplinary rehabilitation. Although variation in several genes has been identified as contributing to NSCLP, most of the genetic susceptibility loci have yet to be defined. To identify additional contributory genes, a high-throughput genomic scan was performed using the Illumina Linkage IVb Panel platform. We genotyped 6008 SNPs in nine non-Hispanic white NSCLP multiplex families and a single large African-American NSCLP multiplex family. Fourteen chromosomal regions were identified with LOD>1.5, including six regions not previously reported. Analysis of the data from the African-American and non-Hispanic white families revealed two likely chromosomal regions: 8q21.3-24.12 and 22q12.2-12.3 with LOD scores of 2.98 and 2.66, respectively. On the basis of biological function, syndecan 2 (SDC2) and growth differentiation factor 6 (GDF6) in 8q21.3-24.12 and myosin heavy-chain 9, non-muscle (MYH9) in 22q12.2-12.3 were selected as candidate genes. Association analyses from these genes yielded marginally significant P-values for SNPs in SDC2 and GDF6 (0.01

Familial component of epilepsy in cleft lip and palate

It is well recognized that offspring of women with epilepsy who are taking anticonvulsant medications have an increased incidence of clefting abnormalities. This increase has been attributed to the teratogenic effects of anticonvulsant medications but an alternative explanation involving a genetic association of epilepsy and clefting has also been proposed. Five family studies attempting to resolve this controversy have been inconclusive either because of study design or analytic limitations. This family study was designed to determine whether epilepsy aggregates in families ascertained by an individual with a clefting disorder. The Mayo Clinic medical linkage registry was used to identify individuals with cleft lip with or without cleft palate and cleft palate in southeast Minnesota from 1935-1986. Only those cases who were 15 years or younger during this period were included in the study. The proband's parents and descendants of their parents, including the proband's sibs, children, grandchildren, niece/nephews, grandnieces/nephews, halfsibs and spouses were also identified and all of their medical records were reviewed for seizure disorders. The standardized morbidity ratios for epilepsy of 0.9 (95% CI 0.2-2.6) observed for first degree relatives (excluding parents) and 0.0 for second degree relatives were not increased. The SMRs ranged from 0.7-2.2 for the individual relative types (parents 1.5, sibs 0.7, children 2.2, probands 1.1, spouses 2.0) and were also not increased. These results do not support the suggestions of some that clefting and epilepsy aggregate together in families. ^

FZD6, MATN2 AND SLC25A32, POSSIBLE CANDIDATE GENES IN NONSYNDROMIC CLEFT LIP AND PALATE

Nonsyndromic cleft lip with or without cleft palate (NSCLP) is a common birth defect with a multifactorial etiology. Despite decades of research, the genetic underpinnings of NSCLP still remain largely unexplained. A genome wide association study (GWAS) of a large NSCLP African American family with seven affected individuals across three generations found evidence for linkage at 8q21.3-24.12 (LOD = 2.98). This region contained three biologically relevant candidate genes: Frizzled-6 (FZD6) (LOD = 2.8), Matrilin-2 (MATN2) (LOD = 2.3), and Solute Carrier Family 25, Member 32 (SLC26A32) (LOD = 1.6). Sequencing of the coding regions and the 5’ and 3’ UTRs of these genes in two affected family members identified a rare intronic variant, rs138557689 (c.-153+432A>C), in FZD6. The rs138557689/C allele segregated with the NSCLP phenotype; in silico analysis predicted and EMSA analysis showed that the 138557689/C allele creates new DNA binding sites. FZD6 is part of the WNT pathway, which is involved in craniofacial development, including midface development and upper lip fusion. Our novel findings suggest that an alteration in FZD6 gene regulation may perturb this tightly controlled biological pathway and in turn contribute to the development of NSCLP in this family. Studies are underway to further define how the rs138557689/C variant affects expression of FZD6.

Maturation of oral feeding skills in preterm infants.

AIM: Safe and successful oral feeding requires proper maturation of sucking, swallowing and respiration. We hypothesized that oral feeding difficulties result from different temporal development of the musculatures implicated in these functions. METHODS: Sixteen medically stable preterm infants (26 to 29 weeks gestation, GA) were recruited. Specific feeding skills were monitored as indirect markers for the maturational process of oral feeding musculatures: rate of milk intake (mL/min); percent milk leakage (lip seal); sucking stage, rate (#/s) and suction/expression ratio; suction amplitude (mmHg), rate and slope (mmHg/s); sucking/swallowing ratio; percent occurrence of swallows at specific phases of respiration. Coefficients of variation (COV) were used as indices of functional stability. Infants, born at 26/27- and 28/29-week GA, were at similar postmenstrual ages (PMA) when taking 1-2 and 6-8 oral feedings per day. RESULTS: Over time, feeding efficiency and several skills improved, some decreased and others remained unchanged. Differences in COVs between the two GA groups demonstrated that, despite similar oral feeding outcomes, maturation levels of certain skills differed. CONCLUSIONS: Components of sucking, swallowing, respiration and their coordinated activity matured at different times and rates. Differences in functional stability of particular outcomes confirm that maturation levels depend on infants' gestational rather than PMA.

In vitro analog of classical conditioning of feeding behavior in aplysia.

The feeding behavior of Aplysia californica can be classically conditioned using tactile stimulation of the lips as a conditioned stimulus (CS) and food as an unconditioned stimulus (US). Moreover, several neural correlates of classical conditioning have been identified. The present study extended previous work by developing an in vitro analog of classical conditioning and by investigating pairing-specific changes in neuronal and synaptic properties. The preparation consisted of the isolated cerebral and buccal ganglia. Electrical stimulation of a lip nerve (AT4) and a branch of the esophageal nerve (En2) served as the CS and US, respectively. Three protocols were used: paired, unpaired, and US alone. Only the paired protocol produced a significant increase in CS-evoked fictive feeding. At the cellular level, classical conditioning enhanced the magnitude of the CS-evoked synaptic input to pattern-initiating neuron B31/32. In addition, paired training enhanced both the magnitude of the CS-evoked synaptic input and the CS-evoked spike activity in command-like neuron CBI-2. The in vitro analog of classical conditioning reproduced all of the cellular changes that previously were identified following behavioral conditioning and has led to the identification of several new learning-related neural changes. In addition, the pairing-specific enhancement of the CS response in CBI-2 indicates that some aspects of associative plasticity may occur at the level of the cerebral sensory neurons.

BMP-SIGNALING REGULATES A COMMON TRANSCRIPTIONAL PROGRAM TO CONTROL FACIAL FORM AND SKELETAL MORPHOGENESIS

Much of the craniofacial skeleton, such as the skull vault, mandible and midface, develops through direct, intramembranous ossification of the cranial neural crest (CNC) derived progenitor cells. Bmp-signaling plays critical roles in normal craniofacial development, and Bmp4 deficiency results in craniofacial abnormalities, such as cleft lip and palate. We performed an in depth analysis of Bmp4, a critical regulator of development, disease, and evolution, in the CNC. Conditional Bmp4 overexpression, using a tetracycline regulated Bmp4 gain of function allele, resulted in facial form changes that were most dramatic after an E10.5 Bmp4 induction. Expression profiling uncovered a signature of Bmp4 induced genes (BIG) composed predominantly of transcriptional regulators controlling self-renewal, osteoblast differentiation, and negative Bmp autoregulation. The complimentary experiment, CNC inactivation of Bmp2, Bmp4, and Bmp7, resulted in complete or partial loss of multiple CNC derived skeletal elements revealing a critical requirement for Bmp-signaling in membranous bone and cartilage development. Importantly, the BIG signature was reduced in Bmp loss of function mutants indicating similar Bmp-regulated target genes underlying facial form modulation and normal skeletal morphogenesis. Chromatin immunoprecipitation (ChIP) revealed a subset of the BIG signature, including Satb2, Smad6, Hand1, Gadd45g and Gata3 that was bound by Smad1/5 in the developing mandible revealing direct, Smad-mediated regulation. These data indicate that Bmp-signaling regulates craniofacial skeletal development and facial form by balancing self-renewal and differentiation pathways in CNC progenitors.

Birth defects among surviving children under seven years of age in Tianjin, China

Birth defects are a leading cause of infant mortality in the developed countries. They are also of increasing concern in many developing countries, such as China. However, prevalence and causes of birth defects in China are inadequately understood.^ The purpose of the present study was to estimated prevalence of birth defects in surviving children under seven years of age in Tianjin, China and investigate determinants of birth defects in the study area.^ The present study took place in Tianjin, China in 1986, involving 22,081 surviving children under seven years of age. Children with birth defects were ascertained through physical examinations by physicians during household visits and ascertainment of birth defects was verified through multiple sources. Of 22,081 surviving children, 524 had birth defects (23.7 per 1,000). The study noted a striking discrepancy in the prevalence of birth defects between urban and rural area. The prevalence of birth defects was 16.3 per 1,000 in the urban and 33.2 per 1,000 in the rural area.^ Using cases of birth defects ascertained from surviving children, a case-control study was carried out. The study observed that first-trimester maternal flu was associated with increased risk of both major and minor birth defects in children after controlling for other maternal factors (adjusted odds ratio (OR) = 8.7, 95% confidence interval (CI) = 4.3-17.3; OR = 3.6, 95% CI = 1.7-7.5). This association could be biased by different reporting of exposure between mothers of children with birth defects and mothers of children without defects. This study indicated that maternal flu was also associated with congenital heart defects and polydactyly after controlling for other maternal factors (adjusted OR = 32.3, 95% CI = 13.3-78.3; adjusted OR = 5.5, 95% CI = 1.1-27.7). The associations remained when affected controls (children with similar birth defects other than congenital heart defects or polydactyly) were used (adjusted OR = 4.3, 95% CI = 1.2-15.3; OR = 1.4, 95% CI = 1.4-7.9). A weak association between first-trimester vaginal bleeding and selected groups of birth defects was found in this study, but the association may be confounded by other factors. Maternal smoking during pregnancy was modestly associated with cleft lip with or without cleft palate (OR = 1.4, 95% = 0.4-4.9), but the association may be due to chance. Some major limitations in this study warrant caution in interpretation of the findings, especially the causal relation. ^

THE EQUILIBRIUM CONSTANT FOR THE GLYCINE SYNTHASE REACTION (ONE-CARBON, METABOLISM, FOLATE)

The equilibrium constant (K(,c)) under physiological conditions (38(DEGREES)C, 0.25 M ionic strength (I), pH 7.0) for the glycine synthase (GS) reaction (E C 2.1.2.1.0) (Equation 1) has been determined. (UNFORMATTED TABLE FOLLOWS)^ 5,10-CH(,2)-H(,4)Folate NADH NH (,4)+ CO(,2) ^ K(,c) = Eq. 1^ H(,4)Folate NAD('+) GLY ^(TABLE ENDS)^ The enzymatic instability of the GS enzyme complex itself has made it necessary to determine the overall K(,c) from the product of constants for the partial reactions of GS determined separately under the same conditions. The partial reactions are the H(,4)Folate-formaldehyde (CH(,2)(OH)(,2)) condensation reaction (Reaction 1) the K(,c) for which has been reported by this laboratory (3.0 x 10('4)), the lipoate (LipS(,2)) dehydrogenase reaction (LipDH) (Reaction 2) and the Gly-Lip^ decarboxylase reaction (Reaction 3) forming reduced lipoate (Lip(SH)(,2)), NH(,4)('+), CO(,2) and CH(,2)(OH)(,2.) (UNFORMATTED TABLE FOLLOWS)(,)^ H(,4)Fote + CH(,2)(OH)(,2) 5,10-CH(,2)-H(,4)Folate (1)^ Lip(SH)(,2) + NAD('+) LipS(,2) + NADH + H('+) (2)^ H('+) + Gly + LipS(,2) Lip(SH)(,2) + NH(,4)('+) CO(,2) + CH(,2)(OH)(,2) (3)^(TABLE ENDS)^ In this work the K(,c) for Reactions 2 and 3 are reported.^ The K(,c)' for the LipDH reaction described by other authors was reported with unexplainable conclusions regarding the pH depend- ence for the reaction. These conclusions would imply otherwise unexpected acid dissociation constants for reduced and oxidized lipoate. The pK(,a)',s for these compounds have been determined to resolve discrepancy. The conclusions are as follows: (1) The K(,c) for the LipDH reaction is 2.08 x 10('-8); (2) The pK(,a)',s for Lip(SH)(,2) are 4.77(-COOH), 9.91(-SH), 11.59(-SH); for LipS(,2) the carboxyl pK(,a)' is 4.77; (3) Contrary to previous literature, the log K(,c)' for the LipDH reaction is a linear function of the pH, a conclusion supported by the values for the dissociation constants.^ The K(,c) for Reaction 3 is the product of constants for Reactions 4-7. (UNFORMATTED TABLE FOLLOWS)^ LipSHSCH(,2)OH + H(,2)O Lip(SH)(,2) + CH(,2)(OH)(,2) (4)^ H(,2)O + LipSHSCH(,2)NH(,3)('+) LipSHSCH(,2)OH + NH(,4)('+) (5)^ LipSHSCH(,2)NH(,2) + H('+) LipSHSCH(,2)NH(,3)('+) (6)^ Gly + LipS(,2) LipSHSCH(,2)NH(,2) + CO(,2) (7)^(TABLE ENDS)^ Reactions 4-6 are non-enzymatic reactions whose constants were determined spectrophotometrically. Reaction 7 was catalyzed by the partially purified P-protein of GS with equilibrium approached from both directions. The value for K(,c) for this reaction is 8.15 x 10('-3). The combined K(,c) for Reactions 4-7 or Reaction 3 is 2.4 M.^ The overall K(,c) for the GS reaction determined by combination of values for Reactions 1-3 is 1.56 x 10('-3). ^

The functions of pitx2 and Bmp signaling in craniofacial development

The formation of the vertebrate face is an extremely complex developmental process, which needs to coordinate the outgrowth of several facial primordia. Facial primordia are small buds made up of mesenchymal masses enclosed by an epithelial layer that surrounds the primitive mouth. The upper jaw is formed by the maxillary process, the lateral nasal process, and the frontonasal process while the mandibular process forms the lower jaw. Recent experiments using genetics in mice and bead implantation approaches have shown that the pitx2 homeobox gene and Bmp signaling play important roles in this complex developmental process. However, the molecular mechanisms underlying the function of pitx2 and Bmp in these events are still unclear. Here, we show that pitx2 is required for oral epithelium maintenance, and branchial arch signaling is pitx2 dosage sensitive by using pitx2 allelic combinations that encode varying levels of pitx2. Maintenance of fgf8 signaling requires only low pitx2 dosage while repression of Bmp signaling requires high pitx2 levels. Different incisor and molar phenotypes in low level pitx2 mutant embryos suggest a distinct requirement for pitx2 in tooth-type development. The results show that pitx2 is required for craniofacial muscle formation and expanded Bmp signaling results in excess bone formation in pitx2 mutant embryos. Fate-mapping studies show that ectopic bone results from excessive bone growth, instead of muscle transformation. Moreover, by using cre/loxp system we show that partial loss of Bmpr-IA in the facial primordia results in cleft lip/palate, abnormal teeth, ectopic teeth and tooth transformation. These phenotypes suggest that Bmp signaling has multiple functions during craniofacial development. The mutant palate shelves can fuse with each other when cultured in vitro, suggesting that cleft palate is secondary to the partial loss of Bmpr-IA. Furthermore, we prove that Bmp4, one of the ligands of Bmpr-IA, plays a role during lip fusion developmental process and partial loss of Bmp4 in the facial primordia results in the lip fusion delay. These results have provided insight to understand the complex signaling cascades that regulate craniofacial development. ^

The role of MEKK3 in T cell homeostasis and IFNgamma production

MEKK3, a member of the MAP3K family, is involved in regulating multiple MAPK and NF-κB pathways. The MAPK and NF-κB signaling pathways are important in regulating T cell functions. MEKK3 is expressed through the development of T cell and also in subsets of T cell in the peripheral. However, the specific role of MEKK3 in T cell function is unknown. To reveal the in vivo function of MEKK3 in T cells, I have generated MEKK3 T cell conditional knock-out mice. Despite a normal thymus development in the conditional knock-out mice, I observed a decrease in the number of peripheral T-cells and impaired T-cell function in response to antigen stimulation. T cells undergo homeostatic proliferation under lymphopenia condition, a process called lymphopenia-induced proliferation (LIP). Using a LIP model, I demonstrated that the reduction of peripheral T cell number is largely due to a severe impairment of the self-antigen/MHC mediated T cell homeostasis. Upon anti-CD3 stimulation, the proliferation of MEKK3-deficient T cell is not significantly affected, but the production of IFNγ by naïve and effector CD4 T cells are markedly decreased. Interestingly, the IL-12/IL-18 driven IFNγ production and MAPK activation in MEKK3-deficient T cells is not affected, suggesting that MEKK3 selectively mediates the TCR induced MAPK signaling. Furthermore, I found that MEKK3 is activated by TCR stimulation in a RAC1/2 dependent manner, but not by IL-12/IL-18 stimulation. Finally, I showed that basal level of ERK and JNK activation is defective under LIP condition. I showed that the TCR induced ERK, JNK and p38 MAPK activation is also defective in MEKK3 deficient CD4 T cells. Taken together, my data demonstrate a crucial role of MEKK3 in T cell homeostasis and IFNγ production through regulating the TCR mediated MAPK pathway. ^

A Bayesian analysis for spatial processes with application to mapping

In geographical epidemiology, maps of disease rates and disease risk provide a spatial perspective for researching disease etiology. For rare diseases or when the population base is small, the rate and risk estimates may be unstable. Empirical Bayesian (EB) methods have been used to spatially smooth the estimates by permitting an area estimate to "borrow strength" from its neighbors. Such EB methods include the use of a Gamma model, of a James-Stein estimator, and of a conditional autoregressive (CAR) process. A fully Bayesian analysis of the CAR process is proposed. One advantage of this fully Bayesian analysis is that it can be implemented simply by using repeated sampling from the posterior densities. Use of a Markov chain Monte Carlo technique such as Gibbs sampler was not necessary. Direct resampling from the posterior densities provides exact small sample inferences instead of the approximate asymptotic analyses of maximum likelihood methods (Clayton & Kaldor, 1987). Further, the proposed CAR model provides for covariates to be included in the model. A simulation demonstrates the effect of sample size on the fully Bayesian analysis of the CAR process. The methods are applied to lip cancer data from Scotland, and the results are compared. ^

Maternal exposure to hazardous air pollutants and oral clefts among offspring: Texas, 1999-2008

Birth defects are a leading cause of infant mortality in the United States. About one in 33 births in the United States is diagnosed with birth defects. Common birth defects include neural tube defects, Down syndrome and oral clefts. The present study focused on oral clefts. ^ Oral clefts refer to the malformation of lip, mouth or both. Birth prevalence of oral clefts in Texas is about 11 per 10,000 births. Etiologically, oral clefts have been classified into two groups, cleft lip with or without cleft palate (CL±P) and isolated cleft palate (CP). In spite of their high prevalence and clinical significance, the etiology of oral clefts in humans has not been well understood. Though a number of risk factors have been identified in epidemiological studies, most of them do not explain the majority of the cases. The need to identify novel risk factors associated with oral clefts provided the motivation for this study. The present study focused on maternal exposure to several hazardous air pollutants. A common subgroup of hazardous air pollutants is the volatile organic compounds found in petroleum derivatives. Four important hydrocarbons in this group are benzene, toluene, ethyl benzene and xylenes (BTEX). ^ The specific aim of this study was to evaluate the association between maternal exposure to environmental levels of BTEX and oral clefts among offspring in Texas for the period 1999-2008. ^ A case-control study design was used to assess if maternal exposure to BTEX increased the risk of oral clefts. The Texas Birth Defects Registry provided data on cases of non-syndromic oral clefts delivered in Texas during the period 1999-2008. Census tract level maternal exposure to BTEX concentrations were obtained from the Hazardous Air Pollutant Exposure Model (HAPEM) developed by the U.S. Environmental Protection Agency. Unconditional logistic regression was used to assess the relationship between maternal exposure to BTEX levels and risk of oral clefts in offspring. ^ In the selected population, mothers who had high estimated exposure to any of the BTEX compounds were not more likely to deliver an offspring with oral clefts. Future research efforts will focus on additional birth defects and thorough assessment of additional potential confounders.^