Systematic analysis of copy number variants of a large cohort of orofacial cleft patients identifies candidate genes for orofacial clefts.

Orofacial clefts (OFCs) represent a large fraction of human birth defects and are one of the most common phenotypes affected by large copy number variants (CNVs). Due to the limited number of CNV patients in individual centers, CNV analyses of a large number of OFC patients are challenging. The present study analyzed 249 genomic deletions and 226 duplications from a cohort of 312 OFC patients reported in two publicly accessible databases of chromosome imbalance and phenotype in humans, DECIPHER and ECARUCA. Genomic regions deleted or duplicated in multiple patients were identified, and genes in these overlapping CNVs were prioritized based on the number of genes encompassed by the region and gene expression in embryonic mouse palate. Our analyses of these overlapping CNVs identified two genes known to be causative for human OFCs, SATB2 and MEIS2, and 12 genes (DGCR6, FGF2, FRZB, LETM1, MAPK3, SPRY1, THBS1, TSHZ1, TTC28, TULP4, WHSC1, WHSC2) that are associated with OFC or orofacial development. Additionally, we report 34 deleted and 24 duplicated genes that have not previously been associated with OFCs but are associated with the BMP, MAPK and RAC1 pathways. Statistical analyses show that the high number of overlapping CNVs is not due to random occurrence. The identified genes are not located in highly variable genomic regions in healthy populations and are significantly enriched for genes that are involved in orofacial development. In summary, we report a CNV analysis pipeline of a large cohort of OFC patients and identify novel candidate OFC genes.