Brazilian multicenter study of association between polymorphisms in CRISPLD2 and JARID2 and non-syndromic oral clefts

Ethnic Differences in the Brazilian Population Influence the Impact of BMP4 Genetic Variants on Susceptibility of Nonsyndromic Orofacial Clefts

OBJECTIVE

The study evaluated the association of BMP4 tag-SNPs and SNP-SNP interactions involving genes active by BMP4 pathway during craniofacial development in the susceptibility of nonsyndromic orofacial clefts (NSOC) in the Brazilian population.

DESIGN

Case-control study.

SETTING

Brazilian Oral Cleft Group.

PARTICIPANTS

The study included 881 healthy controls and 800 patients with different types of NSOC: 232 with cleft lip only (NSCLO), 568 with cleft lip and palate (NSCLP), and 274 with cleft palate only (NSCPO).

INTERVENTIONS

The genomic DNA was genotyped with allelic discrimination assays for five BMP4 tag-SNPs (rs11623717, rs17563, rs2071047, rs2761887 and rs4898820), and analyzed their allelic and genotypic associations using multiple logistic regression. The interactions of these variants with genes involved in the BMP4 signaling pathway, including FGFR1, GREM1, NOG, VAX1 and the 4p16.2 locus, were explored.

MAIN OUTCOME MEASURES

BMP4 variants in the NSOC risk.

RESULTS

Although only nominal p values were identified when the whole sample was considered, subgroup analysis including the patients with high African genomic ancestry showed significant associations of rs2761887 with risk for nonsyndromic cleft lip with or without cleft palate (NSCL  ±  P)[(ORhom: 2.16; 95% CI: 1.21-3.85; p  =  0.01) and (ORrec: 2.05; 95% CI: 1.21-3.47; p  =  0.006)]. Thirteen significant SNP-SNP interactions involving BMP4 and the SNPs at FGFR1, GREM1, NOG and VAX1 and at locus 4p16.2 for increased risk of NSCL  ±  P were identified.

CONCLUSIONS

Our results demonstrate an increased risk of NSCL  ±  P in Brazilian individuals with enrichment of African ancestry in the presence of the BMP4 rs2762887 polymorphism, and reveal relevant genetic contribution of SNP-SNP epistatic interactions involving BMP4 variants to NSCL  ±  P risk.

Concurrent de novo ZFHX4 variant and 16q24.1 deletion in a patient with orofacial clefting; a potential role of ZFHX4 and USP10

A girl with a unilateral cleft lip, alveolus and palate, tooth agenesis, and mild dysmorphic features, without a specific underlying syndrome diagnosis, was genotypically characterized and phenotypically described. Cleft gene panel analysis, single-nucleotide polymorphism (SNP) array, whole genome sequencing (WGS), whole exome sequencing, and quantitative PCR (Q-PCR) analysis were used as diagnostic tests. SNP array revealed a maternal deletion at 16q24.1, encompassing the cleft candidate gene USP10. WES revealed an additional de novo Loss-of-Function variant (p.(Asn838fs)) in the Zinc-Finger-Homeobox-4 (ZFHX4) gene. Q-PCR was performed to explore the effect of the ZFHX4 variant and the deletion in 16q24.1. The mRNA expression of a selection of putative target genes involved in orofacial clefting showed a lowered expression of USP10 (52%), CRISPLD2 (31%), and CRISPLD1 (1%) compared to the control. IRF6 showed no difference in gene expression. This case supports ZFHX4 as a novel cleft gene and suggests USP10 may contribute to the etiology of orofacial clefts in humans.

Brazilian Multiethnic Association Study of Genetic Variant Interactions among FOS, CASP8, MMP2 and CRISPLD2 in the Risk of Nonsyndromic Cleft Lip with or without Cleft Palate

