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 u...

Full description

Bibliographic Details
Main Authors: Yao Xiao, Shengbo Jiao, Miao He, Da Lin, Huanyan Zuo, Jiahao Han, Yonghua Sun, Gang Cao, Zhi Chen, Huan Liu
Format: Article
Language:English
Published: Nature Publishing Group 2022-08-01
Series:International Journal of Oral Science
Online Access:https://doi.org/10.1038/s41368-022-00194-0
_version_ 1811285861823676416
author Yao Xiao
Shengbo Jiao
Miao He
Da Lin
Huanyan Zuo
Jiahao Han
Yonghua Sun
Gang Cao
Zhi Chen
Huan Liu
author_facet Yao Xiao
Shengbo Jiao
Miao He
Da Lin
Huanyan Zuo
Jiahao Han
Yonghua Sun
Gang Cao
Zhi Chen
Huan Liu
author_sort Yao Xiao
collection DOAJ
description 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.
first_indexed 2024-04-13T02:51:13Z
format Article
id doaj.art-55e7ce681fac4bdf868bfbe37daf8017
institution Directory Open Access Journal
issn 1674-2818
2049-3169
language English
last_indexed 2024-04-13T02:51:13Z
publishDate 2022-08-01
publisher Nature Publishing Group
record_format Article
series International Journal of Oral Science
spelling doaj.art-55e7ce681fac4bdf868bfbe37daf80172022-12-22T03:05:50ZengNature Publishing GroupInternational Journal of Oral Science1674-28182049-31692022-08-0114111210.1038/s41368-022-00194-0Chromatin conformation of human oral epithelium can identify orofacial cleft missing functional variantsYao Xiao0Shengbo Jiao1Miao He2Da Lin3Huanyan Zuo4Jiahao Han5Yonghua Sun6Gang Cao7Zhi Chen8Huan Liu9The State Key Laboratory Breeding Base of Basic Science of Stomatology & Key Laboratory for Oral Biomedicine of Ministry of Education, School and Hospital of Stomatology, Wuhan UniversityState Key Laboratory of Freshwater Ecology and Biotechnology, Hubei Hongshan Laboratory, Institute of Hydrobiology, Innovation Academy for Seed Design, Chinese Academy of SciencesThe State Key Laboratory Breeding Base of Basic Science of Stomatology & Key Laboratory for Oral Biomedicine of Ministry of Education, School and Hospital of Stomatology, Wuhan UniversityState Key Laboratory of Agricultural Microbiology, Huazhong Agricultural UniversityThe State Key Laboratory Breeding Base of Basic Science of Stomatology & Key Laboratory for Oral Biomedicine of Ministry of Education, School and Hospital of Stomatology, Wuhan UniversityThe State Key Laboratory Breeding Base of Basic Science of Stomatology & Key Laboratory for Oral Biomedicine of Ministry of Education, School and Hospital of Stomatology, Wuhan UniversityState Key Laboratory of Freshwater Ecology and Biotechnology, Hubei Hongshan Laboratory, Institute of Hydrobiology, Innovation Academy for Seed Design, Chinese Academy of SciencesState Key Laboratory of Agricultural Microbiology, Huazhong Agricultural UniversityThe State Key Laboratory Breeding Base of Basic Science of Stomatology & Key Laboratory for Oral Biomedicine of Ministry of Education, School and Hospital of Stomatology, Wuhan UniversityThe State Key Laboratory Breeding Base of Basic Science of Stomatology & Key Laboratory for Oral Biomedicine of Ministry of Education, School and Hospital of Stomatology, Wuhan UniversityAbstract 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.https://doi.org/10.1038/s41368-022-00194-0
spellingShingle Yao Xiao
Shengbo Jiao
Miao He
Da Lin
Huanyan Zuo
Jiahao Han
Yonghua Sun
Gang Cao
Zhi Chen
Huan Liu
Chromatin conformation of human oral epithelium can identify orofacial cleft missing functional variants
International Journal of Oral Science
title Chromatin conformation of human oral epithelium can identify orofacial cleft missing functional variants
title_full Chromatin conformation of human oral epithelium can identify orofacial cleft missing functional variants
title_fullStr Chromatin conformation of human oral epithelium can identify orofacial cleft missing functional variants
title_full_unstemmed Chromatin conformation of human oral epithelium can identify orofacial cleft missing functional variants
title_short Chromatin conformation of human oral epithelium can identify orofacial cleft missing functional variants
title_sort chromatin conformation of human oral epithelium can identify orofacial cleft missing functional variants
url https://doi.org/10.1038/s41368-022-00194-0
work_keys_str_mv AT yaoxiao chromatinconformationofhumanoralepitheliumcanidentifyorofacialcleftmissingfunctionalvariants
AT shengbojiao chromatinconformationofhumanoralepitheliumcanidentifyorofacialcleftmissingfunctionalvariants
AT miaohe chromatinconformationofhumanoralepitheliumcanidentifyorofacialcleftmissingfunctionalvariants
AT dalin chromatinconformationofhumanoralepitheliumcanidentifyorofacialcleftmissingfunctionalvariants
AT huanyanzuo chromatinconformationofhumanoralepitheliumcanidentifyorofacialcleftmissingfunctionalvariants
AT jiahaohan chromatinconformationofhumanoralepitheliumcanidentifyorofacialcleftmissingfunctionalvariants
AT yonghuasun chromatinconformationofhumanoralepitheliumcanidentifyorofacialcleftmissingfunctionalvariants
AT gangcao chromatinconformationofhumanoralepitheliumcanidentifyorofacialcleftmissingfunctionalvariants
AT zhichen chromatinconformationofhumanoralepitheliumcanidentifyorofacialcleftmissingfunctionalvariants
AT huanliu chromatinconformationofhumanoralepitheliumcanidentifyorofacialcleftmissingfunctionalvariants