Exploring GRHL3 polymorphisms and SNP-SNP interactions in the risk of non-syndromic oral clefts in the Brazilian population

Molecular investigation in individuals with orofacial clefts and microphthalmia-anophthalmia-coloboma spectrum

This study describes genomic findings among individuals with both orofacial clefts (OC) and microphthalmia/anophthalmia/coloboma (MAC) recorded in the Brazilian Database on Craniofacial Anomalies (BDCA). Chromosomal microarray analysis (CMA) and Whole Exome Sequencing (WES) were performed in 17 individuals with OC-MAC. Clinical interpretation of molecular findings was based on data available at the BDCA and on re-examination. No copy number variants (CNVs) classified as likely pathogenic or pathogenic were detected by CMA. WES allowed a conclusive diagnosis in six individuals (35.29%), two of them with variants in the CHD7 gene, and the others with variants in the TFAP2A, POMT1, PTPN11, and TP63 genes with the following syndromes: CHARGE, CHD7-spectrum, Branchiooculofacial, POMT1-spectrum, LEOPARD, and ADULT. Variants of uncertain significance (VUS) possibly associated to the phenotypes were found in six other individuals. Among the individuals with VUSes, three individuals presented variants in genes associated to defects of cilia structure and/or function, including DYNC2H1, KIAA0586, WDR34, INTU, RPGRIP1L, KIF7, and LMNA. These results show that WES was the most effective molecular approach for OC-MAC in this cohort. This study also reinforces the genetic heterogeneity of OC-MAC, and the importance of genes related to ciliopathies in this phenotype.

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.

Genetic markers for non-syndromic orofacial clefts in populations of European ancestry: a meta-analysis

To date, the involvement of various genetic markers in the aetiopathogenesis of non-syndromic orofacial cleft (nsOFC) has been extensively studied. In the present study, we focused on studies performed on populations of European ancestry to systematically review the available literature to define relevant genetic risk factors for nsOFC. Eligible studies were obtained by searching Ovid Medline and Ovid Embase. We gathered the genetic markers from population-based case–control studies on nsOFC, and conducted meta-analysis on the repeatedly reported markers. Whenever possible, we performed stratified analysis based on different nsOFC phenotypes, using allelic, dominant, recessive and overdominant genetic models. Effect sizes were expressed as pooled odds ratios (ORs) with 95% confidence intervals (CIs), and p ≤ 0.05 were considered statistically significant. A total of 84 studies were eligible for this systematic review, with > 700 markers included. Of these, 43 studies were included in the meta-analysis. We analysed 47 genetic variants in 30 genes/loci, which resulted in 226 forest plots. There were statistically significant associations between at least one of the nsOFC phenotypes and 19 genetic variants in 13 genes/loci. These data suggest that IRF6, GRHL3, 8q24, VAX1, TGFA, FOXE1, ABCA4, NOG, GREM1, AXIN2, DVL2, WNT3A and WNT5A have high potential as biomarkers of nsOFC in populations of European descent. Although other meta-analyses that included European samples have been performed on a limited number of genetic variants, this study represents the first meta-analysis of all genetic markers that have been studied in connection with nsOFC in populations of European ancestry.

CYRI-A limits invasive migration through macropinosome formation and integrin uptake regulation

The Scar/WAVE complex drives actin nucleation during cell migration. Interestingly, the same complex is important in forming membrane ruffles during macropinocytosis, a process mediating nutrient uptake and membrane receptor trafficking. Mammalian CYRI-B is a recently described negative regulator of the Scar/WAVE complex by RAC1 sequestration, but its other paralogue, CYRI-A, has not been characterized. Here, we implicate CYRI-A as a key regulator of macropinosome formation and integrin internalization. We find that CYRI-A is transiently recruited to nascent macropinosomes, dependent on PI3K and RAC1 activity. CYRI-A recruitment precedes RAB5A recruitment but follows sharply after RAC1 and actin signaling, consistent with it being a local inhibitor of actin polymerization. Depletion of both CYRI-A and -B results in enhanced surface expression of the α5β1 integrin via reduced internalization. CYRI depletion enhanced migration, invasion, and anchorage-independent growth in 3D. Thus, CYRI-A is a dynamic regulator of macropinocytosis, functioning together with CYRI-B to regulate integrin trafficking.

Association of BMP4 rs17563 Polymorphism with Nonsyndromic Cleft Lip with or without Cleft Palate Risk: Literature Review and Comprehensive Meta-Analysis

Although published individual studies have reported associations between BMP4 rs17563 polymorphism and nonsyndromic cleft lip with or without cleft palate (NSCLP) risk, the results are conflicting. This meta-analysis was conducted to assess the association based on multiple studies. Methods: A comprehensive literature search up to October 1st, 2019 was performed using PubMed, Science Direct, China National Knowledge Infrastructure (CNKI), and Wanfang databases. Results: Fourteen case-control studies with 2,058 NSCLP cases and 2,557 controls were selected. There was no significant association between BMP4 rs17563 polymorphism and risk of NSCLP overall. Subgroup analysis revealed that BMP4 rs17563 polymorphism was associated with NSCLP risk in Chinese and Brazilian populations. Conclusions: This meta-analysis suggests that BMP4 rs17563 polymorphism was not associated with NSCLP risk in overall population. However, BMP4 rs17563 polymorphism may be a risk factor for development of NSCLP in Chinese and Brazilians.

