Ethnic heterogeneity of IRF6 AP-2a binding site promoter SNP association with nonsyndromic cleft lip and palate

Association Between SNPs in 1q32.2 and NSCL ± P in Han Chinese Population

OBJECTIVES

Non-syndromic cleft lip with or without cleft palate (NSCL ± P) is one of the most common birth malformations. Currently, numerous susceptibility SNPs have been reported by GWA studies, however, the replications of them among NSCL ± P from Han Chinese were very limited.

DESIGN

In this study, we selected 16 SNPs around 1q32.2 based on the published GWA studies and replicated them among 302 trios with NSCL ± P from Han Chinese Population. The genotypic data was analyzed with FBAT, PLINK and R package.

SETTING

The study was conducted in a tertiary medical center.

PATIENTS, PARTICIPANTS

302 patients with CL ± P and their parents.

MAIN OUTCOME MEASURES

To ascertain the genetic variants in 1q32.2 in patients with CL ± P in Han Chinese Population.

INTERVENTIONS

Blood samples were collected.

RESULTS

We found T allele (Z = 4.26, p = 0.00002) and T/T homozygotes (Z = 4.4, p = 0.000011) at rs12063989 was significantly over-transmitted among non-syndromic cleft lip with or without cleft palate (NSCL ± P).

CONCLUSIONS

We found rs12063989 exhibited significant association with the occurrence of NSCL ± P, which would provide new evidence for the future study in the etiology of NSCL ± P.

PBX-WNT-P63-IRF6 pathway in nonsyndromic cleft lip and palate

Nonsyndromic cleft lip and palate (NSCLP) is one of the most common craniofacial anomalies in humans, affecting more than 135,000 newborns worldwide. NSCLP has a multifactorial etiology with more than 50 genes postulated to play an etiologic role. The genetic pathway comprised of Pbx-Wnt-p63-Irf6 genes was shown to control facial morphogenesis in mice and proposed as a regulatory pathway for NSCLP. Based on these findings, we investigated whether variation in PBX1, PBX2, and TP63, and their proposed interactions were associated with NSCLP. Fourteen single nucleotide variants (SNVs) in/nearby PBX1, PBX2, and TP63 were genotyped in 780 NSCLP families of nonHispanic white (NHW) and Hispanic ethnicities. Family-based association tests were performed for individual SNVs stratified by ethnicity and family history of NSCLP. Gene-gene interactions were also tested. A significant association was found for PBX2 rs3131300 and NSCLP in combined Hispanic families (p = .003) while nominal association was found for TP63 rs9332461 in multiplex Hispanic families (p = .005). Significant haplotype associations were observed for PBX2 in NHW (p = .0002) and Hispanic families (p = .003), and for TP63 in multiplex Hispanic families (.003). An independent case-control group was used to validate findings, and significant associations were found with PBX1 rs6426870 (p = .007) and TP63 rs9332461 (p = .03). Gene-gene interactions were detected between PBX1/PBX2/TP63 with IRF6 in NHW families, and between PBX1 with WNT9B in both NHW and Hispanic families (p < .0018). This study provides the first evidence for a role of PBX1 and PBX2, additional evidence for the role of TP63, and support for the proposed PBX-WNT-TP63-IRF6 regulatory pathway in the etiology of NSCLP.

Association Studies Between Regulatory Regions of IRF6/TP63 Genes and Nonsyndromic Oral Clefts

OBJECTIVE

To evaluate genetic variants within the regulatory regions of interferon regulatory factor 6 (IRF6) and TP63 for the etiology of nonsyndromic oral clefts risk factors.

DESIGN

We performed allelic transmission disequilibrium test analysis on 5 eligible single-nucleotide polymorphisms (SNPs) and SNP haplotypes using the Family-Based Association Test.

PARTICIPANTS

The study sample consisted of 334 case-parent trios of nonsyndromic oral clefts from Taiwanese population, separated into nonsyndromic cleft lip/palate (NSCL/P) and nonsyndromic cleft palate only (NSCPO) groups.

RESULTS

