Myelomeningocele genotype-phenotype correlation findings in cilia, HH, PCP, and WNT signaling pathways.

Abstract

Myelomeningocele (MMC) is the most severe and frequent type of spina bifida. Its etiology remains poorly understood. The Hedgehog (Hh), Wnt, and planar cell polarity (PCP) signaling pathways are essential for normal tube closure, needing a structural-functional cilium for its adequate function. The present study aimed to investigate the impact of different gene variants (GV) from those pathways on MMC genotype-subphenotype correlations. The study comprised 500 MMC trios and 500 controls, from 16 Telethon centers of 16 Mexican states. Thirty-four GVs of 29 genes from cilia, Hh, PCP, and Wnt pathways, were analyzed, by an Illumina on design microarray. The total sample (T-MMC) was stratified in High-MMC (H-MMC) when thoracic and Low-MMC (L-MMC) when lumbar-sacral vertebrae affected. STATA/SE-12.1 and PLINK software were used for allelic association, TDT, and gene-gene interaction (GGI) analyses, considering p value <.01 as statistically significant differences (SSD). Association analysis showed SSD for COBL-rs10230120, DVL2-rs2074216, PLCB4-rs6077510 GVs in T-MMC and L-MMC, and VANGL2-rs120886448 in T-MMC and H-MMC, and INVS-rs7024375 exclusively in L-MMC. TDT assay showed SSD preferential transmissions of C2CD3-rs826058 in H-MMC, and LRP5-rs3736228, and BBS2-rs1373 in L-MMC. Statistically significant GGI was observed in four in T-MMC, four completely different in L-MMC, and one in H-MMC. Interestingly, no one repeated in subphenotypes. Our results support an association of GVs in Hh, Wnt, PCP, and cilia pathways, with MMC occurrence location, although further validation is needed. Furthermore, present results show a distinctive panel of gene-variants in H-MMC and LMMC subphenotypes, suggesting a feasible genotype-phenotype correlation.