Gene-environment interactions: aligning birth defects research with complex etiology

TRIO RVEMVS: A Bayesian framework for rare variant association analysis with expectation-maximization variable selection using family trio data

It is commonly reported that rare variants may be more functionally related to complex diseases than common variants. However, individual rare variant association tests remain challenging due to low minor allele frequency in the available samples. This paper proposes an expectation maximization variable selection (EMVS) method to simultaneously detect common and rare variants at the individual variant level using family trio data. TRIO_RVEMVS was assessed in both large (1500 families) and small (350 families) datasets based on simulation. The performance of TRIO_RVEMVS was compared with gene-level kernel and burden association tests that use pedigree data (PedGene) and rare-variant extensions of the transmission disequilibrium test (RV-TDT). At the region level, TRIO_RVEMVS outperformed PedGene and RV-TDT when common variants were included. TRIO_RVEMVS performed competitively with PedGene and outperformed RV-TDT when the analysis was only restricted to rare variants. At the individual variants level, with 1,500 trios, the average true positive rate of individual rare variants that were polymorphic across 500 datasets was 12.20%, and the average false positive rate was 0.74%. In the datasets with 350 trios, the average true and false positive rates of individual rare variants were 13.10% and 1.30%, respectively. When applying TRIO_RVEMVS to real data from the Gabriella Miller Kids First Pediatric Research Program, it identified 3 rare variants in q24.21 and q24.22 associated with the risk of orofacial clefts in the Kids First European population.

Examination of piperonyl butoxide developmental toxicity as a Sonic hedgehog pathway inhibitor targeting limb and palate morphogenesis

Piperonyl butoxide (PBO) is a pesticide synergist with widespread use and human exposure that was discovered to inhibit Sonic hedgehog (Shh) signaling, a pathway required for numerous developmental processes. Previous examinations of PBO's potential for developmental toxicity have generated seemingly conflicting results. We investigated the impact of acute PBO exposure targeting Shh pathway activity during palate and limb morphogenesis. Timed-pregnant C57BL/6 J mice were exposed to a single PBO dose (67-1800 mg/kg) at gestational day (GD) 9.75, and litters were collected at GD10.25 and GD10.75 to examine Shh pathway activity or GD17 for phenotypic assessment. PBO exposure induced dose-dependent limb malformations and cleft palate in the highest dose group. Following PBO exposure, reduced expression of the Shh pathway activity markers Gli1 and Ptch1 was observed in the embryonic limb buds and craniofacial processes. These findings provide additional evidence that prenatal PBO exposure targeting Shh pathway activity can result in malformations in mice that parallel common etiologically complex human birth defects.

Use of an Orthodontic and Otolaryngological Approach in an Infant with Holoprosencephaly

Holoprosencephaly is a complex human brain malformation resulting from incomplete cleavage of the prosencephalon into both hemispheres. Congenital nasal pyriform aperture stenosis (CNPAS) is sometimes found in patients with mild forms of holoprosencephaly. Surgical treatment is required. Low-invasive surgical approaches involve balloon dilation of the pyriform opening. We present the case of an 8-day-old girl diagnosed with holoprosencephaly, CNPAS, and the presence of a solitary median maxillary central incisor. Once examined by neonatologist, geneticist, pneumologist, otolaryngologist, and pediatric dentist, a combined otolaryngological-orthodontic approach was used. The obstruction of the right nasal cavity was treated by widening the nasal cavities and stabilizing them with a balloon dilation technique. After surgery, the respiratory space was increased by applying a neonatal palatal expander plate (NPEP) considering the palatal deformity: ogival shaped, anterior vertex growth direction, reduction of transverse diameters. The NPEP promoted distraction of the median palatine suture and assisted the nasal dilation. Therefore, after the insertion of NPEP, the physiological sucking-swallowing mechanism was activated. In infants with CNPAS, NPEP can be useful to ensure the safe stability of nasal dilation. A multidisciplinary approach is fundamental. In our experience, the close collaboration between an otolaryngologist and orthodontist is essential for the management of the patient with CNPAS.

Organ transformation by environmental disruption of protein integrity and epigenetic memory in Drosophila