We found all 3 selected SNPs of the IRF6 gene show significant association with nonsyndromic oral clefts (rs2235371, P = 5.10E-07; rs642961, P = .00194; and rs77542756, P = 9.08E-07). Haplotype analyses identified 3 possible SNP combination haplotypes in the IRF6 gene and found that C-G-G showed significant undertransmission (P = .058), whereas 2 other haplotypes, T-G-A and C-A-G (P = 2.71E-06 and P = 5.00E-04, respectively), were significantly overtransmitted to the NSCL/P children but not to the NSCPO children. For the TP63 gene, we failed to detect evidence of nonsyndromic oral cleft association in the 2 SNPs within the TP63 large intron 1 region.

CONCLUSIONS

We used a family-based analysis in 334 Taiwanese case-parent trios to evaluate selected SNPs of IRF6 genes and TP63 genes for a risk of orofacial clefting. This study provides additional evidence for an association between IRF6 and NSCL/P, including the genetic variants within the 5'-noncoding region of the gene. We also confirmed that NSCL/P and NSCPO individuals belong to different groups. For the TP63, our data did not favor the direct involvement of TAp63 isoforms during orofacial development.

SNP@lincTFBS: an integrated database of polymorphisms in human LincRNA transcription factor binding sites

Large intergenic non-coding RNAs (lincRNAs) are a new class of functional transcripts, and aberrant expression of lincRNAs was associated with several human diseases. The genetic variants in lincRNA transcription factor binding sites (TFBSs) can change lincRNA expression, thereby affecting the susceptibility to human diseases. To identify and annotate these functional candidates, we have developed a database SNP@lincTFBS, which is devoted to the exploration and annotation of single nucleotide polymorphisms (SNPs) in potential TFBSs of human lincRNAs. We identified 6,665 SNPs in 6,614 conserved TFBSs of 2,423 human lincRNAs. In addition, with ChIPSeq dataset, we identified 139,576 SNPs in 304,517 transcription factor peaks of 4,813 lincRNAs. We also performed comprehensive annotation for these SNPs using 1000 Genomes Project datasets across 11 populations. Moreover, one of the distinctive features of SNP@lincTFBS is the collection of disease-associated SNPs in the lincRNA TFBSs and SNPs in the TFBSs of disease-associated lincRNAs. The web interface enables both flexible data searches and downloads. Quick search can be query of lincRNA name, SNP identifier, or transcription factor name. SNP@lincTFBS provides significant advances in identification of disease-associated lincRNA variants and improved convenience to interpret the discrepant expression of lincRNAs. The SNP@lincTFBS database is available at http://bioinfo.hrbmu.edu.cn/SNP_lincTFBS.

Skin Tissue Surface Morphology and Quality of RNA and Protein Extracted from Fresh and Stabilized Human Cleft Lip and Palate Tissue

Cleft lip palate is a human congenital disorder worldwide and the study of this genetic disease requires molecular genetic analysis. This analysis required the use of nucleic acid and protein, thus good quality and quantity of its extraction is important. We are comparing RNA and protein extractions from tissue biopsy of cleft lip palate in both fresh and stabilized condition. Tissue morphology was also captured using Scanning Electron Microscopy (SEM) for any morphology differences. Tissue homogenization may destroy nucleic acid stability but not its morphology. Low RNA concentration from stabilized tissue was found. However, there was no crucial issue of protein extraction, degradation or concentration. Tissue morphology was slight different between normal and CL/P tissue.

Association and Mutation Analyses of the IRF6 Gene in Families with Nonsyndromic and Syndromic Cleft Lip and/or Cleft Palate

Objectives

(1) To detect interferon regulatory factor 6 gene ( IRF6) mutations in newly recruited Van der Woude syndrome (VWS) and popliteal pterygium syndrome (PPS) families. (2) To test for association, in nonsyndromic cleft lip and/or cleft palate (NSCL/P) and in VWS/PPS families, the single nucleotide polymorphism (SNP) rs642961, from the IRF6 enhancer AP-2α region, alone or as haplotype with rs2235371, a coding SNP (Val274Ile).

Design

IRF6 mutation screening was performed by direct sequencing and genotyping of rs642961 and rs2235371 by TaqMan technology.

Patients

Seventy-one Swedish NSCL/P families, 24 Finnish cleft palate (CP) families, and 24 VWS/PPS families (seven newly recruited) were studied.

Results

