Replication of genome wide association identified candidate genes confirm the role of common and rare variants in PAX7 and VAX1 in the etiology of nonsyndromic CL(P)

Identification of functional non-coding variants associated with orofacial cleft

Oral facial cleft (OFC) is a multifactorial disorder that can present as a cleft lip with or without cleft palate (CL/P) or a cleft palate only. Genome wide association studies (GWAS) of isolated OFC have identified common single nucleotide polymorphisms (SNPs) at the 1q32/ IRF6 locus and many other loci where, like IRF6 , the presumed OFC-relevant gene is expressed in embryonic oral epithelium. To identify the functional subset of SNPs at eight such loci we conducted a massively parallel reporter assay in a cell line derived from fetal oral epithelium, revealing SNPs with allele-specific effects on enhancer activity. We filtered these against chromatin-mark evidence of enhancers in relevant cell types or tissues, and then tested a subset in traditional reporter assays, yielding six candidates for functional SNPs in five loci (1q32/ IRF6 , 3q28/ TP63 , 6p24.3/ TFAP2A , 20q12/ MAFB , and 9q22.33/ FOXE1 ). We further tested two SNPs near IRF6 and one near FOXE1 by engineering the genome of induced pluripotent stem cells, differentiating the cells into embryonic oral epithelium, and measuring expression of IRF6 or FOXE1 and binding of transcription factors; the results strongly supported their candidacy. Conditional analyses of a meta-analysis of GWAS suggest that the two functional SNPs near IRF6 account for the majority of risk for CL/P associated with variation at this locus. This study connects genetic variation associated with orofacial cleft to mechanisms of pathogenesis.

Genetic-epigenetic interactions (meQTLs) in orofacial clefts etiology

Objectives

Nonsyndromic orofacial clefts (OFCs) etiology involves multiple genetic and environmental factors with over 60 identified risk loci; however, they account for only a minority of the estimated risk. Epigenetic factors such as differential DNA methylation (DNAm) are also associated with OFCs risk and can alter risk for different cleft types and modify OFCs penetrance. DNAm is a covalent addition of a methyl (CH3) group to the nucleotide cytosine that can lead to changes in expression of the targeted gene. DNAm can be affected by environmental influences and genetic variation via methylation quantitative loci (meQTLs). We hypothesize that aberrant DNAm and the resulting alterations in gene expression play a key role in the etiology of OFCs, and that certain common genetic variants that affect OFCs risk do so by influencing DNAm.

Methods

We used genotype from 10 cleft-associated SNPs and genome-wide DNA methylation data (Illumina 450K array) for 409 cases with OFCs and 456 controls and identified 23 cleft-associated meQTLs. We then used an independent cohort of 362 cleft-discordant sib pairs for replication. We used methylation-specific qPCR to measure methylation levels of each CpG site and combined genotypic and methylation data for an interaction analysis of each SNP-CpG pair using the R package MatrixeQTL in a linear model. We also performed a Paired T-test to analyze differences in DNA methylation between each member of the sibling pairs.

Results

We replicated 9 meQTLs, showing interactions between rs13041247 (MAFB) - cg18347630 (PLCG1) (P=0.04); rs227731 (NOG) - cg08592707 (PPM1E) (P=0.01); rs227731 (NOG) - cg10303698 (CUEDC1) (P=0.001); rs3758249 (FOXE1) - cg20308679 (FRZB) (P=0.04); rs8001641 (SPRY2) - cg19191560 (LGR4) (P=0.04); rs987525(8q24) - cg16561172(MYC) (P=0.00000963); rs7590268(THADA) - cg06873343 (TTYH3) (P=0.04); rs7078160 (VAX1) - cg09487139 (P=0.05); rs560426 (ABCA4/ARHGAP29) - cg25196715 (ABCA4/ARHGAP29) (P=0,03). Paired T-test showed significant differences for cg06873343 (TTYH3) (P=0.04); cg17103269 (LPIN3) (P=0.002), and cg19191560 (LGR4) (P=0.05).

Conclusions

Our results confirm previous evidence that some of the common non-coding variants detected through GWAS studies can influence the risk of OFCs via epigenetic mechanisms, such as DNAm, which can ultimately affect and regulate gene expression. Given the large prevalence of non-coding SNPs in most OFCs genome wide association studies, our findings can potentially address major knowledge gaps, like missing heritability, reduced penetrance, and variable expressivity associated with OFCs phenotypes.

The potential up-regulation risk of 3' UTR SNP (rs10787760 G > A) for the VAX1 gene is associated with NSCLP in the northwest Chinese population

AIMS

To investigate the association between single nucleotide polymorphisms (SNPs) in 3'UTR region of VAX1, SYT14 and PAX7 genes and the risk of non-syndromic cleft palate (NSCLP) in a northwest Chinese population.

MAIN METHODS

A case-control study was conducted in 406 normal controls and 399 NSCLP patients. Using iMLDRTM genotyping technology, eight SNPs of three genes ((rs10787760, rs7086344 at VAX1), (rs1010113, rs851114, and rs485874 at PAX7), and (rs61820397, rs4609425, rs12133399 at SYT14)) were genotyped to investigate the differences in alleles and genotype distribution frequencies between NSCLP patients and healthy controls. RNA Folding Form software was used to predict RNA secondary structure and expression vectors were constructed to explore the function of the relevant SNP. The effect of SNP polymorphism of gene transcription and translation was assessed using qPCR and Western blot analysis.

KEY FINDINGS

Among the eight SNPs of three genes, rs10787760 of VAX1 gene was found to be associated with an increased risk of NSCLP (OR = 1.341, CI = 1.004-1.790) and the GA genotype of rs10787760 increased the risk of cleft lip and/or palate (CL/P) about 1.42 times (p < 0.05), and carrying the A allele might increase the risk of NSCL/P in male (OR = 1.356, 95 % CI = 1.010-1.823). But there was no association observed with cleft palate only (CPO). Cell function experiments revealed that the G to A mutation in rs10787760 up-regulated GFP-VAX1 transcriptional level by 2.39 and 3.13 times in two cell lines respectively, and enhance the protein expression of the VAX1 gene further. RNA secondary structure study showed that the rs10787760 (G > A) had two different secondary structures in 3'UTR region.

SIGNIFICANCE

The rs10787760 variant in the 3'UTR region of VAX1 gene is associated with CL/P in northwest Chinese population. We hypothesize that the machanism of it might be caused by the RNA differenct fold in the 3'UTR region caused by the polymorphism of the gene.

LEVEL OF EVIDENCE

Original Reports.

Ezh2-dependent methylation in oral epithelia promotes secondary palatogenesis

BACKGROUND

In addition to genomic risk variants and environmental influences, increasing evidence suggests epigenetic modifications are important for orofacial development and their alterations can contribute to orofacial clefts. Ezh2 encodes a core catalytic component of the Polycomb repressive complex responsible for addition of methyl marks to Histone H3 as a mechanism of repressing target genes. The role of Ezh2 in orofacial clefts remains unknown.

AIMS

To investigate the epithelial role of Ezh2-dependent methylation in secondary palatogenesis.

METHODS

