Evidence that TGFA influences risk to cleft lip with/without cleft palate through unconventional genetic mechanisms

IRF6 and FGF1 polymorphisms in non-syndromic cleft lip with or without cleft palate in the Polish population

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is the most common developmental defect that significantly affects the morphology and function of the stomatognathic system in children. The etiology of these birth defects is multifactorial, and single nucleotide polymorphisms (SNPs) in IRF6 and FGF1 have been associated with NSCL/P. This study aimed to evaluate whether SNPs in IRF6, namely rs2013162, rs642961, rs2235373, and rs34010 in FGF1, are associated with NSCL/P occurrence in the Polish population. The study included 627 participants: 209 children with NSCL/P and 418 healthy controls. DNA was isolated from saliva in the study group and from umbilical cord blood in controls. Genotyping of polymorphisms was performed using quantitative PCR. There was no statistically significant association of IRF6 gene variants with NSCL/P occurrence, although for rs2013162, AA genotype, odds ratio (OR) = 1.16 and for AC genotype, OR = 0.83; for rs642961, AA genotype, OR = 0.84 and for AG genotype, OR = 1.41; and for rs2235373, AA genotype, OR = 0.79 and for AG, OR = 0.85. In the instance of rs34010 polymorphism in FGF1, the presence of the AA genotype was statistically significant in reducing the risk of NSCL/P (OR = 0.31, p = 0.001). Genetic variation in FGF1 is an important risk marker of NSCL/P in the Polish population, which cannot be stated for the polymorphisms in the IRF6 gene.

Genome-Wide Scan for Parent-of-Origin Effects in a sub-Saharan African Cohort With Nonsyndromic Cleft Lip and/or Cleft Palate (CL/P)

OBJECTIVE

Nonsyndromic cleft lip and/or cleft palate (NSCL/P) have multifactorial etiology where genetic factors, gene-environment interactions, stochastic factors, gene-gene interactions, and parent-of-origin effects (POEs) play cardinal roles. POEs arise when the parental origin of alleles differentially impacts the phenotype of the offspring. The aim of this study was to identify POEs that can increase risk for NSCL/P in humans using a genome-wide dataset.

METHODS

The samples (174 case-parent trios from Ghana, Ethiopia, and Nigeria) included in this study were from the African only genome wide association studies (GWAS) that was published in 2019. Genotyping of individual DNA using over 2 million multiethnic and African ancestry-specific single-nucleotide polymorphisms from the Illumina Multi-Ethnic Genotyping Array v2 15070954 A2 (genome build GRCh37/hg19) was done at the Center for Inherited Diseases Research. After quality control checks, PLINK was employed to carry out POE analysis employing the pooled subphenotypes of NSCL/P.

RESULTS

We observed possible hints of POEs at a cluster of genes at a 1 mega base pair window at the major histocompatibility complex class 1 locus on chromosome 6, as well as at other loci encompassing candidate genes such as ASB18, ANKEF1, AGAP1, GABRD, HHAT, CCT7, DNMT3A, EPHA7, FOXO3, lncRNAs, microRNA, antisense RNAs, ZNRD1, ZFAT, and ZBTB16.

CONCLUSION

Findings from our study suggest that some loci may increase the risk for NSCL/P through POEs. Additional studies are required to confirm these suggestive loci in NSCL/P etiology.

Family based and case-control designs reveal an association of TFAP2A in nonsyndromic cleft lip only among Vietnamese population

Aims

Dozens of causative genes and their mechanisms of nonsyndromic cleft lip with or without cleft palate (NSCL/P) were revealed through genome‐wide association and linkage studies. Results were, however, not always replicated in different populations or methodologies. This study used case–control and family based approaches to investigate the etiology of NSCL/P and its two subtypes: nonsyndromic cleft lip only (NSCLO) and nonsyndromic cleft lip and palate (NSCLP) among the Vietnamese population.

Methods

Two hundred and seventeen NSCL/P case‐parent trios (one affected child and two parents), including 105 NSCLO and 112 NSCLP were involved for a family based design; and 273 ethnic and region‐matched healthy controls with no cleft history in their families were recruited for a case–control design. Three SNPs consisting of TFAP2A (rs1675414 and rs303048) and 8q24 (rs987525) were genotyped using the TaqMan SNP genotyping assay.

Results

TFAP2A rs1675414 was associated with NSCLO, replicated by both case‐control and family based tests. Other SNPs yielded no evidence of susceptibility to NSCL/P or two subtypes.

Conclusion

The current investigation suggests an intriguing role of TFAP2A in the etiology of NSCLO among the Vietnamese population.

This study used case‐control and family‐based approaches to investigate the etiology of NSCL/P and its two subtypes: nonsyndromic cleft lip only (NSCLO), nonsyndromic cleft lip and palate (NSCLP) among Vietnamese population. TFAP2A rs1675414 was associated with NSCLO, replicated by both case‐control and family‐based tests.

Paternal genetic variants and risk of obstructive heart defects: A parent-of-origin approach

Previous research on risk factors for obstructive heart defects (OHDs) focused on maternal and infant genetic variants, prenatal environmental exposures, and their potential interaction effects. Less is known about the role of paternal genetic variants or environmental exposures and risk of OHDs. We examined parent-of-origin effects in transmission of alleles in the folate, homocysteine, or transsulfuration pathway genes on OHD occurrence in offspring. We used data on 569 families of liveborn infants with OHDs born between October 1997 and August 2008 from the National Birth Defects Prevention Study to conduct a family-based case-only study. Maternal, paternal, and infant DNA were genotyped using an Illumina Golden Gate custom single nucleotide polymorphism (SNP) panel. Relative risks (RR), 95% confidence interval (CI), and likelihood ratio tests from log-linear models were used to estimate the parent-of-origin effect of 877 SNPs in 60 candidate genes in the folate, homocysteine, and transsulfuration pathways on the risk of OHDs. Bonferroni correction was applied for multiple testing. We identified 3 SNPs in the transsulfuration pathway and 1 SNP in the folate pathway that were statistically significant after Bonferroni correction. Among infants who inherited paternally-derived copies of the G allele for rs6812588 in the RFC1 gene, the G allele for rs1762430 in the MGMT gene, and the A allele for rs9296695 and rs4712023 in the GSTA3 gene, RRs for OHD were 0.11 (95% CI: 0.04, 0.29, P = 9.16x10^-7), 0.30 (95% CI: 0.17, 0.53, P = 9.80x10^-6), 0.34 (95% CI: 0.20, 0.57, P = 2.28x10^-5), and 0.34 (95% CI: 0.20, 0.58, P = 3.77x10^-5), respectively, compared to infants who inherited maternally-derived copies of the same alleles. We observed statistically significant decreased risk of OHDs among infants who inherited paternal gene variants involved in folate and transsulfuration pathways.

