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

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.

References

[1] Leslie EJ, Marazita ML. Genetics of cleft lip and cleft palate. Am. J. Med. Genet. C. Semin. Med. Genet., 2013, 163C: 246-258.

[2] Dixon MJ, Marazita ML, Beaty TH, Murray JC. Cleft lip and palate: understanding genetic and environmental influences. Nat. Rev. Genet., 2011, 12: 167-178.

[3] Maarse W, . A systematic review of associated structural and chromosomal defects in oral clefts: when is prenatal genetic analysis indicated. J. Med. Genet., 2012, 49: 490-498.

[4] Beaty TH, . A genome-wide association study of cleft lip with and without cleft palate identifies risk variants near MAFB and ABCA4. Nat. Genet., 2010, 42: 525-529.

[5] Birnbaum S, . Key susceptibility locus for nonsyndromic cleft lip with or without cleft palate on chromosome 8q24. Nat. Genet., 2009, 41: 473-477.

[6] Mangold E, . Genome-wide association study identifies two susceptibility loci for nonsyndromic cleft lip with or without cleft palate. Nat. Genet., 2010, 42: 24-26.

[7] Sun Y, . Genome-wide association study identifies a new susceptibility locus for cleft lip with or without a cleft palate. Nat. Commun., 2015, 6

[8] Wolf ZT, . Genome-wide association studies in dogs and humans identify ADAMTS20 as a risk variant for cleft lip and palate. PLoS Genet., 2015, 11: e1005059.

[9] Yu Y, . Genome-wide analyses of non-syndromic cleft lip with palate identify 14 novel loci and genetic heterogeneity. Nat. Commun., 2017, 8

[10] Manolio TA, . Finding the missing heritability of complex diseases. Nature, 2009, 461: 747-753.

[11] Tam V, . Benefits and limitations of genome-wide association studies. Nat. Rev. Genet., 2019, 20: 467-484.

[12] Dudbridge F, Gusnanto A. Estimation of significance thresholds for genomewide association scans. Genet. Epidemiol., 2008, 32: 227-234.

[13] Gusev A, . Integrative approaches for large-scale transcriptome-wide association studies. Nat. Genet., 2016, 48: 245-252.

[14] Gallagher MD, Chen-Plotkin AS. The post-GWAS era: from association to function. Am. J. Hum. Genet., 2018, 102: 717-730.

[15] Mbuyi-Musanzayi S, . Non-syndromic cleft lip and/or cleft palate: epidemiology and risk factors in Lubumbashi (DR Congo), a case-control study. J. Craniomaxillofac Surg., 2018, 46: 1051-1058.

[16] Lane J, Kaartinen V. Signaling networks in palate development. Wiley Interdiscip. Rev. Syst. Biol. Med., 2014, 6: 271-278.

[17] Li C, Lan Y, Jiang R. Molecular and cellular mechanisms of palate development. J. Dent. Res., 2017, 96: 1184-1191.

[18] Liu, H. et al. Analysis of zebrafish periderm enhancers facilitates identification of a regulatory variant near human KRT8/18. Elife 9, https://doi.org/10.7554/eLife.51325 (2020).

[19] Xiao Y, Liu H. Intrinsic differences between the open chromatin regions of oral and epidermal keratinocytes. Chin. J. Dent. Res., 2020, 23: 119-130.

[20] Wang, X. et al. Discovery and validation of sub-threshold genome-wide association study loci using epigenomic signatures. Elife 5, https://doi.org/10.7554/eLife.10557 (2016).

[21] Zhang Z, . [Establishment of human immortalized oral epithelial cell line HIO615 induced by HPV16 E6 and E7]. Zhonghua Kou Qiang Yi Xue Za Zhi, 2002, 37: 12-14.

[22] He M, . Genome-wide analyses identify a novel risk locus for nonsyndromic cleft palate. J. Dent. Res., 2020, 99: 1461-1468.