Associations of CRISPLD2 (cysteine-rich secretory protein LCCL domain containing 2) and genes belonging to its activation pathway, including FOS (Fos proto-oncogene), CASP8 (caspase 8) and MMP2 (matrix metalloproteinase 2), with nonsyndromic orofacial cleft risk, have been reported, but the results are yet unclear. The aim of this study was to evaluate single nucleotide polymorphisms (SNPs) in FOS, CASP8 and MMP2 and to determine their SNP-SNP interactions with CRISPLD2 variants in the risk of nonsyndromic cleft lip with or without cleft palate (NSCL±P) in the Brazilian population. The SNPs rs1046117 (FOS), rs3769825 (CASP8) and rs243836 (MMP2) were genotyped using TaqMan allelic discrimination assays in a case-control sample containing 801 NSCL±P patients (233 nonsyndromic cleft lip only (NSCLO) and 568 nonsyndromic cleft lip and palate (NSCLP)) and 881 healthy controls via logistic regression analysis adjusted for the effects of sex and genomic ancestry proportions with a multiple comparison p value set at ≤0.01. SNP-SNP interactions with rs1546124, rs8061351, rs2326398 and rs4783099 in CRISPLD2 were performed with the model-based multifactor dimensionality reduction test complemented with a 1000 permutation-based strategy. Although the association between FOS rs1046117 and risk of NSCL±P reached only nominal p values, NSCLO risk was significantly higher in carriers of the FOS rs1046117 C allele (OR: 1.28, 95% CI: 1.10-1.64, p = 0.004), TC heterozygous genotype (OR: 1.59, 95% CI: 1.16-2.18, p = 0.003), and in the dominant model (OR: 1.50, 95% CI: 1.10-2.02, p = 0.007). Individually, no significant associations between cleft risk and the SNPs in CASP8 and MMP2 were observed. SNP-SNP interactions involving CRISPLD2 variants and rs1046117 (FOS), rs3769825 (CASP8) and rs243836 (MMP2) yielded several significant p values, mostly driven by FOS rs1046117 and CASP8 rs3769825 in NSCL±P, FOS rs1046117 in NSCLO and CRISPLD2 rs8061351 in NSCLP. Our study is the first in the Brazilian population to reveal the association of FOS rs1046117 with NSCLO risk, and to support that CRISPLD2, CASP8, FOS and MMP2 interactions may be related to the pathogenesis of this common craniofacial malformation.

Phenotypic expansion and variable expressivity in individuals with JARID2-related intellectual disability: A case series

Loss of function variants in JARID2 were recently reported in 16 patients with a neurodevelopmental disorder characterized by delays, intellectual and learning disability, autism, behavioral abnormalities, and dysmorphic features. Most cases were de novo, with only one variant inherited from an affected parent. Here, we present 7 additional individuals from 5 families with pathogenic or likely pathogenic JARID2 variants, confirming this gene-disease association and highlighting palatal abnormalities and heart defect as part of the phenotype. In addition, we report inheritance of JARID2 variants from mildly affected parents, demonstrating the variable expressivity of the disease. We also note the high prevalence of intragenic JARID2 copy number variants, emphasizing the importance of exon-level analysis.

Evaluation of genome-wide association signals for nonsyndromic cleft lip with or without cleft palate in a multiethnic Brazilian population

OBJECTIVE

To evaluate previous nonsyndromic cleft lip with or without cleft palate (NSCL±P) associated signals in 4p16.2, 8p11.23, 12q13.13, 12q13.2 and 17q21.32 in a multiethnic Brazilian cohort.

DESIGN

The single nucleotide polymorphisms (SNPs) rs34246903 in 4p16.2, rs13317 in 8p11.23 (FGFR1, fibroblast growth factor receptor 1), rs3741442 in 12q13.13, rs705704 in 12q13.2 and rs4968247 in 17q21.32 were genotyped with TaqMan allelic discrimination assays in a case-control sample including 801 NSCL±P patients [233 nonsyndromic cleft lip (NSCLO) and 568 nonsyndromic cleft lip and palate (NSCLP)] and 881 healthy controls. Multiple logistic regression analyses, considering sex and genomic ancestry as covariates, were conducted, and the p value was adjusted with Bonferroni multiple correction testing (p ≤ 0.01).

RESULTS

Although several associations were identified, those that resisted the multiple correction testing involved the alleles and genotypes of rs34246903 and rs13317. The NSCLO group had a lower frequency of the minor C allele of rs34246903 compared to controls, giving an odds ratio (OR) of 0.74 [95% confidence interval (CI): 0.59-0.93, p = 0.01]. The rs34246903 CC genotype (homozygous) and the recessive model revealed significant protective associations with NSCLO, yielding ORs of 0.50 (95% CI: 0.29-0.85, p = 0.005) and 0.55 (95% CI: 0.33-0.93, p = 0.01) respectively. The presence of C variant allele of rs13317 (OR: 0.81, 95% CI: 0.69-0.96, p = 0.01) as well the TC genotype (OR: 0.77, 95% CI: 0.62-0.94, p = 0.01) and the dominant model (OR: 0.77, 95% CI: 0.63-0.94, p = 0.009) showed significant associations with reduced risk of NSCL±P.