We used conditional gene-targeting methods to ablate Ezh2 in the surface ectoderm-derived oral epithelium of mouse embryos. We then performed single-cell RNA sequencing combined with immunofluorescence and RT-qPCR to investigate gene expression in conditional mutant palate. We also employed double knockout analyses of Ezh1 and Ezh2 to address if they have synergistic roles in palatogenesis.

RESULTS

We found that conditional inactivation of Ezh2 in oral epithelia results in partially penetrant cleft palate. Double knockout analyses revealed that another family member Ezh1 is dispensable in orofacial development, and it does not have synergistic roles with Ezh2 in palatogenesis. Histochemistry and single-cell RNA-seq analyses revealed dysregulation of cell cycle regulators in the palatal epithelia of Ezh2 mutant mouse embryos disrupts palatogenesis.

CONCLUSION

Ezh2-dependent histone H3K27 methylation represses expression of cell cycle regulator Cdkn1a and promotes proliferation in the epithelium of the developing palatal shelves. Loss of this regulation may perturb movement of the palatal shelves, causing a delay in palate elevation which may result in failure of the secondary palate to close altogether.

Genetic Analysis of the Single-Nucleotide Polymorphisms rs880810, rs545793, rs80094639, and rs13251901 in Nonsyndromic Oral Clefts: A Case–Parent Trio Study

Oral clefts, including cleft lip (CL), cleft palate (CP), and cleft lip and palate (CLP), are the most common types of congenital anomalies of the human face. Various genetic and environmental factors play a role in developing oral clefts. Several studies have shown the association of the PAX7 gene and the 8q24 region with these oral clefts in different populations worldwide. However, there are no reported studies on the possible connection between the PAX7 gene and the 8q24 region nucleotide variants and the risk of developing nonsyndromic oral clefts (NSOC) in the Indian population. Hence, this study aimed to test the possible association between PAX7 gene single-nucleotide polymorphisms (SNPs) rs880810, rs545793,rs80094639, and rs13251901 of the 8q24 region using a case-parent trio design. Forty case-parent trios were selected from the CLP center. Genomic DNA was isolated from the cases and their parents. The rs880810, rs545793, rs80094639, and rs13251901 were genotyped by the MassARRAY technique. PLINK software was used for statistical analysis. All the SNPs were tested for Hardy-Weinberg equilibrium. No statistical significance was found with any SNPs, as none of the genotyped SNPs showed a p -value of less than 0.05. Hence, the rs880810, rs545793, and rs80094639 of the PAX7 gene, and rs13251901 of the 8q24 region are not associated with NSOC in the Indian population.

Identification of Novel Risk Variants of Non-Syndromic Cleft Palate by Targeted Gene Panel Sequencing

Non-syndromic cleft palate (ns-CP) has a genetically heterogeneous aetiology. Numerous studies have suggested a crucial role of rare coding variants in characterizing the unrevealed component of genetic variation in ns-CP called the "missing heritability". Therefore, this study aimed to detect low-frequency variants that are implicated in ns-CP aetiology in the Polish population. For this purpose, coding regions of 423 genes associated with orofacial cleft anomalies and/or involved with facial development were screened in 38 ns-CP patients using the next-generation sequencing technology. After multistage selection and prioritisation, eight novel and four known rare variants that may influence an individual's risk of ns-CP were identified. Among detected alternations, seven were located in novel candidate genes for ns-CP, including COL17A1 (c.2435-1G>A), DLG1 (c.1586G>C, p.Glu562Asp), NHS (c.568G>C, p.Val190Leu-de novo variant), NOTCH2 (c.1997A>G, p.Tyr666Cys), TBX18 (c.647A>T, p.His225Leu), VAX1 (c.400G>A, p.Ala134Thr) and WNT5B (c.716G>T, p.Arg239Leu). The remaining risk variants were identified within genes previously linked to ns-CP, confirming their contribution to this anomaly. This list included ARHGAP29 (c.1706G>A, p.Arg569Gln), FLNB (c.3605A>G, Tyr1202Cys), IRF6 (224A>G, p.Asp75Gly-de novo variant), LRP6 (c.481C>A, p.Pro161Thr) and TP63 (c.353A>T, p.Asn118Ile). In summary, this study provides further insights into the genetic components contributing to ns-CP aetiology and identifies novel susceptibility genes for this craniofacial anomaly.

A Comprehensive Investigation on Potential Risk Factors for NSCL/P in a Rural District of Hebei Province, China

Non syndromic cleft lip with or without palate (NSCL/P), one of the most common birth defects, is closely related to various risk factors. However, information regarding risk factors for NSCL/P in rural districts in China is very limited thus far. The objective of this study was designed to identify the potential risk factors for NSCL/P in rural districts.

A comprehensive retrospective investigation including 435 NSCL/P patients and 402 healthy children was carried out in Hebei Province, China. Multiple logistic regression analysis and transmission disequilibrium test (TDT) were respectively used to identify non-genetic and genetic risk factors for NSCL/P, and then PLINK was used to explore the relationship between non-genetic and genetic risk factors.

The results showed that maternal periconceptional exposure to pesticides and herbicides, as well as low parental education level were involved in the increased risk of NSCL/P, whereas maternal folic acid and multivitamin supplementation use during preconception period were associated with the reduced risk of NSCL/P. TDT analysis identified 2 single nucleotide polymorphisms (SNPs) (rs7078160 and rs4752028) in VAX1 and one SNP (rs17563) in BMP4 as the genetic risk factors for NSCL/P. Further analysis showed that the genetic risk factors were closely related with the negative non-genetic risk factors.

Our study identified the potential risk factors for NSCL/P in rural districts, thus providing a theoretical basis for the prevention of NSCL/P occurrence.

Maternal Transmission of the PAX7 Single Nucleotide Polymorphisms among Indian Cleft Trios

Cleft lip and/or cleft palate (CL/P) is one of the most common congenital anomalies of the human face with a complex etiology involving multiple genetic and environmental factors. Several studies have shown the association of the paired box 7 ( PAX7 ) gene with CL/P in different populations worldwide. However, the current literature reveals no reported case-parent trio studies to evaluate the association between the PAX7 gene and the risk of nonsyndromic cleft lip and/or palate (NSCL/P) in the Indian population. Hence, the purpose of this study was to assess the PAX7 gene single nucleotide polymorphisms (SNPs) in the etiology of NSCL/P among the Indian cleft trios. Forty Indian case-parent trios of NSCL/P were included. The cases and their parents' genomic DNA were extracted. The SNPs rs9439714, rs1339062, rs6695765, rs742071, and rs618941of the PAX7 gene were genotyped using the Agena Bio MassARRAY analysis. The allelic transmission disequilibrium test was performed using PLINK software while pair-wise linkage disequilibrium by the Haploview program. The SNP rs9439714 showed evidence of association ( p -value = 0.02, odds ratio = 3) with NSCL/P. Considering the parent-of-origin effects, the SNPs rs9439714 and rs618941 showed an excess maternal transmission of allele C at rs9439714 ( p -value = 0.05) and G allele at rs618941 ( p -value = 0.04). The results of the present study suggested that the SNPs rs9439714 and rs618941 showed an excess maternal transmission of alleles suggestive of the possible role of the PAX7 gene involvement in the etiology of NSCL/P in the Indian population.

