Joint testing of genotypic and gene-environment interaction identified novel association for BMP4 with non-syndromic CL/P in an Asian population using data from an International Cleft Consortium

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

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

Functional Diversification of SRSF Protein Kinase to Control Ubiquitin-Dependent Neurodevelopmental Signaling

Conserved protein kinases with core cellular functions have been frequently redeployed during metazoan evolution to regulate specialized developmental processes. The Ser/Arg (SR)-rich splicing factor (SRSF) protein kinase (SRPK), which is implicated in splicing regulation, is one such conserved eukaryotic kinase. Surprisingly, we show that SRPK has acquired the capacity to control a neurodevelopmental ubiquitin signaling pathway. In mammalian embryonic stem cells and cultured neurons, SRPK phosphorylates Ser-Arg motifs in RNF12/RLIM, a key developmental E3 ubiquitin ligase that is mutated in an intellectual disability syndrome. Processive phosphorylation by SRPK stimulates RNF12-dependent ubiquitylation of nuclear transcription factor substrates, thereby acting to restrain a neural gene expression program that is aberrantly expressed in intellectual disability. SRPK family genes are also mutated in intellectual disability disorders, and patient-derived SRPK point mutations impair RNF12 phosphorylation. Our data reveal unappreciated functional diversification of SRPK to regulate ubiquitin signaling that ensures correct regulation of neurodevelopmental gene expression.

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.

Accumulation of rare coding variants in genes implicated in risk of human cleft lip with or without cleft palate

Cleft lip with/without cleft palate (CLP) is a common craniofacial malformation with complex etiologies, reflecting both genetic and environmental factors. Most of the suspected genetic risk for CLP has yet to be identified. To further classify risk loci and estimate the contribution of rare variants, we sequenced the exons in 49 candidate genes in 323 CLP cases and 211 nonmalformed controls. Our findings indicated that rare, protein-altering variants displayed markedly higher burdens in CLP cases at relevant loci. First, putative loss-of-function mutations (nonsense, frameshift) were significantly enriched among cases: 13 of 323 cases (~4%) harbored such alleles within these 49 genes, versus one such change in controls (p = 0.01). Second, in gene-level analyses, the burden of rare alleles showed greater case-association for several genes previously implicated in cleft risk. For example, BHMT displayed a 10-fold increase in protein-altering variants in CLP cases (p = .03), including multiple case occurrences of a rare frameshift mutation (K400 fs). Other loci with greater rare, coding allele burdens in cases were in signaling pathways relevant to craniofacial development (WNT9B, BMP4, BMPR1B) as well as the methionine cycle (MTRR). We conclude that rare coding variants may confer risk for isolated CLP.

Evaluation of the presence of MMP-2, TIMP-2, BMP2/4, and TGFβ3 in the facial tissue of children with cleft lip and palate

Cleft lip and palate (CLP) is the most common defect affecting the face. The treatment consists of surgical reconstruction of the anatomical structures of the cleft. Part of the surgical treatment is reconstruction of the alveolar bone by means of autogenic bone grafting (osteoplasty). This study aimed to evaluate the levels of expression of extracellular matrix remodeling factors in the facial tissue of children with a complete unilateral (CU) and a complete bilateral (CB) CLP to assess whether the wound healing process is adequate. Twenty-two CLP patients were enrolled in this study. Tissue samples were collected during alveolar osteoplasty for unilateral (n = 12) or bilateral (n = 10) cleft palate, (age range from 6 years 8 months to 12 years 2 months). Control material was obtained in the case of tooth extraction (age range from 6 years 9 months to 14 years 5 months). Immunohistochemistry was used to assess the levels of matrix metalloproteinase-2 (MMP-2), tissue inhibitor of metalloproteinase-2 (TIMP-2), bone morphogenetic proteins 2 and 4 (BMP2/4), and transforming growth factor β3 (TGFβ3). Numbers of positively stained cells were graded semi-quantitatively. Data were analysed using the Kraskel-Wallis rank test and the Bonferroni correction. The total number of MMP2-positive cells was significantly lower in the CBCLP and in the control group than in the CUCLP (p < 0.001 after the Bonferroni correction). The total number of TIMP2-positive cells was significantly higher in the CUCLP than in the CBCLP and in the control group (p < 0.001; p < 0.003 after the Bonferroni correction). The overall number of BMP2/4, TGFβ3-positive cells was significantly higher in the CUCLP than in the CBCLP and in the control group (p < 0.001 after the Bonferroni correction). The decrease of the relative amount of statistically significant BMP2/4, TGFβ3, MMP-2, TIMP-2 containing bone cells in CBCLP patients identifies affected alveolar bone regeneration and remodeling process.

