MRPL53, a New Candidate Gene for Orofacial Clefting, Identified Using an eQTL Approach.

A valuable approach to understand how individual and population genetic differences can predispose to disease is to assess the impact of genetic variants on cellular functions (e.g., gene expression) of cell and tissue types related to pathological states. To understand the genetic basis of nonsyndromic cleft lip with or without cleft palate (NSCL/P) susceptibility, a complex and highly prevalent congenital malformation, we searched for genetic variants with a regulatory role in a disease-related tissue, the lip muscle (orbicularis oris muscle [OOM]), of affected individuals. From 46 OOM samples, which are frequently discarded during routine corrective surgeries on patients with orofacial clefts, we derived mesenchymal stem cells and correlated the individual genetic variants with gene expression from these cultured cells. Through this strategy, we detected significant cis-eQTLs (i.e., DNA variants affecting gene expression) and selected a few candidates to conduct an association study in a large Brazilian cohort (624 patients and 668 controls). This resulted in the discovery of a novel susceptibility locus for NSCL/P, rs1063588, the best eQTL for the MRPL53 gene, where evidence for association was mostly driven by the Native American ancestry component of our Brazilian sample. MRPL53 (2p13.1) encodes a 39S protein subunit of mitochondrial ribosomes and interacts with MYC, a transcription factor required for normal facial morphogenesis. Our study illustrates not only the importance of sampling admixed populations but also the relevance of measuring the functional effects of genetic variants over gene expression to dissect the complexity of disease phenotypes.

Impact of rare variants in ARHGAP29 to the etiology of oral clefts: role of loss-of-function vs missense variants.

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a prevalent, complex congenital malformation. Genome-wide association studies (GWAS) on NSCL/P have consistently identified association for the 1p22 region, in which ARHGAP29 has emerged as the main candidate gene. ARHGAP29 re-sequencing studies in NSCL/P patients have identified rare variants; however, their clinical impact is still unclear. In this study we identified 10 rare variants in ARHGAP29, including five missense, one in-frame deletion, and four loss-of-function (LoF) variants, in a cohort of 188 familial NSCL/P cases. A significant mutational burden was found for LoF (Sequence Kernel Association Test, p = 0.0005) but not for missense variants in ARHGAP29, suggesting that only LoF variants contribute to the etiology of NSCL/P. Penetrance was estimated as 59%, indicating that heterozygous LoF variants in ARHGAP29 confer a moderate risk to NSCL/P. The GWAS hits in IRF6 (rs642961) and 1p22 (rs560426 and rs4147811) do not seem to contribute to the penetrance of the phenotype, based on co-segregation analysis. Our data show that rare variants leading to haploinsufficiency of ARHGAP29 represent an important etiological clefting mechanism, and genetic testing for this gene might be taken into consideration in genetic counseling of familial cases.