Felicia Cheah

Evidence for gene-environment interaction in a genome wide study of nonsyndromic cleft palate.

Nonsyndromic cleft palate (CP) is a common birth defect with a complex and heterogeneous etiology involving both genetic and environmental risk factors. We conducted a genome‐wide association study (GWAS) using 550 case‐parent trios, ascertained through a CP case collected in an international consortium. Family‐based association tests of single nucleotide polymorphisms (SNP) and three common maternal exposures (maternal smoking, alcohol consumption, and multivitamin supplementation) were used in a combined 2 df test for gene (G) and gene‐environment (G × E) interaction simultaneously, plus a separate 1 df test for G × E interaction alone. Conditional logistic regression models were used to estimate effects on risk to exposed and unexposed children. While no SNP achieved genome‐wide significance when considered alone, markers in several genes attained or approached genome‐wide significance when G × E interaction was included. Among these, MLLT3 and SMC2 on chromosome 9 showed multiple SNPs resulting in an increased risk if the mother consumed alcohol during the peri‐conceptual period (3 months prior to conception through the first trimester). TBK1 on chr. 12 and ZNF236 on chr. 18 showed multiple SNPs associated with higher risk of CP in the presence of maternal smoking. Additional evidence of reduced risk due to G × E interaction in the presence of multivitamin supplementation was observed for SNPs in BAALC on chr. 8. These results emphasize the need to consider G × E interaction when searching for genes influencing risk to complex and heterogeneous disorders, such as nonsyndromic CP. Genet. Epidemiol. 2011.  © 2011 Wiley‐Liss, Inc. 35: 469‐478, 2011

Analysis of candidate genes on chromosome 2 in oral cleft case-parent trios from three populations

Isolated oral clefts, including cleft lip with/without cleft palate (CL/P) and cleft palate (CP), have a complex and heterogeneous etiology. Case-parent trios from three populations were used to study genes spanning chromosome 2, where single nucleotide polymorphic (SNP) markers were analyzed individually and as haplotypes. Case-parent trios from three populations (74 from Maryland, 64 from Singapore and 95 from Taiwan) were genotyped for 962 SNPs in 104 genes on chromosome 2, including two well-recognized candidate genes: TGFA and SATB2. Individual SNPs and haplotypes (in sliding windows of 2–5 SNPs) were used to test for linkage and disequilibrium separately in CL/P and CP trios. A novel candidate gene (ZNF533) showed consistent evidence of linkage and disequilibrium in all three populations for both CL/P and CP. SNPs in key regions of ZNF533 showed considerable variability in estimated genotypic odds ratios and their significance, suggesting allelic heterogeneity. Haplotype frequencies for regions of ZNF533 were estimated and used to partition genetic variance into among-and within-population components. Wright’s fixation index, a measure of genetic diversity, showed little difference between Singapore and Taiwan compared with Maryland. The tensin-1 gene (TNS1) also showed evidence of linkage and disequilibrium among both CL/P and CP trios in all three populations, albeit at a lower level of significance. Additional genes (VAX2, GLI2, ZHFX1B on 2p; WNT6–WNT10A and COL4A3–COL4A4 on 2q) showed consistent evidence of linkage and disequilibrium only among CL/P trios in all three populations, and TGFA showed significant evidence in two of three populations.

Differential parental transmission of markers in RUNX2 among cleft case-parent trios from four populations