Role of PAX7 Gene rs766325 and rs4920520 Polymorphisms in the Etiology of Non-syndromic Cleft Lip and Palate: A Genetic Study

Non-syndromic cleft lip and palate (NSCLP) is one of the most common birth defects in humans with an overall prevalence of ∼1 in 700 live births around the world. The etiology of NSCLP is complex involving multiple genes, environmental factors, and gene-to-gene interactions. Several genome-wide associations (GWA) studies have shown the association of the paired box 7 ( PAX7 ) gene in the etiology of cleft lip and palate in different populations worldwide. However, there are no reported studies on the association between the rs766325 and rs4920520 polymorphisms and the risk of developing NSCLP in the Indian population. Hence, the present study aimed to test for the probable association between rs766325 and rs4920520 polymorphisms among NSCLP Indian population using a case-parent trio design. Forty case-parent trios were selected from the cleft lip and palate center based on the inclusion and exclusion criteria. Genomic DNA was isolated from the cases and their parents. The rs766325 and rs4920520 polymorphisms of the PAX7 gene were analyzed for their association using the MassARRAY analysis. The statistical analysis was done using the PLINK software. The rs766325 and rs4920520 polymorphisms were tested for the Hardy–Weinberg equilibrium. None of the polymorphisms showed any statistical significance. Hence, the rs766325 and rs4920520 polymorphisms of the PAX7 gene were found to be not associated with NSCLP in the Indian case-parent trios.

Genome-Wide Scan for Parent-of-Origin Effects in a sub-Saharan African Cohort With Nonsyndromic Cleft Lip and/or Cleft Palate (CL/P)

OBJECTIVE

Nonsyndromic cleft lip and/or cleft palate (NSCL/P) have multifactorial etiology where genetic factors, gene-environment interactions, stochastic factors, gene-gene interactions, and parent-of-origin effects (POEs) play cardinal roles. POEs arise when the parental origin of alleles differentially impacts the phenotype of the offspring. The aim of this study was to identify POEs that can increase risk for NSCL/P in humans using a genome-wide dataset.

METHODS

The samples (174 case-parent trios from Ghana, Ethiopia, and Nigeria) included in this study were from the African only genome wide association studies (GWAS) that was published in 2019. Genotyping of individual DNA using over 2 million multiethnic and African ancestry-specific single-nucleotide polymorphisms from the Illumina Multi-Ethnic Genotyping Array v2 15070954 A2 (genome build GRCh37/hg19) was done at the Center for Inherited Diseases Research. After quality control checks, PLINK was employed to carry out POE analysis employing the pooled subphenotypes of NSCL/P.

RESULTS

We observed possible hints of POEs at a cluster of genes at a 1 mega base pair window at the major histocompatibility complex class 1 locus on chromosome 6, as well as at other loci encompassing candidate genes such as ASB18, ANKEF1, AGAP1, GABRD, HHAT, CCT7, DNMT3A, EPHA7, FOXO3, lncRNAs, microRNA, antisense RNAs, ZNRD1, ZFAT, and ZBTB16.

CONCLUSION

Findings from our study suggest that some loci may increase the risk for NSCL/P through POEs. Additional studies are required to confirm these suggestive loci in NSCL/P etiology.

Association of genetic polymorphisms of VAX1, MAFB, and NTN1 with nonsyndromic cleft lip with or without cleft palate in Chinese population

Nonsyndromic cleft lip with or without palate (NSCL/P) is a common birth defect involving genetic factors. We conducted this case-control study to verify the association of ten single-nucleotide polymorphisms (SNPs) of six genes (VAX1, MAFB, PAX7, ABCA4, NTN1, and NOG) with NSCL/P in the Chinese population. The study included 249 NSCL/P patients, 62 nonsyndromic cleft palate only (NSCPO) patients and 480 controls. Three loci, namely, VAX1 rs7078160, MAFB rs11696257, and NTN1 rs4791774, were associated with NSCL/P (Bonferroni method adjusted p values were 0.020, 0.00031, and 0.030, respectively). We also found that the disease risk of individuals carrying both VAX1 rs7078160 and NTN1 rs4791774 was higher than those carrying only one of them (p = 4.50 × 10^-4 and 6.03 × 10^-3, respectively). SNPs of genes VAX1 rs7078160, MAFB rs11696257, and NTN1 rs4791774 increased NSCL/P risk in the Chinese population.

PAX7, PAX9 and RYK Expression in Cleft Affected Tissue

Background and Objectives: Cleft lip with or without cleft palate is one of the most common types of congenital malformations. Transcription factors paired box 7 and 9 (PAX7, PAX9) and receptor-like tyrosine kinase (RYK) have been previously associated with the formation of orofacial clefts but their exact possible involvement and interactions in the tissue of specific cleft types remains uncertain. There is a limited number of morphological studies analyzing these specific factors in cleft affected tissue due to ethical aspects and the limited amount of available tissue material. This study analyses the presence of PAX7, PAX9, and RYK immunopositive structures within different cleft affected tissue to assess their possible involvement in cleft morphopathogenesis. Materials and Methods: Cleft affected tissue was collected from non-syndromic orofacial cleft patients during cleft correcting surgery (36 patients with unilateral cleft lip, 13 patients with bilateral cleft lip, 26 patients with isolated cleft palate). Control group oral cavity tissue was obtained from 7 patients without cleft lip and palate. To evaluate the number of immunopositive structures in the cleft affected tissue and the control group, a semiquantitative counting method was used. Non-parametric statistical methods (Kruskal–Wallis H test, Mann–Whitney U test, and Spearman’s rank correlation) were used. Results: Statistically significant differences for the number of PAX7, PAX9, and RYK-positive cells were notified between the controls and the patient groups. Multiple statistically significant correlations between the factors were found in each cleft affected tissue group. Conclusions: PAX7, PAX9, and RYK have a variable involvement and interaction in postnatal morphopathogenesis of orofacial clefts. PAX7 is more associated with the formation of unilateral cleft lip, while PAX9 relates more towards the isolated cleft palate. The stable presence of RYK in all cleft types indicates its possible participation in different facial cleft formations.

The Mafb cleft-associated variant H131Q is not required for palatogenesis in the mouse

BACKGROUND

Orofacial clefts (OFCs) are common birth defects with complex etiology. Genome wide association studies for OFC have identified SNPs in and near MAFB. MAFB is a transcription factor critical for structural development of digits, kidneys, skin, and brain. MAFB is also expressed in the craniofacial region. Previous sequencing of MAFB in a Filipino population revealed a novel missense variant significantly associated with an increased risk for OFC. This MAFB variant, leading to the amino acid change H131Q, was knocked into the mouse Mafb, resulting in the MafbH131Q allele. The MafbH131Q construct was engineered to allow for deletion of Mafb ("Mafbdel ").

RESULTS

Mafbdel/del animals died shortly after birth. Conversely, MafbH131Q/H131Q mice survived into adulthood at Mendelian ratios. Mafbdel/del and MafbH131Q/H131Q heads exhibited normal macroscopic and histological appearance at all embryonic time points evaluated. The periderm was intact based on expression of keratin 6, p63, and E-cadherin. Despite no effect on craniofacial morphogenesis, H131Q inhibited the Mafb-dependent promoter activation of Arhgap29 in palatal mesenchymal, but not ectodermal-derived epithelial cells in a luciferase assay.

