Irf6-Related Gene Regulatory Network Involved in Palate and Lip Development

Uncovering the Pathogenesis of Orofacial Clefts Using Bioinformatics Analysis

OBJECTIVE

Many genes have been found to be associated with the occurrence of the orofacial clefts (OFC). The links between these pathogenic genes are rarely studied. In this study, bioinformatics analysis were performed in order to find associations between OFC-related genes and provide new ideas for etiology study of OFCs.

METHODS

Orofacial clefts-related genes were searched and identified from the Online Mendelian Inheritance of Man (OMIM.org). These genes were then analyzed by bioinformatics methods, including protein-protein interaction network, functional enrichment analysis, module analysis, and hub genes analysis.

RESULTS

After searching the database of OMIM.org and removing duplicate results, 279 genes were finally obtained. These genes were involved to 369 pathways in biological process, 56 in cell component, 64 in molecular function, and 45 in the Kyoto Encyclopedia of Genes and Genomes. Most identified genes were significantly enriched in embryonic appendage morphogenesis (29.17%), embryonic limb morphogenesis (6.06%), and limb development (4.33%) for biological process (Fig. 5A); ciliary tip (42.86%), MKS complex (28.57%), ciliary basal body (14.29%), and ciliary membrane (14.29%) for cell component. The top 10 hub genes were identified, including SHH, GLI2, PTCH1, SMAD4, FGFR1, BMP4, SOX9, SOX2, RUNX2, and CDH1.

CONCLUSIONS

Bioinformatics methods were used to analyze OFC-related genes in this study, including hub gene identifying and analysis, protein-protein interaction network construction, and functional enrichment analysis. Several potential mechanisms related to occurrence of OFCs were also discussed. These results may be helpful for further studies of the etiology of OFC.

Comprehensive analysis of differentially expressed profiles of non‑coding RNAs in peripheral blood and ceRNA regulatory networks in non‑syndromic orofacial clefts

Non-syndromic orofacial clefts (NSOC), which include cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO), are common congenital birth defects in humans. Accumulating evidence indicates that long non-coding RNAs (lncRNAs) and microRNAs (miRNAs or miRs) play important roles in NSOC; however, the potential regulatory associations between them remain largely unknown. In this study, we performed next-generation RNA sequencing (RNA-seq) to identify transcriptome profiles, including mRNAs, lncRNAs and miRNAs, in patients with CL/P and CPO. A total of 36 lncRNAs, 1,341 mRNAs and 60 miRNAs were found to be differentially expressed in the CL/P group compared to the control group, and 57 lncRNAs, 1,255 mRNAs and 162 miRNAs were found to be differentially expressed in the CPO group compared to the control group. Subsequently, reverse transcription-quantitative polymerase chain reaction (RT-qPCR) was performed to validate the expression of selected lncRNAs, miRNAs and mRNAs. In addition, bioinformatics methods were employed to explore the potential functions of ncRNAs and to construct lncRNA-miRNA-mRNA regulatory networks. To the best of our knowledge, this is the first study to comprehensively analyze regulated non-coding RNAs (ncRNAs) in CL/P and CPO, providing a novel perspective on the etiology of NSOC and laying the foundation for future research into the potential regulatory mechanisms of ncRNAs and mRNAs in NSOC.

Systems genetics of nonsyndromic orofacial clefting provides insights into its complex aetiology

Nonsyndromic oral clefting (NSOC) is although one of the most common congenital disorders worldwide, its underlying molecular basis remains elusive. This process has been hindered by the overwhelmingly high level of heterogeneity observed. Given that hitherto multiple loci and genes have been associated with NSOC, and that complex diseases are usually polygenic and show a considerable level of missing heritability, we used a systems genetics approach to reconstruct the NSOC network by integrating human-based physical and regulatory interactome with whole-transcriptome microarray data. We show that the network component contains 53% (23/43) of the curated NSOC-implicated gene set and displays a highly significant propinquity (P < 0.0001) between genes implicated at the genomic level and those differentially expressed at the transcriptome level. In addition, we identified bona fide candidate genes based on topological features and dysregulation (e.g. ANGPTL4), and similarly prioritised genes at GWA loci (e.g. MYC and CREBBP), thus providing further insight into the underlying heterogeneity of NSOC. Gene ontology analysis results were consistent with the NSOC network being associated with embryonic organ morphogenesis and also hinted at an aetiological overlap between NSOC and cancer. We therefore recommend this approach to be applied to other heterogeneous complex diseases to not only provide a molecular framework to unify genes which may seem as disparate entities linked to the same disease, but to also predict and prioritise candidate genes for further validation, thus addressing the missing heritability.