Despite significant progress in understanding epigenetic reprogramming of cells, the mechanistic basis of "organ reprogramming" by (epi-)gene-environment interactions remained largely obscure. Here, we use the ether-induced haltere-to-wing transformations in Drosophila as a model for epigenetic "reprogramming" at the whole organism level. Our findings support a mechanistic chain of events explaining why and how brief embryonic exposure to ether leads to haltere-to-wing transformations manifested at the larval stage and on. We show that ether interferes with protein integrity in the egg, leading to altered deployment of Hsp90 and widespread repression of Trithorax-mediated establishment of active H3K4me3 chromatin marks throughout the genome. Despite this global reduction, Ubx targets and wing development genes preferentially retain higher levels of H3K4me3 that predispose these genes for later up-regulation in the larval haltere disc, hence the wing-like outcome. Consistent with compromised protein integrity during the exposure, the penetrance of bithorax transformations increases by genetic or chemical reduction of Hsp90 function. Moreover, joint reduction in Hsp90 and trx gene dosage can cause bithorax transformations without exposure to ether, supporting an underlying epistasis between Hsp90 and trx loss-of-functions. These findings implicate environmental disruption of protein integrity at the onset of histone methylation with altered epigenetic regulation of developmental patterning genes. The emerging picture provides a unique example wherein the alleviation of the Hsp90 "capacitor function" by the environment drives a morphogenetic shift towards an ancestral-like body plan. The morphogenetic impact of chaperone response during a major setup of epigenetic patterns may be a general scheme for organ transformation by environmental cues.

Genetic Analysis and Functional Assessment of a TGFBR2 Variant in Micrognathia and Cleft Palate

Cleft lip and cleft palate are among the most common congenital anomalies and are the result of incomplete fusion of embryonic craniofacial processes or palatal shelves, respectively. We know that genetics play a large role in these anomalies but the list of known causal genes is far from complete. As part of a larger sequencing effort of patients with micrognathia and cleft palate we identified a candidate variant in transforming growth factor beta receptor 2 (TGFBR2) which is rare, changing a highly conserved amino acid, and predicted to be pathogenic by a number of metrics. The family history and population genetics would suggest this specific variant would be incompletely penetrant, but this gene has been convincingly implicated in craniofacial development. In order to test the hypothesis this might be a causal variant, we used genome editing to create the orthologous variant in a new mouse model. Surprisingly, Tgfbr2V387M mice did not exhibit craniofacial anomalies or have reduced survival suggesting this is, in fact, not a causal variant for cleft palate/ micrognathia. The discrepancy between in silico predictions and mouse phenotypes highlights the complexity of translating human genetic findings to mouse models. We expect these findings will aid in interpretation of future variants seen in TGFBR2 from ongoing sequencing of patients with congenital craniofacial anomalies.

The detection efficacy of noninvasive prenatal genetic testing (NIPT) for sex chromosome abnormalities and copy number variation and its differentiation in pregnant women of different ages

Objective

To analyze the efficacy of noninvasive prenatal genetic testing (NIPT) in detecting fetal sex chromosome abnormalities and copy number variation (CNV), compare the efficacy between NIPT and serological screening alone, and further analyze the fetal sex chromosome abnormalities and CNV differentiation in pregnant women of different ages, so as to provide a reference for the prevention and control of fetal birth defects.

Methods

Clinical data from 22,692 pregnant women admitted to our hospital from January 2013 to December 2022 were retrospectively analyzed. All participants underwent serological screening and NIPT screening to compare fetal chromosomal abnormalities between the two screening modalities. 145 women whose fetus were diagnosed as sex chromosome abnormalities and 36 women whose fetus were diagnosed as CNV abnormalities based on NIPT screening were selected for prenatal diagnosis by amniocentesis or karyotyping. Taking prenatal diagnosis as the standard, the four-grid table method was used to detect the positive predictive value of NIPT screening for fetal sex chromosomal abnormalities and CNV. According to the age, pregnant women were divided into 18-30 years old (n = 9844), 31-35 years old (n = 7612), >35 years old (n = 5236), and then the detection rates of sexual fetal chromosomal abnormalities, CNV and total chromosomal abnormalities were compared in pregnant women.

Results

Among the 22,692 pregnant women in this study, the high-risk proportion of serologic screening with 4.38% was higher than that of NIPT screening with 1.93% (P < 0.05). Among the 145 women with fetal sex chromosome abnormalities screened by NIPT, 122 cases of fetal sex chromosome abnormalities were diagnosed prenatally, including 45, X/47, XXX/47, XYY/47, XXY. The positive predictive values of NIPT screening were 25.00%, 58.82%, 85.71%, and 85.71%, respectively, with an overall predictive value of 44.26%. The positive predictive value of fetal sex chromosome abnormalities in NIPT screening was higher than that of serological screening (P < 0.05). Among the 36 pregnant women with fetal CNV, NIPT screening showed that CNVs≤10 Mb and CNVs>10 Mb were 33.33% and 66.67%, respectively. There were 12 cases of prenatal diagnosis of fetal CNV, among which the NIPT-screened positive predictive values of fetal copy number deletion, duplicate, deletion and duplicate were 50.00%, 57.14% and 100.00%, respectively, with an overall predictive value of 58.33%. The positive predictive value of CNV in NIPT screening was higher than that of serological screening without statistically significant difference (P > 0.05). The results of NIPT screening showed that the detection rate of fetal sex chromosome abnormalities and total abnormalities of pregnant women over 35 years of age was significantly higher than that of pregnant women aged 18-30 and 31-35 years (P < 0.05).