Obstructive heart defects are birth defects that cause obstruction to the blood flow of the developing heart. Common OHDs include coarctation of the aorta, aortic stenosis and pulmonary stenosis. While there is a fair amount of literature indicating an association between maternal genetic variants and OHDs, less is known about the role of paternal genetic variants in the etiology of OHDs. We used a genotype clustering algorithm, SNPMClust, that was developed in-house at the Arkansas Center for Birth Defects Research and Prevention to study the role of paternal genetic variants in the folate, homocysteine and transsulfuration pathways. Maternal, paternal, and infant DNA specimens were collected from participants of the National Birth Defects Prevention Study, a large population-based case-control study in the United States, and were genotyped using an Illumina Golden Gate custom single nucleotide polymorphism (SNP) panel. We identified 4 SNPs in the folate and transsulfuration pathways, rs6812588, rs1762430, rs9296695, and rs4712023, that were associated with a statistically significant decreased risk of OHDs for infants who inherited a paternally-derived copy of the variant allele compared to infants who inherited a maternal copy of the variant allele. In conclusion, we observed a significantly decreased risk and less epigenetic influence of paternal genetic variants on OHDs compared to maternally-derived variants.

Gene-environment interaction in molar-incisor hypomineralization

Molar incisor hypomineralization (MIH) is an enamel condition characterized by lesions ranging in color from white to brown which present rapid caries progression, and mainly affects permanent first molars and incisors. These enamel defects usually occur when there are disturbances during the mineralization or maturation stage of amelogenesis. Both genetic and environmental factors have been suggested to play roles in MIH’s development, but no conclusive risk factors have shown the source of the disease. During head and neck development, the interferon regulatory factor 6 (IRF6) gene is involved in the structure formation of the oral and maxillofacial regions, and the transforming growth factor alpha (TGFA) is an essential cell regulator, acting during proliferation, differentiation, migration and apoptosis. In this present study, it was hypothesized that these genes interact and contribute to predisposition of MIH. Environmental factors affecting children that were 3 years of age or older were also hypothesized to play a role in the disease etiology. Those factors included respiratory issues, malnutrition, food intolerance, infection of any sort and medication intake. A total of 1,065 salivary samples from four different cohorts were obtained, and DNA was extracted from each sample and genotyped for nine different single nucleotide polymorphisms. Association tests and logistic regression implemented in PLINK were used for analyses. A potential interaction between TGFA rs930655 with all markers tested in the cohort from Turkey was identified. These interactions were not identified in the remaining cohorts. Associations (p<0.05) between the use of medication after three years of age and MIH were also found, suggesting that conditions acquired at the age children start to socialize might contribute to the development of MIH.

Environmental mechanisms of orofacial clefts

Orofacial clefts (OFCs) are among the most common birth defects and impart a significant burden on afflicted individuals and their families. It is increasingly understood that many nonsyndromic OFCs are a consequence of extrinsic factors, genetic susceptibilities, and interactions of the two. Therefore, understanding the environmental mechanisms of OFCs is important in the prevention of future cases. This review examines the molecular mechanisms associated with environmental factors that either protect against or increase the risk of OFCs. We focus on essential metabolic pathways, environmental signaling mechanisms, detoxification pathways, behavioral risk factors, and biological hazards that may disrupt orofacial development.

[Exploring parent-of-origin effects for non-syndromic cleft lip with or without cleft palate on PTCH1, PTCH2, SHH, SMO genes in Chinese case-parent trios]

OBJECTIVE

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect, affecting 1.4 per 1 000 live births, and multiple genetic and environmental risk factors influencing its risk. All the known genetic risk factors accounted for a small proportion of the heritability. Several authors have suggested parent-of-origin effects (PoO) may play an important role in the etiology of this complex and heterogeneous malformation. To clarify the genetic association between PTCH1, PTCH2, SHH and SMO in hedgehog (HH) pathway and NSCL/P, as well as testing for potential PoO effects in Chinese case-parent trios.

METHODS