[23] Blanton SH, . Variation in IRF6 contributes to nonsyndromic cleft lip and palate. Am. J. Med. Genet. A, 2005, 137A: 259-262.

[24] Jakobsen LP, . Suggestive linkage to a neighboring region of IRF6 in a cleft lip and palate multiplex family. Am. J. Med. Genet. A, 2007, 143A: 2716-2721.

[25] Kerameddin S, Namipashaki A, Ebrahimi S, Ansari-Pour N. IRF6 is a marker of severity in nonsyndromic cleft lip/palate. J. Dent. Res., 2015, 94: 226S-232S.

[26] Creyghton MP, . Histone H3K27ac separates active from poised enhancers and predicts developmental state. Proc. Natl Acad. Sci. USA, 2010, 107: 21931-21936.

[27] Juriloff DM, Harris MJ. Mouse genetic models of cleft lip with or without cleft palate. Birth Defects Res. A Clin. Mol. Teratol., 2008, 82: 63-77.

[28] Richardson R, . p63 exerts spatio-temporal control of palatal epithelial cell fate to prevent cleft palate. PLoS Genet., 2017, 13: e1006828.

[29] Wessagowit V, Mellerio JE, Pembroke AC, McGrath JA. Heterozygous germline missense mutation in the p63 gene underlying EEC syndrome. Clin. Exp. Dermatol., 2000, 25: 441-443.

[30] Martinelli M, . Evidence of LEF1 fetal-maternal interaction in cleft lip with or without cleft palate in a consistent Italian sample study. Int. J. Immunopathol. Pharmacol., 2011, 24: 15-19.

[31] Rifai L, . Ectodermal dysplasia-like syndrome with mental retardation due to contiguous gene deletion: further clinical and molecular delineation of del(2q32) syndrome. Am. J. Med. Genet. A, 2010, 152A: 111-117.

[32] Lin D, . Digestion-ligation-only Hi-C is an efficient and cost-effective method for chromosome conformation capture. Nat. Genet., 2018, 50: 754-763.

[33] Nagano T, . Comparison of Hi-C results using in-solution versus in-nucleus ligation. Genome Biol., 2015, 16

[34] Schmitt AD, . A compendium of chromatin contact maps reveals spatially active regions in the human genome. Cell Rep., 2016, 17: 2042-2059.

[35] Smith EM, Lajoie BR, Jain G, Dekker J. Invariant TAD boundaries constrain cell-type-specific looping interactions between promoters and distal elements around the CFTR locus. Am. J. Hum. Genet., 2016, 98: 185-201.

[36] Leslie EJ, . A genome-wide association study of nonsyndromic cleft palate identifies an etiologic missense variant in GRHL3. Am. J. Hum. Genet., 2016, 98: 744-754.

[37] Ay F, Bailey TL, Noble WS. Statistical confidence estimation for Hi-C data reveals regulatory chromatin contacts. Genome Res, 2014, 24: 999-1011.

[38] Ghandi M, . gkmSVM: an R package for gapped-kmer SVM. Bioinformatics, 2016, 32: 2205-2207.

[39] Langer L, Sulik K, Pevny L. Cleft palate in a mouse model of SOX2 haploinsufficiency. Cleft Palate Craniofac. J., 2014, 51: 110-114.

[40] Sweat YY, . Sox2 controls periderm and rugae development to inhibit oral adhesions. J. Dent. Res., 2020, 99: 1397-1405.

[41] Schoenfelder S, Fraser P. Long-range enhancer-promoter contacts in gene expression control. Nat. Rev. Genet., 2019, 20: 437-455.

[42] Quillien A, . Robust identification of developmentally active endothelial enhancers in zebrafish using FANS-assisted ATAC-Seq. Cell Rep., 2017, 20: 709-720.

[43] Wilkerson BA, . Open chromatin dynamics in prosensory cells of the embryonic mouse cochlea. Sci. Rep., 2019, 9