Allelic and genotypic frequencies in each phenotype were compared to those of the controls, and no significant difference could be observed. IRF6 gene mutation was detected in six of the seven new VWS/PPS families. Association analysis of the entire VWS/PPS sample set revealed the A allele from rs642961 to be a risk allele. Significant association was detected in the Swedish CP subset of our NSCL/P collection where the G-C haplotype for rs642961-rs2235371 were at risk ( P= .013).

Conclusions

Our results do not support the previously reported association between the A allele of rs642961 and the NSCL phenotype. However, in the VWS/PPS families, the A allele was a risk allele and was, in a large majority (>80%), transmitted on the same chromosome as the IRF6 mutation.

Polymorphisms in FGF12, VCL, CX43 and VAX1 in Brazilian patients with nonsyndromic cleft lip with or without cleft palate

BACKGROUND

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is the most common orofacial birth defect with a wide range prevalence among different populations. Previous association studies with populations from Europe and Asia have identified putative susceptibility markers for NSCL/P in fibroblast growth factor 12 (FGF12), vinculin (VCL), connexin 43 (CX43) and in a region close to the ventral anterior homeobox 1 (VAX1) gene. However, there have thus far been no studies of these markers in NSCL/P Brazilian patients, and as the genetic ancestry of the Brazilian population is highly varied, the predisposition to those disease markers can be different.

METHODS

Herein we conducted a structured association study conditioned on the individual ancestry proportions to determine the role of 16 polymorphic markers within those genes in 300 patients with NSCL/P and 385 unaffected controls.

RESULTS

None of the alleles and genotypes showed association with NSCL/P, though there was a significant association of the haplotype formed by VAX1 rs10787760, rs6585429 and rs1871345 polymorphisms with NSCL/P that did not persist Bonferroni correction for multiple tests.

CONCLUSIONS

Our results are consistent with a lack of involvement of FGF12, VCL and CX43 variants with NSCL/P pathogenesis in Brazilian patients. Furthermore, the higher frequency of a haplotype of VAX1 with NSCL/P patients suggests a low penetrant gene for oral cleft, and warrants further studies.

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.

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.

Region 8q24 is a susceptibility locus for nonsyndromic oral clefting in Brazil

BACKGROUND

Nonsyndromic cleft lip with or without cleft palate is a relatively common craniofacial defect with multifactorial inheritance. The association of the rs987525 single nucleotide variant, located in a gene desert at 8q24.21 region, has been consistently replicated in European populations. We performed a structured association approach combined with transcriptional analysis of the MYC gene to dissect the role of rs987525 in oral clefting susceptibility in the ethnically admixed Brazilian population.

METHODS

We performed the association study conditioned on the individual ancestry proportions in a sample of 563 patients and 336 controls, and in an independent sample of 221 patients and 261 controls. The correlation between rs987525 genotypes and MYC transcriptional levels in orbicularis oris muscle mesenchymal stem cells was also investigated in 42 patients and 4 controls.

RESULTS

We found a significant association in the larger sample (p = 0.0016; OR = 1.80 [95% confidence interval {CI}, 1.21-2.69], for heterozygous genotype, and 2.71 [95% CI, 1.47-4.96] for homozygous genotype). We did not find a significant correlation between rs987525 genotypes and MYC transcriptional levels (p = 0.14; r = -0.22, Spearman Correlation).

CONCLUSIONS

We present a positive association of rs987525 in the Brazilian population for the first time, and it is likely that the European contribution to our population is driving this association. We also cannot discard a role of rs987515 in MYC regulation, because this locus behaves as an expression quantitative locus of MYC in another tissue.

Cleft lip and palate: understanding genetic and environmental influences

Clefts of the lip and/or palate (CLP) are common birth defects of complex aetiology. CLP can occur in isolation or as part of a broad range of chromosomal, Mendelian or teratogenic syndromes. Although there has been marked progress in identifying genetic and environmental triggers for syndromic CLP, the aetiology of the more common non-syndromic (isolated) forms remains poorly characterized. Recently, using a combination of epidemiology, careful phenotyping, genome-wide association studies and analysis of animal models, several distinct genetic and environmental risk factors have been identified and confirmed for non-syndromic CLP. These findings have advanced our understanding of developmental biology and created new opportunities for clinical translational research.