CONCLUSIONS

Mafb is dispensable for murine palatogenesis in vivo, and the cleft-associated variant H131Q, despite its lack of morphogenic effect, altered the expression of Arhgap29 in a cell-dependent context.

Cleft Candidate Genes and Their Products in Human Unilateral Cleft Lip Tissue

Cleft lip and palate are common congenital pathologies that affect the human population worldwide. The formation of cleft lip is associated with multiple genes and their coded proteins, which regulate the development of craniofacial region, but the exact role of these factors is not always clear. The use of morphological studies for evaluation of human cleft-affected tissue has been limited because of insufficiency of available pathological material. The aim of this study was to detect and compare the immunohistochemical expression of cleft candidate gene coded proteins (DLX4, MSX2, HOXB3, SHH, PAX7, SOX3, WNT3A, and FOXE1) in the non-syndromic unilateral cleft lip patient tissue and control group tissue. A semiquantitative counting method was used to evaluate the tissue in biotin-streptavidin-stained slides. Statistically significant differences between the patient and control groups were found for the number of immunoreactive structures for SHH (p = 0.019) and FOXE1 (p = 0.011) in the connective tissue and SOX3 (p = 0.012) in the epithelium. Multiple statistically significant very strong and strong correlations were found between the immunoreactives in cleft-affected tissue. These significant differences and various correlations indicate that multiple morphopathogenetic pathways are possibly involved in unilateral cleft lip pathogenesis. Therefore, we further discuss these possible interactions.

Genetics and signaling mechanisms of orofacial clefts

Craniofacial development involves several complex tissue movements including several fusion processes to form the frontonasal and maxillary structures, including the upper lip and palate. Each of these movements are controlled by many different factors that are tightly regulated by several integral morphogenetic signaling pathways. Subject to both genetic and environmental influences, interruption at nearly any stage can disrupt lip, nasal, or palate fusion and result in a cleft. Here, we discuss many of the genetic risk factors that may contribute to the presentation of orofacial clefts in patients, and several of the key signaling pathways and underlying cellular mechanisms that control lip and palate formation, as identified primarily through investigating equivalent processes in animal models, are examined.

Non-syndromic Cleft Palate: An Overview on Human Genetic and Environmental Risk Factors

The epithelial and mesenchymal cells involved in early embryonic facial development are guided by complex regulatory mechanisms. Any factor perturbing the growth, approach and fusion of the frontonasal and maxillary processes could result in orofacial clefts that represent the most common craniofacial malformations in humans. The rarest and, probably for this reason, the least studied form of cleft involves only the secondary palate, which is posterior to the incisive foramen. The etiology of cleft palate only is multifactorial and involves both genetic and environmental risk factors. The intention of this review is to give the reader an overview of the efforts made by researchers to shed light on the underlying causes of this birth defect. Most of the scientific papers suggesting potential environmental and genetic causes of non-syndromic cleft palate are summarized in this review, including genome-wide association and gene–environment interaction studies.

The risk of nonsyndromic cleft lip with or without cleft palate and Vax1 rs7078160 polymorphisms in southern Han Chinese

Introduction

Non-syndromic cleft lip with or without cleft palate is a common worldwide birth defect due to a combination of environmental and genetic factors. Genome-wide association studies reported the rs7078160 of Vax1 is closely related to non-syndromic cleft lip with or without cleft palate in European populations. The following studies showed the same results in Mongolian, Japanese, Filipino, Vietnamese populations etc. However, conflicting research had been reported in Chinese population,

Objective

The aim of this study was to investigate the association between the rs7078160 polymorphism and non-syndromic cleft lip with or without cleft palate in Southern Chinese patients.

Methods

In this study, we investigated the polymorphism distribution of rs7078160 in 100 complete patient trios (39 patients with non-syndromic cleft lip and palate; 36 patients with non-syndromic cleft lip only; 25 had non-syndromic cleft palate only; and their parents) from Southern ethnic Han Chinese. 60 healthy trios were selected as control. Polymerase chain reaction and Sanger sequencing were used to genotype rs7078160 in Vax1; both case–control and family-based associations were analyzed.

Results

The case–control analyses revealed the rs7078160 polymorphism was significant, associated with non-syndromic cleft lip with or without cleft palate (p = 0.04) and non-syndromic cleft lip and palate (p = 0.01), but not associated with non-syndromic cleft lip only and non-syndromic cleft palate only patients. The genotype composition of rs7078160 comprises mutated homozygous AA, heterozygous AG and wild homozygous GG. Cases with AG + AA genotypes compared with GG homozygotes showed an increased risk of non-syndromic cleft lip with or without cleft palate (p = 0.04, OR = 2.05, 95% CI: 1.01–4.16) and non-syndromic cleft lip and palate (p = 0.01, OR = 3.94, 95% CI: 1.34–11.54). In addition, we did not detect any transmission-disequilibrium in rs7078160 (p = 0.68).

Conclusion

This study suggests that rs7078160 polymorphism is a risk factor of non-syndromic cleft lip with or without cleft palate, and Vax1 is strongly associated with non-syndromic cleft lip with or without cleft palate in Southern Chinese Han populations.

Whole genome sequencing of familial isolated oesophagus atresia uncover shared structural variants

Background

Oesophageal atresia (OA) is a life-threatening developmental defect characterized by a lost continuity between the upper and lower oesophagus. The most common form is a distal connection between the trachea and the oesophagus, i.e. a tracheoesophageal fistula (TEF). The condition may be part of a syndrome or occurs as an isolated feature. The recurrence risk in affected families is increased compared to the population-based incidence suggesting contributing genetic factors.

Methods

To gain insight into gene variants and genes associated with isolated OA we conducted whole genome sequencing on samples from three families with recurrent cases affected by congenital and isolated TEF.

Results

We identified a combination of single nucleotide variants (SNVs), splice site variants (SSV) and structural variants (SV) annotated to altogether 100 coding genes in the six affected individuals.

Conclusion

This study highlights rare SVs among candidate gene variants in our individuals with OA and provides a gene framework for further investigations of genetic factors behind this malformation.

Inherited Eye Diseases with Retinal Manifestations through the Eyes of Homeobox Genes

Retinal development is under the coordinated control of overlapping networks of signaling pathways and transcription factors. The paper was conceived as a review of the data and ideas that have been formed to date on homeobox genes mutations that lead to the disruption of eye organogenesis and result in inherited eye/retinal diseases. Many of these diseases are part of the same clinical spectrum and have high genetic heterogeneity with already identified associated genes. We summarize the known key regulators of eye development, with a focus on the homeobox genes associated with monogenic eye diseases showing retinal manifestations. Recent advances in the field of genetics and high-throughput next-generation sequencing technologies, including single-cell transcriptome analysis have allowed for deepening of knowledge of the genetic basis of inherited retinal diseases (IRDs), as well as improve their diagnostics. We highlight some promising avenues of research involving molecular-genetic and cell-technology approaches that can be effective for IRDs therapy. The most promising neuroprotective strategies are aimed at mobilizing the endogenous cellular reserve of the retina.