Conclusion

NIPT screening could greatly improve the detection efficacy of fetal sex chromosome abnormalities, CNV and other chromosome abnormalities, and decline the false positive rate. However, the positive predictive value of NIPT screening was relatively low, and further prenatal testing and genetic counseling are still required. In addition, NIPT screening for fetal sex chromosome abnormalities, and the detection rate of total abnormalities in pregnant women older than 35 years old were increased significantly, and pregnancy at an advanced age may be one of the risk factors for fetal chromosomal abnormalities.

Rare variants found in multiplex families with orofacial clefts: Does expanding the phenotype make a difference?

Exome sequencing (ES) is now a relatively straightforward process to identify causal variants in Mendelian disorders. However, the same is not true for ES 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 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, SHROOM3, SMC3, TBX3, and TP63 in six 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.

Residential proximity to unconventional oil and gas development and birth defects in Ohio

BACKGROUND

Chemicals used or emitted by unconventional oil and gas development (UOGD) include reproductive/developmental toxicants. Associations between UOGD and certain birth defects were reported in a few studies, with none conducted in Ohio, which experienced a thirty-fold increase in natural gas production between 2010 and 2020.

METHODS

We conducted a registry-based cohort study of 965,236 live births in Ohio from 2010 to 2017. Birth defects were identified in 4653 individuals using state birth records and a state surveillance system. We assigned UOGD exposure based on maternal residential proximity at birth to active UOG wells and a metric specific to the drinking-water exposure pathway that identified UOG wells hydrologically connected to a residence ("upgradient UOG wells"). We estimated odds ratios (ORs) and 95% confidence intervals (CIs) for all structural birth defects combined and specific birth defect types using binary exposure metrics (presence/absence of any UOG well and presence/absence of an upgradient UOG well within 10 km), adjusting for confounders. Additionally, we conducted analyses stratified by urbanicity, infant sex, and social vulnerability.

RESULTS

The odds of any structural defect were 1.13 times higher in children born to mothers living within 10 km of UOGD than those born to unexposed mothers (95%CI: 0.98-1.30). Odds were elevated for neural tube defects (OR: 1.57, 95%CI: 1.12-2.19), limb reduction defects (OR: 1.99, 95%CI: 1.18-3.35), and spina bifida (OR 1.93; 95%CI 1.25-2.98). Hypospadias (males only) was inversely related to UOGD exposure (OR: 0.62, 95%CI: 0.43-0.91). Odds of any structural defect were greater in magnitude but less precise in analyses using the hydrological-specific metric (OR: 1.30; 95%CI: 0.85-1.90), in areas with high social vulnerability (OR: 1.27, 95%CI: 0.99-1.60), and among female offspring (OR: 1.28, 95%CI: 1.06-1.53).

CONCLUSIONS

Our results suggest a positive association between UOGD and certain birth defects, and findings for neural tube defects corroborate results from prior studies.

Rare variants found in multiplex families with orofacial clefts: Does expanding the phenotype make a difference?

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.

Exposure to an environmentally relevant concentration of 17α-ethinylestradiol disrupts craniofacial development of juvenile zebrafish

Endocrine disrupting chemicals (EDCs) can interact with native hormone receptors to interfere with and disrupt hormone signalling that is necessary for a broad range of developmental pathways. EDCs are pervasive in our environment, in particular in our waterways, making aquatic wildlife especially vulnerable to their effects. Many of these EDCs are able to bind to and activate oestrogen receptors, causing aberrant oestrogen signalling. Craniofacial development is an oestrogen-sensitive process, with oestrogen receptors expressed in chondrocytes during critical periods of development. Previous studies have demonstrated a negative effect of high concentrations of oestrogen on early craniofacial patterning in the aquatic model organism, the zebrafish (Danio rerio). In order to determine the impacts of exposure to an oestrogenic EDC, we exposed zebrafish larvae and juveniles to either a high concentration to replicate previous studies, or a low, environmentally relevant concentration of the oestrogenic contaminant, 17α-ethinylestradiol. The prolonged / chronic exposure regimen was used to replicate that seen by many animals in natural waterways. We observed changes to craniofacial morphology in all treatments, and most strikingly in the larvae-juveniles exposed to a low concentration of EE2. In the present study, we have demonstrated that the developmental stage at which exposure occurs can greatly impact phenotypic outcomes, and these results allow us to understand the widespread impact of oestrogenic endocrine disruptors. Given the conservation of key craniofacial development pathways across vertebrates, our model can further be applied in defining the risks of EDCs on mammalian organisms.