We tested for transmission disequilibrium tests (TDT) and PoO effects using 83 common single nucleotide polymorphic (SNP) markers of HH pathway genes from 806 NSCL/P case-parent trios. These trios were drawn from an international consortium established for a genome-wide association studies (GWAS) of non-syndromic oral clefts of multiple ethnicities. DNA samples were collected from each trio. Single marker and haplotype based analysis were performed both in TDT tests and PoO effects. SNPs were excluded if they (ⅰ) had a call rate of < 95%, (ⅱ) had a minor allele frequency (MAF) of < 0.05, (ⅲ) had Mendelian errors over all trios of >5%, (ⅳ) had a genotype distribution in the parents that deviated from the Hardy-Weinberg equilibrium (HWE) (P < 0.000 1). The process was done using Plink (version 1.07, <a href="http://pngu.mgh.harvard.edu/~purcell/plink/data.shtml" target="_blank">http://pngu.mgh.harvard.edu/~purcell/plink/data.shtml</a>). TDT test was performed in Plink v1.07. A log-linear model was used to explore PoO effects using Haplin v6.2.1 as implemented in R package v3.4.2. Significance level was assessed using the Bonferroni correction.

RESULTS

A total of 18 SNPs were dropped due to low MAF, thus leaving 65 SNPs available for the analysis. Thus the Bonferroni threshold was 7.7×10^-4 (0.05/65). Nominal significant association with NSCL/P was found at a SNP (rs4448343 in PTCH1, P=0.023) and six haplotypes (rs10512249-rs4448343, rs1461208-rs7786445, rs10512249-rs4448343, rs16909865-rs10512249-rs4448343, rs1461208-rs7786445-rs12698335, and rs288756-rs288758-rs1151790, P < 0.05). A total of six haplotypes (rs288765-rs1233563, rs12537550-rs11765352, rs872723-rs288765-rs1233563, rs288765-rs1233563-rs288756, rs6459952-rs12537550-rs11765352, and rs12537550-rs11765352-rs6971211) showed PoO effect (P < 0.05). None of the results remained significant after the Bonferroni correction (P>7.7×10^-4).

CONCLUSION

Neither significant association between SNPs within HH pathway and the risk of NSCL/P nor PoO effects was seen in this study.

Evidence for germline non-genetic inheritance of human phenotypes and diseases

It is becoming increasingly apparent that certain phenotypes are inherited across generations independent of the information contained in the DNA sequence, by factors in germ cells that remain largely uncharacterized. As evidence for germline non-genetic inheritance of phenotypes and diseases continues to grow in model organisms, there are fewer reports of this phenomenon in humans, due to a variety of complications in evaluating this mechanism of inheritance in humans. This review summarizes the evidence for germline-based non-genetic inheritance in humans, as well as the significant challenges and important caveats that must be considered when evaluating this process in human populations. Most reports of this process evaluate the association of a lifetime exposure in ancestors with changes in DNA methylation or small RNA expression in germ cells, as well as the association between ancestral experiences and the inheritance of a phenotype in descendants, down to great-grandchildren in some cases. Collectively, these studies provide evidence that phenotypes can be inherited in a DNA-independent manner; the extent to which this process contributes to disease development, as well as the cellular and molecular regulation of this process, remain largely undefined.

Association between methylenetetrahydrofolate reductase polymorphisms and non-syndromic cleft lip with or without palate susceptibility: an updated systematic review and meta-analysis

Methylenetetrahydrofolate reductase (MTHFR) polymorphisms are thought to be involved in the development of cleft lip with or without cleft palate (NSCL/P), but published results are contradictory. We therefore designed an updated meta-analysis to pool eligible studies and to evaluate further the possible relations between MTHFR polymorphisms (c.677C>T and c.1298A>C) and susceptibility to NSCL/P. A comprehensive search based on PubMed, Medline, Web of Science, and Embase databases was made up to February 2018. Twenty-three case-control and 10 case-parent trio studies (including 1149 cases and 1161 controls) were retrieved. Odds ratio (OR) with 95% CI were used to estimate the pooled strength of association under different genetic models. The Q test and I² test were used to estimate heterogeneity among studies, the quality of which was assessed using the Newcastle-Ottawa scale. In the MTHFR c.677C>T polymorphism group, there were significant overall results for the recessive (OR 1.231, 95%CI 1.092 to 1.387) and homozygote (OR 1.252, 95%CI 1.078 to 1.456) models. Subgroup analysis by subjects and ethnicity identified only associations in European mothers for the recessive model and the homozygote model. For the c.1298A>C group, there were no significant results for either European or Asian patients for all genetic models. The MTHFR c.677C>T polymorphism might increase susceptibility to NSCL/P in European mothers, but was negatively associated in Asian patients, and the MTHFR c.1298A>C polymorphism is not involved in the development of NSCL/P in either European or Asian patients.

A parent-of-origin analysis of paternal genetic variants and increased risk of conotruncal heart defects

The association between conotruncal heart defects (CTHDs) and maternal genetic and environmental exposures is well studied. However, little is known about paternal genetic or environmental exposures and risk of CTHDs. We assessed the effect of paternal genetic variants in the folate, homocysteine, and transsulfuration pathways on risk of CTHDs in offspring. We utilized National Birth Defects Prevention Study data to conduct a family-based case only study using 616 live-born infants with CTHDs, born October 1997-August 2008. Maternal, paternal and infant DNA was genotyped using an Illumina® Golden Gate custom single nucleotide polymorphism (SNP) panel. Relative risks (RR) and 95% confidence intervals (CI) from log-linear models determined parent of origin effects for 921 SNPs in 60 candidate genes involved in the folate, homocysteine, and transsulfuration pathways on risk of CTHDs. The risk of CTHD among children who inherited a paternally derived copy of the A allele on GLRX (rs17085159) or the T allele of GLRX (rs12109442) was 0.23 (95%CI: 0.12, 0.42; p = 1.09 × 10^-6 ) and 0.27 (95%CI: 0.14, 0.50; p = 2.06 × 10^-5 ) times the risk among children who inherited a maternal copy of the same allele. The paternally inherited copy of the GSR (rs7818511) A allele had a 0.31 (95%CI: 0.18, 0.53; p = 9.94 × 10^-6 ] risk of CTHD compared to children with the maternal copy of the same allele. The risk of CTHD is less influenced by variants in paternal genes involved in the folate, homocysteine, or transsulfuration pathways than variants in maternal genes in those pathways.

Gene-gene interaction between MSX1 and TP63 in Asian case-parent trios with nonsyndromic cleft lip with or without cleft palate

BACKGROUND

Small ubiquitin-like modification, also known as sumoylation, is a crucial post-translational regulatory mechanisms involved in development of the lip and palate. Recent studies reported two sumoylation target genes, MSX1 and TP63, to have achieved genome-wide level significance in tests of association with nonsyndromic clefts. Here, we performed a candidate gene analysis considering gene-gene and gene-environment interaction for SUMO1, MSX1, and TP63 to further explore the etiology of nonsyndromic cleft lip with or without cleft palate (NSCL/P).

METHODS

A total of 130 single-nucleotide polymorphisms (SNPs) in or near SUMO1, MSX1, and TP63 was analyzed among 1,038 Asian NSCL/P trios ascertained through an international consortium. Conditional logistic regression models were used to explore gene-gene (G × G) and gene-environment (G × E) interaction involving maternal environmental tobacco smoke and multivitamin supplementation. Bonferroni correction was used for G × E analysis and permutation tests were used for G × G analysis.

RESULTS

While transmission disequilibrium tests and gene-environment interaction analysis showed no significant results, we did find signals of gene-gene interaction between SNPs near MSX1 and TP63. Three pairwise interactions yielded significant p values in permutation tests (rs884690 and rs9290890 with p = 9.34 × 10^-5 and empirical p = 1.00 × 10^-4 , rs1022136 and rs4687098 with p = 2.41 × 10^-4 and empirical p = 2.95 × 10^-4 , rs6819546 and rs9681004 with p = 5.15 × 10^-4 and empirical p = 3.02 × 10^-4 ).

CONCLUSION

Gene-gene interaction between MSX1 and TP63 may influence the risk of NSCL/P in Asian populations. Our study provided additional understanding of the genetic etiology of NSCL/P and underlined the importance of considering gene-gene interaction in the etiology of this common craniofacial malformation.

Interaction Effect of RsaI and BamHI Polymorphisms of TGFα, BMP2 and BMP4 on the Occurrence of Non-Syndromic Cleft Lip and Palate in Iranian Patients

BACKGROUND

Orofacial cleft is the most common congenital defect of the maxillofacial region. Its non-syndromic type is multi-factorial, and several genes are involved in its occurrence. This study aimed to assess the interaction effect of Rsal and BamHI polymorphisms of Transforming Growth Factor-alpha (TGFα) gene and Bone Morphogenetic Protein-2 (BMP2) and BMP4 variants on the occurrence of Non-Syndromic Cleft Lip and Palate (NSCLP) in the Iranian population.

METHODS

This case-control study was conducted on 120 children with NSCLP and 215 healthy children. Genotyping of the TGFA/BamHI (rs11466297), TGFA/RsaI (rs37322-48), BMP4 (rs17563) and BMP2 (rs235768) was performed by Polymerase Chain Reaction (PCR) and Restriction Fragment Length Polymorphism (RFLP) methods. Logistic regression was applied to determine the effective factors and the interaction effect of different variants on the occurrence of NSCLP.

RESULTS

Gender of patients had no significant association with the occurrence of NSCLP (p=0.335). Multiple logistic regression showed that the interaction effect of the aforementioned polymorphisms on the occurrence of NSCLP was not statistically significant (p=1.000).

CONCLUSION

Although the individual effect of each of the BMP4, BMP2, RsaI and BamHI variants on the occurrence of NSCLP in the Iranian population has been previously confirmed, their interaction does not play a role in this respect.

Gene-gene interaction of single nucleotide polymorphisms in 16p13.3 may contribute to the risk of non-syndromic cleft lip with or without cleft palate in Chinese case-parent trios

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect with a complex and heterogeneous etiology. A recent genome-wide association study (GWAS) among Chinese populations has identified a new region at 16p13.3 as being associated with NSCL/P, which requires further replication. Here, we attempted to replicate and further clarify the genetic association between this region and NSCL/P, as well as testing for potential gene-gene (G × G) and gene-environment (G × E) interactions. We conducted transmission disequilibrium tests on 69 single nucleotide polymorphisms (SNPs) mapping to 16p13.3 among 806 Chinese case-parent trios ascertained through an international consortium where a GWAS of oral clefts was conducted. G × G, as well as G × E interactions involving maternal environmental tobacco smoke (ETS) and multivitamin supplementation, were explored using conditional logistic regression model. We applied Cordell's method as implemented in the R package TRIO to test for possible interactions. While no SNPs showed evidence of linkage and association with NSCL/P after Bonferroni correction, we found signals of G × G interactions between SNPs in 16p13.3. Nine pairs of SNP-SNP interactions attained significance after Bonferroni correction, among which the most significant interaction was found between rs2072346 (ADCY9) and rs11646137 (intergenic region, P = 7.2 × 10^-5 ). Linkage disequilibrium (LD) analysis revealed only low level of LD between these SNPs. This study failed to confirm the significant association between SNPs within 16p13.3 and the risk of NSCL/P, but underlined the importance of taking into account potential G × G interactions for the genetic association analysis of NSCL/P.

Associations Between TGFA/TGFB3/MSX1 Gene Polymorphisms and Congenital Non-Syndromic Hearing Impairment in a Chinese Population

BACKGROUND The aim of this study was to investigate whether the TGFA/TGFB3/MSX1 gene polymorphisms and haplotypes lead to individual differences between congenital non-syndromic hearing impairment (NSHI) patients and normal people in a Chinese population and to analyze the risk factors for NSHI. MATERIAL AND METHODS Between December 2010 and September 2014, 343 congenital NSHI patients were recruited as cases, and 272 healthy subjects were recruited as controls. Denaturing high-performance liquid chromatography (DHPLC) was used to identify genotypes, SHEsis software was used to conduct gene linkage disequilibrium and haplotype analyses, and regression analysis was performed to identify risk factors for congenital NSHI. RESULTS The distribution of genotype frequencies and allele frequencies of TGFA rs3771494, TGFB3 rs3917201 and rs2268626, and MSX1 rs3821949 and rs62636562 were significantly different between the case and the control groups (all P<0.05). TGFA/TGFB3/MSX1 gene rs3771494, rs1058213, rs3917201, rs2268626, rs3821949, and rs62636562 haplotype analysis showed that haplotype CCGTAC and TTACGT might be protective factors (both P<0.001), while TTGCGC might be a risk factor for the normal population (P<0.001). The other risk factors include paternal smoking, advanced maternal age, maternal sickness history, maternal contact with pesticides or similar drugs, maternal abortion history, maternal medication history, maternal passive smoking history during pregnancy, rs3771494 CT, rs2268626 CC and TC, and rs3821949 GG and AG genotypes were risk factors (all P<0.05), while maternal vitamin supplements during pregnancy, rs3917201 GA, rs62636562 TT and CT genotypes were protective factors for congenital NSHI (all P<0.05). CONCLUSIONS rs3771494, rs3917201, rs2268626, rs3821949 and rs62636562 might be associated with congenital NSHI.

Genomic expression in non syndromic cleft lip and palate patients: A review

Cleft lip and palate are common congenital anomalies with significant medical, psychological, social, and economic ramifications, affecting one in seven hundred live births. Genetic causes of non syndromic cleft lip and/or palate (NSCLP) include chromosomal rearrangements, genetic susceptibility to teratogenic exposures, and complex genetic contributions of multiple genes. Development of the orofacial clefts in an individual will depend on the interaction of several moderately effecting genes with environmental factors. Several candidate genes have been genotyped in different population types, using case parent trio or case control design; also genes have been sequenced and SNPs have been reported. Quantitative and molecular analysis have shown linkage and association studies to be more relevant. Recent literature search shows genome wide association studies using microarray. The aim of this paper was to review the approaches to identify genes associated with NSCLP and to analyze their differential expressions. Although no major gene has been confirmed, a lot of research is ongoing to provide an understanding of the pathophysiology of the orofacial clefts.

Genetic Effect of Transforming Growth Factor Alpha Gene Variants on the Risk of Nonsyndromic Cleft Lip with or without Palate in Korean Populations

Objective

To identify the contribution of TGFA gene variants to the risk of nonsyndromic cleft lip with or without palate (NS-CL±P).

Design

The samples were from 142 Korean NS-CL±P families and 119 control parents having nonaffected children. Minor allele frequency, heterozygosity, and χ2 test for Hardy-Weinberg equilibrium were calculated for each of 10 selected single-nucleotide polymorphisms (SNPs). Ten SNPs were used to examine the association of case-parent trios with the transmission disequilibrium test (TDT) and conditional logistic regression models (CLRMs). Both allelic and genotypic TDTs for individual SNPs and sliding windows of haplotypes consisting of two to five SNPs were tested using family- and haplotype-based association test programs. Genotypic odd ratios (GORs) were obtained from CLRMs using STATA software. The parent-of-origin effect was evaluated for 10 SNPs, and a comparison between 218 case parents and 119 control parents was performed to investigate paternal and maternal ORs.

Results

Family-based TDT and haplotype analysis exhibited no statistical significance, but a relatively meaningful association was shown with rs3771497 (all P < .05; two SNPs, rs3771497 and rs3755377; five SNPs, rs3771497, rs3755377, rs3771485, rs11466212, and rs3771475). G/G homozygotes at rs3771497 have a significant decreased risk of NS-CL±P (GOR = 0.30, P < .01). No SNPs showed parent-of-origin effects. However, in the comparison between case parents and control parents, a single-marker analysis of maternal line showed a significant association with NS-CL±P in rs3771497 ( P < .001, recessive model).

Conclusion

The association of the TGFA gene with NS-CL±P in Korean populations was not clearly found. However, the etiologic effect of the TGFA gene on NS-CL±P patients should be investigated in terms of maternal genotype influence.

Role of angiogenesis-related genes in cleft lip/palate: review of the literature

OBJECTIVES

Cleft lip and cleft palate (CLP) are the most common congenital craniofacial anomalies. They have a multifactorial etiology and result from an incomplete fusion of the facial buds. Two main mechanisms, acting alone or interacting with each other, were evidenced in this fusion defect responsible for CLP: defective tissue development and/or defective apoptosis in normal or defective tissues. The objective of this work was to study the implication and role of angiogenesis-related genes in the etiology of CL/P.

METHODS

Our methodological approach included a systematic and thorough analysis of the genes involved in CL/P (syndromic and non-syndromic forms) including previously identified genes but also genes that could potentially be angiogenesis-related (OMIM, Pub Med).We studied the interactions of these different genes and their relationships with potential environmental factors.

RESULTS

TGFβ, FGA, PDGFc, PDGFRa, FGF, FGFR1, FGFR2 growth factors as well as MMP and TIMP2 proteolytic enzymes are involved in the genesis of CLP (P>L). Furthermore, 18 genes involved in CLP also interact with angiogenesis-related genes.

DISCUSSION

Even if the main angiogenesis-related genes involved in CLP formation are genes participating in several biological activities and their implication might not be always related to angiogenesis defects, they nevertheless remain an undeniably important research pathway. Furthermore, their interactions with environmental factors make them good candidates in the field of CLP prevention.

Association of JARID2 polymorphisms with non-syndromic orofacial clefts in northern Chinese Han population

OBJECTIVES

Non-syndromic orofacial clefts (NSOC) are the most common human craniofacial malformation in all worldwide populations. Recently, the jumoji AT-rich interaction domain 2 (JARID2) had been reported to be a novel candidate gene for non-syndromic cleft lip with or without cleft palate (CL/P). The SNPs rs2076056, rs2237138 and rs2299043 in JARID2 were highly significant in Italian families.

MATERIAL AND METHODS

In the current research, a case-control study was conducted to examine the association between these three SNPs and NSOC in a northern Chinese Han population. Genotyping of the three SNPs were performed using SNaPshot minisequencing technique.

RESULTS

Distribution of rs2237138 genotypes in CL/P group was different from those in the control group (P = 0.04), but significant results did not persist after Benjamini and Hochberg false discovery rate (FDR) correction for multiple tests. Further logistic regression analysis showed that rs2237138 GG genotypes were associated with decreased CL/P susceptibility (OR = 0.34, 95% CI = 0.13-0.84), compared with the AA wild-type homozygote. For the haplotype CGT, a statistically difference was identified between the CL/P group and controls (P = 0.04). And carriers of GAT haplotype were considered to be less frequent among cleft palate only group as compared to controls (P = 0.02). However, both of the haplotypes association did not remain statistically significant after Benjamini and Hochberg FDR correction.

CONCLUSION

We got a weak association between these polymorphisms and NSOC in both single-marker and haplotype analyses. Our data further strengthen the conclusion that JARID2 polymorphisms are associated with NSOC susceptibility.

Review of statistical methodologies for the detection of parent-of-origin effects in family trio genome-wide association data with binary disease traits

The detection of parent-of-origin effects aims to identify whether the functionality of alleles, and in turn associated phenotypic traits, depends on the parental origin of the alleles. Different parent-of-origin effects have been identified through a variety of mechanisms and a number of statistical methodologies for their detection have been proposed, in particular for genome-wide association studies (GWAS). GWAS have had limited success in explaining the heritability of many complex disorders and traits, but successful identification of parent-of-origin effects using trio (mother, father and offspring) GWAS may help shed light on this missing heritability. However, it is important to choose the most appropriate parent-of-origin test or methodology, given knowledge of the phenotype, amount of available data and the type of parent-of-origin effect(s) being considered. This review brings together the parent-of-origin detection methodologies available, comparing them in terms of power and type I error for a number of different simulated data scenarios, and finally offering guidance as to the most appropriate choice for the different scenarios.

Investigation of Parental Transmission of RUNX2 Single Nucleotide Polymorphism and Its Association with Nonsyndromic Cleft Lip with or Without Palate

Objective

To investigate the association and parental transmission of RUNX2 single nucleotide polymorphisms (SNPs) with risk of nonsyndromic cleft lip with or without cleft palate (NS-CL±P).

Design

Four RUNX2 SNPs in 142 Korean NS-CL±P families (nine cleft lip, 26 cleft lip and alveolus, and 107 cleft lip and palate; 76 trios and 66 dyads) were genotyped. The minor allele frequency, heterozygosity, and chi-square test for Hardy-Weinberg equilibrium at each SNP were computed between parents. Pairwise linkage disequilibrium was computed as D′ and r2 for all SNPs. Both allelic and genotypic transmission disequilibrium tests (TDTs) were performed for individual SNPs using a family-based association test program. Sliding windows of haplotypes consisting of two to four SNPs were tested using a haplotype-based association test program. Genotypic odds ratios (GORs) were calculated from conditional logistic regression models. Parent-of-origin effects were assessed using transmission asymmetry test and parent-of-origin likelihood ratio test.

Results

The family-based TDT showed significant evidence of linkage and association at rs1934328 ( P = .001). In the haplotype analysis, two, three, and four haplotypes containing rs1934328 revealed significant associations ( P = .0017, P = .0022, and P = .0020, respectively). The genotypes A/T and T/T at rs1934328 were significantly associated with NS-CL±P compared with the genotype A/A (GOR = 2.75, 95% confidence interval [CI] = 1.39–5.45, P =0.0019 in the dominant model; GOR = 5.38, 95% CI = 1.34–21.68, P = .0046 in the additive model). However, no parent-of origin effect was observed.

Conclusion

These findings suggest possible involvement of RUNX2-rs 194328 in the etiology of NS-CL±P in Korean cleft-parent trios without excess parental transmission.

TGF-α gene variations and increased susceptibility of gastric cancer in an Eastern Chinese Han population

Transforming growth factor-alpha (TGF-α) correlates with deep invasion, advanced stage and poor prognosis in gastric cancer. Genetic variants in the 3' untranslated region (UTR) of TGF-α gene may influence the stability and post-transcriptional regulation of mRNA and contribute to gastric cancer predisposition. To test this hypothesis, we genotyped five polymorphisms in 3'UTR (rs3771527, rs503314, rs473698, rs3732253 and rs538118) and one in 3' near region (rs11466306) of the TGF-α gene by polymerase chain reaction-ligation detection reaction methods (PCR-LDR). We found that GA/AA genotype of rs11466306 in the 3' near gene could increase the risk of overall gastric cancer (adjusted OR = 1.499, 95%CI: 1.101-2.041), compared to the wild homozygous GG genotype. Meanwhile, the risk effect was more obvious in the intestinal gastric cancer and gastric noncardia cancer (adjusted OR = 1.682, 95%CI: 1.188-2.380; adjusted OR = 1.495, 95%CI: 1.072-2.086, respectively), but not for the diffuse type and gastric cardia cancer (p > 0.05). CT/TT genotype for rs3732253 in the 3' UTR was associated with increased risk of intestinal gastric cancer (adjusted OR = 1.464, 95%CI: 1.036-2.069), compared to their wild homozygous genotypes. These findings indicate that potentially functional TGF-α gene variant may contribute to the risk of intestinal gastric cancer and/or gastric noncardia cancer and could be used as molecular markers for detecting intestinal gastric cancer and/or gastric noncardia cancer in Chinese Han population.

TGFA and IRF6 contribute to the risk of nonsyndromic cleft lip with or without cleft palate in northeast China

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) are common birth defects with a complex etiology. Multiple interacting loci and possible environmental factors influence the risk of NSCL/P. 12 single nucleotide polymorphisms (SNPs) in 7 candidate genes were tested using an allele-specific primer extension for case-control and case-parent analyses in northeast China in 236 unrelated patients, 185 mothers and 154 fathers, including 128 complete trios, and 400 control individuals. TGFA and IRF6 genes showed a significant association with NSCL/P. In IRF6, statistical evidence of an association between rs2235371 (p = 0.003), rs2013162 (p<0.0001) and NSCL/P was observed in case-control analyses. Family based association tests (FBATs) showed over-transmission of the C allele at the rs2235371 polymorphism (p = 0.007). In TGFA, associations between rs3771494, rs3771523 (G3822A), rs11466285 (T3851C) and NSCL/P were observed in case-control and FBAT analyses. Associations between other genes (BCL3, TGFB3, MTHFR, PVRL1 and SUMO1) and NSCL/P were not detected.