PAX7 nucleotide variants and the risk of non-syndromic orofacial clefts in the Polish population

OBJECTIVE

The etiology of non-syndromic cleft lip with or without cleft palate (nsCL/P) is multifactorial, heterogeneous, and still not completely understood. The aim of the present study was to examine the associations between common and rare PAX7 nucleotide variants and the risk of this common congenital anomaly in a Polish population.

SUBJECTS AND METHODS

Eight top nsCL/P-associated PAX7 variants identified in our cleft genome-wide association study (GWAS) were selected for replication analysis in an independent group of patients and controls (n = 247 and n = 445, respectively). In addition, mutation screening of the PAX7 protein-coding region was conducted.

RESULTS

Analysis of the pooled data from the GWAS and replication study confirmed that common PAX7 nucleotide variants are significantly associated with the increased risk of nsCL/P. The strongest individual variant was rs1339062 (c.586 + 15617T > C) with a p-value = 2.47E-05 (OR = 1.4, 95%CI: 1.20-1.64). Sequencing analysis identified a novel synonymous PAX7 substitution (c.87G > A, p.Val29Val) in a single patient with nsCLP. This transition located in the early exonic position was predicted to disrupt potential splice enhancer elements.

CONCLUSION

Our study confirmed that PAX7 is a strong candidate gene for nsCL/P. Nucleotide variants of this gene contribute to the etiology of nsCL/P in the homogenous Polish population.

MicroRNA-655-3p and microRNA-497-5p inhibit cell proliferation in cultured human lip cells through the regulation of genes related to human cleft lip

Background

The etiology of cleft lip with or without cleft palate (CL/P), a common congenital birth defect, is complex and involves the contribution of genetic and environmental factors. Although many candidate genes have been identified, the regulation and interaction of these genes in CL/P remain unclear. In addition, the contribution of microRNAs (miRNAs), non-coding RNAs that regulate the expression of multiple genes, to the etiology of CL/P is largely unknown.

Methods

To identify the signatures of causative biological pathways for human CL/P, we conducted a systematic literature review for human CL/P candidate genes and subsequent bioinformatics analyses. Functional enrichment analyses of the candidate CL/P genes were conducted using the pathway databases GO and KEGG. The miRNA-mediated post-transcriptional regulation of the CL/P candidate genes was analyzed with miRanda, PITA, and TargetScan, and miRTarbase. Genotype-phenotype association analysis was conducted using GWAS. The functional significance of the candidate miRNAs was evaluated experimentally in cell proliferation and target gene regulation assays in human lip fibroblasts.

Results

Through an extensive search of the main biomedical databases, we mined 177 genes with mutations or association/linkage reported in individuals with CL/P, and considered them as candidate genes for human CL/P. The genotype-phenotype association study revealed that mutations in 12 genes (ABCA4, ADAM3A, FOXE1, IRF6, MSX2, MTHFR, NTN1, PAX7, TP63, TPM1, VAX1, and WNT9B) were significantly associated with CL/P. In addition, our bioinformatics analysis predicted 16 microRNAs (miRNAs) to be post-transcriptional regulators of CL/P genes. To validate the bioinformatics results, the top six candidate miRNAs (miR-124-3p, miR-369-3p, miR-374a-5p, miR-374b-5p, miR-497-5p, and miR-655-3p) were evaluated by cell proliferation/survival assays and miRNA-gene regulation assays in cultured human lip fibroblasts. We found that miR-497-5p and miR-655-3p significantly suppressed cell proliferation in these cells. Furthermore, the expression of the predicted miRNA-target genes was significantly downregulated by either miR-497-5p or miR-655-3p mimic.

Conclusion

Expression of miR-497-5p and miR-655-3p suppresses cell proliferation through the regulation of human CL/P-candidate genes. This study provides insights into the role of miRNAs in the etiology of CL/P and suggests possible strategies for the diagnosis of CL/P.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0535-2) contains supplementary material, which is available to authorized users.

Evidence for DNA methylation mediating genetic liability to non-syndromic cleft lip/palate

AIM

To determine if nonsyndromic cleft lip with or without cleft palate (nsCL/P) genetic risk variants influence liability to nsCL/P through gene regulation pathways, such as those involving DNA methylation.

MATERIALS & METHODS

nsCL/P genetic summary data and methylation data from four studies were used in conjunction with Mendelian randomization and joint likelihood mapping to investigate potential mediation of nsCL/P genetic variants.

RESULTS & CONCLUSION

Evidence was found at VAX1 (10q25.3), LOC146880 (17q23.3) and NTN1 (17p13.1), that liability to nsCL/P and variation in DNA methylation might be driven by the same genetic variant, suggesting that genetic variation at these loci may increase liability to nsCL/P by influencing DNA methylation. Follow-up analyses using different tissues and gene expression data provided further insight into possible biological mechanisms.

Haploinsufficiency of Homeodomain Proteins Six3, Vax1, and Otx2 Causes Subfertility in Mice via Distinct Mechanisms

Haploinsufficiency occurs when loss of one copy of a diploid gene (hemizygosity) causes a phenotype. It is relatively rare, in that most genes can produce sufficient mRNA and protein from a single copy to prevent any loss of normal activity and function. Reproduction is a complex process relying on migration of GnRH neurons from the olfactory placode to the hypothalamus during development. We have studied 3 different homeodomain genes Otx2, Vax1, and Six3 and found that the deletion of one allele for any of these genes in mice produces subfertility or infertility in one or both sexes, despite the presence of one intact allele. All 3 heterozygous mice have reduced numbers of GnRH neurons, but the mechanisms of subfertility differ significantly. This review compares the subfertility phenotypes and their mechanisms.

Medical genetics and genomic medicine in Nigeria

Medical genetics and genomic medicine in Nigeria.

Vax1 Plays an Indirect Role in the Etiology of Murine Cleft Palate

Cleft lip with or without palate (CLP) and isolated cleft palate (CP) are common human developmental malformations with a complex etiology that reflects a failure of normal facial development. VAX1 encodes a homeobox-containing transcription factor identified as a candidate gene for CLP in human populations, with targeted deletion in mice associated with multiple anomalies, including disruption of the visual apparatus and basal forebrain, lobar holoprosencephaly, and CP. We have investigated Vax1 function during murine palatogenesis but found no evidence for a direct role in this process. Vax1 is not expressed in the developing palate and mutant palatal shelves elevate above the tongue, demonstrating morphology and proliferation indices indistinguishable from wild type. However, mutant mice did have a large midline cavity originating from the embryonic forebrain situated beneath the floor of the hypothalamus and extending through the nasal cavity to expand this region and prevent approximation of the palatal shelves. Interestingly, despite strong expression of Vax1 in ectoderm of the medial nasal processes, the upper lip remained intact in mutant mice. We found further evidence of disrupted craniofacial morphology in Vax1 mutants, including truncation of the midface associated with reduced cell proliferation in forebrain neuroectoderm and frontonasal mesenchyme. Sonic hedgehog (Shh) signal transduction was downregulated in the mutant forebrain, consistent with a role for Vax1 in mediating transduction of this pathway. However, Shh was also reduced in this region, suggestive of a Shh-Vax1 feedback loop during early development of the forebrain and a likely mechanism for the underlying lobar holoprosencephaly. Despite significant associations between VAX1 and human forms of CLP, we find no evidence of a direct role for this transcription factor in development of this region in a mutant mouse model.

