Moreover, in our present results, all variants in chromatin modification genes and CYP26B1 gene are concurrent with other variants in neural tube-related genes (Fig. 3c), so we hypothesize that the craniorachischisis phenotype is also resulted from the synergistic effect of NTC-related polygenic functional variants. The gene discussed is CYP26B1; the disease is craniorachischisis.