Maternal transmission effect of a PDGF-C SNP on nonsyndromic cleft lip with or without palate from a Chinese population

Cleft lip with or without palate (CL/P) is a common congenital anomaly with a high birth prevalence in China. Based on a previous linkage signal of nonsyndromic CL/P (NSCL/P) on the chromosomal region 4q31–q32 from the Chinese populations, we screened the 4q31–q32 region for susceptibility genes in 214 trios of Han Chinese. PDGF-C, an important developmental factor, resides in the region and has been implicated in NSCL/P. However, in our family-based association test (transmission disequilibrium test; TDT), we could not conclude an association between PDGF-C and NSCL/P as previously suggested. Instead, we found strong evidence for parent-of-origin effect at a PDGF-C SNP, rs17035464, by a likelihood ratio test (unadjusted p-value = 0.0018; I_m = 2.46). The location of rs17035464 is 13 kb downstream of a previously reported, NSCL/P-associated SNP, rs28999109. Furthermore, a patient from our sample trios was observed with a maternal segmental uniparental isodisomy (UPD) in a region containing rs17035464. Our findings support the involvement of PDGF-C in the development of oral clefts; moreover, the UPD case report contributes to the collective knowledge of rare variants in the human genome.

Interaction between IRF6 and TGFA genes contribute to the risk of nonsyndromic cleft lip/palate