Imputation of orofacial clefting data identifies novel risk loci and sheds light on the genetic background of cleft lip ± cleft palate and cleft palate only

Nonsyndromic cleft lip with or without cleft palate (nsCL/P) is among the most common human birth defects with multifactorial etiology. Here, we present results from a genome-wide imputation study of nsCL/P in which, after adding replication cohort data, four novel risk loci for nsCL/P are identified (at chromosomal regions 2p21, 14q22, 15q24 and 19p13). On a systematic level, we show that the association signals within this high-density dataset are enriched in functionally-relevant genomic regions that are active in both human neural crest cells (hNCC) and mouse embryonic craniofacial tissue. This enrichment is also detectable in hNCC regions primed for later activity. Using GCTA analyses, we suggest that 30% of the estimated variance in risk for nsCL/P in the European population can be attributed to common variants, with 25.5% contributed to by the 24 risk loci known to date. For each of these, we identify credible SNPs using a Bayesian refinement approach, with two loci harbouring only one probable causal variant. Finally, we demonstrate that there is no polygenic component of nsCL/P detectable that is shared with nonsyndromic cleft palate only (nsCPO). Our data suggest that, while common variants are strongly contributing to risk for nsCL/P, they do not seem to be involved in nsCPO which might be more often caused by rare deleterious variants. Our study generates novel insights into both nsCL/P and nsCPO etiology and provides a systematic framework for research into craniofacial development and malformation.

Genetics and genomics etiology of nonsyndromic orofacial clefts

Orofacial clefts (OFC) are complex birth defects. Studies using contemporary genomic techniques, bioinformatics, and statistical analyses have led to appreciable advances in identifying the causes of syndromic forms of clefts. This commentary gives an overview of the important cleft gene discoveries found using various genomic methods and tools.

Association Studies and Direct DNA Sequencing Implicate Genetic Susceptibility Loci in the Etiology of Nonsyndromic Orofacial Clefts in Sub-Saharan African Populations

Orofacial clefts (OFCs) are congenital dysmorphologies of the human face and oral cavity, with a global incidence of 1 per 700 live births. These anomalies exhibit a multifactorial pattern of inheritance, with genetic and environmental factors both playing crucial roles. Many loci have been implicated in the etiology of nonsyndromic cleft lip with or without cleft palate (NSCL/P) in populations of Asian and European ancestries, through genome-wide association studies and candidate gene studies. However, few populations of African descent have been studied to date. Here, the authors show evidence of an association of some loci with NSCL/P and nonsyndromic cleft palate only (NSCPO) in cohorts from Africa (Ghana, Ethiopia, and Nigeria). The authors genotyped 48 single-nucleotide polymorphisms that were selected from previous genome-wide association studies and candidate gene studies. These markers were successfully genotyped on 701 NSCL/P and 163 NSCPO cases, 1,070 unaffected relatives, and 1,078 unrelated controls. The authors also directly sequenced 7 genes in 184 nonsyndromic OFC (NSOFC) cases and 96 controls from Ghana. Population-specific associations were observed in the case-control analyses of the subpopulations, with West African subpopulations (Ghana and Nigeria) showing a similar pattern of associations. In meta-analyses of the case-control cohort, PAX7 (rs742071, P = 5.10 × 10(-3)), 8q24 (rs987525, P = 1.22 × 10(-3)), and VAX1 (rs7078160, P = 0.04) were nominally associated with NSCL/P, and MSX1 (rs115200552, P = 0.01), TULP4 (rs651333, P = 0.04), CRISPLD2 (rs4783099, P = 0.02), and NOG1 (rs17760296, P = 0.04) were nominally associated with NSCPO. Moreover, 7 loci exhibited evidence of threshold overtransmission in NSOFC cases through the transmission disequilibrium test and through analyses of the family-based association for disease traits. Through DNA sequencing, the authors also identified 2 novel, rare, potentially pathogenic variants (p.Asn323Asp and p.Lys426IlefsTer6) in ARHGAP29 In conclusion, the authors have shown evidence for the association of many loci with NSCL/P and NSCPO. To the best of this knowledge, this study is the first to demonstrate any of these association signals in any African population.

Conditioning adaptive combination of P-values method to analyze case-parent trios with or without population controls

Detection of rare causal variants can help uncover the etiology of complex diseases. Recruiting case-parent trios is a popular study design in family-based studies. If researchers can obtain data from population controls, utilizing them in trio analyses can improve the power of methods. The transmission disequilibrium test (TDT) is a well-known method to analyze case-parent trio data. It has been extended to rare-variant association testing (abbreviated as "rvTDT"), with the flexibility to incorporate population controls. The rvTDT method is robust to population stratification. However, power loss may occur in the conditioning process. Here we propose a "conditioning adaptive combination of P-values method" (abbreviated as "conADA"), to analyze trios with/without unrelated controls. By first truncating the variants with larger P-values, we decrease the vulnerability of conADA to the inclusion of neutral variants. Moreover, because the test statistic is developed by conditioning on parental genotypes, conADA generates valid statistical inference in the presence of population stratification. With regard to statistical methods for next-generation sequencing data analyses, validity may be hampered by population stratification, whereas power may be affected by the inclusion of neutral variants. We recommend conADA for its robustness to these two factors (population stratification and the inclusion of neutral variants).

A multi-ethnic genome-wide association study identifies novel loci for non-syndromic cleft lip with or without cleft palate on 2p24.2, 17q23 and 19q13

Orofacial clefts (OFCs), which include non-syndromic cleft lip with or without cleft palate (CL/P), are among the most common birth defects in humans, affecting approximately 1 in 700 newborns. CL/P is phenotypically heterogeneous and has a complex etiology caused by genetic and environmental factors. Previous genome-wide association studies (GWASs) have identified at least 15 risk loci for CL/P. As these loci do not account for all of the genetic variance of CL/P, we hypothesized the existence of additional risk loci. We conducted a multiethnic GWAS in 6480 participants (823 unrelated cases, 1700 unrelated controls and 1319 case-parent trios) with European, Asian, African and Central and South American ancestry. Our GWAS revealed novel associations on 2p24 near FAM49A, a gene of unknown function (P = 4.22 × 10^-8), and 19q13 near RHPN2, a gene involved in organizing the actin cytoskeleton (P = 4.17 × 10^-8). Other regions reaching genome-wide significance were 1p36 (PAX7), 1p22 (ARHGAP29), 1q32 (IRF6), 8q24 and 17p13 (NTN1), all reported in previous GWASs. Stratification by ancestry group revealed a novel association with a region on 17q23 (P = 2.92 × 10^-8) among individuals with European ancestry. This region included several promising candidates including TANC2, an oncogene required for development, and DCAF7, a scaffolding protein required for craniofacial development. In the Central and South American ancestry group, significant associations with loci previously identified in Asian or European ancestry groups reflected their admixed ancestry. In summary, we have identified novel CL/P risk loci and suggest new genes involved in craniofacial development, confirming the highly heterogeneous etiology of OFCs.