Association between BMP4 gene polymorphisms and cleft lip with or without cleft palate in a population from South China

OBJECTIVE

Previous studies have suggested an association between several polymorphisms of the BMP4 gene and susceptibility to non-syndromic cleft lip with or without cleft palate (NSCL/P) in various populations. However, this association may vary according to ethnic group and the form of NSCL/P. This study analyzed the association between the BMP4 gene polymorphisms rs762642, rs17563, and rs10130587 with the risk of cleft lip only (CLO), cleft palate only (CPO), and cleft lip with palate (CLP) in a population from South China.

METHODS

This case-control study included 165 patients with NSCL/P (53 patients with CPO, 52 with CLO, and 60 with CLP) and 52 healthy volunteers. Peripheral blood samples were collected from all subjects to genotype the rs762642, rs17563, and rs10130587 polymorphisms by direct sequencing. Genotype and allelic frequencies of these polymorphisms were compared between healthy volunteers and patients with various forms of NSCL/P.

RESULTS

The genotype and allelic frequencies of rs762642 differed significantly between subgroups (CPO and CLP) and normal controls, whereas a significant difference was observed only in the CLO subgroup for the rs17563 polymorphism and in the CLO and CLP groups for the rs10130587 polymorphism. In addition, we identified a novel association of a BMP4 gene polymorphism, which was in linkage disequilibrium with the rs10130587 polymorphism, with CLO and CLP.

CONCLUSION

The BMP4 gene polymorphisms rs762642, rs17563, and rs10130587 exhibit different associations with different forms of NSCL/P, suggesting that different forms of NSCL/P may have different etiologies.

Suicidal ideation in systemic lupus erythematosus: NR2A gene polymorphism, clinical and psychosocial factors

Background Systemic lupus erythematosus (SLE) patients are a high-risk population for suicide. Glutamatergic neurosystem genes have been implicated in the neurobiology of depression in SLE and suicidal behaviour in general. However, the role of glutamate receptor gene polymorphisms in suicidal behaviour among SLE patients remains unclear in the context of established clinical and psychosocial factors. We aimed to investigate the association of NR2A gene polymorphism with suicidal ideation in SLE while accounting for the interaction between clinical and psychosocial factors. Methods A total of 130 SLE patients were assessed for mood disorders (MINI International Neuropsychiatric Interview), severity of depression (Patient Health Questionnaire-9), suicidal behaviour (Columbia-Suicide Severity Rating Scale), socio-occupational functioning (Work and Social Adjustment Scale), recent life events (Social Readjustment Rating Scale) and lupus disease activity (SELENA-SLE Disease Activity Index). Eighty-six out of the 130 study participants consented for NR2A genotyping. Results Multivariable logistic regression showed nominal significance for the interaction effect between the NR2A rs2072450 AC genotype and higher severity of socio-occupational impairment with lifetime suicidal ideation in SLE patients ( p = 0.038, odds ratio = 1.364, 95% confidence interval = 1.018-1.827). However, only the association between lifetime mood disorder and lifetime suicidal ideation remained significant after Bonferroni correction ( p < 0.001, odds ratio = 33.834, 95% confidence interval = 7.624-150.138). Conclusions Lifetime mood disorder emerged as a more significant factor for suicidal ideation in SLE compared with NR2A gene polymorphism main and interaction effects. Clinical implications include identification and treatment of mood disorders as an early intervention for suicidal behaviour in SLE. More adequately-powered gene-environment interaction studies are required in the future to clarify the role of glutamate receptor gene polymorphisms in the risk stratification of suicidal behaviour among SLE patients.

Methylome analysis for spina bifida shows SOX18 hypomethylation as a risk factor with evidence for a complex (epi)genetic interplay to affect neural tube development

Background

Neural tube defects (NTDs) are severe congenital malformations that arise from failure of neurulation during early embryonic development. The molecular basis underlying most human NTDs still remains largely unknown. Based on the hypothesis that folic acid prevents NTDs by stimulating methylation reactions, DNA methylation changes could play a role in NTDs. We performed a methylome analysis for patients with myelomeningocele (MMC). Using a candidate CpG analysis for HOX genes, a significant association between HOXB7 hypomethylation and MMC was found.

