Abstract:Background: Disturbances in the intricate processes that control craniofacial morphogenesis can result in birth defects, most common of which are orofacial clefts (OFCs). Non-syndromic cleft lip (nsCL), one of the phenotypic forms amongst OFCs, has a non-random laterality presentation with the left-side being affected twice as often compared to the right-side. This study investigates the etiology of nsCL and the factors contributing to its laterality using a pair of monozygotic twins with mirror-image cleft lip. Methods: We conducted whole-genome sequencing (WGS) analyses in a female twin pair with mirror image nsCL, their affected mother and unaffected father to identify etiopathogenic variants. Additionally, to identify possible cleft lip laterality modifiers, DNA-methylome analysis was conducted to test for differential methylation patterns between the mirror twins. Lastly, DNA methylation patterns were also analyzed on an independent cohort of female cases with unilateral cleft lip (left=22; right=17) for replication purposes. Results: We identified a protein-altering variant in FGF20 (p.Ile79Val) within the fibroblast growth factor interacting family domain segregating with the nsCL in this family. Concurrently, DNA-methylome analysis identified differential methylation regions (DMRs) upstream of Zinc-finger transcription factor ZFP57 (Δβ > 5%). Replication of these results on an independent cohort, confirmed these DMRs, emphasizing their biological significance (p<0.05). Enrichment analysis indicated that these DMRs are involved in DNA methylation during early embryo development (FDR adjusted p-value = 1.3241E-13). Further bioinformatics analyses showed one of these DMRs acting as a binding site for transcription factor AP2A (TFAP2A), a key player in craniofacial development. Interactome analysis also suggested a potential role for ZFP57 in left/right axis specification, thus emphasizing its significance in cleft laterality. Conclusion: This study provides novel insights into the etiology of nsCL and its laterality, suggesting an interplay between etiopathogenic variants and DNA methylation in cleft laterality. Our findings elucidate the intricate mechanisms underlying OFCs development. Understanding these factors may offer new tools for prevention and management of OFCs, alleviating the burden on affected individuals, their families and global health.

Two non-coding variants associated with isolated orofacial cleft promote binding of transcriptional repressors FOXE1 or ETS2 and reduce expression of IRF6

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

Functional analysis of ESRP1/2 gene variants and CTNND1 isoforms in orofacial cleft pathogenesis

Testing Reported Associations of Gene Variants with Non-Syndromic Orofacial Clefts in the Polish Population

Association Between Interferon Regulatory Factor 6 Gene Polymorphisms and Nonsyndromic Cleft Lip with or Without Cleft Palate in a Chinese Population.

Novel Evidence of Association with Nonsyndromic Cleft Lip with or Without Cleft Palate Was Shown for Single Nucleotide Polymorphisms in FOXF2 Gene in an Asian Population.

Shared genetic risk between major orofacial cleft phenotypes in an African population

Case report and functional verification of a novel mutation in the interferon regulatory transcription factor 6 gene in a family with orofacial clefts

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

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

A Multi-Ethnic Genome-Wide Association Study Identifies Novel Loci for Non-Syndromic Cleft Lip with or Without Cleft Palate on 2p24.2, 17q23 and 19Q13

Genome Wide Study of Maternal and Parent-of-origin Effects on the Etiology of Orofacial Clefts.

Association Between Polymorphism of IRF6 Rs2235371 Locus and Nonsyndromic Cleft Lip with or Without Cleft Palate

Identification of rare loss-of-function variants in FAM3B associated with non-syndromic orofacial clefts

High incidence and geographic distribution of cleft palate in Finland are associated with the IRF6 gene

Genetic Variants in METTL16 Affect the Risk of Non‐Syndromic Orofacial Clefts

Multi-omic analyses of a twin pair with mirror image cleft lip identifies pathogenic variant in modified by differential methylation upstream of

[Relationship Between Genetic Polymorphisms of IRF6 Rs642961 and Nonsysdromic Cleft Lip with or Without Cleft Palate].

High incidence and geographic distribution of cleft palate cases in Finland are associated with a regulatory variant in

Association between IRF6, TP63, GREM1 Gene Polymorphisms and Non-Syndromic Orofacial Cleft Phenotypes in Vietnamese Population: A Case–Control and Family-Based Study