Epidemiology, Etiology, and Treatment of Isolated Cleft Palate

Isolated cleft palate (CPO) is the rarest form of oral clefting. The incidence of CPO varies substantially by geography from 1.3 to 25.3 per 10,000 live births, with the highest rates in British Columbia, Canada and the lowest rates in Nigeria, Africa. Stratified by ethnicity/race, the highest rates of CPO are observed in non-Hispanic Whites and the lowest in Africans; nevertheless, rates of CPO are consistently higher in females compared to males. Approximately fifty percent of cases born with cleft palate occur as part of a known genetic syndrome or with another malformation (e.g., congenital heart defects) and the other half occur as solitary defects, referred to often as non-syndromic clefts. The etiology of CPO is multifactorial involving genetic and environmental risk factors. Several animal models have yielded insight into the molecular pathways responsible for proper closure of the palate, including the BMP, TGF-β, and SHH signaling pathways. In terms of environmental exposures, only maternal tobacco smoke has been found to be strongly associated with CPO. Some studies have suggested that maternal glucocorticoid exposure may also be important. Clearly, there is a need for larger epidemiologic studies to further investigate both genetic and environmental risk factors and gene-environment interactions. In terms of treatment, there is a need for long-term comprehensive care including surgical, dental and speech pathology. Overall, five main themes emerge as critical in advancing research: (1) monitoring of the occurrence of CPO (capacity building); (2) detailed phenotyping of the severity (biology); (3) understanding of the genetic and environmental risk factors (primary prevention); (4) access to early detection and multidisciplinary treatment (clinical services); and (5) understanding predictors of recurrence and possible interventions among families with a child with CPO (secondary prevention).

New insights into craniofacial malformations

Development of the human skull and face is a highly orchestrated and complex three-dimensional morphogenetic process, involving hundreds of genes controlling the coordinated patterning, proliferation and differentiation of tissues having multiple embryological origins. Craniofacial malformations that occur because of abnormal development (including cleft lip and/or palate, craniosynostosis and facial dysostoses), comprise over one-third of all congenital birth defects. High-throughput sequencing has recently led to the identification of many new causative disease genes and functional studies have clarified their mechanisms of action. We present recent findings in craniofacial genetics and discuss how this information together with developmental studies in animal models is helping to increase understanding of normal craniofacial development.

Developments in our understanding of the genetic basis of birth defects

Birth defects are a major cause of morbidity and mortality worldwide. There has been much progress in understanding the genetic basis of familial and syndromic forms of birth defects. However, the etiology of nonsydromic birth defects is not well-understood. Although there is still much work to be done, we have many of the tools needed to accomplish the task. Advances in next-generation sequencing have introduced a sea of possibilities, from disease-gene discovery to clinical screening and diagnosis. These advances have been fruitful in identifying a host of candidate disease genes, spanning the spectrum of birth defects. With the advent of CRISPR-Cas9 gene editing, researchers now have a precise tool for characterizing this genetic variation in model systems. Work in model organisms has also illustrated the importance of epigenetics in human development and birth defects etiology. Here we review past and current knowledge in birth defects genetics. We describe genotyping and sequencing methods for the detection and analysis of rare and common variants. We remark on the utility of model organisms and explore epigenetics in the context of structural malformation. We conclude by highlighting approaches that may provide insight into the complex genetics of birth defects.

Palatogenesis and cutaneous repair: A two-headed coin

BACKGROUND

The reparative mechanism that operates following post-natal cutaneous injury is a fundamental survival function that requires a well-orchestrated series of molecular and cellular events. At the end, the body will have closed the hole using processes like cellular proliferation, migration, differentiation and fusion.

RESULTS

These processes are similar to those occurring during embryogenesis and tissue morphogenesis. Palatogenesis, the formation of the palate from two independent palatal shelves growing towards each other and fusing, intuitively, shares many similarities with the closure of a cutaneous wound from the two migrating epithelial fronts.

CONCLUSIONS

In this review, we summarize the current information on cutaneous development, wound healing, palatogenesis and orofacial clefting and propose that orofacial clefting and wound healing are conserved processes that share common pathways and gene regulatory networks.

Genome-wide association studies of suicidal behaviors: a review

Suicidal behaviors represent a fatal dimension of mental ill-health, involving both environmental and heritable (genetic) influences. The putative genetic components of suicidal behaviors have until recent years been mainly investigated by hypothesis-driven research (of "candidate genes"). But technological progress in genotyping has opened the possibilities towards (hypothesis-generating) genomic screens and novel opportunities to explore polygenetic perspectives, now spanning a wide array of possible analyses falling under the term Genome-Wide Association Study (GWAS). Here we introduce and discuss broadly some apparent limitations but also certain developing opportunities of GWAS. We summarize the results from all the eight GWAS conducted up to date focused on suicidality outcomes; treatment emergent suicidal ideation (3 studies), suicide attempts (4 studies) and completed suicides (1 study). Clearly, there are few (if any) genome-wide significant and reproducible findings yet to be demonstrated. We then discuss and pinpoint certain future considerations in relation to sample sizes, the units of genetic associations used, study designs and outcome definitions, psychiatric diagnoses or biological measures, as well as the use of genomic sequencing. We conclude that GWAS should have a lot more potential to show in the case of suicidal outcomes, than what has yet been realized.

Rare functional variants in genome-wide association identified candidate genes for nonsyndromic clefts in the African population

Nonsyndromic clefts of the lip and palate (NSCLP) are complex genetic traits. Together, they are classified as one of the most common birth defects with a prevalence of 1/700 live births. Genome-wide association studies (GWAS) for nonsyndromic cleft lip with or without cleft palate (NSCL[P]) revealed significant association for common single nucleotide polymorphisms near genes involved in craniofacial development i.e., MAFB, PAX7, VAX1, ARHGAP29 (ABCA4 locus), and IRF6. Sequencing of protein coding regions of the NSCL(P) GWAS candidate genes or adjacent genes suggest a role for rare functional variants. Replication studies in the African population did not observe any significant association with the GWAS candidate genes. On the other hand, the role of rare functional variants in GWAS candidate genes has not been evaluated in the African population. We obtained saliva samples from case triads in Nigeria and Ethiopia for Sanger sequencing of the GWAS candidate genes (MAFB, PAX7, VAX1, ARHGAP29, and IRF6) in order to identify rare functional variants. A total of 220 African samples (140 Nigerians and 80 Ethiopians) were sequenced and we found the following new rare variants- p.His165Asn in the MAFB gene, p.Asp428Asn in the PAX7, a splice-site variant that creates a new donor splice-site in PAX7. We also found three previously reported missense variants p.Gly466Ser in PAX7; p.Leu913Ser and Arg955His in ARHGAP29. No de novo mutations were found. Future genome-wide association and sequencing studies should be conducted using samples from Africa in order to identify new molecular genetic factors that contribute to the etiology of NSCLP.

Heterozygous deletion of ventral anterior homeobox (vax1) causes subfertility in mice