Isolated cleft lip with or without cleft palate (CL/P) is among the most common human birth defects, with a prevalence around 1 in 700 live births. The Runt‐related transcription factor 2 (RUNX2) gene has been suggested as a candidate gene for CL/P based largely on mouse models; however, no human studies have focused on RUNX2 as a risk factor for CL/P. This study examines the association between markers in RUNX2 and isolated, nonsyndromic CL/P using a case‐parent trio design, while considering parent‐of‐origin effects. Case‐parent trios from four populations (77 from Maryland, 146 from Taiwan, 35 from Singapore, and 40 from Korea) were genotyped for 24 single nucleotide polymorphisms (SNPs) in the RUNX2 gene. We performed the transmission disequilibrium test on individual SNPs. Parent‐of‐origin effects were assessed using the transmission asymmetry test and the parent‐of‐origin likelihood ratio test (PO‐LRT). When all trios were combined, the transmission asymmetry test revealed a block of 11 SNPs showing excess maternal transmission significant at the P<0.01 level, plus one SNP (rs1934328) showing excess paternal transmission (P=0.002). For the 11 SNPs showing excess maternal transmission, odds ratios of being transmitted to the case from the mother ranged between 3.00 and 4.00. The parent‐of‐origin likelihood ratio tests for equality of maternal and paternal transmission were significant for three individual SNPs (rs910586, rs2819861, and rs1934328). Thus, RUNX2 appears to influence risk of CL/P through a parent‐of‐origin effect with excess maternal transmission. Genet. Epidemiol. 2008.

The FGF and FGFR Gene Family and Risk of Cleft Lip With or Without Cleft Palate

Background Isolated, nonsyndromic cleft lip with or without cleft palate is a common human congenital malformation with a complex and heterogeneous etiology. Genes coding for fibroblast growth factors and their receptors (FGF/FGFR genes) are excellent candidate genes. Methods We tested single-nucleotide polymorphic markers in 10 FGF/FGFR genes (including FGFBP1, FGF2, FGF10, FGF18, FGFR1, FGFR2, FGF19, FGF4, FGF3, and FGF9) for genotypic effects, interactions with one another, and with common maternal environmental exposures in 221 Asian and 76 Maryland case-parent trios ascertained through a child with isolated, nonsyndromic cleft lip with or without cleft palate. Results Both FGFR1 and FGF19 yielded evidence of linkage and association in the transmission disequilibrium test, confirming previous evidence. Haplotypes of three single-nucleotide polymorphisms in FGFR1 were nominally significant among Asian trios. Estimated odds ratios for individual single-nucleotide polymorphic markers and haplotypes of multiple markers in FGF19 ranged from 1.31 to 1.87. We also found suggestive evidence of maternal genotypic effects for markers in FGF2 and FGF10 among Asian trios. Tests for gene-environment (G x E) interaction between markers in FGFR2 and maternal smoking or multivitamin supplementation yielded significant evidence of G x E interaction separately. Tests of gene-gene (G x G) interaction using Cordells method yielded significant evidence between single-nucleotide polymorphisms in FGF9 and FGF18, which was confirmed in an independent sample of trios from an international consortium. Conclusion Our results suggest several genes in the FGF/FGFR family may influence risk for isolated, nonsyndromic cleft lip with or without cleft palate through distinct biological mechanisms.

Maternal transmission effects of the PAX genes among cleft case–parent trios from four populations

Isolated cleft lip with or without cleft palate (CL/P) is among the most common human birth defects, with a prevalence of 1 in 700 live births. The paired box (PAX) genes have been suggested as candidate genes for CL/P based largely on mouse models; however, few human studies have focused on this gene family. This study tests for association between markers in four PAX genes and CL/P using a case-parent trio design considering parent-of-origin effects. Trios from four populations (76 from Maryland, 146 from Taiwan, 35 from Singapore, and 40 from Korea) were genotyped for 34 single nucleotide polymorphisms (SNPs) in the PAX3, PAX6, PAX7, and PAX9 genes. We performed the transmission disequilibrium test (TDT) on individual SNPs. Parent-of-origin effects were assessed using the transmission asymmetry test (TAT) and the parent-of-origin likelihood ratio test (PO-LRT). TDT analysis showed one SNP (rs766325) in PAX7 yielding evidence of linkage and association when parent-of-origin was not considered, with an OR(transmission)=1.62 (P=0.003), and five SNPs in PAX6 (including two pairs in near perfect linkage disequilibrium). TAT analysis of all trios revealed two SNPs in PAX7 and four SNPs in PAX3 showing significant excess maternal transmission. For these six SNPs, the maternal OR(transmission) ranged between 1.74 and 2.40, and PO-LRT was also significant (P-values=0.035–0.012). When this analysis was limited to trios with male cases, SNPs in PAX7 showed higher maternal OR(transmission) and greater significance. PAX genes may influence the risk of CL/P through maternal effects, possibly imprinting, which seems to be stronger among male cases.