High-throughput detection of craniofacial defects in fluorescent zebrafish

Losses and malformations of cranial neural crest cell (cNCC) derivatives are a hallmark of several common brain and face malformations. Nevertheless, the etiology of these cNCC defects remains unknown for many cases, suggesting a complex basis involving interactions between genetic and/or environmental factors. However, the sheer number of possible factors (thousands of genes and hundreds of thousands of toxicants) has hindered identification of specific interactions. Here, we develop a high-throughput analysis that will enable faster identification of multifactorial interactions in the genesis of craniofacial defects. Zebrafish embryos expressing a fluorescent marker of cNCCs (fli1:EGFP) were exposed to a pathway inhibitor standard or environmental toxicant, and resulting changes in fluorescence were measured in high-throughput using a fluorescent microplate reader to approximate cNCC losses. Embryos exposed to the environmental Hedgehog pathway inhibitor piperonyl butoxide (PBO), a Hedgehog pathway inhibitor standard, or alcohol (ethanol) exhibited reduced fli1:EGFP fluorescence at one day post fertilization, which corresponded with craniofacial defects at five days post fertilization. Combining PBO and alcohol in a co-exposure paradigm synergistically reduced fluorescence, demonstrating a multifactorial interaction. Using pathway reporter transgenics, we show that the plate reader assay is sensitive at detecting alterations in Hedgehog signaling, a critical regulator of craniofacial development. We go on to demonstrate that this technique readily detects defects in other important cell types, namely neurons. Together, these findings demonstrate this novel in vivo platform can predict developmental abnormalities and multifactorial interactions in high-throughput.

Gene-environment interactions in birth defect etiology: Challenges and opportunities

Birth defects are relatively common congenital outcomes that significantly impact affected individuals, their families, and communities. Effective development and deployment of prevention and therapeutic strategies for these conditions requires sufficient understanding of etiology, including underlying genetic and environmental causes. Tremendous progress has been made in defining the genetic basis of familial and syndromic forms of birth defects. However, the majority of birth defect cases are considered nonsyndromic and thought to result from multifactorial gene-environment interactions. While substantial advances have been made in elucidating the genetic landscape of these etiologically complex conditions, significant biological and technical constraints have stymied progress toward a refined knowledge of environmental risk factors. Defining specific gene-environment interactions in birth defect etiology is even more challenging. However, progress has been made, including demonstration of critical proofs of concept and development of new conceptual and technical approaches for resolving complex gene-environment interactions. In this review, we discuss current views of multifactorial birth defect etiology, comparing them with other diseases that also involve gene-environment interactions, including primary immunodeficiency and cancer. We describe how various model systems have illuminated mechanisms of multifactorial etiology and these models' individual strengths and weaknesses. Finally, suggestions for areas of future emphasis are proposed.

Gene-environment interactions in the pathogenesis of common craniofacial anomalies

Craniofacial anomalies often exhibit phenotype variability and non-mendelian inheritance due to their multifactorial origin, involving both genetic and environmental factors. A combination of epidemiologic studies, genome-wide association, and analysis of animal models have provided insight into the effects of gene-environment interactions on craniofacial and brain development and the pathogenesis of congenital disorders. In this chapter, we briefly summarize the etiology and pathogenesis of common craniofacial anomalies, focusing on orofacial clefts, hemifacial microsomia, and microcephaly. We then discuss how environmental risk factors interact with genes to modulate the incidence and phenotype severity of craniofacial anomalies. Identifying environmental risk factors and dissecting their interaction with different genes and modifiers is central to improved strategies for preventing craniofacial anomalies.

Gene-alcohol interactions in birth defects