Previous evidence from tooth agenesis studies suggested IRF6 and TGFA interact. Since tooth agenesis is commonly found in individuals with cleft lip/palate (CL/P), we used four large cohorts to evaluate if IRF6 and TGFA interaction contributes to CL/P. Markers within and flanking IRF6 and TGFA genes were tested using Taqman or SYBR green chemistries for case-control analyses in 1,000 Brazilian individuals. We looked for evidence of gene-gene interaction between IRF6 and TGFA by testing if markers associated with CL/P were overtransmitted together in the case-control Brazilian dataset and in the additional family datasets. Genotypes for an additional 142 case-parent trios from South America drawn from the Latin American Collaborative Study of Congenital Malformations (ECLAMC), 154 cases from Latvia, and 8,717 individuals from several cohorts were available for replication of tests for interaction. Tgfa and Irf6 expression at critical stages during palatogenesis was analyzed in wild type and Irf6 knockout mice. Markers in and near IRF6 and TGFA were associated with CL/P in the Brazilian cohort (p<10⁻⁶). IRF6 was also associated with cleft palate (CP) with impaction of permanent teeth (p<10⁻⁶). Statistical evidence of interaction between IRF6 and TGFA was found in all data sets (p = 0.013 for Brazilians; p = 0.046 for ECLAMC; p = 10⁻⁶ for Latvians, and p = 0.003 for the 8,717 individuals). Tgfa was not expressed in the palatal tissues of Irf6 knockout mice. IRF6 and TGFA contribute to subsets of CL/P with specific dental anomalies. Moreover, this potential IRF6-TGFA interaction may account for as much as 1% to 10% of CL/P cases. The Irf6-knockout model further supports the evidence of IRF6-TGFA interaction found in humans.

