Search

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.

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.

Orofacial clefts (OFCs) are the most common craniofacial birth defects and are often categorized into two etiologically distinct groups: cleft lip with or without cleft palate (CL/P) and isolated cleft palate (CP). CP is highly heritable, but there are still relatively few established genetic risk factors associated with its occurrence compared to CL/P. Historically, CP has been studied as a single phenotype despite manifesting across a spectrum of defects involving the hard and/or soft palate. We performed GWAS using transmission disequilibrium tests using 435 case-parent trios to evaluate broad risks for any cleft palate (ACP, n=435), as well as subtype-specific risks for any cleft soft palate (CSP, n=259) and any cleft hard palate (CHP, n=125). We identified a single genome-wide significant locus at 9q33.3 (lead SNP rs7035976, p=4.24×10 -8 ) associated with CHP. One gene at this locus, angiopoietin-like 2 ( ANGPTL2 ), plays a role in osteoblast differentiation. It is expressed in craniofacial tissue of human embryos, as well as in the developing mouse palatal shelves. We found 19 additional loci reaching suggestive significance (p<5×10 -6 ), of which only one overlapped between groups (chromosome 17q24.2, ACP and CSP). Odds ratios (ORs) for each of the 20 loci were most similar across all three groups for SNPs associated with the ACP group, but more distinct when comparing SNPs associated with either the CSP or CHP groups. We also found nominal evidence of replication (p<0.05) for 22 SNPs previously associated with cleft palate (including CL/P). Interestingly, most SNPs associated with CL/P cases were found to convey the opposite effect in those replicated in our dataset for CP only. Ours is the first study to evaluate CP risks in the context of its subtypes and we provide newly reported associations affecting the broad risk for CP as well as evidence of subtype-specific risks.

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p=6.86×10 -7 ). Of the 15 variants within SEC24D , 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p=0.0009), and created binding sites for 23 TFs, including Pax6 (p=6.12×10 -5 ) and Pax9 (p= 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study demonstrates that rare genetic variation contributes to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC., Competing Interests: Conflicts of Interest The authors have no conflicts to disclose.

Whole-exome sequencing (WES) is now a relatively straightforward process to identify causal variants in Mendelian disorders. However, the same is not true for WES in families where the inheritance patterns are less clear, and a complex etiology is suspected. Orofacial clefts (OFCs) are highly heritable birth defects with both Mendelian and complex etiologies. The phenotypic spectrum of OFCs may include overt clefts and several subclinical phenotypes, such as discontinuities in the orbicularis oris muscle (OOM) in the upper lip, velopharyngeal insufficiency (VPI), microform clefts or bifid uvulas. We hypothesize that expanding the OFC phenotype to include these phenotypes can clarify inheritance patterns in multiplex families, making them appear more Mendelian. We performed whole-exome sequencing to find rare, likely causal genetic variants in 31 multiplex OFC families, which included families with multiple individuals with OFCs and individuals with subclinical phenotypes. We identified likely causal variants in COL11A2, IRF6, KLF4, SHROOM3, SMC3, TP63 , and TBX3 in seven families. Although we did not find clear evidence supporting the subclinical phenotype hypothesis, our findings support a role for rare variants in the etiology of OFCs.