Most human birth defects are thought to result from complex interactions between combinations of genetic and environmental factors. This is true even for conditions that, at face value, may appear simple and straightforward, like fetal alcohol spectrum disorders (FASD). FASD describe the full range of structural and neurological disruptions that result from prenatal alcohol exposure. While FASD require alcohol exposure, evidence from human and animal model studies demonstrate that additional genetic and/or environmental factors can influence the embryo's susceptibility to alcohol. Only a limited number of alcohol interactions in birth defects have been identified, with many sensitizing genetic and environmental factors likely yet to be identified. Because of this, while unsatisfying, there is no definitively "safe" dose of alcohol for all pregnancies. Determining these other factors, as well as mechanistically characterizing known interactions, is critical for better understanding and preventing FASD and requires combined scrutiny of human and model organism studies.

Brain Organization and Human Diseases

The cortex is a highly organized structure that develops from the caudal regions of the segmented neural tube. Its spatial organization sets the stage for future functional arealization. Here, we suggest using a developmental perspective to describe and understand the etiology of common cortical malformations and their manifestation in the human brain.

Subphenotypes in Non-Syndromic Orofacial Cleft Patients Based on the Tooth Agenesis Code (TAC)

Background: It was the aim of this study to investigate tooth agenesis patterns, which are expressed to different subphenotypes according to the TAC method in a spectrum of non-syndromic orofacial cleft patients. Methods: A total of 183 orofacial cleft patient records were assessed for tooth agenesis and TAC patterns. The association between TAC and sex, and cleft type was examined, and logistic regression models were additionally applied. Additionally, the distribution of missing teeth by cleft type and the tooth agenesis inter-quadrant association were examined. Results: The most frequent cleft type was CLPL (n = 72; 39.3%), while the maxillary left lateral incisor was the most frequently missing tooth that was strongly dependent on the cleft type (29.5%, p < 0.001). Of the 31 TAC patterns identified, four were the most prevalent and occurred in 80.8% of the sample, while 20 TAC patterns were unique. Cleft type contrary to sex (p = 0.405) was found to play a significant role in TAC distribution (p = 0.001). The logistic regression’s results suggested that overall, neither sex nor cleft type were associated with tooth agenesis. Prevalence of tooth agenesis in each quadrant clearly depended on cleft type; and there was a strong association found between tooth agenesis in different quadrants. Conclusions: Thirty-one different subphenotypes were identified in TAC patterns. The first four TAC patterns accounted for the 80.8% of the sample’s variability while twenty of the patterns were unique. A strong association was present between TAC pattern and cleft type. No association was found between the sex of the patient, tooth agenesis and TAC patterns. Tooth agenesis depended strongly on the cleft type, and the most frequently missing tooth was the maxillary left lateral incisor. The interquadrant association for tooth agenesis found suggests a genetic link in the etiology of clefts.

Analysis of Gene-Environment Interactions Related to Developmental Disorders

Various genetic and environmental factors are associated with developmental disorders (DDs). It has been suggested that interaction between genetic and environmental factors (G × E) is involved in the etiology of DDs. There are two major approaches to analyze the interaction: genome-wide and candidate gene-based approaches. In this mini-review, we demonstrate how these approaches can be applied to reveal the G × E related to DDs focusing on zebrafish and mouse models. We also discuss novel approaches to analyze the G × E associated with DDs.

Genome-wide Interaction Study Implicates VGLL2 and Alcohol Exposure and PRL and Smoking in Orofacial Cleft Risk

Non-syndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect, affecting approximately 1 in 700 births. NSCL/P has complex etiology including several known genes and environmental factors; however, known genetic risk variants only account for a small fraction of the heritability of NSCL/P. It is commonly suggested that gene-by-environment (G×E) interactions may help explain some of the “missing” heritability of NSCL/P. We conducted a genome-wide G×E interaction study in cases and controls of European ancestry with three common maternal exposures during pregnancy: alcohol, smoking, and vitamin use using a two-stage design. After selecting 127 loci with suggestive 2df tests for gene and G x E effects, 40 loci showed significant G x E effects after correcting for multiple tests. Notable interactions included SNPs of 6q22 near VGLL2 with alcohol and 6p22.3 near PRL with smoking. These interactions could provide new insights into the etiology of CL/P and new opportunities to modify risk through behavioral changes.

Concepts in Multifactorial Etiology of Developmental Disorders: Gene-Gene and Gene-Environment Interactions in Holoprosencephaly

Many common developmental disorders are thought to arise from a complex set of genetic and environmental risk factors. These factors interact with each other to affect the strength and duration of key developmental signaling pathways, thereby increasing the possibility that they fail to achieve the thresholds required for normal embryonic patterning. One such disorder, holoprosencephaly (HPE), serves as a useful model system in understanding various forms of multifactorial etiology. Genomic analysis of HPE cases, epidemiology, and mechanistic studies of animal models have illuminated multiple potential ways that risk factors interact to produce adverse developmental outcomes. Among these are: 1) interactions between driver and modifier genes; 2) oligogenic inheritance, wherein each parent provides predisposing variants in one or multiple distinct loci; 3) interactions between genetic susceptibilities and environmental risk factors that may be insufficient on their own; and 4) interactions of multiple genetic variants with multiple non-genetic risk factors. These studies combine to provide concepts that illuminate HPE and are also applicable to additional disorders with complex etiology, including neural tube defects, congenital heart defects, and oro-facial clefting.