CONCLUSION

Our study is the first to support the association of rs34246903 (4p16.2) with NSCLO and rs13317 within FGFR1 with NSCL±P in the highly admixed Brazilian population. Further studies are needed to determine the functionality of those SNPs or to identify the causal markers in linkage disequilibrium with those susceptibility markers.

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.

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.

JARID2 haploinsufficiency is associated with a clinically distinct neurodevelopmental syndrome

PURPOSE

JARID2, located on chromosome 6p22.3, is a regulator of histone methyltransferase complexes that is expressed in human neurons. So far, 13 individuals sharing clinical features including intellectual disability (ID) were reported with de novo heterozygous deletions in 6p22-p24 encompassing the full length JARID2 gene (OMIM 601594). However, all published individuals to date have a deletion of at least one other adjoining gene, making it difficult to determine if JARID2 is the critical gene responsible for the shared features. We aim to confirm JARID2 as a human disease gene and further elucidate the associated clinical phenotype.

METHODS

Chromosome microarray analysis, exome sequencing, and an online matching platform (GeneMatcher) were used to identify individuals with single-nucleotide variants or deletions involving JARID2.

RESULTS

We report 16 individuals in 15 families with a deletion or single-nucleotide variant in JARID2. Several of these variants are likely to result in haploinsufficiency due to nonsense-mediated messenger RNA (mRNA) decay. All individuals have developmental delay and/or ID and share some overlapping clinical characteristics such as facial features with those who have larger deletions involving JARID2.

CONCLUSION

We report that JARID2 haploinsufficiency leads to a clinically distinct neurodevelopmental syndrome, thus establishing gene-disease validity for the purpose of diagnostic reporting.

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

OBJECTIVE

To investigate the association of single-nucleotide polymorphisms (SNP) in grainyhead-like 3 (GRHL3) and to verify its possible interactions with others genes responsible for craniofacial development in the risk of non-syndromic oral cleft (NSOC).

METHODS

Applying TaqMan allelic discrimination assays, we evaluated GRHL3 SNPs (rs10903078, rs41268753, and rs4648975) in an ancestry-structured case-control sample composed of 1,127 Brazilian participants [272 non-syndromic cleft palate only (NSCPO), 242 non-syndromic cleft lip only (NSCLO), 319 non-syndromic cleft lip and palate (NSCLP), and 294 healthy controls]. Additionally, SNP-SNP interactions of GRHL3 and previously reported variants in FAM49A, FOXE1, NTN1, and VAX1 were verified in non-syndromic cleft lip with or without cleft palate (NSCL ± P). To eliminate false-positive associations, Bonferroni correction or 1,000 permutation method was applied.

RESULTS

The multiple logistic regression analysis showed that the CC genotype of rs10903078 (p = .03) and the haplotype C-C formed by the SNPs rs10903078 and rs41268753 (p = .04) were associated with NSCLO, but the p-values did not withstand Bonferroni correction. However, SNP-SNP test revealed significant interactions between GRHL3 SNPs and FAM49A (rs7552), FOXE1 (rs3758249), VAX1 (rs7078160 and rs751231), and NTN1 (rs9891446).

CONCLUSIONS

Our results confirm the importance of GRHL3 and its interactions with previously NSOC-associated genes, including FAM49A, FOXE1, NTN1, and VAX1, in the pathogenesis of NSOC in the Brazilian population.

A systematic literature review on the European, African and Amerindian genetic ancestry components on Brazilian health outcomes