Association between genes on chromosome 4p16 and non-syndromic oral clefts in four populations

Isolated cleft lip with or without cleft palate and cleft palate are among the most common human birth defects. Several candidate gene studies on MSX1 have shown significant association between markers in MSX1 and risk of oral clefts, and re-sequencing studies have identified multiple mutations in MSX1 in a small minority of cases, which may account for 1–2% of all isolated oral clefts cases. We explored the 2-Mb region around MSX1, using a marker map of 393 single nucleotide polymorphisms (SNPs) in 297 cleft lip, with or without cleft palate, case–parent trios and 84 cleft palate trios from Maryland, Taiwan, Singapore, and Korea. Both individual markers and haplotypes of two to five SNPs showed several regions yielding statistical evidence for linkage and disequilibrium. Two genes (STK32B and EVC) yielded consistent evidence from cleft lip, with or without cleft palate, trios in all four populations. These two genes plus EVC2 also yielded suggestive evidence for linkage and disequilibrium among cleft palate trios. This analysis suggests that several genes, not just MSX1, in this region may influence risk of oral clefts.

Bmp4 was associated with nscl/p in an asian population

Background The Bone Morphogenetic Protein 4 gene (BMP4) is located in chromosome 14q22-q23 which has shown evidence of linkage for isolated nonsyndromic cleft lip with or without cleft palate (NSCL/P) in a genome wide linkage analysis of human multiplex families. BMP4 has been shown to play crucial roles in lip and palatal development in animal models. Several candidate gene association analyses also supported its potential risk for NSCL/P, however, results across these association studies have been inconsistent. The aim of the current study was to test for possible association between markers in and around the BMP4 gene and NSCL/P in Asian and Maryland trios. Methodology/Principal Findings Family Based Association Test was used to test for deviation from Mendelian assortment for 12 SNPs in and around BMP4. Nominal significant evidence of linkage and association was seen for three SNPs (rs10130587, rs2738265 and rs2761887) in 221 Asian trios and for one SNP (rs762642) in 76 Maryland trios. Statistical significance still held for rs10130587 after Bonferroni correction (corrected p = 0.019) among the Asian group. Estimated odds ratio for carrying the apparent high risk allele at this SNP was 1.61 (95%CI = 1.20, 2.18). Conclusions Our results provided further evidence of association between BMP4 and NSCL/P.

Genome wide study of maternal and parent-of-origin effects on the etiology of orofacial clefts.

We performed a genome wide association analysis of maternally‐mediated genetic effects and parent‐of‐origin (POO) effects on risk of orofacial clefting (OC) using over 2,000 case‐parent triads collected through an international cleft consortium. We used log‐linear regression models to test individual SNPs. For SNPs with a P‐value <10−5 for maternal genotypic effects, we also applied a haplotype‐based method, TRIMM, to extract potential information from clusters of correlated SNPs. None of the SNPs were significant at the genome wide level. Our results suggest neither maternal genome nor POO effects play major roles in the etiology of OC in our sample. This finding is consistent with previous genetic studies and recent population‐based cohort studies in Norway and Denmark, which showed no apparent difference between mother‐to‐offspring and father‐to‐offspring recurrence of clefting. We, however, cannot completely rule out maternal genome or POO effects as risk factors because very small effects might not be detectable with our sample size, they may influence risk through interactions with environmental exposures or may act through a more complex network of interacting genes. Thus, the most promising SNPs identified by this study may still be worth further investigation. Published 2012. This article is a U.S. Government work and is in the public domain in the USA.

