Abstract:Genome-wide association studies (GWASs) are the most widely used method to identify genetic risk loci associated with orofacial clefts (OFC). However, despite the increasing size of cohort, GWASs are still insufficient to detect all the heritability, suggesting there are more associations under the current stringent statistical threshold. In this study, we obtained an integrated epigenomic dataset based on the chromatin conformation of a human oral epithelial cell line (HIOEC) using RNA-seq, ATAC-seq, H3K27ac ChIP-seq, and DLO Hi-C. Presumably, this epigenomic dataset could reveal the missing functional variants located in the oral epithelial cell active enhancers/promoters along with their risk target genes, despite relatively less-stringent statistical association with OFC. Taken a non-syndromic cleft palate only (NSCPO) GWAS data of the Chinese Han population as an example, 3664 SNPs that cannot reach the strict significance threshold were subjected to this functional identification pipeline. In total, 254 potential risk SNPs residing in active cis-regulatory elements interacting with 1 718 promoters of oral epithelium-expressed genes were screened. Gapped k-mer machine learning based on enhancers interacting with epithelium-expressed genes along with in vivo and in vitro reporter assays were employed as functional validation. Among all the potential SNPs, we chose and confirmed that the risk alleles of rs560789 and rs174570 reduced the epithelial-specific enhancer activity by preventing the binding of transcription factors related to epithelial development. In summary, we established chromatin conformation datasets of human oral epithelial cells and provided a framework for testing and understanding how regulatory variants impart risk for clefts.

Genetic association and functional validation of ZFP36L2 in non-syndromic orofacial cleft subtypes

Integrating GWAS and eQTL to predict genes and pathways for non-syndromic cleft lip with or without palate

Genome-wide analyses of non-syndromic cleft lip with palate identify 14 novel loci and genetic heterogeneity

Integrating Transcriptomics and Genomics to Identify Fibroblast Growth Factor/receptor Candidate Genes for Non‐syndromic Orofacial Cleft in Chinese

Associations between microRNA binding site SNPs in FGFs and FGFRs and the risk of non-syndromic orofacial cleft

Targeted Re-Sequencing of the 2p21 Locus Identifies Non-Syndromic Cleft Lip Only Novel Susceptibility Gene ZFP36L2

Multi-ancestry Genome Wide Association Study Meta-analysis of Non-syndromic Orofacial Clefts

Chromatin conformation of human oral epithelium can identify orofacial cleft missing functional variants

Identification of Common Non-Coding Variants at 1p22 That Are Functional for Non-Syndromic Orofacial Clefting

An Association Study of Candidate Gene Variants in Chinese Nonsyndromic Cleft Lip with or without Palate Subjects

Unveiling dysregulated lncRNAs and networks in non-syndromic cleft lip with or without cleft palate pathogenesis

Genetic Variant in Lncrna RP11-462G12.2 Confers the Risk of Non-Syndromic Orofacial Clefts

Variants in miRNA regulome and their association with the risk of nonsyndromic orofacial clefts.

Genome-Wide Analyses Identify A Novel Risk Locus For Nonsyndromic Cleft Palate

Author Response for "integrating GWAS and Eqtl to Predict Genes and Pathways for Non‐syndromic Cleft Lip with or Without Palate"

Role of Epigenetics and Mirnas in Orofacial Clefts

Genome-wide association study identifies a new susceptibility locus for cleft lip with or without a cleft palate.

Identification of Functional Variants for Cleft Lip with or Without Cleft Palate in or Near PAX7, FGFR2, and NOG by Targeted Sequencing of GWAS Loci.

Combining genetic and single-cell expression data reveals cell types and novel candidate genes for orofacial clefting

Association Analysis of SNPs in GRHL3, FAF1, and KCNJ2 with NSCPO Sub-Phenotypes in Han Chinese

Identification of rare variants in PTCH2 associated with non-syndromic orofacial clefts