The variables such as race, skin colour and ethnicity have become intensely discussed in medicine research, as a response to the rising debate over the importance of the ethnic-racial dimension in the scope of health-disease processes. The aim of this study was to identify the European (EUR), African (AFR) and Amerindian (AMR) ancestries on Brazilian health outcomes through a systematic literature review. This study was carried out by searching in three electronic databases, for studies published between 2005 and 2017. A total of 13 papers were eligible. The search identified the following health outcomes: visceral leishmaniosis, malaria, Alzheimer’s disease, neuromyelitis optica, multiple sclerosis, prostate cancer, non-syndromic cleft lip/palate, chronic heart failure, sickle cell disease, primary congenital glaucoma, preterm labour, preterm premature rupture of membranes, systemic lupus erythematosus and type 1 diabetes mellitus. Research paper assessments were guided by the STROBE instrument, and agreements between results were determined by comparing the points attributed by two authors. Increased EUR ancestry was identified from preterm labour (PTL), type 1 diabetes (T1D) and non-syndromic cleft lip with or without cleft palate (NSCL), as well as in patients presenting aggressive prostate cancer prognoses. On the other hand, the highest AFR ancestral component was verified from systemic lupus erythematosus (SLE) and primary congenital glaucoma (PCG) cases, presenting worse prognoses. AMR ancestry may be a protective factor in the development of Alzheimer’s disease (AD). The worst hemodynamic parameters in cases of heart failure (HF) were identified among individuals with greater AMR and AFR ancestry indices.

Interactions between superoxide dismutase and paraoxonase polymorphic variants in nonsyndromic cleft lip with or without cleft palate in the Brazilian population

During development, oxidative stress is hypothesized to mediate embryotoxicity, which may be intensified by exposition to environmental factors and by genetic variations in the enzymes involved in protecting cells from these damaging effects, including superoxide dismutase (SOD) and paraoxonase (PON). The aim of this study was to evaluate the influence of single-nucleotide polymorphisms (SNP) in genes associated with the neutralization of oxidative stress (SOD and PON family members) in the risk of nonsyndromic oral cleft in the Brazilian population. Initially, we tested for association between 28 SNP in SOD1, SOD2, SOD3, PON1, PON2, and PON3 among 325 nonsyndromic cleft lip with or without cleft palate (NSCL±P) case-parent trios. Multiple logistic regression analyses were used to explore gene, GxG and GxE, involving factors that induce oxidative stress accumulation during pregnancy. Signals still significant after both Bonferroni correction and in permutation test were subsequently confirmed in an ancestry-structured case-control analysis with 722 NSCL±P and 866 controls from the same population. In the trio sample, transmission disequilibrium test (TDT) (allele and haplotype) and GxE analysis showed no significant associations, but multiple pairwise GxG interactions involving 10 SNP in PON1, PON2, and PON3 were detected and further examined in the case-control sample. The PON1 rs2237583 and PON2 rs17166879 yielded significant evidence of SNP-SNP interactions after adjustment for multiple tests (both Bonferroni correction and 10,000 permutation test). The C allele and the CT genotype of PON1 rs2237583 were associated with significant protective effects against NSCL±P, while rs3917490 showed a significant association only in the sample composed of patients displaying high African ancestry. Our results reveal associations between rs2237583 and rs3917490 in PON1 and GxG interactions containing rs2237583 and rs17166879 with the susceptibility of NSCL±P in the Brazilian population. Furthermore, this study underlines the recent tendency of taking into account potential GxG interactions to clarify the underlying mechanisms associated with the etiology of this common malformation. Environ. Mol. Mutagen. 60: 185-196, 2019. © 2018 Wiley Periodicals, Inc.

Understanding the participation of GREM1 polymorphisms in nonsyndromic cleft lip with or without cleft palate in the Brazilian population

BACKGROUND

GREM1, which encodes Gremlin 1, an antagonist of bone morphogenic proteins with effects on proliferation and apoptosis, has been considered a candidate gene for nonsyndromic cleft lip with or without cleft palate (NSCL±P). In this study, we investigated potential associations of single nucleotide polymorphisms (SNP) in GREM1 and NSCL±P risk in the Brazilian population. Additionally, SNP-SNP interactions of GREM1 with previously reported rs1880646 variant in NTN1 (netrin 1), a gene also responsible for apoptotic phenotypes were verified.

METHODS

Applying Taqman allelic discrimination assays, we evaluated the variants rs16969681, rs16969816, rs16969862, and rs1258763 in 325 case-parent trios and in 1,588 isolated samples in a case-control study. Allelic and genotypic analyses, as well as interaction tests assessing gene-environmental factor (GxE) and SNP-SNP interaction with rs1880646 variant in NTN1, were performed based on logistic regression analysis adjusted for the effects of gender and genomic ancestry proportions.