The known genetic causes of idiopathic hypogonadotropic hypogonadism (IHH) are often associated with the loss of GnRH neurons, leading to the disruption of the hypothalamic pituitary gonadal axis and subfertility. The majority of IHH cases have unknown origins and likely arise from compound mutations in more than one gene. Here we identify the homeodomain transcription factor ventral anterior homeobox1 (Vax1) as a potential genetic contributor to polygenic IHH. Although otherwise healthy, male and female Vax1 heterozygous (HET) mice are subfertile, indicating dosage sensitivity for the Vax1 allele. Although Vax1 mRNA is expressed in the pituitary, hypothalamus, and testis, we did not detect Vax1 mRNA in the sperm, ovary, or isolated pituitary gonadotropes. Whereas Vax1 HET females produced normal numbers of superovulated oocytes, corpora lutea numbers were reduced along with a slight increase in circulating basal LH and estrogen. The subfertility originated in the hypothalamus in which kisspeptin and GnRH transcripts were altered along with a substantial reduction of GnRH neuron number. Although the pituitary responded normally to a GnRH challenge, diestrus females had reduced LHβ and FSHβ in diestrus. Furthermore, Vax1 HET males had reduced GnRH mRNA and neuron numbers, whereas the pituitary had normal transcript levels and response to GnRH. Interestingly, the Vax1 HET males had an 88% reduction of motile sperm. Taken together, our data suggest that Vax1 HET subfertility originates in the hypothalamus by disrupting the hypothalamic-pituitary-gonadal axis. In addition, male subfertility may also be due to an unknown effect of Vax1 in the testis.

Nonsyndromic cleft lip with or without cleft palate: Increased burden of rare variants within Gremlin-1, a component of the bone morphogenetic protein 4 pathway

BACKGROUND

The genes Gremlin-1 (GREM1) and Noggin (NOG) are components of the bone morphogenetic protein 4 pathway, which has been implicated in craniofacial development. Both genes map to recently identified susceptibility loci (chromosomal region 15q13, 17q22) for nonsyndromic cleft lip with or without cleft palate (nsCL/P). The aim of the present study was to determine whether rare variants in either gene are implicated in nsCL/P etiology.

METHODS

The complete coding regions, untranslated regions, and splice sites of GREM1 and NOG were sequenced in 96 nsCL/P patients and 96 controls of Central European ethnicity. Three burden and four nonburden tests were performed. Statistically significant results were followed up in a second case-control sample (n = 96, respectively). For rare variants observed in cases, segregation analyses were performed.

RESULTS

In NOG, four rare sequence variants (minor allele frequency < 1%) were identified. Here, burden and nonburden analyses generated nonsignificant results. In GREM1, 33 variants were identified, 15 of which were rare. Of these, five were novel. Significant p-values were generated in three nonburden analyses. Segregation analyses revealed incomplete penetrance for all variants investigated.

CONCLUSION

Our study did not provide support for NOG being the causal gene at 17q22. However, the observation of a significant excess of rare variants in GREM1 supports the hypothesis that this is the causal gene at chr. 15q13. Because no single causal variant was identified, future sequencing analyses of GREM1 should involve larger samples and the investigation of regulatory elements.

Disrupting hedgehog and WNT signaling interactions promotes cleft lip pathogenesis

Cleft lip, which results from impaired facial process growth and fusion, is one of the most common craniofacial birth defects. Many genes are known to be involved in the etiology of this disorder; however, our understanding of cleft lip pathogenesis remains incomplete. In the present study, we uncovered a role for sonic hedgehog (SHH) signaling during lip fusion. Mice carrying compound mutations in hedgehog acyltransferase (Hhat) and patched1 (Ptch1) exhibited perturbations in the SHH gradient during frontonasal development, which led to hypoplastic nasal process outgrowth, epithelial seam persistence, and cleft lip. Further investigation revealed that enhanced SHH signaling restricts canonical WNT signaling in the lambdoidal region by promoting expression of genes encoding WNT inhibitors. Moreover, reduction of canonical WNT signaling perturbed p63/interferon regulatory factor 6 (p63/IRF6) signaling, resulting in increased proliferation and decreased cell death, which was followed by persistence of the epithelial seam and cleft lip. Consistent with our results, mutations in genes that disrupt SHH and WNT signaling have been identified in both mice and humans with cleft lip. Collectively, our data illustrate that altered SHH signaling contributes to the etiology and pathogenesis of cleft lip through antagonistic interactions with other gene regulatory networks, including the canonical WNT and p63/IRF6 signaling pathways.

Novel IRF6 mutations in families with Van Der Woude syndrome and popliteal pterygium syndrome from sub-Saharan Africa

Orofacial clefts (OFC) are complex genetic traits that are often classified as syndromic or nonsyndromic clefts. Currently, there are over 500 types of syndromic clefts in the Online Mendelian Inheritance in Man (OMIM) database, of which Van der Woude syndrome (VWS) is one of the most common (accounting for 2% of all OFC). Popliteal pterygium syndrome (PPS) is considered to be a more severe form of VWS. Mutations in the IRF6 gene have been reported worldwide to cause VWS and PPS. Here, we report studies of families with VWS and PPS in sub-Saharan Africa. We screened the DNA of eight families with VWS and one family with PPS from Nigeria and Ethiopia by Sanger sequencing of the most commonly affected exons in IRF6 (exons 3, 4, 7, and 9). For the VWS families, we found a novel nonsense variant in exon 4 (p.Lys66X), a novel splice-site variant in exon 4 (p.Pro126Pro), a novel missense variant in exon 4 (p.Phe230Leu), a previously reported splice-site variant in exon 7 that changes the acceptor splice site, and a known missense variant in exon 7 (p.Leu251Pro). A previously known missense variant was found in exon 4 (p.Arg84His) in the PPS family. All the mutations segregate in the families. Our data confirm the presence of IRF6-related VWS and PPS in sub-Saharan Africa and highlights the importance of screening for novel mutations in known genes when studying diverse global populations. This is important for counseling and prenatal diagnosis for high-risk families.

Pax7 is regulated by cMyb during early neural crest development through a novel enhancer

The neural crest (NC) is a migratory population of cells unique to vertebrates that generates many diverse derivatives. NC cells arise during gastrulation at the neural plate border (NPB), which is later elevated as the neural folds (NFs) form and fuse in the dorsal region of the closed neural tube, from where NC cells emigrate. In chick embryos, Pax7 is an early marker, and necessary component of NC development. Unlike other early NPB markers, which are co-expressed in lateral ectoderm, medial neural plate or posterior-lateral mesoderm, Pax7 early expression seems more restricted to the NPB. However, the molecular mechanisms controlling early Pax7 expression remain poorly understood. Here, we identify a novel enhancer of Pax7 in avian embryos that replicates the expression of Pax7 associated with early NC development. Expression from this enhancer is found in early NPB, NFs and early emigrating NC, but unlike Pax7, which is also expressed in mesodermal derivatives, this enhancer is not active in somites. Further analysis demonstrates that cMyb is able to interact with this enhancer and modulates reporter and endogenous early Pax7 expression; thus, cMyb is identified as a novel regulator of Pax7 in early NC development.