[44] Wilderman A, VanOudenhove J, Kron J, Noonan JP, Cotney J. High-resolution epigenomic atlas of human embryonic craniofacial development. Cell Rep., 2018, 23: 1581-1597.

[45] Welzenbach J, . Integrative approaches generate insights into the architecture of non-syndromic cleft lip with or without cleft palate. HGG Adv., 2021, 2: 100038.

[46] Calo E, Wysocka J. Modification of enhancer chromatin: what, how, and why. Mol. Cell, 2013, 49: 825-837.

[47] Liu H, . Analysis of zebrafish periderm enhancers facilitates identification of a regulatory variant near human KRT8/18. eLife, 2020, 9: e51325.

[48] Rivera CM, Ren B. Mapping human epigenomes. Cell, 2013, 155: 39-55.

[49] Dixon JR, . Topological domains in mammalian genomes identified by analysis of chromatin interactions. Nature, 2012, 485: 376-380.

[50] Beiraghi S, . Identification and characterization of a novel gene disrupted by a pericentric inversion inv(4)(p13.1q21.1) in a family with cleft lip. Gene, 2003, 309: 11-21.

[51] Igal RA, Sinner DI. Stearoyl-CoA desaturase 5 (SCD5), a Delta-9 fatty acyl desaturase in search of a function. Biochim. Biophys. Acta Mol. Cell Biol. Lipids, 2021, 1866: 158840.

[52] Fan YY, . Characterization of an arachidonic acid-deficient (Fads1 knockout) mouse model. J. Lipid Res., 2012, 53: 1287-1295.

[53] Bikle DD, Xie Z, Tu CL. Calcium regulation of keratinocyte differentiation. Expert Rev. Endocrinol. Metab., 2012, 7: 461-472.

[54] Zhang Z, . Robust capturing chromosome conformation using the DLO Hi-C 2.0 method. J. Genet. Genomics, 2020, 47: 655-658.

[55] Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat. Biotechnol., 2016, 34: 525-527.

[56] Hong, P. et al. The DLO Hi-C tool for digestion-ligation-only Hi-C chromosome conformation capture data analysis. Genes (Basel) 11, https://doi.org/10.3390/genes11030289 (2020).

[57] Durand NC, . Juicer provides a one-click system for analyzing loop-resolution Hi-C experiments. Cell Syst., 2016, 3: 95-98.

[58] Miura H, Poonperm R, Takahashi S, Hiratani I. Practical analysis of Hi-C data: generating A/B compartment profiles. Methods Mol. Biol., 2018, 1861: 221-245.

[59] Kaul A, Bhattacharyya S, Ay F. Identifying statistically significant chromatin contacts from Hi-C data with FitHiC2. Nat. Protoc., 2020, 15: 991-1012.

[60] Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics, 2014, 30: 2114-2120.

[61] Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat. Methods, 2012, 9: 357-359.

[62] Li H, . The sequence alignment/map format and SAMtools. Bioinformatics, 2009, 25: 2078-2079.

[63] Zhang Y, . Model-based analysis of ChIP-Seq (MACS). Genome Biol., 2008, 9

[64] Ramirez F, . deepTools2: a next generation web server for deep-sequencing data analysis. Nucleic Acids Res., 2016, 44: W160-W165.

[65] Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics, 2010, 26: 841-842.

[66] Fisher S, . Evaluating the biological relevance of putative enhancers using Tol2 transposon-mediated transgenesis in zebrafish. Nat. Protoc., 2006, 1: 1297-1305.

[67] Liu H, . Identification of common non-coding variants at 1p22 that are functional for non-syndromic orofacial clefting. Nat. Commun., 2017, 8

Funding

National Natural Science Foundation of China (National Science Foundation of China)(82071077, 81771057, 81400477, 82181340279, 82071110, 81771066)

The Fundamental Research Funds for the Central Universities; The Young Top-notch Talent Cultivation Program of Hubei Province