The heterogeneous genetic architectures of orofacial clefts.

Orofacial clefts (OFCs) are common, affecting 1:1000 live births. OFCs occur across a phenotypic spectrum - including cleft lip (CL), cleft lip and palate (CLP), or cleft palate (CP) - and can be further subdivided based on laterality, severity, or specific structures affected. Herein we review what is known about the genetic architecture underlying each of these subtypes, considering both shared and subtype-specific risks. While there are more known genetic similarities between CL and CLP than CP, recent research supports both shared and subtype-specific genetic risk factors within and between phenotypic classifications of OFCs. Larger sample sizes and deeper phenotyping data will be of increasing importance for the discovery of novel genetic risk factors for OFCs and various subtypes going forward.

Investigating gene functions and single-cell expression profiles of de novo variants in orofacial clefts.

Orofacial clefts (OFCs) are common congenital birth defects with various etiologies, including genetic variants. Online Mendelian Inheritance in Man (OMIM) annotated several hundred genes involving OFCs. Furthermore, several hundreds of de novo variants (DNVs) have been identified from individuals with OFCs. Some DNVs are related to known OFC genes or pathways, but there are still many DNVs whose relevance to OFC development is unknown. To explore novel gene functions and their cellular expression profiles, we focused on DNVs in genes that were not listed in OMIM. We collected 960 DNVs in 853 genes from published studies and curated these genes, based on the DNVs' deleteriousness, into 230 and 23 genes related to cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO), respectively. For comparison, we curated 178 CL/P and 277 CPO genes from OMIM. In CL/P, the pathways enriched in DNV and OMIM genes were significantly overlapped (p = 0.002). Single-cell RNA sequencing (scRNA-seq) analysis of mouse lip development revealed that both gene sets had abundant expression in the ectoderm (DNV genes: adjusted p = 0.032, OMIM genes: adjusted p < 0.0002), while only DNV genes were enriched in the endothelium (adjusted p = 0.032). Although we did not achieve significant findings using CPO gene sets, which was mainly due to the limited number of DNV genes, scRNA-seq analysis implicated various expression patterns among DNV and OMIM genes. Our results suggest that combinatory pathway and scRNA-seq data analyses are helpful for contextualizing genes in OFC development.

Genome-wide study of gene-by-sex interactions identifies risks for cleft palate.

Structural birth defects affect 3-4% of all live births and, depending on the type, tend to manifest in a sex-biased manner. Orofacial clefts (OFCs) are the most common craniofacial structural birth defects and are often divided into cleft lip with or without cleft palate (CL/P) and cleft palate only (CP). Previous studies have found sex-specific risks for CL/P, but these risks have yet to be evaluated in CP. CL/P is more common in males and CP is more frequently observed in females, so we hypothesized there would also be sex-specific differences for CP. Using a trio-based cohort, we performed sex-stratified genome-wide association studies (GWAS) based on proband sex followed by a genome-wide gene-by-sex (GxS) interaction testing. There were 13 loci significant for GxS interactions, with the top finding in LTBP1 (RR=3.37 [2.04 - 5.56], p=1.93x10 -6 ). LTBP1 plays a role in regulating TGF-B bioavailability, and knockdown in both mice and zebrafish lead to craniofacial anomalies. Further, there is evidence for differential expression of LTBP1 between males and females in both mice and humans. Therefore, we tested the association between the imputed genetically regulated gene expression of genes with significant GxS interactions and the CP phenotype. We found significant association for LTBP1 in cell cultured fibroblasts in female probands (p=0.0013) but not in males. Taken altogether, we show there are sex-specific risks for CP that are otherwise undetectable in a combined sex cohort, and LTBP1 is a candidate risk gene, particularly in females.

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

Orofacial cleft (OFC) is a common human congenital anomaly. Epithelial-specific RNA splicing regulators ESRP1 and ESRP2 regulate craniofacial morphogenesis and their disruption result in OFC in zebrafish, mouse and humans. Using esrp1/2 mutant zebrafish and murine Py2T cell line models, we functionally tested the pathogenicity of human ESRP1/2 gene variants. We found that many variants predicted by in silico methods to be pathogenic were functionally benign. Esrp1 also regulates the alternative splicing of Ctnnd1 and these genes are co-expressed in the embryonic and oral epithelium. In fact, over-expression of ctnnd1 is sufficient to rescue morphogenesis of epithelial-derived structures in esrp1/2 zebrafish mutants. Additionally, we identified 13 CTNND1 variants from genome sequencing of OFC cohorts, confirming CTNND1 as a key gene in human OFC. This work highlights the importance of functional assessment of human gene variants and demonstrates the critical requirement of Esrp-Ctnnd1 acting in the embryonic epithelium to regulate palatogenesis.

Building a growing genomic data repository for maternal and fetal health through the PING Consortium.

Background: Prenatally transmitted viruses can cause severe damage to the developing brain. There is unexplained variability in prenatal brain injury and postnatal neurodevelopmental outcomes, suggesting disease modifiers. Discordant outcomes among dizygotic twins could be explained by genetic susceptibly or protection. Among several well-recognized threats to the developing brain, Zika is a mosquito-borne, positive-stranded RNA virus that was originally isolated in Uganda and spread to cause epidemics in Africa, Asia, and the Americas. In the Americas, the virus caused congenital Zika syndrome and a multitude of neurodevelopmental disorders. As of now, there is no preventative treatment or cure for the adverse outcomes caused by prenatal Zika infection. The Prenatal Infection and Neurodevelopmental Genetics (PING) Consortium was initiated in 2016 to identify factors modulating prenatal brain injury and postnatal neurodevelopmental outcomes for Zika and other prenatal viral infections. Methods: The Consortium has pooled information from eight multi-site studies conducted at 23 research centers in six countries to build a growing clinical and genomic data repository. This repository is being mined to search for modifiers of virally induced brain injury and developmental outcomes. Multilateral partnerships include commitments with Children's National Hospital (USA), Instituto Nacional de Salud (Colombia), the Natural History of Zika Virus Infection in Gestation program (Brazil), and Zika Instituto Fernandes Figueira (Brazil), in addition to the Centers for Disease Control and Prevention and the National Institutes of Health. Discussion: Our goal in bringing together these sets of patient data was to test the hypothesis that personal and populational genetic differences affect the severity of brain injury after a prenatal viral infection and modify neurodevelopmental outcomes. We have enrolled 4,102 mothers and 3,877 infants with 3,063 biological samples and clinical data covering over 80 phenotypic fields and 5,000 variables. There were several notable challenges in bringing together cohorts enrolled in different studies, including variability in the timepoints evaluated and the collected clinical data and biospecimens. Thus far, we have performed whole exome sequencing on 1,226 participants. Here, we present the Consortium's formation and the overarching study design. We began our investigation with prenatal Zika infection with the goal of applying this knowledge to other prenatal infections and exposures that can affect brain development.