Identification of Novel Variants in Cleft Palate-Associated Genes in Brazilian Patients With Non-syndromic Cleft Palate Only

The identification of genetic risk factors for non-syndromic oral clefts is of great importance for better understanding the biological processes related to this heterogeneous and complex group of diseases. Herein we applied whole-exome sequencing to identify potential variants related to non-syndromic cleft palate only (NSCPO) in the multiethnic Brazilian population. Thirty NSCPO samples and 30 sex- and genetic ancestry-matched healthy controls were pooled (3 pools with 10 samples for each group) and subjected to whole-exome sequencing. After filtering, the functional affects, individually and through interactions, of the selected variants and genes were assessed by bioinformatic analyses. As a group, 399 variants in 216 genes related to palatogenesis/cleft palate, corresponding to 6.43%, were exclusively identified in the NSCPO pools. Among those genes are 99 associated with syndromes displaying cleft palate in their clinical spectrum and 92 previously related to cleft lip palate. The most significantly biological processes and pathways overrepresented in the NSCPO-identified genes were associated with the folic acid metabolism, highlighting the interaction between LDL receptor-related protein 6 (LRP6) and 5-methyltetrahydrofolate-homocysteine methyltransferase (MTR) that interconnect two large networks. This study yields novel data on characterization of specific variants and complex processes and pathways related to NSCPO, including many variants in genes of the folate/homocysteine pathway, and confirms that variants in genes related to syndromic cleft palate and cleft lip-palate may cause NSCPO.

Potential interactions among single nucleotide polymorphisms in bone- and cartilage-related genes in skeletal malocclusions

OBJECTIVE

To investigate SNPs in bone- and cartilage-related genes and their interaction in the aetiology of sagittal and vertical skeletal malocclusions.

SETTINGS AND SAMPLE POPULATION

This study included 143 patients and classified as follows: skeletal class I (n = 77), class II (n = 47) and class III (n = 19); maxillary retrusion (n = 39), protrusion (n = 52) and well-positioned maxilla (n = 52); mandibular retrognathism (n = 50), prognathism (n = 50) and well-positioned mandible (n = 43); normofacial (n = 72), dolichofacial (n = 55) and brachyfacial (n = 16).

MATERIALS AND METHODS

Steiner's ANB, SNA, SNB angles and Ricketts' NBa-PtGn angle were measured to determine the skeletal malocclusion and the vertical pattern. Nine SNPs in BMP2, BMP4, SMAD6, RUNX2, WNT3A and WNT11 were genotyped. Chi-squared test was used to compare genotypes among the groups. Multifactor dimensionality reduction (MDR) and binary logistic regression analysis, both using gender and age as co-variables, were also used. We performed Bonferroni correction for multiple testing.

RESULTS

Significant associations at P < .05 were observed for SNPs rs1005464 (P = .042) and rs235768 (P = .021) in BMP2 with mandibular retrognathism and for rs59983488 (RUNX2) with maxillary protrusion (P = .04) as well as for rs708111 (WNT3A) with skeletal class III (P = .02; dominant model), rs1533767 (WNT11) with a brachyfacial skeletal pattern (P = .01, OR = 0.10; dominant model) and for rs3934908 (SMAD6) with prognathism (P = .02; recessive model). After the Bonferroni correction, none of the SNPs remained associated. The MDR predicted some interaction for skeletal class II, dolichofacial and brachyfacial phenotypes.

CONCLUSION

Our results suggest that SNPs in BMP2, BMP4, SMAD6, RUNX2, WNT3A and WNT11 could be involved in the aetiology of sagittal and vertical malocclusions.

Association between forkhead box E1 polymorphisms and risk of non-syndromic cleft lip with or without cleft palate: A meta-analysis

OBJECTIVE

The purpose of the present work was to investigate the association between forkhead box E1 (FOXE1) and the risk of non-syndromic cleft lip with or without cleft palate (NSCL/P).

MATERIALS AND METHODS

Relevant studies were searched in several professional databases up to 31 July 2019. The pooled odds ratios (ORs) and 95% confidence intervals (95% CIs) were calculated using a fixed-effect model or a random-effect model to analyse the relationship between FOXE1 polymorphisms and NSCL/P.

RESULTS

A total of four single nucleotide polymorphisms (SNPs), including rs3758249, rs4460498, rs1443434 and rs10217225, were analysed. The overall findings showed that FOXE1 rs4460498 was statistically associated with NSCL/P (including cleft lip with or without cleft palate (CL/P) and cleft palate only (CPO)). Genotypes CC and CT of rs4460498 were significantly more closely correlated with NSCL/P (including CL/P and CPO) than genotype TT (NSCL/P: TT vs CC, OR = 0.630, P = .000; TT vs TC + CC, OR = 0.775, P = .020; CL/P: TT vs CC, OR = 0.664, P = .000; TT vs TC + CC, OR = 0.738, P = .006. CPO: TT vs CC, OR = 0.761, P = .027; TT vs TC + CC, OR = 0.792, P = .045). For rs10217225, only the TT genotype might have contributed to the elevated risk of CL/P (TT vs CC OR = 2.236, P = .000). The other FOXE1 polymorphisms were not associated with NSCLP, CL/P or CPO.

CONCLUSION

The meta-analysis provided confirmation that the polymorphism of FOXE1 rs10217225 was correlated with an increased risk of CL/P, and the polymorphism of FOXE1 rs4460498 was a protective factor for NSCL/P, including CLP and CPO.