Detection of parent-of-origin specific expression quantitative trait loci by cis-association analysis of gene expression in trios

Parent-of-origin (PofO) effects, such as imprinting are a phenomenon in which homologous chromosomes exhibit differential gene expression and epigenetic modifications according to their parental origin. Such non-Mendelian inheritance patterns are generally ignored by conventional association studies, as these tests consider the maternal and paternal alleles as equivalent. To identify regulatory regions that show PofO effects on gene expression (imprinted expression Quantitative Trait Loci, ieQTLs), here we have developed a novel method in which we associate SNP genotypes of defined parental origin with gene expression levels. We applied this method to study 59 HapMap phase II parent-offspring trios. By analyzing mother/father/child trios, rules of Mendelian inheritance allowed the parental origin to be defined for ∼95% of SNPs in each child. We used 680,475 informative SNPs and corresponding expression data for 92,167 probe sets from Affymetrix GeneChip Human Exon 1.0 ST arrays and performed four independent cis-association analyses with the expression level of RefSeq genes within 1 Mb using PLINK. Independent analyses of maternal and paternal genotypes identified two significant cis-ieQTLs (p<10⁻⁷) at which expression of genes SFT2D2 and SRRT associated exclusively with maternally inherited SNPs rs3753292 and rs6945374, respectively. 28 additional suggestive cis-associations with only maternal or paternal SNPs were found at a lower stringency threshold of p<10⁻⁶, including associations with two known imprinted genes PEG10 and TRAPPC9, demonstrating the efficacy of our method. Furthermore, comparison of our method that utilizes independent analyses of maternal and paternal genotypes with the Likelihood Ratio Test (LRT) showed it to be more effective for detecting imprinting effects than the LRT. Our method represents a novel approach that can identify imprinted regulatory elements that control gene expression, suggesting novel PofO effects in the human genome.