Engineering Epithelial-Mesenchymal Microtissues to Study Cell-Cell Interactions in Development

Intercellular signaling drives human development, but there is a paucity of in vitro models that recapitulate important tissue architecture while remaining operationally simple and scalable. As an example, formation of the upper lip and palate requires the orchestrated proliferation and fusion of embryonic facial growth centers and is dependent on paracrine epithelial-mesenchymal signaling through multiple pathways including the Sonic Hedgehog (SHH), transforming growth factor-beta (Tgf-β), bone morphogenic protein (BMP), and epidermal growth factor (EGF) pathways. We have developed a robust, throughput-compatible microphysiological system to model intercellular signaling including epithelial-mesenchymal interactions that is useful for studying both normal and abnormal orofacial development. We describe the construction and operation of an engineered microplate created using CNC micromilling of 96-well microtiter plates capable of containing up to 20 epithelial-mesenchymal microtissues. A dense three-dimensional mesenchyme is created by embedding cells (O9-1, 3T3) in a biomimetic hydrogel. An epithelial layer is then overlayed on the microtissue by loading cells in engineered microchannels that flank the microtissue. The result is an engineering epithelial-mesenchymal interface that is both on and perpendicular to the imaging plane making it suitable for high-content imaging and analysis. The resulting microtissues and device are compatible with diverse analytical techniques including fluorescent and luminescent cell health and enzymatic reporter assays, gene expression analyses, and protein staining. This tractable model and approach promise to shed light on critical processes in intercellular signaling events in orofacial development and beyond.

Systems biology analysis of human genomes points to key pathways conferring spina bifida risk

Spina bifida (SB) is a debilitating birth defect caused by multiple gene and environment interactions. Though SB shows non-Mendelian inheritance, genetic factors contribute to an estimated 70% of cases. Nevertheless, identifying human mutations conferring SB risk is challenging due to its relative rarity, genetic heterogeneity, incomplete penetrance, and environmental influences that hamper genome-wide association studies approaches to untargeted discovery. Thus, SB genetic studies may suffer from population substructure and/or selection bias introduced by typical candidate gene searches. We report a population based, ancestry-matched whole-genome sequence analysis of SB genetic predisposition using a systems biology strategy to interrogate 298 case-control subject genomes (149 pairs). Genes that were enriched in likely gene disrupting (LGD), rare protein-coding variants were subjected to machine learning analysis to identify genes in which LGD variants occur with a different frequency in cases versus controls and so discriminate between these groups. Those genes with high discriminatory potential for SB significantly enriched pathways pertaining to carbon metabolism, inflammation, innate immunity, cytoskeletal regulation, and essential transcriptional regulation consistent with their having impact on the pathogenesis of human SB. Additionally, an interrogation of conserved noncoding sequences identified robust variant enrichment in regulatory regions of several transcription factors critical to embryonic development. This genome-wide perspective offers an effective approach to the interrogation of coding and noncoding sequence variant contributions to rare complex genetic disorders.

Genetic models and approaches to study orofacial clefts

INTRODUCTION

Orofacial clefts (OFCs) are common craniofacial birth defects with heterogeneous phenotype and etiology. Geneticists have applied nearly every available method and technology to further our understanding of the genetic architectures of OFCs.

OBJECTIVE

This review describes the evidence for a genetic etiology in OFCs, statistical genetic approaches employed to identify genetic causes, and how the results have shaped our current understanding of the genetic architectures of syndromic and nonsyndromic OFCs.

CONCLUSION

There has been rapid progress towards elucidating the genetic architectures of OFCs due to the availability of large collections of DNA samples from cases, controls, and families with OFCs and the consistent adoption of new methodologies and novel statistical approaches as they are developed. Genetic studies have identified rare and common variants influencing risk of OFCs in both Mendelian and complex forms of OFCs, blurring the distinctions traditional categories used in genetic studies and clinical medicine.

Δ9-Tetrahydrocannabinol inhibits Hedgehog-dependent patterning during development

