Background: Nonsyndromic cleft lip with/without cleft palate (nsCL/P) is a congenital malformation of multifactorial etiology. Research has identified >40 genome-wide significant risk loci, which explain less than 40% of nsCL/P heritability. Studies show that some of the hidden heritability is explained by rare penetrant variants.
Methods: To identify new candidate genes, we searched for highly penetrant de novo variants (DNVs) in 50 nsCL/P patient/parent-trios with a low polygenic risk for the phenotype (discovery). We prioritized DNV-carrying candidate genes from the discovery for resequencing in independent cohorts of 1010 nsCL/P patients of diverse ethnicities and 1574 population-matched controls (replication). Segregation analyses and rare variant association in the replication cohort, in combination with additional data (genome-wide association data, expression, protein-protein-interactions), were used for final prioritization.
Conclusion: In the discovery step, 60 DNVs were identified in 60 genes, including a variant in the established nsCL/P risk gene CDH1. Re-sequencing of 32 prioritized genes led to the identification of 373 rare, likely pathogenic variants. Finally, MDN1 and PAXIP1 were prioritized as top candidates. Our findings demonstrate that DNV detection, including polygenic risk score analysis, is a powerful tool for identifying nsCL/P candidate genes, which can also be applied to other multifactorial congenital malformations.Funding information: The present study was supported by the German Research Foundation (DFG)-Grants BE 3828/8-1, LU 1944/2-1, MA 2546/5-1, and LU1944/3-1.
ACKNOWLEDGMENTS: The authors thank all patients, relatives, and control individuals for their participation. We thank the German support group for individuals with cleft lip and/or palate (Wolfgang Rosenthal Gesellschaft) for assistance with recruitment.We acknowledge the invaluable assistance of all clinical, laboratory, and bioinformatic personnel. The authors thank the Next Generation Sequencing Core Facility of the Medical Faculty of the University of Bonn for sequencing the samples that were used in this study. DbGaP datasets were accessed through dbGaP accession number phs000094.v1.p1 (Supplemental Acknowledgments). Finally, the authors thank the Genome Aggregation Database (gnomAD), and all groups that provided exome and genome variant data to this resource. A full list of gnomAD contributors is provided in the gnomAD flagship paper (Karczewski et al., 2020). Open Access funding enabled and organized by Projekt DEAL