Gene-gene interaction for nonsyndromic cleft lip with or without cleft palate in Chilean case-parent trios

OBJECTIVE

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a birth defect for which several genes susceptibility genes been proposed. Consequently, it has been suggested that many of these genes belong to common inter-related pathways during craniofacial development gene-gene interaction. We evaluated the presence of gene-gene interaction for single nucleotide polymorphisms within interferon regulatory factor 6 (IRF6), muscle segment homeobox 1 (MSX1), bone morphogenetic protein 4 (BMP4) and transforming growth factor 3 (TGFB3) genes in NSCL/P risk in Chilean case-parent trios.

DESIGN

From previous studies, we retrieved genotypes for 13 polymorphic variants within these four genes in 152 case-parent trios. Using the trio package (R) we evaluate the gene-gen interaction in genetic markers pairs applying a 1°-of-freedom test (1df) and a confirmatory 4°-of-freedom (4df) test for epistasis followed by both a permutation test and a Benjamini-Hochberg test for multiple comparisons adjustment.

RESULTS

We found evidence of gene-gene interaction for rs6446693 (MSX1) and rs2268625 (TGFB3) (4df p = 0.024; permutation p = 0.015, Benjamini-Hochberg p = 0.001).

CONCLUSIONS

A significant gene-gene interaction was detected for rs6446693 (MSX1) and rs2268625 (TGFB3). This finding is concordant with research in animal models showing that MSX1 and TGFB3 are expressed in common molecular pathways acting in an epistatic manner during maxillofacial development.

Gene Interactions Provide Evidence for Signaling Pathways Involved in Cleft Lip/Palate in Humans

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a common craniofacial birth defect that has a complex etiology. Genome-wide association studies have recently identified new loci associated with NSCL/P, but these loci have not been analyzed in a Mexican Mestizo population. A complex etiology implies the presence of genetic interactions, but there is little available information regarding this in NSCL/P, and no signaling pathway has been clearly implicated in humans. Here, we analyzed the associations of 24 single nucleotide polymorphisms (SNPs) with NSCL/P in a Mexican Mestizo population (133 cases, 263 controls). The multifactorial dimensionality reduction method was used to examine gene-gene and gene-folic acid consumption interactions for the 24 SNPs analyzed in this study and for 2 additional SNPs that had previously been genotyped in the same study population. Six SNPs located in paired box 7, ventral anterior homeobox 1, sprouty RTK signaling antagonist 2, bone morphogenetic protein 4, and tropomyosin 1 genes were associated with higher risks of NSCL/P (P = 0.0001 to 0.04); 2 SNPs, 1 each in netrin 1 and V-maf avian musculoaponeurotic fibrosarcoma oncogene homolog B, were associated with a lower risk of NSCL/P (P = 0.013 to 0.03); and 2 SNPs, 1 each in ATP binding cassette subfamily A member 4 (ABCA4) and noggin, showed associations with NSCL/P that approached the threshold of significance (P = 0.056 to 0.07). In addition, 6 gene-gene interactions (P = 0.0001 to 0.001) and an ABCA4-folic acid consumption interaction (P < 0.0001) were identified. On the basis of these results, combined with those of previous association studies in the literature and biological characterizations of murine models, we propose an interaction network in which interferon regulatory factor 6 plays a central role in the etiology of NSCL/P.

Toward an orofacial gene regulatory network

Orofacial clefting is a common birth defect with significant morbidity. A panoply of candidate genes have been discovered through synergy of animal models and human genetics. Among these, variants in interferon regulatory factor 6 (IRF6) cause syndromic orofacial clefting and contribute risk toward isolated cleft lip and palate (1/700 live births). Rare variants in IRF6 can lead to Van der Woude syndrome (1/35,000 live births) and popliteal pterygium syndrome (1/300,000 live births). Furthermore, IRF6 regulates GRHL3 and rare variants in this downstream target can also lead to Van der Woude syndrome. In addition, a common variant (rs642961) in the IRF6 locus is found in 30% of the world's population and contributes risk for isolated orofacial clefting. Biochemical studies revealed that rs642961 abrogates one of four AP-2alpha binding sites. Like IRF6 and GRHL3, rare variants in TFAP2A can also lead to syndromic orofacial clefting with lip pits (branchio-oculo-facial syndrome). The literature suggests that AP-2alpha, IRF6 and GRHL3 are part of a pathway that is essential for lip and palate development. In addition to updating the pathways, players and pursuits, this review will highlight some of the current questions in the study of orofacial clefting.