RESULTS

The risk alleles of all SNP were undertransmitted in NSCL±P trios, though the case-control analysis confirmed only the association with rs16969862 alleles (OR: 0.78, 95% CI: 0.63-0.96, p = .02). The GxE interaction analysis revealed a significant interaction between maternal environmental contact with agrotoxics and rs16969816 (OR: 0.25, 95% CI: 0.08-0.74, p = .01), and pairwise interaction test with NTN1 rs1880646 yielded significant p values in the 1,000 permutation test for rs16969681, rs16969816, and rs16969862.

CONCLUSION

The GREM1 is involved in the etiology of NSCL±P in the Brazilian population and reveal that the interaction between GREM1 and NTN1 may be related with the pathogenesis of this common craniofacial malformation.

Dental genetics in Brazil: Where we are

Dentistry constitutes the basic nucleus of professionals of higher level of health in Brazil with one of the largest concentrations of dentists per capita in the world. However, the genetic in dentistry in Brazil is explored, basically, in research field. Future actions need to be performed in order to deep the whole knowledge about diagnosis and treatment of diseases with genetic basis in dentistry, in Brazil.

Association Between CRISPLD2 Polymorphisms and the Risk of Nonsyndromic Clefts of the Lip and/or Palate: A Meta-analysis

OBJECTIVE

Nonsyndromic clefts of the lip and/or palate (NSCL/P) are one of the most common polygenic diseases. Recently, many studies focused on the association between CRISPLD2 polymorphisms and NSCL/P risk. However, some studies have shown opposite results. In this study, meta-analysis was used to confirm whether CRISPLD2 polymorphism was associated with NSCL/P, and the possible mechanism between CRISPLD2 and NSCL/P was explored.

METHODS

Relevant studies were conducted on PubMed, Ovid, EBSCO, CINAHL, FMRS, Web of Science, CNKI, and Wanfang databases from their inception up to June 31, 2016. Review Manager 5.0.24 was used to analyze whether CRISPLD2 polymorphism was involved in NSCL/P by pooling odds ratios (ORs) and 95% confidence intervals (CIs). Potential publication bias was evaluated by visual inspection of the funnel plot.

RESULTS

CRISPLD2 rs4783099 was associated with cleft lip and/or palate (CL/P) statistically (OR = 3.18, P < .01). Compared to genotype TT, genotypes CC and CT were correlated significantly (OR = 2.04, P = .04) with CL/P. No evidence showed an association between genetic variation at the CRISPLD2 locus and cleft palate only (CP).

CONCLUSION

The polymorphism of CRISPLD2 rs4783099 is correlated with an increased risk of CL/P.

Possible effect of SNAIL family transcriptional repressor 1 polymorphisms in non-syndromic cleft lip with or without cleft palate

OBJECTIVE

Orofacial development is a complex process subjected to failure impairing. Indeed, the cleft of the lip and/or of the palate is among the most frequent inborn malformations. The JARID2 gene has been suggested to be involved in non-syndromic cleft lip with or without cleft palate (nsCL/P) etiology. JARID2 interacts with the polycomb repressive complex 2 (PRC2) in regulating the expression patterns of developmental genes by modifying the chromatin state.

MATERIALS AND METHODS

Genes coding for the PRC2 components, as well as other genes active in cell differentiation and embryonic development, were selected for a family-based association study to verify their involvement in nsCL/P. A total of 632 families from Italy and Asia participated to the study.

RESULTS

Evidence of allelic association was found with polymorphisms of SNAI1; in particular, the rs16995010-G allele was undertransmitted to the nsCL/P cases [P = 0.004, odds ratio = 0.69 (95% C.I. 0.54-0.89)]. However, the adjusted significance value corrected for all the performed tests was P = 0.051.

CONCLUSIONS

The findings emerging by the present study suggest for the first time an involvement of SNAI1 in the nsCL/P onset.

CLINICAL RELEVANCE

Interestingly, SNAI1 is known to promote epithelial to mesenchymal transition by repressing E-cadherin expression, but it needs an intact PRC2 to act this function. Alterations of this process could contribute to the complex etiology of nsCL/P.