Genomic, cDNA, and embryonic expression analysis of zebrafish transforming growth factor beta 3 (tgfβ3)

TGFβ3, a member of the transforming growth factor β family, regulates a spectrum of biological processes and is involved in mammalian pulmonary and craniofacial development. Homologs of human TGFβ3 have been identified in several vertebrate species. We sequenced a cDNA clone of zebrafish tgfβ3, consisting of a 271‐bp 5′ untranslated region, a 1,233‐bp open reading frame that encodes a predicted 410 amino acid peptide, and a 527‐bp 3′ untranslated region. Using 5′ rapid amplification of cDNA ends, the transcription start site of this gene was determined to lie an additional 29 nucleotides upstream. The gene is composed of seven exons and maps to a segment of linkage group 17 that is syntenic to the human TGFβ3 locus on chromosome 14q24. One stimulating protein 1 (Sp1) and two (TATA binding protein) (TBP) transcription factor binding sites were identified in the putative promoter segment upstream of the transcription start site. Comparative alignment analysis revealed a high degree of tgfβ3 nucleotide and amino acid identity between zebrafish and other species, including complete conservation of the cysteine knot structure that facilitates protein–protein interaction. Also, 9 of 10 amino acid residues critical for ligand/receptor binding in human TGFβ3 are conserved in zebrafish, suggesting a high degree of functional conservation even in lower vertebrates. Zebrafish tgfβ3 transcripts were first detected in the notochord (10 somite to high‐pec stage), followed by expression in the developing pharyngeal arch and neurocranial cartilage (18 somite to protruding mouth stage), lens and heart (21 somite to protruding mouth stage), and pectoral fins (prim‐25 to protruding mouth stage). The strong expression in the pectoral fins, not reported in the orthologous mammalian forelimb, suggests a modified or novel function of tgfβ3 during early fish development. Developmental Dynamics 232:1021–1030, 2005.

tgfβ3 regulation of chondrogenesis and osteogenesis in zebrafish is mediated through formation and survival of a subpopulation of the cranial neural crest

Zebrafish tgfbeta3 is strongly expressed in a subpopulation of the migrating neural crest cells, developing pharyngeal arches and neurocranial cartilages. To study the regulatory role of tgfbeta3 in head skeletal formation, we knocked down tgfbeta3 in zebrafish and found impaired craniofacial chondrogenesis, evident by malformations in selected neurocranial and pharyngeal arch cartilages. Over-expressing tgfbeta3 in embryos resulted in smaller craniofacial cartilages without any gross malformations. These defects suggest that tgfbeta3 is required for normal chondrogenesis. To address the cellular mechanisms that lead to the observed malformations, we analyzed cranial neural crest development in morphant and tgfbeta3 over-expressing fish. We observed reduced pre-migratory and migratory cranial neural crest, the precursors of the neurocranial cartilage and pharyngeal arches, in tgfbeta3 knockdown embryos. In contrast, only the migratory neural crest was reduced in embryos over-expressing tgfbeta3. This raised the possibility that the reduced number of cranial neural crest cells is a result of increased apoptosis. Consistent with this, markedly elevated TUNEL staining in the midbrain and hindbrain, and developing pharyngeal arch region was observed in morphants, while tgfbeta3 over-expressing embryos showed marginally increased apoptosis in the developing pharyngeal arch region. We propose that both Tgfbeta3 suppression and over-expression result in reduced chondrocyte and osteocyte formation, but to different degrees and through different mechanisms. In Tgfbeta3 suppressed embryos, this is due to impaired formation and survival of a subpopulation of cranial neural crest cells through markedly increased apoptosis in regions containing the cranial neural crest cells, while in Tgfbeta3 over-expressing embryos, the milder phenotype is also due to a slightly elevated apoptosis in these regions. Therefore, proper cranial neural crest formation and survival, and ultimately craniofacial chondrogenesis and osteogenesis, are dependent on tight regulation of Tgfbeta3 protein levels in zebrafish.