Genetic and environmental risk factors for submucous cleft palate

A multifactorial aetiology with genetic and environmental factors is assumed for orofacial clefts. Submucous cleft palate (SMCP), a subgroup of cleft palates with insufficient median fusion of the muscles of the soft palate hidden under the mucosa, has a prevalence of 1:1,250-1:5,000. We described the prevalence of risk factors among 103 German patients with the subtype SMCP and genotyped 24 single nucleotide polymorphisms (SNPs) from 12 candidate genes for orofacial clefts. Analysis of risk factors yielded a positive history for maternal cigarette smoking during pregnancy in 25.2% of the patients, and this was significantly more frequent than in the normal population. The group of patients differed in allele frequencies at SNP rs3917192 of the gene TGFB3 (nominal P = 0.053) and at SNP rs5752638 of the gene MN1 (nominal P = 0.075) compared with 279 control individuals. Our results indicate a potential role of maternal smoking during pregnancy for the formation of SMCP. The analysis of genetic variants hints at the contribution of TGFB3 and MN1 in the aetiology of SMCPs.

The evolution of human genetic studies of cleft lip and cleft palate

Orofacial clefts (OFCs)--primarily cleft lip and cleft palate--are among the most common birth defects in all populations worldwide, and have notable population, ethnicity, and gender differences in birth prevalence. Interest in these birth defects goes back centuries, as does formal scientific interest; scientists often used OFCs as examples or evidence during paradigm shifts in human genetics, and have also used virtually every new method of human genetic analysis to deepen our understanding of OFC. This review traces the evolution of human genetic investigations of OFC, highlights the specific insights gained about OFC through the years, and culminates in a review of recent key OFC genetic findings resulting from the powerful tools of the genomics era. Notably, OFC represents a major success for genome-wide approaches, and the field is poised for further breakthroughs in the near future.

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.

Male and female differential reproductive rate could explain parental transmission asymmetry of mutation origin in Hirschsprung disease

Hirschsprung disease (HSCR, aganglionic megacolon) is a complex and heterogeneous disease with an incidence of 1 in 5000 live births. Despite the multifactorial determination of HSCR in the vast majority of cases, there is a monogenic subgroup for which private rare RET coding sequence mutations with high penetrance are found (45% of HSCR familial cases). An asymmetrical parental origin is observed for RET coding sequence mutations with a higher maternal inheritance. A parent-of-origin effect is usually assumed. Here we show that a differential reproductive rate for males and females also leads to an asymmetrical parental origin, which was never considered as a possible explanation till now. In the case of HSCR, we show a positive association between penetrance of the mutation and parental transmission asymmetry: no parental transmission asymmetry is observed in sporadic RET CDS mutation carrier cases for which penetrance of the mutation is low, whereas a parental transmission asymmetry is observed in affected sib-pairs for which penetrance of the mutation is higher. This allows us to conclude that the explanation for this parental asymmetry is that more severe mutations have resulted in a differential reproductive rate between male and female carriers.

Evidence of gene-environment interaction for the RUNX2 gene and environmental tobacco smoke in controlling the risk of cleft lip with/without cleft palate

This study examined the association between 49 markers in the Runt-related transcription factor 2 (RUNX2) gene and nonsyndromic cleft lip with/without cleft palate (CL/P) among 326 Chinese case-parent trios, while considering gene-environment (GxE) interaction and parent-of-origin effects. Five single-nucleotide polymorphisms (SNPs) showed significant evidence of linkage and association with CL/P and these results were replicated in an independent European sample of 825 case-parent trios. We also report compelling evidence for interaction between markers in RUNX2 and environmental tobacco smoke (ETS). Although most marginal SNP effects (i.e., ignoring maternal exposures) were not statistically significant, eight SNPs were significant when considering possible interaction with ETS when testing for gene (G) and GxE interaction simultaneously or when considering GxE alone. Independent samples from European populations showed consistent evidence of significant GxETS interaction at two SNPs (rs6904353 and rs7748231). Our results suggest genetic variation in RUNX2 may influence susceptibility to CL/P through interacting with ETS.

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 × E) interaction between markers in FGFR2 and maternal smoking or multivitamin supplementation yielded significant evidence of G × E interaction separately. Tests of gene-gene (G × G) interaction using Cordell's 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 bronchodilator use and the risk of orofacial clefts

BACKGROUND

Few epidemiological studies have explored the relationship between orofacial clefts and bronchodilators. We assessed whether mothers who used bronchodilators during early pregnancy were at an increased risk of delivering infants with orofacial clefts.

METHODS

We used National Birth Defects Prevention Study case-control data from mothers of 2711 infants with orofacial clefts and 6482 mothers of live born infants without birth defects, delivered during 1997 through 2005. Information on medication use from 3 months before pregnancy through delivery was collected using a standardized interview. Logistic regression was used to estimate adjusted odds ratios (aOR) and 95% confidence intervals (CIs) for maternal bronchodilator use during the periconceptional period (1 month before pregnancy through the third month of pregnancy) while controlling for other covariates.

RESULTS