Methods

In the current study, we analyzed leukocyte methylome data of ten patients with MMC and six controls using Illumina Methylation Analyzer and WateRmelon R-packages and performed validation studies using larger MMC and control cohorts with Sequenom EpiTYPER.

Results

The methylome analysis showed 75 CpGs in 45 genes that are significantly differentially methylated in MMC patients. CpG-specific methylation differences were next replicated for the top six candidate genes ABAT, CNTNAP1, SLC1A6, SNED1, SOX18, and TEPP but only for the SOX18 locus a significant overall hypomethylation was observed (P value = 0.0003). Chemically induced DNA demethylation in HEK cells resulted in SOX18 hypomethylation and increased expression. Injection of sox18 mRNA in zebrafish resulted in abnormal neural tube formation. Quantification of DNA methylation for the SOX18 locus was also determined for five families where parents had normal methylation values compared to significant lower values for both the MMC as their non-affected child. SOX18 methylation studies were performed for a MMC patient with a paternally inherited chromosomal deletion that includes BMP4. The patient showed extreme SOX18 hypomethylation similar to his healthy mother while his father had normal methylation values.

Conclusions

This is the first genome-wide methylation study in leukocytes for patients with NTDs. We report SOX18 as a novel MMC risk gene but our findings also suggest that SOX18 hypomethylation must interplay with environmental and (epi)genetic factors to cause NTDs. Further studies are needed that combine methylome data with next-generation sequencing approaches to unravel NTD etiology.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-016-0272-8) contains supplementary material, which is available to authorized users.

DNA Methylation in Newborns and Maternal Smoking in Pregnancy: Genome-wide Consortium Meta-analysis

Epigenetic modifications, including DNA methylation, represent a potential mechanism for environmental impacts on human disease. Maternal smoking in pregnancy remains an important public health problem that impacts child health in a myriad of ways and has potential lifelong consequences. The mechanisms are largely unknown, but epigenetics most likely plays a role. We formed the Pregnancy And Childhood Epigenetics (PACE) consortium and meta-analyzed, across 13 cohorts (n = 6,685), the association between maternal smoking in pregnancy and newborn blood DNA methylation at over 450,000 CpG sites (CpGs) by using the Illumina 450K BeadChip. Over 6,000 CpGs were differentially methylated in relation to maternal smoking at genome-wide statistical significance (false discovery rate, 5%), including 2,965 CpGs corresponding to 2,017 genes not previously related to smoking and methylation in either newborns or adults. Several genes are relevant to diseases that can be caused by maternal smoking (e.g., orofacial clefts and asthma) or adult smoking (e.g., certain cancers). A number of differentially methylated CpGs were associated with gene expression. We observed enrichment in pathways and processes critical to development. In older children (5 cohorts, n = 3,187), 100% of CpGs gave at least nominal levels of significance, far more than expected by chance (p value < 2.2 × 10(-16)). Results were robust to different normalization methods used across studies and cell type adjustment. In this large scale meta-analysis of methylation data, we identified numerous loci involved in response to maternal smoking in pregnancy with persistence into later childhood and provide insights into mechanisms underlying effects of this important exposure.

Common mechanisms in development and disease: BMP signaling in craniofacial development

BMP signaling is one of the key pathways regulating craniofacial development. It is involved in the early patterning of the head, the development of cranial neural crest cells, and facial patterning. It regulates development of its mineralized structures, such as cranial bones, maxilla, mandible, palate, and teeth. Targeted mutations in the mouse have been instrumental to delineate the functional involvement of this signaling network in different aspects of craniofacial development. Gene polymorphisms and mutations in BMP pathway genes have been associated with various non-syndromic and syndromic human craniofacial malformations. The identification of intricate cellular interactions and underlying molecular pathways illustrate the importance of local fine-regulation of Bmp signaling to control proliferation, apoptosis, epithelial-mesenchymal interactions, and stem/progenitor differentiation during craniofacial development. Thus, BMP signaling contributes both to shape and functionality of our facial features. BMP signaling also regulates postnatal craniofacial growth and is associated with dental structures life-long. A more detailed understanding of BMP function in growth, homeostasis, and repair of postnatal craniofacial tissues will contribute to our ability to rationally manipulate this signaling network in the context of tissue engineering.

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.