Many developmental disorders are thought to arise from an interaction between genetic and environmental risk factors. The Hedgehog (HH) signaling pathway regulates myriad developmental processes, and pathway inhibition is associated with birth defects, including holoprosencephaly (HPE). Cannabinoids are HH pathway inhibitors, but little is known of their effects on HH-dependent processes in mammalian embryos, and their mechanism of action is unclear. We report that the psychoactive cannabinoid Δ9-tetrahydrocannabinol (THC) induces two hallmark HH loss-of-function phenotypes (HPE and ventral neural tube patterning defects) in Cdon mutant mice, which have a subthreshold deficit in HH signaling. THC therefore acts as a 'conditional teratogen', dependent on a complementary but insufficient genetic insult. In vitro findings indicate that THC is a direct inhibitor of the essential HH signal transducer smoothened. The canonical THC receptor, cannabinoid receptor-type 1, is not required for THC to inhibit HH signaling. Cannabis consumption during pregnancy may contribute to a combination of risk factors underlying specific developmental disorders. These findings therefore have significant public health relevance.

Gene-teratogen interactions influence the penetrance of birth defects by altering Hedgehog signaling strength

Birth defects result from interactions between genetic and environmental factors, but the mechanisms remain poorly understood. We find that mutations and teratogens interact in predictable ways to cause birth defects by changing target cell sensitivity to Hedgehog (Hh) ligands. These interactions converge on a membrane protein complex, the MMM complex, that promotes degradation of the Hh transducer Smoothened (SMO). Deficiency of the MMM component MOSMO results in elevated SMO and increased Hh signaling, causing multiple birth defects. In utero exposure to a teratogen that directly inhibits SMO reduces the penetrance and expressivity of birth defects in Mosmo-/- embryos. Additionally, tissues that develop normally in Mosmo-/- embryos are refractory to the teratogen. Thus, changes in the abundance of the protein target of a teratogen can change birth defect outcomes by quantitative shifts in Hh signaling. Consequently, small molecules that re-calibrate signaling strength could be harnessed to rescue structural birth defects.

The Role of Sonic Hedgehog in Human Holoprosencephaly and Short-Rib Polydactyly Syndromes

The Hedgehog (HH) signalling pathway is one of the major pathways controlling cell differentiation and proliferation during human development. This pathway is complex, with HH function influenced by inhibitors, promotors, interactions with other signalling pathways, and non-genetic and cellular factors. Many aspects of this pathway are not yet clarified. The main features of Sonic Hedgehog (SHH) signalling are discussed in relation to its function in human development. The possible role of SHH will be considered using examples of holoprosencephaly and short-rib polydactyly (SRP) syndromes. In these syndromes, there is wide variability in phenotype even with the same genetic mutation, so that other factors must influence the outcome. SHH mutations were the first identified genetic causes of holoprosencephaly, but many other genes and environmental factors can cause malformations in the holoprosencephaly spectrum. Many patients with SRP have genetic defects affecting primary cilia, structures found on most mammalian cells which are thought to be necessary for canonical HH signal transduction. Although SHH signalling is affected in both these genetic conditions, there is little overlap in phenotype. Possible explanations will be canvassed, using data from published human and animal studies. Implications for the understanding of SHH signalling in humans will be discussed.

Symmetry and fluctuation of cell movements in neural crest-derived facial mesenchyme

In the face, symmetry is established when bilateral streams of neural crest cells leave the neural tube at the same time, follow identical migration routes and then give rise to the facial prominences. However, developmental instability exists, particularly surrounding the steps of lip fusion. The causes of instability are unknown but inability to cope with developmental fluctuations are a likely cause of congenital malformations, such as non-syndromic orofacial clefts. Here, we tracked cell movements over time in the frontonasal mass, which forms the facial midline and participates in lip fusion, using live-cell imaging of chick embryos. Our mathematical examination of cell velocity vectors uncovered temporal fluctuations in several parameters, including order/disorder, symmetry/asymmetry and divergence/convergence. We found that treatment with a Rho GTPase inhibitor completely disrupted the temporal fluctuations in all measures and blocked morphogenesis. Thus, we discovered that genetic control of symmetry extends to mesenchymal cell movements and that these movements are of the type that could be perturbed in asymmetrical malformations, such as non-syndromic cleft lip.

This article has an associated ‘The people behind the papers’ interview.

Highlighted Article: Live imaging of the chick embryo face followed by mathematical analysis of mesenchymal cell tracks captures novel fluctuations between states of order/disorder as well as symmetry/asymmetry, revealing developmental instabilities that are part of normal morphogenesis.

Genetic background but not prostatic epithelial beta-catenin influences susceptibility of male mice to testosterone and estradiol-induced urinary dysfunction