We observed an association between maternal bronchodilator use during the periconceptional period and cleft lip only (CLO) (aOR = 1.77, 95% CI: 1.08-2.88). The risk of cleft palate only (CPO) (aOR = 1.53, 95% CI: 0.99-2.37) was elevated but was not statistically significant. No association was observed for maternal bronchodilator use and the risk of cleft lip with cleft palate (aOR = 0.78, 95% CI: 0.46-1.31). The most commonly used bronchodilator was albuterol (88.7%). Maternal albuterol use was associated with CLO (aOR = 1.79, 95% CI: 1.07-2.99) and CPO (aOR = 1.65, 95% CI: 1.06-2.58).

CONCLUSIONS

We observed a statistically significant association between maternal bronchodilator use during the periconceptional period and the risk of CLO after controlling for other risk factors. It is unclear whether the increased odds ratios observed in this study are due to the bronchodilators, the severity of asthma, or both, or to chance alone. Further studies to disentangle the role of asthma or asthma medications would help clarify these findings.

Genetics of nonsyndromic orofacial clefts

With an average worldwide prevalence of approximately 1.2/1000 live births, orofacial clefts are the most common craniofacial birth defects in humans. Like other complex disorders, these birth defects are thought to result from the complex interplay of multiple genes and environmental factors. Significant progress in the identification of underlying genes and pathways has benefited from large populations available for study, increased international collaboration, rapid advances in genotyping technology, and major improvements in analytic approaches. Here we review recent advances in genetic epidemiological approaches to complex traits and their applications to studies of nonsyndromic orofacial clefts. Our main aim is to bring together a discussion of new and previously identified candidate genes to create a more cohesive picture of interacting pathways that shape the human craniofacial region. In future directions, we highlight the need to search for copy number variants that affect gene dosage and rare variants that are possibly associated with a higher disease penetrance. In addition, sequencing of protein-coding regions in candidate genes and screening for genetic variation in noncoding regulatory elements will help advance this important area of research.

IRF6 polymorphisms are associated with nonsyndromic orofacial clefts in a Chinese Han population

IRF6 plays an important role in orofacial development. In the present study, we genotyped two polymorphisms (rs642961 and rs2235371) within the IRF6 locus and estimated their associations with risk of nonsyndromic orofacial clefts (NSOC), including the subgroups, in a hospital-based case-control study in a Chinese Han population. In the single locus analyses, we found rs642961 AG and AG/AA genotypes were associated with increased risk of NSOC, especially cleft lip with or without cleft palate (CL/P) and cleft lip with cleft palate (CLP), while significantly decreased risks were associated with rs2235371 CT and CT/TT genotypes. When examining the combined effects of these two polymorphisms and using the rs642961 A and rs2235371 C alleles as the risk alleles, we found genotypes containing 2-4 risk alleles conferred high risk to NSOC, CL/P, and CLP. Furthermore, to test whether rs642961 could modulate IRF6 expression in vivo, we surgically collected lip skin tissues within the adjacent region of lip cleft site and found rs642961 genotypes were associated with differential levels of IRF6 mRNA and protein expression in an allele-dosage manner, providing the first evidence that rs642961 affected IRF6 expression in vivo. Taken together, these findings confirm the contribution of IRF6 genetic variants in the etiology of NSOC in a Chinese Han population.

Evidence of gene-environment interaction for the IRF6 gene and maternal multivitamin supplementation in controlling the risk of cleft lip with/without cleft palate

Although multiple genes have been identified as genetic risk factors for isolated, non-syndromic cleft lip with/without cleft palate (CL/P), a complex and heterogeneous birth defect, interferon regulatory factor 6 gene (IRF6) is one of the best documented genetic risk factors. In this study, we tested for association between markers in IRF6 and CL/P in 326 Chinese case-parent trios, considering gene-environment interaction for two common maternal exposures, and parent-of-origin effects. CL/P case-parent trios from three sites in mainland China and Taiwan were genotyped for 22 single nucleotide polymorphisms (SNPs) in IRF6. The transmission disequilibrium test was used to test for marginal effects of individual SNPs. We used PBAT to screen the SNPs and haplotypes for gene-environment (G×E) interaction and conditional logistic regression models to quantify effect sizes for SNP-environment interaction. After Bonferroni correction, 14 SNPs showed statistically significant association with CL/P. Evidence of G×E interaction was found for both maternal exposures, multivitamin supplementation and environmental tobacco smoke (ETS). Two SNPs showed evidence of interaction with multivitamin supplementation in conditional logistic regression models (rs2076153 nominal P=0.019, rs17015218 nominal P=0.012). In addition, rs1044516 yielded evidence for interaction with maternal ETS (nominal P=0.041). Haplotype analysis using PBAT also suggested interaction between SNPs in IRF6 and both multivitamin supplementation and ETS. However, no evidence for maternal genotypic effects or significant parent-of-origin effects was seen in these data. These results suggest IRF6 gene may influence risk of CL/P through interaction with multivitamin supplementation and ETS in the Chinese population.

FOXE1 association with both isolated cleft lip with or without cleft palate, and isolated cleft palate

Nonsyndromic orofacial clefts are a common complex birth defect caused by genetic and environmental factors and/or their interactions. A previous genome-wide linkage scan discovered a novel locus for cleft lip with or without cleft palate (CL/P) at 9q22-q33. To identify the etiologic gene, we undertook an iterative and complementary fine mapping strategy using family-based CL/P samples from Colombia, USA and the Philippines. Candidate genes within 9q22-q33 were sequenced, revealing 32 new variants. Concurrently, 397 SNPs spanning the 9q22-q33 2-LOD-unit interval were tested for association. Significant SNP and haplotype association signals (P = 1.45E - 08) narrowed the interval to a 200 kb region containing: FOXE1, C9ORF156 and HEMGN. Association results were replicated in CL/P families of European descent and when all populations were combined the two most associated SNPs, rs3758249 (P = 5.01E - 13) and rs4460498 (P = 6.51E - 12), were located inside a 70 kb high linkage disequilibrium block containing FOXE1. Association signals for Caucasians and Asians clustered 5' and 3' of FOXE1, respectively. Isolated cleft palate (CP) was also associated, indicating that FOXE1 plays a role in two phenotypes thought to be genetically distinct. Foxe1 expression was found in the epithelium undergoing fusion between the medial nasal and maxillary processes. Mutation screens of FOXE1 identified two family-specific missense mutations at highly conserved amino acids. These data indicate that FOXE1 is a major gene for CL/P and provides new insights for improved counseling and genetic interaction studies.