Urinary voiding dysfunction in aging men can cause bothersome symptoms and irreparable tissue damage. Underlying mechanisms are not fully known. We previously demonstrated that subcutaneous, slow-release testosterone and estradiol implants (T+E2) drive a pattern of urinary voiding dysfunction in male mice that resembles that of aging men. The initial goal of this study was to test the hypothesis that prostatic epithelial beta-catenin (Ctnnb1) is required for T+E2-mediated voiding dysfunction. Targeted Ctnnb1 deletion did not significantly change voiding function in control or T+E2 treated mice but led to the surprising discovery that the C57BL/6J × FVB/NJ × 129S1 mixed genetic background onto which Ctnnb1 loss of function alleles were maintained is profoundly susceptible to voiding dysfunction. The mixed background mice develop a more rapid T+E2-mediated increase in spontaneous urine spotting, are more impaired in ability to initiate bladder contraction, and develop larger and heavier bladders than T+E2 treated C57BL/6J pure bred mice. To better understand mechanisms, we separately evaluated contributions of T and E2 and found that E2 mediates voiding dysfunction. Our findings that genetic factors serve as modifiers of responsiveness to T and E2 demonstrate the need to control for genetic background in studies of male voiding dysfunction. We also show that genetic factors could control severity of voiding dysfunction. We demonstrate the importance of E2 as a key mediator of voiding impairment, and show that the concentration of E2 in subcutaneous implants determines the severity of voiding dysfunction in mice, demonstrating that the mouse model is tunable, a factor which is important for future pharmacological intervention studies.

A Microphysiological Approach to Evaluate Effectors of Intercellular Hedgehog Signaling in Development

Paracrine signaling in the tissue microenvironment is a central mediator of morphogenesis, and modeling this dynamic intercellular activity in vitro is critical to understanding normal and abnormal development. For example, Sonic Hedgehog (Shh) signaling is a conserved mechanism involved in multiple developmental processes and strongly linked to human birth defects including orofacial clefts of the lip and palate. SHH ligand produced, processed, and secreted from the epithelial ectoderm is shuttled through the extracellular matrix where it binds mesenchymal receptors, establishing a gradient of transcriptional response that drives orofacial morphogenesis. In humans, complex interactions of genetic predispositions and environmental insults acting on diverse molecular targets are thought to underlie orofacial cleft etiology. Consequently, there is a need for tractable in vitro approaches that model this complex cellular and environmental interplay and are sensitive to disruption across the multistep signaling cascade. We developed a microplate-based device that supports an epithelium directly overlaid onto an extracellular matrix-embedded mesenchyme, mimicking the basic tissue architecture of developing orofacial tissues. SHH ligand produced from the epithelium generated a gradient of SHH-driven transcription in the adjacent mesenchyme, recapitulating the gradient of pathway activity observed in vivo. Shh pathway activation was antagonized by small molecule inhibitors of epithelial secretory, extracellular matrix transport, and mesenchymal sensing targets, supporting the use of this approach in high-content chemical screening of the complete Shh pathway. Together, these findings demonstrate a novel and practical microphysiological model with broad utility for investigating epithelial-mesenchymal interactions and environmental signaling disruptions in development.

Examining the developmental toxicity of piperonyl butoxide as a Sonic hedgehog pathway inhibitor

Piperonyl butoxide (PBO) is a semisynthetic chemical present in hundreds of pesticide formulations used in agricultural, commercial, and residential settings. PBO acts as a pesticide synergist by inhibiting insect cytochrome P450 enzymes and is often present at much higher concentrations than active insecticidal ingredients. PBO was recently discovered to also inhibit Sonic hedgehog (Shh) signaling, a key molecular pathway in embryonic development and in brain and face morphogenesis. Recent animal model studies have shown that in utero PBO exposure can cause overt craniofacial malformations or more subtle neurodevelopmental abnormalities. Related adverse developmental outcomes in humans are etiologically heterogeneous, and, while studies are limited, PBO exposure during pregnancy has been linked to neurodevelopmental deficits. Contextualized in PBO's newly recognized mechanism as a Shh signaling inhibitor, these findings support more rigorous examination of the developmental toxicity of PBO and its potential contribution to etiologically complex human birth defects. In this review, we highlight environmental sources of human PBO exposure and summarize existing animal studies examining the developmental impact of prenatal PBO exposure. Also presented are critical knowledge gaps in our understanding of PBO's pharmacokinetics and potential role in gene-environment and environment-environment interactions that should be addressed to better understand the human health impact of environmental PBO exposure.